Inheritance diagram for XtalDedxAnalysis:

Public Member Functions
	XtalDedxAnalysis (const edm::ParameterSet &)

	~XtalDedxAnalysis () override

Public Member Functions inherited from edm::one::EDAnalyzer< 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 &, const edm::EventSetup &) override

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

void	beginJob () override

void	endJob () 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 Attributes
const edm::InputTag	caloHitSource_

const double	energyMax_

TH1F *	meE1T0_ [4]

TH1F *	meE1T1_ [4]

TH1F *	meE9T0_ [4]

TH1F *	meE9T1_ [4]

TH1F *	meNHit_ [4]

TH1I *	mType_

const std::string	simTkLabel_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_calo_

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

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

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 XtalDedxAnalysis.cc.

Constructor & Destructor Documentation

◆ XtalDedxAnalysis()

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

explicit

Definition at line 64 of file XtalDedxAnalysis.cc.

References caloHitSource_, energyMax_, mps_fire::i, TFileService::kSharedResource, meE1T0_, meE1T1_, meE9T0_, meE9T1_, meNHit_, mType_, Skims_PA_cff::name, simTkLabel_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, and runGCPTkAlMap::title.

     : caloHitSource_(ps.getParameter<edm::InputTag>("caloHitSource")),
       simTkLabel_(ps.getParameter<std::string>("moduleLabelTk")),
       energyMax_(ps.getParameter<double>("energyMax")),
       tok_calo_(consumes<edm::PCaloHitContainer>(caloHitSource_)),
       tok_tk_(consumes<edm::SimTrackContainer>(edm::InputTag(simTkLabel_))),
       tok_vtx_(consumes<edm::SimVertexContainer>(edm::InputTag(simTkLabel_))) {
   usesResource(TFileService::kSharedResource);
   edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis::Source " << caloHitSource_ << " Track Label "
                                         << simTkLabel_ << " Energy Max " << energyMax_;
 
   // Book histograms
   edm::Service<TFileService> tfile;
 
   if (!tfile.isAvailable())
     throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                       << "please add it to config file";
   // Histograms for Hits
   std::string types[4] = {"total", "by dE/dx", "by delta-ray", "by bremms"};
   char name[20], title[80];
   for (int i = 0; i < 4; i++) {
     sprintf(name, "Hits%d", i);
     sprintf(title, "Number of hits (%s)", types[i].c_str());
     meNHit_[i] = tfile->make<TH1F>(name, title, 5000, 0., 5000.);
     meNHit_[i]->GetXaxis()->SetTitle(title);
     meNHit_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "E1T0%d", i);
     sprintf(title, "E1 (Loss %s) in GeV", types[i].c_str());
     meE1T0_[i] = tfile->make<TH1F>(name, title, 5000, 0, energyMax_);
     meE1T0_[i]->GetXaxis()->SetTitle(title);
     meE1T0_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "E9T0%d", i);
     sprintf(title, "E9 (Loss %s) in GeV", types[i].c_str());
     meE9T0_[i] = tfile->make<TH1F>(name, title, 5000, 0, energyMax_);
     meE9T0_[i]->GetXaxis()->SetTitle(title);
     meE9T0_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "E1T1%d", i);
     sprintf(title, "E1 (Loss %s with t < 400 ns) in GeV", types[i].c_str());
     meE1T1_[i] = tfile->make<TH1F>(name, title, 5000, 0, energyMax_);
     meE1T1_[i]->GetXaxis()->SetTitle(title);
     meE1T1_[i]->GetYaxis()->SetTitle("Events");
     sprintf(name, "E9T1%d", i);
     sprintf(title, "E9 (Loss %s with t < 400 ns) in GeV", types[i].c_str());
     meE9T1_[i] = tfile->make<TH1F>(name, title, 5000, 0, energyMax_);
     meE9T1_[i]->GetXaxis()->SetTitle(title);
     meE9T1_[i]->GetYaxis()->SetTitle("Events");
   }
   sprintf(name, "PDGType");
   sprintf(title, "PDG ID of first level secondary");
   mType_ = tfile->make<TH1I>(name, title, 5000, -2500, 2500);
   mType_->GetXaxis()->SetTitle(title);
   mType_->GetYaxis()->SetTitle("Tracks");
 }

◆ ~XtalDedxAnalysis()

XtalDedxAnalysis::~XtalDedxAnalysis ( )

