Inheritance diagram for HcalGeomCheck:

Classes
struct	hitsinfo

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

	~HcalGeomCheck () 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)

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

Private Member Functions
void	analyzeHits (int, const std::string &, const std::vector< PCaloHit > &)

Private Attributes
const std::string	caloHitSource_

const int	depthMax_

std::vector< TH1D * >	h_E_

std::map< int, TH1D * >	h_eta_

std::map< std::pair< int, int >, TH1D * >	h_phi_

TH2D *	h_RZ_

std::vector< TH1D * >	h_T_

const HcalDDDRecConstants *	hcons_

const int	ietaMax_

const int	ietaMin_

const int	nbinR_

const int	nbinZ_

const double	rmax_

const double	rmin_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hits_

edm::ESGetToken< HcalDDDRecConstants, HcalRecNumberingRecord >	tok_HRNDC_

const int	verbosity_

const double	zmax_

const double	zmin_

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

Constructor & Destructor Documentation

◆ HcalGeomCheck()

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

explicit

Definition at line 76 of file HcalGeomCheck.cc.

     :
  
       caloHitSource_(iConfig.getParameter<std::string>("caloHitSource")),
       ietaMin_(iConfig.getUntrackedParameter<int>("ietaMin", -41)),
       ietaMax_(iConfig.getUntrackedParameter<int>("ietaMax", 41)),
       depthMax_(iConfig.getUntrackedParameter<int>("depthMax", 7)),
       rmin_(iConfig.getUntrackedParameter<double>("rMin", 0.0)),
       rmax_(iConfig.getUntrackedParameter<double>("rMax", 5500.0)),
       zmin_(iConfig.getUntrackedParameter<double>("zMin", -12500.0)),
       zmax_(iConfig.getUntrackedParameter<double>("zMax", 12500.0)),
       nbinR_(iConfig.getUntrackedParameter<int>("nBinR", 550)),
       nbinZ_(iConfig.getUntrackedParameter<int>("nBinZ", 2500)),
       verbosity_(iConfig.getUntrackedParameter<int>("verbosity", 0)) {
   usesResource(TFileService::kSharedResource);
   tok_hits_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", caloHitSource_));
   tok_HRNDC_ = esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>();
 }

References caloHitSource_, TFileService::kSharedResource, tok_hits_, and tok_HRNDC_.

◆ ~HcalGeomCheck()

HcalGeomCheck::~HcalGeomCheck ( )

inlineoverride

Definition at line 48 of file HcalGeomCheck.cc.

