Inheritance diagram for HGcalHitCheck:

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

	~HGcalHitCheck () override

Public Member Functions inherited from edm::one::EDAnalyzer<>
	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	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 std::string	caloHitSource_

const std::string	g4Label_

const edm::ESGetToken< HGCalGeometry, IdealGeometryRecord >	geomToken_

TH1D *	h_hit1_

TH1D *	h_hit2_

std::vector< TH1D * >	h_hitF_

std::vector< TH1D * >	h_hitL_

std::vector< TH1D * >	h_hitP_

TH1D *	h_hits_

bool	histos_

const int	layers_

const std::string	nameDetector_

const std::string	nameSense_

const std::string	tag_

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

const int	verbosity_

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 31 of file HGCalHitCheck.cc.

Constructor & Destructor Documentation

◆ HGcalHitCheck()

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

Definition at line 52 of file HGCalHitCheck.cc.

References nameSense_.

     : g4Label_(ps.getParameter<std::string>("moduleLabel")),
       caloHitSource_(ps.getParameter<std::string>("caloHitSource")),
       nameSense_(ps.getParameter<std::string>("nameSense")),
       nameDetector_(ps.getParameter<std::string>("tag")),
       tag_(ps.getParameter<std::string>("nameDevice")),
       layers_(ps.getParameter<int>("layers")),
       verbosity_(ps.getParameter<int>("verbosity")),
       tok_calo_(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, caloHitSource_))),
       geomToken_(esConsumes<HGCalGeometry, IdealGeometryRecord>(edm::ESInputTag{"", nameSense_})),
       histos_(false) {
   edm::LogVerbatim("HitStudy") << "Test Hit ID for " << nameDetector_ << " using SimHits for " << nameSense_
                                << " with module Label: " << g4Label_ << "   Hits: " << caloHitSource_;
 }

◆ ~HGcalHitCheck()

HGcalHitCheck::~HGcalHitCheck ( )

inlineoverride

Definition at line 34 of file HGCalHitCheck.cc.

34 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 135 of file HGCalHitCheck.cc.

References MillePedeFileConverter_cfg::e, relativeConstraints::geom, geomToken_, edm::EventSetup::getData(), h_hit1_, h_hit2_, h_hitF_, h_hitL_, h_hitP_, h_hits_, histos_, hfClusterShapes_cfi::hits, l1ctLayer2EG_cff::id, hit::id, info(), edm::HandleBase::isValid(), nameSense_, TrackingDataMCValidation_Standalone_cff::nhits, tok_calo_, verbosity_, HGCalTypes::WaferFull, HGCalWaferIndex::waferU(), and HGCalWaferIndex::waferV().

                                                                       {
   if (verbosity_ > 0)
     edm::LogVerbatim("HitStudy") << "Run = " << e.id().run() << " Event = " << e.id().event();
 
   // get hcalGeometry
   const HGCalGeometry* geom = &iS.getData(geomToken_);
   const HGCalDDDConstants& hgc = geom->topology().dddConstants();
   const std::vector<DetId>& validIds = geom->getValidDetIds();
   edm::LogVerbatim("HitStudy") << "Detector " << nameSense_ << " with " << validIds.size() << " valid cells";
 
   const edm::Handle<edm::PCaloHitContainer>& hitsCalo = e.getHandle(tok_calo_);
   bool getHits = (hitsCalo.isValid());
   uint32_t nhits = (getHits) ? hitsCalo->size() : 0;
   uint32_t wafer(0), tiles(0);
   if (verbosity_ > 1)
     edm::LogVerbatim("HitStudy") << "HGcalHitCheck: Input flags Hits " << getHits << " with " << nhits << " hits";
   if (histos_)
     h_hits_->Fill(nhits);
 
   if (getHits) {
     std::vector<PCaloHit> hits;
     hits.insert(hits.end(), hitsCalo->begin(), hitsCalo->end());
     if (!hits.empty()) {
       for (auto hit : hits) {
         if (histos_) {
           if ((nameSense_ == "HGCalEESensitive") || (nameSense_ == "HGCalHESiliconSensitive")) {
             ++wafer;
             HGCSiliconDetId id(hit.id());
             int lay = id.layer();
             h_hitL_[lay - 1]->Fill(nhits);
             HGCalParameters::waferInfo info = hgc.waferInfo(lay, id.waferU(), id.waferV());
             if (info.part == HGCalTypes::WaferFull) {
               h_hit1_->Fill(nhits);
               h_hitF_[lay - 1]->Fill(nhits);
             } else {
               h_hit2_->Fill(nhits);
               h_hitP_[lay - 1]->Fill(nhits);
             }
           } else {
             ++tiles;
             HGCScintillatorDetId id(hit.id());
             int lay = id.layer();
             h_hitL_[lay - 1]->Fill(nhits);
             int sipm = id.sipm();
             if (sipm == 1) {
               h_hit2_->Fill(nhits);
               h_hitP_[lay - 1]->Fill(nhits);
             } else {
               h_hit1_->Fill(nhits);
               h_hitF_[lay - 1]->Fill(nhits);
             }
           }
         }
       }
     }
   }
   edm::LogVerbatim("HitStudy") << "Total hits = " << nhits << " Wafer DetIds = " << wafer << " Tile DetIds = " << tiles;
 }

