Inheritance diagram for HitParentTest:

Public Member Functions
	HitParentTest (const edm::ParameterSet &ps)

	~HitParentTest () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const *> const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c) override

void	beginJob () override

void	beginRun (edm::Run const &, edm::EventSetup const &) override

void	endJob () override

void	endRun (edm::Run const &, edm::EventSetup const &) override

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Private Member Functions
void	analyzeAPDHits (const std::vector< PCaloHit > &, const edm::Handle< edm::SimTrackContainer > &, const edm::Handle< edm::SimVertexContainer > &, int depth)

void	analyzeHits (const std::vector< PCaloHit > &, const edm::Handle< edm::SimTrackContainer > &, const edm::Handle< edm::SimVertexContainer > &, int type)

math::XYZTLorentzVectorD	getOldestParentVertex (const edm::Handle< edm::SimTrackContainer > &, const edm::Handle< edm::SimVertexContainer > &, edm::SimTrackContainer::const_iterator track)

edm::SimTrackContainer::const_iterator	parentSimTrack (const edm::Handle< edm::SimTrackContainer > &, const edm::Handle< edm::SimVertexContainer > &, edm::SimTrackContainer::const_iterator thisTrkItr)

bool	simTrackPresent (const edm::Handle< edm::SimTrackContainer > &, int id)

bool	validSimTrack (const edm::Handle< edm::SimTrackContainer > &, const edm::Handle< edm::SimVertexContainer > &, unsigned int simTkId, edm::SimTrackContainer::const_iterator thisTrkItr)

Private Attributes
std::string	detector [7]

const std::string	g4Label_

bool	histos

const std::string	hitLabEB_

const std::string	hitLabEE_

const std::string	hitLabES_

const std::string	hitLabHC_

TH1F *	hitRho [7]

TH1F *	hitType [7]

TH1F *	hitZ [7]

unsigned int	noGenParticle [7]

unsigned int	noParent [7]

unsigned int	noSimTrack [7]

unsigned int	num_apd_hits_no_gen_particle [2]

unsigned int	num_apd_hits_no_parent [2]

unsigned int	num_apd_hits_no_simtrack [2]

std::map< int, unsigned >	particle_type_count

const std::string	sourceLabel

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_eb_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_ee_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_es_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hc_

const edm::EDGetTokenT< edm::SimTrackContainer >	tok_tk_

const edm::EDGetTokenT< edm::SimVertexContainer >	tok_vtx_

unsigned int	total_num_apd_hits_seen [2]

unsigned int	totalHits [7]

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Detailed Description

Definition at line 31 of file HitParentTest.cc.

Constructor & Destructor Documentation

◆ HitParentTest()

HitParentTest::HitParentTest ( const edm::ParameterSet & ps )

Definition at line 98 of file HitParentTest.cc.

References detector, g4Label_, hitLabEB_, hitLabEE_, hitLabES_, hitLabHC_, mps_fire::i, TFileService::kSharedResource, noGenParticle, noParent, noSimTrack, num_apd_hits_no_gen_particle, num_apd_hits_no_parent, num_apd_hits_no_simtrack, AlCaHLTBitMon_QueryRunRegistry::string, total_num_apd_hits_seen, and totalHits.

     : sourceLabel(ps.getUntrackedParameter<std::string>("SourceLabel", "VtxSmeared")),
       g4Label_(ps.getUntrackedParameter<std::string>("ModuleLabel", "g4SimHits")),
       hitLabEB_(ps.getUntrackedParameter<std::string>("EBCollection", "EcalHitsEB")),
       hitLabEE_(ps.getUntrackedParameter<std::string>("EECollection", "EcalHitsEE")),
       hitLabES_(ps.getUntrackedParameter<std::string>("ESCollection", "EcalHitsES")),
       hitLabHC_(ps.getUntrackedParameter<std::string>("HCCollection", "HcalHits")),
       tok_eb_(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, hitLabEB_))),
       tok_ee_(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, hitLabEE_))),
       tok_es_(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, hitLabES_))),
       tok_hc_(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, hitLabHC_))),
       tok_tk_(consumes<edm::SimTrackContainer>(edm::InputTag(g4Label_))),
       tok_vtx_(consumes<edm::SimVertexContainer>(edm::InputTag(g4Label_))) {
   usesResource(TFileService::kSharedResource);
 
   edm::LogVerbatim("HitParentTest") << "Module Label: " << g4Label_ << "   Hits: " << hitLabEB_ << ", " << hitLabEE_
                                     << ", " << hitLabES_ << ", " << hitLabHC_;
 
   for (unsigned int i = 0; i < 2; ++i) {
     total_num_apd_hits_seen[i] = 0;
     num_apd_hits_no_parent[i] = 0;
     num_apd_hits_no_simtrack[i] = 0;
     num_apd_hits_no_gen_particle[i] = 0;
   }
   std::string dets[7] = {"EB", "EE", "ES", "HB", "HE", "HO", "HF"};
   for (unsigned int i = 0; i < 7; ++i) {
     detector[i] = dets[i];
     totalHits[i] = 0;
     noParent[i] = 0;
     noSimTrack[i] = 0;
     noGenParticle[i] = 0;
   }
 }

