ana::SBNOsc::Chi2Sensitivity::EventSample Class Reference

Public Member Functions
	EventSample (const fhicl::ParameterSet &config)
std::vector< double >	Oscillate (double sinth, double dm2) const
std::vector< double >	Signal () const
std::vector< double >	Background () const

Public Attributes
std::string	fName
double	fBaseline
	baseline (along z) in cm More...
double	fBeamCenterX
	Center of beam in detector coordinates in x-dir in cm. More...
double	fBeamCenterY
	Center of beam in detector coordinates in y-dir in cm. More...
double	fBeamFrontZ
	Front face of detector along z-dir in cm. More...
std::array< double, 2 >	fXlim
	Detector size in cm. More...
std::array< double, 2 >	fYlim
	Detector size in cm. More...
std::array< double, 2 >	fZlim
	Detector size in cm. More...
double	fScalePOT
	Factor for POT (etc.) scaling. More...
double	fPOT
int	fOscType
	Oscilaltion type: 0 == None, 1 == numu -> nue, 2 == numu -> numu. More...
std::vector< double >	fEnergyBinScale
TH1D *	fBkgCounts
	Background count Histogram. More...
TH3D *	fSignalCounts
	Signal Count Histogram. More...
std::vector< double >	fBins
	Reco Energy bin limits. More...
std::vector< double >	fTrueEBins
	True energy bin limits. More...
std::vector< double >	fDistBins
	Distance bin limits. More...

Detailed Description

Definition at line 55 of file Chi2Sensitivity.h.

Constructor & Destructor Documentation

ana::SBNOsc::Chi2Sensitivity::EventSample::EventSample

(

const fhicl::ParameterSet &

config

)

Definition at line 161 of file Chi2Sensitivity.cxx.

                                                                       {
    // scaling stuff
    fName = config.get<std::string>("name", "");
    fScalePOT = config.get<double>("ScalePOT", 0.);
    fPOT = 0.;

    // setup detector stuff
    std::vector<double> xlim = config.get<std::vector<double> >("DetX");
    std::vector<double> ylim = config.get<std::vector<double> >("DetY");
    std::vector<double> zlim = config.get<std::vector<double> >("DetZ");
    fXlim[0] = xlim[0];
    fXlim[1] = xlim[1];
    fYlim[0] = ylim[0];
    fYlim[1] = ylim[1];
    fZlim[0] = zlim[0];
    fZlim[1] = zlim[1];

    // oscillation type
    fOscType = config.get<int>("OscType", 0);

    // bins of things
    //
    // Reco energy
    std::vector<double> bins = config.get<std::vector<double> >("binlims");
    for (auto const& bin: bins) {
        fBins.push_back(bin);
    }

    // True energy
    std::vector<double> truee = config.get<std::vector<double> >("TrueELims");
    double minTrueE = truee.at(0);
    double maxTrueE = truee.at(1);
    int numTrueEBinLimits = config.get<int>("NumTrueEBins", 1) + 1;
    double trueE_binwidth = (maxTrueE - minTrueE) / numTrueEBinLimits;
    for (int i = 0; i < numTrueEBinLimits; i ++) {
        fTrueEBins.push_back(minTrueE + i * trueE_binwidth);
    }

    // distance
    //
    // Take distance along z-axis as limits
    int numBinsPerMeter = config.get<int>("NumDistanceBinsPerMeter", 1);
    double numBinsPerCM = numBinsPerMeter * 0.01; // unit conversion

    // get beam center and baseline in cm
    fBaseline = config.get<double>("Baseline"); 
    fBeamCenterX = config.get<double>("BeamCenterX", 0.);
    fBeamCenterY = config.get<double>("BeamCenterY", 0.);
    fBeamFrontZ = config.get<double>("BeamFrontZ", 0.);

    // get full range of distances across detector
    double detector_distance = sqrt((fZlim[1] - fZlim[0]) * (fZlim[1] - fZlim[0]) \
        + (fXlim[1] - fXlim[0]) * (fXlim[1] - fXlim[0]) \
        + (fYlim[1] - fYlim[0]) * (fYlim[1] - fYlim[0]));
    // get full range of beam
    // taken approximately from eyeballing MCFlux information
    double beam_width = sqrt( 125. * 125. /* x */ + 125. * 125. /* y */ + 5500. * 5500. /* z */);
    unsigned n_bins = (unsigned) ((detector_distance + beam_width) * numBinsPerCM) + 1;
    unsigned n_limits = n_bins + 1;
    double dist_binwidth = 1./numBinsPerCM;

    // minimum distance from decayed vertex to detector
    // taken approximately from MCFlux info
    double min_baseline = fBaseline - 5500.; 
 
    for (unsigned i = 0; i < n_limits; i++) {
        fDistBins.push_back((min_baseline + (i * dist_binwidth)) / 100000. /* cm -> km */);
    }

    // scaling in reco energy bins
    fEnergyBinScale = config.get<std::vector<double>>("energy_bin_scale", {});
    if (fEnergyBinScale.size() == 0) {
      fEnergyBinScale = std::vector<double>(fBins.size() - 1, 1.0);
    }
    else assert(fEnergyBinScale.size() == fBins.size() - 1);

    // setup histograms
    fBkgCounts = new TH1D((fName + " Background").c_str(), (fName + " Background").c_str(), fBins.size() - 1, &fBins[0]);
    fSignalCounts = new TH3D((fName + " Signal").c_str(), (fName + " Signal").c_str(), 
        fBins.size() - 1, &fBins[0], fTrueEBins.size() - 1, &fTrueEBins[0], fDistBins.size() - 1, &fDistBins[0]);

}

Member Function Documentation

std::vector< double > ana::SBNOsc::Chi2Sensitivity::EventSample::Background

(

)

const

Definition at line 101 of file Chi2Sensitivity.cxx.

                                                               {
    std::vector<double> ret(fBkgCounts->GetNbinsX(), 0);
    for (unsigned bin = 0; bin < fBkgCounts->GetNbinsX(); bin++) {
        ret[bin] = fBkgCounts->GetBinContent(bin+1);
    }
    return ret;
}

std::vector< double > ana::SBNOsc::Chi2Sensitivity::EventSample::Oscillate	(	double	sinth,
		double	dm2
	)		const

Definition at line 57 of file Chi2Sensitivity.cxx.

                                                                                      {
    ROOT::Math::IntegratorOneDim integrator;
    // return histogram binned by reco energy
    std::vector<double> ret(fBkgCounts->GetNbinsX(), 0);
    // iterate over signal bins
    for (unsigned Ebin = 0; Ebin < fSignalCounts->GetNbinsY(); Ebin++) {
        double energy = (fTrueEBins[Ebin] + fTrueEBins[Ebin+1]) / 2.;
        for (unsigned Dbin = 0; Dbin < fSignalCounts->GetNbinsZ(); Dbin++) {
            double distance = (fDistBins[Dbin] + fDistBins[Dbin+1]) / 2.;
            for (unsigned bin = 0; bin < fSignalCounts->GetNbinsX(); bin++) {
                double counts = fSignalCounts->GetBinContent(bin+1, Ebin+1, Dbin+1); 

                if (fOscType == 1) {
                    counts *= numu_to_nue(distance/energy, sinth, dm2); 
                }
                else if (fOscType == 2) {
                    counts *= NumuOscillate(integrator, fDistBins[Dbin], fDistBins[Dbin+1], fTrueEBins[Ebin], fTrueEBins[Ebin+1], dm2, sinth);
                }

                ret[bin] += counts;  
            }
        }
    }

    return ret;
}

std::vector< double > ana::SBNOsc::Chi2Sensitivity::EventSample::Signal

(

)

const

