Inheritance diagram for pfpf::PFPAna:

Public Types
using	ModuleType = EDAnalyzer

Public Member Functions
	PFPAna (fhicl::ParameterSet const &pset)

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)

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)

std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)

template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)

template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
std::string const &	processName () const

bool	wantAllEvents () const noexcept

bool	wantEvent (ScheduleID id, Event const &e) const

Handle< TriggerResults >	getTriggerResults (Event const &e) const

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 ()

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

Detailed Description

Definition at line 44 of file PFPAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 108 of file PFPAna_module.cc.

     : EDAnalyzer(pset)
     , fHitsModuleLabel(pset.get<std::string>("HitsModuleLabel"))
     , fClusterModuleLabel(pset.get<std::string>("ClusterModuleLabel"))
     , fTrackModuleLabel(pset.get<std::string>("TrackModuleLabel"))
     , fPFParticleModuleLabel(pset.get<std::string>("PFParticleModuleLabel"))
     , fVertexModuleLabel(pset.get<std::string>("VertexModuleLabel"))
     , fElecKERange(pset.get<std::vector<float>>("ElecKERange"))
     , fMuonKERange(pset.get<std::vector<float>>("MuonKERange"))
     , fPionKERange(pset.get<std::vector<float>>("PionKERange"))
     , fKaonKERange(pset.get<std::vector<float>>("KaonKERange"))
     , fProtKERange(pset.get<std::vector<float>>("ProtKERange"))
     , fTrackWeightOption(pset.get<short>("TrackWeightOption"))
     , fMergeDaughters(pset.get<bool>("MergeDaughters"))
     , fSkipCosmics(pset.get<bool>("SkipCosmics"))
     , fPrintLevel(pset.get<short>("PrintLevel"))
   {}

Member Function Documentation

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

overrideprivate

