SterileSpec Class Reference

#include <SterileSpec.h>

List of all members.

Public Types

enum  whichPred { CC, NC }
enum  whichDet { ND, FD }
enum  whichBinningScheme {
  LE, LEME, rLEME, fine,
  RHC
}

Public Member Functions

 SterileSpec ()
virtual ~SterileSpec ()
void SetNDMCHistograms (TH2D *NDNCTrueNC, TH2D *NDNCNuMu, TH2D *NDNCBeamNue, TH2D *NDNCAppNue, TH2D *NDNCAppNuTau, TH2D *NDCCTrueNC, TH2D *NDCCNuMu, TH2D *NDCCBeamNue, TH2D *NDCCAppNue, TH2D *NDCCAppNuTau)
void SetFDMCHistograms (TH2D *FDNCTrueNC, TH2D *FDNCNuMu, TH2D *FDNCBeamNue, TH2D *FDNCAppNue, TH2D *FDNCAppNuTau, TH2D *FDCCTrueNC, TH2D *FDCCNuMu, TH2D *FDCCBeamNue, TH2D *FDCCAppNue, TH2D *FDCCAppNuTau)
void SetNDMCHistogramsNoRebin (TH2D *NDNCTrueNC, TH2D *NDNCNuMu, TH2D *NDNCBeamNue, TH2D *NDNCAppNue, TH2D *NDNCAppNuTau, TH2D *NDCCTrueNC, TH2D *NDCCNuMu, TH2D *NDCCBeamNue, TH2D *NDCCAppNue, TH2D *NDCCAppNuTau)
void SetFDMCHistogramsNoRebin (TH2D *FDNCTrueNC, TH2D *FDNCNuMu, TH2D *FDNCBeamNue, TH2D *FDNCAppNue, TH2D *FDNCAppNuTau, TH2D *FDCCTrueNC, TH2D *FDCCNuMu, TH2D *FDCCBeamNue, TH2D *FDCCAppNue, TH2D *FDCCAppNuTau)
void SetMatrices (TMatrixD *CCInvert, TMatrixD *NCInvert, TMatrixD *CC, TMatrixD *NC)
void SetNDBaseline (double baseline)
void SetFDBaseline (double baseline)
TH1D * CreateTotalSpectrum (NuMMParameters pars, whichPred Pred, whichDet Det)
std::vector< Double_t > GenerateBinningScheme (whichPred Pred)
TH2D * Rebin4Chi2 (TH2D *orig, whichPred Pred)
TH1D * Rebin4Chi2 (TH1D *orig, whichPred Pred)
void UseDm241 (bool YesorNo)
void SetBinningScheme (whichBinningScheme bs)

Public Attributes

SterileSpecCompComp
TH2D * NDNC_TrueNC
TH2D * NDNC_NuMu
TH2D * NDNC_BeamNue
TH2D * NDNC_AppNue
TH2D * NDNC_AppNuTau
TH2D * NDCC_TrueNC
TH2D * NDCC_NuMu
TH2D * NDCC_BeamNue
TH2D * NDCC_AppNue
TH2D * NDCC_AppNuTau
TH2D * FDNC_TrueNC
TH2D * FDNC_NuMu
TH2D * FDNC_BeamNue
TH2D * FDNC_AppNue
TH2D * FDNC_AppNuTau
TH2D * FDCC_TrueNC
TH2D * FDCC_NuMu
TH2D * FDCC_BeamNue
TH2D * FDCC_AppNue
TH2D * FDCC_AppNuTau
TH1D * vtruenc
TH1D * vnumu
TH1D * vbeamnue
TH1D * vappnue
TH1D * vappnutau

Private Member Functions

 ClassDef (SterileSpec, 0)

Private Attributes

TMatrixD * CoVarCCinvert
TMatrixD * CoVarNCinvert
TMatrixD * CoVarCC
TMatrixD * CoVarNC
Double_t ND_baseline
Double_t FD_baseline
bool use41
whichBinningScheme binningScheme

Detailed Description

Definition at line 19 of file SterileSpec.h.


Member Enumeration Documentation

Enumerator:
LE 
LEME 
rLEME 
fine 
RHC 

Definition at line 24 of file SterileSpec.h.

00024 {LE, LEME, rLEME, fine, RHC};