48 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 111 of file HcalGeomCheck.cc.

                                                                                {
   //Now the hits
   edm::Handle<edm::PCaloHitContainer> theCaloHitContainer;
   iEvent.getByToken(tok_hits_, theCaloHitContainer);
   if (theCaloHitContainer.isValid()) {
     if (verbosity_ > 0)
       edm::LogVerbatim("HcalValidation") << " PcalohitItr = " << theCaloHitContainer->size();
  
     //Merge hits for the same DetID
     std::map<HcalDetId, hitsinfo> map_hits;
     unsigned int nused(0);
     for (auto const& hit : *(theCaloHitContainer.product())) {
       unsigned int id = hit.id();
       HcalDetId detId = HcalHitRelabeller::relabel(id, hcons_);
       double energy = hit.energy();
       double time = hit.time();
       int subdet = detId.subdet();
       int ieta = detId.ieta();
       int iphi = detId.iphi();
       int depth = detId.depth();
       std::pair<double, double> etaphi = hcons_->getEtaPhi(subdet, ieta, iphi);
       double rz = hcons_->getRZ(subdet, ieta, depth);
       if (verbosity_ > 2)
         edm::LogVerbatim("HcalValidation") << "i/p " << subdet << ":" << ieta << ":" << iphi << ":" << depth << " o/p "
                                            << etaphi.first << ":" << etaphi.second << ":" << rz;
       HepGeom::Point3D<float> gcoord = HepGeom::Point3D<float>(rz * cos(etaphi.second) / cosh(etaphi.first),
                                                                rz * sin(etaphi.second) / cosh(etaphi.first),
                                                                rz * tanh(etaphi.first));
       nused++;
       double tof = (gcoord.mag() * CLHEP::mm) / CLHEP::c_light;
       if (verbosity_ > 1)
         edm::LogVerbatim("HcalValidation")
             << "Detector " << subdet << " ieta = " << ieta << " iphi = " << iphi << " depth = " << depth
             << " positon = " << gcoord << " energy = " << energy << " time = " << time << ":" << tof;
       time -= tof;
       if (time < 0)
         time = 0;
       hitsinfo hinfo;
       if (map_hits.count(detId) != 0) {
         hinfo = map_hits[detId];
       } else {
         hinfo.phi = gcoord.getPhi();
         hinfo.eta = gcoord.getEta();
         hinfo.time = time;
       }
       hinfo.energy += energy;
       map_hits[detId] = hinfo;
  
       h_RZ_->Fill(gcoord.z(), gcoord.rho());
     }
  
     //Fill in histograms
     for (auto const& hit : map_hits) {
       hitsinfo hinfo = hit.second;
       int subdet = hit.first.subdet();
       if (subdet > 0 && subdet < static_cast<int>(h_E_.size())) {
         int depth = hit.first.depth();
         int ieta = hit.first.ieta();
         int iphi = hit.first.iphi();
         if (verbosity_ > 1)
           edm::LogVerbatim("HGCalValidation")
               << " ----------------------   eta = " << ieta << ":" << hinfo.eta << " phi = " << iphi << ":" << hinfo.phi
               << " depth = " << depth << " E = " << hinfo.energy << " T = " << hinfo.time;
         h_E_[subdet]->Fill(hinfo.energy);
         h_T_[subdet]->Fill(hinfo.time);
         auto itr1 = h_phi_.find(std::pair<int, int>(ieta, depth));
         if (itr1 != h_phi_.end())
           itr1->second->Fill(iphi);
         auto itr2 = h_eta_.find(depth);
         if (itr2 != h_eta_.end())
           itr2->second->Fill(ieta);
       }
     }
   } else if (verbosity_ > 0) {
     edm::LogVerbatim("HcalValidation") << "PCaloHitContainer does not "
                                        << "exist for HCAL";
   }
 }

References funct::cos(), LEDCalibrationChannels::depth, HcalDetId::depth(), HCALHighEnergyHPDFilter_cfi::energy, fastmath::etaphi(), HcalDDDRecConstants::getEtaPhi(), HcalDDDRecConstants::getRZ(), h_E_, h_eta_, h_phi_, h_RZ_, h_T_, hcons_, hit::id, LEDCalibrationChannels::ieta, HcalDetId::ieta(), iEvent, LEDCalibrationChannels::iphi, HcalDetId::iphi(), edm::HandleBase::isValid(), edm::Handle< T >::product(), HcalHitRelabeller::relabel(), funct::sin(), HcalDetId::subdet(), protons_cff::time, tok_hits_, and verbosity_.

◆ analyzeHits()

void HcalGeomCheck::analyzeHits	(	int	,
		const std::string &	,
		const std::vector< PCaloHit > &
	)

private

◆ beginJob()

