Public Member Functions | Static Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
ana::ProportionalDecomp Class Reference

Splits Data proportionally according to MC. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-05-06/CAFAna/Decomp/ProportionalDecomp.h"

Inheritance diagram for ana::ProportionalDecomp:
ana::IDecomp

Public Member Functions

 ProportionalDecomp (SpectrumLoaderBase &loaderMC, SpectrumLoaderBase &loaderData, const std::string &label, const Binning &bins, const Var &var, const Cut &cut, const SystShifts &shiftMC=kNoShift, const SystShifts &shiftData=kNoShift, const Var &wei=kUnweighted)
 
 ProportionalDecomp (SpectrumLoaderBase &loaderMC, SpectrumLoaderBase &loaderData, const HistAxis &axis, const Cut &cut, const SystShifts &shiftMC=kNoShift, const SystShifts &shiftData=kNoShift, const Var &wei=kUnweighted)
 
 ProportionalDecomp (Loaders &loaders, const HistAxis &axis, const Cut &cut, const SystShifts &shiftMC=kNoShift, const SystShifts &shiftData=kNoShift, const Var &wei=kUnweighted)
 
Spectrum NumuComponent () const override
 
Spectrum AntiNumuComponent () const override
 
Spectrum NCTotalComponent () const override
 
Spectrum NCComponent () const override
 
Spectrum NCAntiComponent () const override
 
Spectrum NueComponent () const override
 
Spectrum AntiNueComponent () const override
 
Spectrum MC_NumuComponent () const override
 
Spectrum MC_AntiNumuComponent () const override
 
Spectrum MC_NueComponent () const override
 
Spectrum MC_AntiNueComponent () const override
 
Spectrum MC_NCTotalComponent () const override
 
Spectrum MC_NCComponent () const override
 
Spectrum MC_NCAntiComponent () const override
 
Spectrum Data_Component () const override
 
void SaveTo (TDirectory *dir, const std::string &name) const override
 

Static Public Member Functions

static std::unique_ptr< ProportionalDecomp > LoadFrom (TDirectory *dir, const std::string &name)
 

Protected Member Functions

 ProportionalDecomp ()
 

Protected Attributes

Spectrum fData
 
Spectrum fNCTot
 
Spectrum fNC
 
Spectrum fNCAnti
 
Spectrum fNue
 
Spectrum fAntiNue
 
Spectrum fNumu
 
Spectrum fAntiNumu
 
Spectrum fTotal
 

Detailed Description

Splits Data proportionally according to MC.

Definition at line 20 of file ProportionalDecomp.h.

Constructor & Destructor Documentation

ana::ProportionalDecomp::ProportionalDecomp ( SpectrumLoaderBase &  loaderMC,
SpectrumLoaderBase &  loaderData,
const std::string &  label,
const Binning &  bins,
const Var &  var,
const Cut &  cut,
const SystShifts &  shiftMC = kNoShift,
const SystShifts &  shiftData = kNoShift,
const Var &  wei = kUnweighted 
)

Definition at line 17 of file ProportionalDecomp.cxx.

26  : ProportionalDecomp(loaderMC,
27  loaderData,
28  HistAxis(label, bins, var),
29  cut,
30  shiftMC,
31  shiftData,
32  wei)
33  {}
stan::math::var
var_value< double > var
Definition: StanTypedefs.h:14
label
const char * label
Definition: cana_time_peak.C:34
ana::bins
const Binning bins
Definition: NumuCC_CPiBin.h:8
cut
const Cut cut
Definition: exporter_fd.C:30
ana::HistAxis
_HistAxis< Var > HistAxis
Definition: HistAxis.h:10
ana::ProportionalDecomp::ProportionalDecomp ( SpectrumLoaderBase &  loaderMC,
SpectrumLoaderBase &  loaderData,
const HistAxis &  axis,
const Cut &  cut,
const SystShifts &  shiftMC = kNoShift,
const SystShifts &  shiftData = kNoShift,
const Var &  wei = kUnweighted 
)

Definition at line 36 of file ProportionalDecomp.cxx.

