#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-24/3FlavorAna/Ana2018/numu/ND_DataMC_Comp_Systs_header.h"

Public Member Functions
	DataMCPair (Selection sel, TangibleAxis tanAxis, SpectrumLoader &loaderData, SpectrumLoader &loaderMC, std::vector< const ISyst * > systs, const Var BinWeight, const Cut &bkg=kNoCut)

std::string	ShortName () const

const char *	CName () const

void	DrawMCSyst (const int iSyst=-1) const

void	OverlayDataMCSyst (const int iSyst=-1) const

void	DrawMCNormSyst (const int iSyst=-1) const

void	OverlayDataMCSystNorm (const int iSyst=-1) const

void	OverlayDataMCSystExtraData (const DataMCPair &extraData, const std::string dataLegendTitle, const std::string extraLegendTitle, const int iSyst=-1, const int dataColor=kBlack, const int extraColor=kGray+1, const int dataMarker=kFullCircle, const int extraMarker=kFullCircle) const

void	OverlayDataMCSystNormExtraData (const DataMCPair &extraData, const std::string dataLegendTitle, const std::string extraLegendTitle, const int iSyst=-1, const int dataColor=kBlack, const int extraColor=kGray+1, const int dataMarker=kFullCircle, const int extraMarker=kFullCircle) const

TLegend *	DrawLegendArea (const int dataColor=kBlack, std::string dataLegendTitle="Data", const int dataMarker=kFullCircle) const

TLegend *	DrawLegendPOT (const int dataColor=kBlack, std::string dataLegendTitle="Data", const int dataMarker=kFullCircle) const

void	AddExposureArea () const

void	AddExposurePOT () const

void	AestheticsArea () const

void	AestheticsPOT () const

void	DrawData (const int color=kBlack, const int marker=kFullCircle) const

void	DrawMCBkg () const

float	Purity () const

	DataMCPair (Selection sel, TangibleAxis tanAxis, SpectrumLoader &loaderData, SpectrumLoader &loaderMC, std::vector< const ISyst * > systs, const Var BinWeight, const Cut &bkg=kNoCut, const Cut &wrs=kNoCut)

std::string	ShortName () const

void	OverlayDataMCSyst (const int iSyst=-1) const

void	OverlayDataMCSystNorm (const int iSyst=-1) const

void	DrawMCSyst (const int iSyst=-1) const

void	DrawMCNormSyst (const int iSyst=-1) const

TLegend *	DrawLegend () const

void	AddExposure (bool IsPOT) const

void	RedrawAxes () const

void	DrawData () const

void	DrawMCBkg () const

void	DrawWrong () const

float	Purity () const

Public Attributes
Selection	fSel

TangibleAxis	fAxis

Spectrum	fData

Spectrum	fMC

Spectrum	fMCBkg

std::vector< Spectrum >	fUps

std::vector< Spectrum >	fDowns

std::string	fShortName

Spectrum	fWrong

Detailed Description

Definition at line 192 of file hadEFrac_nd_data_mc_systs.C.

Constructor & Destructor Documentation

DataMCPair::DataMCPair	(	Selection	sel,
		TangibleAxis	tanAxis,
		SpectrumLoader &	loaderData,
		SpectrumLoader &	loaderMC,
		std::vector< const ISyst * >	systs,
		const Var	BinWeight,
		const Cut &	bkg = `kNoCut`
	)

inline

Definition at line 196 of file hadEFrac_nd_data_mc_systs.C.

                                                                                         :
   fSel(sel),
   fAxis(tanAxis),
   //fData(loaderData, fAxis.fObj, fSel.fObj, kNoShift, kPPFXFluxCVWgt * kXSecCVWgt2017 * kBinWidthWeight),
   //fMC(loaderMC, fAxis.fObj, fSel.fObj, kNoShift, kPPFXFluxCVWgt * kXSecCVWgt2017 * kBinWidthWeight),
   //fMCBkg(loaderMC, fAxis.fObj, fSel.fObj && bkg, kNoShift, kPPFXFluxCVWgt * kXSecCVWgt2017 * kBinWidthWeight),
   fData (loaderData, fAxis.fObj, fSel.fObj       , kNoShift, BinWeight),
   fMC   (loaderMC  , fAxis.fObj, fSel.fObj       , kNoShift, BinWeight),
   fMCBkg(loaderMC  , fAxis.fObj, fSel.fObj && bkg, kNoShift, BinWeight),
   fShortName(fAxis.fShortName + "_" + fSel.fShortName)
   {
     fUps.reserve(systs.size());
     fDowns.reserve(systs.size());
     for(const auto& syst:systs){
       //fUps.emplace_back(loaderMC, fAxis.fObj, fSel.fObj, SystShifts(syst, +1), kPPFXFluxCVWgt * kXSecCVWgt2017 * kBinWidthWeight);
       //fDowns.emplace_back(loaderMC, fAxis.fObj, fSel.fObj, SystShifts(syst, -1), kPPFXFluxCVWgt * kXSecCVWgt2017 * kBinWidthWeight);
       fUps  .emplace_back(loaderMC, fAxis.fObj, fSel.fObj, SystShifts(syst, +1), BinWeight);
       fDowns.emplace_back(loaderMC, fAxis.fObj, fSel.fObj, SystShifts(syst, -1), BinWeight);
     }
   };

