Inheritance diagram for HGCalBHValidation:

Public Member Functions
	HGCalBHValidation (const edm::ParameterSet &ps)

	~HGCalBHValidation () override=default

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 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 (edm::Event const &, edm::EventSetup const &) override

void	beginJob () override

void	beginRun (edm::Run const &, edm::EventSetup const &) 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)

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	digis_

edm::Service< TFileService >	fs_

const std::string	g4Label_

TH2D *	hdi2Oc_

TH2D *	hdi3Oc_

TH1D *	hdigEn_

TH1D *	hdigLn_

TH2D *	hdigOc_

const std::string	hits_

TH2D *	hsi2Oc_

TH2D *	hsi3Oc_

TH1D *	hsimE1_

TH1D *	hsimE2_

TH1D *	hsimLn_

TH2D *	hsimOc_

TH1D *	hsimTm_

const int	iSample_

const int	ringMax_

const double	threshold_

const edm::EDGetTokenT< HGCalDigiCollection >	tok_digi_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hits_

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 32 of file HGCalBHValidation.cc.

Constructor & Destructor Documentation

◆ HGCalBHValidation()

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

explicit

Definition at line 61 of file HGCalBHValidation.cc.

References digis_, g4Label_, hits_, ProducerED_cfi::InputTag, iSample_, TFileService::kSharedResource, and threshold_.

     : g4Label_(ps.getParameter<std::string>("ModuleLabel")),
       hits_((ps.getParameter<std::string>("HitCollection"))),
       digis_(ps.getParameter<edm::InputTag>("DigiCollection")),
       iSample_(ps.getParameter<int>("Sample")),
       threshold_(ps.getParameter<double>("Threshold")),
       tok_hits_(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, hits_))),
       tok_digi_(consumes<HGCalDigiCollection>(digis_)),
       ringMax_(50) {
   usesResource(TFileService::kSharedResource);
 
   edm::LogVerbatim("HGCalValidation") << "HGCalBHValidation::Input for SimHit:" << edm::InputTag(g4Label_, hits_)
                                       << "  Digits:" << digis_ << "  Sample: " << iSample_ << "  Threshold "
                                       << threshold_;
 }

◆ ~HGCalBHValidation()

HGCalBHValidation::~HGCalBHValidation ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 107 of file HGCalBHValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), TauDecayModes::dec, hcalRecHitTable_cff::depth, das-up-to-nevents::df, MillePedeFileConverter_cfg::e, edm::SortedCollection< T, SORT >::end(), hcalRecHitTable_cff::energy, hdi2Oc_, hdi3Oc_, hdigEn_, hdigLn_, hdigOc_, DetId::HGCalHSc, hsi2Oc_, hsi3Oc_, hsimE1_, hsimE2_, hsimLn_, hsimOc_, hsimTm_, mps_fire::i, l1ctLayer2EG_cff::id, HGCScintillatorDetId::iphi(), iSample_, edm::HandleBase::isValid(), ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, HGCScintillatorDetId::layer(), phi, relativeConstraints::ring, HGCScintillatorDetId::ring(), edm::SortedCollection< T, SORT >::size(), threshold_, hcalRecHitTable_cff::time, tok_digi_, tok_hits_, and HGCScintillatorDetId::zside().

                                                                        {
   edm::LogVerbatim("HGCalValidation") << "HGCalBHValidation:Run = " << e.id().run() << " Event = " << e.id().event();
 
   //SimHits
   const edm::Handle<edm::PCaloHitContainer>& hitsHE = e.getHandle(tok_hits_);
   edm::LogVerbatim("HGCalValidation") << "HGCalBHValidation.: PCaloHitContainer"
                                       << " obtained with flag " << hitsHE.isValid();
   if (hitsHE.isValid()) {
     edm::LogVerbatim("HGCalValidation") << "HGCalBHValidation: PCaloHit buffer " << hitsHE->size();
     unsigned i(0);
     std::map<unsigned int, double> map_try;
     for (edm::PCaloHitContainer::const_iterator it = hitsHE->begin(); it != hitsHE->end(); ++it) {
       double energy = it->energy();
       double time = it->time();
       unsigned int id = it->id();
       if (DetId(id).det() == DetId::HGCalHSc) {
         double ring = 0.01 + HGCScintillatorDetId(id).ring();
         double phi = HGCScintillatorDetId(id).iphi() - 0.01;
         int lay = HGCScintillatorDetId(id).layer();
         double ring1 = ring * HGCScintillatorDetId(id).zside();
         hsi2Oc_->Fill(ring1, phi, energy);
         hsimE1_->Fill(energy);
         hsimTm_->Fill(time, energy);
         hsimOc_->Fill(ring1, phi, energy);
         hsimLn_->Fill(lay, energy);
         hsi3Oc_->Fill(ring1, lay, energy);
         double ensum(0);
         if (map_try.count(id) != 0)
           ensum = map_try[id];
         ensum += energy;
         map_try[id] = ensum;
         edm::LogVerbatim("HGCalValidation") << "HGCalBHHit[" << ++i << "] ID " << std::hex << " " << id << std::dec
                                             << " " << HGCScintillatorDetId(id) << " E " << energy << " time " << time;
       }
     }
     for (std::map<unsigned int, double>::iterator itr = map_try.begin(); itr != map_try.end(); ++itr) {
       hsimE2_->Fill((*itr).second);
     }
   }
 
   //Digits
   unsigned int kount(0);
   const edm::Handle<HGCalDigiCollection>& hecoll = e.getHandle(tok_digi_);
   edm::LogVerbatim("HGCalValidation") << "HGCalBHValidation.: "
                                       << "HGCalDigiCollection obtained with"
                                       << " flag " << hecoll.isValid();
   if (hecoll.isValid()) {
     edm::LogVerbatim("HGCalValidation") << "HGCalBHValidation: HGCalDigi "
                                         << "buffer " << hecoll->size();
     for (HGCalDigiCollection::const_iterator it = hecoll->begin(); it != hecoll->end(); ++it) {
       const HGCalDataFrame& df(*it);
       double energy = df[iSample_].data();
       if (DetId(df.id()).det() == DetId::HGCalHSc) {
         HGCScintillatorDetId cell(df.id());
         int depth = cell.layer();
         if (energy > threshold_) {
           double ring = 0.01 + cell.ring();
           double phi = cell.iphi() - 0.01;
           double ring1 = cell.zside() * ring;
           hdi2Oc_->Fill(ring1, phi);
           hdigEn_->Fill(energy);
           hdigOc_->Fill(ring1, phi);
           hdigLn_->Fill(depth);
           hdi3Oc_->Fill(ring1, depth);
           edm::LogVerbatim("HGCalValidation")
               << "HGCalBHDigit[" << ++kount << "] ID " << cell << " E " << energy << ":" << (energy > threshold_);
         }
       }
     }
   }
 }