◆ beginJob()

void HGcalHitCheck::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 79 of file HGCalHitCheck.cc.

References compareTotals::fs, h_hit1_, h_hit2_, h_hitF_, h_hitL_, h_hitP_, h_hits_, histos_, dqmdumpme::k, layers_, Skims_PA_cff::name, nameSense_, tag_, and runGCPTkAlMap::title.

                              {
   edm::Service<TFileService> fs;
   if (!fs.isAvailable()) {
     edm::LogVerbatim("HitStudy") << "TFileService unavailable: no histograms";
   } else {
     histos_ = true;
     char name[100], title[200];
     sprintf(name, "HitsL");
     sprintf(title, "Number of hits in %s for %s", nameSense_.c_str(), tag_.c_str());
     h_hits_ = fs->make<TH1D>(name, title, 1000, 0, 5000.);
     h_hits_->GetXaxis()->SetTitle(title);
     h_hits_->GetYaxis()->SetTitle("Hits");
     h_hits_->Sumw2();
     sprintf(name, "HitsF");
     sprintf(title, "Number of hits in %s for %s in Full Wafers or SiPM 2", nameSense_.c_str(), tag_.c_str());
     h_hit1_ = fs->make<TH1D>(name, title, 1000, 0, 5000.);
     h_hit1_->GetXaxis()->SetTitle(title);
     h_hit1_->GetYaxis()->SetTitle("Hits");
     h_hit1_->Sumw2();
     sprintf(name, "HitsP");
     sprintf(title, "Number of hits in %s for %s in Partial Wafers or SiPM 4", nameSense_.c_str(), tag_.c_str());
     h_hit2_ = fs->make<TH1D>(name, title, 1000, 0, 5000.);
     h_hit2_->GetXaxis()->SetTitle(title);
     h_hit2_->GetYaxis()->SetTitle("Hits");
     h_hit2_->Sumw2();
     for (int k = 0; k < layers_; ++k) {
       sprintf(name, "HitsL%d", k + 1);
       sprintf(title, "Number of hits in %s for %s in Layer %d", nameSense_.c_str(), tag_.c_str(), k + 1);
       h_hitL_.emplace_back(fs->make<TH1D>(name, title, 1000, 0, 5000.));
       h_hitL_.back()->GetXaxis()->SetTitle(title);
       h_hitL_.back()->GetYaxis()->SetTitle("Hits");
       h_hitL_.back()->Sumw2();
       sprintf(name, "HitsF%d", k + 1);
       sprintf(title,
               "Number of hits in %s for %s in Full Wafers or SiPM 2 of Layer %d",
               nameSense_.c_str(),
               tag_.c_str(),
               k + 1);
       h_hitF_.emplace_back(fs->make<TH1D>(name, title, 1000, 0, 5000.));
       h_hitF_.back()->GetXaxis()->SetTitle(title);
       h_hitF_.back()->GetYaxis()->SetTitle("Hits");
       h_hitF_.back()->Sumw2();
       sprintf(name, "HitsP%d", k + 1);
       sprintf(title,
               "Number of hits in %s for %s in Partial Wafers or SiPM 4  of Layer %d",
               nameSense_.c_str(),
               tag_.c_str(),
               k + 1);
       h_hitP_.emplace_back(fs->make<TH1D>(name, title, 1000, 0, 5000.));
       h_hitP_.back()->GetXaxis()->SetTitle(title);
       h_hitP_.back()->GetYaxis()->SetTitle("Hits");
       h_hitP_.back()->Sumw2();
     }
   }
 }