DataMCPair::DataMCPair	(	Selection	sel,
		TangibleAxis	tanAxis,
		SpectrumLoader &	loaderData,
		SpectrumLoader &	loaderMC,
		std::vector< const ISyst * >	systs,
		const Var	BinWeight,
		const Cut &	bkg = `kNoCut`,
		const Cut &	wrs = `kNoCut`
	)

inline

Definition at line 95 of file ND_DataMC_Comp_Systs_header.h.

                                                                                                                :
   fSel(sel),
     fAxis(tanAxis),
     fData (loaderData, fAxis.fObj, fSel.fObj       , kNoShift, BinWeight),
     fMC   (loaderMC  , fAxis.fObj, fSel.fObj       , kNoShift, BinWeight),
     fWrong(loaderMC  , fAxis.fObj, fSel.fObj && wrs, kNoShift, BinWeight),
     fMCBkg(loaderMC  , fAxis.fObj, fSel.fObj && bkg, kNoShift, BinWeight),
     fShortName(fAxis.fShortName + "_" + fSel.fShortName)
       {
         fUps.reserve(systs.size());
         fDowns.reserve(systs.size());
         for(const auto& syst:systs){
           fUps  .emplace_back(loaderMC, fAxis.fObj, fSel.fObj, SystShifts(syst, +1), BinWeight);
           fDowns.emplace_back(loaderMC, fAxis.fObj, fSel.fObj, SystShifts(syst, -1), BinWeight);
         }
       };

Member Function Documentation

void DataMCPair::AddExposure ( bool IsPOT ) const

inline

Definition at line 284 of file ND_DataMC_Comp_Systs_header.h.

References ana::AutoPlaceLegend(), Bless, om::cout, allTimeWatchdog::endl, Integral(), IsFHC, place(), pot, sFHC, string, and MakeMiniprodValidationCuts::text.

   {
     // Format Latex and text placement
     TLatex text;
     text.SetTextAlign(22);
     text.SetTextSize(0.04);
     // Figure out means and integrals.
     std::stringstream expo;
     std::stringstream expoData;
     std::stringstream expoMC;
     float pot = fData.POT();
     float dataMean      = (fData.ToTH1(pot)) -> GetMean();
     float dataMeanError = (fData.ToTH1(pot)) -> GetMeanError();
     float dataIntegral  = (fData.ToTH1(pot)) -> Integral( 0, (fData.ToTH1(pot))->GetNbinsX() );
     float MCMean        = (fMC  .ToTH1(pot)) -> GetMean();
     float MCMeanError   = (fMC  .ToTH1(pot)) -> GetMeanError();
     float MCIntegral    = (fMC  .ToTH1(pot)) -> Integral();
     // Decide my stringsteams.
     expoData << fixed << setprecision(3) << "Data mean: " << dataMean << ", Integral: " << dataIntegral;
     expoMC   << fixed << setprecision(3) << "Mont mean: " << MCMean   << ", Integral: " << MCIntegral;
 
     // Figure out my character strings.
     std::string sFHC = (IsFHC == true ? "Neutrino beam":"Antineutrino beam");
     std::string sPOT = (IsPOT == true ? "POT":"Area");
     // Figure out which quartile I'm looking at.
     std::string sQnt = "";
     if      ( ShortName().find("Quant1") != std::string::npos) sQnt = "Quartile 1";
     else if ( ShortName().find("Quant2") != std::string::npos) sQnt = "Quartile 2";
     else if ( ShortName().find("Quant3") != std::string::npos) sQnt = "Quartile 3";
     else if ( ShortName().find("Quant4") != std::string::npos) sQnt = "Quartile 4";
 
     TLegend* place = AutoPlaceLegend(0.0001, 0.0001, 0.54);
     if        ( !Bless || sQnt != "" ) {
       place = AutoPlaceLegend(0.0001, 0.0001, 0.525);
     } else if ( !Bless && sQnt != "" ) {
       place = AutoPlaceLegend(0.0001, 0.0001, 0.50);
     }
     // Print stuff out to the screen...
     std::cout << "\n\t Looking at " << ShortName() << ", " << sFHC << ", " << sPOT 
               << "\n\t Data Mean / Integral are; " << dataMean << ", " << dataIntegral
               << "\n\t Mont Mean / Integral are; " << MCMean   << ", " << MCIntegral
               << std::endl;
 
     // Decide what I'm putting in the legend
     std::string       Descrip = sPOT+" Normalized";
     if ( sQnt != "" ) Descrip = "#splitline{"+sPOT+" Normalized}{"+sQnt+"}";
     if ( !Bless ) {
       Descrip   = "#splitline{"+sPOT+" Normalized}{#splitline{"+expoData.str()+"}{"+expoMC.str()+"}}";
       if ( sQnt != "" )
         Descrip = "#splitline{"+sPOT+" Normalized}{#splitline{"+sQnt+"}{#splitline{"+expoData.str()+"}{"+expoMC.str()+"}}}";
     }
     // Print the legend
     text.DrawLatexNDC( place->GetX1(), place->GetY1(), Descrip.c_str() );
   }

