NuMMRunCPT Class Reference

#include <NuMMRunCPT.h>

Inheritance diagram for NuMMRunCPT:
NuMMRun

List of all members.

Public Member Functions

 NuMMRunCPT (NuMMHelperCPT *helper, NuMatrixSpectrum *ndNuData, NuMatrixSpectrum *ndBarData, NuMatrixSpectrum *fdNuData, NuMatrixSpectrum *fdBarData)
virtual Double_t ComparePredWithData (const NuMMParameters &pars)
virtual std::pair
< NuMatrixSpectrum,
NuMatrixSpectrum
MakeFDPred (const NuMMParameters &pars)
virtual std::vector< TH1D > WriteFDPredHistos (const NuMMParameters &pars) const
void PredictNus (Bool_t predict)
 Will the neutrino prediction be calculated?
void PredictNuBars (Bool_t predict)
 Will the antineutrino prediction be calculated?

Private Attributes

NuMatrixSpectrumfndNuData
NuMatrixSpectrumfndBarData
NuMatrixSpectrumffdNuData
NuMatrixSpectrumffdBarData
NuMMHelperCPTfHelper
Bool_t fPredictNeutrinos
 Controls neutrino component of prediction.
Bool_t fPredictAntiNeutrinos
 Controls antineutrino component of prediction.

Detailed Description

Definition at line 15 of file NuMMRunCPT.h.


Constructor & Destructor Documentation

NuMMRunCPT::NuMMRunCPT ( NuMMHelperCPT helper,
NuMatrixSpectrum ndNuData,
NuMatrixSpectrum ndBarData,
NuMatrixSpectrum fdNuData,
NuMatrixSpectrum fdBarData 
)

Definition at line 17 of file NuMMRunCPT.cxx.

00022   : NuMMRun(),
00023     fPredictNeutrinos(true),
00024     fPredictAntiNeutrinos(true)
00025 {
00026   fndNuData = ndNuData;
00027   fndBarData = ndBarData;
00028   ffdNuData = fdNuData;
00029   ffdBarData = fdBarData;
00030   fHelper = helper;
00031 }


Member Function Documentation

Double_t NuMMRunCPT::ComparePredWithData ( const NuMMParameters pars  )  [virtual]

Implements NuMMRun.

Definition at line 466 of file NuMMRunCPT.cxx.

00467 {
00468   const pair<NuMatrixSpectrum,NuMatrixSpectrum> predictions =
00469     this->MakeFDPred(pars);
00470   Double_t like = 0.0;
00471 
00472   if (fPredictNeutrinos) {
00473     like += this->StatsLikelihood(predictions.first.Spectrum(),
00474                                         ffdNuData->Spectrum());
00475   }
00476   if (fPredictAntiNeutrinos) {
00477     like += this->StatsLikelihood(predictions.second.Spectrum(),
00478                                 ffdBarData->Spectrum());
00479   }
00480   return like;
00481 }

pair< NuMatrixSpectrum, NuMatrixSpectrum > NuMMRunCPT::MakeFDPred ( const NuMMParameters pars  )  [virtual]

Definition at line 225 of file NuMMRunCPT.cxx.

