cluster::ClusterAna Class Reference

Inheritance diagram for cluster::ClusterAna:

Public Types
using	ModuleType = EDAnalyzer

Public Member Functions
	ClusterAna (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
void	endJob () 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	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
short	moduleID
std::ofstream	outFile

Detailed Description

Definition at line 44 of file ClusterAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

cluster::ClusterAna::ClusterAna ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 116 of file ClusterAna_module.cc.

References fClusterModuleLabel, fPrintLevel, moduleID, and outFile.

    : EDAnalyzer(pset)
    , fHitsModuleLabel(pset.get<std::string>("HitsModuleLabel"))
    , fClusterModuleLabel(pset.get<std::string>("ClusterModuleLabel"))
    , 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"))
  {

    if (fPrintLevel == -1) {
      // encode the clustermodule label into an integer
      moduleID = 0;
      size_t found = fClusterModuleLabel.find("traj");
      if (found != std::string::npos) moduleID = 1;
      found = fClusterModuleLabel.find("line");
      if (found != std::string::npos) moduleID = 2;
      found = fClusterModuleLabel.find("fuzz");
      if (found != std::string::npos) moduleID = 3;
      found = fClusterModuleLabel.find("pand");
      if (found != std::string::npos) moduleID = 4;
      std::cout << "fClusterModuleLabel " << fClusterModuleLabel << " ID " << moduleID << "\n";

      std::string fileName = fClusterModuleLabel + ".tru";
      outFile.open(fileName);
    }
  }

Member Function Documentation

void cluster::ClusterAna::analyze ( const art::Event &  evt )

overrideprivate

Definition at line 208 of file ClusterAna_module.cc.

References util::abs(), sim::ParticleList::begin(), simb::MCNeutrino::CCNC(), DEFINE_ART_MODULE, simb::MCParticle::E(), sim::ParticleList::end(), art::EventID::event(), fClusterModuleLabel, fCRE, fCREP2, fCRP, fElecKERange, fHitsModuleLabel, art::fill_ptr_vector(), fKaonKERange, fMergeDaughters, fMuonKERange, 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, fPionKERange, fPrintLevel, fProtKERange, fSkipCosmics, fTrackWeightOption, fVertexModuleLabel, art::ProductRetriever::getByLabel(), simb::MCTruth::GetNeutrino(), art::Event::id(), simb::MCNeutrino::InteractionType(), ipart, art::Handle< T >::isValid(), geo::GeometryCore::Iterate(), simb::kBeamNeutrino, simb::kCC, simb::kCosmicRay, simb::kNC, simb::kSingleParticle, simb::MCParticle::Mass(), mfp, moduleID, simb::MCParticle::Mother(), geo::GeometryCore::Nplanes(), simb::MCNeutrino::Nu(), simb::MCTruth::Origin(), outFile, part, cheat::ParticleInventoryService::ParticleList(), simb::MCParticle::PdgCode(), recob::Hit::PeakTime(), simb::MCParticle::Process(), art::PtrVector< T >::push_back(), art::right(), art::PtrVector< T >::size(), util::size(), simb::MCParticle::TrackId(), cheat::BackTrackerService::TrackIdsToHits_Ps(), cheat::ParticleInventoryService::TrackIdToMCTruth_P(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), simb::MCParticle::Vz(), geo::WireID::Wire, recob::Hit::WireID(), and recob::Vertex::XYZ().

  {
    // code stolen from TrackAna_module.cc
    art::ServiceHandle<geo::Geometry const> geom;
    if (geom->Nplanes() > 3) {
      mf::LogError("ClusterAna") << "Too many planes for this code...";
      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 (allhits.size() == 0) 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) {
      std::cout << "No clusters in event. Hits size " << allhits.size() << "\n";
      // don't return or the statistics will be off
    }

    // get 3D vertices
    art::Handle<std::vector<recob::Vertex>> vertexListHandle;
    double xyz[3] = {0, 0, 0};
    art::PtrVector<recob::Vertex> recoVtxList;
    if (vertexListHandle.isValid()) {
      evt.getByLabel(fVertexModuleLabel, vertexListHandle);
      for (unsigned int ii = 0; ii < vertexListHandle->size(); ++ii) {
        art::Ptr<recob::Vertex> vertex(vertexListHandle, ii);
        recoVtxList.push_back(vertex);
        vertex->XYZ(xyz);
      } // ii
    }   // vertexListHandle

    // 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 (or single particle)
    int neutTrackID = -1;
    std::vector<int> tidlist;
    float neutEnergy = -1.;
    int neutIntType = -1;
    int neutCCNC = -1;
    bool isNeutrino = false;
    bool isSingleParticle = false;

    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 (theTruth->Origin() == simb::kBeamNeutrino) isNeutrino = true;
      if (theTruth->Origin() == simb::kSingleParticle) isSingleParticle = true;
      float KE = 1000 * (part->E() - part->Mass());

      // Get the neutrino track ID. Assume that there is only one neutrino
      // interaction and it is first in the list of BeamNeutrino particles
      if ((isNeutrino || isSingleParticle) && neutTrackID < 0) {
        neutTrackID = trackID;
        if (isNeutrino) {
          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 <

      if (fPrintLevel > 3)
        mf::LogVerbatim("ClusterAna") << "Origin " << theTruth->Origin() << " trackID " << trackID
                                      << " PDG " << part->PdgCode() << " KE " << (int)KE << " MeV "
                                      << " neutTrackID " << neutTrackID << " Mother "
                                      << part->Mother() << " Proc " << part->Process();

      bool isCharged = (pdg == 11) || (pdg == 13) || (pdg == 211) || (pdg == 321) || (pdg == 2212);

      if (!isCharged) continue;

      // 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("ClusterAna")
          << "plist2[" << plist2.size() - 1 << "]"
          << " Origin " << theTruth->Origin() << " trackID " << trackID << " PDG "
          << part->PdgCode() << " KE " << (int)KE << " Mother " << part->Mother() + neutTrackID - 1
          << " Proc " << part->Process();
    }

    if (plist2.size() == 0) 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("ClusterAna") << "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
      unsigned short ii, dpl, jj, jpl, kpl;
      for (ii = 0; ii < plist2.size(); ++ii) {
        dpl = plist2.size() - 1 - ii;
        // 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 or single particle daughters
        if (motherID != neutTrackID) continue;
        // count the number of daughters
        int ndtr = 0;
        for (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 (jj = 0; jj < plist2.size(); ++jj) {
          jpl = plist2.size() - 1 - jj;
          if (plist2[jpl]->TrackId() == motherID) {
            mpl = jpl;
            break;
          }
        } // jpl
        // mother not found for some reason
        if (mpl < 0) {
          mf::LogVerbatim("ClusterAna")
            << " mother of daughter " << dpl << " not found. mpl = " << mpl;
          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 : geom->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(geom->Nplanes());
      std::vector<short> nRec(geom->Nplanes());
      std::vector<short> nTruRec(geom->Nplanes());
      std::vector<float> eff(geom->Nplanes());
      std::vector<float> pur(geom->Nplanes());
      std::vector<float> ep(geom->Nplanes());
      for (unsigned int plane = 0; plane < geom->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
            mom = (*rit).second;
          } // (*rit).first == mom
          ++rit;
        } // rit

        if (nTru[plane] == 0) continue;
        // Ignore short truth clusters
        if (nTru[plane] < 3) continue;
        short icl = truToCl[ipl][plane];
        if (icl < 0) { nRec[plane] = 0; }
        else {
          nRec[plane] = nRecHitInCl[icl];
        }
        nTruRec[plane] = nTruHitInCl[ipl][plane];
        eff[plane] = 0;
        pur[plane] = 0;
        if (nTru[plane] > 0) eff[plane] = (float)nTruRec[plane] / (float)nTru[plane];
        if (nRec[plane] > 0) pur[plane] = (float)nTruRec[plane] / (float)nRec[plane];
        // do this to prevent histogram under/over flow
        if (eff[plane] == 0) { eff[plane] = 0.001; }
        if (pur[plane] == 0) { pur[plane] = 0.001; }
        if (eff[plane] >= 1) { eff[plane] = 0.999; }
        if (pur[plane] >= 1) { pur[plane] = 0.999; }
        ep[plane] = eff[plane] * pur[plane];
        if (fPrintLevel == -1)
          outFile << moduleID << " " << evt.id().event() << " " << trackID << " " << PDG << " "
                  << (int)KE << " " << plane << " " << nTru[plane] << " " << std::setprecision(3)
                  << eff[plane] << " " << pur[plane] << "\n";
      } // 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("ClusterAna")
            << ">>>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("ClusterAna")
            << ">>>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("ClusterAna");
          mfp << " Truth P:W:T ";
          for (unsigned int plane = 0; plane < geom->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("ClusterAna")
        << "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 cluster::ClusterAna::beginJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 151 of file ClusterAna_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.

  {

    // get access to the TFile service
    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", 200, 0, 2000);
    fNuEP2_KE_muon = tfs->make<TProfile>("NuEP2_KE_muon", "NuEP2 muon vs KE", 200, 0, 2000);
    fNuEP2_KE_pion = tfs->make<TProfile>("NuEP2_KE_pion", "NuEP2 pion vs KE", 200, 0, 2000);
    fNuEP2_KE_kaon = tfs->make<TProfile>("NuEP2_KE_kaon", "NuEP2 kaon vs KE", 200, 0, 2000);
    fNuEP2_KE_prot = tfs->make<TProfile>("NuEP2_KE_prot", "NuEP2 proton vs KE", 200, 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);
  }

void cluster::ClusterAna::endJob ( )

overrideprivatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 203 of file ClusterAna_module.cc.

References fPrintLevel, and outFile.

  {
    if (fPrintLevel == -1) outFile.close();
  }

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 cluster::ClusterAna::fClusterModuleLabel

private

Definition at line 93 of file ClusterAna_module.cc.

Referenced by analyze(), and ClusterAna().

TH1F* cluster::ClusterAna::fCRE

private

Definition at line 57 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fCREP2

private

Definition at line 56 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fCRP

private

Definition at line 58 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

std::vector<float> cluster::ClusterAna::fElecKERange

private

Definition at line 95 of file ClusterAna_module.cc.

Referenced by analyze().

std::string cluster::ClusterAna::fHitsModuleLabel

private

Definition at line 92 of file ClusterAna_module.cc.

Referenced by analyze().

std::vector<float> cluster::ClusterAna::fKaonKERange

private

Definition at line 98 of file ClusterAna_module.cc.

Referenced by analyze().

bool cluster::ClusterAna::fMergeDaughters

private

Definition at line 101 of file ClusterAna_module.cc.

Referenced by analyze().

std::vector<float> cluster::ClusterAna::fMuonKERange

private

Definition at line 96 of file ClusterAna_module.cc.

Referenced by analyze().

TH1F* cluster::ClusterAna::fNClusters

private

Definition at line 53 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNHitInCluster

private

Definition at line 54 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuE_elec

private

Definition at line 71 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuE_kaon

private

Definition at line 74 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuE_muon

private

Definition at line 72 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuE_pion

private

Definition at line 73 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuE_prot

private

Definition at line 75 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuEP2_elec

private

Definition at line 66 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuEP2_kaon

private

Definition at line 69 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* cluster::ClusterAna::fNuEP2_KE_elec

private

Definition at line 86 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* cluster::ClusterAna::fNuEP2_KE_kaon

private

Definition at line 89 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* cluster::ClusterAna::fNuEP2_KE_muon

private

Definition at line 87 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* cluster::ClusterAna::fNuEP2_KE_pion

private

Definition at line 88 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TProfile* cluster::ClusterAna::fNuEP2_KE_prot

private

Definition at line 90 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuEP2_muon

private

Definition at line 67 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuEP2_pion

private

Definition at line 68 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuEP2_prot

private

Definition at line 70 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuKE_elec

private

Definition at line 60 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuKE_kaon

private

Definition at line 63 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuKE_muon

private

Definition at line 61 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuKE_pion

private

Definition at line 62 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuKE_prot

private

Definition at line 64 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuP_elec

private

Definition at line 76 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuP_kaon

private

Definition at line 79 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuP_muon

private

Definition at line 77 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuP_pion

private

Definition at line 78 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuP_prot

private

Definition at line 80 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuVtx_dx

private

Definition at line 82 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuVtx_dy

private

Definition at line 83 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::ClusterAna::fNuVtx_dz

private

Definition at line 84 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

std::vector<float> cluster::ClusterAna::fPionKERange

private

Definition at line 97 of file ClusterAna_module.cc.

Referenced by analyze().

short cluster::ClusterAna::fPrintLevel

private

Definition at line 104 of file ClusterAna_module.cc.

Referenced by analyze(), ClusterAna(), and endJob().

std::vector<float> cluster::ClusterAna::fProtKERange

private

Definition at line 99 of file ClusterAna_module.cc.

Referenced by analyze().

bool cluster::ClusterAna::fSkipCosmics

private

Definition at line 103 of file ClusterAna_module.cc.

Referenced by analyze(), and beginJob().

short cluster::ClusterAna::fTrackWeightOption

private

Definition at line 100 of file ClusterAna_module.cc.

Referenced by analyze().

std::string cluster::ClusterAna::fVertexModuleLabel

private

Definition at line 94 of file ClusterAna_module.cc.

Referenced by analyze().

short cluster::ClusterAna::moduleID

private

Definition at line 105 of file ClusterAna_module.cc.

Referenced by analyze(), and ClusterAna().

std::ofstream cluster::ClusterAna::outFile

private

Definition at line 107 of file ClusterAna_module.cc.

Referenced by analyze(), ClusterAna(), and endJob().

---

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

• larreco/v09_25_00/source/larreco/ClusterFinder/ClusterAna_module.cc