Inheritance diagram for HcalSimHitAnalysis:

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

	~HcalSimHitAnalysis () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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)

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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...

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
std::string	g4Label_

std::string	hcalHits_

TH2F *	poszn_ [ndets_]

TH2F *	poszp_ [ndets_]

bool	testNumber_

edm::EDGetTokenT < edm::PCaloHitContainer >	tok_calo_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	tok_geom_

edm::ESGetToken < HcalDDDRecConstants, HcalRecNumberingRecord >	tok_HRNDC_

bool	verbose_

Static Private Attributes
static const int	ndets_ = 4

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 40 of file HcalSimHitAnalysis.cc.

Constructor & Destructor Documentation

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

Definition at line 62 of file HcalSimHitAnalysis.cc.

References g4Label_, edm::ParameterSet::getUntrackedParameter(), hcalHits_, HLT_FULL_cff::InputTag, TFileService::kSharedResource, AlCaHLTBitMon_QueryRunRegistry::string, testNumber_, tok_calo_, tok_geom_, tok_HRNDC_, and verbose_.

                                                                 {
   usesResource(TFileService::kSharedResource);
 
   g4Label_ = ps.getUntrackedParameter<std::string>("moduleLabel", "g4SimHits");
   hcalHits_ = ps.getUntrackedParameter<std::string>("HitCollection", "HcalHits");
   verbose_ = ps.getUntrackedParameter<bool>("Verbose", false);
   testNumber_ = ps.getUntrackedParameter<bool>("TestNumber", true);
 
   tok_calo_ = consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label_, hcalHits_));
   tok_geom_ = esConsumes<CaloGeometry, CaloGeometryRecord>();
   tok_HRNDC_ = esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord>();
 
   edm::LogVerbatim("HitStudy") << "Module Label: " << g4Label_ << "   Hits: " << hcalHits_ << " testNumber "
                                << testNumber_;
 }

HcalSimHitAnalysis::~HcalSimHitAnalysis ( )

inlineoverride

Definition at line 43 of file HcalSimHitAnalysis.cc.

43 {}

Member Function Documentation

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 120 of file HcalSimHitAnalysis.cc.

References funct::abs(), edm::EventID::event(), edm::Event::getByToken(), edm::EventSetup::getData(), HcalGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, gpuClustering::id, edm::EventBase::id(), hit::id, edm::HandleBase::isValid(), nhits, PV3DBase< T, PVType, FrameType >::perp(), poszn_, poszp_, HcalHitRelabeller::relabel(), edm::EventID::run(), testNumber_, tok_calo_, tok_geom_, tok_HRNDC_, verbose_, x, PV3DBase< T, PVType, FrameType >::x(), y, PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and ecaldqm::zside().

                                                                            {
   if (verbose_)
     edm::LogVerbatim("HitStudy") << "Run = " << e.id().run() << " Event = " << e.id().event();
 
   // get hcalGeometry
   const CaloGeometry* geo = &iS.getData(tok_geom_);
   const HcalGeometry* hgeom = static_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal, HcalBarrel));
   const auto& pHRNDC = iS.getData(tok_HRNDC_);
   const HcalDDDRecConstants* hcons = &pHRNDC;
 
   edm::Handle<edm::PCaloHitContainer> hitsCalo;
   e.getByToken(tok_calo_, hitsCalo);
   bool getHits = (hitsCalo.isValid());
   uint32_t nhits = (getHits) ? hitsCalo->size() : 0;
   if (verbose_)
     edm::LogVerbatim("HitStudy") << "HcalSimHitAnalysis: Input flags Hits " << getHits << " with " << nhits << " hits";
   if (getHits) {
     std::vector<PCaloHit> hits;
     hits.insert(hits.end(), hitsCalo->begin(), hitsCalo->end());
     if (!hits.empty()) {
       std::map<HcalDetId, double> hitMap;
       for (auto hit : hits) {
         double edep = hit.energy();
         uint32_t id = hit.id();
         HcalDetId hid = (testNumber_) ? HcalHitRelabeller::relabel(id, hcons) : HcalDetId(id);
         auto it = hitMap.find(hid);
         if (it == hitMap.end()) {
           hitMap[hid] = edep;
         } else {
           (it->second) += edep;
         }
       }
 
       for (auto it : hitMap) {
         HcalDetId id(it.first);
         GlobalPoint gpos = hgeom->getPosition(id);
         HcalSubdetector subdet = (id).subdet();
         int indx =
             ((subdet == HcalBarrel)
                  ? 0
                  : ((subdet == HcalEndcap) ? 1 : ((subdet == HcalOuter) ? 2 : ((subdet == HcalForward) ? 3 : -1))));
         if (verbose_)
           edm::LogVerbatim("HitStudy") << "HcalSimHitAnalysis: " << id << ":" << it.second << " at " << gpos
                                        << " subdet " << subdet << ":" << indx;
         if (indx >= 0) {
           double x = ((indx == 0) || (indx == 2)) ? std::abs(gpos.z()) : gpos.x();
           double y = ((indx == 0) || (indx == 2)) ? (gpos.perp()) : gpos.y();
           if (id.zside() >= 0)
             poszp_[indx]->Fill(x, y);
           else
             poszn_[indx]->Fill(x, y);
         }
       }
     }
   }
 }