◆ beginJob()

void HGCalBHValidation::beginJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 40 of file HGCalBHValidation.cc.

40 {}

◆ beginRun()

void HGCalBHValidation::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	es
	)

overrideprotected

Definition at line 87 of file HGCalBHValidation.cc.

References fs_, hdi2Oc_, hdi3Oc_, hdigEn_, hdigLn_, hdigOc_, hsi2Oc_, hsi3Oc_, hsimE1_, hsimE2_, hsimLn_, hsimOc_, hsimTm_, TFileService::make(), and ringMax_.

                                                                        {
   edm::LogVerbatim("HGCalValidation") << "HGCalBHValidation::Maximum Number of"
                                       << " ring sectors:" << ringMax_ << "\nHitsValidationHcal::Booking the Histograms";
 
   //Histograms for Sim Hits
   hsimE1_ = fs_->make<TH1D>("SimHitEn1", "Sim Hit Energy", 1000, 0.0, 1.0);
   hsimE2_ = fs_->make<TH1D>("SimHitEn2", "Sim Hit Energy", 1000, 0.0, 1.0);
   hsimTm_ = fs_->make<TH1D>("SimHitTime", "Sim Hit Time", 1000, 0.0, 500.0);
   hsimLn_ = fs_->make<TH1D>("SimHitLong", "Sim Hit Long. Profile", 50, 0.0, 25.0);
   hsimOc_ = fs_->make<TH2D>("SimHitOccup", "Sim Hit Occupnacy", 2 * ringMax_ + 1, -ringMax_, ringMax_ + 1, 360, 0, 360);
   hsi2Oc_ = fs_->make<TH2D>("SimHitOccu2", "Sim Hit Occupnacy", 2 * ringMax_ + 1, -ringMax_, ringMax_ + 1, 360, 0, 360);
   hsi3Oc_ = fs_->make<TH2D>("SimHitOccu3", "Sim Hit Occupnacy", 2 * ringMax_ + 1, -ringMax_, ringMax_ + 1, 50, 0, 25);
   //Histograms for Digis
   hdigEn_ = fs_->make<TH1D>("DigiEnergy", "Digi ADC Sample", 1000, 0.0, 1000.0);
   hdigLn_ = fs_->make<TH1D>("DigiLong", "Digi Long. Profile", 50, 0.0, 25.0);
   hdigOc_ = fs_->make<TH2D>("DigiOccup", "Digi Occupnacy", 2 * ringMax_ + 1, -ringMax_, ringMax_ + 1, 360, 0, 360);
   hdi2Oc_ = fs_->make<TH2D>("DigiOccu2", "Digi Occupnacy", 2 * ringMax_ + 1, -ringMax_, ringMax_ + 1, 360, 0, 360);
   hdi3Oc_ = fs_->make<TH2D>("DigiOccu3", "Digi Occupnacy", 2 * ringMax_ + 1, -ringMax_, ringMax_ + 1, 50, 0, 25);
 }

◆ endRun()

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

inlineoverrideprotected

Definition at line 42 of file HGCalBHValidation.cc.

42 {}

◆ fillDescriptions()

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

static

Definition at line 77 of file HGCalBHValidation.cc.

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

                                                                                    {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("ModuleLabel", "g4SimHits");
   desc.add<std::string>("HitCollection", "HGCHitsHEback");
   desc.add<edm::InputTag>("DigiCollection", edm::InputTag("simHGCalUnsuppressedDigis", "HEback"));
   desc.add<int>("Sample", 5);
   desc.add<double>("Threshold", 15.0);
   descriptions.add("hgcalBHAnalysis", desc);
 }

Member Data Documentation

◆ digis_

const edm::InputTag HGCalBHValidation::digis_

private

Definition at line 48 of file HGCalBHValidation.cc.

Referenced by HGCalBHValidation().

◆ fs_

edm::Service<TFileService> HGCalBHValidation::fs_

private

Definition at line 46 of file HGCalBHValidation.cc.

Referenced by beginRun().

◆ g4Label_

const std::string HGCalBHValidation::g4Label_

private

Definition at line 47 of file HGCalBHValidation.cc.

Referenced by HGCalBHValidation().

◆ hdi2Oc_

TH2D * HGCalBHValidation::hdi2Oc_

private

Definition at line 58 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hdi3Oc_

TH2D * HGCalBHValidation::hdi3Oc_

private

Definition at line 58 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hdigEn_

TH1D * HGCalBHValidation::hdigEn_

private

Definition at line 56 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hdigLn_

TH1D * HGCalBHValidation::hdigLn_

private

Definition at line 56 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hdigOc_

TH2D* HGCalBHValidation::hdigOc_

private

Definition at line 58 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hits_

const std::string HGCalBHValidation::hits_

private

Definition at line 47 of file HGCalBHValidation.cc.

Referenced by HGCalBHValidation().

◆ hsi2Oc_

TH2D * HGCalBHValidation::hsi2Oc_

private

Definition at line 57 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hsi3Oc_

TH2D * HGCalBHValidation::hsi3Oc_

private

Definition at line 57 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hsimE1_

TH1D* HGCalBHValidation::hsimE1_

private

Definition at line 55 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hsimE2_

TH1D * HGCalBHValidation::hsimE2_

private

Definition at line 55 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hsimLn_

TH1D* HGCalBHValidation::hsimLn_

private

Definition at line 56 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hsimOc_

TH2D* HGCalBHValidation::hsimOc_

private

Definition at line 57 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ hsimTm_

TH1D * HGCalBHValidation::hsimTm_

private

Definition at line 55 of file HGCalBHValidation.cc.

Referenced by analyze(), and beginRun().

◆ iSample_

const int HGCalBHValidation::iSample_

private

Definition at line 49 of file HGCalBHValidation.cc.

Referenced by analyze(), and HGCalBHValidation().

◆ ringMax_

const int HGCalBHValidation::ringMax_

private

Definition at line 53 of file HGCalBHValidation.cc.

Referenced by beginRun().

◆ threshold_

const double HGCalBHValidation::threshold_

private

Definition at line 50 of file HGCalBHValidation.cc.

Referenced by analyze(), and HGCalBHValidation().

◆ tok_digi_

const edm::EDGetTokenT<HGCalDigiCollection> HGCalBHValidation::tok_digi_

private

Definition at line 52 of file HGCalBHValidation.cc.

Referenced by analyze().

◆ tok_hits_

const edm::EDGetTokenT<edm::PCaloHitContainer> HGCalBHValidation::tok_hits_

private

Definition at line 51 of file HGCalBHValidation.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGCalBHValidation()

◆ ~HGCalBHValidation()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ beginRun()

◆ endRun()

◆ fillDescriptions()

Member Data Documentation

◆ digis_

◆ fs_

◆ g4Label_

◆ hdi2Oc_

◆ hdi3Oc_

◆ hdigEn_

◆ hdigLn_

◆ hdigOc_

◆ hits_

◆ hsi2Oc_

◆ hsi3Oc_

◆ hsimE1_

◆ hsimE2_

◆ hsimLn_

◆ hsimOc_

◆ hsimTm_

◆ iSample_

◆ ringMax_

◆ threshold_

◆ tok_digi_

◆ tok_hits_