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 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 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_, 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_