void DataMCPair::AddExposureArea ( ) const

inline

Definition at line 447 of file hadEFrac_nd_data_mc_systs.C.

References ana::AutoPlaceLegend(), allTimeWatchdog::endl, genie::utils::style::Format(), place(), pot, and MakeMiniprodValidationCuts::text.

   {
     std::stringstream expo;
     std::stringstream expoData;
     std::stringstream expoMC;
     float pot = fData.POT();
     float dataMean = (fData.ToTH1(pot)) -> GetMean();
     float dataMeanError = (fData.ToTH1(pot)) -> GetMeanError();
     float MCMean = (fMC.ToTH1(pot)) -> GetMean();
     float MCMeanError = (fMC.ToTH1(pot)) -> GetMeanError();
     expo.precision(2);
     expoData.precision(2);
     expoMC.precision(2);
     expo  << "ND area norm., "
           << TString::Format("%.02lf #times 10^{20}", pot/1e20) << " POT " <<endl;
     //expoData  << TString::Format("Data mean: %.02lf #pm %.02lf GeV  ", dataMean, dataMeanError) <<endl;
     //expoMC    << TString::Format("MC mean:   %.02lf #pm %.02lf GeV  ", MCMean, MCMeanError) <<endl;
     expoData  << TString::Format("Data mean: %.02lf GeV  ", dataMean) <<endl;
     expoMC    << TString::Format("MC mean:   %.02lf GeV  ", MCMean  ) <<endl;
 
     TLatex text;
     TLegend* place = AutoPlaceLegend(0.0001, 0.0001, 0.52);
     text.SetTextAlign(22);
     text.SetTextSize(0.04);
 
     // text.DrawLatexNDC(place->GetX1(), place->GetY1(), 
     //                (expo.str()).c_str() );
 
     text.DrawLatexNDC(place->GetX1(), place->GetY1(), 
                       ("#splitline{" + expo.str() + "}{ #splitline{" + 
                        expoData.str() + "}{" + expoMC.str() + "}}").c_str()
                       );
 
   }

void DataMCPair::AddExposurePOT ( ) const

inline

Definition at line 482 of file hadEFrac_nd_data_mc_systs.C.

References ana::AutoPlaceLegend(), om::cout, genie::utils::style::Format(), place(), pot, and MakeMiniprodValidationCuts::text.

   {
     std::stringstream expo;
     float pot = fData.POT();
     expo.precision(2);
     expo  << "ND POT norm., "
           << TString::Format("%.02lf #times 10^{20}", pot/1e20) << " POT";
     TLatex text;
     TLegend* place = AutoPlaceLegend(0.0001, 0.0001, 0.56);
     text.SetTextAlign(22);
     text.SetTextSize(0.04);
 
     text.DrawLatexNDC(place->GetX1(), place->GetY1(), expo.str().c_str());
 
     if (ShortName().find("numuE") != std::string::npos) {
       std::cout << "\n" << ShortName() << "\t Data: " << fData.ToTH1(pot)->Integral()
                 << "\t MC: " << fMC.ToTH1(pot)->Integral() << "\n";
     }
   }

void DataMCPair::AestheticsArea ( ) const

inline

Definition at line 502 of file hadEFrac_nd_data_mc_systs.C.

References ana::CenterTitles(), and RedrawAxes().

   {
     CenterTitles();
     RedrawAxes();
     AddExposureArea();
   }

void DataMCPair::AestheticsPOT ( ) const

inline

Definition at line 509 of file hadEFrac_nd_data_mc_systs.C.

References ana::CenterTitles(), and RedrawAxes().

   {
     CenterTitles();
     RedrawAxes();
     AddExposurePOT();
   }

const char* DataMCPair::CName ( ) const

inline

Definition at line 221 of file hadEFrac_nd_data_mc_systs.C.

222 {return ShortName().c_str();};

DataMCPair::ShortName

std::string ShortName() const

Definition: hadEFrac_nd_data_mc_systs.C:219

void DataMCPair::DrawData ( ) const

inline

Definition at line 362 of file ND_DataMC_Comp_Systs_header.h.

References om::cout, DataCol, DataMar, e, allTimeWatchdog::endl, and Scale().

                         {
     
     TH1D* hData = fData.ToTH1(fData.POT());
     std::cout << "\t Data integral is " << hData->Integral() << ", mean is " << hData->GetMean() << std::endl;
 
     if (ShortName().find("NuMuE") != std::string::npos) {
       hData -> Scale(0.1,"width");
     }
     hData -> Scale( 1e-3 );
     hData->SetMarkerColor( DataCol );
     hData->SetMarkerStyle( DataMar );
     hData->SetLineColor  ( DataCol );
     hData->Draw("same");
   };

void DataMCPair::DrawData	(	const int	color = `kBlack`,
		const int	marker = `kFullCircle`
	)		const

inline

