Inheritance diagram for cluster::DBclusterAna:

Public Types
using	ModuleType = EDAnalyzer

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

virtual	~DBclusterAna ()

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)
	read access to event More...

void	beginJob ()

Private Attributes
TH1F *	fNoParticles_pdg

TH1F *	fNoParticles_trackid

TH1F *	fNoParticles_trackid_mother

TH1F *	fNoParticles_trackid_per_event

TH1F *	fNoParticles_pdg_per_event

TH1F *	fCl_for_Muon

TH1F *	fNoClustersInEvent

TH1F *	fPercentNoise

TH1F *	fno_of_clusters_per_track

TH1F *	fPercent_lost_muon_hits

TH1F *	fPercent_lost_electron_hits

TH1F *	fPercent_lost_positron_hits

TH1F *	fPercent_lost_111_hits

TH1F *	fPercent_lost_211_hits

TH1F *	fPercent_lost_m211_hits

TH1F *	fPercent_lost_2212_hits

TH1F *	fPercent_lost_2112_hits

TH1F *	fPercent_lost_muon_energy

TH1F *	fPercent_lost_electron_energy

TH1F *	fPercent_lost_positron_energy

TH1F *	fPercent_lost_111_energy

TH1F *	fPercent_lost_211_energy

TH1F *	fPercent_lost_m211_energy

TH1F *	fPercent_lost_2212_energy

TH1F *	fPercent_lost_2112_energy

TH1F *	fEnergy

TH2F *	fbrian_in

TH2F *	fbrian_coll

std::string	fDigitModuleLabel

std::string	fHitsModuleLabel

std::string	fLArG4ModuleLabel

std::string	fClusterFinderModuleLabel

std::string	fCalDataModuleLabel

std::string	fGenieGenModuleLabel

Detailed Description

Definition at line 45 of file DBclusterAna_module.cc.

Member Typedef Documentation

using art::EDAnalyzer::ModuleType = EDAnalyzer

inherited

Definition at line 22 of file EDAnalyzer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 100 of file DBclusterAna_module.cc.

     : EDAnalyzer(pset)
     , fDigitModuleLabel(pset.get<std::string>("DigitModuleLabel"))
     , fHitsModuleLabel(pset.get<std::string>("HitsModuleLabel"))
     , fLArG4ModuleLabel(pset.get<std::string>("LArGeantModuleLabel"))
     , fClusterFinderModuleLabel(pset.get<std::string>("ClusterFinderModuleLabel"))
     , fCalDataModuleLabel(pset.get<std::string>("CalDataModuleLabel"))
     , fGenieGenModuleLabel(pset.get<std::string>("GenieGenModuleLabel"))
   {}

cluster::DBclusterAna::~DBclusterAna ( )

virtual

Definition at line 111 of file DBclusterAna_module.cc.

111 {}

Member Function Documentation

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

private

read access to event

Definition at line 266 of file DBclusterAna_module.cc.

References util::abs(), sim::ParticleList::at(), sim::ParticleList::begin(), DEFINE_ART_MODULE, sim::ParticleList::end(), art::EventID::event(), fCalDataModuleLabel, fCl_for_Muon, fClusterFinderModuleLabel, fDigitModuleLabel, fGenieGenModuleLabel, fHitsModuleLabel, art::fill_ptr_vector(), fno_of_clusters_per_track, fNoClustersInEvent, fNoParticles_pdg, fNoParticles_pdg_per_event, fNoParticles_trackid, fNoParticles_trackid_mother, fNoParticles_trackid_per_event, fPercent_lost_111_energy, fPercent_lost_111_hits, fPercent_lost_2112_energy, fPercent_lost_2112_hits, fPercent_lost_211_energy, fPercent_lost_211_hits, fPercent_lost_2212_energy, fPercent_lost_2212_hits, fPercent_lost_electron_energy, fPercent_lost_electron_hits, fPercent_lost_m211_energy, fPercent_lost_m211_hits, fPercent_lost_muon_energy, fPercent_lost_muon_hits, fPercent_lost_positron_energy, fPercent_lost_positron_hits, fPercentNoise, art::ProductRetriever::getByLabel(), simb::MCTruth::GetParticle(), hits(), cheat::BackTrackerService::HitToEveTrackIDEs(), cheat::BackTrackerService::HitToTrackIDEs(), art::Event::id(), art::Event::isRealData(), geo::GeometryCore::Iterate(), simb::MCTruth::NParticles(), part, cheat::ParticleInventoryService::ParticleList(), simb::MCParticle::PdgCode(), art::PtrVector< T >::push_back(), art::Event::run(), cheat::ParticleInventoryService::SetEveIdCalculator(), and art::PtrVector< T >::size().

   {
     std::cout << "run    : " << evt.run() << std::endl;
     std::cout << "event  : " << evt.id().event() << std::endl;
     //----------------------------------------------------------------
 
     /* This is basically a module for studying MC efficiency/purity. Kick out now if not MC. EC, 8-Oct-2010 */
     if (evt.isRealData()) {
       std::cout << "**** DBclusterAna: Bailing. Don't call this module if you're not MC. "
                 << std::endl;
       exit(1);
     }
 
     art::ServiceHandle<geo::Geometry const> geom;
     art::ServiceHandle<cheat::BackTrackerService const> bt_serv;
     art::ServiceHandle<cheat::ParticleInventoryService> pi_serv;
 
     art::Handle<std::vector<raw::RawDigit>> rdListHandle;
     evt.getByLabel(fDigitModuleLabel, rdListHandle);
     art::Handle<std::vector<recob::Hit>> hitListHandle;
     evt.getByLabel(fHitsModuleLabel, hitListHandle);
     art::Handle<std::vector<simb::MCTruth>> mctruthListHandle;
     evt.getByLabel(fGenieGenModuleLabel, mctruthListHandle);
     art::Handle<std::vector<recob::Cluster>> clusterListHandle;
     evt.getByLabel(fClusterFinderModuleLabel, clusterListHandle);
     art::Handle<std::vector<recob::Wire>> wireListHandle;
     evt.getByLabel(fCalDataModuleLabel, wireListHandle);
 
     art::FindManyP<recob::Hit> fmh(clusterListHandle, evt, fClusterFinderModuleLabel);
 
     //----------------------------------------------------------------
 
     //----------------------------------------------------------------
 
     art::PtrVector<raw::RawDigit> rawdigits;
 
     for (size_t ii = 0; ii < rdListHandle->size(); ++ii) {
       art::Ptr<raw::RawDigit> rawdigit(rdListHandle, ii);
       rawdigits.push_back(rawdigit);
     }
 
     //get the simb::MCParticle collection from the art::Event and then use the
     //Simulation/SimListUtils object to create a sim::ParticleList from the art::Event.
     pi_serv->SetEveIdCalculator(new sim::EmEveIdCalculator);
 
     sim::ParticleList const& _particleList = pi_serv->ParticleList();
 
     std::vector<int> mc_trackids;
 
     for (auto i = _particleList.begin(); i != _particleList.end(); ++i) {
       int trackID = (*i).first;
       mc_trackids.push_back(trackID);
     }
 
     //----------------------------------------------------------------
     art::PtrVector<simb::MCTruth> mclist;
     for (size_t ii = 0; ii < mctruthListHandle->size(); ++ii) {
       art::Ptr<simb::MCTruth> mctparticle(mctruthListHandle, ii);
       mclist.push_back(mctparticle);
     }
 
     //--------------------------------------------------
     std::vector<art::Ptr<recob::Hit>> hits_vec;
 
     //--------------------------------------------------
     std::vector<art::Ptr<recob::Hit>> hits;
     art::fill_ptr_vector(hits, hitListHandle);
     //---------------------------------------------------
 
     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);
     }
 
     std::cout << "in Efficiency, clusters.size()= " << clusters.size() << std::endl;
 
     //---------------------------------------------------------------
     art::PtrVector<recob::Wire> wirelist;
 
     for (size_t ii = 0; ii < wireListHandle->size(); ++ii) {
       art::Ptr<recob::Wire> wireHolder(wireListHandle, ii);
 
       wirelist.push_back(wireHolder);
     }
 
     //...........................................................................
     // How many different particles do we have in a cluster???
     //  -- I will answer that by 3 methods
     //  1) count different TrackIDs in each cluster
     //  2)  count different TrackIDs in each cluster with the use of "Mother" to get rid of the ones that just randomly got their TrackIDs changed by Geant4
     //  3) count different pdg codes in each cluster <--probably not very good b/c if you have 2 "different" electrons in a cluster it's going to count them as one.
     //..........................................................................
     // How many clusters it takes to contain a particle???
     // take each TrackID and count in how many clusters it appears
     //.........................................................................
 
     double no_of_particles_in_cluster = 0;
     double sum_vec_trackid = 0;
     double no_of_clusters = 0;
     double total_no_hits_in_clusters = 0;
     std::vector<int> vec_pdg;
     std::vector<int> vec_trackid, vec_trackid_mother, vec_trackid_mother_en;
     std::vector<int> all_trackids;
     std::vector<int> ids;
     std::vector<int>::iterator it, it2, it3, it4, it5, it6, it7, it8;
     int no_cl_for_muon = 0;
     int no_cl_for_electron = 0;
     int no_cl_for_positron = 0;
     int no_cl_for_pion_111 = 0;
     int no_cl_for_pion_211 = 0;
     int no_cl_for_pion_m211 = 0;
     int no_cl_for_proton = 0;
     double noCluster = 0;
     double _hit_13 = 0, _hit_11 = 0, _hit_m_11 = 0, _hit_111 = 0, _hit_211 = 0, _hit_m211 = 0,
            _hit_2212 = 0, _hit_2112 = 0;
     double _en_13 = 0, _en_11 = 0, _en_m11 = 0, _en_111 = 0, _en_211 = 0, _en_m211 = 0,
            _en_2212 = 0, _en_2112 = 0;
     std::vector<double> diff_vec;
 
     double hit_energy = 0;
 
     auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(evt);
     if (clusters.size() != 0 && hits.size() != 0) {
       for (auto const& plane : geom->Iterate<geo::PlaneGeo>(geo::TPCID{0, 0})) {
         geo::View_t view = plane.View();
         for (size_t j = 0; j < clusters.size(); ++j) {
 
           if (clusters[j]->View() == view) {
 
             std::vector<art::Ptr<recob::Hit>> _hits = fmh.at(j);
             art::Ptr<recob::Hit> _hits_ptr;
 
             for (size_t p = 0; p < _hits.size(); ++p) {
               _hits_ptr = _hits[p];
               hits_vec.push_back(_hits_ptr);
             }
 
             std::vector<art::Ptr<recob::Hit>>::iterator itr = hits_vec.begin();
 
             while (itr != hits_vec.end()) {
               diff_vec.clear();
 
               hit_energy = _hits[itr - hits_vec.begin()]->Integral();
 
               std::vector<sim::TrackIDE> trackides = bt_serv->HitToTrackIDEs(clockData, *itr);
               std::vector<sim::TrackIDE> eveides = bt_serv->HitToEveTrackIDEs(clockData, *itr);
 
               std::vector<sim::TrackIDE>::iterator idesitr = trackides.begin();
 
               while (idesitr != trackides.end()) {
 
                 vec_trackid.push_back((*idesitr).trackID);
 
                 const simb::MCParticle* particle = _particleList.at((*idesitr).trackID);
                 int pdg = particle->PdgCode();
                 if (pdg == 13 || pdg == -13) {
                   _hit_13++;
                   _en_13 += hit_energy * ((*idesitr).energyFrac);
                 }
                 idesitr++;
               }
 
               itr++;
 
             } // loop thru hits
 
             it = find(vec_pdg.begin(), vec_pdg.end(), 13);
             if (it != vec_pdg.end()) { no_cl_for_muon++; }
 
             it2 = find(vec_pdg.begin(), vec_pdg.end(), 11);
             if (it2 != vec_pdg.end()) { no_cl_for_electron++; }
 
             it3 = find(vec_pdg.begin(), vec_pdg.end(), -11);
             if (it3 != vec_pdg.end()) { no_cl_for_positron++; }
 
             it4 = find(vec_pdg.begin(), vec_pdg.end(), 111);
             if (it4 != vec_pdg.end()) { no_cl_for_pion_111++; }
             it6 = find(vec_pdg.begin(), vec_pdg.end(), 211);
             if (it6 != vec_pdg.end()) { no_cl_for_pion_211++; }
             it7 = find(vec_pdg.begin(), vec_pdg.end(), -211);
             if (it7 != vec_pdg.end()) { no_cl_for_pion_m211++; }
             it8 = find(vec_pdg.begin(), vec_pdg.end(), 2212);
             if (it8 != vec_pdg.end()) { no_cl_for_proton++; }
 
             sort(vec_pdg.begin(), vec_pdg.end());
             vec_pdg.erase(unique(vec_pdg.begin(), vec_pdg.end()), vec_pdg.end());
 
             //same for vec_trackid:
 
             sort(vec_trackid.begin(), vec_trackid.end());
             vec_trackid.erase(unique(vec_trackid.begin(), vec_trackid.end()), vec_trackid.end());
             for (unsigned int ii = 0; ii < vec_trackid.size(); ++ii) {
               all_trackids.push_back(vec_trackid[ii]);
             }
 
             //...............................................................
             //Also Make vec_trackid_mother unique:
 
             sort(vec_trackid_mother.begin(), vec_trackid_mother.end());
             vec_trackid_mother.erase(unique(vec_trackid_mother.begin(), vec_trackid_mother.end()),
                                      vec_trackid_mother.end());
 
             // Also Make vec_trackid_mother_en unique:
 
             sort(vec_trackid_mother_en.begin(), vec_trackid_mother_en.end());
             vec_trackid_mother_en.erase(
               unique(vec_trackid_mother_en.begin(), vec_trackid_mother_en.end()),
               vec_trackid_mother_en.end());
 
             //........................................................................
 
             // Q: How many clusters it takes to contain a certain particle?
             for (unsigned int ii = 0; ii < mc_trackids.size(); ++ii) {
               it5 = find(vec_trackid.begin(), vec_trackid.end(), mc_trackids[ii]);
               if (it5 != vec_trackid.end()) {
                 ids.push_back(mc_trackids[ii]); // then make it unique
                 noCluster++;
               }
             }
 
             fNoParticles_pdg->Fill(vec_pdg.size());
 
             fNoParticles_trackid->Fill(vec_trackid.size());
 
             fNoParticles_trackid_mother->Fill(vec_trackid_mother.size());
 
             no_of_clusters++;
 
             no_of_particles_in_cluster += vec_pdg.size();
             sum_vec_trackid += vec_trackid_mother.size();
             total_no_hits_in_clusters += _hits.size();
 
             vec_pdg.clear();
 
             vec_trackid.clear();
 
             vec_trackid_mother.clear();
 
             vec_trackid_mother_en.clear();
 
           } // end if cluster is in correct view
 
           hits_vec.clear();
         } // for each cluster
 
         sort(all_trackids.begin(), all_trackids.end());
         all_trackids.erase(unique(all_trackids.begin(), all_trackids.end()), all_trackids.end());
 
         // now I have a vector(all_trackids) that only contains unique trackids
         // go to it and search for each trackid in every cluster
 
         sort(ids.begin(), ids.end());
         ids.erase(unique(ids.begin(), ids.end()), ids.end());
         double no_of_clusters_per_track = noCluster / ids.size();
 
         // Filling Histograms:
 
         if (no_cl_for_muon != 0) { fCl_for_Muon->Fill(no_cl_for_muon); }
 
         fno_of_clusters_per_track->Fill(no_of_clusters_per_track);
 
         no_cl_for_muon = 0;
         no_cl_for_electron = 0;
         no_cl_for_positron = 0;
         no_cl_for_pion_111 = 0;
         no_cl_for_pion_211 = 0;
         no_cl_for_pion_m211 = 0;
         no_cl_for_proton = 0;
         noCluster = 0;
         ids.clear();
 
       } // for each plane
       double result = no_of_particles_in_cluster / no_of_clusters;
 
       double no_noise_hits = hits.size() - total_no_hits_in_clusters;
       double percent_noise = double(no_noise_hits / hits.size()) * 100;
 
       double no_trackid_per_cl_per_event = sum_vec_trackid / no_of_clusters;
       fNoParticles_trackid_per_event->Fill(no_trackid_per_cl_per_event);
       fNoParticles_pdg_per_event->Fill(result);
       fNoClustersInEvent->Fill(clusters.size());
 
       fPercentNoise->Fill(percent_noise);
 
       //-----------------------------------------------------------------------
       for (unsigned int i = 0; i < mclist.size(); ++i) {
         art::Ptr<simb::MCTruth> mc(mclist[i]);
 
         for (int ii = 0; ii < mc->NParticles(); ++ii) {
           simb::MCParticle part(mc->GetParticle(ii));
           std::cout << "FROM MC TRUTH,the particle's pdg code is: " << part.PdgCode() << std::endl;
           std::cout << "with energy= " << part.E();
           if (abs(part.PdgCode()) == 13) { std::cout << " I have a muon!!!" << std::endl; }
           if (abs(part.PdgCode()) == 111) { std::cout << " I have a pi zero!!!" << std::endl; }
         }
       }
     }
 
     // Now I can load just hits (not clusters) and see how many of what kind are
     // lost due to dbscan marking them as noise:
 
     // Load hits and copy most of the code from above:
 
     double hit_13 = 0, hit_11 = 0, hit_m_11 = 0, hit_111 = 0, hit_211 = 0, hit_m211 = 0,
            hit_2212 = 0, hit_2112 = 0;
     double en_13 = 0, en_11 = 0, en_m11 = 0, en_111 = 0, en_211 = 0, en_m211 = 0, en_2212 = 0,
            en_2112 = 0;
 
     std::vector<art::Ptr<recob::Hit>>::iterator itr = hits.begin();
     while (itr != hits.end()) {
 
       std::vector<sim::TrackIDE> trackides = bt_serv->HitToTrackIDEs(clockData, *itr);
       std::vector<sim::TrackIDE> eveides = bt_serv->HitToEveTrackIDEs(clockData, *itr);
 
       std::vector<sim::TrackIDE>::iterator idesitr = trackides.begin();
 
       hit_energy = hits[itr - hits.begin()]->Integral();
 
       while (idesitr != trackides.end()) {
 
         const simb::MCParticle* particle = _particleList.at((*idesitr).trackID);
 
         int pdg = particle->PdgCode();
         diff_vec.clear();
 
         if (pdg == 13 || pdg == -13) {
           hit_13++;
           en_13 += hit_energy * ((*idesitr).energyFrac);
         }
         idesitr++;
 
       } // trackIDs
 
       itr++;
     } // hits
 
     if (hit_13 != 0) { fPercent_lost_muon_hits->Fill(100 - ((_hit_13 / hit_13) * 100)); }
     if (hit_11 != 0) { fPercent_lost_electron_hits->Fill(100 - ((_hit_11 / hit_11) * 100)); }
     if (hit_m_11 != 0) { fPercent_lost_positron_hits->Fill(100 - ((_hit_m_11 / hit_m_11) * 100)); }
     if (hit_111 != 0) { fPercent_lost_111_hits->Fill(100 - ((_hit_111 / hit_111) * 100)); }
 
     if (hit_211 != 0) { fPercent_lost_211_hits->Fill(100 - ((_hit_211 / hit_211) * 100)); }
 
     if (hit_m211 != 0) { fPercent_lost_m211_hits->Fill(100 - ((_hit_m211 / hit_m211) * 100)); }
     if (hit_2212 != 0) { fPercent_lost_2212_hits->Fill(100 - ((_hit_2212 / hit_2212) * 100)); }
     if (hit_2112 != 0) { fPercent_lost_2112_hits->Fill(100 - ((_hit_2112 / hit_2112) * 100)); }
 
     std::cout << "****** mu E from clusters = " << _en_13 << std::endl;
     std::cout << "****** mu E from hits = " << en_13 << std::endl;
 
     if (en_13 != 0) { fPercent_lost_muon_energy->Fill(100 - ((_en_13 / en_13) * 100)); }
     if (en_11 != 0) { fPercent_lost_electron_energy->Fill(100 - ((_en_11 / en_11) * 100)); }
     if (en_m11 != 0) { fPercent_lost_positron_energy->Fill(100 - ((_en_m11 / en_m11) * 100)); }
     if (en_111 != 0) { fPercent_lost_111_energy->Fill(100 - ((_en_111 / en_111) * 100)); }
     if (en_211 != 0) { fPercent_lost_211_energy->Fill(100 - ((_en_211 / en_211) * 100)); }
     if (en_m211 != 0) { fPercent_lost_m211_energy->Fill(100 - ((_en_m211 / en_m211) * 100)); }
     if (en_2212 != 0) { fPercent_lost_2212_energy->Fill(100 - ((_en_2212 / en_2212) * 100)); }
     if (en_2112 != 0) { fPercent_lost_2112_energy->Fill(100 - ((_en_2112 / en_2112) * 100)); }
   }

void cluster::DBclusterAna::beginJob ( )

privatevirtual

Reimplemented from art::EDAnalyzer.

Definition at line 113 of file DBclusterAna_module.cc.

References fbrian_coll, fbrian_in, fCl_for_Muon, fEnergy, fno_of_clusters_per_track, fNoClustersInEvent, fNoParticles_pdg, fNoParticles_pdg_per_event, fNoParticles_trackid, fNoParticles_trackid_mother, fNoParticles_trackid_per_event, fPercent_lost_111_energy, fPercent_lost_111_hits, fPercent_lost_2112_energy, fPercent_lost_2112_hits, fPercent_lost_211_energy, fPercent_lost_211_hits, fPercent_lost_2212_energy, fPercent_lost_2212_hits, fPercent_lost_electron_energy, fPercent_lost_electron_hits, fPercent_lost_m211_energy, fPercent_lost_m211_hits, fPercent_lost_muon_energy, fPercent_lost_muon_hits, fPercent_lost_positron_energy, fPercent_lost_positron_hits, and fPercentNoise.

   {
     // get access to the TFile service
     art::ServiceHandle<art::TFileService const> tfs;
 
     fNoParticles_pdg_per_event = tfs->make<TH1F>(
       "fNoParticles_pdg_per_event", "Average # of Particles per cluster for each event", 500, 0, 5);
     fNoParticles_pdg = tfs->make<TH1F>(
       "fNoParticles_pdg", "Number of Particles in a Cluster for each cluster", 500, 0, 5);
     fNoParticles_trackid = tfs->make<TH1F>(
       "fNoParticles_trackid", "Number of different TrackIDs in a Cluster", 300, 0, 30);
 
     fNoParticles_trackid_mother =
       tfs->make<TH1F>("fNoParticles_trackid_mother",
                       "Number of different TrackIDs in a Cluster(using mother)for each cluster",
                       300,
                       0,
                       30);
 
     fNoParticles_trackid_per_event =
       tfs->make<TH1F>("fNoParticles_trackid_per_event",
                       "Avg Number of different TrackIDs per Cluster per event",
                       300,
                       0,
                       30);
     fCl_for_Muon =
       tfs->make<TH1F>("fCl_for_Muon", "Number of Clusters for Muon per plane (pdg)", 1500, 0, 15);
 
     fNoClustersInEvent =
       tfs->make<TH1F>("fNoClustersInEvent", "Number of Clusters in an Event", 50, 0, 50);
 
     fPercentNoise =
       tfs->make<TH1F>("fPercentNoise", "% of hits that were marked as Noise by DBSCAN", 250, 0, 25);
 
     fno_of_clusters_per_track = tfs->make<TH1F>(
       "fno_of_clusters_per_track", "Number of Clusters per TrackID per plane", 1500, 0, 15);
 
     fPercent_lost_muon_hits =
       tfs->make<TH1F>("fPercent_lost_muon_hits",
                       "Number of muon hits excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_electron_hits =
       tfs->make<TH1F>("fPercent_lost_electron_hits",
                       "Number of electron hits excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_positron_hits =
       tfs->make<TH1F>("fPercent_lost_positron_hits",
                       "Number of positron hits excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_111_hits =
       tfs->make<TH1F>("fPercent_lost_111_hits",
                       "Number of pion(111) hits excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_211_hits =
       tfs->make<TH1F>("fPercent_lost_211_hits",
                       "Number of pion(211) hits excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_m211_hits =
       tfs->make<TH1F>("fPercent_lost_m211_hits",
                       "Number of pion(-211) hits excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_2212_hits =
       tfs->make<TH1F>("fPercent_lost_2212_hits",
                       "Number of proton hits excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_2112_hits =
       tfs->make<TH1F>("fPercent_lost_2112_hits",
                       "Number of neutron hits excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
 
     fPercent_lost_muon_energy = tfs->make<TH1F>("fPercent_lost_muon_energy",
                                                 " muon energy excluded by dbscan in % (per Event)",
                                                 10000,
                                                 0,
                                                 100);
     fPercent_lost_electron_energy =
       tfs->make<TH1F>("fPercent_lost_electron_energy",
                       "electron energy excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_positron_energy =
       tfs->make<TH1F>("fPercent_lost_positron_energy",
                       " positron energy excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_111_energy =
       tfs->make<TH1F>("fPercent_lost_111_energy",
                       "pion(111) energy excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_211_energy =
       tfs->make<TH1F>("fPercent_lost_211_energy",
                       "pion(211) energy excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_m211_energy =
       tfs->make<TH1F>("fPercent_lost_m211_energy",
                       " pion(-211) energy excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
     fPercent_lost_2212_energy = tfs->make<TH1F>("fPercent_lost_2212_energy",
                                                 "proton energy excluded by dbscan in % (per Event)",
                                                 10000,
                                                 0,
                                                 100);
     fPercent_lost_2112_energy =
       tfs->make<TH1F>("fPercent_lost_2112_energy",
                       "neutron energy excluded by dbscan in % (per Event)",
                       10000,
                       0,
                       100);
 
     fEnergy = tfs->make<TH1F>("fEnergy", "energy for each voxel", 100000, 0, 0.0005);
 
     fbrian_in = tfs->make<TH2F>("fbrian_in",
                                 ";# Electrons deposited; # Electrons detected by hitfinder",
                                 1000,
                                 0,
                                 10000000,
                                 1000,
                                 0,
                                 10000000);
     fbrian_coll = tfs->make<TH2F>("fbrian_coll",
                                   ";# Electrons deposited; # Electrons detected by hitfinder",
                                   1000,
                                   0,
                                   10000000,
                                   1000,
                                   0,
                                   10000000);
   }

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

TH2F* cluster::DBclusterAna::fbrian_coll

private

Definition at line 84 of file DBclusterAna_module.cc.

Referenced by beginJob().

TH2F* cluster::DBclusterAna::fbrian_in

private

Definition at line 83 of file DBclusterAna_module.cc.

Referenced by beginJob().

std::string cluster::DBclusterAna::fCalDataModuleLabel

private

Definition at line 90 of file DBclusterAna_module.cc.

Referenced by analyze().

TH1F* cluster::DBclusterAna::fCl_for_Muon

private

Definition at line 61 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

std::string cluster::DBclusterAna::fClusterFinderModuleLabel

private

Definition at line 89 of file DBclusterAna_module.cc.

Referenced by analyze().

std::string cluster::DBclusterAna::fDigitModuleLabel

private

Definition at line 86 of file DBclusterAna_module.cc.

Referenced by analyze().

TH1F* cluster::DBclusterAna::fEnergy

private

Definition at line 82 of file DBclusterAna_module.cc.

Referenced by beginJob().

std::string cluster::DBclusterAna::fGenieGenModuleLabel

private

Definition at line 91 of file DBclusterAna_module.cc.

Referenced by analyze().

std::string cluster::DBclusterAna::fHitsModuleLabel

private

Definition at line 87 of file DBclusterAna_module.cc.

Referenced by analyze().

std::string cluster::DBclusterAna::fLArG4ModuleLabel

private

Definition at line 88 of file DBclusterAna_module.cc.

TH1F* cluster::DBclusterAna::fno_of_clusters_per_track

private

Definition at line 64 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fNoClustersInEvent

private

Definition at line 62 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fNoParticles_pdg

private

Definition at line 56 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fNoParticles_pdg_per_event

private

Definition at line 60 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fNoParticles_trackid

private

Definition at line 57 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fNoParticles_trackid_mother

private

Definition at line 58 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fNoParticles_trackid_per_event

private

Definition at line 59 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_111_energy

private

Definition at line 77 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_111_hits

private

Definition at line 68 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_2112_energy

private

Definition at line 81 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_2112_hits

private

Definition at line 72 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_211_energy

private

Definition at line 78 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_211_hits

private

Definition at line 69 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_2212_energy

private

Definition at line 80 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_2212_hits

private

Definition at line 71 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_electron_energy

private

Definition at line 75 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_electron_hits

private

Definition at line 66 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_m211_energy

private

Definition at line 79 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_m211_hits

private

Definition at line 70 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_muon_energy

private

Definition at line 74 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_muon_hits

private

Definition at line 65 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_positron_energy

private

Definition at line 76 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercent_lost_positron_hits

private

Definition at line 67 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

TH1F* cluster::DBclusterAna::fPercentNoise

private

Definition at line 63 of file DBclusterAna_module.cc.

Referenced by analyze(), and beginJob().

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

larreco/v09_25_00/source/larreco/ClusterFinder/DBclusterAna_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