Inheritance diagram for HFPMTHitAnalyzer:

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

	~HFPMTHitAnalyzer () 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 final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () 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 &&)=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

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

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)

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

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

void	beginJob () override

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

void	endJob () override

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

Private Attributes
int	event_no

std::string	g4Label

TH1F *	h_HFDepHit

TH1F *	h_HFEta [3]

TH1F *	h_HFPhi [3]

std::string	hcalHits

TH1F *	hHF12_time [3]

TH1F *	hHF12_time_Ewt [3]

TH1F *	hHF1_time [3]

TH1F *	hHF1_time_Ewt [3]

TH1F *	hHF2_time [3]

TH1F *	hHF2_time_Ewt [3]

TH1F *	hHF_e_1 [3]

TH1F *	hHF_e_12 [3]

TH1F *	hHF_e_2 [3]

TH1F *	hHF_em_1 [3]

TH1F *	hHF_em_12 [3]

TH1F *	hHF_em_2 [3]

TH1F *	hHF_had_1 [3]

TH1F *	hHF_had_12 [3]

TH1F *	hHF_had_2 [3]

TH1F *	hHF_MC_e

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_calo_

edm::EDGetTokenT< edm::HepMCProduct >	tok_evt_

edm::EDGetTokenT< edm::SimTrackContainer >	tok_track_

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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)

Detailed Description

Definition at line 34 of file HFPMTHitAnalyzer.cc.

Constructor & Destructor Documentation

◆ HFPMTHitAnalyzer()

HFPMTHitAnalyzer::HFPMTHitAnalyzer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 68 of file HFPMTHitAnalyzer.cc.

                                                                  {
   usesResource(TFileService::kSharedResource);
  
   tok_evt_ = consumes<edm::HepMCProduct>(
       edm::InputTag(iConfig.getUntrackedParameter<std::string>("SourceLabel", "VtxSmeared")));
   g4Label = iConfig.getUntrackedParameter<std::string>("ModuleLabel", "g4SimHits");
   hcalHits = iConfig.getUntrackedParameter<std::string>("HitCollection", "HcalHits");
  
   tok_calo_ = consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label, hcalHits));
   tok_track_ = consumes<edm::SimTrackContainer>(edm::InputTag(g4Label));
 }

References g4Label, edm::ParameterSet::getUntrackedParameter(), hcalHits, TFileService::kSharedResource, AlCaHLTBitMon_QueryRunRegistry::string, tok_calo_, tok_evt_, and tok_track_.

◆ ~HFPMTHitAnalyzer()

HFPMTHitAnalyzer::~HFPMTHitAnalyzer ( )

inlineoverride

Definition at line 37 of file HFPMTHitAnalyzer.cc.

37 {}

Member Function Documentation

◆ analyze()

void HFPMTHitAnalyzer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 203 of file HFPMTHitAnalyzer.cc.

                                                                                   {
   ++event_no;
   if (event_no % 500 == 0)
     std::cout << "Event # " << event_no << " processed.\n";
  
   std::vector<PCaloHit> caloHits;
   edm::Handle<edm::PCaloHitContainer> hitsHcal;
   iEvent.getByToken(tok_calo_, hitsHcal);
   std::vector<SimTrack> simTracks;
   edm::Handle<edm::SimTrackContainer> Tracks;
   iEvent.getByToken(tok_track_, Tracks);
  
   edm::Handle<edm::HepMCProduct> EvtHandle;
   iEvent.getByToken(tok_evt_, EvtHandle);
   const HepMC::GenEvent *myGenEvent = EvtHandle->GetEvent();
  
   float orig_energy = 0;
   for (HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end();
        ++p) {
     orig_energy = (*p)->momentum().e();
     hHF_MC_e->Fill(orig_energy);
   }
  
   caloHits.insert(caloHits.end(), hitsHcal->begin(), hitsHcal->end());
   analyzeHits(caloHits, *Tracks);
 }

References analyzeHits(), gather_cfg::cout, event_no, edm::HepMCProduct::GetEvent(), hHF_MC_e, iEvent, AlCaHLTBitMon_ParallelJobs::p, TrackCandidateProducer_cfi::simTracks, tok_calo_, tok_evt_, tok_track_, and gather_cfg::Tracks.

◆ analyzeHits()

void HFPMTHitAnalyzer::analyzeHits	(	std::vector< PCaloHit > &	hits,
		const std::vector< SimTrack > &	tracks1
	)

private

Definition at line 230 of file HFPMTHitAnalyzer.cc.

                                                                                                 {
   int nHit = hits.size();
   float energy1[3], energy2[3], energy12[3];
   float em1[3], had1[3], em2[3], had2[3], em12[3], had12[3];
   for (int i = 0; i < 3; i++) {
     energy1[i] = energy2[i] = energy12[i] = 0;
     em1[i] = em2[i] = em12[i] = 0;
     had1[i] = had2[i] = had12[i] = 0;
   }
   edm::LogVerbatim("HFShower") << "HFPMTHitAnalyser::Entry " << event_no << " with " << nHit << " hits";
   for (int i = 0; i < nHit; i++) {
     double energy = hits[i].energy();
     double em = hits[i].energyEM();
     double had = hits[i].energyHad();
     double time = hits[i].time();
     uint32_t id_ = hits[i].id();
     uint16_t pmtHit = hits[i].depth();
     uint16_t depthX = pmtHit;
     edm::LogVerbatim("HFShower") << "HFPMTHitAnalyser::Hit " << i << " ID " << std::hex << id_ << std::dec << " PMT "
                                  << pmtHit << " E (e|h|t) " << em << " " << had << " " << energy << " Time " << time;
  
     HcalDetId id = HcalDetId(id_);
     int det = id.det();
     int subdet = id.subdet();
     int depth = id.depth();
     if (pmtHit != 0)
       pmtHit = 1;
  
     if (det == 4) {
       if (subdet == static_cast<int>(HcalForward)) {
         h_HFDepHit->Fill(double(depth + 10 * depthX));
         if (depthX > 0) {
           int ieta = id.ietaAbs();
           int iphi = id.iphi();
           if (depth != 1) {
             ieta += 50;
             iphi += 50;
           }
           if (depthX <= 3) {
             h_HFEta[depthX - 1]->Fill(double(ieta));
             h_HFPhi[depthX - 1]->Fill(double(iphi));
           }
         }
         if (depth == 1) {
           energy1[pmtHit] += energy;
           energy12[pmtHit] += energy;
           em1[pmtHit] += em;
           em12[pmtHit] += em;
           had1[pmtHit] += had;
           had12[pmtHit] += had;
           energy1[2] += energy;
           energy12[2] += energy;
           em1[2] += em;
           em12[2] += em;
           had1[2] += had;
           had12[2] += had;
           hHF1_time[pmtHit]->Fill(time);
           hHF1_time[2]->Fill(time);
           hHF1_time_Ewt[pmtHit]->Fill(time, energy);
           hHF1_time_Ewt[2]->Fill(time, energy);
         }
         if (depth == 2) {
           energy2[pmtHit] += energy;
           energy12[pmtHit] += energy;
           em2[pmtHit] += em;
           em12[pmtHit] += em;
           had2[pmtHit] += had;
           had12[pmtHit] += had;
           energy2[2] += energy;
           energy12[2] += energy;
           em2[2] += em;
           em12[2] += em;
           had2[2] += had;
           had12[2] += had;
           hHF2_time[pmtHit]->Fill(time);
           hHF2_time[2]->Fill(time);
           hHF2_time_Ewt[pmtHit]->Fill(time, energy);
           hHF2_time_Ewt[2]->Fill(time, energy);
         }
       }
     }
   }
   for (int i = 0; i < 3; i++) {
     hHF_e_1[i]->Fill(energy1[i]);
     hHF_e_2[i]->Fill(energy2[i]);
     hHF_e_12[i]->Fill(energy12[i]);
     hHF_em_1[i]->Fill(em1[i]);
     hHF_em_2[i]->Fill(em2[i]);
     hHF_em_12[i]->Fill(em12[i]);
     hHF_had_1[i]->Fill(had1[i]);
     hHF_had_2[i]->Fill(had2[i]);
     hHF_had_12[i]->Fill(had12[i]);
     edm::LogVerbatim("HFShower") << "HFPMTHitAnalyser:: Type " << i << " Energy 1|2| " << energy1[i] << " "
                                  << energy2[i] << " " << energy12[i] << " EM Energy 1|2| " << em1[i] << " " << em2[i]
                                  << " " << em12[i] << " Had Energy 1|2| " << had1[i] << " " << had2[i] << " "
                                  << had12[i];
   }
 }