Definition at line 516 of file hadEFrac_nd_data_mc_systs.C.

References plot_validation_datamc_2018::color.

   {
     TH1D* hData = fData.ToTH1(fData.POT());
     hData->GetYaxis()->SetTitle("Events / 0.1 GeV"); // only for cust bin neutrino plot
     hData->SetMarkerColor(color);
     hData->SetMarkerStyle(marker);
     hData->SetLineColor(color);
 
     hData->Draw("SAME");
   };

TLegend* DataMCPair::DrawLegend ( ) const

inline

Definition at line 241 of file ND_DataMC_Comp_Systs_header.h.

References ana::AutoPlaceLegend(), DataCol, DataMar, FullPred, IsFHC, MECModelEnuComparisons::leg, PredErr, TotalBkg, and WrongSi.

                               {
 
     TLegend* leg = AutoPlaceLegend(0.25, 0.19, 0.68);
 
     TH1F* colMC = new TH1F();
     colMC   ->SetLineColor ( FullPred );
     TH1F* colWrong = new TH1F();
     colWrong->SetLineColor ( WrongSi  );
     colWrong->SetFillColor ( WrongSi  );
     TH1F* colMCBkg = new TH1F();
     colMCBkg->SetLineColor ( TotalBkg );
     colMCBkg->SetFillColor ( TotalBkg );
     TH1F* colData = new TH1F();
     colData->SetLineColor  ( DataCol  );
     colData->SetMarkerColor( DataCol  );
     colData->SetMarkerStyle( DataMar ) ;
 
     leg->AddEntry((TObject*)colData , "Data"      , "ple");
     leg->AddEntry((TObject*)colMC   , "Simulation", "l"  );
     TLegendEntry *entryns1=leg->AddEntry("error","1-#sigma syst. range","f");
     if (IsFHC) leg->AddEntry((TObject*)colWrong, "Wrong Sign: #bar{#nu_{#mu}} CC", "f"  );
     else       leg->AddEntry((TObject*)colWrong, "Wrong Sign: #nu_{#mu} CC"      , "f"  );
     leg->AddEntry((TObject*)colMCBkg, "Total Background", "f"  );
     entryns1->SetLineColor( PredErr );
     entryns1->SetFillStyle( 1001    );
     entryns1->SetFillColor( PredErr );
     leg->SetTextSize(0.04);
 
     leg->SetFillColor(0);
     leg->SetFillStyle(0);
 
     for(const auto& obj:*leg->GetListOfPrimitives())
       {
         if(obj->InheritsFrom("TAttFill")){
           ((TAttFill*)obj)->SetFillStyle(0);
           //((TAttFill*)obj)->SetFillColor(0);
         }
 
       }
     leg->Draw();
     return leg;
   }

TLegend* DataMCPair::DrawLegendArea	(	const int	dataColor = `kBlack`,
		std::string	dataLegendTitle = `"Data"`,
		const int	dataMarker = `kFullCircle`
	)		const

inline

Definition at line 369 of file hadEFrac_nd_data_mc_systs.C.

References ana::AutoPlaceLegend(), kBlue, kRed, ana::kTotalMCColor, and MECModelEnuComparisons::leg.

   {
 
     TLegend* leg = AutoPlaceLegend(0.36, 0.19, 0.68);
 
     TH1F* colMC = new TH1F();
     colMC->SetLineColor(kTotalMCColor);
     TH1F* colMCBkg = new TH1F();
     colMCBkg->SetLineColor(kBlue);
     TH1F* colData = new TH1F();
     colData->SetLineColor(dataColor);
     colData->SetMarkerColor(dataColor);
     colData->SetMarkerStyle(dataMarker);
 
     leg->AddEntry((TObject*)colMC, "Simulated Selected Events", "l");
     leg->AddEntry((TObject*)colMCBkg, "Simulated Background", "l");
     leg->AddEntry((TObject*)colData, dataLegendTitle.c_str(), "ple");
     TLegendEntry *entryns1=leg->AddEntry("error","Shape-only 1-#sigma syst. range","f");
     entryns1->SetFillStyle(1001);
     entryns1->SetFillColor(kRed-10);
     leg->SetTextSize(0.04);
 
     leg->SetFillColor(0);
     leg->SetFillStyle(0);
 
     for(const auto& obj:*leg->GetListOfPrimitives())
     {
       if(obj->InheritsFrom("TAttFill")){
         ((TAttFill*)obj)->SetFillStyle(0);
         //((TAttFill*)obj)->SetFillColor(0);
       }
 
     }
     leg->Draw();
     return leg;
   }

TLegend* DataMCPair::DrawLegendPOT	(	const int	dataColor = `kBlack`,
		std::string	dataLegendTitle = `"Data"`,
		const int	dataMarker = `kFullCircle`
	)		const

inline

Definition at line 408 of file hadEFrac_nd_data_mc_systs.C.

References ana::AutoPlaceLegend(), kBlue, kRed, ana::kTotalMCColor, and MECModelEnuComparisons::leg.

   {
 
     TLegend* leg = AutoPlaceLegend(0.36, 0.19, 0.68);
 
     TH1F* colMC = new TH1F();
     colMC->SetLineColor(kTotalMCColor);
     TH1F* colMCBkg = new TH1F();
     colMCBkg->SetLineColor(kBlue);
     TH1F* colData = new TH1F();
     colData->SetLineColor(dataColor);
     colData->SetMarkerColor(dataColor);
     colData->SetMarkerStyle(dataMarker);
 
     leg->AddEntry((TObject*)colMC, "Simulated Selected Events", "l");
     leg->AddEntry((TObject*)colMCBkg, "Simulated Background", "l");
     leg->AddEntry((TObject*)colData, dataLegendTitle.c_str(), "ple");
     TLegendEntry *entryns1=leg->AddEntry("error","Full 1-#sigma syst. range","f");
     entryns1->SetFillStyle(1001);
     entryns1->SetFillColor(kRed-10);
     leg->SetTextSize(0.04);
 
     leg->SetFillColor(0);
     leg->SetFillStyle(0);
 
     for(const auto& obj:*leg->GetListOfPrimitives())
     {
       if(obj->InheritsFrom("TAttFill")){
         ((TAttFill*)obj)->SetFillStyle(0);
         //((TAttFill*)obj)->SetFillColor(0);
       }
 
     }
     leg->Draw();
     return leg;
   }

void DataMCPair::DrawMCBkg ( ) const

inline

Definition at line 377 of file ND_DataMC_Comp_Systs_header.h.

References om::cout, e, allTimeWatchdog::endl, Scale(), and TotalBkg.

   {
     TH1D* hMCBkg = fMCBkg.ToTH1(fData.POT());
     std::cout << "\t MCBkg integral is " << hMCBkg->Integral() << ", mean is " << hMCBkg->GetMean() << std::endl;
 
     if (ShortName().find("NuMuE") != std::string::npos) {
       hMCBkg -> Scale(0.1,"width");
     }
     hMCBkg->Scale( 1e-3 );
     hMCBkg->SetLineColor( TotalBkg );
     hMCBkg->SetFillColor( TotalBkg );
     hMCBkg->Draw("hist same");
   };

void DataMCPair::DrawMCBkg ( ) const

inline

Definition at line 527 of file hadEFrac_nd_data_mc_systs.C.

References kBlue.

   {
     TH1D* hMCBkg = fMCBkg.ToTH1(fData.POT());
     hMCBkg->GetYaxis()->SetTitle("Events / 0.1 GeV"); // only for cust bin neutrino plot
     hMCBkg->SetLineColor(kBlue);
     hMCBkg->Draw("HIST SAME");
   };

void DataMCPair::DrawMCNormSyst ( const int iSyst = -1 ) const

inline

Definition at line 184 of file ND_DataMC_Comp_Systs_header.h.

References e, FullPred, MECModelEnuComparisons::i, ana::Spectrum::Integral(), ana::kPOT, norm, ana::Spectrum::OverridePOT(), ana::PlotWithSystErrorBand(), ana::Spectrum::POT(), PredErr, and Scale().

                                                   {
     // New vectors for area normalized systs
     std::vector<Spectrum> normUps;
     std::vector<Spectrum> normDowns;
     
     Spectrum tempMC(fMC.ToTH1(fMC.POT()),fMC.POT(),fMC.Livetime());
     tempMC.OverridePOT( fMC.POT() * fMC.ToTH1(fData.POT())->GetSumOfWeights() / fData.ToTH1(fData.POT())->GetSumOfWeights() );
     
     TH1* hNom = tempMC.ToTH1( fData.POT(), kBlack, kSolid, kPOT );
     if ( ShortName().find("NuMuE") != std::string::npos ) hNom -> Scale(0.1,"width");
     hNom -> Scale( 1e-3 );
 
     if (iSyst >= 0) {
       // We'll have one entry if there's just one syst to plot.
       normUps.reserve(1);
       normDowns.reserve(1);
     } else {
       // Otherwise we need all systs.
       normUps.reserve(fUps.size());
       normDowns.reserve(fUps.size());
     }
     
     for(size_t i = 0; i < fUps.size(); ++i) {
       // bail if we want a particular syst and it's not the one we want.
       if(iSyst >= 0 && (int)i != iSyst) continue;
       normUps.emplace_back(fUps[i]);
       Spectrum& norm = normUps.back();
       norm.OverridePOT(norm.POT() * norm.Integral(1) / tempMC.Integral(1));
     }
     
     for ( size_t i = 0; i < fDowns.size(); ++i)
       {
         // bail if we want a particular syst and it's not the one we want.
         if(iSyst >= 0 && (int)i != iSyst) continue;
         normDowns.emplace_back(fDowns[i]);
         Spectrum& norm = normDowns.back();
         norm.OverridePOT(norm.POT() * norm.Integral(1) / tempMC.Integral(1));
       }
 
     std::vector <TH1*> hUps, hDos;
     for (const auto& Up:normUps) {
       // Ups.
       TH1* hTempUp = Up.ToTH1( fData.POT(), kBlack, kSolid, kPOT );
       if ( ShortName().find("NuMuE") != std::string::npos ) hTempUp -> Scale(0.1,"width");
       hTempUp -> Scale( 1e-3 );
       hUps.push_back( hTempUp );
     }
     for (const auto& Do:normDowns) {
       // Dos.
       TH1* hTempDo = Do.ToTH1( fData.POT(), kBlack, kSolid, kPOT );
       if ( ShortName().find("NuMuE") != std::string::npos ) hTempDo -> Scale(0.1,"width");
       hTempDo -> Scale( 1e-3 );
       hDos.push_back( hTempDo );
     }
     PlotWithSystErrorBand( hNom, hUps, hDos, FullPred, PredErr, 1.5, true );
   }