void HcalGeomCheck::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 198 of file HcalGeomCheck.cc.

                              {
   if (verbosity_ > 2)
     edm::LogVerbatim("HcalValidation") << "HcalGeomCheck:: Enter beginJob";
   edm::Service<TFileService> fs;
   h_RZ_ = fs->make<TH2D>("RZ", "R vs Z", nbinZ_, zmin_, zmax_, nbinR_, rmin_, rmax_);
   if (verbosity_ > 2)
     edm::LogVerbatim("HcalValidation") << "HcalGeomCheck: booked scatterplot RZ";
   char name[20], title[100];
   for (int depth = 1; depth < depthMax_; ++depth) {
     for (int ieta = ietaMin_; ieta <= ietaMax_; ++ieta) {
       sprintf(name, "phi%d%d", ieta, depth);
       sprintf(title, "i#phi (i#eta = %d, depth = %d)", ieta, depth);
       h_phi_[std::pair<int, int>(ieta, depth)] = fs->make<TH1D>(name, title, 400, -20, 380);
       if (verbosity_ > 2)
         edm::LogVerbatim("HcalValidation") << "HcalGeomCheck: books " << title;
     }
     sprintf(name, "eta%d", depth);
     sprintf(title, "i#eta (depth = %d)", depth);
     h_eta_[depth] = fs->make<TH1D>(name, title, 100, -50, 50);
     if (verbosity_ > 2)
       edm::LogVerbatim("HcalValidation") << "HcalGeomCheck: books " << title;
   }
  
   std::vector<std::string> dets = {"HCAL", "HB", "HE", "HF", "HO"};
   for (unsigned int ih = 0; ih < dets.size(); ++ih) {
     sprintf(name, "E_%s", dets[ih].c_str());
     sprintf(title, "Energy deposit in %s (MeV)", dets[ih].c_str());
     h_E_.emplace_back(fs->make<TH1D>(name, title, 1000, 0.0, 1.0));
     if (verbosity_ > 2)
       edm::LogVerbatim("HcalValidation") << "HcalGeomCheck: books " << title;
     sprintf(name, "T_%s", dets[ih].c_str());
     sprintf(title, "Time of hit in %s (ns)", dets[ih].c_str());
     h_T_.emplace_back(fs->make<TH1D>(name, title, 1000, 0.0, 200.0));
     if (verbosity_ > 2)
       edm::LogVerbatim("HcalValidation") << "HcalGeomCheck: books " << title;
   }
 }

References LEDCalibrationChannels::depth, depthMax_, h_E_, h_eta_, h_phi_, h_RZ_, h_T_, LEDCalibrationChannels::ieta, ietaMax_, ietaMin_, TFileService::make(), Skims_PA_cff::name, nbinR_, nbinZ_, rmax_, rmin_, runGCPTkAlMap::title, verbosity_, zmax_, and zmin_.

◆ beginRun()

void HcalGeomCheck::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	iSetup
	)

overrideprotected

Definition at line 191 of file HcalGeomCheck.cc.

                                                                        {
   const auto& pHRNDC = iSetup.getData(tok_HRNDC_);
   hcons_ = &pHRNDC;
   if (verbosity_ > 0)
     edm::LogVerbatim("HcalValidation") << "Obtain HcalDDDRecConstants from Event Setup";
 }

References edm::EventSetup::getData(), hcons_, tok_HRNDC_, and verbosity_.

◆ endRun()

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

inlineoverrideprotected

Definition at line 55 of file HcalGeomCheck.cc.

55 {}

◆ fillDescriptions()

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

static

Definition at line 95 of file HcalGeomCheck.cc.

                                                                                {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("caloHitSource", "HcalHits");
   desc.addUntracked<int>("ietaMin", -41);
   desc.addUntracked<int>("ietaMax", 41);
   desc.addUntracked<int>("depthMax", 7);
   desc.addUntracked<double>("rMin", 0.0);
   desc.addUntracked<double>("rMax", 5500.0);
   desc.addUntracked<double>("zMin", -12500.0);
   desc.addUntracked<double>("zMax", 12500.0);
   desc.addUntracked<int>("nBinR", 550);
   desc.addUntracked<int>("nBinZ", 250);
   desc.addUntracked<int>("verbosity", 0);
   descriptions.add("hcalGeomCheck", desc);
 }

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

Member Data Documentation

◆ caloHitSource_

const std::string HcalGeomCheck::caloHitSource_

private

Definition at line 61 of file HcalGeomCheck.cc.

Referenced by HcalGeomCheck().

◆ depthMax_