Enumerator:
ND 
FD 

Definition at line 23 of file SterileSpec.h.

00023 {ND, FD};

Enumerator:
CC 
NC 

Definition at line 22 of file SterileSpec.h.

00022 {CC, NC};


Constructor & Destructor Documentation

SterileSpec::SterileSpec (  ) 

Definition at line 22 of file SterileSpec.cxx.

00022                         : ND_baseline(0.0), FD_baseline(0.0), binningScheme(LE)
00023 {
00024   // Will create the comppnents of the FD and ND spectra.
00025   Comp = new SterileSpecComp();
00026   use41 = false;
00027 }

SterileSpec::~SterileSpec (  )  [virtual]

Definition at line 29 of file SterileSpec.cxx.

00030 {
00031 }


Member Function Documentation

SterileSpec::ClassDef ( SterileSpec  ,
 
) [private]
TH1D * SterileSpec::CreateTotalSpectrum ( NuMMParameters  pars,
whichPred  Pred,
whichDet  Det 
)

Definition at line 124 of file SterileSpec.cxx.

References CC, Comp, SterileSpecComp::CreateSpectrumComponent(), FD, FD_baseline, FDCC_AppNue, FDCC_AppNuTau, FDCC_BeamNue, FDCC_NuMu, FDCC_TrueNC, FDNC_AppNue, FDNC_AppNuTau, FDNC_BeamNue, FDNC_NuMu, FDNC_TrueNC, SterileSpecComp::kAppNue, SterileSpecComp::kAppNuTau, SterileSpecComp::kBeamNue, SterileSpecComp::kNuMu, SterileSpecComp::kTrueNC, SterileSpecComp::LoadGenericHelpers(), NC, ND, ND_baseline, NDCC_AppNue, NDCC_AppNuTau, NDCC_BeamNue, NDCC_NuMu, NDCC_TrueNC, NDNC_AppNue, NDNC_AppNuTau, NDNC_BeamNue, NDNC_NuMu, NDNC_TrueNC, vappnue, vappnutau, vbeamnue, vnumu, and vtruenc.

Referenced by FarOverNearFit::CalculateNDIntegrals(), FarOverNearFitPQ::CalculateNDIntegrals(), FarOverNearFit::ComparePredWithData(), FarOverNearFitPQ::ComparePredWithData(), FarOverNearFitPQ::CreateFakeHighStatsData(), FarOverNearFit::CreateFakeHighStatsData(), and FarOverNearFit::CreatePrediction().

00125 {
00126   Double_t baseline = 0.0;
00127 
00128   if(Det == FD){
00129     
00130     baseline = FD_baseline;
00131 
00132     if(Pred == CC){
00133       Comp->LoadGenericHelpers(FDCC_TrueNC, FDCC_NuMu, FDCC_BeamNue, FDCC_AppNue, FDCC_AppNuTau);
00134     }
00135     else if(Pred == NC){
00136       Comp->LoadGenericHelpers(FDNC_TrueNC, FDNC_NuMu, FDNC_BeamNue, FDNC_AppNue, FDNC_AppNuTau);
00137     }
00138     else{
00139       std::cout<< "You didnt specify which event type or detector... asserting." << std::endl;
00140       assert(false);
00141     }
00142   }
00143   if(Det == ND){
00144 
00145     baseline = ND_baseline;
00146 
00147     if(Pred == CC){
00148       Comp->LoadGenericHelpers(NDCC_TrueNC, NDCC_NuMu, NDCC_BeamNue, NDCC_AppNue, NDCC_AppNuTau);
00149     }
00150     else if(Pred == NC){
00151       Comp->LoadGenericHelpers(NDNC_TrueNC, NDNC_NuMu, NDNC_BeamNue, NDNC_AppNue, NDNC_AppNuTau);
00152     }
00153     else{ 
00154       std::cout << "You didnt specify which event type or detector... asserting." << std::endl;
00155       assert(false);
00156     }
00157   }
00158 
00159   // Now the right helpers are loaded, create the components
00160   vtruenc   = Comp->CreateSpectrumComponent(pars, SterileSpecComp::kTrueNC,   baseline);
00161   vnumu     = Comp->CreateSpectrumComponent(pars, SterileSpecComp::kNuMu,     baseline);
00162   vbeamnue  = Comp->CreateSpectrumComponent(pars, SterileSpecComp::kBeamNue,  baseline);
00163   vappnue   = Comp->CreateSpectrumComponent(pars, SterileSpecComp::kAppNue,   baseline);
00164   vappnutau = Comp->CreateSpectrumComponent(pars, SterileSpecComp::kAppNuTau, baseline);
00165 
00166   TH1D* hTotal = new TH1D(*vtruenc);
00167   hTotal->Add(vnumu);
00168   hTotal->Add(vbeamnue);
00169   hTotal->Add(vappnue);
00170   hTotal->Add(vappnutau);
00171 
00172   return hTotal;
00173 }