References TauDecayModes::dec, LEDCalibrationChannels::depth, HCALHighEnergyHPDFilter_cfi::energy, event_no, h_HFDepHit, h_HFEta, h_HFPhi, HcalForward, hHF1_time, hHF1_time_Ewt, hHF2_time, hHF2_time_Ewt, hHF_e_1, hHF_e_12, hHF_e_2, hHF_em_1, hHF_em_12, hHF_em_2, hHF_had_1, hHF_had_12, hHF_had_2, hfClusterShapes_cfi::hits, mps_fire::i, LEDCalibrationChannels::ieta, LEDCalibrationChannels::iphi, and protons_cff::time.

Referenced by analyze().

◆ beginJob()

void HFPMTHitAnalyzer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 88 of file HFPMTHitAnalyzer.cc.

                                 {
   event_no = 0;
   char name[20], title[120], sub[11];
  
   edm::Service<TFileService> fs;
   if (!fs.isAvailable())
     throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                       << "please add it to config file";
  
   TFileDirectory HFHitsDir = fs->mkdir("HFPMTHits");
   h_HFDepHit = HFHitsDir.make<TH1F>("Hit20", "Depths in HF", 50, 0., 50.);
   h_HFDepHit->GetXaxis()->SetTitle("Depths in HF");
   for (int i = 0; i < 3; i++) {
     if (i == 0)
       sprintf(sub, "(PMT)");
     else if (i == 1)
       sprintf(sub, "(Bundle)");
     else
       sprintf(sub, "(Jungle)");
     sprintf(name, "Eta%d", i);
     sprintf(title, "Eta Index of hits in %s", sub);
     h_HFEta[i] = HFHitsDir.make<TH1F>(name, title, 100, 0, 100.);
     h_HFEta[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Phi%d", i);
     sprintf(title, "Phi Index of hits in %s", sub);
     h_HFPhi[i] = HFHitsDir.make<TH1F>(name, title, 100, 0, 100.);
     h_HFPhi[i]->GetXaxis()->SetTitle(title);
   }
   TFileDirectory HFSourcePart = fs->mkdir("HFMCinfo");
   hHF_MC_e = HFSourcePart.make<TH1F>("MCEnergy", "Energy of Generated Particle", 1000, 0., 500.);
   hHF_MC_e->GetXaxis()->SetTitle("Energy of Generated Particle");
  
   //energy Histograms
   TFileDirectory HFPCaloHitEnergyDir = fs->mkdir("HFPCaloHitEnergy");
   for (int i = 0; i < 3; i++) {
     if (i == 0)
       sprintf(sub, "(Absorber)");
     else if (i == 1)
       sprintf(sub, "(PMT)");
     else
       sprintf(sub, "(All)");
     sprintf(name, "Energy1%d", i);
     sprintf(title, "Energy in depth 1 %s", sub);
     hHF_e_1[i] = HFPCaloHitEnergyDir.make<TH1F>(name, title, 1000, 0., 500.);
     hHF_e_1[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Energy2%d", i);
     sprintf(title, "Energy in depth 2 %s", sub);
     hHF_e_2[i] = HFPCaloHitEnergyDir.make<TH1F>(name, title, 1000, 0., 500.);
     hHF_e_2[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Energy12%d", i);
     sprintf(title, "Energy in depths 1,2 %s", sub);
     hHF_e_12[i] = HFPCaloHitEnergyDir.make<TH1F>(name, title, 1000, 0., 500.);
     hHF_e_12[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Em1%d", i);
     sprintf(title, "EM energy in depth 1 %s", sub);
     hHF_em_1[i] = HFPCaloHitEnergyDir.make<TH1F>(name, title, 1000, 0., 500.);
     hHF_em_1[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Em2%d", i);
     sprintf(title, "EM energy in depth 2 %s", sub);
     hHF_em_2[i] = HFPCaloHitEnergyDir.make<TH1F>(name, title, 1000, 0., 500.);
     hHF_em_2[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Em12%d", i);
     sprintf(title, "EM energy in depths 1,2 %s", sub);
     hHF_em_12[i] = HFPCaloHitEnergyDir.make<TH1F>(name, title, 1000, 0., 500.);
     hHF_em_12[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Had1%d", i);
     sprintf(title, "Had energy in depth 1 %s", sub);
     hHF_had_1[i] = HFPCaloHitEnergyDir.make<TH1F>(name, title, 1000, 0., 0.1);
     hHF_had_1[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Had2%d", i);
     sprintf(title, "Had energy in depth 2 %s", sub);
     hHF_had_2[i] = HFPCaloHitEnergyDir.make<TH1F>(name, title, 1000, 0., 0.1);
     hHF_had_2[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Had12%d", i);
     sprintf(title, "Had energy in depths 1,2 %s", sub);
     hHF_had_12[i] = HFPCaloHitEnergyDir.make<TH1F>(name, title, 1000, 0., 0.1);
     hHF_had_12[i]->GetXaxis()->SetTitle(title);
   }
  
   //Timing Histograms
   TFileDirectory HFPCaloHitTimeDir = fs->mkdir("HFPCaloHitTime");
   for (int i = 0; i < 3; i++) {
     if (i == 0)
       sprintf(sub, "(Absorber)");
     else if (i == 1)
       sprintf(sub, "(PMT)");
     else
       sprintf(sub, "(All)");
     sprintf(name, "Time1Ewt%d", i);
     sprintf(title, "Time (energy weighted) in depth 1 %s", sub);
     hHF1_time_Ewt[i] = HFPCaloHitTimeDir.make<TH1F>(name, title, 400, 0., 400.);
     hHF1_time_Ewt[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Time2Ewt%d", i);
     sprintf(title, "Time (energy weighted) in depth 2 %s", sub);
     hHF2_time_Ewt[i] = HFPCaloHitTimeDir.make<TH1F>(name, title, 400, 0., 400.);
     hHF2_time_Ewt[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Time12Ewt%d", i);
     sprintf(title, "Time (energy weighted) in depths 1,2 %s", sub);
     hHF12_time_Ewt[i] = HFPCaloHitTimeDir.make<TH1F>(name, title, 400, 0., 400.);
     hHF12_time_Ewt[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Time1%d", i);
     sprintf(title, "Time in depth 1 %s", sub);
     hHF1_time[i] = HFPCaloHitTimeDir.make<TH1F>(name, title, 400, 0., 400.);
     hHF1_time[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Time2%d", i);
     sprintf(title, "Time in depth 2 %s", sub);
     hHF2_time[i] = HFPCaloHitTimeDir.make<TH1F>(name, title, 400, 0., 400.);
     hHF2_time[i]->GetXaxis()->SetTitle(title);
     sprintf(name, "Time12%d", i);
     sprintf(title, "Time in depths 1,2 %s", sub);
     hHF12_time[i] = HFPCaloHitTimeDir.make<TH1F>(name, title, 400, 0., 400.);
     hHF12_time[i]->GetXaxis()->SetTitle(title);
   }
 }

References event_no, h_HFDepHit, h_HFEta, h_HFPhi, hHF12_time, hHF12_time_Ewt, hHF1_time, hHF1_time_Ewt, hHF2_time, hHF2_time_Ewt, hHF_e_1, hHF_e_12, hHF_e_2, hHF_em_1, hHF_em_12, hHF_em_2, hHF_had_1, hHF_had_12, hHF_had_2, hHF_MC_e, mps_fire::i, edm::Service< T >::isAvailable(), TFileDirectory::make(), TFileService::mkdir(), Skims_PA_cff::name, and runGCPTkAlMap::title.

◆ beginRun()

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

inlineoverrideprivate

Definition at line 42 of file HFPMTHitAnalyzer.cc.

42 {}

◆ endJob()

void HFPMTHitAnalyzer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 329 of file HFPMTHitAnalyzer.cc.

329 {}

◆ endRun()

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

inlineoverrideprivate

Definition at line 45 of file HFPMTHitAnalyzer.cc.

45 {}

◆ fillDescriptions()

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

static

Definition at line 80 of file HFPMTHitAnalyzer.cc.

                                                                                   {
   edm::ParameterSetDescription desc;
   desc.addUntracked<std::string>("SourceLabel", "generatorSmeared");
   desc.addUntracked<std::string>("ModuleLabel", "g4SimHits");
   desc.addUntracked<std::string>("HitCollection", "HcalHits");
   descriptions.add("HFPMTHitAnalyzer", desc);
 }