const int HcalGeomCheck::depthMax_

private

Definition at line 62 of file HcalGeomCheck.cc.

Referenced by beginJob().

◆ h_E_

std::vector<TH1D*> HcalGeomCheck::h_E_

private

Definition at line 73 of file HcalGeomCheck.cc.

Referenced by analyze(), and beginJob().

◆ h_eta_

std::map<int, TH1D*> HcalGeomCheck::h_eta_

private

Definition at line 72 of file HcalGeomCheck.cc.

Referenced by analyze(), and beginJob().

◆ h_phi_

std::map<std::pair<int, int>, TH1D*> HcalGeomCheck::h_phi_

private

Definition at line 71 of file HcalGeomCheck.cc.

Referenced by analyze(), and beginJob().

◆ h_RZ_

TH2D* HcalGeomCheck::h_RZ_

private

Definition at line 70 of file HcalGeomCheck.cc.

Referenced by analyze(), and beginJob().

◆ h_T_

std::vector<TH1D*> HcalGeomCheck::h_T_

private

Definition at line 73 of file HcalGeomCheck.cc.

Referenced by analyze(), and beginJob().

◆ hcons_

const HcalDDDRecConstants* HcalGeomCheck::hcons_

private

Definition at line 65 of file HcalGeomCheck.cc.

Referenced by analyze(), and beginRun().

◆ ietaMax_

const int HcalGeomCheck::ietaMax_

private

Definition at line 62 of file HcalGeomCheck.cc.

Referenced by beginJob().

◆ ietaMin_

const int HcalGeomCheck::ietaMin_

private

Definition at line 62 of file HcalGeomCheck.cc.

Referenced by beginJob().

◆ nbinR_

const int HcalGeomCheck::nbinR_

private

Definition at line 64 of file HcalGeomCheck.cc.

Referenced by beginJob().

◆ nbinZ_

const int HcalGeomCheck::nbinZ_

private

Definition at line 64 of file HcalGeomCheck.cc.

Referenced by beginJob().

◆ rmax_

const double HcalGeomCheck::rmax_

private

Definition at line 63 of file HcalGeomCheck.cc.

Referenced by beginJob().

◆ rmin_

const double HcalGeomCheck::rmin_

private

Definition at line 63 of file HcalGeomCheck.cc.

Referenced by beginJob().

◆ tok_hits_

edm::EDGetTokenT<edm::PCaloHitContainer> HcalGeomCheck::tok_hits_

private

Definition at line 66 of file HcalGeomCheck.cc.

Referenced by analyze(), and HcalGeomCheck().

◆ tok_HRNDC_

edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> HcalGeomCheck::tok_HRNDC_

private

Definition at line 67 of file HcalGeomCheck.cc.

Referenced by beginRun(), and HcalGeomCheck().

◆ verbosity_

const int HcalGeomCheck::verbosity_

private

Definition at line 64 of file HcalGeomCheck.cc.

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

◆ zmax_

const double HcalGeomCheck::zmax_

private

Definition at line 63 of file HcalGeomCheck.cc.

Referenced by beginJob().

◆ zmin_

const double HcalGeomCheck::zmin_

private

Definition at line 63 of file HcalGeomCheck.cc.

Referenced by beginJob().

Classes

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HcalGeomCheck()

◆ ~HcalGeomCheck()

Member Function Documentation

◆ analyze()

◆ analyzeHits()

◆ beginJob()

◆ beginRun()

◆ endRun()

◆ fillDescriptions()

Member Data Documentation

◆ caloHitSource_

◆ depthMax_

◆ h_E_

◆ h_eta_

◆ h_phi_

◆ h_RZ_

◆ h_T_

◆ hcons_

◆ ietaMax_

◆ ietaMin_

◆ nbinR_

◆ nbinZ_

◆ rmax_

◆ rmin_

◆ tok_hits_

◆ tok_HRNDC_

◆ verbosity_

◆ zmax_

◆ zmin_