Public Types | Public Member Functions | Static Public Member Functions | Protected Types | Protected Member Functions | Protected Attributes | List of all members
lem::MakePID Class Reference

Calculate final PID variables from match details. More...

Inheritance diagram for lem::MakePID:
art::EDProducer art::ProducerBase art::Consumer art::EngineCreator art::ProductRegistryHelper

Public Types

using ModuleType = EDProducer
 
using WorkerType = WorkerT< EDProducer >
 
template<typename UserConfig , typename KeysToIgnore = void>
using Table = ProducerBase::Table< UserConfig, KeysToIgnore >
 

Public Member Functions

 MakePID (const fhicl::ParameterSet &pset)
 
 ~MakePID ()
 
virtual void beginJob ()
 
virtual void beginRun (art::Run &run)
 
virtual void produce (art::Event &evt)
 
virtual void reconfigure (const fhicl::ParameterSet &pset)
 
template<typename PROD , BranchType B = InEvent>
ProductID getProductID (std::string const &instanceName={}) const
 
template<typename PROD , BranchType B>
ProductID getProductID (ModuleDescription const &moduleDescription, std::string const &instanceName) const
 
bool modifiesEvent () const
 
template<typename T , BranchType = InEvent>
ProductToken< T > consumes (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< T > consumes (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< T > consumesView (InputTag const &it)
 
template<typename T , BranchType = InEvent>
ProductToken< T > mayConsume (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ProductToken< T > mayConsume (InputTag const &it)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , art::BranchType BT>
art::ViewToken< T > mayConsumeView (InputTag const &it)
 
base_engine_t & createEngine (seed_t seed)
 
base_engine_t & createEngine (seed_t seed, std::string const &kind_of_engine_to_make)
 
base_engine_t & createEngine (seed_t seed, std::string const &kind_of_engine_to_make, label_t const &engine_label)
 
seed_t get_seed_value (fhicl::ParameterSet const &pset, char const key[]="seed", seed_t const implicit_seed=-1)
 

Static Public Member Functions

static cet::exempt_ptr< Consumer > non_module_context ()
 

Protected Types

enum  WeightScheme {
  kUnweighted, kSimple, kExp, kSigOnly,
  kSigOnlyExp, kBkgOnly, kBkgOnlyExp, kSig50Only,
  kSig90Only, kPi0Only
}
 
enum  MatchVar {
  kSigFrac, kBkgFrac, kNumuFrac, kY,
  kQFrac, kSigFrac50, kSigFrac90, kL0,
  kL1, kTrueEVis, kQESigFrac, kResSigFrac,
  kDISSigFrac, kQENCFrac, kResNCFrac, kDISNCFrac,
  kQENumuFrac, kResNumuFrac, kDISNumuFrac
}
 

Protected Member Functions

void AvgAndFit (const std::vector< MatchSummary > &matches, WeightScheme weightScheme, MatchVar matchVar, double &mean, double &fit) const
 
CurrentProcessingContext const * currentContext () const
 
void validateConsumedProduct (BranchType const bt, ProductInfo const &pi)
 
void prepareForJob (fhicl::ParameterSet const &pset)
 
void showMissingConsumes () const
 

Protected Attributes

bool fScaleCalE
 
std::vector< double > fExpTable
 
double fNCScale
 
double fSigScale
 
double fEnrichScale
 
double fExpLength
 
double fExpPow
 
bool fTrimMatchLists
 
bool fUseFitInput
 
bool fPositionalWeight
 
std::string fMatchesLabel
 
std::string fSummarizeLabel
 
TMVA::Reader * fReader
 
dec::Forest * fForest
 
TH1 * fFlatHist
 

Detailed Description

Calculate final PID variables from match details.

Definition at line 48 of file MakePID_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer
inherited

Definition at line 34 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>
using art::EDProducer::Table = ProducerBase::Table<UserConfig, KeysToIgnore>
inherited

Definition at line 43 of file EDProducer.h.

using art::EDProducer::WorkerType = WorkerT<EDProducer>
inherited

Definition at line 35 of file EDProducer.h.

Member Enumeration Documentation

enum lem::MakePID::MatchVar
protected
Enumerator
kSigFrac 
kBkgFrac 
kNumuFrac 
kY 
kQFrac 
kSigFrac50 
kSigFrac90 
kL0 
kL1 
kTrueEVis 
kQESigFrac 
kResSigFrac 
kDISSigFrac 
kQENCFrac 
kResNCFrac 
kDISNCFrac 
kQENumuFrac 
kResNumuFrac 
kDISNumuFrac 

Definition at line 75 of file MakePID_module.cc.

75  {
76  kSigFrac,
77  kBkgFrac,
78  kNumuFrac,
79  kY,
80  kQFrac,
81  kSigFrac50,
82  kSigFrac90,
83  kL0,
84  kL1,
85  kTrueEVis,
86  kQESigFrac,
87  kResSigFrac,
88  kDISSigFrac,
89  kQENCFrac,
90  kResNCFrac,
91  kDISNCFrac,
92  kQENumuFrac,
93  kResNumuFrac,
94  kDISNumuFrac
95  };
enum lem::MakePID::WeightScheme
protected
Enumerator
kUnweighted 
kSimple 
kExp 
kSigOnly 
kSigOnlyExp 
kBkgOnly 
kBkgOnlyExp 
kSig50Only 
kSig90Only 
kPi0Only 

Definition at line 62 of file MakePID_module.cc.

62  {
63  kUnweighted,
64  kSimple,
65  kExp,
66  kSigOnly,
67  kSigOnlyExp,
68  kBkgOnly,
69  kBkgOnlyExp,
70  kSig50Only,
71  kSig90Only,
72  kPi0Only
73  };

Constructor & Destructor Documentation

lem::MakePID::MakePID ( const fhicl::ParameterSet &  pset)
explicit

Definition at line 124 of file MakePID_module.cc.

References ana::assert(), om::cout, allTimeWatchdog::endl, util::EnvExpansion(), fFlatHist, fForest, fReader, lem::dec::Forest::FromFile(), fhicl::ParameterSet::get(), gJunk, gTMVAVar, path, reconfigure(), and string.

125  {
126  reconfigure(pset);
127 
128  produces<std::vector<lem::PIDDetails>>();
129  produces<art::Assns<lem::PIDDetails, rb::Cluster>>();
130 
131  const std::string path = pset.get<std::string>("LibraryPath");
132  const std::string libPath = util::EnvExpansion(path) + "/";
133 
134  fReader = 0;
135  if(pset.get<bool>("FillTMVAs")){
136  fReader = new TMVA::Reader;
137 
138  fReader->AddVariable("pidExtrap", &gTMVAVar[0]); // vars.fPIDFit
139  fReader->AddVariable("yexp", &gTMVAVar[1]); // vars.fMeanYExp
140  fReader->AddVariable("qfracexp", &gTMVAVar[2]); // vars.fMeanQFracExp
141  fReader->AddVariable("ediff", &gTMVAVar[3]); // vars.fEnergyDiffExp
142  fReader->AddVariable("calE", &gTMVAVar[4]); // 2*totalRecoGeV
143  fReader->AddVariable("enrichExp", &gTMVAVar[5]); // vars.fEnrichFracExp
144 
145  fReader->AddSpectator("isSig", &gJunk);
146  fReader->AddSpectator("myW", &gJunk);
147 
148  const std::string tmvaPath = libPath+"weights/";
149 
150  fReader->BookMVA("BDTG", tmvaPath+"TMVAClassification_BDTG.weights.xml");
151  fReader->BookMVA("BDT", tmvaPath+"TMVAClassification_BDT.weights.xml");
152  fReader->BookMVA("BDTD", tmvaPath+"TMVAClassification_BDTD.weights.xml");
153  fReader->BookMVA("BDTMitFisher", tmvaPath+"TMVAClassification_BDTMitFisher.weights.xml");
154  fReader->BookMVA("MLP", tmvaPath+"TMVAClassification_MLP.weights.xml");
155 
156  fReader->BookMVA("BDTG_CC", tmvaPath+"TMVAClassification_BDTG_CC.weights.xml");
157  fReader->BookMVA("MLP_CC", tmvaPath+"TMVAClassification_MLP_CC.weights.xml");
158  }
159 
160  // fANN = fANNEnrich = 0;
161  // if(pset.get<bool>("FillANNs")){
162  // fANN = fann_create_from_file((libPath+"fann.net").c_str());
163  // assert(fANN);
164  // fann_print_connections(fANN);
165 
166  // fANNEnrich = fann_create_from_file((libPath+"fann_enrich.net").c_str());
167  // assert(fANNEnrich);
168  // fann_print_connections(fANNEnrich);
169  // }
170 
171  fForest = 0;
172  if(pset.get<bool>("FillDecTree")){
173  std::cout << "Loading LEM dectree from " << libPath << std::endl;
174  fForest = dec::Forest::FromFile(libPath+"dectree.root");
175  }
176 
177  fFlatHist = 0;
178  if(pset.get<bool>("Flatten")){
179  TFile* fFlat = new TFile((libPath+"flatten.root").c_str());
180  assert(!fFlat->IsZombie());
181  fFlatHist = (TH1*)fFlat->Get("conv");
182  assert(fFlatHist);
183  }
184  }
lem::MakePID::reconfigure
virtual void reconfigure(const fhicl::ParameterSet &pset)
Definition: MakePID_module.cc:194
lem::dec::Forest::FromFile
static Forest * FromFile(const std::string &fname)
Load PID from a file.
Definition: DecisionTree.cxx:59
util::EnvExpansion
std::string EnvExpansion(const std::string &inString)
Function to expand environment variables.
Definition: EnvExpand.cxx:8
lem::MakePID::fForest
dec::Forest * fForest
Definition: MakePID_module.cc:118
gTMVAVar
float gTMVAVar[6]
Definition: MakePID_module.cc:39
fhicl::ParameterSet::get
T get(std::string const &key) const
Definition: ParameterSet.h:231
om::cout
OStream cout
Definition: OStream.cxx:6
path
const std::string path
Definition: plot_BEN.C:43
ana::assert
assert(nhit_max >=nhit_nbins)
gJunk
float gJunk
Definition: MakePID_module.cc:40
lem::MakePID::fReader
TMVA::Reader * fReader
Definition: MakePID_module.cc:117
string
enum BeamMode string
lem::MakePID::~MakePID ( )

Definition at line 187 of file MakePID_module.cc.

References fForest, and fReader.

188  {
189  delete fReader;
190  delete fForest;
191  }
lem::MakePID::fForest
dec::Forest * fForest
Definition: MakePID_module.cc:118
lem::MakePID::fReader
TMVA::Reader * fReader
Definition: MakePID_module.cc:117

Member Function Documentation

void lem::MakePID::AvgAndFit ( const std::vector< MatchSummary > &  matches,
WeightScheme  weightScheme,
MatchVar  matchVar,
double &  mean,
double &  fit 
) const
protected

Definition at line 235 of file MakePID_module.cc.

References abs(), plot_validation_datamc::c, lem::MatchSummary::ccnc, stan::math::exp(), fExpLength, fExpPow, fExpTable, fPositionalWeight, MECModelEnuComparisons::i, makeTrainCVSamples::int, lem::MatchSummary::isSig, kBkgFrac, kBkgOnly, kBkgOnlyExp, simb::kCC, simb::kDIS, kDISNCFrac, kDISNumuFrac, kDISSigFrac, kExp, kL0, kL1, simb::kNC, kNumuFrac, kPi0Only, simb::kQE, kQENCFrac, kQENumuFrac, kQESigFrac, kQFrac, simb::kRes, kResNCFrac, kResNumuFrac, kResSigFrac, kSig50Only, kSig90Only, kSigFrac, kSigFrac50, kSigFrac90, kSigOnly, kSigOnlyExp, kTrueEVis, kUnweighted, kY, util::LinFit(), util::LinFitUnweighted(), m, cafExposure::match, lem::MatchSummary::mode, lem::MatchSummary::pdg, lem::MatchSummary::photE0, lem::MatchSummary::photE1, lem::MatchSummary::photL0, lem::MatchSummary::photL1, lem::MatchSummary::potential, cet::pow(), lem::MatchSummary::qFrac, PandAna.Demos.pi0_spectra::tot, lem::MatchSummary::trueEVis, PandAna.Demos.tute_pid_validation::var, w, lem::MatchSummary::weight, and lem::MatchSummary::y.

Referenced by produce().

238  {
239  std::vector<double> xs, ys, ws;
240  double tot = 0;
241  mean = 0;
242 
243  const double maxPotential = matches.empty() ? 1 : matches[matches.size()-1].potential;
244 
245  for(unsigned int i = 0; i < matches.size(); ++i){
246  const MatchSummary& match = matches[i];
247 
248  double w = match.weight;
249 
250  if(weightScheme == kUnweighted) w = 1;
251 
252  if(weightScheme == kExp ||
253  weightScheme == kSigOnlyExp ||
254  weightScheme == kBkgOnlyExp){
255  if(fPositionalWeight){
256  w *= fExpTable[i];
257  }
258  else{
259  w *= exp(-pow(match.potential/maxPotential, fExpPow)/fExpLength);
260  }
261  }
262 
263  if(weightScheme == kSigOnly || weightScheme == kSigOnlyExp){
264  if(!match.isSig) w = 0;
265  }
266 
267  if(weightScheme == kBkgOnly || weightScheme == kBkgOnlyExp){
268  if(match.isSig) w = 0;
269  }
270 
271  if(weightScheme == kSig50Only){
272  if(!match.isSig || match.y > 0.50) w = 0;
273  }
274  if(weightScheme == kSig90Only){
275  if(!match.isSig || match.y > 0.90) w = 0;
276  }
277  if(weightScheme == kPi0Only) w = (match.photL0 > 0 || match.photL1 > 0 || match.photE0 > 0 || match.photE1 > 0);
278 
279  double var = 0;
280  switch(matchVar){
281  case kSigFrac: var = int(match.isSig); break;
282  case kBkgFrac: var = 1-int(match.isSig); break;
283  case kNumuFrac: var = int(match.ccnc == simb::kCC && abs(match.pdg) == 14); break;
284  case kY: var = match.y; break;
285  case kQFrac: var = match.qFrac; break;
286  case kSigFrac50: var = int(match.isSig && match.y < .50); break;
287  case kSigFrac90: var = int(match.isSig && match.y < .90); break;
288  case kL0: var = match.photL0; break;
289  case kL1: var = match.photL1; break;
290  case kTrueEVis: var = match.trueEVis; break;
291  case kQESigFrac: var = (match.isSig && match.mode == simb::kQE); break;
292  case kResSigFrac: var = (match.isSig && match.mode == simb::kRes); break;
293  case kDISSigFrac: var = (match.isSig && match.mode == simb::kDIS); break;
294  case kQENCFrac: var = (match.ccnc == simb::kNC && match.mode == simb::kQE); break;
295  case kResNCFrac: var = (match.ccnc == simb::kNC && match.mode == simb::kRes); break;
296  case kDISNCFrac: var = (match.ccnc == simb::kNC && match.mode == simb::kDIS); break;
297  case kQENumuFrac: var = (match.ccnc == simb::kCC && abs(match.pdg) == 14 && match.mode == simb::kQE); break;
298  case kResNumuFrac: var = (match.ccnc == simb::kCC && abs(match.pdg) == 14 && match.mode == simb::kRes); break;
299  case kDISNumuFrac: var = (match.ccnc == simb::kCC && abs(match.pdg) == 14 && match.mode == simb::kDIS); break;
300  }
301 
302  xs.push_back(match.potential);
303  ys.push_back(var);
304  ws.push_back(w);
305 
306  tot += w;
307  mean += w*var;
308  } // end for i
309 
310  mean /= tot;
311 
312  if(ws.size() < 2){
313  fit = -1;
314  }
315  else{
316  double m, c;
317  if(weightScheme == kUnweighted)
318  util::LinFitUnweighted(xs, ys, m, c);
319  else
320  util::LinFit(xs, ys, ws, m, c);
321  fit = m*xs[0]+c;
322  }
323  }
PandAna.Demos.tute_pid_validation.var
list var
Definition: tute_pid_validation.py:17
PandAna.Demos.pi0_spectra.tot
tot
Definition: pi0_spectra.py:117
cet::pow
constexpr T pow(T x)
Definition: pow.h:75
abs
void abs(TH1 *hist)
Definition: wrong_sign_uncertainty.C:716
lem::MakePID::fExpTable
std::vector< double > fExpTable
Definition: MakePID_module.cc:101
overlay_prestage_def.matches
matches
Definition: overlay_prestage_def.py:135
m
#define m
Definition: MakeGENIEReweightTable_module.cc:407
stan::math::exp
fvar< T > exp(const fvar< T > &x)
Definition: exp.hpp:10
plot_validation_datamc.c
c
Definition: plot_validation_datamc.py:198
util::LinFitUnweighted
void LinFitUnweighted(const std::vector< double > &x, const std::vector< double > &y, double &m, double &c)
Simplified version of LinFit.
Definition: MathUtil.cxx:8
util::LinFit
double LinFit(const std::vector< double > &x, const std::vector< double > &y, const std::vector< double > &w, double &m, double &c)
Find the best-fit line to a collection of points in 2-D by minimizing the squared vertical distance f...
Definition: MathUtil.cxx:36
w
Float_t w
Definition: plot.C:20
void lem::MakePID::beginJob ( )
virtual

Reimplemented from art::EDProducer.

Definition at line 208 of file MakePID_module.cc.

References stan::math::exp(), fExpLength, fExpPow, fExpTable, MECModelEnuComparisons::i, lem::kMaxNumMatches, and cet::pow().

209  {
210  fExpTable.resize(kMaxNumMatches);
211  for(unsigned int i = 0; i < kMaxNumMatches; ++i){
212  fExpTable[i] = exp(-pow(i/fExpLength, fExpPow));
213  }
214  }
lem::kMaxNumMatches
const unsigned int kMaxNumMatches
Definition: FindMatches.h:29
cet::pow
constexpr T pow(T x)
Definition: pow.h:75
lem::MakePID::fExpTable
std::vector< double > fExpTable
Definition: MakePID_module.cc:101
stan::math::exp
fvar< T > exp(const fvar< T > &x)
Definition: exp.hpp:10
void lem::MakePID::beginRun ( art::Run &  run)
virtual

Reimplemented from art::EDProducer.

Definition at line 217 of file MakePID_module.cc.

References ana::assert(), om::cout, art::Handle< T >::failedToGet(), fEnrichScale, fMatchesLabel, fNCScale, fSigScale, art::DataViewImpl::getByLabel(), lem::LibrarySummary::nSwapNC, lem::LibrarySummary::nTrueNC, lem::LibrarySummary::totBkg, lem::LibrarySummary::totEnrich, and lem::LibrarySummary::totSig.

218  {
219  art::Handle<lem::LibrarySummary> libsum;
220  run.getByLabel(fMatchesLabel, libsum);
221 
222  assert(!libsum.failedToGet());
223 
224  fNCScale = libsum->nTrueNC/double(libsum->nTrueNC+libsum->nSwapNC);
225  fSigScale = double(libsum->totBkg-libsum->nSwapNC)/libsum->totSig;
226  fEnrichScale = double(libsum->totBkg-libsum->nSwapNC)/libsum->totEnrich;
227 
228  std::cout << "Reweighted to sig, bkg, enrich = "
229  << libsum->totSig*fSigScale << ", "
230  << libsum->totBkg+(fNCScale-1)*libsum->nTrueNC-libsum->nSwapNC << ", "
231  << libsum->totEnrich*fEnrichScale;
232  }
om::cout
OStream cout
Definition: OStream.cxx:6
lem::MakePID::fMatchesLabel
std::string fMatchesLabel
Definition: MakePID_module.cc:114
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
Definition: DataViewImpl.h:344
ana::assert
assert(nhit_max >=nhit_nbins)
art::Handle::failedToGet
bool failedToGet() const
Definition: Handle.h:196
template<typename T , BranchType = InEvent>
ProductToken<T> art::Consumer::consumes ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ProductToken<T> art::Consumer::consumes ( InputTag const &  it)
inherited

Definition at line 146 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

147 {
148  if (!moduleContext_)
149  return ProductToken<T>::invalid();
150 
151  consumables_[BT].emplace_back(ConsumableType::Product,
152  TypeID{typeid(T)},
153  it.label(),
154  it.instance(),
155  it.process());
156  return ProductToken<T>{it};
157 }
it
set< int >::iterator it
Definition: ProduceBadRunsByFailureMode.C:101
art::ProductToken::invalid
static ProductToken< T > invalid()
Definition: ProductToken.h:47
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:137
T
double T
Definition: Xdiff_gwt.C:5
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:135
template<typename T , art::BranchType BT>
void art::Consumer::consumesMany ( )
inherited

Definition at line 161 of file Consumer.h.

References T.

162 {
163  if (!moduleContext_)
164  return;
165 
166  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
167 }
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:137
T
double T
Definition: Xdiff_gwt.C:5
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:135
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::Consumer::consumesView ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ViewToken<T> art::Consumer::consumesView ( InputTag const &  it)
inherited

Definition at line 171 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

172 {
173  if (!moduleContext_)
174  return ViewToken<T>::invalid();
175 
176  consumables_[BT].emplace_back(ConsumableType::ViewElement,
177  TypeID{typeid(T)},
178  it.label(),
179  it.instance(),
180  it.process());
181  return ViewToken<T>{it};
182 }
it
set< int >::iterator it
Definition: ProduceBadRunsByFailureMode.C:101
art::ViewToken::invalid
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:137
T
double T
Definition: Xdiff_gwt.C:5
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:135
base_engine_t& art::EngineCreator::createEngine ( seed_t  seed)
inherited

Referenced by evgen::CosmicPionGen::CosmicPionGen(), evgen::CosmicsGen::CosmicsGen(), fuzz::FuzzyKVertex::FuzzyKVertex(), g4n::G4Gen::G4Gen(), photrans::ImprovedTransport::ImprovedTransport(), art::MixFilter< T >::initEngine_(), evgen::IsoGen::IsoGen(), rsim::MakeNoiseSpectrumFile::MakeNoiseSpectrumFile(), photrans::PhotonSystRewriter::PhotonSystRewriter(), photrans::PhotonTransport::PhotonTransport(), rsim::ReadoutSim::ReadoutSim(), evgen::SingleGen::SingleGen(), and evgen::SupernovaGen::SupernovaGen().

base_engine_t& art::EngineCreator::createEngine ( seed_t  seed,
std::string const &  kind_of_engine_to_make 
)
inherited
base_engine_t& art::EngineCreator::createEngine ( seed_t  seed,
std::string const &  kind_of_engine_to_make,
label_t const &  engine_label 
)
inherited
CurrentProcessingContext const* art::EDProducer::currentContext ( ) const
protectedinherited
seed_t art::EngineCreator::get_seed_value ( fhicl::ParameterSet const &  pset,
char const  key[] = "seed",
seed_t const  implicit_seed = -1 
)
inherited

Referenced by art::MixFilter< T >::initEngine_().

template<typename PROD , BranchType B>
ProductID art::EDProducer::getProductID ( std::string const &  instanceName = {}) const
inlineinherited

Definition at line 123 of file EDProducer.h.

References art::EDProducer::moduleDescription_.

Referenced by skim::NueSkimmer::CopyMichelSlice(), and skim::NueSkimmer::CopyMichelTrack().

124  {
125  return ProducerBase::getProductID<PROD, B>(moduleDescription_,
126  instanceName);
127  }
art::EDProducer::moduleDescription_
ModuleDescription moduleDescription_
Definition: EDProducer.h:115
template<typename PROD , BranchType B>
ProductID art::ProducerBase::getProductID ( ModuleDescription const &  moduleDescription,
std::string const &  instanceName 
) const
inherited

Definition at line 56 of file ProducerBase.h.

References art::ModuleDescription::moduleLabel().

Referenced by art::ProducerBase::modifiesEvent().

58  {
59  auto const& pd =
60  get_ProductDescription<PROD>(B, md.moduleLabel(), instanceName);
61  return pd.productID();
62  }
fetch_tb_beamline_files.md
md
Definition: fetch_tb_beamline_files.py:11
template<typename T , BranchType = InEvent>
ProductToken<T> art::Consumer::mayConsume ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ProductToken<T> art::Consumer::mayConsume ( InputTag const &  it)
inherited

Definition at line 189 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

190 {
191  if (!moduleContext_)
192  return ProductToken<T>::invalid();
193 
194  consumables_[BT].emplace_back(ConsumableType::Product,
195  TypeID{typeid(T)},
196  it.label(),
197  it.instance(),
198  it.process());
199  return ProductToken<T>{it};
200 }
it
set< int >::iterator it
Definition: ProduceBadRunsByFailureMode.C:101
art::ProductToken::invalid
static ProductToken< T > invalid()
Definition: ProductToken.h:47
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:137
T
double T
Definition: Xdiff_gwt.C:5
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:135
template<typename T , art::BranchType BT>
void art::Consumer::mayConsumeMany ( )
inherited

Definition at line 204 of file Consumer.h.

References T.

205 {
206  if (!moduleContext_)
207  return;
208 
209  consumables_[BT].emplace_back(ConsumableType::Many, TypeID{typeid(T)});
210 }
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:137
T
double T
Definition: Xdiff_gwt.C:5
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:135
template<typename Element , BranchType = InEvent>
ViewToken<Element> art::Consumer::mayConsumeView ( InputTag const &  )
inherited
template<typename T , art::BranchType BT>
art::ViewToken<T> art::Consumer::mayConsumeView ( InputTag const &  it)
inherited

Definition at line 214 of file Consumer.h.

References art::InputTag::instance(), PandAna.reco_validation.prod5_pid_validation::invalid, art::InputTag::label(), art::InputTag::process(), and T.

215 {
216  if (!moduleContext_)
217  return ViewToken<T>::invalid();
218 
219  consumables_[BT].emplace_back(ConsumableType::ViewElement,
220  TypeID{typeid(T)},
221  it.label(),
222  it.instance(),
223  it.process());
224  return ViewToken<T>{it};
225 }
it
set< int >::iterator it
Definition: ProduceBadRunsByFailureMode.C:101
art::ViewToken::invalid
static ViewToken< Element > invalid()
Definition: ProductToken.h:75
art::Consumer::consumables_
ConsumableProducts consumables_
Definition: Consumer.h:137
T
double T
Definition: Xdiff_gwt.C:5
art::Consumer::moduleContext_
bool moduleContext_
Definition: Consumer.h:135
bool art::ProducerBase::modifiesEvent ( ) const
inlineinherited

Definition at line 40 of file ProducerBase.h.

References art::ProducerBase::getProductID(), and string.

41  {
42  return true;
43  }
static cet::exempt_ptr<Consumer> art::Consumer::non_module_context ( )
staticinherited

Referenced by art::RootOutput::beginSubRun(), art::ProcessPackage< L >::postScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::Run >::End::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::postScheduleSignal(), art::BeginEndPackage< Level::SubRun >::End::postScheduleSignal(), art::ProcessPackage< L >::preScheduleSignal(), art::BeginEndPackage< Level::Run >::Begin::preScheduleSignal(), art::BeginEndPackage< Level::SubRun >::Begin::preScheduleSignal(), art::EmptyEvent::readRun_(), and art::EmptyEvent::readSubRun_().

void art::Consumer::prepareForJob ( fhicl::ParameterSet const &  pset)
protectedinherited
void lem::MakePID::produce ( art::Event &  evt)
virtual

Implements art::EDProducer.

Definition at line 326 of file MakePID_module.cc.

References ana::assert(), AvgAndFit(), plot_validation_datamc::c, lem::dec::Forest::Classify(), util::CreateAssn(), DEFINE_ART_MODULE(), e, stan::math::exp(), lem::PIDExtraVars::fAvgExpE, lem::PIDExtraVars::fAvgInvE, lem::PIDExtraVars::fAvgSigE, lem::PIDExtraVars::fBestBkgPotential, lem::PIDExtraVars::fBestEnrichPotential, lem::PIDExtraVars::fBestSigPotential, lem::PIDExtraVars::fChiBkg, lem::PIDExtraVars::fChiBkgExp, lem::PIDExtraVars::fChiSig, lem::PIDExtraVars::fChiSigExp, lem::PIDExtraVars::fDecTreePID, lem::PIDExtraVars::fDISNCFracExp, lem::PIDExtraVars::fDISNumuFracExp, lem::PIDExtraVars::fDISSigFracExp, lem::PIDExtraVars::fEnergyDiff, lem::PIDExtraVars::fEnergyDiffExp, lem::PIDExtraVars::fEnrichFrac, lem::PIDExtraVars::fEnrichFracExp, lem::PIDExtraVars::fEnrichFracFit, lem::PIDExtraVars::fEnrichQFrac, lem::PIDExtraVars::fEnrichQFracExp, lem::PIDExtraVars::fEnrichQFracFit, fEnrichScale, fExpLength, fExpPow, fExpTable, fFlatHist, fForest, fMatchesLabel, lem::PIDExtraVars::fMeanL0Exp, lem::PIDExtraVars::fMeanL0Fit, lem::PIDExtraVars::fMeanL0PiOnlyExp, lem::PIDExtraVars::fMeanL0PiOnlyFit, lem::PIDExtraVars::fMeanL1Exp, lem::PIDExtraVars::fMeanL1Fit, lem::PIDExtraVars::fMeanL1PiOnlyExp, lem::PIDExtraVars::fMeanL1PiOnlyFit, lem::PIDExtraVars::fMeanMatchEExp, lem::PIDExtraVars::fMeanQFrac, lem::PIDExtraVars::fMeanQFracBkg, lem::PIDExtraVars::fMeanQFracBkgExp, lem::PIDExtraVars::fMeanQFracBkgFit, lem::PIDExtraVars::fMeanQFracExp, lem::PIDExtraVars::fMeanQFracFit, lem::PIDExtraVars::fMeanQFracSig, lem::PIDExtraVars::fMeanQFracSigExp, lem::PIDExtraVars::fMeanQFracSigFit, lem::PIDExtraVars::fMeanQFracSigY50Fit, lem::PIDExtraVars::fMeanQFracSigY90Fit, lem::PIDExtraVars::fMeanY, lem::PIDExtraVars::fMeanYExp, lem::PIDExtraVars::fMeanYFit, lem::PIDExtraVars::fMeanYSig, lem::PIDExtraVars::fMeanYSigExp, lem::PIDExtraVars::fMeanYSigFit, lem::PIDExtraVars::fMeanYSigY50Fit, lem::PIDExtraVars::fMeanYSigY90Fit, PandAna.Demos.pi0_spectra::fmt, fNCScale, lem::PIDExtraVars::fPIDExp, lem::PIDExtraVars::fPIDExpFit, lem::PIDExtraVars::fPIDFit, lem::PIDExtraVars::fPIDFitUnweight, lem::PIDExtraVars::fPIDFracNumuFit, lem::PIDExtraVars::fPIDy50Fit, lem::PIDExtraVars::fPIDy90Fit, fPositionalWeight, lem::PIDExtraVars::fQENCFracExp, lem::PIDExtraVars::fQENumuFracExp, lem::PIDExtraVars::fQESigFracExp, fReader, lem::PIDExtraVars::fResNCFracExp, lem::PIDExtraVars::fResNumuFracExp, lem::PIDExtraVars::fResSigFracExp, fScaleCalE, fSigScale, fSummarizeLabel, lem::PIDExtraVars::fTMVABDT, lem::PIDExtraVars::fTMVABDTD, lem::PIDExtraVars::fTMVABDTG, lem::PIDExtraVars::fTMVABDTG_CC, lem::PIDExtraVars::fTMVABDTMF, lem::PIDExtraVars::fTMVAMLP, lem::PIDExtraVars::fTMVAMLP_CC, fTrimMatchLists, fUseFitInput, art::DataViewImpl::getByLabel(), gTMVAVar, MECModelEnuComparisons::i, dumpEventsToText::inputs, kBkgFrac, kBkgOnly, kBkgOnlyExp, kDISNCFrac, kDISNumuFrac, kDISSigFrac, kExp, kL0, kL1, kNumuFrac, kPi0Only, kQENCFrac, kQENumuFrac, kQESigFrac, kQFrac, kResNCFrac, kResNumuFrac, kResSigFrac, kSig50Only, kSig90Only, kSigFrac, kSigFrac50, kSigFrac90, kSigOnly, kSigOnlyExp, kSimple, kTrueEVis, kUnweighted, kY, m, lem::MatchList::matches, overlay_prestage_def::matches, lem::PIDDetails::matchList, getGoodRuns4SAM::n, num, BlessedPlots::pid, cet::pow(), art::Event::put(), util::sqr(), lem::dec::Evt::vars, lem::PIDDetails::vars, vars, and w.

327  {
328  std::unique_ptr<std::vector<lem::PIDDetails> > pidcol(new std::vector<lem::PIDDetails>);
329  std::unique_ptr<art::Assns<lem::PIDDetails, rb::Cluster>> assns(new art::Assns<lem::PIDDetails, rb::Cluster>);
330 
331  art::Handle<std::vector<lem::MatchList> > matchlists;
332  evt.getByLabel(fMatchesLabel, matchlists);
333 
334  art::Handle<std::vector<LEMInput> > inputs;
335  evt.getByLabel(fSummarizeLabel, inputs);
336 
337  assert(inputs->size() == matchlists->size());
338 
339  art::FindManyP<rb::Cluster> fmt(matchlists, evt, fMatchesLabel);
340 
341  const unsigned int matchMax = matchlists->size();
342  for(unsigned int matchIdx = 0; matchIdx < matchMax; ++matchIdx){
343  MatchList matchlist = (*matchlists)[matchIdx];
344 
345  // Get all the slices associated with this match list
346  std::vector<art::Ptr<rb::Cluster> > slices = fmt.at(matchIdx);
347  // Associated slices can be empty when we're running over lemsum files,
348  // where the slices have been stripped out. In that case, don't make
349  // associations, LEMAssociator will sort it out later.
350  assert(slices.empty() || slices.size() == 1);
351 
352  const double totalRecoGeV = (*inputs)[matchIdx].totalGeV;
353 
354  std::vector<MatchSummary> matches, enrichedMatches;
355  for(unsigned int i = 0; i < matchlist.matches.size(); ++i){
356  if(matchlist.matches[i].enrich){
357  enrichedMatches.push_back(matchlist.matches[i]);
358  }
359  else{
360  matches.push_back(matchlist.matches[i]);
361  }
362  }
363 
364  if(matches.empty()){
365  lem::PIDDetails pid(-1);
366  pidcol->push_back(pid);
367  if(!slices.empty())
368  util::CreateAssn(*this, evt, *pidcol, slices[0], *assns);
369  continue;
370  }
371 
372  // Equalize stats
373  for(unsigned int n = 0; n < matches.size(); ++n){
374  if(matches[n].isSig) matches[n].weight *= fSigScale;
375  if(matches[n].ccnc == 1) matches[n].weight *= fNCScale;
376  }
377  for(unsigned int n = 0; n < enrichedMatches.size(); ++n){
378  enrichedMatches[n].weight *= fEnrichScale;
379  }
380 
381 
382  PIDExtraVars vars;
383 
384  std::vector<double> xs, ys, ws, ws2;
385 
386  double meanESig = 0;
387  double meanEBkg = 0;
388  double meanESigExp = 0;
389  double meanEBkgExp = 0;
390  double totSig = 0;
391  double totBkg = 0;
392  double totSigExp = 0;
393  double totBkgExp = 0;
394  double totW = 0;
395 
396  double invE = 0, invEDenom = 0;
397  double expE = 0, expEDenom = 0;
398  double sigE = 0, sigEDenom = 0;
399 
400  const double maxPotential = matches.empty() ? 1 : matches[matches.size()-1].potential;
401  for(unsigned int i = 0; i < matches.size(); ++i){
402  const double w = matches[i].weight;
403 
404  const double invEVal = w/(matches[i].potential+1e-5);
405  if(matches[i].isSig) invE += invEVal; else invE -= invEVal;
406  invEDenom += invEVal;
407 
408  const double expEVal = w*exp(-matches[i].potential);
409  if(matches[i].isSig) expE += expEVal; else expE -= expEVal;
410  expEDenom += expEVal;
411 
412  const double sigEVal = w/(1+matches[i].potential);
413  if(matches[i].isSig) sigE += sigEVal; else sigE -= sigEVal;
414  sigEDenom += sigEVal;
415 
416  xs.push_back(matches[i].potential);
417  ys.push_back(matches[i].isSig);
418  ws.push_back(w);
419  double expVal = fExpTable[i];
420  if(!fPositionalWeight){
421  expVal = exp(-pow(matches[i].potential/maxPotential, fExpPow)/fExpLength);
422  }
423  ws2.push_back(w*expVal);
424 
425  if(matches[i].isSig){
426  totSig += w;
427  totSigExp += expVal*w;
428  meanESig += w*matches[i].potential;
429  meanESigExp += expVal*w*matches[i].potential;
430  }
431  else{
432  totBkg += w;
433  totBkgExp += expVal*w;
434  meanEBkg += w*matches[i].potential;
435  meanEBkgExp += expVal*w*matches[i].potential;
436  }
437 
438  totW += w;
439  } // end for i
440 
441  if(ws.size() < 2){
442  vars.fChiSig = 0;
443  vars.fChiBkg = 0;
444  vars.fChiSigExp = 0;
445  vars.fChiBkgExp = 0;
446  }
447  else{
448  for(int c = 0; c <= 1; ++c){
449  double num = 0, denom = 0;
450  for(unsigned int i = 0; i < xs.size(); ++i){
451  num += ws[i]*(xs[i]*ys[i]-c*xs[i]);
452  denom += xs[i]*xs[i]*ws[i];
453  }
454  const double m = num/denom;
455  double chi = 0;
456  double totW = 0;
457  for(unsigned int i = 0; i < xs.size(); ++i){
458  chi += ws[i]*util::sqr(ys[i]-m*xs[i]-c);
459  totW += ws[i];
460  }
461  chi /= totW;
462  if(c == 0) vars.fChiBkg = chi;
463  if(c == 1) vars.fChiSig = chi;
464  } // end for c
465 
466  for(int c = 0; c <= 1; ++c){
467  double num = 0, denom = 0;
468  for(unsigned int i = 0; i < xs.size(); ++i){
469  num += ws2[i]*(xs[i]*ys[i]-c*xs[i]);
470  denom += xs[i]*xs[i]*ws2[i];
471  }
472  const double m = num/denom;
473  double chi = 0;
474  double totW = 0;
475  for(unsigned int i = 0; i < xs.size(); ++i){
476  chi += ws2[i]*util::sqr(ys[i]-m*xs[i]-c);
477  totW += ws2[i];
478  }
479  chi /= totW;
480  if(c == 0) vars.fChiBkgExp = chi;
481  if(c == 1) vars.fChiSigExp = chi;
482  } // end for c
483  }
484 
485  if(totSig) meanESig /= totSig;
486  if(totBkg) meanEBkg /= totBkg;
487  if(totSigExp) meanESigExp /= totSigExp;
488  if(totBkgExp) meanEBkgExp /= totBkgExp;
489 
490 
491  double junk;
492 
493  AvgAndFit(matches, kSimple, kSigFrac, junk, vars.fPIDFit);
494  AvgAndFit(matches, kUnweighted, kSigFrac, junk, vars.fPIDFitUnweight);
495  AvgAndFit(matches, kExp, kSigFrac, vars.fPIDExp, vars.fPIDExpFit);
496 
497  AvgAndFit(matches, kSimple, kY, vars.fMeanY, vars.fMeanYFit);
498  AvgAndFit(matches, kExp, kY, vars.fMeanYExp, junk);
499 
500  AvgAndFit(matches, kSigOnly, kY, vars.fMeanYSig, vars.fMeanYSigFit);
501  AvgAndFit(matches, kSigOnlyExp, kY, vars.fMeanYSigExp, junk);
502 
503  AvgAndFit(matches, kSimple, kQFrac, vars.fMeanQFrac, vars.fMeanQFracFit);
504  AvgAndFit(matches, kExp, kQFrac, vars.fMeanQFracExp, junk);
505 
506  AvgAndFit(matches, kSigOnly, kQFrac, vars.fMeanQFracSig, vars.fMeanQFracSigFit);
507  AvgAndFit(matches, kSigOnlyExp, kQFrac, vars.fMeanQFracSigExp, junk);
508 
509  AvgAndFit(matches, kBkgOnly, kQFrac, vars.fMeanQFracBkg, vars.fMeanQFracBkgFit);
510  AvgAndFit(matches, kBkgOnlyExp, kQFrac, vars.fMeanQFracBkgExp, junk);
511 
512  AvgAndFit(matches, kExp, kTrueEVis, vars.fMeanMatchEExp, junk);
513 
514  AvgAndFit(matches, kExp, kQESigFrac, vars.fQESigFracExp, junk);
515  AvgAndFit(matches, kExp, kResSigFrac, vars.fResSigFracExp, junk);
516  AvgAndFit(matches, kExp, kDISSigFrac, vars.fDISSigFracExp, junk);
517  AvgAndFit(matches, kExp, kQENCFrac, vars.fQENCFracExp, junk);
518  AvgAndFit(matches, kExp, kResNCFrac, vars.fResNCFracExp, junk);
519  AvgAndFit(matches, kExp, kDISNCFrac, vars.fDISNCFracExp, junk);
520  AvgAndFit(matches, kExp, kQENumuFrac, vars.fQENumuFracExp, junk);
521  AvgAndFit(matches, kExp, kResNumuFrac, vars.fResNumuFracExp, junk);
522  AvgAndFit(matches, kExp, kDISNumuFrac, vars.fDISNumuFracExp, junk);
523 
524  // Mix in the usual signal matches
525  const int emOrig = enrichedMatches.size();
526  for(unsigned int n = 0; n < matches.size(); ++n){
527  if(matches[n].isSig) enrichedMatches.push_back(matches[n]);
528  }
529  std::sort(enrichedMatches.begin(), enrichedMatches.end());
530  enrichedMatches.resize(emOrig);
531 
532  AvgAndFit(enrichedMatches, kSimple, kBkgFrac, vars.fEnrichFrac, vars.fEnrichFracFit);
533  AvgAndFit(enrichedMatches, kExp, kBkgFrac, vars.fEnrichFracExp, junk);
534 
535  AvgAndFit(enrichedMatches, kBkgOnly, kQFrac, vars.fEnrichQFrac, vars.fEnrichQFracFit);
536  AvgAndFit(enrichedMatches, kBkgOnlyExp, kQFrac, vars.fEnrichQFracExp, junk);
537 
538  AvgAndFit(enrichedMatches, kSimple, kL0, junk, vars.fMeanL0Fit);
539  AvgAndFit(enrichedMatches, kSimple, kL1, junk, vars.fMeanL1Fit);
540  AvgAndFit(enrichedMatches, kExp, kL0, vars.fMeanL0Exp, junk);
541  AvgAndFit(enrichedMatches, kExp, kL1, vars.fMeanL1Exp, junk);
542 
543  AvgAndFit(enrichedMatches, kPi0Only, kL0, junk, vars.fMeanL0PiOnlyFit);
544  AvgAndFit(enrichedMatches, kPi0Only, kL1, junk, vars.fMeanL1PiOnlyFit);
545  AvgAndFit(enrichedMatches, kExp, kL0, vars.fMeanL0PiOnlyExp, junk);
546  AvgAndFit(enrichedMatches, kExp, kL1, vars.fMeanL1PiOnlyExp, junk);
547 
548  AvgAndFit(matches, kSimple, kSigFrac50, junk, vars.fPIDy50Fit);
549  AvgAndFit(matches, kSimple, kSigFrac90, junk, vars.fPIDy90Fit);
550 
551  AvgAndFit(matches, kSig50Only, kQFrac, junk, vars.fMeanQFracSigY50Fit);
552  AvgAndFit(matches, kSig90Only, kQFrac, junk, vars.fMeanQFracSigY90Fit);
553 
554  AvgAndFit(matches, kSig50Only, kY, junk, vars.fMeanYSigY50Fit);
555  AvgAndFit(matches, kSig90Only, kY, junk, vars.fMeanYSigY90Fit);
556 
557  AvgAndFit(matches, kSimple, kNumuFrac, junk, vars.fPIDFracNumuFit);
558 
559  vars.fEnergyDiff = meanESig-meanEBkg;
560  vars.fEnergyDiffExp = meanESigExp-meanEBkgExp;
561 
562  vars.fBestSigPotential = vars.fBestBkgPotential = vars.fBestEnrichPotential = 2;
563  for(unsigned int n = 0; n < matches.size(); ++n){
564  if(matches[n].isSig){
565  vars.fBestSigPotential = matches[n].potential;
566  break;
567  }
568  }
569  for(unsigned int n = 0; n < matches.size(); ++n){
570  if(!matches[n].isSig){
571  vars.fBestBkgPotential = matches[n].potential;
572  break;
573  }
574  }
575  for(unsigned int n = 0; n < enrichedMatches.size(); ++n){
576  if(!enrichedMatches[n].isSig){
577  vars.fBestEnrichPotential = enrichedMatches[n].potential;
578  break;
579  }
580  }
581 
582  vars.fAvgInvE = invE/invEDenom;
583  vars.fAvgExpE = expE/expEDenom;
584  vars.fAvgSigE = sigE/sigEDenom;
585 
586  // if(fANN || fANNEnrich){
587  // fann_type input[6];
588  // input[0] = fUseFitInput ? vars.fPIDFit : vars.fPIDExp;
589  // input[1] = vars.fMeanYExp;
590  // input[2] = vars.fMeanQFracExp;
591  // input[3] = vars.fEnergyDiffExp;
592 
593  // if(fANN){
594  // fann_type* calc_out = fann_run(fANN, input);
595  // vars.fFannId = calc_out[0];
596  // }
597 
598  // input[4] = 2*totalRecoGeV;
599  // input[5] = vars.fEnrichFracExp;
600 
601  // if(fANNEnrich){
602  // fann_type* calc_out = fann_run(fANNEnrich, input);
603  // vars.fFannIdEnrich = calc_out[0];
604  // }
605  // }
606 
607  if(fReader){
608  gTMVAVar[0] = fUseFitInput ? vars.fPIDFit : vars.fPIDExp;
609  gTMVAVar[1] = vars.fMeanYExp;
610  gTMVAVar[2] = vars.fMeanQFracExp;
611  gTMVAVar[3] = vars.fEnergyDiffExp;
612  if(fScaleCalE)
613  gTMVAVar[4] = 2*totalRecoGeV;
614  else
615  gTMVAVar[4] = totalRecoGeV;
616  gTMVAVar[5] = vars.fEnrichFracExp;
617 
618  vars.fTMVABDTG = fReader->EvaluateMVA("BDTG");
619  vars.fTMVABDT = fReader->EvaluateMVA("BDT");
620  vars.fTMVABDTD = fReader->EvaluateMVA("BDTD");
621  vars.fTMVABDTMF = fReader->EvaluateMVA("BDTMitFisher");
622  vars.fTMVAMLP = fReader->EvaluateMVA("MLP");
623 
624  vars.fTMVABDTG_CC = fReader->EvaluateMVA("BDTG_CC");
625  vars.fTMVAMLP_CC = fReader->EvaluateMVA("MLP_CC");
626  }
627 
628  if(fForest){
629  if(vars.fPIDExp == -1){
630  vars.fDecTreePID = -1;
631  }
632  else{
633  dec::Evt decevt;
634  decevt.vars[0] = fUseFitInput ? vars.fPIDFit : vars.fPIDExp;
635  decevt.vars[1] = vars.fMeanYExp;
636  decevt.vars[2] = vars.fMeanQFracExp;
637  decevt.vars[3] = vars.fEnergyDiffExp;
638  if(fScaleCalE)
639  decevt.vars[4] = 2*totalRecoGeV;
640  else
641  decevt.vars[4] = totalRecoGeV;
642  decevt.vars[5] = vars.fEnrichFracExp;
643  vars.fDecTreePID = fForest->Classify(decevt);
644  }
645  }
646 
647  double defaultPID = vars.fDecTreePID;
648 
649  // Flattened version of the Decision Tree variable
650  if(fFlatHist) defaultPID = fFlatHist->GetBinContent(fFlatHist->FindBin(vars.fDecTreePID));
651 
652 
653  lem::PIDDetails pid(defaultPID);
654 
655  if(fTrimMatchLists) matchlist.matches.resize(1);
656 
657  // TODO: sample size equalization doesn't survive this
658  pid.matchList = matchlist;
659  pid.vars = vars;
660 
661  pidcol->push_back(pid);
662 
663  if(!slices.empty()){
664  util::CreateAssn(*this, evt, *pidcol, slices[0], *assns);
665  }
666  } // end for matchIdx
667 
668  evt.put(std::move(pidcol));
669  evt.put(std::move(assns));
670  }
util::CreateAssn
static bool CreateAssn(art::EDProducer const &prod, art::Event &evt, std::vector< T > &a, art::Ptr< U > b, art::Assns< T, U > &assn, size_t indx=UINT_MAX, std::string const &instance=std::string())
Create a 1 to 1 association between a new product and one already in the event.
cet::pow
constexpr T pow(T x)
Definition: pow.h:75
lem::MakePID::fExpTable
std::vector< double > fExpTable
Definition: MakePID_module.cc:101
util::sqr
T sqr(T x)
More efficient square function than pow(x,2)
Definition: MathUtil.h:23
overlay_prestage_def.matches
matches
Definition: overlay_prestage_def.py:135
m
#define m
Definition: MakeGENIEReweightTable_module.cc:407
art::Event::put
ProductID put(std::unique_ptr< PROD > &&product)
Definition: Event.h:102
lem::MakePID::fForest
dec::Forest * fForest
Definition: MakePID_module.cc:118
lem::dec::Forest::Classify
double Classify(const Evt &evt) const
Calculate the PID value for evt.
Definition: DecisionTree.cxx:49
gTMVAVar
float gTMVAVar[6]
Definition: MakePID_module.cc:39
stan::math::exp
fvar< T > exp(const fvar< T > &x)
Definition: exp.hpp:10
plot_validation_datamc.c
c
Definition: plot_validation_datamc.py:198
vars
const std::map< std::pair< std::string, std::string >, Variable > vars
Definition: DataMCComp_header.h:233
lem::MakePID::fSummarizeLabel
std::string fSummarizeLabel
Definition: MakePID_module.cc:115
lem::MakePID::fMatchesLabel
std::string fMatchesLabel
Definition: MakePID_module.cc:114
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &productInstanceName, Handle< PROD > &result) const
Definition: DataViewImpl.h:344
num
int num
Definition: f2_nu.C:119
ana::assert
assert(nhit_max >=nhit_nbins)
PandAna.Demos.pi0_spectra.fmt
fmt
Definition: pi0_spectra.py:172
lem::PIDDetails
Attach LEM-specific info to the base PID object.
Definition: PIDDetails.h:20
e
Float_t e
Definition: plot.C:35
lem::MakePID::AvgAndFit
void AvgAndFit(const std::vector< MatchSummary > &matches, WeightScheme weightScheme, MatchVar matchVar, double &mean, double &fit) const
Definition: MakePID_module.cc:235
w
Float_t w
Definition: plot.C:20
lem::MakePID::fReader
TMVA::Reader * fReader
Definition: MakePID_module.cc:117
void lem::MakePID::reconfigure ( const fhicl::ParameterSet &  pset)
virtual

Definition at line 194 of file MakePID_module.cc.

References fExpLength, fExpPow, fMatchesLabel, fPositionalWeight, fScaleCalE, fSummarizeLabel, fTrimMatchLists, fUseFitInput, fhicl::ParameterSet::get(), and string.

Referenced by MakePID().

195  {
196  fScaleCalE = pset.get<bool>("ScaleCalE");
197 
198  fTrimMatchLists = pset.get<bool>("TrimMatchLists");
199  fUseFitInput = pset.get<bool>("UseFitInput");
200  fPositionalWeight = pset.get<bool>("PositionalWeight");
201  fMatchesLabel = pset.get<std::string>("MatchesLabel");
202  fSummarizeLabel = pset.get<std::string>("SummarizeLabel");
203  fExpLength = pset.get<double>("ExpLength");
204  fExpPow = pset.get<double>("ExpPow");
205  }
fhicl::ParameterSet::get
T get(std::string const &key) const
Definition: ParameterSet.h:231
lem::MakePID::fSummarizeLabel
std::string fSummarizeLabel
Definition: MakePID_module.cc:115
lem::MakePID::fMatchesLabel
std::string fMatchesLabel
Definition: MakePID_module.cc:114
string
enum BeamMode string
void art::Consumer::showMissingConsumes ( ) const
protectedinherited

Referenced by art::RootOutput::endJob().

void art::Consumer::validateConsumedProduct ( BranchType const  bt,
ProductInfo const &  pi 
)
protectedinherited

Member Data Documentation

double lem::MakePID::fEnrichScale
protected

Definition at line 105 of file MakePID_module.cc.

Referenced by beginRun(), and produce().

double lem::MakePID::fExpLength
protected

Definition at line 107 of file MakePID_module.cc.

Referenced by AvgAndFit(), beginJob(), produce(), and reconfigure().

double lem::MakePID::fExpPow
protected

Definition at line 107 of file MakePID_module.cc.

Referenced by AvgAndFit(), beginJob(), produce(), and reconfigure().

std::vector<double> lem::MakePID::fExpTable
protected

Definition at line 101 of file MakePID_module.cc.

Referenced by AvgAndFit(), beginJob(), and produce().

TH1* lem::MakePID::fFlatHist
protected

Definition at line 120 of file MakePID_module.cc.

Referenced by MakePID(), and produce().

dec::Forest* lem::MakePID::fForest
protected

Definition at line 118 of file MakePID_module.cc.

Referenced by MakePID(), produce(), and ~MakePID().

std::string lem::MakePID::fMatchesLabel
protected

Definition at line 114 of file MakePID_module.cc.

Referenced by beginRun(), produce(), and reconfigure().

double lem::MakePID::fNCScale
protected

Definition at line 103 of file MakePID_module.cc.

Referenced by beginRun(), and produce().

bool lem::MakePID::fPositionalWeight
protected

Definition at line 113 of file MakePID_module.cc.

Referenced by AvgAndFit(), produce(), and reconfigure().

TMVA::Reader* lem::MakePID::fReader
protected

Definition at line 117 of file MakePID_module.cc.

Referenced by MakePID(), produce(), and ~MakePID().

bool lem::MakePID::fScaleCalE
protected

Definition at line 60 of file MakePID_module.cc.

Referenced by produce(), and reconfigure().

double lem::MakePID::fSigScale
protected

Definition at line 104 of file MakePID_module.cc.

Referenced by beginRun(), and produce().

std::string lem::MakePID::fSummarizeLabel
protected

Definition at line 115 of file MakePID_module.cc.

Referenced by produce(), and reconfigure().

bool lem::MakePID::fTrimMatchLists
protected

Definition at line 111 of file MakePID_module.cc.

Referenced by produce(), and reconfigure().

bool lem::MakePID::fUseFitInput
protected

Definition at line 112 of file MakePID_module.cc.

Referenced by produce(), and reconfigure().

The documentation for this class was generated from the following file:
Generated on Sat Nov 28 2020 05:19:51 for NOvA by   doxygen 1.8.11