References NuMatrix1D::Add(), NuMMHelperCPT::BarBeamMatrix(), NuMMHelperCPT::BarXSecGraph(), NuMatrixSpectrum::Divide(), NuMMParameters::Dm2(), NuMMParameters::Dm2Bar(), NuMatrixSpectrum::ExtrapolateNDToFD(), NuMMHelperCPT::FDBarEfficiency(), NuMMHelperCPT::FDBarNCContamination(), NuMMHelperCPT::FDBarPurity(), NuMMHelperCPT::FDBarRecoVsTrue(), NuMMHelperCPT::FDNuEfficiency(), NuMMHelperCPT::FDNuNCContamination(), NuMMHelperCPT::FDNuPurity(), NuMMHelperCPT::FDNuRecoVsTrue(), NuMMHelperCPT::FDNuTauEfficiency(), NuMMHelperCPT::FDNuTauRecoVsTrue(), NuMMHelperCPT::FDTauBarEfficiency(), NuMMHelperCPT::FDTauBarRecoVsTrue(), NuMMHelperCPT::FDWrongSignBarEfficiency(), NuMMHelperCPT::FDWrongSignBarRecoVsTrue(), NuMMHelperCPT::FDWrongSignNuEfficiency(), NuMMHelperCPT::FDWrongSignNuRecoVsTrue(), ffdBarData, NuMMRun::fFDFidMass, ffdNuData, fHelper, fndBarData, NuMMRun::fNDFidMass, fndNuData, NuMMRun::fQuietMode, NuMatrix::GetPOT(), NuMatrixSpectrum::InverseOscillateLinearInterp(), NuMatrixSpectrum::Multiply(), NuMMParameters::NCBackgroundScale(), NuMMHelperCPT::NDBarEfficiency(), NuMMHelperCPT::NDBarPurity(), NuMMHelperCPT::NDBarRecoVsTrue(), NuMMHelperCPT::NDNuEfficiency(), NuMMHelperCPT::NDNuPurity(), NuMMHelperCPT::NDNuRecoVsTrue(), NuMMParameters::Normalisation(), NuMMHelperCPT::NuBeamMatrix(), NuMMHelperCPT::NuXSecGraph(), NuMatrixSpectrum::OscillateLinearInterp(), NuMatrixSpectrum::RecoToTrue(), NuMatrix::ResetPOT(), NuMMParameters::ShwEnScale(), NuMMParameters::Sn2(), NuMMParameters::Sn2Bar(), NuMatrixSpectrum::TrueToReco(), NuMMHelperCPT::XSecGraphNuTaus(), and NuMMHelperCPT::XSecGraphTauBars().

