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

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

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	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, overlapproblemtsosanalyzer_cfi::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;
           }
         }
       }
     }
   }
 }