void HcalSimHitAnalysis::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 87 of file HcalSimHitAnalysis.cc.

References mps_fire::i, edm::Service< T >::isAvailable(), TFileService::make(), mergeVDriftHistosByStation::name, ndets_, poszn_, poszp_, AlCaHLTBitMon_QueryRunRegistry::string, compare::tfile, and runGCPTkAlMap::title.

                                   {
   edm::Service<TFileService> tfile;
   if (!tfile.isAvailable())
     throw cms::Exception("BadConfig") << "TFileService unavailable: "
                                       << "please add it to config file";
   char name[20], title[120];
   std::string dets[ndets_] = {"HB", "HE", "HO", "HF"};
   int nx[ndets_] = {100, 100, 350, 160};
   double xlo[ndets_] = {0, -300, 0, -160};
   double xhi[ndets_] = {500, 300, 3500, 160};
   int ny[ndets_] = {100, 100, 50, 160};
   double ylo[ndets_] = {170, -300, 375, -160};
   double yhi[ndets_] = {370, 300, 425, 160};
   std::string xttl[ndets_] = {"|z| (cm)", "x (cm)", "|z| (cm)", "x (cm)"};
   std::string yttl[ndets_] = {"#rho (cm)", "y (cm)", "#rho (cm)", "y (cm)"};
   for (int i = 0; i < ndets_; i++) {
     sprintf(name, "poszp%d", i);
     sprintf(title, "%s+", dets[i].c_str());
     poszp_[i] = tfile->make<TH2F>(name, title, nx[i], xlo[i], xhi[i], ny[i], ylo[i], yhi[i]);
     poszp_[i]->GetXaxis()->SetTitle(xttl[i].c_str());
     poszp_[i]->GetYaxis()->SetTitle(yttl[i].c_str());
     sprintf(title, "%s-", dets[i].c_str());
     poszp_[i]->GetYaxis()->SetTitleOffset(1.2);
     poszp_[i]->Sumw2();
     sprintf(name, "poszn%d", i);
     poszn_[i] = tfile->make<TH2F>(name, title, nx[i], xlo[i], xhi[i], ny[i], ylo[i], yhi[i]);
     poszn_[i]->GetXaxis()->SetTitle(xttl[i].c_str());
     poszn_[i]->GetYaxis()->SetTitle(yttl[i].c_str());
     poszn_[i]->GetYaxis()->SetTitleOffset(1.2);
     poszn_[i]->Sumw2();
   }
 }

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

inlineoverrideprotected

Definition at line 49 of file HcalSimHitAnalysis.cc.

49 {}

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

inlineoverrideprotected

Definition at line 50 of file HcalSimHitAnalysis.cc.

50 {}

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

static

Definition at line 78 of file HcalSimHitAnalysis.cc.

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

                                                                                     {
   edm::ParameterSetDescription desc;
   desc.addUntracked<std::string>("ModuleLabel", "g4SimHits");
   desc.addUntracked<std::string>("HitCollection", "HcalHits");
   desc.addUntracked<bool>("Verbose", false);
   desc.addUntracked<bool>("TestNumber", true);
   descriptions.add("hcalSimHitAnalysis", desc);
 }