00226 {
00227   if (!fQuietMode) {
00228     cout << "sn2: " << pars.Sn2() << "; dm2: " << pars.Dm2() 
00229     << "; sn2bar: " << pars.Sn2Bar() << "; dm2bar: " << pars.Dm2Bar()
00230     << endl
00231     << "; norm: " << pars.Normalisation() 
00232     << "; NCBack: " << pars.NCBackgroundScale()
00233     << "; ShwEn: " << pars.ShwEnScale() << endl;
00234   }
00235   
00236   //Get the neutrinos to the FD
00237   NuMatrixSpectrum nuPrediction(*fndNuData);
00238   nuPrediction.Multiply(fHelper->NDNuPurity());
00239   nuPrediction.RecoToTrue(fHelper->NDNuRecoVsTrue());
00240   nuPrediction.Divide(fHelper->NDNuEfficiency());
00241   nuPrediction.Divide(fHelper->NuXSecGraph());
00242   nuPrediction.Divide(fndNuData->GetPOT());
00243   nuPrediction.Divide(fNDFidMass);
00244   nuPrediction.ExtrapolateNDToFD(fHelper->NuBeamMatrix());
00245   nuPrediction.Multiply(fFDFidMass);
00246   nuPrediction.Multiply(ffdNuData->GetPOT());
00247   nuPrediction.ResetPOT(ffdNuData->GetPOT());
00248 
00249   //Get the antineutrinos to the FD
00250   NuMatrixSpectrum barPrediction(*fndBarData);
00251   barPrediction.Multiply(fHelper->NDBarPurity());
00252   barPrediction.RecoToTrue(fHelper->NDBarRecoVsTrue());
00253   barPrediction.Divide(fHelper->NDBarEfficiency());
00254   barPrediction.Divide(fHelper->BarXSecGraph());
00255   barPrediction.Divide(fndBarData->GetPOT());
00256   barPrediction.Divide(fNDFidMass);
00257   barPrediction.ExtrapolateNDToFD(fHelper->BarBeamMatrix());
00258   barPrediction.Multiply(fFDFidMass);
00259   barPrediction.Multiply(ffdBarData->GetPOT());
00260   barPrediction.ResetPOT(ffdBarData->GetPOT());
00261 
00262   //Get the taus (ignoring any wrong-sign taubars)
00263   NuMatrixSpectrum potentialNuTaus(nuPrediction);
00264   potentialNuTaus.Multiply(fHelper->XSecGraphNuTaus());
00265   potentialNuTaus.Multiply(fHelper->FDNuTauEfficiency());
00266   potentialNuTaus.InverseOscillateLinearInterp(pars.Dm2(),pars.Sn2());
00267   potentialNuTaus.TrueToReco(fHelper->FDNuTauRecoVsTrue());
00268   //Get the taubars (ignoring any wrong-sign taus)
00269   NuMatrixSpectrum potentialTauBars(barPrediction);
00270   potentialTauBars.Multiply(fHelper->XSecGraphTauBars());
00271   potentialTauBars.Multiply(fHelper->FDTauBarEfficiency());
00272   potentialTauBars.InverseOscillateLinearInterp(pars.Dm2Bar(),pars.Sn2Bar());
00273   potentialTauBars.TrueToReco(fHelper->FDTauBarRecoVsTrue());
00274 
00275   //Cross sections
00276   nuPrediction.Multiply(fHelper->NuXSecGraph());
00277   barPrediction.Multiply(fHelper->BarXSecGraph());
00278 
00279   //Get the neutrino NC background
00280   NuMatrixSpectrum nuNCBackground(nuPrediction);
00281   nuNCBackground.Multiply(fHelper->FDNuEfficiency());
00282   nuNCBackground.TrueToReco(fHelper->FDNuRecoVsTrue());
00283   nuNCBackground.Divide(fHelper->FDNuPurity());
00284   nuNCBackground.Multiply(fHelper->FDNuNCContamination());
00285 
00286   //Get the antineutrino NC background
00287   NuMatrixSpectrum barNCBackground(barPrediction);
00288   barNCBackground.Multiply(fHelper->FDBarEfficiency());
00289   barNCBackground.TrueToReco(fHelper->FDBarRecoVsTrue());
00290   barNCBackground.Divide(fHelper->FDBarPurity());
00291   barNCBackground.Multiply(fHelper->FDBarNCContamination());
00292   
00293   //Get the wrong-sign neutrino background for the antineutrino
00294   //prediction
00295   NuMatrixSpectrum wrongSignNuMus(nuPrediction);
00296   wrongSignNuMus.Multiply(fHelper->FDWrongSignNuEfficiency());
00297   wrongSignNuMus.OscillateLinearInterp(pars.Dm2(),pars.Sn2());
00298   wrongSignNuMus.TrueToReco(fHelper->FDWrongSignNuRecoVsTrue());
00299 
00300   //Get the wrong-sign antineutrino background for the neutrino
00301   //prediction
00302   NuMatrixSpectrum wrongSignNuBars(barPrediction);
00303   wrongSignNuBars.Multiply(fHelper->FDWrongSignBarEfficiency());
00304   wrongSignNuBars.OscillateLinearInterp(pars.Dm2Bar(),pars.Sn2Bar());
00305   wrongSignNuBars.TrueToReco(fHelper->FDWrongSignBarRecoVsTrue());
00306 
00307   //Oscillations
00308   nuPrediction.OscillateLinearInterp(pars.Dm2(),pars.Sn2());
00309   barPrediction.OscillateLinearInterp(pars.Dm2Bar(),pars.Sn2Bar());
00310   //Efficiencies
00311   nuPrediction.Multiply(fHelper->FDNuEfficiency());
00312   barPrediction.Multiply(fHelper->FDBarEfficiency());
00313 
00314   //True to reco
00315   nuPrediction.TrueToReco(fHelper->FDNuRecoVsTrue());
00316   barPrediction.TrueToReco(fHelper->FDBarRecoVsTrue());
00317 
00318   //Add in backgrounds
00319   nuPrediction.Add(wrongSignNuBars);
00320   barPrediction.Add(wrongSignNuMus);
00321   nuPrediction.Add(nuNCBackground);
00322   nuPrediction.Add(potentialNuTaus);
00323   barPrediction.Add(barNCBackground);
00324   barPrediction.Add(potentialTauBars);
00325 
00326   pair<NuMatrixSpectrum,NuMatrixSpectrum>
00327     predictions(nuPrediction,barPrediction);
00328   return predictions;
00329 }