43  : fData (loaderData, axis, cut, shiftData, wei),
44  // nc
45  fNCTot (loaderMC, axis, cut && kIsNC, shiftMC, wei),
46  fNC (loaderMC, axis, cut && kIsNC && !kIsAntiNu, shiftMC, wei),
47  fNCAnti (loaderMC, axis, cut && kIsNC && kIsAntiNu, shiftMC, wei),
48  //end nc
49  fNue (loaderMC, axis, cut && kIsBeamNue&&!kIsAntiNu, shiftMC, wei),
50  fAntiNue (loaderMC, axis, cut && kIsBeamNue&& kIsAntiNu, shiftMC, wei),
51  fNumu (loaderMC, axis, cut && kIsNumuCC &&!kIsAntiNu, shiftMC, wei),
52  fAntiNumu(loaderMC, axis, cut && kIsNumuCC && kIsAntiNu, shiftMC, wei),
53  fTotal (loaderMC, axis, cut && kHasNu, shiftMC, wei)
54  {}
ana::kIsNumuCC
const Cut kIsNumuCC(CCFlavSel(14, 14))
Select CC .
ana::kIsAntiNu
const Cut kIsAntiNu([](const caf::SRProxy *sr){if(sr->mc.nnu==0) return false;assert(sr->mc.nnu==1);return sr->mc.nu[0].pdg< 0;})
Is this truly an antineutrino?
Definition: TruthCuts.h:53
ana::kIsBeamNue
const Cut kIsBeamNue(CCFlavSel(12, 12))
Select CC .
ana::kIsNC
const Cut kIsNC([](const caf::SRProxy *sr){if(sr->mc.nnu==0) return false;assert(sr->mc.nnu==1);return!sr->mc.nu[0].iscc;})
Is this a Neutral Current event?
Definition: TruthCuts.h:8
cut
const Cut cut
Definition: exporter_fd.C:30
ana::kHasNu
const Cut kHasNu([](const caf::SRProxy *sr){if(sr->mc.nnu==0) return false;assert(sr->mc.nnu==1);return true;})
Definition: TruthCuts.h:56
ana::ProportionalDecomp::ProportionalDecomp ( Loaders &  loaders,
const HistAxis &  axis,
const Cut &  cut,
const SystShifts &  shiftMC = kNoShift,
const SystShifts &  shiftData = kNoShift,
const Var &  wei = kUnweighted 
)

Definition at line 57 of file ProportionalDecomp.cxx.

63  : ProportionalDecomp(loaders.GetLoader(caf::kNEARDET, Loaders::kMC),
64  loaders.GetLoader(caf::kNEARDET, Loaders::kData),
65  axis,
66  cut,
67  shiftMC,
68  shiftData,
69  wei)
70  {}
caf::kNEARDET
Near Detector underground.
Definition: SREnums.h:10
cut
const Cut cut
Definition: exporter_fd.C:30
loaders
std::vector< Loaders * > loaders
Definition: syst_header.h:386
ana::ProportionalDecomp::ProportionalDecomp ( )
inlineprotected

Definition at line 72 of file ProportionalDecomp.h.

73  : fData(Spectrum::Uninitialized()),
74  //nc
75  fNCTot (Spectrum::Uninitialized()),
76  fNC (Spectrum::Uninitialized()),
77  fNCAnti (Spectrum::Uninitialized()),
78  //end nc
79  fNue (Spectrum::Uninitialized()),
80  fAntiNue(Spectrum::Uninitialized()),
81  fNumu(Spectrum::Uninitialized()),
82  fAntiNumu(Spectrum::Uninitialized()),
83  fTotal(Spectrum::Uninitialized())
84  {};
ana::Spectrum::Uninitialized
static Spectrum Uninitialized()
Definition: Spectrum.h:144

Member Function Documentation

Spectrum ana::ProportionalDecomp::AntiNueComponent ( ) const
overridevirtual

Implements ana::IDecomp.

Definition at line 108 of file ProportionalDecomp.cxx.

References fAntiNue, fData, and fTotal.

Referenced by nue_decomp_scales_for_make_decomp().

109  {
110  return (fAntiNue/fTotal)*fData;
111  }
Spectrum ana::ProportionalDecomp::AntiNumuComponent ( ) const
overridevirtual

Implements ana::IDecomp.

Definition at line 79 of file ProportionalDecomp.cxx.

References fAntiNumu, fData, and fTotal.

Referenced by nue_decomp_scales_for_make_decomp().

80  {
81  return (fAntiNumu/fTotal)*fData;
82  }
Spectrum ana::ProportionalDecomp::Data_Component ( ) const
inlineoverridevirtual

Reimplemented from ana::IDecomp.

Definition at line 65 of file ProportionalDecomp.h.

References dir, fData, LoadFrom(), SaveTo(), and string.

Referenced by nue_decomp_scales_for_make_decomp().

65 {return fData;}
std::unique_ptr< ProportionalDecomp > ana::ProportionalDecomp::LoadFrom ( TDirectory *  dir,
const std::string &  name 
)
static

Definition at line 143 of file ProportionalDecomp.cxx.

References ana::assert(), dir, ana::Spectrum::LoadFrom(), and runNovaSAM::ret.

Referenced by BlessedPlotsAna(), BlessedPlotsAnaByPeriod(), CutTableAna(), Data_Component(), plot_nominal_spectra(), and PlotNus17Prediction().

144  {
145  dir = dir->GetDirectory(name.c_str()); // switch to subdir
146  assert(dir);
147 
148  std::unique_ptr<ProportionalDecomp> ret(new ProportionalDecomp);
149 
150  // nc
151  ret->fNCTot = *Spectrum::LoadFrom(dir, "nc_tot_comp");
152  ret->fNC = *Spectrum::LoadFrom(dir, "nc_comp");
153  ret->fNCAnti = *Spectrum::LoadFrom(dir, "nc_anti_comp");
154  //end nc
155  ret->fData = *Spectrum::LoadFrom(dir, "data_comp");
156  ret->fNue = *Spectrum::LoadFrom(dir, "nue_comp");
157  ret->fAntiNue = *Spectrum::LoadFrom(dir, "antinue_comp");
158  ret->fNumu = *Spectrum::LoadFrom(dir, "numu_comp");
159  ret->fAntiNumu = *Spectrum::LoadFrom(dir, "antinumu_comp");
160  ret->fTotal = *Spectrum::LoadFrom(dir, "total_comp");
161 
162  delete dir;
163 
164  return ret;
165  }
name
const XML_Char * name
Definition: expat.h:151
ana::Spectrum::LoadFrom
static std::unique_ptr< Spectrum > LoadFrom(TDirectory *dir, const std::string &name)
Definition: Spectrum.cxx:563
dir
TDirectory * dir
Definition: macro.C:5
ana::assert
assert(nhit_max >=nhit_nbins)
Spectrum ana::ProportionalDecomp::MC_AntiNueComponent ( ) const
inlineoverridevirtual

Reimplemented from ana::IDecomp.

Definition at line 60 of file ProportionalDecomp.h.

References fAntiNue.

Referenced by nue_decomp_scales_for_make_decomp().

60 {return fAntiNue;}
Spectrum ana::ProportionalDecomp::MC_AntiNumuComponent ( ) const
inlineoverridevirtual

Reimplemented from ana::IDecomp.

Definition at line 58 of file ProportionalDecomp.h.

References fAntiNumu.

Referenced by nue_decomp_scales_for_make_decomp().

58 {return fAntiNumu;}
Spectrum ana::ProportionalDecomp::MC_NCAntiComponent ( ) const
inlineoverridevirtual

Reimplemented from ana::IDecomp.

Definition at line 63 of file ProportionalDecomp.h.

References fNCAnti.

63 {return fNCAnti;}
Spectrum ana::ProportionalDecomp::MC_NCComponent ( ) const
inlineoverridevirtual

Reimplemented from ana::IDecomp.

Definition at line 62 of file ProportionalDecomp.h.

References fNC.

Referenced by nue_decomp_scales_for_make_decomp().

62 {return fNC;}
Spectrum ana::ProportionalDecomp::MC_NCTotalComponent ( ) const
inlineoverridevirtual

Reimplemented from ana::IDecomp.

Definition at line 61 of file ProportionalDecomp.h.

References fNCTot.

61 {return fNCTot;}
Spectrum ana::ProportionalDecomp::MC_NueComponent ( ) const
inlineoverridevirtual

Reimplemented from ana::IDecomp.

Definition at line 59 of file ProportionalDecomp.h.

References fNue.

Referenced by nue_decomp_scales_for_make_decomp().

59 {return fNue;}
Spectrum ana::ProportionalDecomp::MC_NumuComponent ( ) const
inlineoverridevirtual

Reimplemented from ana::IDecomp.

Definition at line 57 of file ProportionalDecomp.h.

References fNumu.

Referenced by nue_decomp_scales_for_make_decomp().

57 {return fNumu;}
Spectrum ana::ProportionalDecomp::NCAntiComponent ( ) const
overridevirtual

Implements ana::IDecomp.

Definition at line 96 of file ProportionalDecomp.cxx.

References fData, fNCAnti, and fTotal.

97  {
98  return (fNCAnti/fTotal)*fData;
99  }
Spectrum ana::ProportionalDecomp::NCComponent ( ) const
overridevirtual

Implements ana::IDecomp.

Definition at line 91 of file ProportionalDecomp.cxx.

References fData, fNC, and fTotal.

Referenced by nue_decomp_scales_for_make_decomp().

92  {
93  return (fNC/fTotal)*fData;
94  }
Spectrum ana::ProportionalDecomp::NCTotalComponent ( ) const
overridevirtual

Reimplemented from ana::IDecomp.

Definition at line 86 of file ProportionalDecomp.cxx.

References fData, fNCTot, and fTotal.

87  {
88  return (fNCTot/fTotal)*fData;
89  }
Spectrum ana::ProportionalDecomp::NueComponent ( ) const
overridevirtual

Implements ana::IDecomp.

Definition at line 102 of file ProportionalDecomp.cxx.

References fData, fNue, and fTotal.

Referenced by nue_decomp_scales_for_make_decomp().

103  {
104  return (fNue/fTotal)*fData;
105  }
Spectrum ana::ProportionalDecomp::NumuComponent ( ) const
overridevirtual

Implements ana::IDecomp.

Definition at line 73 of file ProportionalDecomp.cxx.

References fData, fNumu, and fTotal.

Referenced by nue_decomp_scales_for_make_decomp().

74  {
75  return (fNumu/fTotal)*fData;
76  }
void ana::ProportionalDecomp::SaveTo ( TDirectory *  dir,
const std::string &  name 
) const
overridevirtual

Implements ana::IDecomp.

Definition at line 114 of file ProportionalDecomp.cxx.

References dir, fAntiNue, fAntiNumu, fData, fNC, fNCAnti, fNCTot, fNue, fNumu, fTotal, ana::Spectrum::SaveTo(), and tmp.

Referenced by BlessedPlotsLoad(), CutTableLoad(), Data_Component(), make_decomp(), make_michel(), MakeNus17Prediction(), MakeNusPrediction(), and nue_decomp_scales().

115  {
116  TDirectory* tmp = gDirectory;
117 
118  dir = dir->mkdir(name.c_str()); // switch to subdir
119  dir->cd();
120 
121  TObjString("ProportionalDecomp").Write("type");
122 
123  fData.SaveTo(dir, "data_comp");
124  //nc
125  fNCTot.SaveTo(dir, "nc_tot_comp");
126  fNC.SaveTo(dir, "nc_comp");
127  fNCAnti.SaveTo(dir, "nc_anti_comp");
128  //end nc
129  fNue.SaveTo(dir, "nue_comp");
130  fAntiNue.SaveTo(dir, "antinue_comp");
131  fNumu.SaveTo(dir, "numu_comp");
132  fAntiNumu.SaveTo(dir, "antinumu_comp");
133  fTotal.SaveTo(dir, "total_comp");
134 
135  dir->Write();
136  delete dir;
137 
138  tmp->cd();
139  }
name
const XML_Char * name
Definition: expat.h:151
tmp
Float_t tmp
Definition: plot.C:36
ana::Spectrum::SaveTo
void SaveTo(TDirectory *dir, const std::string &name) const
Definition: Spectrum.cxx:534
dir
TDirectory * dir
Definition: macro.C:5

Member Data Documentation

Spectrum ana::ProportionalDecomp::fAntiNue
protected

Definition at line 93 of file ProportionalDecomp.h.

Referenced by AntiNueComponent(), MC_AntiNueComponent(), and SaveTo().

Spectrum ana::ProportionalDecomp::fAntiNumu
protected

Definition at line 95 of file ProportionalDecomp.h.

Referenced by AntiNumuComponent(), MC_AntiNumuComponent(), and SaveTo().

Spectrum ana::ProportionalDecomp::fData
protected

Definition at line 84 of file ProportionalDecomp.h.

Referenced by AntiNueComponent(), AntiNumuComponent(), Data_Component(), NCAntiComponent(), NCComponent(), NCTotalComponent(), NueComponent(), NumuComponent(), and SaveTo().

Spectrum ana::ProportionalDecomp::fNC
protected

Definition at line 89 of file ProportionalDecomp.h.

Referenced by MC_NCComponent(), NCComponent(), and SaveTo().

Spectrum ana::ProportionalDecomp::fNCAnti
protected

Definition at line 90 of file ProportionalDecomp.h.

Referenced by MC_NCAntiComponent(), NCAntiComponent(), and SaveTo().

Spectrum ana::ProportionalDecomp::fNCTot
protected

Definition at line 88 of file ProportionalDecomp.h.

Referenced by MC_NCTotalComponent(), NCTotalComponent(), and SaveTo().

Spectrum ana::ProportionalDecomp::fNue
protected

Definition at line 92 of file ProportionalDecomp.h.

Referenced by MC_NueComponent(), NueComponent(), and SaveTo().

Spectrum ana::ProportionalDecomp::fNumu
protected

Definition at line 94 of file ProportionalDecomp.h.

Referenced by MC_NumuComponent(), NumuComponent(), and SaveTo().

Spectrum ana::ProportionalDecomp::fTotal
protected

Definition at line 96 of file ProportionalDecomp.h.

Referenced by AntiNueComponent(), AntiNumuComponent(), NCAntiComponent(), NCComponent(), NCTotalComponent(), NueComponent(), NumuComponent(), and SaveTo().

The documentation for this class was generated from the following files:
Generated on Fri May 7 2021 14:13:37 for NOvA by   doxygen 1.8.11