std::vector< Double_t > SterileSpec::GenerateBinningScheme ( whichPred  Pred  ) 

Definition at line 242 of file SterileSpec.cxx.

References binningScheme, CC, MuELoss::e, fine, LE, LEME, NC, RHC, and rLEME.

Referenced by Rebin4Chi2().

00243 {
00244   // Define the Binning
00245   std::vector<Double_t> binscheme;
00246 
00247   // Set min Energy
00248   binscheme.push_back(0.0);
00249 
00250   // Start the Bin Count
00251   Int_t nbins = 1;
00252 
00253   // binning scheme for MINOS
00254   if(binningScheme == LE){
00255     if(Pred == CC){
00256       binscheme.push_back(2.0); nbins++; Double_t binWidth = 0.25;
00257       while(binscheme[nbins-1] < 5-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00258       binWidth = 0.5;
00259       while(binscheme[nbins-1] < 15-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00260       binWidth = 1.00;
00261       while(binscheme[nbins-1] < 20-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00262       binWidth = 2.0;
00263       while(binscheme[nbins-1] < 30-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00264       binWidth = 2.0;
00265       while(binscheme[nbins-1] < 40-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00266     }
00267 
00268     if(Pred == NC){
00269       binscheme.push_back(0.5); nbins++; Double_t binWidth = 0.25;
00270       while(binscheme[nbins-1] < 5-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00271       binWidth = 4;
00272       while(binscheme[nbins-1] < 20-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00273       binWidth = 19.00;
00274       while(binscheme[nbins-1] < 40-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00275     }
00276   }
00277   // binning scheme for overlayed MINOS and MINOS+ 
00278   else if(binningScheme == LEME){
00279     if(Pred == CC){
00280       binscheme.push_back(1.75); nbins++; Double_t binWidth = 0.25;
00281       while(binscheme[nbins-1] < 15-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00282       binWidth = 1.;
00283       while(binscheme[nbins-1] < 30-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00284       binWidth = 2.;
00285       while(binscheme[nbins-1] < 40-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00286     }
00287 
00288     if(Pred == NC){
00289       binscheme.push_back(0.75); nbins++; Double_t binWidth = 0.25;
00290       while(binscheme[nbins-1] < 8-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00291       binWidth = 2.;
00292       while(binscheme[nbins-1] < 20-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00293       binWidth = 10.;
00294       while(binscheme[nbins-1] < 40-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00295     }
00296   }
00297   // reduced binning scheme for overlayed MINOS and MINOS+ 
00298   else if(binningScheme == rLEME){
00299     if(Pred == CC){
00300       binscheme.push_back(2.0); nbins++; Double_t binWidth = 1.0;
00301       while(binscheme[nbins-1] < 15-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00302       binWidth = 3.;
00303       while(binscheme[nbins-1] < 30-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00304       binWidth = 5.;
00305       while(binscheme[nbins-1] < 40-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00306     }
00307     if(Pred == NC){
00308       binscheme.push_back(1.0); nbins++; Double_t binWidth = 1.0;
00309       while(binscheme[nbins-1] < 8-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00310       binWidth = 4.;
00311       while(binscheme[nbins-1] < 20-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00312       binWidth = 10.;
00313       while(binscheme[nbins-1] < 40-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00314     }
00315   }
00316   // fine binning scheme, binning used in helper files generation 
00317   else if(binningScheme == fine){
00318     if(Pred == CC){
00319       binscheme.push_back(0.05); nbins++; Double_t binWidth = 0.05;
00320       while(binscheme[nbins-1] < 5-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00321       binWidth = 0.25;
00322       while(binscheme[nbins-1] < 20-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00323       binWidth = 1.0;
00324       while(binscheme[nbins-1] < 30-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00325       binWidth = 2.0;
00326       while(binscheme[nbins-1] < 40-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00327     }
00328     if(Pred == NC){
00329       binscheme.push_back(0.05); nbins++; Double_t binWidth = 0.05;
00330       while(binscheme[nbins-1] < 5-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00331       binWidth = 0.25;
00332       while(binscheme[nbins-1] < 20-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00333       binWidth = 1.0;
00334       while(binscheme[nbins-1] < 30-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00335       binWidth = 2.0;
00336       while(binscheme[nbins-1] < 40-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00337     }
00338   }
00339   // RHC MINOS binning scheme 
00340   else if(binningScheme == RHC){
00341     if(Pred == CC){
00342       binscheme.push_back(3.0); nbins++; Double_t binWidth = .5;
00343       while(binscheme[nbins-1] < 4.5-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00344       binWidth = 1.5;
00345       while(binscheme[nbins-1] < 7.5-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00346       binWidth = 2.5;
00347       while(binscheme[nbins-1] < 10-1e-3  ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00348       binWidth = 4.0;
00349       while(binscheme[nbins-1] < 14-1e-3  ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00350       binWidth = 6.0;
00351       while(binscheme[nbins-1] < 20-1e-3  ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00352       binWidth = 20.0;
00353       while(binscheme[nbins-1] < 40-1e-3  ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00354     }
00355     if(Pred == NC){
00356       binscheme.push_back(0.75); nbins++; Double_t binWidth = 0.5;
00357       while(binscheme[nbins-1] < 1.25-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00358       binWidth = 0.75;
00359       while(binscheme[nbins-1] < 3.5-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00360       binWidth = 1.00;
00361       while(binscheme[nbins-1] < 4.5-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00362       binWidth = 1.50;
00363       while(binscheme[nbins-1] < 6.0-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00364       binWidth = 3.00;
00365       while(binscheme[nbins-1] < 9.0-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00366       binWidth = 5.00;
00367       while(binscheme[nbins-1] < 14.0-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00368       binWidth = 6.00;
00369       while(binscheme[nbins-1] < 20.0-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00370       binWidth = 20.00;
00371       while(binscheme[nbins-1] < 40.0-1e-3 ){ binscheme.push_back(binscheme[nbins-1] + binWidth); nbins++;}
00372     }
00373   }
00374   else{
00375     cout << "binning scheme error" << endl;
00376     abort();
00377   }
00378 
00379   binscheme.push_back(200);
00380 
00381   return binscheme;
00382 }

TH1D * SterileSpec::Rebin4Chi2 ( TH1D *  orig,
whichPred  Pred 
)

Definition at line 220 of file SterileSpec.cxx.

References GenerateBinningScheme().

00221 {
00222   std::vector<Double_t> edges = GenerateBinningScheme(Pred);
00223 
00224   TString hName = orig->GetName();
00225   hName += "_RecoRebinned";
00226   TString hTitle = orig->GetTitle();
00227 
00228   TH1D* hRebinned = new TH1D(hName, hTitle, edges.size()-1, &edges[0]);
00229   hRebinned->SetDirectory(0);
00230 
00231   for(Int_t iX = 1; iX <= orig->GetNbinsX(); ++iX){
00232 
00233     Double_t binX = orig->GetXaxis()->GetBinCenter(iX);
00234     Double_t content = orig->GetBinContent(iX);
00235 
00236     hRebinned->Fill(binX, content);
00237   }
00238 
00239   return hRebinned;
00240 }

TH2D * SterileSpec::Rebin4Chi2 ( TH2D *  orig,
whichPred  Pred 
)

Definition at line 181 of file SterileSpec.cxx.

References GenerateBinningScheme().

Referenced by SetFDMCHistograms(), and SetNDMCHistograms().

00182 {
00183   // This was designed for rebinning the reco to true histograms - this only rebins the reco
00184   // because its a waste of CPU to loop over so many reco bins only to rebin after.
00185 
00186   std::vector<Double_t> edges = GenerateBinningScheme(Pred);
00187 
00188   std::vector<Double_t> trueBins;
00189   Int_t nbinsTrue = orig->GetXaxis()->GetNbins();
00190   trueBins.push_back( orig->GetXaxis()->GetBinLowEdge(1) );
00191 
00192   for(Int_t ibin = 1; ibin <= nbinsTrue; ++ibin){
00193     trueBins.push_back( orig->GetXaxis()->GetBinUpEdge(ibin) );
00194   }
00195 
00196   TString hName = orig->GetName();
00197   hName += "_RecoRebinned";
00198   TString hTitle = orig->GetTitle();
00199 
00200   TH2D* hRebinned = new TH2D(hName, hTitle, nbinsTrue, &trueBins[0], edges.size()-1, &edges[0]);
00201   hRebinned->SetDirectory(0);
00202 
00203   // Now refill the contents into the new histogram
00204   for (Int_t iX = 1; iX <= orig->GetXaxis()->GetNbins(); ++iX){
00205 
00206     Double_t  binX = orig->GetXaxis()->GetBinCenter(iX);
00207 
00208     for (Int_t iY = 1; iY <= orig->GetYaxis()->GetNbins(); ++iY){
00209       
00210       Double_t binY = orig->GetYaxis()->GetBinCenter(iY);
00211       Double_t content = orig->GetBinContent(iX, iY);
00212 
00213       hRebinned->Fill(binX, binY, content);
00214     }
00215   }
00216 
00217   return hRebinned;
00218 }

void SterileSpec::SetBinningScheme ( whichBinningScheme  bs  )  [inline]

Definition at line 53 of file SterileSpec.h.

References binningScheme.

00053 {binningScheme = bs;};

void SterileSpec::SetFDBaseline ( double  baseline  ) 

Definition at line 119 of file SterileSpec.cxx.

References FD_baseline.

Referenced by FarOverNearFitPQ::Initalise(), and FarOverNearFit::Initalise().

00120 {
00121   FD_baseline = baseline;
00122 }

void SterileSpec::SetFDMCHistograms ( TH2D *  FDNCTrueNC,
TH2D *  FDNCNuMu,
TH2D *  FDNCBeamNue,
TH2D *  FDNCAppNue,
TH2D *  FDNCAppNuTau,
TH2D *  FDCCTrueNC,
TH2D *  FDCCNuMu,
TH2D *  FDCCBeamNue,
TH2D *  FDCCAppNue,
TH2D *  FDCCAppNuTau 
)

Definition at line 83 of file SterileSpec.cxx.

References CC, FDCC_AppNue, FDCC_AppNuTau, FDCC_BeamNue, FDCC_NuMu, FDCC_TrueNC, FDNC_AppNue, FDNC_AppNuTau, FDNC_BeamNue, FDNC_NuMu, FDNC_TrueNC, NC, and Rebin4Chi2().

Referenced by FarOverNearFit::Add(), FarOverNearFitPQ::Initalise(), and FarOverNearFit::Initalise().

00085 {
00086   // FDMC, Nue and Nutau histograms (rebin for the sterile analysis)
00087   FDNC_TrueNC   = Rebin4Chi2(FDNCTrueNC,   NC);
00088   FDNC_NuMu     = Rebin4Chi2(FDNCNuMu,     NC);
00089   FDNC_BeamNue  = Rebin4Chi2(FDNCBeamNue,  NC);
00090   FDNC_AppNue   = Rebin4Chi2(FDNCAppNue,   NC);
00091   FDNC_AppNuTau = Rebin4Chi2(FDNCAppNuTau, NC);
00092 
00093   FDCC_TrueNC   = Rebin4Chi2(FDCCTrueNC,   CC);
00094   FDCC_NuMu     = Rebin4Chi2(FDCCNuMu,     CC);
00095   FDCC_BeamNue  = Rebin4Chi2(FDCCBeamNue,  CC);
00096   FDCC_AppNue   = Rebin4Chi2(FDCCAppNue,   CC);
00097   FDCC_AppNuTau = Rebin4Chi2(FDCCAppNuTau, CC);
00098 
00099 }

void SterileSpec::SetFDMCHistogramsNoRebin ( TH2D *  FDNCTrueNC,
TH2D *  FDNCNuMu,
TH2D *  FDNCBeamNue,
TH2D *  FDNCAppNue,
TH2D *  FDNCAppNuTau,
TH2D *  FDCCTrueNC,
TH2D *  FDCCNuMu,
TH2D *  FDCCBeamNue,
TH2D *  FDCCAppNue,
TH2D *  FDCCAppNuTau 
)

Definition at line 67 of file SterileSpec.cxx.

References FDCC_AppNue, FDCC_AppNuTau, FDCC_BeamNue, FDCC_NuMu, FDCC_TrueNC, FDNC_AppNue, FDNC_AppNuTau, FDNC_BeamNue, FDNC_NuMu, and FDNC_TrueNC.

00069 {
00070   FDNC_TrueNC   = FDNCTrueNC;
00071   FDNC_NuMu     = FDNCNuMu;
00072   FDNC_BeamNue  = FDNCBeamNue;
00073   FDNC_AppNue   = FDNCAppNue;
00074   FDNC_AppNuTau = FDNCAppNuTau;
00075 
00076   FDCC_TrueNC   = FDCCTrueNC;
00077   FDCC_NuMu     = FDCCNuMu;
00078   FDCC_BeamNue  = FDCCBeamNue;
00079   FDCC_AppNue   = FDCCAppNue;
00080   FDCC_AppNuTau = FDCCAppNuTau;
00081 }

void SterileSpec::SetMatrices ( TMatrixD *  CCInvert,
TMatrixD *  NCInvert,
TMatrixD *  CC,
TMatrixD *  NC 
)

Definition at line 101 of file SterileSpec.cxx.

References CoVarCC, CoVarCCinvert, CoVarNC, and CoVarNCinvert.

Referenced by FarOverNearFitPQ::Initalise(), and FarOverNearFit::Initalise().

00102 {
00103   // This class needs to matrices to create fake data
00104   //These are the matrices used in the fit
00105   CoVarCCinvert = CCInvert;
00106   CoVarNCinvert = NCInvert;
00107 
00108   // These are the oringial matrices before inversion
00109   // these are used for the fluc of systematics
00110   CoVarCC = CC;
00111   CoVarNC = NC;
00112 }

void SterileSpec::SetNDBaseline ( double  baseline  ) 

Definition at line 114 of file SterileSpec.cxx.

References ND_baseline.

Referenced by FarOverNearFitPQ::Initalise(), and FarOverNearFit::Initalise().

00115 {
00116   ND_baseline = baseline;
00117 }

void SterileSpec::SetNDMCHistograms ( TH2D *  NDNCTrueNC,
TH2D *  NDNCNuMu,
TH2D *  NDNCBeamNue,
TH2D *  NDNCAppNue,
TH2D *  NDNCAppNuTau,
TH2D *  NDCCTrueNC,
TH2D *  NDCCNuMu,
TH2D *  NDCCBeamNue,
TH2D *  NDCCAppNue,
TH2D *  NDCCAppNuTau 
)

Definition at line 49 of file SterileSpec.cxx.

References CC, NC, NDCC_AppNue, NDCC_AppNuTau, NDCC_BeamNue, NDCC_NuMu, NDCC_TrueNC, NDNC_AppNue, NDNC_AppNuTau, NDNC_BeamNue, NDNC_NuMu, NDNC_TrueNC, and Rebin4Chi2().

Referenced by FarOverNearFit::Add(), FarOverNearFitPQ::Initalise(), and FarOverNearFit::Initalise().

00051 {
00052   //NDMC, Nue and Nutau histograms  (rebin for the sterile analysis) 
00053   NDNC_TrueNC   = Rebin4Chi2(NDNCTrueNC,   NC);
00054   NDNC_NuMu     = Rebin4Chi2(NDNCNuMu,     NC);
00055   NDNC_BeamNue  = Rebin4Chi2(NDNCBeamNue,  NC);
00056   NDNC_AppNue   = Rebin4Chi2(NDNCAppNue,   NC);
00057   NDNC_AppNuTau = Rebin4Chi2(NDNCAppNuTau, NC);
00058 
00059   NDCC_TrueNC   = Rebin4Chi2(NDCCTrueNC,   CC);
00060   NDCC_NuMu     = Rebin4Chi2(NDCCNuMu,     CC);
00061   NDCC_BeamNue  = Rebin4Chi2(NDCCBeamNue,  CC);
00062   NDCC_AppNue   = Rebin4Chi2(NDCCAppNue,   CC);
00063   NDCC_AppNuTau = Rebin4Chi2(NDCCAppNuTau, CC);
00064 
00065 }

void SterileSpec::SetNDMCHistogramsNoRebin ( TH2D *  NDNCTrueNC,
TH2D *  NDNCNuMu,
TH2D *  NDNCBeamNue,
TH2D *  NDNCAppNue,
TH2D *  NDNCAppNuTau,
TH2D *  NDCCTrueNC,
TH2D *  NDCCNuMu,
TH2D *  NDCCBeamNue,
TH2D *  NDCCAppNue,
TH2D *  NDCCAppNuTau 
)

Definition at line 33 of file SterileSpec.cxx.

References NDCC_AppNue, NDCC_AppNuTau, NDCC_BeamNue, NDCC_NuMu, NDCC_TrueNC, NDNC_AppNue, NDNC_AppNuTau, NDNC_BeamNue, NDNC_NuMu, and NDNC_TrueNC.

00035 {
00036   NDNC_TrueNC   = NDNCTrueNC;
00037   NDNC_NuMu     = NDNCNuMu;
00038   NDNC_BeamNue  = NDNCBeamNue;
00039   NDNC_AppNue   = NDNCAppNue;
00040   NDNC_AppNuTau = NDNCAppNuTau;
00041 
00042   NDCC_TrueNC   = NDCCTrueNC;
00043   NDCC_NuMu     = NDCCNuMu;
00044   NDCC_BeamNue  = NDCCBeamNue;
00045   NDCC_AppNue   = NDCCAppNue;
00046   NDCC_AppNuTau = NDCCAppNuTau;
00047 }

void SterileSpec::UseDm241 ( bool  YesorNo  ) 

Definition at line 175 of file SterileSpec.cxx.

References Comp, use41, and SterileSpecComp::UseDm241().

00176 {
00177   use41 = YesorNo;
00178   Comp->UseDm241(use41);
00179 }


Member Data Documentation

Definition at line 78 of file SterileSpec.h.

Referenced by GenerateBinningScheme(), and SetBinningScheme().

Definition at line 30 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), and UseDm241().

TMatrixD* SterileSpec::CoVarCC [private]

Definition at line 71 of file SterileSpec.h.

Referenced by SetMatrices().

TMatrixD* SterileSpec::CoVarCCinvert [private]

Definition at line 69 of file SterileSpec.h.

Referenced by SetMatrices().

TMatrixD* SterileSpec::CoVarNC [private]

Definition at line 72 of file SterileSpec.h.

Referenced by SetMatrices().

TMatrixD* SterileSpec::CoVarNCinvert [private]

Definition at line 70 of file SterileSpec.h.

Referenced by SetMatrices().

Double_t SterileSpec::FD_baseline [private]

Definition at line 75 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), and SetFDBaseline().

Definition at line 61 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Definition at line 61 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Definition at line 61 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Definition at line 61 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Definition at line 61 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Definition at line 60 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Definition at line 60 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Definition at line 60 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Definition at line 60 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Definition at line 60 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetFDMCHistograms(), and SetFDMCHistogramsNoRebin().

Double_t SterileSpec::ND_baseline [private]

Definition at line 75 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), and SetNDBaseline().

Definition at line 57 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

Definition at line 57 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

Definition at line 57 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

Definition at line 57 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

Definition at line 57 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

Definition at line 53 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

Definition at line 53 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

Definition at line 53 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

Definition at line 53 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

Definition at line 53 of file SterileSpec.h.

Referenced by CreateTotalSpectrum(), SetNDMCHistograms(), and SetNDMCHistogramsNoRebin().

bool SterileSpec::use41 [private]

Definition at line 76 of file SterileSpec.h.

Referenced by UseDm241().

Definition at line 64 of file SterileSpec.h.

Referenced by CreateTotalSpectrum().

Definition at line 64 of file SterileSpec.h.

Referenced by CreateTotalSpectrum().

Definition at line 64 of file SterileSpec.h.

Referenced by CreateTotalSpectrum().

Definition at line 64 of file SterileSpec.h.

Referenced by CreateTotalSpectrum().

Definition at line 64 of file SterileSpec.h.

Referenced by CreateTotalSpectrum().


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

Generated on 22 Nov 2017 for loon by  doxygen 1.6.1