void NuMMRunCPT::PredictNuBars ( Bool_t  predict  )  [inline]

Will the antineutrino prediction be calculated?

Definition at line 35 of file NuMMRunCPT.h.

References fPredictAntiNeutrinos.

00035                                      {
00036       fPredictAntiNeutrinos = predict;
00037   }

void NuMMRunCPT::PredictNus ( Bool_t  predict  )  [inline]

Will the neutrino prediction be calculated?

Definition at line 31 of file NuMMRunCPT.h.

References fPredictNeutrinos.

00031                                   {
00032       fPredictNeutrinos = predict;
00033   }

std::vector< TH1D > NuMMRunCPT::WriteFDPredHistos ( const NuMMParameters pars  )  const [virtual]

Implements NuMMRun.

Definition at line 35 of file NuMMRunCPT.cxx.

References NuMatrix1D::Add(), NuMMHelperCPT::BarBeamMatrix(), NuMMHelperCPT::BarXSecGraph(), NuMatrixSpectrum::Divide(), NuMMParameters::Dm2(), NuMMParameters::Dm2Bar(), NuMatrixSpectrum::ExtrapolateNDToFD(), NuMMHelperCPT::FDBarEfficiency(), NuMMHelperCPT::FDBarNCContamination(), NuMMHelperCPT::FDBarPurity(), NuMMHelperCPT::FDBarRecoVsTrue(), NuMMHelperCPT::FDNuEfficiency(), NuMMHelperCPT::FDNuNCContamination(), NuMMHelperCPT::FDNuPurity(), NuMMHelperCPT::FDNuRecoVsTrue(), NuMMHelperCPT::FDNuTauEfficiency(), NuMMHelperCPT::FDNuTauRecoVsTrue(), NuMMHelperCPT::FDTauBarEfficiency(), NuMMHelperCPT::FDTauBarRecoVsTrue(), NuMMHelperCPT::FDWrongSignBarEfficiency(), NuMMHelperCPT::FDWrongSignBarRecoVsTrue(), NuMMHelperCPT::FDWrongSignNuEfficiency(), NuMMHelperCPT::FDWrongSignNuRecoVsTrue(), ffdBarData, NuMMRun::fFDFidMass, ffdNuData, fHelper, fndBarData, NuMMRun::fNDFidMass, fndNuData, NuMatrix::GetPOT(), NuMatrixSpectrum::InverseOscillateLinearInterp(), NuMatrixSpectrum::Multiply(), NuMMParameters::NCBackgroundScale(), NuMMHelperCPT::NDBarEfficiency(), NuMMHelperCPT::NDBarPurity(), NuMMHelperCPT::NDBarRecoVsTrue(), NuMMHelperCPT::NDNuEfficiency(), NuMMHelperCPT::NDNuPurity(), NuMMHelperCPT::NDNuRecoVsTrue(), NuMMParameters::Normalisation(), NuMMHelperCPT::NuBeamMatrix(), NuMMHelperCPT::NuXSecGraph(), NuMatrixSpectrum::OscillateLinearInterp(), NuMatrixSpectrum::RecoToTrue(), NuMatrix::ResetPOT(), NuMMParameters::ShwEnScale(), NuMMParameters::Sn2(), NuMMParameters::Sn2Bar(), NuMatrixSpectrum::Spectrum(), NuMatrixSpectrum::TrueToReco(), NuMMHelperCPT::XSecGraphNuTaus(), and NuMMHelperCPT::XSecGraphTauBars().

00036 {
00037   //Set up a vector to push the histograms into.
00038   vector<TH1D> vHistos;
00039 
00040   //Put the ND data in the vector:
00041   TH1D hNDNuData(*(fndNuData->Spectrum()));
00042   hNDNuData.SetName("ndDataNQ");
00043   hNDNuData.SetTitle("ndDataNQ");
00044   vHistos.push_back(*(new TH1D(hNDNuData)));
00045   TH1D hNDBarData(*(fndBarData->Spectrum()));
00046   hNDBarData.SetName("ndDataPQ");
00047   hNDBarData.SetTitle("ndDataPQ");
00048   vHistos.push_back(*(new TH1D(hNDBarData)));
00049 
00050 cout << "sn2: " << pars.Sn2() << "; dm2: " << pars.Dm2() 
00051      << "; sn2bar: " << pars.Sn2Bar() << "; dm2bar: " << pars.Dm2Bar()
00052      << endl
00053      << "; norm: " << pars.Normalisation() 
00054      << "; NCBack: " << pars.NCBackgroundScale()
00055      << "; ShwEn: " << pars.ShwEnScale() << endl;
00056 
00057   //Get the neutrinos to the FD
00058   NuMatrixSpectrum nuPrediction(*fndNuData);
00059   nuPrediction.Multiply(fHelper->NDNuPurity());
00060   nuPrediction.RecoToTrue(fHelper->NDNuRecoVsTrue());
00061   nuPrediction.Divide(fHelper->NDNuEfficiency());
00062   nuPrediction.Divide(fHelper->NuXSecGraph());
00063   nuPrediction.Divide(fndNuData->GetPOT());
00064   nuPrediction.Divide(fNDFidMass);
00065   nuPrediction.ExtrapolateNDToFD(fHelper->NuBeamMatrix());
00066   nuPrediction.Multiply(fFDFidMass);
00067   nuPrediction.Multiply(ffdNuData->GetPOT());
00068   nuPrediction.ResetPOT(ffdNuData->GetPOT());
00069   
00070   //Get the antineutrinos to the FD
00071   NuMatrixSpectrum barPrediction(*fndBarData);
00072   barPrediction.Multiply(fHelper->NDBarPurity());
00073   barPrediction.RecoToTrue(fHelper->NDBarRecoVsTrue());
00074   barPrediction.Divide(fHelper->NDBarEfficiency());
00075   barPrediction.Divide(fHelper->BarXSecGraph());
00076   barPrediction.Divide(fndBarData->GetPOT());
00077   barPrediction.Divide(fNDFidMass);
00078   barPrediction.ExtrapolateNDToFD(fHelper->BarBeamMatrix());
00079   barPrediction.Multiply(fFDFidMass);
00080   barPrediction.Multiply(ffdBarData->GetPOT());
00081   barPrediction.ResetPOT(ffdBarData->GetPOT());
00082   
00083   //Get the taus (ignoring any wrong-sign taubars)
00084   NuMatrixSpectrum potentialNuTaus(nuPrediction);
00085   potentialNuTaus.Multiply(fHelper->XSecGraphNuTaus());
00086   potentialNuTaus.Multiply(fHelper->FDNuTauEfficiency());
00087   potentialNuTaus.InverseOscillateLinearInterp(pars.Dm2(),pars.Sn2());
00088   //Get the taubars (ignoring any wrong-sign taus)
00089   NuMatrixSpectrum potentialTauBars(barPrediction);
00090   potentialTauBars.Multiply(fHelper->XSecGraphTauBars());
00091   potentialTauBars.Multiply(fHelper->FDTauBarEfficiency());
00092   potentialTauBars.InverseOscillateLinearInterp(pars.Dm2Bar(),pars.Sn2Bar());
00093   
00094   // Write out the true tau spectra
00095   TH1D hFDTausNuMu(*(potentialNuTaus.Spectrum()));
00096   hFDTausNuMu.SetName("fdTausNuMu");
00097   vHistos.push_back(hFDTausNuMu);
00098   TH1D hFDTausNuBar(*(potentialTauBars.Spectrum()));
00099   hFDTausNuBar.SetName("fdTausNuBar");
00100   vHistos.push_back(hFDTausNuBar);
00101 
00102   // Do the tau true to reco speerately so that we can write out the spectra
00103   potentialNuTaus.TrueToReco(fHelper->FDNuTauRecoVsTrue());
00104   potentialTauBars.TrueToReco(fHelper->FDTauBarRecoVsTrue());
00105 
00106   //Cross sections
00107   nuPrediction.Multiply(fHelper->NuXSecGraph());
00108   barPrediction.Multiply(fHelper->BarXSecGraph());
00109 
00110   //Get the neutrino NC background
00111   NuMatrixSpectrum nuNCBackground(nuPrediction);
00112   nuNCBackground.Multiply(fHelper->FDNuEfficiency());
00113   nuNCBackground.TrueToReco(fHelper->FDNuRecoVsTrue());
00114   nuNCBackground.Divide(fHelper->FDNuPurity());
00115   nuNCBackground.Multiply(fHelper->FDNuNCContamination());
00116 
00117   //Get the antineutrino NC background
00118   NuMatrixSpectrum barNCBackground(barPrediction);
00119   barNCBackground.Multiply(fHelper->FDBarEfficiency());
00120   barNCBackground.TrueToReco(fHelper->FDBarRecoVsTrue());
00121   barNCBackground.Divide(fHelper->FDBarPurity());
00122   barNCBackground.Multiply(fHelper->FDBarNCContamination());
00123   
00124   //Get the wrong-sign neutrino background for the antineutrino
00125   //prediction
00126   NuMatrixSpectrum wrongSignNuMus(nuPrediction);
00127   wrongSignNuMus.Multiply(fHelper->FDWrongSignNuEfficiency());
00128   wrongSignNuMus.OscillateLinearInterp(pars.Dm2(),pars.Sn2());
00129 
00130   //Get the wrong-sign antineutrino background for the neutrino
00131   //prediction
00132   NuMatrixSpectrum wrongSignNuBars(barPrediction);
00133   wrongSignNuBars.Multiply(fHelper->FDWrongSignBarEfficiency());
00134   wrongSignNuBars.OscillateLinearInterp(pars.Dm2Bar(),pars.Sn2Bar());
00135 
00136 
00137   TH1D hFDWrongSignNQ(*(wrongSignNuBars.Spectrum()));
00138   hFDWrongSignNQ.SetName("fdWrongSignNQ");
00139   vHistos.push_back(hFDWrongSignNQ);
00140 
00141   TH1D hFDWrongSignPQ(*(wrongSignNuMus.Spectrum()));
00142   hFDWrongSignPQ.SetName("fdWrongSignPQ");
00143   vHistos.push_back(hFDWrongSignPQ);
00144 
00145   // Do the truetoreco seperately so that we can write out the histos first
00146   wrongSignNuBars.TrueToReco(fHelper->FDWrongSignBarRecoVsTrue());
00147   wrongSignNuMus.TrueToReco(fHelper->FDWrongSignNuRecoVsTrue());
00148 
00149   //Oscillations
00150   nuPrediction.OscillateLinearInterp(pars.Dm2(),pars.Sn2());
00151   barPrediction.OscillateLinearInterp(pars.Dm2Bar(),pars.Sn2Bar());
00152   //Efficiencies
00153   nuPrediction.Multiply(fHelper->FDNuEfficiency());
00154   barPrediction.Multiply(fHelper->FDBarEfficiency());
00155 
00156   // Write out the true spectra for the prediction components
00157   TH1D hFDNoBackTrueNQ(*(nuPrediction.Spectrum()));
00158   hFDNoBackTrueNQ.SetName("fdBasePredictionNuMus");
00159   vHistos.push_back(hFDNoBackTrueNQ);
00160 
00161   TH1D hFDNoBackTruePQ(*(barPrediction.Spectrum()));
00162   hFDNoBackTruePQ.SetName("fdBasePredictionNuBars");
00163   vHistos.push_back(hFDNoBackTruePQ);
00164   
00165 //True to reco
00166   nuPrediction.TrueToReco(fHelper->FDNuRecoVsTrue());
00167   barPrediction.TrueToReco(fHelper->FDBarRecoVsTrue());
00168 
00169   // Dump out the prediction spectrum before adding in the backgrounds
00170   TH1D hFDNoBackNQ(*(nuPrediction.Spectrum()));
00171   hFDNoBackNQ.SetName("fdBasePredictionNQ");
00172   vHistos.push_back(hFDNoBackNQ);
00173 
00174   TH1D hFDNoBackPQ(*(barPrediction.Spectrum()));
00175   hFDNoBackPQ.SetName("fdBasePredictionPQ");
00176   vHistos.push_back(hFDNoBackPQ);
00177 
00178   //Add in backgrounds
00179   nuPrediction.Add(wrongSignNuBars);
00180   barPrediction.Add(wrongSignNuMus);
00181   nuPrediction.Add(nuNCBackground);
00182   nuPrediction.Add(potentialNuTaus);
00183   barPrediction.Add(barNCBackground);
00184   barPrediction.Add(potentialTauBars);
00185 
00186   //Put the componants of the predictions in the vector
00187 
00188   TH1D hFDPredictionNQ(*(nuPrediction.Spectrum()));
00189   hFDPredictionNQ.SetName("fdPredictionNQ");
00190   vHistos.push_back(hFDPredictionNQ);
00191 
00192   TH1D hFDPredictionPQ(*(barPrediction.Spectrum()));
00193   hFDPredictionPQ.SetName("fdPredictionPQ");
00194   vHistos.push_back(hFDPredictionPQ);
00195 
00196   TH1D hFDNCNQ(*(nuNCBackground.Spectrum()));
00197   hFDNCNQ.SetName("fdNCNQ");
00198   vHistos.push_back(hFDNCNQ);
00199 
00200   TH1D hFDNCPQ(*(barNCBackground.Spectrum()));
00201   hFDNCPQ.SetName("fdNCPQ");
00202   vHistos.push_back(hFDNCPQ);
00203 
00204   TH1D hFDTausNQ(*(potentialNuTaus.Spectrum()));
00205   hFDTausNQ.SetName("fdTausNQ");
00206   vHistos.push_back(hFDTausNQ);
00207 
00208   TH1D hFDTausPQ(*(potentialTauBars.Spectrum()));
00209   hFDTausPQ.SetName("fdTausPQ");
00210   vHistos.push_back(hFDTausPQ);
00211 
00212   TH1D hFDDataNQ(*(ffdNuData->Spectrum()));
00213   hFDDataNQ.SetName("fdDataNQ");
00214   vHistos.push_back(hFDDataNQ);
00215 
00216   TH1D hFDDataPQ(*(ffdBarData->Spectrum()));
00217   hFDDataPQ.SetName("fdDataPQ");
00218   vHistos.push_back(hFDDataPQ);
00219 
00220   return vHistos;
00221 }


Member Data Documentation

Definition at line 43 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

Definition at line 42 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

Reimplemented from NuMMRun.

Definition at line 44 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

Definition at line 41 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

Definition at line 40 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

Controls antineutrino component of prediction.

Definition at line 49 of file NuMMRunCPT.h.

Referenced by PredictNuBars().

Controls neutrino component of prediction.

Definition at line 47 of file NuMMRunCPT.h.

Referenced by PredictNus().


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

Generated on 13 Sep 2017 for loon by  doxygen 1.6.1