inlineoverride

Definition at line 39 of file XtalDedxAnalysis.cc.

39 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 126 of file XtalDedxAnalysis.cc.

References analyzeHits(), MillePedeFileConverter_cfg::e, edm::HandleBase::isValid(), tok_calo_, tok_tk_, and tok_vtx_.

                                                                        {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis::Run = " << e.id().run() << " Event = " << e.id().event();
 #endif
   std::vector<PCaloHit> caloHits;
   const edm::Handle<edm::PCaloHitContainer> &pCaloHits = e.getHandle(tok_calo_);
 
   const edm::Handle<edm::SimTrackContainer> &simTk = e.getHandle(tok_tk_);
   const edm::Handle<edm::SimVertexContainer> &simVtx = e.getHandle(tok_vtx_);
 
   if (pCaloHits.isValid()) {
     caloHits.insert(caloHits.end(), pCaloHits->begin(), pCaloHits->end());
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis: Hit buffer " << caloHits.size();
 #endif
     analyzeHits(caloHits, simTk, simVtx);
   }
 }

◆ analyzeHits()

void XtalDedxAnalysis::analyzeHits	(	std::vector< PCaloHit > &	hits,
		const edm::Handle< edm::SimTrackContainer > &	SimTk,
		const edm::Handle< edm::SimVertexContainer > &	SimVtx
	)

protected

Definition at line 145 of file XtalDedxAnalysis.cc.

References hcalRecHitTable_cff::energy, hfClusterShapes_cfi::hits, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::iv, meE1T0_, meE1T1_, meE9T0_, meE9T1_, meNHit_, mType_, class-composition::parent, and hcalRecHitTable_cff::time.

Referenced by analyze().

                                                                                    {
   edm::SimTrackContainer::const_iterator simTrkItr;
   int nHit = hits.size();
   double e10[4], e90[4], e11[4], e91[4], hit[4];
   for (int i = 0; i < 4; i++)
     e10[i] = e90[i] = e11[i] = e91[i] = hit[i] = 0;
   for (int i = 0; i < nHit; i++) {
     double energy = hits[i].energy();
     double time = hits[i].time();
     unsigned int id_ = hits[i].id();
     int trackID = hits[i].geantTrackId();
     int type = 1;
     for (simTrkItr = SimTk->begin(); simTrkItr != SimTk->end(); simTrkItr++) {
       if (trackID == (int)(simTrkItr->trackId())) {
         int thePID = simTrkItr->type();
         if (thePID == 11)
           type = 2;
         else if (thePID != -13 && thePID != 13)
           type = 3;
         break;
       }
     }
     hit[0]++;
     hit[type]++;
     e90[0] += energy;
     e90[type] += energy;
     if (time < 400) {
       e91[0] += energy;
       e91[type] += energy;
     }
     if (id_ == 22) {
       e10[0] += energy;
       e10[type] += energy;
       if (time < 400) {
         e11[0] += energy;
         e11[type] += energy;
       }
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis:Hit[" << i << "] ID " << id_ << " E " << energy
                                           << " time " << time << " track " << trackID << " type " << type;
 #endif
   }
   for (int i = 0; i < 4; i++) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis:Type(" << i << ") Hit " << hit[i] << " E10 " << e10[i]
                                           << " E11 " << e11[i] << " E90 " << e90[i] << " E91 " << e91[i];
 #endif
     meNHit_[i]->Fill(hit[i]);
     meE1T0_[i]->Fill(e10[i]);
     meE9T0_[i]->Fill(e90[i]);
     meE1T1_[i]->Fill(e11[i]);
     meE9T1_[i]->Fill(e91[i]);
   }
 
   // Type of the secondary (coming directly from a generator level track)
   int nvtx = 0, k1 = 0;
   edm::SimVertexContainer::const_iterator simVtxItr;
   for (simVtxItr = SimVtx->begin(); simVtxItr != SimVtx->end(); simVtxItr++)
     nvtx++;
 #ifdef EDM_ML_DEBUG
   int ntrk = 0;
   for (simTrkItr = SimTk->begin(); simTrkItr != SimTk->end(); simTrkItr++)
     ntrk++;
   edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis: " << ntrk << " tracks and " << nvtx << " vertices";
 #endif
   for (simTrkItr = SimTk->begin(); simTrkItr != SimTk->end(); simTrkItr++, ++k1) {
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("CherenkovAnalysis") << "Track " << k1 << " PDGId " << simTrkItr->type() << " Vertex ID "
                                           << simTrkItr->vertIndex() << " Generator " << simTrkItr->noGenpart();
 #endif
     if (simTrkItr->noGenpart()) {              // This is a secondary
       int vertIndex = simTrkItr->vertIndex();  // Vertex index of origin
       if (vertIndex >= 0 && vertIndex < nvtx) {
         simVtxItr = SimVtx->begin();
         for (int iv = 0; iv < vertIndex; iv++)
           simVtxItr++;
         int parent = simVtxItr->parentIndex(), k2 = 0;
         for (edm::SimTrackContainer::const_iterator trkItr = SimTk->begin(); trkItr != SimTk->end(); trkItr++, ++k2) {
 #ifdef EDM_ML_DEBUG
           edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis::Track " << k2 << " ID " << trkItr->trackId()
                                                 << " (" << parent << ")  Generator " << trkItr->noGenpart();
 #endif
           if ((int)trkItr->trackId() == parent) {  // Parent track
             if (!trkItr->noGenpart()) {            // Generator level
 #ifdef EDM_ML_DEBUG
               edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis::Track found with ID " << simTrkItr->type();
 #endif
               mType_->Fill(simTrkItr->type());
             }
             break;
           }
         }
       }
     }
   }
 }