◆ ~HitParentTest()

HitParentTest::~HitParentTest ( )

inlineoverride

Definition at line 34 of file HitParentTest.cc.

34 {}

Member Function Documentation

◆ analyze()

void HitParentTest::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 171 of file HitParentTest.cc.

References analyzeAPDHits(), analyzeHits(), hcalRecHitTable_cff::depth, MillePedeFileConverter_cfg::e, edm::HandleBase::isValid(), tok_eb_, tok_ee_, tok_es_, tok_hc_, tok_tk_, and tok_vtx_.

                                                                    {
   edm::LogVerbatim("HitParentTest") << "HitParentTes::Run = " << e.id().run() << " Event = " << e.id().event();
 
   // get PCaloHits for ecal barrel
   const edm::Handle<edm::PCaloHitContainer>& caloHitEB = e.getHandle(tok_eb_);
 
   // get PCaloHits for ecal endcap
   const edm::Handle<edm::PCaloHitContainer>& caloHitEE = e.getHandle(tok_ee_);
 
   // get PCaloHits for preshower
   const edm::Handle<edm::PCaloHitContainer>& caloHitES = e.getHandle(tok_es_);
 
   // get PCaloHits for hcal
   const edm::Handle<edm::PCaloHitContainer>& caloHitHC = e.getHandle(tok_hc_);
 
   // get sim tracks
   const edm::Handle<edm::SimTrackContainer>& simTk = e.getHandle(tok_tk_);
 
   // get sim vertices
   const edm::Handle<edm::SimVertexContainer>& simVtx = e.getHandle(tok_vtx_);
 
   edm::LogVerbatim("HitParentTest") << "HitParentTest: hits valid[EB]: " << caloHitEB.isValid()
                                     << " valid[EE]: " << caloHitEE.isValid() << " valid[ES]: " << caloHitES.isValid()
                                     << " valid[HC]: " << caloHitHC.isValid();
 
   if (caloHitEB.isValid()) {
     for (int depth = 1; depth <= 2; ++depth)
       analyzeAPDHits(*caloHitEB, simTk, simVtx, depth);
     analyzeHits(*caloHitEB, simTk, simVtx, 0);
   }
   if (caloHitEE.isValid())
     analyzeHits(*caloHitEE, simTk, simVtx, 1);
   if (caloHitES.isValid())
     analyzeHits(*caloHitES, simTk, simVtx, 2);
   if (caloHitHC.isValid())
     analyzeHits(*caloHitHC, simTk, simVtx, 3);
 }

◆ analyzeAPDHits()

void HitParentTest::analyzeAPDHits	(	const std::vector< PCaloHit > &	hits,
		const edm::Handle< edm::SimTrackContainer > &	simTk,
		const edm::Handle< edm::SimVertexContainer > &	simVtx,
		int	depth
	)

private

Definition at line 337 of file HitParentTest.cc.

References hcalRecHitTable_cff::depth, newFWLiteAna::found, getOldestParentVertex(), hfClusterShapes_cfi::hits, match(), num_apd_hits_no_gen_particle, num_apd_hits_no_parent, num_apd_hits_no_simtrack, parentSimTrack(), particle_type_count, simTrackPresent(), total_num_apd_hits_seen, and validSimTrack().

Referenced by analyze().

                                               {
   for (std::vector<PCaloHit>::const_iterator hit_it = hits.begin(); hit_it != hits.end(); ++hit_it) {
     if (hit_it->depth() == depth) {
       ++total_num_apd_hits_seen[depth - 1];
 
       // get the geant track id
       int hit_geant_track_id = hit_it->geantTrackId();
 
       if (hit_geant_track_id <= 0) {
         ++num_apd_hits_no_parent[depth - 1];
       } else {
         bool found = false;
         // check whether this id is actually there in the list of simtracks
         for (edm::SimTrackContainer::const_iterator simTrkItr = simTk->begin(); simTrkItr != simTk->end() && !found;
              ++simTrkItr) {
           if (hit_geant_track_id == (int)(simTrkItr->trackId())) {
             found = true;
             bool match = validSimTrack(simTk, simVtx, hit_geant_track_id, simTrkItr);
 
             edm::LogVerbatim("HitParentTest")
                 << "APDHIT Match = " << match << " hit_geant_track_id = " << hit_geant_track_id
                 << " particle id=" << simTrkItr->type();
 
             if (!match) {
               edm::LogVerbatim("HitParentTest") << "NO MATCH FOUND !";
               ++num_apd_hits_no_gen_particle[depth - 1];
             }
 
             int apd_pid = simTrkItr->type();
             std::map<int, unsigned>::iterator count_it = particle_type_count.find(apd_pid);
             if (count_it == particle_type_count.end())
               // first occurence of this particle pid
               particle_type_count[apd_pid] = 1;
             else
               ++count_it->second;
 
             //--------------------
             // check where the oldest parent has its vertex. Should be close to the
             // beam pipe...
             math::XYZTLorentzVectorD oldest_parent_vertex = getOldestParentVertex(simTk, simVtx, simTrkItr);
 
             edm::SimTrackContainer::const_iterator oldest_parent_track = parentSimTrack(simTk, simVtx, simTrkItr);
 
             edm::LogVerbatim("HitParentTest")
                 << "APD hit pid = " << apd_pid << " APD hit track id = " << hit_geant_track_id
                 << " depth = " << hit_it->depth() << " OLDEST PARENT'S VERTEX: " << oldest_parent_vertex
                 << " rho = " << oldest_parent_vertex.Rho() << " OLDEST PARENT'S PID: " << oldest_parent_track->type()
                 << " OLDEST PARENT'S track id: " << oldest_parent_track->trackId() << "\n"
                 << "APD hit vertex index: " << simTrkItr->vertIndex() << " (tot #vertices: " << simVtx->size() << ")"
                 << "\n"
                 << "APD hit vertex parent track: " << (*simVtx)[simTrkItr->vertIndex()].parentIndex()
                 << " present=" << simTrackPresent(simTk, (*simVtx)[simTrkItr->vertIndex()].parentIndex());
 
           }  // a match was found
         }    // geant track id found in simtracks
 
         if (!found)
           ++num_apd_hits_no_simtrack[depth - 1];
       }  // hits with a valid SimTrack Id
     }    // right depth index
   }      // loop over all calohits (of the given depth)
 }

◆ analyzeHits()

void HitParentTest::analyzeHits	(	const std::vector< PCaloHit > &	hits,
		const edm::Handle< edm::SimTrackContainer > &	simTk,
		const edm::Handle< edm::SimVertexContainer > &	simVtx,
		int	type
	)

private

performs some checks on hits

Definition at line 264 of file HitParentTest.cc.

References detector, RemoveAddSevLevel::flag, newFWLiteAna::found, getOldestParentVertex(), HcalEndcap, HcalForward, HcalOuter, histos, hitRho, hfClusterShapes_cfi::hits, hitType, hitZ, l1ctLayer2EG_cff::id, match(), noGenParticle, noParent, noSimTrack, parentSimTrack(), simTrackPresent(), HcalDetId::subdet(), totalHits, and validSimTrack().

Referenced by analyze().

                                           {
   for (std::vector<PCaloHit>::const_iterator hit_it = hits.begin(); hit_it != hits.end(); ++hit_it) {
     int id(type), flag(0);
     if (type == 3) {
       HcalDetId id_ = HcalDetId(hit_it->id());
       int subdet = id_.subdet();
       if (subdet == static_cast<int>(HcalEndcap))
         id = type + 1;
       else if (subdet == static_cast<int>(HcalOuter))
         id = type + 2;
       else if (subdet == static_cast<int>(HcalForward))
         id = type + 3;
     }
     ++totalHits[id];
 
     // get the geant track id
     int hit_geant_track_id = hit_it->geantTrackId();
 
     if (hit_geant_track_id <= 0) {
       ++noParent[id];
       flag = 1;
     } else {
       bool found = false;
       flag = 2;
       // check whether this id is actually there in the list of simtracks
       for (edm::SimTrackContainer::const_iterator simTrkItr = simTk->begin(); simTrkItr != simTk->end() && !found;
            ++simTrkItr) {
         if (hit_geant_track_id == (int)(simTrkItr->trackId())) {
           found = true;
           flag = 3;
           bool match = validSimTrack(simTk, simVtx, hit_geant_track_id, simTrkItr);
 
           edm::LogVerbatim("HitParentTest")
               << "[" << detector[type] << "] Match = " << match << " hit_geant_track_id=" << hit_geant_track_id
               << " particle id=" << simTrkItr->type();
 
           if (!match) {
             edm::LogVerbatim("HitParentTest") << "NO MATCH FOUND !";
             ++noGenParticle[id];
           }
 
           // beam pipe...
           int pid = simTrkItr->type();
           math::XYZTLorentzVectorD oldest_parent_vertex = getOldestParentVertex(simTk, simVtx, simTrkItr);
 
           edm::SimTrackContainer::const_iterator oldest_parent_track = parentSimTrack(simTk, simVtx, simTrkItr);
 
           edm::LogVerbatim("HitParentTest")
               << "[" << detector[type] << "] Hit pid = " << pid << "  hit track id = " << hit_geant_track_id
               << " Oldest Parent's Vertex: " << oldest_parent_vertex << " rho = " << oldest_parent_vertex.Rho()
               << " Oldest Parent's pid: " << oldest_parent_track->type()
               << " Oldest Parent's track id: " << oldest_parent_track->trackId()
               << "\nHit vertex index: " << simTrkItr->vertIndex() << " (tot #vertices: " << simVtx->size() << ")"
               << "\nHit vertex parent track: " << (*simVtx)[simTrkItr->vertIndex()].parentIndex()
               << " present=" << simTrackPresent(simTk, (*simVtx)[simTrkItr->vertIndex()].parentIndex());
           if (histos) {
             hitRho[id]->Fill(oldest_parent_vertex.Rho());
             hitZ[id]->Fill(oldest_parent_vertex.Z());
           }
         }  // a match was found
       }    // geant track id found in simtracks
 
       if (!found)
         ++noSimTrack[id];
     }  // hits with a valid SimTrack Id
     if (histos)
       hitType[id]->Fill((double)(flag));
   }  // loop over all calohits (of the given depth)
 }

◆ beginJob()

void HitParentTest::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 143 of file HitParentTest.cc.

References detector, histos, hitRho, hitType, hitZ, mps_fire::i, Skims_PA_cff::name, compare::tfile, and runGCPTkAlMap::title.

                              {
   edm::Service<TFileService> tfile;
   if (!tfile.isAvailable()) {
     edm::LogVerbatim("HitParentTest") << "TFileService unavailable: no histograms";
     histos = false;
   } else {
     char name[20], title[100];
     histos = true;
     for (unsigned int i = 0; i < 7; ++i) {
       sprintf(name, "HitType%d", i);
       sprintf(title, "Hit types for %s", detector[i].c_str());
       hitType[i] = tfile->make<TH1F>(name, title, 10, 0., 10.);
       hitType[i]->GetXaxis()->SetTitle(title);
       hitType[i]->GetYaxis()->SetTitle("Hits");
       sprintf(name, "RhoVertex%d", i);
       sprintf(title, "#rho of the vertex for %s Hits", detector[i].c_str());
       hitRho[i] = tfile->make<TH1F>(name, title, 10000, 0., 100.);
       hitRho[i]->GetXaxis()->SetTitle(title);
       hitRho[i]->GetYaxis()->SetTitle("Hits");
       sprintf(name, "ZVertex%d", i);
       sprintf(title, "z of the vertex for %s Hits", detector[i].c_str());
       hitZ[i] = tfile->make<TH1F>(name, title, 2000, -100., 100.);
       hitZ[i]->GetXaxis()->SetTitle(title);
       hitZ[i]->GetYaxis()->SetTitle("Hits");
     }
   }
 }

◆ beginRun()

void HitParentTest::beginRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

inlineoverrideprotected

Definition at line 39 of file HitParentTest.cc.

39 {}

◆ endJob()

void HitParentTest::endJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 221 of file HitParentTest.cc.

References hcalRecHitTable_cff::depth, detector, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, noGenParticle, noParent, noSimTrack, num_apd_hits_no_gen_particle, num_apd_hits_no_parent, num_apd_hits_no_simtrack, particle_type_count, jetUpdater_cfi::sort, total_num_apd_hits_seen, and totalHits.

                            {
   edm::LogVerbatim("HitParentTest") << "HitParentTest::Total number of APD hits seen: "
                                     << total_num_apd_hits_seen[0] + total_num_apd_hits_seen[1];
   for (unsigned int depth = 0; depth < 2; ++depth) {
     edm::LogVerbatim("HitParentTest") << "APD Hits in depth = " << depth + 1
                                       << " Total = " << total_num_apd_hits_seen[depth];
     edm::LogVerbatim("HitParentTest") << "summary of errors:\n"
                                       << " number of APD hits with zero geant track id: "
                                       << num_apd_hits_no_parent[depth] << "\n"
                                       << "number of APD hits for which the parent simtrack was not found "
                                       << "in the simtrack collection: " << num_apd_hits_no_simtrack[depth] << "\n"
                                       << "number of APD hits for which no generator particle was "
                                       << "found: " << num_apd_hits_no_gen_particle[depth] << "\n";
   }
 
   for (unsigned int det = 0; det < 7; ++det) {
     edm::LogVerbatim("HitParentTest") << "Total number of hits seen in " << detector[det] << ": " << totalHits[det];
     edm::LogVerbatim("HitParentTest") << "summary of errors:\n"
                                       << "number of hits with zero geant track id: " << noParent[det] << "\n"
                                       << "number of hits for which the parent simtrack was not found in "
                                       << "the simtrack collection: " << noSimTrack[det] << "\n"
                                       << "number of hits for which no generator particle was found: "
                                       << noGenParticle[det] << "\n";
   }
 
   // sort in decreasing order of occurence
   std::vector<int> sorted_pids;
   for (std::map<int, unsigned>::const_iterator it = particle_type_count.begin(); it != particle_type_count.end(); ++it)
     sorted_pids.push_back(it->first);
 
   // now sort the pids
 
   std::sort(sorted_pids.begin(), sorted_pids.end(), HitParentTestComparison(particle_type_count));
 
   edm::LogVerbatim("HitParentTest") << "HitParentTest::Frequency particle types through the APD "
                                     << "(pid/frequency):";
   for (unsigned i = 0; i < sorted_pids.size(); ++i) {
     int pid = sorted_pids[i];
     edm::LogVerbatim("HitParentTest") << "  pid " << std::setw(6) << pid << ": count " << std::setw(6)
                                       << particle_type_count[pid];
   }
 }

◆ endRun()

void HitParentTest::endRun	(	edm::Run const &	,
		edm::EventSetup const &
	)

inlineoverrideprotected

Definition at line 41 of file HitParentTest.cc.

41 {}

◆ fillDescriptions()

void HitParentTest::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 132 of file HitParentTest.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                {
   edm::ParameterSetDescription desc;
   desc.addUntracked<std::string>("SourceLabel", "generatorSmeared");
   desc.addUntracked<std::string>("ModuleLabel", "g4SimHits");
   desc.addUntracked<std::string>("EBCollection", "EcalHitsEB");
   desc.addUntracked<std::string>("EECollection", "EcalHitsEE");
   desc.addUntracked<std::string>("ESCollection", "EcalHitsES");
   desc.addUntracked<std::string>("HCCollection", "HcalHits");
   descriptions.add("HitParentTest", desc);
 }

◆ getOldestParentVertex()

math::XYZTLorentzVectorD HitParentTest::getOldestParentVertex	(	const edm::Handle< edm::SimTrackContainer > &	simTk,
		const edm::Handle< edm::SimVertexContainer > &	simVtx,
		edm::SimTrackContainer::const_iterator	track
	)

private

Definition at line 411 of file HitParentTest.cc.

References parentSimTrack(), and HLT_2024v14_cff::track.

Referenced by analyzeAPDHits(), and analyzeHits().

                                                                                                         {
   static const math::XYZTLorentzVectorD invalid_vertex(
       10000, 10000, 10000, 10000);  // default value if no valid vertex found
 
   edm::SimTrackContainer::const_iterator oldest_parent_track = parentSimTrack(simTk, simVtx, track);
 
   int vertex_index = oldest_parent_track->vertIndex();
 
   // sanity checks
   if (vertex_index < 0 || vertex_index >= (int)(simVtx->size()))
     return invalid_vertex;
 
   return (*simVtx)[vertex_index].position();
 }

◆ parentSimTrack()

edm::SimTrackContainer::const_iterator HitParentTest::parentSimTrack	(	const edm::Handle< edm::SimTrackContainer > &	simTk,
		const edm::Handle< edm::SimVertexContainer > &	simVtx,
		edm::SimTrackContainer::const_iterator	thisTrkItr
	)

private

Definition at line 428 of file HitParentTest.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, createfilelist::int, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, and class-composition::parent.

Referenced by analyzeAPDHits(), analyzeHits(), and getOldestParentVertex().

                                                                                                                     {
   edm::SimTrackContainer::const_iterator itr = simTk->end();
 
   int vertIndex = thisTrkItr->vertIndex();
   int type = thisTrkItr->type();
   int charge = (int)thisTrkItr->charge();
   edm::LogVerbatim("HitParentTest") << "SimTrackParent:: " << thisTrkItr->trackId() << " Vertex " << vertIndex
                                     << " Type " << type << " Charge " << charge;
 
   if (vertIndex == -1)
     return thisTrkItr;
   else if (vertIndex >= (int)simVtx->size())
     return itr;
 
   edm::SimVertexContainer::const_iterator simVtxItr = simVtx->begin();
   for (int iv = 0; iv < vertIndex; iv++)
     simVtxItr++;
   int parent = simVtxItr->parentIndex();
 
   if (parent < 0 && simVtxItr != simVtx->begin()) {
     const math::XYZTLorentzVectorD pos1 = simVtxItr->position();
     for (simVtxItr = simVtx->begin(); simVtxItr != simVtx->end(); ++simVtxItr) {
       if (simVtxItr->parentIndex() > 0) {
         const math::XYZTLorentzVectorD pos2 = pos1 - simVtxItr->position();
         double dist = pos2.P();
         if (dist < 0.001) {
           parent = simVtxItr->parentIndex();
           break;
         }
       }
     }
   }
   for (edm::SimTrackContainer::const_iterator simTrkItr = simTk->begin(); simTrkItr != simTk->end(); simTrkItr++) {
     if ((int)simTrkItr->trackId() == parent && simTrkItr != thisTrkItr)
       return parentSimTrack(simTk, simVtx, simTrkItr);
   }
 
   return thisTrkItr;
 }

◆ simTrackPresent()

bool HitParentTest::simTrackPresent	(	const edm::Handle< edm::SimTrackContainer > &	simTk,
		int	id
	)

private

Definition at line 403 of file HitParentTest.cc.

Referenced by analyzeAPDHits(), and analyzeHits().

                                                                                         {
   for (edm::SimTrackContainer::const_iterator simTrkItr = simTk->begin(); simTrkItr != simTk->end(); ++simTrkItr) {
     if ((int)(simTrkItr->trackId()) == id)
       return true;
   }
   return false;
 }

◆ validSimTrack()

bool HitParentTest::validSimTrack	(	const edm::Handle< edm::SimTrackContainer > &	simTk,
		const edm::Handle< edm::SimVertexContainer > &	simVtx,
		unsigned int	simTkId,
		edm::SimTrackContainer::const_iterator	thisTrkItr
	)

private

Definition at line 470 of file HitParentTest.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, and class-composition::parent.

Referenced by analyzeAPDHits(), and analyzeHits().

                                                                                  {
   edm::LogVerbatim("HitParentTest") << "Inside validSimTrack: trackId " << thisTrkItr->trackId() << " vtxIndex "
                                     << thisTrkItr->vertIndex() << " to be matched to " << simTkId;
 
   //This track originates from simTkId
   if (thisTrkItr->trackId() == simTkId)
     return true;
 
   //Otherwise trace back the history using SimTracks and SimVertices
   int vertIndex = thisTrkItr->vertIndex();
   if (vertIndex == -1 || vertIndex >= (int)simVtx->size())
     return false;
   edm::SimVertexContainer::const_iterator simVtxItr = simVtx->begin();
   for (int iv = 0; iv < vertIndex; iv++)
     simVtxItr++;
   int parent = simVtxItr->parentIndex();
   edm::LogVerbatim("HitParentTest") << "validSimTrack:: parent index " << parent << " ";
   if (parent < 0 && simVtxItr != simVtx->begin()) {
     const math::XYZTLorentzVectorD pos1 = simVtxItr->position();
     for (simVtxItr = simVtx->begin(); simVtxItr != simVtx->end(); ++simVtxItr) {
       if (simVtxItr->parentIndex() > 0) {
         const math::XYZTLorentzVectorD pos2 = pos1 - simVtxItr->position();
         double dist = pos2.P();
         if (dist < 0.001) {
           parent = simVtxItr->parentIndex();
           break;
         }
       }
     }
   }
   edm::LogVerbatim("HitParentTest") << "HitParentTes::final index " << parent;
   for (edm::SimTrackContainer::const_iterator simTrkItr = simTk->begin(); simTrkItr != simTk->end(); simTrkItr++) {
     if ((int)simTrkItr->trackId() == parent && simTrkItr != thisTrkItr)
       return validSimTrack(simTk, simVtx, simTkId, simTrkItr);
   }
 
   return false;
 }

Member Data Documentation

◆ detector

std::string HitParentTest::detector[7]

private

Definition at line 91 of file HitParentTest.cc.

Referenced by analyzeHits(), beginJob(), endJob(), and HitParentTest().

◆ g4Label_

const std::string HitParentTest::g4Label_

private

Definition at line 68 of file HitParentTest.cc.

Referenced by HitParentTest().

◆ histos

bool HitParentTest::histos

private

Definition at line 92 of file HitParentTest.cc.

Referenced by analyzeHits(), and beginJob().

◆ hitLabEB_

const std::string HitParentTest::hitLabEB_

private

Definition at line 68 of file HitParentTest.cc.

Referenced by HitParentTest().

◆ hitLabEE_

const std::string HitParentTest::hitLabEE_

private

Definition at line 68 of file HitParentTest.cc.

Referenced by HitParentTest().

◆ hitLabES_

const std::string HitParentTest::hitLabES_

private

Definition at line 68 of file HitParentTest.cc.

Referenced by HitParentTest().

◆ hitLabHC_

const std::string HitParentTest::hitLabHC_

private

Definition at line 68 of file HitParentTest.cc.

Referenced by HitParentTest().

◆ hitRho

TH1F * HitParentTest::hitRho[7]

private

Definition at line 95 of file HitParentTest.cc.

Referenced by analyzeHits(), and beginJob().

◆ hitType

TH1F* HitParentTest::hitType[7]

private

Definition at line 95 of file HitParentTest.cc.

Referenced by ntupleDataFormat._RecoHitAdaptor::_hits(), analyzeHits(), and beginJob().

◆ hitZ

TH1F * HitParentTest::hitZ[7]

private

Definition at line 95 of file HitParentTest.cc.

Referenced by analyzeHits(), and beginJob().

◆ noGenParticle

unsigned int HitParentTest::noGenParticle[7]

private

Definition at line 94 of file HitParentTest.cc.

Referenced by analyzeHits(), endJob(), and HitParentTest().

◆ noParent

unsigned int HitParentTest::noParent[7]

private

Definition at line 93 of file HitParentTest.cc.

Referenced by analyzeHits(), endJob(), and HitParentTest().

◆ noSimTrack

unsigned int HitParentTest::noSimTrack[7]

private

Definition at line 94 of file HitParentTest.cc.

Referenced by analyzeHits(), endJob(), and HitParentTest().

◆ num_apd_hits_no_gen_particle

unsigned int HitParentTest::num_apd_hits_no_gen_particle[2]

private

number of APD hits for which no generator particle was found

Definition at line 85 of file HitParentTest.cc.

Referenced by analyzeAPDHits(), endJob(), and HitParentTest().

◆ num_apd_hits_no_parent

unsigned int HitParentTest::num_apd_hits_no_parent[2]

private

Definition at line 78 of file HitParentTest.cc.

Referenced by analyzeAPDHits(), endJob(), and HitParentTest().

◆ num_apd_hits_no_simtrack

unsigned int HitParentTest::num_apd_hits_no_simtrack[2]

private

number of apd hits for which the parent sim track id was not found in the simtrack collection.

Definition at line 82 of file HitParentTest.cc.

Referenced by analyzeAPDHits(), endJob(), and HitParentTest().

◆ particle_type_count

std::map<int, unsigned> HitParentTest::particle_type_count

private

'histogram' of types of particles going through the APD. Maps from numeric particle code to the number of occurences.

Definition at line 89 of file HitParentTest.cc.

Referenced by analyzeAPDHits(), and endJob().

◆ sourceLabel

const std::string HitParentTest::sourceLabel

private

Definition at line 68 of file HitParentTest.cc.

◆ tok_eb_

const edm::EDGetTokenT<edm::PCaloHitContainer> HitParentTest::tok_eb_

private

Definition at line 69 of file HitParentTest.cc.

Referenced by analyze().

◆ tok_ee_

const edm::EDGetTokenT<edm::PCaloHitContainer> HitParentTest::tok_ee_

private

Definition at line 70 of file HitParentTest.cc.

Referenced by analyze().

◆ tok_es_

const edm::EDGetTokenT<edm::PCaloHitContainer> HitParentTest::tok_es_

private

Definition at line 71 of file HitParentTest.cc.

Referenced by analyze().

◆ tok_hc_

const edm::EDGetTokenT<edm::PCaloHitContainer> HitParentTest::tok_hc_

private

Definition at line 72 of file HitParentTest.cc.

Referenced by analyze().

◆ tok_tk_

const edm::EDGetTokenT<edm::SimTrackContainer> HitParentTest::tok_tk_

private

Definition at line 73 of file HitParentTest.cc.

Referenced by analyze().

◆ tok_vtx_

const edm::EDGetTokenT<edm::SimVertexContainer> HitParentTest::tok_vtx_

private

Definition at line 74 of file HitParentTest.cc.

Referenced by analyze().

◆ total_num_apd_hits_seen

unsigned int HitParentTest::total_num_apd_hits_seen[2]

private

error and other counters

Definition at line 77 of file HitParentTest.cc.

Referenced by analyzeAPDHits(), endJob(), and HitParentTest().

◆ totalHits

unsigned int HitParentTest::totalHits[7]

private

Definition at line 93 of file HitParentTest.cc.

Referenced by analyzeHits(), endJob(), and HitParentTest().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HitParentTest()

◆ ~HitParentTest()

Member Function Documentation

◆ analyze()

◆ analyzeAPDHits()

◆ analyzeHits()

◆ beginJob()

◆ beginRun()

◆ endJob()

◆ endRun()

◆ fillDescriptions()

◆ getOldestParentVertex()

◆ parentSimTrack()

◆ simTrackPresent()

◆ validSimTrack()

Member Data Documentation

◆ detector

◆ g4Label_

◆ histos

◆ hitLabEB_

◆ hitLabEE_

◆ hitLabES_

◆ hitLabHC_

◆ hitRho

◆ hitType

◆ hitZ

◆ noGenParticle

◆ noParent

◆ noSimTrack

◆ num_apd_hits_no_gen_particle

◆ num_apd_hits_no_parent

◆ num_apd_hits_no_simtrack

◆ particle_type_count

◆ sourceLabel

◆ tok_eb_

◆ tok_ee_

◆ tok_es_

◆ tok_hc_

◆ tok_tk_

◆ tok_vtx_

◆ total_num_apd_hits_seen

◆ totalHits