◆ endJob()

void HGcalHitCheck::endJob ( void )

inlineoverrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 40 of file HGCalHitCheck.cc.

40 {}

◆ fillDescriptions()

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

static

Definition at line 67 of file HGCalHitCheck.cc.

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

                                                                                {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("moduleLabel", "g4SimHits");
   desc.add<std::string>("caloHitSource", "HGCHitsEE");
   desc.add<std::string>("nameSense", "HGCalEESensitive");
   desc.add<std::string>("nameDevice", "HGCal EE");
   desc.add<std::string>("tag", "DDD");
   desc.add<int>("layers", 26);
   desc.add<int>("verbosity", 0);
   descriptions.add("hgcalHitCheckEE", desc);
 }

Member Data Documentation

◆ caloHitSource_

const std::string HGcalHitCheck::caloHitSource_

private

Definition at line 43 of file HGCalHitCheck.cc.

◆ g4Label_

const std::string HGcalHitCheck::g4Label_

private

Definition at line 43 of file HGCalHitCheck.cc.

◆ geomToken_

const edm::ESGetToken<HGCalGeometry, IdealGeometryRecord> HGcalHitCheck::geomToken_

private

Definition at line 46 of file HGCalHitCheck.cc.

Referenced by analyze().

◆ h_hit1_

TH1D * HGcalHitCheck::h_hit1_

private

Definition at line 48 of file HGCalHitCheck.cc.

Referenced by analyze(), and beginJob().

◆ h_hit2_

TH1D * HGcalHitCheck::h_hit2_

private

Definition at line 48 of file HGCalHitCheck.cc.

Referenced by analyze(), and beginJob().

◆ h_hitF_

std::vector<TH1D*> HGcalHitCheck::h_hitF_

private

Definition at line 49 of file HGCalHitCheck.cc.

Referenced by analyze(), and beginJob().

◆ h_hitL_

std::vector<TH1D*> HGcalHitCheck::h_hitL_

private

Definition at line 49 of file HGCalHitCheck.cc.

Referenced by analyze(), and beginJob().

◆ h_hitP_

std::vector<TH1D*> HGcalHitCheck::h_hitP_

private

Definition at line 49 of file HGCalHitCheck.cc.

Referenced by analyze(), and beginJob().

◆ h_hits_

TH1D* HGcalHitCheck::h_hits_

private

Definition at line 48 of file HGCalHitCheck.cc.

Referenced by analyze(), and beginJob().

◆ histos_

bool HGcalHitCheck::histos_

private

Definition at line 47 of file HGCalHitCheck.cc.

Referenced by analyze(), and beginJob().

◆ layers_

const int HGcalHitCheck::layers_

private

Definition at line 44 of file HGCalHitCheck.cc.

Referenced by beginJob().

◆ nameDetector_

const std::string HGcalHitCheck::nameDetector_

private

Definition at line 43 of file HGCalHitCheck.cc.

◆ nameSense_

const std::string HGcalHitCheck::nameSense_

private

Definition at line 43 of file HGCalHitCheck.cc.

Referenced by analyze(), beginJob(), and HGcalHitCheck().

◆ tag_

const std::string HGcalHitCheck::tag_

private

Definition at line 43 of file HGCalHitCheck.cc.

Referenced by beginJob().

◆ tok_calo_

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

private

Definition at line 45 of file HGCalHitCheck.cc.

Referenced by analyze().

◆ verbosity_

const int HGcalHitCheck::verbosity_

private

Definition at line 44 of file HGCalHitCheck.cc.

Referenced by analyze().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGcalHitCheck()

◆ ~HGcalHitCheck()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ fillDescriptions()

Member Data Documentation

◆ caloHitSource_

◆ g4Label_

◆ geomToken_

◆ h_hit1_

◆ h_hit2_

◆ h_hitF_

◆ h_hitL_

◆ h_hitP_

◆ h_hits_

◆ histos_

◆ layers_

◆ nameDetector_

◆ nameSense_

◆ tag_

◆ tok_calo_

◆ verbosity_