Inheritance diagram for XtalDedxAnalysis:

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

	~XtalDedxAnalysis () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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

std::vector< ESProxyIndex > 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 &	operator= (EDConsumerBase &&)=default

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

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)

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

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

void	beginJob () override

void	endJob () override

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

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

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

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

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

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

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

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

Private Attributes
edm::InputTag	caloHitSource_

TH1F *	meE1T0_ [4]

TH1F *	meE1T1_ [4]

TH1F *	meE9T0_ [4]

TH1F *	meE9T1_ [4]

TH1F *	meNHit_ [4]

TH1I *	mType_

std::string	simTkLabel_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_calo_

edm::EDGetTokenT< edm::SimTrackContainer >	tok_tk_

edm::EDGetTokenT< edm::SimVertexContainer >	tok_vtx_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Detailed Description

Definition at line 37 of file XtalDedxAnalysis.cc.

Constructor & Destructor Documentation

◆ XtalDedxAnalysis()

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

explicit

Definition at line 63 of file XtalDedxAnalysis.cc.

                                                             {
   caloHitSource_ = ps.getParameter<edm::InputTag>("caloHitSource");
   simTkLabel_ = ps.getUntrackedParameter<std::string>("moduleLabelTk", "g4SimHits");
   double energyMax = ps.getParameter<double>("EnergyMax");
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("CherenkovAnalysis") << "XtalDedxAnalysis::Source " << caloHitSource_ << " Track Label " << simTkLabel_
                                     << " Energy Max " << energyMax;
 #endif
   // register for data access
   tok_calo_ = consumes<edm::PCaloHitContainer>(caloHitSource_);
   tok_tk_ = consumes<edm::SimTrackContainer>(edm::InputTag(simTkLabel_));
   tok_vtx_ = consumes<edm::SimVertexContainer>(edm::InputTag(simTkLabel_));
  
   // 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");
 }

References caloHitSource_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mps_fire::i, meE1T0_, meE1T1_, meE9T0_, meE9T1_, meNHit_, mType_, Skims_PA_cff::name, simTkLabel_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, runGCPTkAlMap::title, tok_calo_, tok_tk_, and tok_vtx_.

◆ ~XtalDedxAnalysis()

XtalDedxAnalysis::~XtalDedxAnalysis ( )

inlineoverride

Definition at line 40 of file XtalDedxAnalysis.cc.

40 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Implements edm::EDAnalyzer.

Definition at line 119 of file XtalDedxAnalysis.cc.

                                                                        {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("CherenkovAnalysis") << "XtalDedxAnalysis::Run = " << e.id().run() << " Event = " << e.id().event();
 #endif
   std::vector<PCaloHit> caloHits;
   edm::Handle<edm::PCaloHitContainer> pCaloHits;
   e.getByToken(tok_calo_, pCaloHits);
  
   std::vector<SimTrack> theSimTracks;
   edm::Handle<edm::SimTrackContainer> simTk;
   e.getByToken(tok_tk_, simTk);
  
   std::vector<SimVertex> theSimVertex;
   edm::Handle<edm::SimVertexContainer> simVtx;
   e.getByToken(tok_vtx_, simVtx);
  
   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);
   }
 }

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

◆ analyzeHits()

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

protected

Definition at line 144 of file XtalDedxAnalysis.cc.

                                                                              {
   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, ntrk = 0, k1 = 0;
   edm::SimVertexContainer::const_iterator simVtxItr;
   for (simTrkItr = SimTk->begin(); simTrkItr != SimTk->end(); simTrkItr++)
     ntrk++;
   for (simVtxItr = SimVtx->begin(); simVtxItr != SimVtx->end(); simVtxItr++)
     nvtx++;
 #ifdef EDM_ML_DEBUG
   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;
           }
         }
       }
     }
   }
 }