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

NuMatrixSpectrum * fndNuData
NuMatrixSpectrum * fndBarData
NuMatrixSpectrum * ffdNuData
NuMatrixSpectrum * ffdBarData
NuMMHelperCPT * fHelper
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

NuMatrixSpectrum* NuMMRunCPT::ffdBarData [private]

Definition at line 43 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

NuMatrixSpectrum* NuMMRunCPT::ffdNuData [private]

Definition at line 42 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

NuMMHelperCPT* NuMMRunCPT::fHelper [private]

Reimplemented from NuMMRun.

Definition at line 44 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

NuMatrixSpectrum* NuMMRunCPT::fndBarData [private]

Definition at line 41 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

NuMatrixSpectrum* NuMMRunCPT::fndNuData [private]

Definition at line 40 of file NuMMRunCPT.h.

Referenced by MakeFDPred(), and WriteFDPredHistos().

Bool_t NuMMRunCPT::fPredictAntiNeutrinos [private]

Controls antineutrino component of prediction.

Definition at line 49 of file NuMMRunCPT.h.

Referenced by PredictNuBars().

Bool_t NuMMRunCPT::fPredictNeutrinos [private]

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 Fri May 24 22:04:16 2013 for loon by  doxygen 1.4.7