void DataMCPair::DrawMCNormSyst ( const int iSyst = -1 ) const

inline

Definition at line 251 of file hadEFrac_nd_data_mc_systs.C.

References MECModelEnuComparisons::i, ana::Spectrum::Integral(), ana::kTotalMCColor, ana::kTotalMCErrorBandColor, norm, ana::Spectrum::OverridePOT(), ana::PlotWithSystErrorBand(), and ana::Spectrum::POT().

   {
         // New vectors for area normalized systs
     std::vector<Spectrum> normUps;
     std::vector<Spectrum> normDowns;
 
     Spectrum tempMC(fMC.ToTH1(fMC.POT()),fMC.POT(),fMC.Livetime());
     tempMC.OverridePOT(fMC.POT() * fMC.ToTH1(fData.POT())->GetSumOfWeights() / fData.ToTH1(fData.POT())->GetSumOfWeights());
 
     if (iSyst >= 0)
     {
       // We'll have one entry if there's just one syst to plot.
       normUps.reserve(1);
       normDowns.reserve(1);
     }
     else
     {
       // Otherwise we need all systs.
       normUps.reserve(fUps.size());
       normDowns.reserve(fUps.size());
     }
     for(size_t i = 0; i < fUps.size(); ++i)
     {
       // bail if we want a particular syst and it's not the one we want.
       if(iSyst >= 0 && (int)i != iSyst) continue;
       normUps.emplace_back(fUps[i]);
       Spectrum& norm = normUps.back();
       norm.OverridePOT(norm.POT() * norm.Integral(1) / tempMC.Integral(1));
     }
     for(size_t i = 0; i < fDowns.size(); ++i)
     {
       // bail if we want a particular syst and it's not the one we want.
       if(iSyst >= 0 && (int)i != iSyst) continue;
       normDowns.emplace_back(fDowns[i]);
       Spectrum& norm = normDowns.back();
       norm.OverridePOT(norm.POT() * norm.Integral(1) / tempMC.Integral(1));
     }
     PlotWithSystErrorBand(tempMC, normUps, normDowns, fData.POT(),
               kTotalMCColor, kTotalMCErrorBandColor, 1.5);
 
   }

void DataMCPair::DrawMCSyst ( const int iSyst = -1 ) const

inline

Definition at line 159 of file ND_DataMC_Comp_Systs_header.h.

References e, FullPred, ana::kPOT, ana::PlotWithSystErrorBand(), PredErr, and Scale().

   {
     TH1* hNom =  fMC.ToTH1(fData.POT(), kBlack, kSolid, kPOT);
     if ( ShortName().find("NuMuE") != std::string::npos ) hNom -> Scale(0.1,"width");
     hNom -> Scale( 1e-3 );
     std::vector <TH1*> hUps, hDos;
     for (const auto& Up:fUps) {
       // Ups.
       TH1* hTempUp = Up.ToTH1( fData.POT(), kBlack, kSolid, kPOT );
       if ( ShortName().find("NuMuE") != std::string::npos ) hTempUp -> Scale(0.1,"width");
       hTempUp -> Scale( 1e-3 );
       hUps.push_back( hTempUp );
     }
     for (const auto& Do:fDowns) {
       // Dos.
       TH1* hTempDo = Do.ToTH1( fData.POT(), kBlack, kSolid, kPOT );
       if ( ShortName().find("NuMuE") != std::string::npos ) hTempDo -> Scale(0.1,"width");
       hTempDo -> Scale( 1e-3 );
       hDos.push_back( hTempDo );
     }
     
     PlotWithSystErrorBand( hNom, hUps, hDos, FullPred, PredErr, 1.5, true );    
   }

void DataMCPair::DrawMCSyst ( const int iSyst = -1 ) const

inline

Definition at line 227 of file hadEFrac_nd_data_mc_systs.C.

References ana::assert(), ana::kTotalMCColor, ana::kTotalMCErrorBandColor, and ana::PlotWithSystErrorBand().

   {
     if(iSyst < 0)
       PlotWithSystErrorBand(fMC, fUps, fDowns, fData.POT(),
                 kTotalMCColor, kTotalMCErrorBandColor, 1.5);
     else
     {
       assert(iSyst < fUps.size());
       PlotWithSystErrorBand(fMC, {fUps[iSyst]}, {fDowns[iSyst]}, fData.POT(),
                 kTotalMCColor, kTotalMCErrorBandColor, 1.5);
     }
 
   }

