Public Member Functions | Private Member Functions | Private Attributes | List of all members
pfpf::PFPAna Class Reference
Inheritance diagram for pfpf::PFPAna:
art::EDAnalyzer art::detail::Analyzer art::detail::LegacyModule art::Observer art::ModuleBase

Public Member Functions

 PFPAna (fhicl::ParameterSet const &pset)
 
- Public Member Functions inherited from art::EDAnalyzer
 EDAnalyzer (fhicl::ParameterSet const &pset)
 
template<typename Config >
 EDAnalyzer (Table< Config > const &config)
 
std::string workerType () const
 
- Public Member Functions inherited from art::detail::Analyzer
virtual ~Analyzer () noexcept
 
 Analyzer (fhicl::ParameterSet const &pset)
 
template<typename Config >
 Analyzer (Table< Config > const &config)
 
void doBeginJob (SharedResources const &resources)
 
void doEndJob ()
 
void doRespondToOpenInputFile (FileBlock const &fb)
 
void doRespondToCloseInputFile (FileBlock const &fb)
 
void doRespondToOpenOutputFiles (FileBlock const &fb)
 
void doRespondToCloseOutputFiles (FileBlock const &fb)
 
bool doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
 
bool doEndRun (RunPrincipal &rp, ModuleContext const &mc)
 
bool doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
 
bool doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
 
bool doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
 
- Public Member Functions inherited from art::Observer
 ~Observer () noexcept
 
 Observer (Observer const &)=delete
 
 Observer (Observer &&)=delete
 
Observer & operator= (Observer const &)=delete
 
Observer & operator= (Observer &&)=delete
 
void registerProducts (ProductDescriptions &, ModuleDescription const &)
 
void fillDescriptions (ModuleDescription const &)
 
fhicl::ParameterSetID selectorConfig () const
 
- Public Member Functions inherited from art::ModuleBase
virtual ~ModuleBase () noexcept
 
 ModuleBase ()
 
ModuleDescription const & moduleDescription () const
 
void setModuleDescription (ModuleDescription const &)
 
std::array< std::vector< ProductInfo >, NumBranchTypes > const & getConsumables () const
 
void sortConsumables (std::string const &current_process_name)
 
template<typename T , BranchType BT>
ViewToken< T > consumesView (InputTag const &tag)
 
template<typename T , BranchType BT>
ViewToken< T > mayConsumeView (InputTag const &tag)
 

Private Member Functions

void analyze (const art::Event &evt) override
 
void beginJob () override
 

Private Attributes

TH1F * fNClusters
 
TH1F * fNHitInCluster
 
TH1F * fCREP2
 
TH1F * fCRE
 
TH1F * fCRP
 
TH1F * fNuKE_elec
 
TH1F * fNuKE_muon
 
TH1F * fNuKE_pion
 
TH1F * fNuKE_kaon
 
TH1F * fNuKE_prot
 
TH1F * fNuEP2_elec
 
TH1F * fNuEP2_muon
 
TH1F * fNuEP2_pion
 
TH1F * fNuEP2_kaon
 
TH1F * fNuEP2_prot
 
TH1F * fNuE_elec
 
TH1F * fNuE_muon
 
TH1F * fNuE_pion
 
TH1F * fNuE_kaon
 
TH1F * fNuE_prot
 
TH1F * fNuP_elec
 
TH1F * fNuP_muon
 
TH1F * fNuP_pion
 
TH1F * fNuP_kaon
 
TH1F * fNuP_prot
 
TH1F * fNuVtx_dx
 
TH1F * fNuVtx_dy
 
TH1F * fNuVtx_dz
 
TProfile * fNuEP2_KE_elec
 
TProfile * fNuEP2_KE_muon
 
TProfile * fNuEP2_KE_pion
 
TProfile * fNuEP2_KE_kaon
 
TProfile * fNuEP2_KE_prot
 
std::string fHitsModuleLabel
 
std::string fClusterModuleLabel
 
std::string fTrackModuleLabel
 
std::string fPFParticleModuleLabel
 
std::string fVertexModuleLabel
 
std::vector< float > fElecKERange
 
std::vector< float > fMuonKERange
 
std::vector< float > fPionKERange
 
std::vector< float > fKaonKERange
 
std::vector< float > fProtKERange
 
short fTrackWeightOption
 
bool fMergeDaughters
 
bool fSkipCosmics
 
short fPrintLevel
 

Additional Inherited Members

- Public Types inherited from art::EDAnalyzer
using WorkerType = WorkerT< EDAnalyzer >
 
using ModuleType = EDAnalyzer
 
- Protected Member Functions inherited from art::Observer
std::string const & processName () const
 
bool wantAllEvents () const noexcept
 
bool wantEvent (ScheduleID id, Event const &e) const
 
Handle< TriggerResults > getTriggerResults (Event const &e) const
 
 Observer (fhicl::ParameterSet const &config)
 
 Observer (std::vector< std::string > const &select_paths, std::vector< std::string > const &reject_paths, fhicl::ParameterSet const &config)
 
- Protected Member Functions inherited from art::ModuleBase
ConsumesCollector & consumesCollector ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > consumes (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > consumesView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void consumesMany ()
 
template<typename T , BranchType = InEvent>
ProductToken< T > mayConsume (InputTag const &)
 
template<typename Element , BranchType = InEvent>
ViewToken< Element > mayConsumeView (InputTag const &)
 
template<typename T , BranchType = InEvent>
void mayConsumeMany ()
 

Detailed Description

Definition at line 43 of file PFPAna_module.cc.

Constructor & Destructor Documentation

pfpf::PFPAna::PFPAna ( fhicl::ParameterSet const &  pset)
explicit

Definition at line 107 of file PFPAna_module.cc.

108  : EDAnalyzer(pset)
109  , fHitsModuleLabel(pset.get<std::string>("HitsModuleLabel"))
110  , fClusterModuleLabel(pset.get<std::string>("ClusterModuleLabel"))
111  , fTrackModuleLabel(pset.get<std::string>("TrackModuleLabel"))
112  , fPFParticleModuleLabel(pset.get<std::string>("PFParticleModuleLabel"))
113  , fVertexModuleLabel(pset.get<std::string>("VertexModuleLabel"))
114  , fElecKERange(pset.get<std::vector<float>>("ElecKERange"))
115  , fMuonKERange(pset.get<std::vector<float>>("MuonKERange"))
116  , fPionKERange(pset.get<std::vector<float>>("PionKERange"))
117  , fKaonKERange(pset.get<std::vector<float>>("KaonKERange"))
118  , fProtKERange(pset.get<std::vector<float>>("ProtKERange"))
119  , fTrackWeightOption(pset.get<short>("TrackWeightOption"))
120  , fMergeDaughters(pset.get<bool>("MergeDaughters"))
121  , fSkipCosmics(pset.get<bool>("SkipCosmics"))
122  , fPrintLevel(pset.get<short>("PrintLevel"))
123  {}
pfpf::PFPAna::fElecKERange
std::vector< float > fElecKERange
Definition: PFPAna_module.cc:94
pfpf::PFPAna::fVertexModuleLabel
std::string fVertexModuleLabel
Definition: PFPAna_module.cc:93
string
std::string string
Definition: nybbler.cc:12
pfpf::PFPAna::fHitsModuleLabel
std::string fHitsModuleLabel
Definition: PFPAna_module.cc:89
pfpf::PFPAna::fProtKERange
std::vector< float > fProtKERange
Definition: PFPAna_module.cc:98
art::EDAnalyzer::EDAnalyzer
EDAnalyzer(fhicl::ParameterSet const &pset)
Definition: EDAnalyzer.h:25
pfpf::PFPAna::fClusterModuleLabel
std::string fClusterModuleLabel
Definition: PFPAna_module.cc:90
pfpf::PFPAna::fPionKERange
std::vector< float > fPionKERange
Definition: PFPAna_module.cc:96
pfpf::PFPAna::fTrackWeightOption
short fTrackWeightOption
Definition: PFPAna_module.cc:99
pfpf::PFPAna::fPFParticleModuleLabel
std::string fPFParticleModuleLabel
Definition: PFPAna_module.cc:92
pfpf::PFPAna::fMuonKERange
std::vector< float > fMuonKERange
Definition: PFPAna_module.cc:95
pfpf::PFPAna::fTrackModuleLabel
std::string fTrackModuleLabel
Definition: PFPAna_module.cc:91
pfpf::PFPAna::fKaonKERange
std::vector< float > fKaonKERange
Definition: PFPAna_module.cc:97

Member Function Documentation

void pfpf::PFPAna::analyze ( const art::Event &  evt)
overrideprivate

Definition at line 178 of file PFPAna_module.cc.

179  {
180  // code stolen from TrackAna_module.cc
181  art::ServiceHandle<geo::Geometry const> geom;
182  if (geom->Nplanes() > 3) return;
183 
184  // get all hits in the event
185  art::Handle<std::vector<recob::Hit>> hitListHandle;
186  evt.getByLabel(fHitsModuleLabel, hitListHandle);
187  std::vector<art::Ptr<recob::Hit>> allhits;
188  art::fill_ptr_vector(allhits, hitListHandle);
189  if (empty(allhits)) return;
190 
191  // get clusters and cluster-hit associations
192  art::Handle<std::vector<recob::Cluster>> clusterListHandle;
193  evt.getByLabel(fClusterModuleLabel, clusterListHandle);
194  art::FindManyP<recob::Hit> fmh(clusterListHandle, evt, fClusterModuleLabel);
195  if (clusterListHandle->size() == 0) return;
196 
197  // get 3D vertices
198  art::Handle<std::vector<recob::Vertex>> vertexListHandle;
199  evt.getByLabel(fVertexModuleLabel, vertexListHandle);
200  art::PtrVector<recob::Vertex> recoVtxList;
201  double xyz[3] = {0, 0, 0};
202  for (unsigned int ii = 0; ii < vertexListHandle->size(); ++ii) {
203  art::Ptr<recob::Vertex> vertex(vertexListHandle, ii);
204  recoVtxList.push_back(vertex);
205  vertex->XYZ(xyz);
206  }
207 
208  // get PFParticles
209  art::Handle<std::vector<recob::PFParticle>> PFPListHandle;
210  evt.getByLabel(fPFParticleModuleLabel, PFPListHandle);
211  art::PtrVector<recob::PFParticle> recoPFPList;
212  for (unsigned int ii = 0; ii < PFPListHandle->size(); ++ii) {
213  art::Ptr<recob::PFParticle> pfp(PFPListHandle, ii);
214  recoPFPList.push_back(pfp);
215  mf::LogVerbatim("PFPAna") << "PFParticle PDG " << pfp->PdgCode();
216  }
217 
218  // list of all true particles
219  art::ServiceHandle<cheat::BackTrackerService const> bt_serv;
220  art::ServiceHandle<cheat::ParticleInventoryService const> pi_serv;
221  sim::ParticleList const& plist = pi_serv->ParticleList();
222  // list of all true particles that will be considered
223  std::vector<const simb::MCParticle*> plist2;
224  // true (reconstructed) hits for each particle in plist2
225  std::vector<std::vector<art::Ptr<recob::Hit>>> hlist2;
226  // index of cluster matched to each particle in plist2 in each plane
227  std::vector<std::vector<short>> truToCl;
228  // number of true hits in each plane and cluster
229  std::vector<std::vector<unsigned short>> nTruHitInCl;
230  //number of reconstructed hits in all clusters
231  std::vector<unsigned short> nRecHitInCl;
232 
233  // calculate average EP2 for every event to facilitate code development
234  // Beam Neutrinos - muons and not-muons
235 
236  float aveNuEP2mu = 0.;
237  float numNuEP2mu = 0.;
238  float aveNuEP2nm = 0.;
239  float numNuEP2nm = 0.;
240  // Cosmic Rays
241  float aveCREP2 = 0.;
242  float numCREP2 = 0.;
243 
244  // track ID of the neutrino
245  int neutTrackID = -1;
246  std::vector<int> tidlist;
247  float neutEnergy = -1.;
248  int neutIntType = -1;
249  int neutCCNC = -1;
250 
251  auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
252 
253  for (sim::ParticleList::const_iterator ipart = plist.begin(); ipart != plist.end(); ++ipart) {
254  const simb::MCParticle* part = (*ipart).second;
255  assert(part != 0);
256  int pdg = abs(part->PdgCode());
257  int trackID = part->TrackId();
258  art::Ptr<simb::MCTruth> theTruth = pi_serv->TrackIdToMCTruth_P(trackID);
259  if (fSkipCosmics && theTruth->Origin() == simb::kCosmicRay) continue;
260 
261  if (fPrintLevel > 3)
262  mf::LogVerbatim("PFPAna") << "Pre-Cuts origin " << theTruth->Origin() << " trackID "
263  << trackID << " PDG " << part->PdgCode() << " E " << part->E()
264  << " mass " << part->Mass() << " Mother "
265  << part->Mother() + neutTrackID << " Proc " << part->Process();
266 
267  // Get the neutrino track ID. Assume that there is only one neutrino
268  // interaction and it is first in the list of BeamNeutrino particles
269  if (theTruth->Origin() == simb::kBeamNeutrino && neutTrackID < 0) {
270  neutTrackID = trackID;
271  simb::MCNeutrino theNeutrino = theTruth->GetNeutrino();
272  neutEnergy = 1000. * theNeutrino.Nu().E();
273  neutIntType = theNeutrino.InteractionType();
274  neutCCNC = theNeutrino.CCNC();
275  for (unsigned short iv = 0; iv < recoVtxList.size(); ++iv) {
276  recoVtxList[iv]->XYZ(xyz);
277  fNuVtx_dx->Fill(part->Vx() - xyz[0]);
278  fNuVtx_dy->Fill(part->Vy() - xyz[1]);
279  fNuVtx_dz->Fill(part->Vz() - xyz[2]);
280  } // iv
281  } // theTruth->Origin() == simb::kBeamNeutrino && neutTrackID <
282 
283  bool isCharged = (pdg == 11) || (pdg == 13) || (pdg == 211) || (pdg == 321) || (pdg == 2212);
284 
285  if (!isCharged) continue;
286 
287  float KE = 1000 * (part->E() - part->Mass());
288  // KE (MeV) cuts
289  if (pdg == 11) {
290  if (fElecKERange[0] < 0) continue;
291  // only allow primary electrons
292  if (part->Process() != "primary") continue;
293  if (KE < fElecKERange[0] || KE > fElecKERange[1]) continue;
294  }
295  if (pdg == 13) {
296  if (fMuonKERange[0] < 0) continue;
297  if (KE < fMuonKERange[0] || KE > fMuonKERange[1]) continue;
298  }
299  if (pdg == 211) {
300  if (fPionKERange[0] < 0) continue;
301  if (KE < fPionKERange[0] || KE > fPionKERange[1]) continue;
302  }
303  if (pdg == 321) {
304  if (fKaonKERange[0] < 0) continue;
305  if (KE < fKaonKERange[0] || KE > fKaonKERange[1]) continue;
306  }
307  if (pdg == 2212) {
308  if (fProtKERange[0] < 0) continue;
309  if (KE < fProtKERange[0] || KE > fProtKERange[1]) continue;
310  }
311  // ignore secondaries from neutron interactions
312  if (part->Process() == "NeutronInelastic") continue;
313  plist2.push_back(part);
314  tidlist.push_back(trackID);
315  // initialize the true->(cluster,plane) association
316  std::vector<short> temp{-1, -1, -1};
317  truToCl.push_back(temp);
318  // initialize the true hit count
319  std::vector<unsigned short> temp2(3);
320  nTruHitInCl.push_back(temp2);
321 
322  if (fPrintLevel > 2)
323  mf::LogVerbatim("PFPAna") << plist2.size() - 1 << " Origin " << theTruth->Origin()
324  << " trackID " << trackID << " PDG " << part->PdgCode() << " KE "
325  << (int)KE << " Mother " << part->Mother() + neutTrackID
326  << " Proc " << part->Process();
327  }
328 
329  if (empty(plist2)) return;
330 
331  // get the hits (in all planes) that are matched to the true tracks
332  hlist2 = bt_serv->TrackIdsToHits_Ps(clockData, tidlist, allhits);
333  if (hlist2.size() != plist2.size()) {
334  mf::LogError("PFPAna") << "MC particle list size " << plist2.size()
335  << " != size of MC particle true hits lists " << hlist2.size();
336  return;
337  }
338  tidlist.clear();
339 
340  // vector of (mother, daughter) pairs
341  std::vector<std::pair<unsigned short, unsigned short>> moda;
342  // Deal with mother-daughter tracks
343  if (fMergeDaughters && neutTrackID >= 0) {
344  // Assume that daughters appear later in the list. Step backwards
345  // to accumulate all generations of daughters
346  for (unsigned short dpl = plist2.size() - 1; dpl > 0; --dpl) {
347  // no mother
348  if (plist2[dpl]->Mother() == 0) continue;
349  // electron
350  if (abs(plist2[dpl]->PdgCode()) == 11) continue;
351  // the actual mother trackID is offset from the neutrino trackID
352  int motherID = neutTrackID + plist2[dpl]->Mother() - 1;
353  // ensure that we are only looking at BeamNeutrino daughters
354  if (motherID < 0) continue;
355  // count the number of daughters
356  int ndtr = 0;
357  for (unsigned short kpl = 0; kpl < plist2.size(); ++kpl) {
358  if (plist2[kpl]->Mother() == motherID) ++ndtr;
359  }
360  // require only one daughter
361  if (ndtr > 1) continue;
362  // find the mother in the list
363  int mpl = -1;
364  for (unsigned short jpl = dpl - 1; jpl > 0; --jpl) {
365  if (plist2[jpl]->TrackId() == motherID) {
366  mpl = jpl;
367  break;
368  }
369  } // jpl
370  // mother not found for some reason
371  if (mpl < 0) continue;
372  // ensure that PDG code for mother and daughter are the same
373  if (plist2[dpl]->PdgCode() != plist2[mpl]->PdgCode()) continue;
374  moda.push_back(std::make_pair(mpl, dpl));
375  } // dpl
376  } // MergeDaughters
377 
378  // Now match reconstructed clusters to true particles.
379  art::PtrVector<recob::Cluster> clusters;
380  for (unsigned int ii = 0; ii < clusterListHandle->size(); ++ii) {
381  art::Ptr<recob::Cluster> clusterHolder(clusterListHandle, ii);
382  clusters.push_back(clusterHolder);
383  }
384 
385  fNClusters->Fill(clusterListHandle->size());
386  nRecHitInCl.resize(clusters.size());
387 
388  // get the plane from the view. Perhaps there is a method that does
389  // this somewhere...
390  std::map<geo::View_t, unsigned int> ViewToPlane;
391  for (unsigned int plane = 0; plane < geom->Nplanes(); ++plane) {
392  geo::View_t view = geom->Plane(plane).View();
393  ViewToPlane[view] = plane;
394  }
395  for (size_t icl = 0; icl < clusters.size(); ++icl) {
396  unsigned int plane = ViewToPlane[clusters[icl]->View()];
397  std::vector<art::Ptr<recob::Hit>> cluhits = fmh.at(icl);
398  fNHitInCluster->Fill(cluhits.size());
399  nRecHitInCl[icl] = cluhits.size();
400  // count the number of hits matched to each true particle in plist2
401  std::vector<unsigned short> nHitInPl2(plist2.size());
402  for (size_t iht = 0; iht < cluhits.size(); ++iht) {
403  /*
404  mf::LogVerbatim("PFPAna")
405  <<"Clus Hit "<<cluhits[iht]->View()
406  <<":"<<cluhits[iht]->WireID().Wire
407  <<":"<<(int)cluhits[iht]->PeakTime();
408 */
409  // look for this hit in all of the truth hit lists
410  short hitInPl2 = -1;
411  for (unsigned short ipl = 0; ipl < plist2.size(); ++ipl) {
412  unsigned short imat = 0;
413  for (imat = 0; imat < hlist2[ipl].size(); ++imat) {
414  if (cluhits[iht] == hlist2[ipl][imat]) break;
415  } // imat
416  if (imat < hlist2[ipl].size()) {
417  hitInPl2 = ipl;
418  break;
419  }
420  } // ipl
421  if (hitInPl2 < 0) continue;
422  // Assign the hit count to the mother if this is a daughter.
423  // Mother-daughter pairs are entered in the moda vector in reverse
424  // order, so assign daughter hits to the highest generation mother.
425  for (unsigned short imd = 0; imd < moda.size(); ++imd) {
426  if (moda[imd].second == hitInPl2) hitInPl2 = moda[imd].first;
427  }
428  // count
429  ++nHitInPl2[hitInPl2];
430  } // iht
431  // Associate the cluster with the truth particle that has the highest
432  // number of cluster hits
433  unsigned short nhit = 0;
434  short imtru = -1;
435  for (unsigned int ipl = 0; ipl < nHitInPl2.size(); ++ipl) {
436  if (nHitInPl2[ipl] > nhit) {
437  nhit = nHitInPl2[ipl];
438  imtru = ipl;
439  }
440  } // ipl
441  // make the cluster->(true,plane) association and save the
442  // number of true hits in the cluster
443  if (imtru != -1) {
444  // clobber a previously made association?
445  if (nhit > nTruHitInCl[imtru][plane]) {
446  truToCl[imtru][plane] = icl;
447  nTruHitInCl[imtru][plane] = nhit;
448  }
449  } // imtru != 1
450  } // icl
451 
452  // ready to calculate Efficiency, Purity in each plane and EP2
453  for (unsigned short ipl = 0; ipl < plist2.size(); ++ipl) {
454  // ignore daughters
455  bool skipit = false;
456  for (unsigned short ii = 0; ii < moda.size(); ++ii) {
457  if (moda[ii].second == ipl) {
458  skipit = true;
459  break;
460  }
461  } // ii
462  if (skipit) continue;
463  // ignore true particles with few true hits. Outside the detector
464  // or not reconstructable?
465  if (hlist2[ipl].size() < 3) continue;
466 
467  int trackID = plist2[ipl]->TrackId();
468  art::Ptr<simb::MCTruth> theTruth = pi_serv->TrackIdToMCTruth_P(trackID);
469  bool isCosmic = (theTruth->Origin() == simb::kCosmicRay);
470  float KE = 1000 * (plist2[ipl]->E() - plist2[ipl]->Mass());
471  int PDG = abs(plist2[ipl]->PdgCode());
472 
473  std::vector<short> nTru(geom->Nplanes());
474  std::vector<short> nRec(geom->Nplanes());
475  std::vector<short> nTruRec(geom->Nplanes());
476  std::vector<float> eff(geom->Nplanes());
477  std::vector<float> pur(geom->Nplanes());
478  std::vector<float> ep(geom->Nplanes());
479  for (unsigned int plane = 0; plane < geom->Nplanes(); ++plane) {
480  // count the number of true hits in this plane for the true particle.
481  // First count the mother hits
482  for (unsigned short ii = 0; ii < hlist2[ipl].size(); ++ii) {
483  if (ViewToPlane[hlist2[ipl][ii]->View()] == plane) ++nTru[plane];
484  } // ii
485  // mf::LogVerbatim("PFPAna")
486  // <<"Chk mom "<<ipl<<" plane "<<plane<<" nTru "<<nTru[plane];
487  // next look for daughters and count those hits in all generations
488  unsigned short mom = ipl;
489  std::vector<std::pair<unsigned short, unsigned short>>::reverse_iterator rit =
490  moda.rbegin();
491  while (rit != moda.rend()) {
492  if ((*rit).first == mom) {
493  unsigned short dau = (*rit).second;
494  for (unsigned short jj = 0; jj < hlist2[dau].size(); ++jj) {
495  if (ViewToPlane[hlist2[dau][jj]->View()] == plane) ++nTru[plane];
496  } // jj
497  // It is likely that one hit appears in the mother list
498  // as well as the daughter list, so subtract one from the count
499  --nTru[plane];
500  mom = (*rit).second;
501  // mf::LogVerbatim("PFPAna")<<"new mom "<<mom<<" nTru
502  // "<<nTru[plane];
503  } // (*rit).first == mom
504  ++rit;
505  } // rit
506  // mf::LogVerbatim("PFPAna")<<"Chk dau "<<nTru[plane];
507  if (nTru[plane] == 0) {
508  // mf::LogVerbatim("PFPAna")<<"No true hits in plane "<<plane
509  // <<" for truth particle "<<ipl;
510  continue;
511  }
512  short icl = truToCl[ipl][plane];
513  nRec[plane] = nRecHitInCl[icl];
514  nTruRec[plane] = nTruHitInCl[ipl][plane];
515  // mf::LogVerbatim("PFPAna")<<"icl "<<icl<<" nRec "<<nRec[plane]
516  // <<" nTruRec "<<nTruRec[plane];
517  if (nTru[plane] > 0) eff[plane] = (float)nTruRec[plane] / (float)nTru[plane];
518  if (nRec[plane] > 0) pur[plane] = (float)nTruRec[plane] / (float)nRec[plane];
519  ep[plane] = eff[plane] * pur[plane];
520  } // plane
521  // sort the ep values in ascending order
522  std::vector<float> temp;
523  temp = ep;
524  std::sort(temp.begin(), temp.end());
525  // EP2 is the second highest value
526  unsigned short ii = temp.size() - 2;
527  float ep2 = temp[ii];
528  // find the plane that defined EP2
529  short ep2Plane = 0;
530  short ep2Cluster = 0;
531  for (unsigned short jj = 0; jj < temp.size(); ++jj) {
532  if (ep[jj] == ep2) {
533  ep2Plane = jj;
534  ep2Cluster = truToCl[ipl][ep2Plane];
535  break;
536  }
537  } // jj
538  // find the US and DS ends of the cluster for printing
539  std::array<double, 2> clBeg, clEnd;
540  if (ep2Cluster >= 0) {
541  clBeg[0] = clusters[ep2Cluster]->StartWire();
542  clBeg[1] = clusters[ep2Cluster]->StartTick();
543  clEnd[0] = clusters[ep2Cluster]->EndWire();
544  clEnd[1] = clusters[ep2Cluster]->EndTick();
545  }
546  else {
547  clBeg.fill(0.);
548  clEnd.fill(0.);
549  }
550  // fill histograms
551  if (isCosmic) {
552  fCREP2->Fill(ep2);
553  fCRE->Fill(eff[ep2Plane]);
554  fCRP->Fill(pur[ep2Plane]);
555  aveCREP2 += ep2;
556  numCREP2 += 1.;
557  if (fPrintLevel > 1)
558  mf::LogVerbatim("PFPAna")
559  << ">>>CREP2 " << std::fixed << std::setprecision(2) << ep2 << " E " << eff[ep2Plane]
560  << std::setprecision(2) << " P " << pur[ep2Plane] << " P:W:T " << ep2Plane << ":"
561  << (int)clBeg[0] << ":" << (int)clBeg[1] << "-" << ep2Plane << ":" << (int)clEnd[0]
562  << ":" << (int)clEnd[1] << " PDG " << PDG << " KE " << (int)KE << " MeV";
563  } // isCosmic
564  else {
565  float wght = 1.;
566  if (fTrackWeightOption == 1) wght = KE;
567  // accumulate statistics for muons and not-muons
568  if (PDG == 13) {
569  aveNuEP2mu += ep2 * wght;
570  numNuEP2mu += wght;
571  }
572  else {
573  aveNuEP2nm += ep2 * wght;
574  numNuEP2nm += wght;
575  }
576  if (PDG == 11) {
577  fNuKE_elec->Fill(KE, wght);
578  fNuE_elec->Fill(eff[ep2Plane], wght);
579  fNuP_elec->Fill(pur[ep2Plane], wght);
580  fNuEP2_elec->Fill(ep2, wght);
581  fNuEP2_KE_elec->Fill(KE, ep2, wght);
582  }
583  else if (PDG == 13) {
584  fNuKE_muon->Fill(KE, wght);
585  fNuE_muon->Fill(eff[ep2Plane], wght);
586  fNuP_muon->Fill(pur[ep2Plane], wght);
587  fNuEP2_muon->Fill(ep2, wght);
588  fNuEP2_KE_muon->Fill(KE, ep2, wght);
589  }
590  else if (PDG == 211) {
591  fNuKE_pion->Fill(KE, wght);
592  fNuE_pion->Fill(eff[ep2Plane], wght);
593  fNuP_pion->Fill(pur[ep2Plane], wght);
594  fNuEP2_pion->Fill(ep2, wght);
595  fNuEP2_KE_pion->Fill(KE, ep2, wght);
596  }
597  else if (PDG == 321) {
598  fNuKE_kaon->Fill(KE, wght);
599  fNuE_kaon->Fill(eff[ep2Plane], wght);
600  fNuP_kaon->Fill(pur[ep2Plane], wght);
601  fNuEP2_kaon->Fill(ep2, wght);
602  fNuEP2_KE_kaon->Fill(KE, ep2, wght);
603  }
604  else if (PDG == 2212) {
605  fNuKE_prot->Fill(KE, wght);
606  fNuE_prot->Fill(eff[ep2Plane], wght);
607  fNuP_prot->Fill(pur[ep2Plane], wght);
608  fNuEP2_prot->Fill(ep2, wght);
609  fNuEP2_KE_prot->Fill(KE, ep2, wght);
610  }
611  if (fPrintLevel > 1)
612  mf::LogVerbatim("PFPAna")
613  << ">>>NuEP2 " << std::fixed << std::setprecision(2) << ep2 << " E " << eff[ep2Plane]
614  << std::setprecision(2) << " P " << pur[ep2Plane] << " P:W:T " << ep2Plane << ":"
615  << (int)clBeg[0] << ":" << (int)clBeg[1] << "-" << ep2Plane << ":" << (int)clEnd[0]
616  << ":" << (int)clEnd[1] << " PDG " << PDG << " KE " << (int)KE << " MeV ";
617  if (fPrintLevel > 2) {
618  // print out the begin/end true hits
619  mf::LogVerbatim mfp("PFPAna");
620  mfp << " Truth P:W:T ";
621  for (unsigned int plane = 0; plane < geom->Nplanes(); ++plane) {
622  unsigned short loW = 9999;
623  unsigned short loT = 0;
624  unsigned short hiW = 0;
625  unsigned short hiT = 0;
626  for (unsigned short ii = 0; ii < hlist2[ipl].size(); ++ii) {
627  if (ViewToPlane[hlist2[ipl][ii]->View()] == plane) {
628  art::Ptr<recob::Hit> theHit = hlist2[ipl][ii];
629  if (theHit->WireID().Wire < loW) {
630  loW = theHit->WireID().Wire;
631  loT = theHit->PeakTime();
632  }
633  if (theHit->WireID().Wire > hiW) {
634  hiW = theHit->WireID().Wire;
635  hiT = theHit->PeakTime();
636  }
637  } // correct view
638  } // ii
639  mfp << plane << ":" << loW << ":" << loT << "-" << plane << ":" << hiW << ":" << hiT
640  << " ";
641  } // plane
642  } // fPrintLevel > 2
643  } // !isCosmic
644  } // ipl
645 
646  float ave1 = -1.;
647  if (numNuEP2mu > 0.) ave1 = aveNuEP2mu / numNuEP2mu;
648 
649  float ave2 = -1.;
650  if (numNuEP2nm > 0.) ave2 = aveNuEP2nm / numNuEP2nm;
651 
652  float ave3 = -1.;
653  if (numCREP2 > 0.) ave3 = aveCREP2 / numCREP2;
654 
655  if (fPrintLevel > 0) {
656  std::string nuType = "Other";
657  if (neutCCNC == simb::kCC) {
658  if (neutIntType == 1001) nuType = "CCQE";
659  if (neutIntType == 1091) nuType = "DIS";
660  if (neutIntType == 1097) nuType = "COH";
661  if (neutIntType > 1002 && neutIntType < 1091) nuType = "RES";
662  }
663  else if (neutCCNC == simb::kNC) {
664  nuType = "NC";
665  }
666  else {
667  nuType = "Unknown";
668  }
669  mf::LogVerbatim("PFPAna") << "EvtEP2 " << evt.id().event() << " NuType " << nuType << " Enu "
670  << std::fixed << std::setprecision(0) << neutEnergy << std::right
671  << std::fixed << std::setprecision(2) << " NuMuons " << ave1
672  << " NuPiKp " << ave2 << " CosmicRays " << ave3 << " CCNC "
673  << neutCCNC << " IntType " << neutIntType;
674  }
675  } // analyze
simb::MCParticle::E
double E(const int i=0) const
Definition: MCParticle.h:233
mf::LogVerbatim
MaybeLogger_< ELseverityLevel::ELsev_info, true > LogVerbatim
Definition: MessageLogger.h:216
pfpf::PFPAna::fElecKERange
std::vector< float > fElecKERange
Definition: PFPAna_module.cc:94
simb::MCParticle::PdgCode
int PdgCode() const
Definition: MCParticle.h:212
simb::MCNeutrino::CCNC
int CCNC() const
Definition: MCNeutrino.h:148
simb::MCTruth::GetNeutrino
const simb::MCNeutrino & GetNeutrino() const
Definition: MCTruth.h:77
geo::GeometryCore::Plane
PlaneGeo const & Plane(unsigned int const p, unsigned int const tpc=0, unsigned int const cstat=0) const
Returns the specified wire.
Definition: GeometryCore.h:2994
pfpf::PFPAna::fNuEP2_prot
TH1F * fNuEP2_prot
Definition: PFPAna_module.cc:67
mf::MaybeLogger_
Definition: MessageLogger.h:113
View
AdcChannelData::View View
Definition: test_AdcChannelData.cxx:23
pfpf::PFPAna::fVertexModuleLabel
std::string fVertexModuleLabel
Definition: PFPAna_module.cc:93
reco_momentum_tuples.float
float
Definition: reco_momentum_tuples.py:12
geo::View_t
enum geo::_plane_proj View_t
Enumerate the possible plane projections.
string
std::string string
Definition: nybbler.cc:12
recob::Hit::WireID
geo::WireID WireID() const
Definition: Hit.h:233
simb::MCParticle::Mother
int Mother() const
Definition: MCParticle.h:213
simb::MCNeutrino::Nu
const simb::MCParticle & Nu() const
Definition: MCNeutrino.h:146
simb::MCTruth::Origin
simb::Origin_t Origin() const
Definition: MCTruth.h:74
simb::MCParticle::Mass
double Mass() const
Definition: MCParticle.h:239
pfpf::PFPAna::fNuP_kaon
TH1F * fNuP_kaon
Definition: PFPAna_module.cc:76
pfpf::PFPAna::fHitsModuleLabel
std::string fHitsModuleLabel
Definition: PFPAna_module.cc:89
pfpf::PFPAna::fNuKE_prot
TH1F * fNuKE_prot
Definition: PFPAna_module.cc:62
pfpf::PFPAna::fProtKERange
std::vector< float > fProtKERange
Definition: PFPAna_module.cc:98
geo::WireID::Wire
WireID_t Wire
Index of the wire within its plane.
Definition: geo_types.h:580
const_iterator
intermediate_table::const_iterator const_iterator
Definition: intermediate_table.cc:28
genie::units::MeV
static constexpr double MeV
Definition: Units.h:129
mf::LogError
MaybeLogger_< ELseverityLevel::ELsev_error, false > LogError
Definition: MessageLogger.h:211
simb::MCParticle::Process
std::string Process() const
Definition: MCParticle.h:215
setprecision
Q_EXPORT QTSManip setprecision(int p)
Definition: qtextstream.h:343
pfpf::PFPAna::fClusterModuleLabel
std::string fClusterModuleLabel
Definition: PFPAna_module.cc:90
simb::MCParticle::TrackId
int TrackId() const
Definition: MCParticle.h:210
pfpf::PFPAna::fPionKERange
std::vector< float > fPionKERange
Definition: PFPAna_module.cc:96
pfpf::PFPAna::fNuEP2_KE_kaon
TProfile * fNuEP2_KE_kaon
Definition: PFPAna_module.cc:86
PDG
const uint PDG
Definition: qregexp.cpp:140
util::size
decltype(auto) constexpr size(T &&obj)
ADL-aware version of std::size.
Definition: StdUtils.h:92
pfpf::PFPAna::fNuVtx_dy
TH1F * fNuVtx_dy
Definition: PFPAna_module.cc:80
pfpf::PFPAna::fNuKE_elec
TH1F * fNuKE_elec
Definition: PFPAna_module.cc:58
simb::MCNeutrino::InteractionType
int InteractionType() const
Definition: MCNeutrino.h:150
abs
T abs(T value)
Definition: sparse_vector_test.cc:30
geo::PlaneGeo::View
View_t View() const
Which coordinate does this plane measure.
Definition: PlaneGeo.h:184
pfpf::PFPAna::fNuEP2_KE_elec
TProfile * fNuEP2_KE_elec
Definition: PFPAna_module.cc:83
art::DataViewImpl::getByLabel
bool getByLabel(std::string const &label, std::string const &instance, Handle< PROD > &result) const
Definition: DataViewImpl.h:633
geo::GeometryCore::Nplanes
unsigned int Nplanes(unsigned int tpc=0, unsigned int cstat=0) const
Returns the total number of wire planes in the specified TPC.
Definition: GeometryCore.h:2848
pfpf::PFPAna::fNuVtx_dz
TH1F * fNuVtx_dz
Definition: PFPAna_module.cc:81
genie::utils::res::PdgCode
int PdgCode(Resonance_t res, int Q)
(resonance id, charge) -> PDG code
Definition: BaryonResUtils.cxx:181
pfpf::PFPAna::fNuEP2_muon
TH1F * fNuEP2_muon
Definition: PFPAna_module.cc:64
pfpf::PFPAna::fNuE_kaon
TH1F * fNuE_kaon
Definition: PFPAna_module.cc:71
art::PtrVector::push_back
void push_back(Ptr< U > const &p)
Definition: PtrVector.h:435
pfpf::PFPAna::fNuP_muon
TH1F * fNuP_muon
Definition: PFPAna_module.cc:74
pfpf::PFPAna::fTrackWeightOption
short fTrackWeightOption
Definition: PFPAna_module.cc:99
cheat::ParticleInventoryService::TrackIdToMCTruth_P
const art::Ptr< simb::MCTruth > & TrackIdToMCTruth_P(int id) const
Definition: ParticleInventoryService_service.cc:142
submit_ppfxjobs.temp
temp
Definition: submit_ppfxjobs.py:110
pfpf::PFPAna::fNuEP2_KE_pion
TProfile * fNuEP2_KE_pion
Definition: PFPAna_module.cc:85
nhit
uint8_t nhit
Definition: CRTFragment.hh:201
pfpf::PFPAna::fNuEP2_kaon
TH1F * fNuEP2_kaon
Definition: PFPAna_module.cc:66
pfpf::PFPAna::fNuP_elec
TH1F * fNuP_elec
Definition: PFPAna_module.cc:73
pfpf::PFPAna::fPFParticleModuleLabel
std::string fPFParticleModuleLabel
Definition: PFPAna_module.cc:92
pfpf::PFPAna::fNHitInCluster
TH1F * fNHitInCluster
Definition: PFPAna_module.cc:52
cheat::BackTrackerService::TrackIdsToHits_Ps
std::vector< std::vector< art::Ptr< recob::Hit > > > TrackIdsToHits_Ps(detinfo::DetectorClocksData const &clockData, std::vector< int > const &tkIds, std::vector< art::Ptr< recob::Hit >> const &hitsIn) const
Definition: BackTrackerService_service.cc:243
pfpf::PFPAna::fNuKE_muon
TH1F * fNuKE_muon
Definition: PFPAna_module.cc:59
art::PtrVector::size
size_type size() const
Definition: PtrVector.h:302
pfpf::PFPAna::fNuP_prot
TH1F * fNuP_prot
Definition: PFPAna_module.cc:77
cheat::ParticleInventoryService::ParticleList
const sim::ParticleList & ParticleList() const
Definition: ParticleInventoryService_service.cc:114
pfpf::PFPAna::fNuVtx_dx
TH1F * fNuVtx_dx
Definition: PFPAna_module.cc:79
keras_to_tensorflow.int
int
Definition: keras_to_tensorflow.py:69
pfpf::PFPAna::fNuE_prot
TH1F * fNuE_prot
Definition: PFPAna_module.cc:72
pfpf::PFPAna::fNuEP2_KE_prot
TProfile * fNuEP2_KE_prot
Definition: PFPAna_module.cc:87
simb::MCParticle::Vx
double Vx(const int i=0) const
Definition: MCParticle.h:221
recob::Hit::PeakTime
float PeakTime() const
Time of the signal peak, in tick units.
Definition: Hit.h:218
pfpf::PFPAna::fNuE_muon
TH1F * fNuE_muon
Definition: PFPAna_module.cc:69
pfpf::PFPAna::fNClusters
TH1F * fNClusters
Definition: PFPAna_module.cc:51
pfpf::PFPAna::fNuEP2_KE_muon
TProfile * fNuEP2_KE_muon
Definition: PFPAna_module.cc:84
pfpf::PFPAna::fNuE_pion
TH1F * fNuE_pion
Definition: PFPAna_module.cc:70
pfpf::PFPAna::fNuE_elec
TH1F * fNuE_elec
Definition: PFPAna_module.cc:68
pfpf::PFPAna::fNuEP2_elec
TH1F * fNuEP2_elec
Definition: PFPAna_module.cc:63
pfpf::PFPAna::fNuKE_pion
TH1F * fNuKE_pion
Definition: PFPAna_module.cc:60
simb::MCParticle::Vz
double Vz(const int i=0) const
Definition: MCParticle.h:223
pfpf::PFPAna::fNuP_pion
TH1F * fNuP_pion
Definition: PFPAna_module.cc:75
run_cnn_1class.pdg
pdg
Definition: run_cnn_1class.py:54
art::EventID::event
EventNumber_t event() const
Definition: EventID.h:116
pfpf::PFPAna::fMuonKERange
std::vector< float > fMuonKERange
Definition: PFPAna_module.cc:95
pfpf::PFPAna::fNuKE_kaon
TH1F * fNuKE_kaon
Definition: PFPAna_module.cc:61
art::fill_ptr_vector
void fill_ptr_vector(std::vector< Ptr< T >> &ptrs, H const &h)
Definition: Ptr.h:297
util::quantities::second
second_as<> second
Type of time stored in seconds, in double precision.
Definition: spacetime.h:85
simb::MCNeutrino
Event generator information.
Definition: MCNeutrino.h:18
if
if(!yymsg) yymsg
art::Event::id
EventID id() const
Definition: Event.cc:34
simb::MCParticle::Vy
double Vy(const int i=0) const
Definition: MCParticle.h:222
pfpf::PFPAna::fKaonKERange
std::vector< float > fKaonKERange
Definition: PFPAna_module.cc:97
util::empty
decltype(auto) constexpr empty(T &&obj)
ADL-aware version of std::empty.
Definition: StdUtils.h:97
simb::kCosmicRay
Cosmic rays.
Definition: MCTruth.h:24
pfpf::PFPAna::fNuEP2_pion
TH1F * fNuEP2_pion
Definition: PFPAna_module.cc:65
simb::kBeamNeutrino
Beam neutrinos.
Definition: MCTruth.h:23
vertex
vertex reconstruction
Definition: AggregateVertexAna_module.cc:32
void pfpf::PFPAna::beginJob ( )
overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 127 of file PFPAna_module.cc.

128  {
129  art::ServiceHandle<art::TFileService const> tfs;
130 
131  fNClusters = tfs->make<TH1F>("fNoClustersInEvent", "Number of Clusters", 40, 0, 400);
132  fNHitInCluster = tfs->make<TH1F>("fNHitInCluster", "NHitInCluster", 100, 0, 100);
133 
134  if (!fSkipCosmics) {
135  // Cosmic ray histos
136  fCREP2 = tfs->make<TH1F>("CREP2", "CREP2", 50, 0, 1);
137  fCRE = tfs->make<TH1F>("CRE", "CR Efficiency", 50, 0, 1);
138  fCRP = tfs->make<TH1F>("CRP", "CR Purity", 50, 0, 1);
139  }
140  // Neutrino Int histos
141  fNuKE_elec = tfs->make<TH1F>("NuKE_elec", "NuKE electron", 100, 0, 4000);
142  fNuKE_muon = tfs->make<TH1F>("NuKE_muon", "NuKE muon", 100, 0, 4000);
143  fNuKE_pion = tfs->make<TH1F>("NuKE_pion", "NuKE pion", 100, 0, 4000);
144  fNuKE_kaon = tfs->make<TH1F>("NuKE_kaon", "NuKE kaon", 100, 0, 4000);
145  fNuKE_prot = tfs->make<TH1F>("NuKE_prot", "NuKE proton", 100, 0, 4000);
146 
147  fNuEP2_elec = tfs->make<TH1F>("NuEP2_elec", "NuEP2 electron", 50, 0, 1);
148  fNuEP2_muon = tfs->make<TH1F>("NuEP2_muon", "NuEP2 muon", 50, 0, 1);
149  fNuEP2_pion = tfs->make<TH1F>("NuEP2_pion", "NuEP2 pion", 50, 0, 1);
150  fNuEP2_kaon = tfs->make<TH1F>("NuEP2_kaon", "NuEP2 kaon", 50, 0, 1);
151  fNuEP2_prot = tfs->make<TH1F>("NuEP2_prot", "NuEP2 proton", 50, 0, 1);
152 
153  fNuE_elec = tfs->make<TH1F>("NuE_elec", "Nu Efficiency electron", 50, 0, 1);
154  fNuE_muon = tfs->make<TH1F>("NuE_muon", "Nu Efficiency muon", 50, 0, 1);
155  fNuE_pion = tfs->make<TH1F>("NuE_pion", "Nu Efficiency pion", 50, 0, 1);
156  fNuE_kaon = tfs->make<TH1F>("NuE_kaon", "Nu Efficiency kaon", 50, 0, 1);
157  fNuE_prot = tfs->make<TH1F>("NuE_prot", "Nu Efficiency proton", 50, 0, 1);
158 
159  fNuP_elec = tfs->make<TH1F>("NuP_elec", "Nu Purity electron", 50, 0, 1);
160  fNuP_muon = tfs->make<TH1F>("NuP_muon", "Nu Purity muon", 50, 0, 1);
161  fNuP_pion = tfs->make<TH1F>("NuP_pion", "Nu Purity pion", 50, 0, 1);
162  fNuP_kaon = tfs->make<TH1F>("NuP_kaon", "Nu Purity kaon", 50, 0, 1);
163  fNuP_prot = tfs->make<TH1F>("NuP_prot", "Nu Purity proton", 50, 0, 1);
164 
165  // True - Reco vertex difference
166  fNuVtx_dx = tfs->make<TH1F>("Vtx dx", "Vtx dx", 80, -10, 10);
167  fNuVtx_dy = tfs->make<TH1F>("Vtx dy", "Vtx dy", 80, -10, 10);
168  fNuVtx_dz = tfs->make<TH1F>("Vtx dz", "Vtx dz", 80, -10, 10);
169 
170  fNuEP2_KE_elec = tfs->make<TProfile>("NuEP2_KE_elec", "NuEP2 electron vs KE", 20, 0, 2000);
171  fNuEP2_KE_muon = tfs->make<TProfile>("NuEP2_KE_muon", "NuEP2 muon vs KE", 20, 0, 2000);
172  fNuEP2_KE_pion = tfs->make<TProfile>("NuEP2_KE_pion", "NuEP2 pion vs KE", 20, 0, 2000);
173  fNuEP2_KE_kaon = tfs->make<TProfile>("NuEP2_KE_kaon", "NuEP2 kaon vs KE", 20, 0, 2000);
174  fNuEP2_KE_prot = tfs->make<TProfile>("NuEP2_KE_prot", "NuEP2 proton vs KE", 20, 0, 2000);
175  }
pfpf::PFPAna::fNuEP2_prot
TH1F * fNuEP2_prot
Definition: PFPAna_module.cc:67
pfpf::PFPAna::fNuP_kaon
TH1F * fNuP_kaon
Definition: PFPAna_module.cc:76
pfpf::PFPAna::fNuKE_prot
TH1F * fNuKE_prot
Definition: PFPAna_module.cc:62
pfpf::PFPAna::fNuEP2_KE_kaon
TProfile * fNuEP2_KE_kaon
Definition: PFPAna_module.cc:86
pfpf::PFPAna::fNuVtx_dy
TH1F * fNuVtx_dy
Definition: PFPAna_module.cc:80
pfpf::PFPAna::fNuKE_elec
TH1F * fNuKE_elec
Definition: PFPAna_module.cc:58
pfpf::PFPAna::fNuEP2_KE_elec
TProfile * fNuEP2_KE_elec
Definition: PFPAna_module.cc:83
pfpf::PFPAna::fNuVtx_dz
TH1F * fNuVtx_dz
Definition: PFPAna_module.cc:81
pfpf::PFPAna::fNuEP2_muon
TH1F * fNuEP2_muon
Definition: PFPAna_module.cc:64
pfpf::PFPAna::fNuE_kaon
TH1F * fNuE_kaon
Definition: PFPAna_module.cc:71
pfpf::PFPAna::fNuP_muon
TH1F * fNuP_muon
Definition: PFPAna_module.cc:74
pfpf::PFPAna::fNuEP2_KE_pion
TProfile * fNuEP2_KE_pion
Definition: PFPAna_module.cc:85
pfpf::PFPAna::fNuEP2_kaon
TH1F * fNuEP2_kaon
Definition: PFPAna_module.cc:66
pfpf::PFPAna::fNuP_elec
TH1F * fNuP_elec
Definition: PFPAna_module.cc:73
pfpf::PFPAna::fNHitInCluster
TH1F * fNHitInCluster
Definition: PFPAna_module.cc:52
pfpf::PFPAna::fNuKE_muon
TH1F * fNuKE_muon
Definition: PFPAna_module.cc:59
pfpf::PFPAna::fNuP_prot
TH1F * fNuP_prot
Definition: PFPAna_module.cc:77
pfpf::PFPAna::fNuVtx_dx
TH1F * fNuVtx_dx
Definition: PFPAna_module.cc:79
pfpf::PFPAna::fNuE_prot
TH1F * fNuE_prot
Definition: PFPAna_module.cc:72
pfpf::PFPAna::fNuEP2_KE_prot
TProfile * fNuEP2_KE_prot
Definition: PFPAna_module.cc:87
pfpf::PFPAna::fNuE_muon
TH1F * fNuE_muon
Definition: PFPAna_module.cc:69
pfpf::PFPAna::fNClusters
TH1F * fNClusters
Definition: PFPAna_module.cc:51
pfpf::PFPAna::fNuEP2_KE_muon
TProfile * fNuEP2_KE_muon
Definition: PFPAna_module.cc:84
pfpf::PFPAna::fNuE_pion
TH1F * fNuE_pion
Definition: PFPAna_module.cc:70
pfpf::PFPAna::fNuE_elec
TH1F * fNuE_elec
Definition: PFPAna_module.cc:68
pfpf::PFPAna::fNuEP2_elec
TH1F * fNuEP2_elec
Definition: PFPAna_module.cc:63
pfpf::PFPAna::fNuKE_pion
TH1F * fNuKE_pion
Definition: PFPAna_module.cc:60
pfpf::PFPAna::fNuP_pion
TH1F * fNuP_pion
Definition: PFPAna_module.cc:75
pfpf::PFPAna::fNuKE_kaon
TH1F * fNuKE_kaon
Definition: PFPAna_module.cc:61
pfpf::PFPAna::fNuEP2_pion
TH1F * fNuEP2_pion
Definition: PFPAna_module.cc:65

Member Data Documentation

std::string pfpf::PFPAna::fClusterModuleLabel
private

Definition at line 90 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fCRE
private

Definition at line 55 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fCREP2
private

Definition at line 54 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fCRP
private

Definition at line 56 of file PFPAna_module.cc.

std::vector<float> pfpf::PFPAna::fElecKERange
private

Definition at line 94 of file PFPAna_module.cc.

std::string pfpf::PFPAna::fHitsModuleLabel
private

Definition at line 89 of file PFPAna_module.cc.

std::vector<float> pfpf::PFPAna::fKaonKERange
private

Definition at line 97 of file PFPAna_module.cc.

bool pfpf::PFPAna::fMergeDaughters
private

Definition at line 100 of file PFPAna_module.cc.

std::vector<float> pfpf::PFPAna::fMuonKERange
private

Definition at line 95 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNClusters
private

Definition at line 51 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNHitInCluster
private

Definition at line 52 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuE_elec
private

Definition at line 68 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuE_kaon
private

Definition at line 71 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuE_muon
private

Definition at line 69 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuE_pion
private

Definition at line 70 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuE_prot
private

Definition at line 72 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuEP2_elec
private

Definition at line 63 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuEP2_kaon
private

Definition at line 66 of file PFPAna_module.cc.

TProfile* pfpf::PFPAna::fNuEP2_KE_elec
private

Definition at line 83 of file PFPAna_module.cc.

TProfile* pfpf::PFPAna::fNuEP2_KE_kaon
private

Definition at line 86 of file PFPAna_module.cc.

TProfile* pfpf::PFPAna::fNuEP2_KE_muon
private

Definition at line 84 of file PFPAna_module.cc.

TProfile* pfpf::PFPAna::fNuEP2_KE_pion
private

Definition at line 85 of file PFPAna_module.cc.

TProfile* pfpf::PFPAna::fNuEP2_KE_prot
private

Definition at line 87 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuEP2_muon
private

Definition at line 64 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuEP2_pion
private

Definition at line 65 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuEP2_prot
private

Definition at line 67 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuKE_elec
private

Definition at line 58 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuKE_kaon
private

Definition at line 61 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuKE_muon
private

Definition at line 59 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuKE_pion
private

Definition at line 60 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuKE_prot
private

Definition at line 62 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuP_elec
private

Definition at line 73 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuP_kaon
private

Definition at line 76 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuP_muon
private

Definition at line 74 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuP_pion
private

Definition at line 75 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuP_prot
private

Definition at line 77 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuVtx_dx
private

Definition at line 79 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuVtx_dy
private

Definition at line 80 of file PFPAna_module.cc.

TH1F* pfpf::PFPAna::fNuVtx_dz
private

Definition at line 81 of file PFPAna_module.cc.

std::string pfpf::PFPAna::fPFParticleModuleLabel
private

Definition at line 92 of file PFPAna_module.cc.

std::vector<float> pfpf::PFPAna::fPionKERange
private

Definition at line 96 of file PFPAna_module.cc.

short pfpf::PFPAna::fPrintLevel
private

Definition at line 102 of file PFPAna_module.cc.

std::vector<float> pfpf::PFPAna::fProtKERange
private

Definition at line 98 of file PFPAna_module.cc.

bool pfpf::PFPAna::fSkipCosmics
private

Definition at line 101 of file PFPAna_module.cc.

std::string pfpf::PFPAna::fTrackModuleLabel
private

Definition at line 91 of file PFPAna_module.cc.

short pfpf::PFPAna::fTrackWeightOption
private

Definition at line 99 of file PFPAna_module.cc.

std::string pfpf::PFPAna::fVertexModuleLabel
private

Definition at line 93 of file PFPAna_module.cc.

The documentation for this class was generated from the following file:
Generated by   doxygen 1.8.11