Definition at line 177 of file PFPAna_module.cc.

   {
     // code stolen from TrackAna_module.cc
     auto const& wireReadoutGeom = art::ServiceHandle<geo::WireReadout>()->Get();
     if (wireReadoutGeom.Nplanes() > 3) return;
 
     // get all hits in the event
     art::Handle<std::vector<recob::Hit>> hitListHandle;
     evt.getByLabel(fHitsModuleLabel, hitListHandle);
     std::vector<art::Ptr<recob::Hit>> allhits;
     art::fill_ptr_vector(allhits, hitListHandle);
     if (empty(allhits)) return;
 
     // get clusters and cluster-hit associations
     art::Handle<std::vector<recob::Cluster>> clusterListHandle;
     evt.getByLabel(fClusterModuleLabel, clusterListHandle);
     art::FindManyP<recob::Hit> fmh(clusterListHandle, evt, fClusterModuleLabel);
     if (clusterListHandle->size() == 0) return;
 
     // get 3D vertices
     art::Handle<std::vector<recob::Vertex>> vertexListHandle;
     evt.getByLabel(fVertexModuleLabel, vertexListHandle);
     art::PtrVector<recob::Vertex> recoVtxList;
     double xyz[3] = {0, 0, 0};
     for (unsigned int ii = 0; ii < vertexListHandle->size(); ++ii) {
       art::Ptr<recob::Vertex> vertex(vertexListHandle, ii);
       recoVtxList.push_back(vertex);
       vertex->XYZ(xyz);
     }
 
     // get PFParticles
     art::Handle<std::vector<recob::PFParticle>> PFPListHandle;
     evt.getByLabel(fPFParticleModuleLabel, PFPListHandle);
     art::PtrVector<recob::PFParticle> recoPFPList;
     for (unsigned int ii = 0; ii < PFPListHandle->size(); ++ii) {
       art::Ptr<recob::PFParticle> pfp(PFPListHandle, ii);
       recoPFPList.push_back(pfp);
       mf::LogVerbatim("PFPAna") << "PFParticle PDG " << pfp->PdgCode();
     }
 
     // list of all true particles
     art::ServiceHandle<cheat::BackTrackerService const> bt_serv;
     art::ServiceHandle<cheat::ParticleInventoryService const> pi_serv;
     sim::ParticleList const& plist = pi_serv->ParticleList();
     // list of all true particles that will be considered
     std::vector<const simb::MCParticle*> plist2;
     // true (reconstructed) hits for each particle in plist2
     std::vector<std::vector<art::Ptr<recob::Hit>>> hlist2;
     // index of cluster matched to each particle in plist2 in each plane
     std::vector<std::vector<short>> truToCl;
     // number of true hits in each plane and cluster
     std::vector<std::vector<unsigned short>> nTruHitInCl;
     //number of reconstructed hits in all clusters
     std::vector<unsigned short> nRecHitInCl;
 
     // calculate average EP2 for every event to facilitate code development
     // Beam Neutrinos - muons and not-muons
 
     float aveNuEP2mu = 0.;
     float numNuEP2mu = 0.;
     float aveNuEP2nm = 0.;
     float numNuEP2nm = 0.;
     // Cosmic Rays
     float aveCREP2 = 0.;
     float numCREP2 = 0.;
 
     // track ID of the neutrino
     int neutTrackID = -1;
     std::vector<int> tidlist;
     float neutEnergy = -1.;
     int neutIntType = -1;
     int neutCCNC = -1;
 
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
 
     for (sim::ParticleList::const_iterator ipart = plist.begin(); ipart != plist.end(); ++ipart) {
       const simb::MCParticle* part = (*ipart).second;
       assert(part != 0);
       int pdg = abs(part->PdgCode());
       int trackID = part->TrackId();
       art::Ptr<simb::MCTruth> theTruth = pi_serv->TrackIdToMCTruth_P(trackID);
       if (fSkipCosmics && theTruth->Origin() == simb::kCosmicRay) continue;
 
       if (fPrintLevel > 3)
         mf::LogVerbatim("PFPAna") << "Pre-Cuts origin " << theTruth->Origin() << " trackID "
                                   << trackID << " PDG " << part->PdgCode() << " E " << part->E()
                                   << " mass " << part->Mass() << " Mother "
                                   << part->Mother() + neutTrackID << " Proc " << part->Process();
 
       // Get the neutrino track ID. Assume that there is only one neutrino
       // interaction and it is first in the list of BeamNeutrino particles
       if (theTruth->Origin() == simb::kBeamNeutrino && neutTrackID < 0) {
         neutTrackID = trackID;
         simb::MCNeutrino theNeutrino = theTruth->GetNeutrino();
         neutEnergy = 1000. * theNeutrino.Nu().E();
         neutIntType = theNeutrino.InteractionType();
         neutCCNC = theNeutrino.CCNC();
         for (unsigned short iv = 0; iv < recoVtxList.size(); ++iv) {
           recoVtxList[iv]->XYZ(xyz);
           fNuVtx_dx->Fill(part->Vx() - xyz[0]);
           fNuVtx_dy->Fill(part->Vy() - xyz[1]);
           fNuVtx_dz->Fill(part->Vz() - xyz[2]);
         } // iv
       }   // theTruth->Origin() == simb::kBeamNeutrino && neutTrackID <
 
       bool isCharged = (pdg == 11) || (pdg == 13) || (pdg == 211) || (pdg == 321) || (pdg == 2212);
 
       if (!isCharged) continue;
 
       float KE = 1000 * (part->E() - part->Mass());
       // KE (MeV) cuts
       if (pdg == 11) {
         if (fElecKERange[0] < 0) continue;
         // only allow primary electrons
         if (part->Process() != "primary") continue;
         if (KE < fElecKERange[0] || KE > fElecKERange[1]) continue;
       }
       if (pdg == 13) {
         if (fMuonKERange[0] < 0) continue;
         if (KE < fMuonKERange[0] || KE > fMuonKERange[1]) continue;
       }
       if (pdg == 211) {
         if (fPionKERange[0] < 0) continue;
         if (KE < fPionKERange[0] || KE > fPionKERange[1]) continue;
       }
       if (pdg == 321) {
         if (fKaonKERange[0] < 0) continue;
         if (KE < fKaonKERange[0] || KE > fKaonKERange[1]) continue;
       }
       if (pdg == 2212) {
         if (fProtKERange[0] < 0) continue;
         if (KE < fProtKERange[0] || KE > fProtKERange[1]) continue;
       }
       // ignore secondaries from neutron interactions
       if (part->Process() == "NeutronInelastic") continue;
       plist2.push_back(part);
       tidlist.push_back(trackID);
       // initialize the true->(cluster,plane) association
       std::vector<short> temp{-1, -1, -1};
       truToCl.push_back(temp);
       // initialize the true hit count
       std::vector<unsigned short> temp2(3);
       nTruHitInCl.push_back(temp2);
 
       if (fPrintLevel > 2)
         mf::LogVerbatim("PFPAna") << plist2.size() - 1 << " Origin " << theTruth->Origin()
                                   << " trackID " << trackID << " PDG " << part->PdgCode() << " KE "
                                   << (int)KE << " Mother " << part->Mother() + neutTrackID
                                   << " Proc " << part->Process();
     }
 
     if (empty(plist2)) return;
 
     // get the hits (in all planes) that are matched to the true tracks
     hlist2 = bt_serv->TrackIdsToHits_Ps(clockData, tidlist, allhits);
     if (hlist2.size() != plist2.size()) {
       mf::LogError("PFPAna") << "MC particle list size " << plist2.size()
                              << " != size of MC particle true hits lists " << hlist2.size();
       return;
     }
     tidlist.clear();
 
     // vector of (mother, daughter) pairs
     std::vector<std::pair<unsigned short, unsigned short>> moda;
     // Deal with mother-daughter tracks
     if (fMergeDaughters && neutTrackID >= 0) {
       // Assume that daughters appear later in the list. Step backwards
       // to accumulate all generations of daughters
       for (unsigned short dpl = plist2.size() - 1; dpl > 0; --dpl) {
         // no mother
         if (plist2[dpl]->Mother() == 0) continue;
         // electron
         if (abs(plist2[dpl]->PdgCode()) == 11) continue;
         // the actual mother trackID is offset from the neutrino trackID
         int motherID = neutTrackID + plist2[dpl]->Mother() - 1;
         // ensure that we are only looking at BeamNeutrino daughters
         if (motherID < 0) continue;
         // count the number of daughters
         int ndtr = 0;
         for (unsigned short kpl = 0; kpl < plist2.size(); ++kpl) {
           if (plist2[kpl]->Mother() == motherID) ++ndtr;
         }
         // require only one daughter
         if (ndtr > 1) continue;
         // find the mother in the list
         int mpl = -1;
         for (unsigned short jpl = dpl - 1; jpl > 0; --jpl) {
           if (plist2[jpl]->TrackId() == motherID) {
             mpl = jpl;
             break;
           }
         } // jpl
         // mother not found for some reason
         if (mpl < 0) continue;
         // ensure that PDG code for mother and daughter are the same
         if (plist2[dpl]->PdgCode() != plist2[mpl]->PdgCode()) continue;
         moda.push_back(std::make_pair(mpl, dpl));
       } //  dpl
     }   // MergeDaughters
 
     // Now match reconstructed clusters to true particles.
     art::PtrVector<recob::Cluster> clusters;
     for (unsigned int ii = 0; ii < clusterListHandle->size(); ++ii) {
       art::Ptr<recob::Cluster> clusterHolder(clusterListHandle, ii);
       clusters.push_back(clusterHolder);
     }
 
     fNClusters->Fill(clusterListHandle->size());
     nRecHitInCl.resize(clusters.size());
 
     // get the plane from the view. Perhaps there is a method that does
     // this somewhere...
     std::map<geo::View_t, unsigned int> ViewToPlane;
     for (auto const& plane : wireReadoutGeom.Iterate<geo::PlaneGeo>(geo::TPCID{0, 0})) {
       geo::View_t view = plane.View();
       ViewToPlane[view] = plane.ID().Plane;
     }
     for (size_t icl = 0; icl < clusters.size(); ++icl) {
       unsigned int plane = ViewToPlane[clusters[icl]->View()];
       std::vector<art::Ptr<recob::Hit>> cluhits = fmh.at(icl);
       fNHitInCluster->Fill(cluhits.size());
       nRecHitInCl[icl] = cluhits.size();
       // count the number of hits matched to each true particle in plist2
       std::vector<unsigned short> nHitInPl2(plist2.size());
       for (size_t iht = 0; iht < cluhits.size(); ++iht) {
         // look for this hit in all of the truth hit lists
         short hitInPl2 = -1;
         for (unsigned short ipl = 0; ipl < plist2.size(); ++ipl) {
           unsigned short imat = 0;
           for (imat = 0; imat < hlist2[ipl].size(); ++imat) {
             if (cluhits[iht] == hlist2[ipl][imat]) break;
           } // imat
           if (imat < hlist2[ipl].size()) {
             hitInPl2 = ipl;
             break;
           }
         } // ipl
         if (hitInPl2 < 0) continue;
         // Assign the hit count to the mother if this is a daughter.
         // Mother-daughter pairs are entered in the moda vector in reverse
         // order, so assign daughter hits to the highest generation mother.
         for (unsigned short imd = 0; imd < moda.size(); ++imd) {
           if (moda[imd].second == hitInPl2) hitInPl2 = moda[imd].first;
         }
         // count
         ++nHitInPl2[hitInPl2];
       } // iht
       // Associate the cluster with the truth particle that has the highest
       // number of cluster hits
       unsigned short nhit = 0;
       short imtru = -1;
       for (unsigned int ipl = 0; ipl < nHitInPl2.size(); ++ipl) {
         if (nHitInPl2[ipl] > nhit) {
           nhit = nHitInPl2[ipl];
           imtru = ipl;
         }
       } // ipl
       // make the cluster->(true,plane) association and save the
       // number of true hits in the cluster
       if (imtru != -1) {
         // clobber a previously made association?
         if (nhit > nTruHitInCl[imtru][plane]) {
           truToCl[imtru][plane] = icl;
           nTruHitInCl[imtru][plane] = nhit;
         }
       } // imtru != 1
     }   // icl
 
     // ready to calculate Efficiency, Purity in each plane and EP2
     for (unsigned short ipl = 0; ipl < plist2.size(); ++ipl) {
       // ignore daughters
       bool skipit = false;
       for (unsigned short ii = 0; ii < moda.size(); ++ii) {
         if (moda[ii].second == ipl) {
           skipit = true;
           break;
         }
       } // ii
       if (skipit) continue;
       // ignore true particles with few true hits. Outside the detector
       // or not reconstructable?
       if (hlist2[ipl].size() < 3) continue;
 
       int trackID = plist2[ipl]->TrackId();
       art::Ptr<simb::MCTruth> theTruth = pi_serv->TrackIdToMCTruth_P(trackID);
       bool isCosmic = (theTruth->Origin() == simb::kCosmicRay);
       float KE = 1000 * (plist2[ipl]->E() - plist2[ipl]->Mass());
       int PDG = abs(plist2[ipl]->PdgCode());
 
       std::vector<short> nTru(wireReadoutGeom.Nplanes());
       std::vector<short> nRec(wireReadoutGeom.Nplanes());
       std::vector<short> nTruRec(wireReadoutGeom.Nplanes());
       std::vector<float> eff(wireReadoutGeom.Nplanes());
       std::vector<float> pur(wireReadoutGeom.Nplanes());
       std::vector<float> ep(wireReadoutGeom.Nplanes());
       for (unsigned int plane = 0; plane < wireReadoutGeom.Nplanes(); ++plane) {
         // count the number of true hits in this plane for the true particle.
         // First count the mother hits
         for (unsigned short ii = 0; ii < hlist2[ipl].size(); ++ii) {
           if (ViewToPlane[hlist2[ipl][ii]->View()] == plane) ++nTru[plane];
         } // ii
         // next look for daughters and count those hits in all generations
         unsigned short mom = ipl;
         std::vector<std::pair<unsigned short, unsigned short>>::reverse_iterator rit =
           moda.rbegin();
         while (rit != moda.rend()) {
           if ((*rit).first == mom) {
             unsigned short dau = (*rit).second;
             for (unsigned short jj = 0; jj < hlist2[dau].size(); ++jj) {
               if (ViewToPlane[hlist2[dau][jj]->View()] == plane) ++nTru[plane];
             } // jj
             // It is likely that one hit appears in the mother list
             // as well as the daughter list, so subtract one from the count
             --nTru[plane];
             mom = (*rit).second;
           } // (*rit).first == mom
           ++rit;
         } // rit
         if (nTru[plane] == 0) { continue; }
         short icl = truToCl[ipl][plane];
         nRec[plane] = nRecHitInCl[icl];
         nTruRec[plane] = nTruHitInCl[ipl][plane];
         if (nTru[plane] > 0) eff[plane] = (float)nTruRec[plane] / (float)nTru[plane];
         if (nRec[plane] > 0) pur[plane] = (float)nTruRec[plane] / (float)nRec[plane];
         ep[plane] = eff[plane] * pur[plane];
       } // plane
       // sort the ep values in ascending order
       std::vector<float> temp;
       temp = ep;
       std::sort(temp.begin(), temp.end());
       // EP2 is the second highest value
       unsigned short ii = temp.size() - 2;
       float ep2 = temp[ii];
       // find the plane that defined EP2
       short ep2Plane = 0;
       short ep2Cluster = 0;
       for (unsigned short jj = 0; jj < temp.size(); ++jj) {
         if (ep[jj] == ep2) {
           ep2Plane = jj;
           ep2Cluster = truToCl[ipl][ep2Plane];
           break;
         }
       } // jj
       // find the US and DS ends of the cluster for printing
       std::array<double, 2> clBeg, clEnd;
       if (ep2Cluster >= 0) {
         clBeg[0] = clusters[ep2Cluster]->StartWire();
         clBeg[1] = clusters[ep2Cluster]->StartTick();
         clEnd[0] = clusters[ep2Cluster]->EndWire();
         clEnd[1] = clusters[ep2Cluster]->EndTick();
       }
       else {
         clBeg.fill(0.);
         clEnd.fill(0.);
       }
       // fill histograms
       if (isCosmic) {
         fCREP2->Fill(ep2);
         fCRE->Fill(eff[ep2Plane]);
         fCRP->Fill(pur[ep2Plane]);
         aveCREP2 += ep2;
         numCREP2 += 1.;
         if (fPrintLevel > 1)
           mf::LogVerbatim("PFPAna")
             << ">>>CREP2 " << std::fixed << std::setprecision(2) << ep2 << " E " << eff[ep2Plane]
             << std::setprecision(2) << " P " << pur[ep2Plane] << " P:W:T " << ep2Plane << ":"
             << (int)clBeg[0] << ":" << (int)clBeg[1] << "-" << ep2Plane << ":" << (int)clEnd[0]
             << ":" << (int)clEnd[1] << " PDG " << PDG << " KE " << (int)KE << " MeV";
       } // isCosmic
       else {
         float wght = 1.;
         if (fTrackWeightOption == 1) wght = KE;
         // accumulate statistics for muons and not-muons
         if (PDG == 13) {
           aveNuEP2mu += ep2 * wght;
           numNuEP2mu += wght;
         }
         else {
           aveNuEP2nm += ep2 * wght;
           numNuEP2nm += wght;
         }
         if (PDG == 11) {
           fNuKE_elec->Fill(KE, wght);
           fNuE_elec->Fill(eff[ep2Plane], wght);
           fNuP_elec->Fill(pur[ep2Plane], wght);
           fNuEP2_elec->Fill(ep2, wght);
           fNuEP2_KE_elec->Fill(KE, ep2, wght);
         }
         else if (PDG == 13) {
           fNuKE_muon->Fill(KE, wght);
           fNuE_muon->Fill(eff[ep2Plane], wght);
           fNuP_muon->Fill(pur[ep2Plane], wght);
           fNuEP2_muon->Fill(ep2, wght);
           fNuEP2_KE_muon->Fill(KE, ep2, wght);
         }
         else if (PDG == 211) {
           fNuKE_pion->Fill(KE, wght);
           fNuE_pion->Fill(eff[ep2Plane], wght);
           fNuP_pion->Fill(pur[ep2Plane], wght);
           fNuEP2_pion->Fill(ep2, wght);
           fNuEP2_KE_pion->Fill(KE, ep2, wght);
         }
         else if (PDG == 321) {
           fNuKE_kaon->Fill(KE, wght);
           fNuE_kaon->Fill(eff[ep2Plane], wght);
           fNuP_kaon->Fill(pur[ep2Plane], wght);
           fNuEP2_kaon->Fill(ep2, wght);
           fNuEP2_KE_kaon->Fill(KE, ep2, wght);
         }
         else if (PDG == 2212) {
           fNuKE_prot->Fill(KE, wght);
           fNuE_prot->Fill(eff[ep2Plane], wght);
           fNuP_prot->Fill(pur[ep2Plane], wght);
           fNuEP2_prot->Fill(ep2, wght);
           fNuEP2_KE_prot->Fill(KE, ep2, wght);
         }
         if (fPrintLevel > 1)
           mf::LogVerbatim("PFPAna")
             << ">>>NuEP2 " << std::fixed << std::setprecision(2) << ep2 << " E " << eff[ep2Plane]
             << std::setprecision(2) << " P " << pur[ep2Plane] << " P:W:T " << ep2Plane << ":"
             << (int)clBeg[0] << ":" << (int)clBeg[1] << "-" << ep2Plane << ":" << (int)clEnd[0]
             << ":" << (int)clEnd[1] << " PDG " << PDG << " KE " << (int)KE << " MeV ";
         if (fPrintLevel > 2) {
           // print out the begin/end true hits
           mf::LogVerbatim mfp("PFPAna");
           mfp << " Truth P:W:T ";
           for (unsigned int plane = 0; plane < wireReadoutGeom.Nplanes(); ++plane) {
             unsigned short loW = 9999;
             unsigned short loT = 0;
             unsigned short hiW = 0;
             unsigned short hiT = 0;
             for (unsigned short ii = 0; ii < hlist2[ipl].size(); ++ii) {
               if (ViewToPlane[hlist2[ipl][ii]->View()] == plane) {
                 art::Ptr<recob::Hit> theHit = hlist2[ipl][ii];
                 if (theHit->WireID().Wire < loW) {
                   loW = theHit->WireID().Wire;
                   loT = theHit->PeakTime();
                 }
                 if (theHit->WireID().Wire > hiW) {
                   hiW = theHit->WireID().Wire;
                   hiT = theHit->PeakTime();
                 }
               } // correct view
             }   // ii
             mfp << plane << ":" << loW << ":" << loT << "-" << plane << ":" << hiW << ":" << hiT
                 << " ";
           } // plane
         }   // fPrintLevel > 2
       }     // !isCosmic
     }       // ipl
 
     float ave1 = -1.;
     if (numNuEP2mu > 0.) ave1 = aveNuEP2mu / numNuEP2mu;
 
     float ave2 = -1.;
     if (numNuEP2nm > 0.) ave2 = aveNuEP2nm / numNuEP2nm;
 
     float ave3 = -1.;
     if (numCREP2 > 0.) ave3 = aveCREP2 / numCREP2;
 
     if (fPrintLevel > 0) {
       std::string nuType = "Other";
       if (neutCCNC == simb::kCC) {
         if (neutIntType == 1001) nuType = "CCQE";
         if (neutIntType == 1091) nuType = "DIS";
         if (neutIntType == 1097) nuType = "COH";
         if (neutIntType > 1002 && neutIntType < 1091) nuType = "RES";
       }
       else if (neutCCNC == simb::kNC) {
         nuType = "NC";
       }
       else {
         nuType = "Unknown";
       }
       mf::LogVerbatim("PFPAna") << "EvtEP2 " << evt.id().event() << " NuType " << nuType << " Enu "
                                 << std::fixed << std::setprecision(0) << neutEnergy << std::right
                                 << std::fixed << std::setprecision(2) << " NuMuons " << ave1
                                 << " NuPiKp " << ave2 << " CosmicRays " << ave3 << " CCNC "
                                 << neutCCNC << " IntType " << neutIntType;
     }
   } // analyze