void DataMCPair::DrawWrong ( ) const

inline

Definition at line 391 of file ND_DataMC_Comp_Systs_header.h.

References om::cout, e, allTimeWatchdog::endl, Scale(), and WrongSi.

   {
     TH1D* hWrong = fWrong.ToTH1(fData.POT());
     std::cout << "\t hWrong integral is " << hWrong->Integral() << ", mean is " << hWrong->GetMean() << std::endl;
 
     if (ShortName().find("NuMuE") != std::string::npos) {
       hWrong -> Scale(0.1,"width");
     }
     hWrong->Scale( 1e-3 );
     hWrong->SetLineColor( WrongSi );
     hWrong->SetFillColor( WrongSi );
     hWrong->Draw("hist same");
   };

void DataMCPair::OverlayDataMCSyst ( const int iSyst = -1 ) const

inline

Definition at line 118 of file ND_DataMC_Comp_Systs_header.h.

References ana::CenterTitles(), DataCol, DataMar, Draw(), DrawLegend(), MECModelEnuComparisons::leg, and RedrawAxes().

                                                      {
     
     TH1D* hDummy = fData.ToTH1(fData.POT());
     hDummy->SetMarkerColor(DataCol);
     hDummy->SetMarkerStyle(DataMar);
     hDummy->SetLineColor  (DataCol);
     hDummy->Draw();
 
     DrawMCSyst(iSyst);
     DrawMCBkg();
     DrawWrong();
     DrawData();
     hDummy -> Draw( "axis same" );
     TLegend* leg = DrawLegend();
     CenterTitles();
     RedrawAxes();
     AddExposure( true );
   };

void DataMCPair::OverlayDataMCSyst ( const int iSyst = -1 ) const

inline

Definition at line 240 of file hadEFrac_nd_data_mc_systs.C.

References MECModelEnuComparisons::leg.

   {
     DrawData();
     DrawMCSyst(iSyst);
     DrawData();
     DrawMCBkg();
     TLegend* leg = DrawLegendPOT();
     AestheticsPOT();
 
   };

void DataMCPair::OverlayDataMCSystExtraData	(	const DataMCPair &	extraData,
		const std::string	dataLegendTitle,
		const std::string	extraLegendTitle,
		const int	iSyst = `-1`,
		const int	dataColor = `kBlack`,
		const int	extraColor = `kGray+1`,
		const int	dataMarker = `kFullCircle`,
		const int	extraMarker = `kFullCircle`
	)		const

inline

Definition at line 316 of file hadEFrac_nd_data_mc_systs.C.

References ana::DataMCPair::DrawData(), and MECModelEnuComparisons::leg.

   {
 
     DrawMCSyst(iSyst);
     DrawData(dataColor, dataMarker);
     extraData.DrawData(extraColor, extraMarker);
     DrawMCBkg();
     TLegend* leg = DrawLegendPOT(dataColor, dataLegendTitle, dataMarker);
     TH1F* colExtraData = new TH1F();
     colExtraData->SetMarkerColor(extraColor);
     colExtraData->SetMarkerStyle(extraMarker);
     colExtraData->SetLineColor(extraColor);
     leg->AddEntry((TObject*)colExtraData, extraLegendTitle.c_str(), "ple");
     AestheticsPOT();
 
   };

void DataMCPair::OverlayDataMCSystNorm ( const int iSyst = -1 ) const

inline

Definition at line 138 of file ND_DataMC_Comp_Systs_header.h.

References ana::CenterTitles(), DataCol, DataMar, DrawLegend(), MECModelEnuComparisons::leg, and RedrawAxes().

                                                          {
 
     TH1D* hDummy = fData.ToTH1(fData.POT());
     hDummy->SetMarkerColor(DataCol);
     hDummy->SetMarkerStyle(DataMar);
     hDummy->SetLineColor  (DataCol);
     hDummy->Draw();
 
     DrawMCNormSyst(iSyst);
     DrawMCBkg();
     DrawWrong();
     DrawData();
     hDummy->Draw( "axis same" );
     TLegend* leg = DrawLegend();
     CenterTitles();
     RedrawAxes();
     AddExposure( false );
   };

void DataMCPair::OverlayDataMCSystNorm ( const int iSyst = -1 ) const

inline

Definition at line 294 of file hadEFrac_nd_data_mc_systs.C.

References MECModelEnuComparisons::leg.

   {
 
     const int tempcolor=kBlack;
     const int tempmarker=kFullCircle;
     TH1D* hDummy = fData.ToTH1(fData.POT());
     hDummy->GetYaxis()->SetTitle("Events / 0.1 GeV"); // only for cust bin neutrino plot
     //hDummy->GetYaxis()->SetTitleOffset(1.1);
     hDummy->SetMarkerColor(tempcolor);
     hDummy->SetMarkerStyle(tempmarker);
     hDummy->SetLineColor(tempcolor);
     hDummy->Draw();
 
     DrawMCNormSyst(iSyst);
     DrawData();
     DrawMCBkg();
     TLegend* leg = DrawLegendArea();
     AestheticsArea();
   };