Definition at line 84 of file Chi2Sensitivity.cxx.

                                                           {
    // return histogram binned by reco energy
    std::vector<double> ret(fBkgCounts->GetNbinsX(), 0);
    // iterate over signal bins
    for (unsigned Ebin = 0; Ebin < fSignalCounts->GetNbinsY(); Ebin++) {
        for (unsigned Dbin = 0; Dbin < fSignalCounts->GetNbinsZ(); Dbin++) {
            for (unsigned bin = 0; bin < fSignalCounts->GetNbinsX(); bin++) {
                double counts = fSignalCounts->GetBinContent(bin+1, Ebin+1, Dbin+1); 
                ret[bin] += counts;  
            }
        }
    }

    return ret;
}

Member Data Documentation

double ana::SBNOsc::Chi2Sensitivity::EventSample::fBaseline

baseline (along z) in cm

Definition at line 68 of file Chi2Sensitivity.h.

double ana::SBNOsc::Chi2Sensitivity::EventSample::fBeamCenterX

Center of beam in detector coordinates in x-dir in cm.

Definition at line 69 of file Chi2Sensitivity.h.

double ana::SBNOsc::Chi2Sensitivity::EventSample::fBeamCenterY

Center of beam in detector coordinates in y-dir in cm.

Definition at line 70 of file Chi2Sensitivity.h.

double ana::SBNOsc::Chi2Sensitivity::EventSample::fBeamFrontZ

Front face of detector along z-dir in cm.

Definition at line 71 of file Chi2Sensitivity.h.

std::vector<double> ana::SBNOsc::Chi2Sensitivity::EventSample::fBins

Reco Energy bin limits.

Definition at line 85 of file Chi2Sensitivity.h.

TH1D* ana::SBNOsc::Chi2Sensitivity::EventSample::fBkgCounts

Background count Histogram.

Definition at line 81 of file Chi2Sensitivity.h.

std::vector<double> ana::SBNOsc::Chi2Sensitivity::EventSample::fDistBins

Distance bin limits.

Definition at line 87 of file Chi2Sensitivity.h.

std::vector<double> ana::SBNOsc::Chi2Sensitivity::EventSample::fEnergyBinScale

Definition at line 78 of file Chi2Sensitivity.h.

std::string ana::SBNOsc::Chi2Sensitivity::EventSample::fName

Definition at line 67 of file Chi2Sensitivity.h.

int ana::SBNOsc::Chi2Sensitivity::EventSample::fOscType

Oscilaltion type: 0 == None, 1 == numu -> nue, 2 == numu -> numu.

Definition at line 77 of file Chi2Sensitivity.h.

double ana::SBNOsc::Chi2Sensitivity::EventSample::fPOT

Definition at line 76 of file Chi2Sensitivity.h.

double ana::SBNOsc::Chi2Sensitivity::EventSample::fScalePOT

Factor for POT (etc.) scaling.

Definition at line 75 of file Chi2Sensitivity.h.

TH3D* ana::SBNOsc::Chi2Sensitivity::EventSample::fSignalCounts

Signal Count Histogram.

Definition at line 82 of file Chi2Sensitivity.h.

std::vector<double> ana::SBNOsc::Chi2Sensitivity::EventSample::fTrueEBins

True energy bin limits.

Definition at line 86 of file Chi2Sensitivity.h.

std::array<double, 2> ana::SBNOsc::Chi2Sensitivity::EventSample::fXlim

Detector size in cm.

Definition at line 72 of file Chi2Sensitivity.h.

std::array<double, 2> ana::SBNOsc::Chi2Sensitivity::EventSample::fYlim

Detector size in cm.

Definition at line 73 of file Chi2Sensitivity.h.

std::array<double, 2> ana::SBNOsc::Chi2Sensitivity::EventSample::fZlim

Detector size in cm.

Definition at line 74 of file Chi2Sensitivity.h.

---

The documentation for this class was generated from the following files:

• Chi2Sensitivity.h
• Chi2Sensitivity.cxx