#include <ThinSliceSample.h>

Public Member Functions
	ThinSliceSample (std::string name, int flux_type, const std::vector< fhicl::ParameterSet > &selections, const std::vector< double > &incident_bins, const std::vector< double > &true_incident_bins, size_t beam_energy_bin, bool is_signal=false, std::pair< double, double > range={0., 0.})

	~ThinSliceSample ()

void	SetFactor (double f)

const std::map< int, TH1 * > &	GetSelectionHists () const

const std::map< int, std::vector< TH1 * > > &	GetShifts (std::string syst_name) const

const std::map< int, std::vector< TSpline3 * > > &	GetSplines (std::string syst_name) const

const std::map< std::string, std::map< int, std::vector< TSpline3 * > > > &	GetAllSplines () const

void	SetSystematicSplines (const std::map< std::string, std::map< int, std::vector< TSpline3 * >>> &input)

double	GetSplineWeight (std::string syst_name, double par_val, int selection_ID, double val) const

void	AddSystematicShift (TH1 *hist, std::string syst_name, int selection_ID)

void	SaveSystematics (std::string syst_name, TDirectory *dir)

const std::map< int, TH1 * > &	GetRebinnedSelectionHists () const

TH1 *	GetSelectionHist (int id)

TH1D &	GetTrueIncidentHist ()

TH1 *	GetRebinnedSelectionHist (int id)

const std::string &	GetName () const

const int &	GetFluxType () const

const double &	GetNominalFlux () const

const double &	GetVariedFlux () const

void	AddFlux (double val=1.)

void	AddVariedFlux (double val=1.)

void	FillSystematicShift (std::string syst_name, int selection_ID, const std::vector< double > &vals)

void	FillSystematicShift (std::string syst_name, int selection_ID, const std::vector< double > &vals, const std::vector< double > &weights)

void	SetSystematicVals (std::string syst_name, std::vector< double > &vals)

void	MakeSystematicSplines (std::string syst_name)

void	AddIncidentEnergies (const std::vector< double > &vals, double weight=1.)

void	AddESliceEnergies (const std::pair< double, double > &vals, double weight=1.)

void	FillTrueIncidentHist (const std::vector< double > &vals, double weight=1.)

void	FillSelectionHist (int id, double val, double weight=1.)

template<size_t N>
void	FillSelectionHist (int id, const double(&vals)[N], double weight=1.)

void	FillHistFromIncidentEnergies (TH1D &hist)

void	FillESliceHist (TH1D &hist)

void	ScaleHists (double val)

void	ScaleIncidentEnergies (double val)

void	ScaleESliceEnergies (double val)

void	Reset ()

void	ScaleVariedFlux (double val)

void	ScaleToDataMC ()

void	SetDataMCScale (double val)

void	SetFactorAndScale (double val)

void	ExtraFactor (double val)

double	GetFactor ()

void	ResetFactor ()

void	SetFactorToBestFit ()

double	GetBestFitFactor ()

void	SetBestFit ()

bool	CheckIsSignal ()

bool	CheckInSignalRange (double val)

double	RangeLowEnd ()

double	RangeHighEnd ()

const std::pair< double, double > &	GetRange () const

void	RefillRebinnedHists ()

void	MakeRebinnedHists ()

Private Member Functions
void	Rebin1D (TH1 sel_hist, TH1 rebinned)

void	Rebin2D (TH1 sel_hist, TH1 rebinned)

void	Rebin3D (TH1 sel_hist, TH1 rebinned)

Private Attributes
double	fFactor = 1.

double	fBestFitFactor = 1.

bool	fBestFitIsSet = false

std::string	fSampleName

int	fFluxType

double	fNominalFlux = 0.

double	fVariedFlux = 0.

double	fDataMCScale = 1.

bool	fIsSignal

std::pair< double, double >	fRange

std::map< int, TH1 * >	fSelectionHists

TH1D	fTrueIncidentHist

std::map< int, TH1 * >	fSelectionHistsRebinned

bool	fMadeRebinned = false

std::map< std::string, std::map< int, std::vector< TSpline3 * > > >	fSystematicSplines

std::map< std::string, std::map< int, std::vector< TH1 * > > >	fSystematicShifts

std::map< std::string, std::vector< double > >	fSystematicVals

std::vector< std::pair< double, double > >	fIncidentEnergies

std::vector< std::pair< std::pair< double, double >, double > >	fESliceEnergies

Detailed Description

Definition at line 23 of file ThinSliceSample.h.

Constructor & Destructor Documentation

protoana::ThinSliceSample::ThinSliceSample	(	std::string	name,
		int	flux_type,
		const std::vector< fhicl::ParameterSet > &	selections,
		const std::vector< double > &	incident_bins,
		const std::vector< double > &	true_incident_bins,
		size_t	beam_energy_bin,
		bool	is_signal = `false`,
		std::pair< double, double >	range = `{0., 0.}`
	)

Definition at line 10 of file ThinSliceSample.cxx.

     : fSampleName(name),
       fFluxType(flux_type),
       fIsSignal(is_signal),
       fRange(range) {
 
   std::string inc_name = "";
   std::string title = name +
                       (is_signal ?
                           ("(" + protoana::PreciseToString(range.first) + " "
                            + protoana::PreciseToString(range.second) + ")") :
                            "") +
                       ";Reconstructed KE (MeV)";
   if (is_signal) {
     inc_name = "sample_" + name + "_" +
                protoana::PreciseToString(range.first) + "_" +
                protoana::PreciseToString(range.second) + "_" +
                "_incident_hist_" +
                std::to_string(beam_energy_bin);
   }
   else {
     inc_name = "sample_" + name + "_incident_hist_" + std::to_string(beam_energy_bin);
   }
   inc_name += ("_" + std::to_string(beam_energy_bin));
   //fIncidentHist = TH1D(inc_name.c_str(), title.c_str(), incident_bins.size() - 1,
   //                     &incident_bins[0]);
 
   inc_name += "_true";
   fTrueIncidentHist = TH1D(inc_name.c_str(), title.c_str(),
                            true_incident_bins.size() - 1,
                            &true_incident_bins[0]);
 
   for (auto it = selections.begin(); it != selections.end(); ++it) {
     std::string sel_name = "";
     if (is_signal) {
       sel_name = "sample_" + name + "_" +
                  protoana::PreciseToString(range.first) + "_" +
                  protoana::PreciseToString(range.second) + "_selected_" +
                  it->get<std::string>("Name") + "_hist";
     }
     else {
       sel_name = "sample_" + name + "_selected_" +
                  it->get<std::string>("Name") + "_hist";
     }
     sel_name += "_" + std::to_string(beam_energy_bin);
 
     std::vector<std::vector<double>> selected_bins
         = it->get<std::vector<std::vector<double>>>("RecoBins");
     if (selected_bins.size() == 1) {
       fSelectionHists[it->get<int>("ID")] = new TH1D(
           sel_name.c_str(), title.c_str(), selected_bins[0].size() - 1,
           &selected_bins[0][0]);
     }
     else if (selected_bins.size() == 2) {
       fSelectionHists[it->get<int>("ID")] = new TH2D(
           sel_name.c_str(), title.c_str(),
           selected_bins[0].size() - 1, &selected_bins[0][0],
           selected_bins[1].size() - 1, &selected_bins[1][0]);
     }
     else if (selected_bins.size() == 3) {
       fSelectionHists[it->get<int>("ID")] = new TH3D(
           sel_name.c_str(), title.c_str(),
           selected_bins[0].size() - 1, &selected_bins[0][0],
           selected_bins[1].size() - 1, &selected_bins[1][0],
           selected_bins[2].size() - 1, &selected_bins[2][0]);
     }
     /*else {
       throw
     }*/
   }
   MakeRebinnedHists();
 }

protoana::ThinSliceSample::~ThinSliceSample ( )

inline

Definition at line 32 of file ThinSliceSample.h.

32 {};

Member Function Documentation

void protoana::ThinSliceSample::AddESliceEnergies	(	const std::pair< double, double > &	vals,
		double	weight = `1.`
	)

inline

Definition at line 213 of file ThinSliceSample.h.

                                                                                    {
     fESliceEnergies.push_back({vals, weight});
   };

void protoana::ThinSliceSample::AddFlux ( double val = 1. )

inline

Definition at line 131 of file ThinSliceSample.h.

                                 {
     fNominalFlux += val;
     fVariedFlux += val;
   };

void protoana::ThinSliceSample::AddIncidentEnergies	(	const std::vector< double > &	vals,
		double	weight = `1.`
	)

inline

Definition at line 208 of file ThinSliceSample.h.

                                                                                {
     for (auto v : vals)
       fIncidentEnergies.push_back({v, weight});
   };

void protoana::ThinSliceSample::AddSystematicShift	(	TH1 *	hist,
		std::string	syst_name,
		int	selection_ID
	)

inline

Definition at line 65 of file ThinSliceSample.h.

                                             {
     fSystematicShifts[syst_name][selection_ID].push_back(hist);
   };

void protoana::ThinSliceSample::AddVariedFlux ( double val = 1. )

inline

Definition at line 136 of file ThinSliceSample.h.

                                       {
     fVariedFlux += val;
   }

bool protoana::ThinSliceSample::CheckInSignalRange ( double val )

inline

Definition at line 363 of file ThinSliceSample.h.

363 {return ((fRange.first < val) &&

364 (val <= fRange.second));};

val

Definition: registry_via_id_test_2.cc:15

protoana::ThinSliceSample::fRange

std::pair< double, double > fRange

Definition: ThinSliceSample.h:382

bool protoana::ThinSliceSample::CheckIsSignal ( )

inline

Definition at line 362 of file ThinSliceSample.h.

362 {return fIsSignal;};

protoana::ThinSliceSample::fIsSignal

bool fIsSignal

Definition: ThinSliceSample.h:381

void protoana::ThinSliceSample::ExtraFactor ( double val )

inline

Definition at line 331 of file ThinSliceSample.h.

                                {
     SetFactorAndScale(val*fFactor);
   };

void protoana::ThinSliceSample::FillESliceHist ( TH1D & hist )

inline

Definition at line 257 of file ThinSliceSample.h.

                                    {
     for (auto e : fESliceEnergies) {
       std::pair<double, double> vals = e.first;
       double w = e.second;
       int last_bin = hist.FindBin(vals.first);
       int first_bin = hist.FindBin(vals.second);
       for (int i = first_bin; i <= last_bin; ++i)
         //hist.AddBinContent(i, fFactor);
         hist.AddBinContent(i, fFactor*w);
     }
   };

void protoana::ThinSliceSample::FillHistFromIncidentEnergies ( TH1D & hist )

inline

Definition at line 250 of file ThinSliceSample.h.

                                                  {
     for (auto vals : fIncidentEnergies) {
       //hist.Fill(vals.first, fFactor/*vals.second*/);
       hist.Fill(vals.first, fFactor*vals.second);
     }
   };

void protoana::ThinSliceSample::FillSelectionHist	(	int	id,
		double	val,
		double	weight = `1.`
	)

inline

Definition at line 223 of file ThinSliceSample.h.

                                                                  {
     if (fSelectionHists.find(id) != fSelectionHists.end()) {
       fSelectionHists.at(id)->Fill(val, weight);
     }
   };

template<size_t N>

void protoana::ThinSliceSample::FillSelectionHist	(	int	id,
		const double(&)	vals[N],
		double	weight = `1.`
	)

inline

Definition at line 229 of file ThinSliceSample.h.

                                                                  {
     if (N < 1 || N > 3) {
       std::string message = "Error: trying to fill hists with too many values";
       message += std::to_string(N);
       throw std::runtime_error(message);
     }
 
     if (fSelectionHists.find(id) != fSelectionHists.end()) {
       if (N == 1) {
         fSelectionHists.at(id)->Fill(vals[0], weight);
       }
       else if (N == 2) {
         ((TH2D*)fSelectionHists.at(id))->Fill(vals[0], vals[1], weight);
       }
       else if (N == 3) {
         ((TH3D*)fSelectionHists.at(id))->Fill(vals[0], vals[1], vals[2], weight);
       }
     }
   }

void protoana::ThinSliceSample::FillSystematicShift	(	std::string	syst_name,
		int	selection_ID,
		const std::vector< double > &	vals
	)

inline

Definition at line 140 of file ThinSliceSample.h.

                                                            {
     if (vals.size() !=
         fSystematicShifts[syst_name][selection_ID].size()) {
       std::string message = "ThinSliceSample: Input systematic shift values and number of shift hists differ"; 
       throw std::runtime_error(message);
     }
 
     for (size_t i = 0; i < vals.size(); ++i) {
       fSystematicShifts[syst_name][selection_ID][i]->Fill(vals[i]);
     }
   };

void protoana::ThinSliceSample::FillSystematicShift	(	std::string	syst_name,
		int	selection_ID,
		const std::vector< double > &	vals,
		const std::vector< double > &	weights
	)

inline

Definition at line 154 of file ThinSliceSample.h.

                                                               {
     if (vals.size() !=
         fSystematicShifts[syst_name][selection_ID].size()) {
       std::string message = "ThinSliceSample: Input systematic shift values and number of shift hists differ"; 
       throw std::runtime_error(message);
     }
 
     for (size_t i = 0; i < vals.size(); ++i) {
       fSystematicShifts[syst_name][selection_ID][i]->Fill(vals[i], weights[i]);
     }
   };

void protoana::ThinSliceSample::FillTrueIncidentHist	(	const std::vector< double > &	vals,
		double	weight = `1.`
	)

inline

Definition at line 217 of file ThinSliceSample.h.

                                                                                 {
     for (size_t i = 0; i < vals.size(); ++i) {
       fTrueIncidentHist.Fill(vals.at(i), weight);
     }
   };

const std::map<std::string, std::map<int, std::vector<TSpline3 *> > >& protoana::ThinSliceSample::GetAllSplines ( ) const

inline

Definition at line 51 of file ThinSliceSample.h.

                             {
     return fSystematicSplines;
   };

double protoana::ThinSliceSample::GetBestFitFactor ( )

inline

Definition at line 349 of file ThinSliceSample.h.

                             {
     return fBestFitFactor;
   };

double protoana::ThinSliceSample::GetFactor ( void )

inline

Definition at line 335 of file ThinSliceSample.h.

                      {
     return fFactor;
   };

const int& protoana::ThinSliceSample::GetFluxType ( ) const

inline

Definition at line 119 of file ThinSliceSample.h.

                                   {
     return fFluxType;
   };

const std::string& protoana::ThinSliceSample::GetName ( void ) const

inline

Definition at line 115 of file ThinSliceSample.h.

                                     {
     return fSampleName;
   };

const double& protoana::ThinSliceSample::GetNominalFlux ( ) const

inline

Definition at line 123 of file ThinSliceSample.h.

                                         {
     return fNominalFlux;
   };

const std::pair<double, double>& protoana::ThinSliceSample::GetRange ( ) const

inline

Definition at line 368 of file ThinSliceSample.h.

368 {return fRange;};

protoana::ThinSliceSample::fRange

std::pair< double, double > fRange

Definition: ThinSliceSample.h:382

TH1* protoana::ThinSliceSample::GetRebinnedSelectionHist ( int id )

inline

Definition at line 111 of file ThinSliceSample.h.

                                          {
     return fSelectionHistsRebinned.at(id);
   };

const std::map<int, TH1 *>& protoana::ThinSliceSample::GetRebinnedSelectionHists ( ) const

inline

Definition at line 89 of file ThinSliceSample.h.

                                                                {
     return fSelectionHistsRebinned;
   };

TH1* protoana::ThinSliceSample::GetSelectionHist ( int id )

inline

Definition at line 93 of file ThinSliceSample.h.

                                  {
     return fSelectionHists.at(id);
   };

const std::map<int, TH1 *>& protoana::ThinSliceSample::GetSelectionHists ( ) const

inline

Definition at line 36 of file ThinSliceSample.h.

                                                        {
     return fSelectionHists;
   };

const std::map<int, std::vector<TH1 *> >& protoana::ThinSliceSample::GetShifts ( std::string syst_name ) const

inline

Definition at line 41 of file ThinSliceSample.h.

                                            {
     return fSystematicShifts.at(syst_name); 
   };

const std::map<int, std::vector<TSpline3 *> >& protoana::ThinSliceSample::GetSplines ( std::string syst_name ) const

inline

Definition at line 46 of file ThinSliceSample.h.

                                             {
     return fSystematicSplines.at(syst_name);
   };

double protoana::ThinSliceSample::GetSplineWeight	(	std::string	syst_name,
		double	par_val,
		int	selection_ID,
		double	val
	)		const

inline

Definition at line 59 of file ThinSliceSample.h.

                                                              {
     int bin = fSelectionHists.at(selection_ID)->FindBin(val);
     return fSystematicSplines.at(syst_name).at(selection_ID).at(bin-1)->Eval(par_val);
   };

TH1D& protoana::ThinSliceSample::GetTrueIncidentHist ( )

inline

Definition at line 102 of file ThinSliceSample.h.

                                {
     return fTrueIncidentHist;
   };

const double& protoana::ThinSliceSample::GetVariedFlux ( ) const

inline

Definition at line 127 of file ThinSliceSample.h.

                                        {
     return fVariedFlux;
   };

void protoana::ThinSliceSample::MakeRebinnedHists ( )

Definition at line 89 of file ThinSliceSample.cxx.

                                                 {
   if (!fMadeRebinned) {
     //std::string inc_name = fIncidentHist.GetName();
     //inc_name += "Rebinned";
     //fIncidentHistRebinned = TH1D(inc_name.c_str(), fIncidentHist.GetTitle(),
     //                             fIncidentHist.GetNbinsX(), 0, fIncidentHist.GetNbinsX());
     //for (int i = 1; i <= fIncidentHist.GetNbinsX(); ++i) {
     //  fIncidentHistRebinned.SetBinContent(i, fIncidentHist.GetBinContent(i));
 
     //  double low_edge = fIncidentHist.GetXaxis()->GetBinLowEdge(i);
     //  double up_edge = fIncidentHist.GetXaxis()->GetBinUpEdge(i);
     //  std::string bin_label = (low_edge < 0. ? "< 0." :
     //                           (protoana::PreciseToString(low_edge, 0) + " - " +
     //                            protoana::PreciseToString(up_edge, 0)));
     //  fIncidentHistRebinned.GetXaxis()->SetBinLabel(i, bin_label.c_str());
     //}
 
     for (auto it = fSelectionHists.begin(); it != fSelectionHists.end(); ++it) {
       TH1 * sel_hist = (TH1 *)it->second;
       std::string name = sel_hist->GetName();
       name += "Rebinned";
       
       size_t nAxes = 1;
       if (sel_hist->GetNbinsY() > 1) ++nAxes;
       if (sel_hist->GetNbinsZ() > 1) ++nAxes;
 
       if (nAxes == 1) {
         fSelectionHistsRebinned[it->first] = new TH1D(
             name.c_str(), sel_hist->GetTitle(),
             sel_hist->GetNbinsX(), 0, sel_hist->GetNbinsX());
         Rebin1D(sel_hist, fSelectionHistsRebinned[it->first]);
       }
       else if (nAxes == 2) {
         fSelectionHistsRebinned[it->first] = new TH2D(
             name.c_str(), sel_hist->GetTitle(),
             sel_hist->GetNbinsX(), 0, sel_hist->GetNbinsX(),
             sel_hist->GetNbinsY(), 0, sel_hist->GetNbinsY());
         Rebin2D(sel_hist, fSelectionHistsRebinned[it->first]);
       }
       else if (nAxes == 3) {
         fSelectionHistsRebinned[it->first] = new TH3D(
             name.c_str(), sel_hist->GetTitle(),
             sel_hist->GetNbinsX(), 0, sel_hist->GetNbinsX(),
             sel_hist->GetNbinsY(), 0, sel_hist->GetNbinsY(),
             sel_hist->GetNbinsZ(), 0, sel_hist->GetNbinsZ());
         Rebin3D(sel_hist, fSelectionHistsRebinned[it->first]);
       }
     }
 
     fMadeRebinned = true;
   }
 }

void protoana::ThinSliceSample::MakeSystematicSplines ( std::string syst_name )

inline

Definition at line 173 of file ThinSliceSample.h.

                                                    {
     for (auto it2 = fSystematicShifts[syst_name].begin();
          it2 != fSystematicShifts[syst_name].end(); ++it2) {
       int selection_ID = it2->first;
       std::vector<TH1*> hists = it2->second;
       TH1D * selection_hist = (TH1D*)fSelectionHists[selection_ID];
       for (int i = 1; i <= selection_hist->GetNbinsX(); ++i) {
         std::vector<double> vars;
         for (size_t j = 0; j < hists.size(); ++j) {
           if (selection_hist->GetBinContent(i) < 1.e-5) {
             vars.push_back(1.);
           }
           else {
             vars.push_back(
                 hists[j]->GetBinContent(i)/selection_hist->GetBinContent(i));
           }
         }
         //if (do_insert)
           vars.insert(vars.begin() + vars.size()/2, 1.);
         std::string spline_name = selection_hist->GetName();
         spline_name += "_" + syst_name + "_Spline_" + std::to_string(i);
         fSystematicSplines[syst_name][selection_ID].push_back(
             new TSpline3(spline_name.c_str(), &fSystematicVals[syst_name][0],
                          &vars[0], vars.size()));
       }
     }
   };

double protoana::ThinSliceSample::RangeHighEnd ( )

inline

Definition at line 366 of file ThinSliceSample.h.

366 {return fRange.second;};

protoana::ThinSliceSample::fRange

std::pair< double, double > fRange

Definition: ThinSliceSample.h:382

double protoana::ThinSliceSample::RangeLowEnd ( )

inline

Definition at line 365 of file ThinSliceSample.h.

365 {return fRange.first;};

protoana::ThinSliceSample::fRange

std::pair< double, double > fRange

Definition: ThinSliceSample.h:382

void protoana::ThinSliceSample::Rebin1D	(	TH1 *	sel_hist,
		TH1 *	rebinned
	)

private

Definition at line 142 of file ThinSliceSample.cxx.

                                                                     {
   for (int i = 1; i <= sel_hist->GetNbinsX(); ++i) {
     double low_x = sel_hist->GetXaxis()->GetBinLowEdge(i);
     double up_x = sel_hist->GetXaxis()->GetBinUpEdge(i);
     std::string bin_label = (low_x < 0. ? "< 0." :
                              (protoana::PreciseToString(low_x, 0) + " - " +
                               protoana::PreciseToString(up_x, 0)));
     rebinned->GetXaxis()->SetBinLabel(i, bin_label.c_str());
 
     rebinned->SetBinContent(i, sel_hist->GetBinContent(i));
   }
 }

void protoana::ThinSliceSample::Rebin2D	(	TH1 *	sel_hist,
		TH1 *	rebinned
	)

private

Definition at line 155 of file ThinSliceSample.cxx.

                                                                     {
   for (int i = 1; i <= sel_hist->GetNbinsX(); ++i) {
     double low_x = sel_hist->GetXaxis()->GetBinLowEdge(i);
     double up_x = sel_hist->GetXaxis()->GetBinUpEdge(i);
     std::string bin_label = (low_x < 0. ? "< 0." :
                              (protoana::PreciseToString(low_x, 0) + " - " +
                               protoana::PreciseToString(up_x, 0)));
     rebinned->GetXaxis()->SetBinLabel(
         i, bin_label.c_str());
     for (int j = 1; j <= sel_hist->GetNbinsY(); ++j) {
       double low_y = sel_hist->GetYaxis()->GetBinLowEdge(j);
       double up_y = sel_hist->GetYaxis()->GetBinUpEdge(j);
       std::string y_label = (low_y < 0. ? "< 0." :
                              (protoana::PreciseToString(low_y, 0) + " - " +
                               protoana::PreciseToString(up_y, 0)));
       rebinned->GetYaxis()->SetBinLabel(j, bin_label.c_str());
       rebinned->SetBinContent(i, j, sel_hist->GetBinContent(i, j));
     }
   }
 }

void protoana::ThinSliceSample::Rebin3D	(	TH1 *	sel_hist,
		TH1 *	rebinned
	)

private

Definition at line 176 of file ThinSliceSample.cxx.

                                                                     {
   for (int i = 1; i <= sel_hist->GetNbinsX(); ++i) {
     double low_x = sel_hist->GetXaxis()->GetBinLowEdge(i);
     double up_x = sel_hist->GetXaxis()->GetBinUpEdge(i);
     std::string bin_label = (low_x < 0. ? "< 0." :
                              (protoana::PreciseToString(low_x, 0) + " - " +
                               protoana::PreciseToString(up_x, 0)));
     rebinned->GetXaxis()->SetBinLabel(i, bin_label.c_str());
     for (int j = 1; j <= sel_hist->GetNbinsY(); ++j) {
       double low_y = sel_hist->GetYaxis()->GetBinLowEdge(j);
       double up_y = sel_hist->GetYaxis()->GetBinUpEdge(j);
       std::string y_label = (low_y < 0. ? "< 0." :
                              (protoana::PreciseToString(low_y, 0) + " - " +
                               protoana::PreciseToString(up_y, 0)));
       rebinned->GetYaxis()->SetBinLabel(j, bin_label.c_str());
 
       for (int k = 1; k <= sel_hist->GetNbinsY(); ++k) {
         double low_z = sel_hist->GetYaxis()->GetBinLowEdge(k);
         double up_z = sel_hist->GetYaxis()->GetBinUpEdge(k);
         std::string y_label = (low_z < 0. ? "< 0." :
                                (protoana::PreciseToString(low_z, 0) + " - " +
                                 protoana::PreciseToString(up_z, 0)));
         rebinned->GetZaxis()->SetBinLabel(k, bin_label.c_str());
 
         rebinned->SetBinContent(i, j, k, sel_hist->GetBinContent(i, j, k));
       }
     }
   }
 }

void protoana::ThinSliceSample::RefillRebinnedHists ( )

Definition at line 206 of file ThinSliceSample.cxx.

                                                   {
   //for (int i = 1; i <= fIncidentHist.GetNbinsX(); ++i) {
   //  fIncidentHistRebinned.SetBinContent(i, fIncidentHist.GetBinContent(i));
   //}
   for (auto it = fSelectionHistsRebinned.begin();
        it != fSelectionHistsRebinned.end(); ++it) {
     for (int i = 1; i <= it->second->GetNbinsX(); ++i) {
       it->second->SetBinContent(i, fSelectionHists[it->first]->GetBinContent(i));
     }
   }
 }

void protoana::ThinSliceSample::Reset ( void )

inline

Definition at line 292 of file ThinSliceSample.h.

                {
     fVariedFlux = 0.;
     for (auto it = fSelectionHists.begin(); it != fSelectionHists.end(); ++it) {
       it->second->Reset();
     }
 
     fTrueIncidentHist.Reset();
 
     fIncidentEnergies.clear();
     fESliceEnergies.clear();
     fFactor = 1.;
   };

void protoana::ThinSliceSample::ResetFactor ( )

inline

Definition at line 339 of file ThinSliceSample.h.

                      {
     ScaleHists(1./fFactor);
     fVariedFlux *= (1./fFactor);
     fFactor = 1.;
   };

void protoana::ThinSliceSample::SaveSystematics	(	std::string	syst_name,
		TDirectory *	dir
	)

inline

Definition at line 70 of file ThinSliceSample.h.

                                                               {
     dir->cd();
     //for (auto it = fSystematicShifts[syst_name].begin();
     //     it != fSystematicShifts[syst_name].end(); ++it) {
     //  for (size_t i = 0; i < it->second.size(); ++i) {
     //    it->second[i]->Write();
     //  }
     //} 
     for (auto it = fSystematicSplines[syst_name].begin();
          it != fSystematicSplines[syst_name].end(); ++it) {
       for (size_t i = 0; i < it->second.size(); ++i) {
         TCanvas c(it->second[i]->GetName(), "");
         it->second[i]->Draw();
         std::string name = syst_name + std::to_string(i);
         c.Write(syst_name.c_str());
       }
     }
   };

void protoana::ThinSliceSample::ScaleESliceEnergies ( double val )

inline

Definition at line 285 of file ThinSliceSample.h.

                                        {
     for (auto it = fESliceEnergies.begin();
          it != fESliceEnergies.end(); ++it) {
       it->second *= val;
     }
   };

void protoana::ThinSliceSample::ScaleHists ( double val )

inline

Definition at line 269 of file ThinSliceSample.h.

                               {
     for (auto it = fSelectionHists.begin(); it != fSelectionHists.end(); ++it) {
       it->second->Scale(val);
     }
     RefillRebinnedHists();
 
     fTrueIncidentHist.Scale(val);
   };

void protoana::ThinSliceSample::ScaleIncidentEnergies ( double val )

inline

Definition at line 278 of file ThinSliceSample.h.

                                          {
     for (auto it = fIncidentEnergies.begin();
          it != fIncidentEnergies.end(); ++it) {
       it->second *= val;
     }
   };

void protoana::ThinSliceSample::ScaleToDataMC ( )

inline

Definition at line 309 of file ThinSliceSample.h.

                        {
     ScaleHists(fDataMCScale);
     ScaleIncidentEnergies(fDataMCScale);
     ScaleESliceEnergies(fDataMCScale);
     ScaleVariedFlux(fDataMCScale);
   }

void protoana::ThinSliceSample::ScaleVariedFlux ( double val )

inline

Definition at line 305 of file ThinSliceSample.h.

                                    {
     fVariedFlux *= val;
   };

void protoana::ThinSliceSample::SetBestFit ( )

inline

Definition at line 353 of file ThinSliceSample.h.

                     {
     if (fBestFitIsSet) {
       return;
     }
 
     fBestFitFactor = fFactor;
     fBestFitIsSet = true;
   };

void protoana::ThinSliceSample::SetDataMCScale ( double val )

inline

Definition at line 316 of file ThinSliceSample.h.

                                   {
     fDataMCScale = val;
     ScaleHists(fDataMCScale);
     ScaleIncidentEnergies(fDataMCScale);
     fNominalFlux *= val;
     fVariedFlux *= val;
   };

void protoana::ThinSliceSample::SetFactor ( double f )

inline

Definition at line 34 of file ThinSliceSample.h.

34 {fFactor = f;};

keras_to_tensorflow.f

f

Definition: keras_to_tensorflow.py:162

protoana::ThinSliceSample::fFactor

double fFactor

Definition: ThinSliceSample.h:374

void protoana::ThinSliceSample::SetFactorAndScale ( double val )

inline

Definition at line 324 of file ThinSliceSample.h.

                                      {
     ResetFactor();
     fFactor = val;
     fVariedFlux *= val;
     ScaleHists(val);
   };

void protoana::ThinSliceSample::SetFactorToBestFit ( )

inline

Definition at line 345 of file ThinSliceSample.h.

                             {
     SetFactorAndScale(fBestFitFactor);
   };

void protoana::ThinSliceSample::SetSystematicSplines ( const std::map< std::string, std::map< int, std::vector< TSpline3 * >>> & input )

inline

Definition at line 55 of file ThinSliceSample.h.

                                                                                                  {
     fSystematicSplines = input;
   };

void protoana::ThinSliceSample::SetSystematicVals	(	std::string	syst_name,
		std::vector< double > &	vals
	)

inline

Definition at line 169 of file ThinSliceSample.h.

                                                                         {
     fSystematicVals[syst_name] = vals;
   };

Member Data Documentation

double protoana::ThinSliceSample::fBestFitFactor = 1.

private

Definition at line 374 of file ThinSliceSample.h.

bool protoana::ThinSliceSample::fBestFitIsSet = false

private

Definition at line 375 of file ThinSliceSample.h.

double protoana::ThinSliceSample::fDataMCScale = 1.

private

Definition at line 380 of file ThinSliceSample.h.

std::vector<std::pair<std::pair<double, double>, double> > protoana::ThinSliceSample::fESliceEnergies

private

Definition at line 408 of file ThinSliceSample.h.

double protoana::ThinSliceSample::fFactor = 1.

private

Definition at line 374 of file ThinSliceSample.h.

int protoana::ThinSliceSample::fFluxType

private

Definition at line 377 of file ThinSliceSample.h.

std::vector<std::pair<double, double> > protoana::ThinSliceSample::fIncidentEnergies

private

Definition at line 407 of file ThinSliceSample.h.

bool protoana::ThinSliceSample::fIsSignal

private

Definition at line 381 of file ThinSliceSample.h.

bool protoana::ThinSliceSample::fMadeRebinned = false

private

Definition at line 393 of file ThinSliceSample.h.

double protoana::ThinSliceSample::fNominalFlux = 0.

private

Definition at line 378 of file ThinSliceSample.h.

std::pair<double, double> protoana::ThinSliceSample::fRange

private

Definition at line 382 of file ThinSliceSample.h.

std::string protoana::ThinSliceSample::fSampleName

private

Definition at line 376 of file ThinSliceSample.h.

std::map<int, TH1 *> protoana::ThinSliceSample::fSelectionHists

private

Definition at line 387 of file ThinSliceSample.h.

std::map<int, TH1 *> protoana::ThinSliceSample::fSelectionHistsRebinned

private

Definition at line 390 of file ThinSliceSample.h.

std::map<std::string, std::map<int, std::vector<TH1 *> > > protoana::ThinSliceSample::fSystematicShifts

private

Definition at line 403 of file ThinSliceSample.h.

std::map<std::string, std::map<int, std::vector<TSpline3 *> > > protoana::ThinSliceSample::fSystematicSplines

private

Definition at line 401 of file ThinSliceSample.h.

std::map<std::string, std::vector<double> > protoana::ThinSliceSample::fSystematicVals

private

Definition at line 404 of file ThinSliceSample.h.

TH1D protoana::ThinSliceSample::fTrueIncidentHist

private

Definition at line 389 of file ThinSliceSample.h.

double protoana::ThinSliceSample::fVariedFlux = 0.

private

Definition at line 379 of file ThinSliceSample.h.

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

protoduneana/protoduneana/Utilities/FitUtils/ThinSliceSample.h
protoduneana/protoduneana/Utilities/FitUtils/ThinSliceSample.cxx

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation