Incorporate into a weight Var like Var wei(BeamWeight(...
Definition at line 55 of file BeamSysts.h.
Definition at line 350 of file BeamSysts.cxx.
References om::cerr, fillBadChanDBTables::det, caf::Proxy< caf::SRNeutrino >::det, caf::Proxy< caf::SRNeutrino >::E, allTimeWatchdog::endl, energy, ana::BeamSystOrWeightBase::fHistName, ana::BeamSystOrWeightBase::fHistos, flux, make_syst_table_plots::h, ana::BeamSystOrWeightBase::InitializeHistograms(), caf::Proxy< caf::SRNeutrino >::isRHC, ana::BeamSystOrWeightBase::kAntiNue, ana::BeamSystOrWeightBase::kAntiNumu, caf::kFARDET, ana::BeamSystOrWeightBase::kFD, ana::BeamSystOrWeightBase::kFHC, ana::BeamSystOrWeightBase::kND, caf::kNEARDET, ana::BeamSystOrWeightBase::kNue, ana::BeamSystOrWeightBase::kNumu, ana::BeamSystOrWeightBase::kRHC, caf::Proxy< caf::SRNeutrino >::pdgorig, and canMan::sign().
Referenced by operator()(), and ana::BeamSyst::TruthShift().
375 case 14: flav =
kNumu;
break;
376 case 12: flav =
kNue;
break;
389 if (energy > h->GetXaxis()->GetXmin() &&
390 energy < h->GetXaxis()->GetXmax() ){
391 return h->Interpolate(energy);
397 if (energy > h->GetXaxis()->GetXmax())
398 return h->GetBinContent(h->GetNbinsX());
Near Detector underground.
Far Detector at Ash River.
caf::Proxy< caf::Det_t > det
caf::Proxy< short int > pdgorig
void InitializeHistograms() const
TH1D * fHistos[kNumFluxType][kNumDets][kNumSigns][kNumFlavors]
void ana::BeamSystOrWeightBase::InitializeHistograms |
( |
| ) |
const |
|
protectedinherited |
Definition at line 281 of file BeamSysts.cxx.
References om::cerr, fillBadChanDBTables::det, allTimeWatchdog::endl, ana::BeamSystOrWeightBase::fFileName, ana::BeamSystOrWeightBase::fHistName, ana::BeamSystOrWeightBase::fHistos, fin, flux, compareCafs::histName, ana::BeamSystOrWeightBase::kAntiNue, ana::BeamSystOrWeightBase::kAntiNumu, ana::BeamSystOrWeightBase::kFD, ana::BeamSystOrWeightBase::kFHC, ana::BeamSystOrWeightBase::kMinus, ana::BeamSystOrWeightBase::kND, ana::BeamSystOrWeightBase::kNue, ana::BeamSystOrWeightBase::kNumu, ana::BeamSystOrWeightBase::kPlus, ana::BeamSystOrWeightBase::kRHC, and canMan::sign().
Referenced by ana::BeamSystOrWeightBase::GetWeight().
284 if(
fHistos[0][0][0][0])
return;
296 TIter iterHist(gDirectory->GetListOfKeys());
299 bool isSeparatedByFlavor =
false;
301 while((keyHist = (TKey*)iterHist())) {
302 TString
histName = keyHist->GetName();
307 if(histName.Contains(
"FHC")) flux =
kFHC;
308 if(histName.Contains(
"RHC")) flux =
kRHC;
311 if(histName.Contains(
"ND")) det =
kND;
312 if(histName.Contains(
"FD")) det =
kFD;
315 if(histName.Contains(
"minus")) sign =
kMinus;
316 if(histName.Contains(
"plus" )) sign =
kPlus;
319 if(histName.Contains(
"numu")) flav =
kNumu;
320 if(histName.Contains(
"nue")) flav =
kNue;
321 if(histName.Contains(
"anumu")) flav =
kAntiNumu;
322 if(histName.Contains(
"anue")) flav =
kAntiNue;
329 if(histName.Contains(
"numu")){
330 isSeparatedByFlavor=
true;
340 if (!isSeparatedByFlavor) {
342 <<
" doesn't have required flavor information; aborting"
TH1D * fHistos[kNumFluxType][kNumDets][kNumSigns][kNumFlavors]