◆ beginJob()

void XtalDedxAnalysis::beginJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 44 of file XtalDedxAnalysis.cc.

44 {}

◆ endJob()

void XtalDedxAnalysis::endJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 46 of file XtalDedxAnalysis.cc.

46 {}

◆ fillDescriptions()

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

static

Definition at line 118 of file XtalDedxAnalysis.cc.

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

                                                                                   {
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("caloHitSource", edm::InputTag("g4SimHits", "EcalHitsEB"));
   desc.add<std::string>("moduleLabelTk", "g4SimHits");
   desc.add<double>("energyMax", 2.0);
   descriptions.add("xtalDedxAnalysis", desc);
 }

Member Data Documentation

◆ caloHitSource_

const edm::InputTag XtalDedxAnalysis::caloHitSource_

private

Definition at line 53 of file XtalDedxAnalysis.cc.

Referenced by XtalDedxAnalysis().

◆ energyMax_

const double XtalDedxAnalysis::energyMax_

private

Definition at line 55 of file XtalDedxAnalysis.cc.

Referenced by XtalDedxAnalysis().

◆ meE1T0_

TH1F * XtalDedxAnalysis::meE1T0_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

◆ meE1T1_

TH1F * XtalDedxAnalysis::meE1T1_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

◆ meE9T0_

TH1F * XtalDedxAnalysis::meE9T0_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

◆ meE9T1_

TH1F * XtalDedxAnalysis::meE9T1_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

◆ meNHit_

TH1F* XtalDedxAnalysis::meNHit_[4]

private

Definition at line 60 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

◆ mType_

TH1I* XtalDedxAnalysis::mType_

private

Definition at line 61 of file XtalDedxAnalysis.cc.

Referenced by analyzeHits(), and XtalDedxAnalysis().

◆ simTkLabel_

const std::string XtalDedxAnalysis::simTkLabel_

private

Definition at line 54 of file XtalDedxAnalysis.cc.

Referenced by XtalDedxAnalysis().

◆ tok_calo_

const edm::EDGetTokenT<edm::PCaloHitContainer> XtalDedxAnalysis::tok_calo_

private

Definition at line 56 of file XtalDedxAnalysis.cc.

Referenced by analyze().

◆ tok_tk_

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

private

Definition at line 57 of file XtalDedxAnalysis.cc.

Referenced by analyze().

◆ tok_vtx_

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

private

Definition at line 58 of file XtalDedxAnalysis.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ XtalDedxAnalysis()

◆ ~XtalDedxAnalysis()

Member Function Documentation

◆ analyze()

◆ analyzeHits()

◆ beginJob()

◆ endJob()

◆ fillDescriptions()

Member Data Documentation

◆ caloHitSource_

◆ energyMax_

◆ meE1T0_

◆ meE1T1_

◆ meE9T0_

◆ meE9T1_

◆ meNHit_

◆ mType_

◆ simTkLabel_

◆ tok_calo_

◆ tok_tk_

◆ tok_vtx_