void DataMCPair::OverlayDataMCSystNormExtraData	(	const DataMCPair &	extraData,
		const std::string	dataLegendTitle,
		const std::string	extraLegendTitle,
		const int	iSyst = `-1`,
		const int	dataColor = `kBlack`,
		const int	extraColor = `kGray+1`,
		const int	dataMarker = `kFullCircle`,
		const int	extraMarker = `kFullCircle`
	)		const

inline

Definition at line 343 of file hadEFrac_nd_data_mc_systs.C.

References ana::DataMCPair::DrawData(), and MECModelEnuComparisons::leg.

   {
 
 
     DrawMCNormSyst(iSyst);
     DrawData(dataColor, dataMarker);
     extraData.DrawData(extraColor);
     DrawMCBkg();
     TLegend* leg = DrawLegendArea(dataColor, dataLegendTitle, dataMarker);
     TH1F* colExtraData = new TH1F();
     colExtraData->SetMarkerColor(extraColor);
     colExtraData->SetMarkerStyle(extraMarker);
     colExtraData->SetLineColor(extraColor);
     leg->AddEntry((TObject*)colExtraData, extraLegendTitle.c_str(), "ple");
     AestheticsArea();
 
   };

float DataMCPair::Purity ( ) const

inline

Definition at line 407 of file ND_DataMC_Comp_Systs_header.h.

References GetEntries().

   {
     return 1 - fMCBkg.ToTH1(fData.POT())->GetEntries() /
       fMC.ToTH1(fData.POT())->GetEntries();
   }

float DataMCPair::Purity ( ) const

inline

Definition at line 535 of file hadEFrac_nd_data_mc_systs.C.

References GetEntries().

   {
     return 1 - fMCBkg.ToTH1(fData.POT())->GetEntries() /
                   fMC.ToTH1(fData.POT())->GetEntries();
         }

void DataMCPair::RedrawAxes ( ) const

inline

Definition at line 339 of file ND_DataMC_Comp_Systs_header.h.

References string.

   {
     // --- Loop through the objects in the pad. 
     for(const auto& obj:*(gPad->GetListOfPrimitives())) {
       // --- Want to repaint all TH1's.
       if(obj->InheritsFrom("TH1")) {
         // --- Check if I have the Reco Energy canvas?      
         std::string XTitle = ((TH1*)obj)->GetXaxis()->GetTitle();
         std::stringstream YAxisStr;
         YAxisStr << "10^{3} Events / " << Form("%.2f#times10^{20} POT",fData.POT()/1e20);
         if ( XTitle.find("Reconstructed Neutrino Energ") != std::string::npos ) {
           YAxisStr << " / 0.1 GeV";
         }
         ((TH1*)obj)->GetYaxis()->SetTitle( YAxisStr.str().c_str() );
         // Repaint the axes.
         ((TH1*)obj)->GetXaxis()->Paint();
         //((TH1*)obj)->GetYaxis()->SetTitleOffset(1.1);
         ((TH1*)obj)->GetYaxis()->Paint();
         ((TH1*)obj)->GetZaxis()->Paint();
       }
     }
   }

std::string DataMCPair::ShortName ( ) const

inline

Definition at line 114 of file ND_DataMC_Comp_Systs_header.h.

References ana::Tangible< T >::fShortName.

115 {return fShortName;};

DataMCPair::fShortName

std::string fShortName

Definition: hadEFrac_nd_data_mc_systs.C:548

std::string DataMCPair::ShortName ( ) const

inline

Definition at line 219 of file hadEFrac_nd_data_mc_systs.C.

References ana::Tangible< T >::fShortName.

220 {return fShortName;};

DataMCPair::fShortName

std::string fShortName

Definition: hadEFrac_nd_data_mc_systs.C:548

Member Data Documentation

TangibleAxis DataMCPair::fAxis

Definition at line 542 of file hadEFrac_nd_data_mc_systs.C.

Spectrum DataMCPair::fData

Definition at line 543 of file hadEFrac_nd_data_mc_systs.C.

std::vector<Spectrum> DataMCPair::fDowns

Definition at line 547 of file hadEFrac_nd_data_mc_systs.C.

Spectrum DataMCPair::fMC

Definition at line 544 of file hadEFrac_nd_data_mc_systs.C.

Spectrum DataMCPair::fMCBkg

Definition at line 545 of file hadEFrac_nd_data_mc_systs.C.

Selection DataMCPair::fSel

Definition at line 541 of file hadEFrac_nd_data_mc_systs.C.

std::string DataMCPair::fShortName

Definition at line 548 of file hadEFrac_nd_data_mc_systs.C.

std::vector<Spectrum> DataMCPair::fUps

Definition at line 546 of file hadEFrac_nd_data_mc_systs.C.

Spectrum DataMCPair::fWrong

Definition at line 418 of file ND_DataMC_Comp_Systs_header.h.

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

/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-24/3FlavorAna/Ana2017/numu/hadEFrac_nd_data_mc_systs.C
/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-24/3FlavorAna/Ana2018/numu/ND_DataMC_Comp_Systs_header.h

Public Member Functions

Public Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation