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

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const 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

bool	wantsStreamRuns () const

	~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

ESProxyIndex const *	esGetTokenIndices (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::vector< ModuleDescription const * > &modules, 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

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices 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)

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 ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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

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

void	analyzeHits (const std::vector< PCaloHit > &, int type)

math::XYZTLorentzVectorD	getOldestParentVertex (edm::SimTrackContainer::const_iterator track)

edm::SimTrackContainer::const_iterator	parentSimTrack (edm::SimTrackContainer::const_iterator thisTrkItr)

bool	simTrackPresent (int id)

bool	validSimTrack (unsigned int simTkId, edm::SimTrackContainer::const_iterator thisTrkItr)

Private Attributes
std::string	detector [7]

std::string	g4Label

bool	histos

std::string	hitLabEB

std::string	hitLabEE

std::string	hitLabES

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

edm::Handle< edm::SimTrackContainer >	SimTk

edm::Handle< edm::SimVertexContainer >	SimVtx

std::string	sourceLabel

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_eb_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_ee_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_es_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hc_

edm::EDGetTokenT< edm::SimTrackContainer >	tok_tk_

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 36 of file HitParentTest.cc.

Constructor & Destructor Documentation

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

Definition at line 93 of file HitParentTest.cc.

References detector, g4Label, edm::ParameterSet::getUntrackedParameter(), hitLabEB, hitLabEE, hitLabES, hitLabHC, mps_fire::i, HLT_2018_cff::InputTag, TFileService::kSharedResource, noGenParticle, noParent, noSimTrack, num_apd_hits_no_gen_particle, num_apd_hits_no_parent, num_apd_hits_no_simtrack, sourceLabel, AlCaHLTBitMon_QueryRunRegistry::string, tok_eb_, tok_ee_, tok_es_, tok_hc_, tok_tk_, tok_vtx_, total_num_apd_hits_seen, and totalHits.

                                                       {
   usesResource(TFileService::kSharedResource);
 
   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");
   edm::LogVerbatim("HitParentTest") << "Module Label: " << g4Label << "   Hits: " << hitLabEB << ", " << hitLabEE
                                     << ", " << hitLabES << ", " << hitLabHC;
 
   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));
 
   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 ( )

inlineoverride

Definition at line 39 of file HitParentTest.cc.

References beginJob(), and fillDescriptions().

39 {}

Member Function Documentation

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

overrideprotected

Definition at line 167 of file HitParentTest.cc.

References analyzeAPDHits(), analyzeHits(), LEDCalibrationChannels::depth, edm::EventID::event(), edm::Event::getByToken(), edm::EventBase::id(), edm::HandleBase::isValid(), edm::EventID::run(), SimTk, SimVtx, tok_eb_, tok_ee_, tok_es_, tok_hc_, tok_tk_, and tok_vtx_.

Referenced by beginRun().

                                                                    {
   edm::LogVerbatim("HitParentTest") << "HitParentTes::Run = " << e.id().run() << " Event = " << e.id().event();
 
   // get PCaloHits for ecal barrel
   edm::Handle<edm::PCaloHitContainer> caloHitEB;
   e.getByToken(tok_eb_, caloHitEB);
 
   // get PCaloHits for ecal endcap
   edm::Handle<edm::PCaloHitContainer> caloHitEE;
   e.getByToken(tok_ee_, caloHitEE);
 
   // get PCaloHits for preshower
   edm::Handle<edm::PCaloHitContainer> caloHitES;
   e.getByToken(tok_es_, caloHitES);
 
   // get PCaloHits for hcal
   edm::Handle<edm::PCaloHitContainer> caloHitHC;
   e.getByToken(tok_hc_, caloHitHC);
 
   // get sim tracks
   e.getByToken(tok_tk_, SimTk);
 
   // get sim vertices
   e.getByToken(tok_vtx_, SimVtx);
 
   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, depth);
     analyzeHits(*caloHitEB, 0);
   }
   if (caloHitEE.isValid())
     analyzeHits(*caloHitEE, 1);
   if (caloHitES.isValid())
     analyzeHits(*caloHitES, 2);
   if (caloHitHC.isValid())
     analyzeHits(*caloHitHC, 3);
 }

void HitParentTest::analyzeAPDHits	(	const std::vector< PCaloHit > &	hits,
		int	depth
	)

private

Definition at line 334 of file HitParentTest.cc.

References LEDCalibrationChannels::depth, newFWLiteAna::found, getOldestParentVertex(), match(), num_apd_hits_no_gen_particle, num_apd_hits_no_parent, num_apd_hits_no_simtrack, parentSimTrack(), particle_type_count, SimTk, simTrackPresent(), SimVtx, total_num_apd_hits_seen, and validSimTrack().

Referenced by analyze(), and endRun().

                                                                              {
   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(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(simTrkItr);
 
             edm::SimTrackContainer::const_iterator oldest_parent_track = parentSimTrack(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((*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)
 }

void HitParentTest::analyzeHits	(	const std::vector< PCaloHit > &	hits,
		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, hitType, hitZ, triggerObjects_cff::id, match(), noGenParticle, noParent, noSimTrack, parentSimTrack(), SimTk, simTrackPresent(), SimVtx, HcalDetId::subdet(), totalHits, and validSimTrack().

Referenced by analyze(), and endRun().

                                                                          {
   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(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(simTrkItr);
 
           edm::SimTrackContainer::const_iterator oldest_parent_track = parentSimTrack(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((*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)
 }

void HitParentTest::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 139 of file HitParentTest.cc.

References detector, histos, hitRho, hitType, hitZ, mps_fire::i, edm::Service< T >::isAvailable(), TFileService::make(), Skims_PA_cff::name, compare::tfile, and overlapproblemtsosanalyzer_cfi::title.

Referenced by ~HitParentTest().

                              {
   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");
     }
   }
 }

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

inlineoverrideprotected

Definition at line 44 of file HitParentTest.cc.

References analyze(), HltBtagPostValidation_cff::c, and MillePedeFileConverter_cfg::e.

44 {}

void HitParentTest::endJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 221 of file HitParentTest.cc.

References LEDCalibrationChannels::depth, detector, dqm-mbProfile::format, mps_fire::i, noGenParticle, noParent, noSimTrack, num_apd_hits_no_gen_particle, num_apd_hits_no_parent, num_apd_hits_no_simtrack, particle_type_count, total_num_apd_hits_seen, and totalHits.

Referenced by endRun().

                            {
   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 " << boost::format("%6d") % pid << ": count "
                                       << boost::format("%6d") % particle_type_count[pid];
   }
 }

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

inlineoverrideprotected

Definition at line 46 of file HitParentTest.cc.

References analyzeAPDHits(), analyzeHits(), LEDCalibrationChannels::depth, endJob(), getOldestParentVertex(), parentSimTrack(), simTrackPresent(), HLT_2018_cff::track, and validSimTrack().

46 {}

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

static

Definition at line 128 of file HitParentTest.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::addUntracked(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~HitParentTest().

                                                                                {
   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);
 }

math::XYZTLorentzVectorD HitParentTest::getOldestParentVertex ( edm::SimTrackContainer::const_iterator track )

private

Definition at line 405 of file HitParentTest.cc.

References parentSimTrack(), and SimVtx.

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

                                                                                                     {
   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(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();
 }

edm::SimTrackContainer::const_iterator HitParentTest::parentSimTrack ( edm::SimTrackContainer::const_iterator thisTrkItr )

private

Definition at line 420 of file HitParentTest.cc.

References begin, ALCARECOTkAlJpsiMuMu_cff::charge, createfilelist::int, class-composition::parent, SimTk, and SimVtx.

Referenced by analyzeAPDHits(), analyzeHits(), endRun(), 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(simTrkItr);
   }
 
   return thisTrkItr;
 }

bool HitParentTest::simTrackPresent ( int id )

private

Definition at line 397 of file HitParentTest.cc.

References SimTk.

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

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

bool HitParentTest::validSimTrack	(	unsigned int	simTkId,
		edm::SimTrackContainer::const_iterator	thisTrkItr
	)

private

Definition at line 460 of file HitParentTest.cc.

References begin, DEFINE_FWK_MODULE, class-composition::parent, SimTk, and SimVtx.

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

                                                                                                      {
   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(simTkId, simTrkItr);
   }
 
   return false;
 }