void pfpf::PFPAna::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 127 of file PFPAna_module.cc.

References fCRE, fCREP2, fCRP, fNClusters, fNHitInCluster, fNuE_elec, fNuE_kaon, fNuE_muon, fNuE_pion, fNuE_prot, fNuEP2_elec, fNuEP2_kaon, fNuEP2_KE_elec, fNuEP2_KE_kaon, fNuEP2_KE_muon, fNuEP2_KE_pion, fNuEP2_KE_prot, fNuEP2_muon, fNuEP2_pion, fNuEP2_prot, fNuKE_elec, fNuKE_kaon, fNuKE_muon, fNuKE_pion, fNuKE_prot, fNuP_elec, fNuP_kaon, fNuP_muon, fNuP_pion, fNuP_prot, fNuVtx_dx, fNuVtx_dy, fNuVtx_dz, and fSkipCosmics.

   {
     art::ServiceHandle<art::TFileService const> tfs;
 
     fNClusters = tfs->make<TH1F>("fNoClustersInEvent", "Number of Clusters", 40, 0, 400);
     fNHitInCluster = tfs->make<TH1F>("fNHitInCluster", "NHitInCluster", 100, 0, 100);
 
     if (!fSkipCosmics) {
       // Cosmic ray histos
       fCREP2 = tfs->make<TH1F>("CREP2", "CREP2", 50, 0, 1);
       fCRE = tfs->make<TH1F>("CRE", "CR Efficiency", 50, 0, 1);
       fCRP = tfs->make<TH1F>("CRP", "CR Purity", 50, 0, 1);
     }
     // Neutrino Int histos
     fNuKE_elec = tfs->make<TH1F>("NuKE_elec", "NuKE electron", 100, 0, 4000);
     fNuKE_muon = tfs->make<TH1F>("NuKE_muon", "NuKE muon", 100, 0, 4000);
     fNuKE_pion = tfs->make<TH1F>("NuKE_pion", "NuKE pion", 100, 0, 4000);
     fNuKE_kaon = tfs->make<TH1F>("NuKE_kaon", "NuKE kaon", 100, 0, 4000);
     fNuKE_prot = tfs->make<TH1F>("NuKE_prot", "NuKE proton", 100, 0, 4000);
 
     fNuEP2_elec = tfs->make<TH1F>("NuEP2_elec", "NuEP2 electron", 50, 0, 1);
     fNuEP2_muon = tfs->make<TH1F>("NuEP2_muon", "NuEP2 muon", 50, 0, 1);
     fNuEP2_pion = tfs->make<TH1F>("NuEP2_pion", "NuEP2 pion", 50, 0, 1);
     fNuEP2_kaon = tfs->make<TH1F>("NuEP2_kaon", "NuEP2 kaon", 50, 0, 1);
     fNuEP2_prot = tfs->make<TH1F>("NuEP2_prot", "NuEP2 proton", 50, 0, 1);
 
     fNuE_elec = tfs->make<TH1F>("NuE_elec", "Nu Efficiency electron", 50, 0, 1);
     fNuE_muon = tfs->make<TH1F>("NuE_muon", "Nu Efficiency muon", 50, 0, 1);
     fNuE_pion = tfs->make<TH1F>("NuE_pion", "Nu Efficiency pion", 50, 0, 1);
     fNuE_kaon = tfs->make<TH1F>("NuE_kaon", "Nu Efficiency kaon", 50, 0, 1);
     fNuE_prot = tfs->make<TH1F>("NuE_prot", "Nu Efficiency proton", 50, 0, 1);
 
     fNuP_elec = tfs->make<TH1F>("NuP_elec", "Nu Purity electron", 50, 0, 1);
     fNuP_muon = tfs->make<TH1F>("NuP_muon", "Nu Purity muon", 50, 0, 1);
     fNuP_pion = tfs->make<TH1F>("NuP_pion", "Nu Purity pion", 50, 0, 1);
     fNuP_kaon = tfs->make<TH1F>("NuP_kaon", "Nu Purity kaon", 50, 0, 1);
     fNuP_prot = tfs->make<TH1F>("NuP_prot", "Nu Purity proton", 50, 0, 1);
 
     // True - Reco vertex difference
     fNuVtx_dx = tfs->make<TH1F>("Vtx dx", "Vtx dx", 80, -10, 10);
     fNuVtx_dy = tfs->make<TH1F>("Vtx dy", "Vtx dy", 80, -10, 10);
     fNuVtx_dz = tfs->make<TH1F>("Vtx dz", "Vtx dz", 80, -10, 10);
 
     fNuEP2_KE_elec = tfs->make<TProfile>("NuEP2_KE_elec", "NuEP2 electron vs KE", 20, 0, 2000);
     fNuEP2_KE_muon = tfs->make<TProfile>("NuEP2_KE_muon", "NuEP2 muon vs KE", 20, 0, 2000);
     fNuEP2_KE_pion = tfs->make<TProfile>("NuEP2_KE_pion", "NuEP2 pion vs KE", 20, 0, 2000);
     fNuEP2_KE_kaon = tfs->make<TProfile>("NuEP2_KE_kaon", "NuEP2 kaon vs KE", 20, 0, 2000);
     fNuEP2_KE_prot = tfs->make<TProfile>("NuEP2_KE_prot", "NuEP2 proton vs KE", 20, 0, 2000);
   }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::consumes ( InputTag const & tag )

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

   {
     return collector_.consumes<T, BT>(tag);
   }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

   {
     return collector_;
   }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

   {
     collector_.consumesMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::consumesView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const & tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

   {
     return collector_.consumesView<T, BT>(tag);
   }

void art::detail::Analyzer::doBeginJob ( SharedResources const & resources )

inherited

Definition at line 25 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     setupQueues(resources);
     ProcessingFrame const frame{ScheduleID{}};
     beginJobWithFrame(frame);
   }

bool art::detail::Analyzer::doBeginRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 68 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::doBeginSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 84 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     beginSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

void art::detail::Analyzer::doEndJob ( )

inherited

Definition at line 33 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     endJobWithFrame(frame);
   }

bool art::detail::Analyzer::doEndRun	(	RunPrincipal &	rp,
		ModuleContext const &	mc
	)

inherited

Definition at line 76 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endRunWithFrame(std::as_const(rp).makeRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::doEndSubRun	(	SubRunPrincipal &	srp,
		ModuleContext const &	mc
	)

inherited

Definition at line 92 of file Analyzer.cc.

References art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{mc.scheduleID()};
     endSubRunWithFrame(std::as_const(srp).makeSubRun(mc), frame);
     return true;
   }

bool art::detail::Analyzer::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
	)

inherited

Definition at line 100 of file Analyzer.cc.

References e, and art::ModuleContext::scheduleID().

Referenced by art::detail::Analyzer::Analyzer().

   {
     auto const e = std::as_const(ep).makeEvent(mc);
     if (wantEvent(mc.scheduleID(), e)) {
       ++counts_run;
       ProcessingFrame const frame{mc.scheduleID()};
       analyzeWithFrame(e, frame);
       ++counts_passed;
     }
     return true;
   }

void art::detail::Analyzer::doRespondToCloseInputFile ( FileBlock const & fb )

inherited

Definition at line 47 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseInputFileWithFrame(fb, frame);
   }

void art::detail::Analyzer::doRespondToCloseOutputFiles ( FileBlock const & fb )

inherited

Definition at line 61 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToCloseOutputFilesWithFrame(fb, frame);
   }

void art::detail::Analyzer::doRespondToOpenInputFile ( FileBlock const & fb )

inherited

Definition at line 40 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenInputFileWithFrame(fb, frame);
   }

void art::detail::Analyzer::doRespondToOpenOutputFiles ( FileBlock const & fb )

inherited

Definition at line 54 of file Analyzer.cc.

Referenced by art::detail::Analyzer::Analyzer().

   {
     ProcessingFrame const frame{ScheduleID{}};
     respondToOpenOutputFilesWithFrame(fb, frame);
   }

std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

   {
     return collector_.getConsumables();
   }

Handle< TriggerResults > art::Observer::getTriggerResults ( Event const & e ) const

protectedinherited

Definition at line 75 of file Observer.cc.

References art::ProductRetriever::get(), and art::Observer::selectors_.

Referenced by art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (selectors_) {
       return selectors_->getOneTriggerResults(e);
     }
 
     // The following applies for cases where no SelectEvents entries
     // exist.
     Handle<TriggerResults> h;
     if (e.get(empty_process_name, h)) {
       return h;
     }
     return Handle<TriggerResults>{};
   }

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const & wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

   {
     return doMakeWorker(wp);
   }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::mayConsume ( InputTag const & tag )

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

   {
     return collector_.mayConsume<T, BT>(tag);
   }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

   {
     collector_.mayConsumeMany<T, BT>();
   }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const & )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const & tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

   {
     return collector_.mayConsumeView<T, BT>(tag);
   }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

   {
     if (md_.has_value()) {
       return *md_;
     }
 
     throw Exception{errors::LogicError,
                     "There was an error while calling moduleDescription().\n"}
       << "The moduleDescription() base-class member function cannot be called\n"
          "during module construction.  To determine which module is "
          "responsible\n"
          "for calling it, find the '<module type>:<module "
          "label>@Construction'\n"
          "tag in the message prefix above.  Please contact artists@fnal.gov\n"
          "for guidance.\n";
   }

string const & art::Observer::processName ( ) const

protectedinherited

Definition at line 57 of file Observer.cc.

References art::Observer::process_name_.

Referenced by art::FileDumperOutput::printPrincipal().

   {
     return process_name_;
   }

void art::ModuleBase::setModuleDescription ( ModuleDescription const & md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

   {
     md_ = md;
   }

void art::ModuleBase::sortConsumables ( std::string const & current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

   {
     // Now that we know we have seen all the consumes declarations,
     // sort the results for fast lookup later.
     collector_.sortConsumables(current_process_name);
   }

bool art::Observer::wantAllEvents ( ) const

inlineprotectednoexceptinherited

Definition at line 31 of file Observer.h.

References e, art::Observer::getTriggerResults(), art::Observer::wantAllEvents_, and art::Observer::wantEvent().

     {
       return wantAllEvents_;
     }

bool art::Observer::wantEvent	(	ScheduleID	id,
		Event const &	e
	)		const

protectedinherited

Definition at line 63 of file Observer.cc.

References art::Observer::rejectors_, art::Observer::selectors_, and art::Observer::wantAllEvents_.

Referenced by art::OutputModule::doEvent(), art::OutputModule::doWriteEvent(), and art::Observer::wantAllEvents().

   {
     if (wantAllEvents_) {
       return true;
     }
     bool const select_event = selectors_ ? selectors_->matchEvent(id, e) : true;
     bool const reject_event =
       rejectors_ ? rejectors_->matchEvent(id, e) : false;
     return select_event and not reject_event;
   }

Member Data Documentation

std::string pfpf::PFPAna::fClusterModuleLabel

private

Definition at line 91 of file PFPAna_module.cc.

Referenced by analyze().

TH1F* pfpf::PFPAna::fCRE

private

Definition at line 56 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fCREP2

private

Definition at line 55 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fCRP

private

Definition at line 57 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

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

private

Definition at line 95 of file PFPAna_module.cc.

Referenced by analyze().

std::string pfpf::PFPAna::fHitsModuleLabel

private

Definition at line 90 of file PFPAna_module.cc.

Referenced by analyze().

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

private

Definition at line 98 of file PFPAna_module.cc.

Referenced by analyze().

bool pfpf::PFPAna::fMergeDaughters

private

Definition at line 101 of file PFPAna_module.cc.

Referenced by analyze().

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

private

Definition at line 96 of file PFPAna_module.cc.

Referenced by analyze().

TH1F* pfpf::PFPAna::fNClusters

private

Definition at line 52 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNHitInCluster

private

Definition at line 53 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuE_elec

private

Definition at line 69 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuE_kaon

private

Definition at line 72 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuE_muon

private

Definition at line 70 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuE_pion

private

Definition at line 71 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuE_prot

private

Definition at line 73 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuEP2_elec

private

Definition at line 64 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuEP2_kaon

private

Definition at line 67 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* pfpf::PFPAna::fNuEP2_KE_elec

private

Definition at line 84 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* pfpf::PFPAna::fNuEP2_KE_kaon

private

Definition at line 87 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* pfpf::PFPAna::fNuEP2_KE_muon

private

Definition at line 85 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* pfpf::PFPAna::fNuEP2_KE_pion

private

Definition at line 86 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* pfpf::PFPAna::fNuEP2_KE_prot

private

Definition at line 88 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuEP2_muon

private

Definition at line 65 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuEP2_pion

private

Definition at line 66 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuEP2_prot

private

Definition at line 68 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuKE_elec

private

Definition at line 59 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuKE_kaon

private

Definition at line 62 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuKE_muon

private

Definition at line 60 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuKE_pion

private

Definition at line 61 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuKE_prot

private

Definition at line 63 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuP_elec

private

Definition at line 74 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuP_kaon

private

Definition at line 77 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuP_muon

private

Definition at line 75 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuP_pion

private

Definition at line 76 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuP_prot

private

Definition at line 78 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuVtx_dx

private

Definition at line 80 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuVtx_dy

private

Definition at line 81 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* pfpf::PFPAna::fNuVtx_dz

private

Definition at line 82 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

std::string pfpf::PFPAna::fPFParticleModuleLabel

private

Definition at line 93 of file PFPAna_module.cc.

Referenced by analyze().

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

private

Definition at line 97 of file PFPAna_module.cc.

Referenced by analyze().

short pfpf::PFPAna::fPrintLevel

private

Definition at line 103 of file PFPAna_module.cc.

Referenced by analyze().

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

private

Definition at line 99 of file PFPAna_module.cc.

Referenced by analyze().

bool pfpf::PFPAna::fSkipCosmics

private

Definition at line 102 of file PFPAna_module.cc.

Referenced by analyze(), and beginJob().

std::string pfpf::PFPAna::fTrackModuleLabel

private

Definition at line 92 of file PFPAna_module.cc.

short pfpf::PFPAna::fTrackWeightOption

private

Definition at line 100 of file PFPAna_module.cc.

Referenced by analyze().

std::string pfpf::PFPAna::fVertexModuleLabel

private

Definition at line 94 of file PFPAna_module.cc.

Referenced by analyze().

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

larreco/v10_01_07/source/larreco/TrackFinder/PFPAna_module.cc

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation