Inheritance diagram for HGCalSimHitStudy:

Classes
struct	hitsinfo

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

	~HGCalSimHitStudy () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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::vector< ModuleDescription const * > &modules, 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

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::vector< std::string >	caloHitSources_

const double	etamax_

const double	etamin_

std::vector< TH1D * >	h_C1_

std::vector< TH1D * >	h_C2_

std::vector< TH1D * >	h_E_

std::vector< TH2D * >	h_EtaPhi_

std::vector< TH2D * >	h_EtFiZm_

std::vector< TH2D * >	h_EtFiZp_

std::vector< TH1D * >	h_LayerZm_

std::vector< TH1D * >	h_LayerZp_

std::vector< TH1D * >	h_Ly_

std::vector< TH2D * >	h_RZ_

std::vector< TH1D * >	h_T_

std::vector< TH1D * >	h_W1_

std::vector< TH1D * >	h_W2_

std::vector< TH2D * >	h_XY_

const HcalDDDRecConstants *	hcons_

std::vector< bool >	heRebuild_

std::vector< const HGCalDDDConstants * >	hgcons_

const bool	ifLayer_

const bool	ifNose_

std::vector< int >	layerFront_

std::vector< int >	layers_

const std::vector< std::string >	nameDetectors_

const int	nbinEta_

const int	nbinR_

const int	nbinZ_

const int	nLayers_

const double	rmax_

const double	rmin_

std::vector< edm::ESGetToken< HGCalDDDConstants, IdealGeometryRecord > >	tok_hgcGeom_

std::vector< 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 48 of file HGCalSimHitStudy.cc.

Constructor & Destructor Documentation

◆ HGCalSimHitStudy()

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

explicit

Definition at line 92 of file HGCalSimHitStudy.cc.

     : nameDetectors_(iConfig.getParameter<std::vector<std::string> >("detectorNames")),
       caloHitSources_(iConfig.getParameter<std::vector<std::string> >("caloHitSources")),
       rmin_(iConfig.getUntrackedParameter<double>("rMin", 0.0)),
       rmax_(iConfig.getUntrackedParameter<double>("rMax", 3000.0)),
       zmin_(iConfig.getUntrackedParameter<double>("zMin", 3000.0)),
       zmax_(iConfig.getUntrackedParameter<double>("zMax", 6000.0)),
       etamin_(iConfig.getUntrackedParameter<double>("etaMin", 1.0)),
       etamax_(iConfig.getUntrackedParameter<double>("etaMax", 3.0)),
       nbinR_(iConfig.getUntrackedParameter<int>("nBinR", 300)),
       nbinZ_(iConfig.getUntrackedParameter<int>("nBinZ", 300)),
       nbinEta_(iConfig.getUntrackedParameter<int>("nBinEta", 200)),
       nLayers_(iConfig.getUntrackedParameter<int>("layers", 50)),
       verbosity_(iConfig.getUntrackedParameter<int>("verbosity", 0)),
       ifNose_(iConfig.getUntrackedParameter<bool>("ifNose", false)),
       ifLayer_(iConfig.getUntrackedParameter<bool>("ifLayer", false)) {
   usesResource(TFileService::kSharedResource);
  
   for (auto const& name : nameDetectors_) {
     if (name == "HCal") {
       heRebuild_.emplace_back(true);
       tok_hrndc_ = esConsumes<HcalDDDRecConstants, HcalRecNumberingRecord, edm::Transition::BeginRun>();
       tok_hgcGeom_.emplace_back(esConsumes<HGCalDDDConstants, IdealGeometryRecord, edm::Transition::BeginRun>(
           edm::ESInputTag{"", "HGCalHEScintillatorSensitive"}));
     } else {
       heRebuild_.emplace_back(false);
       tok_hgcGeom_.emplace_back(
           esConsumes<HGCalDDDConstants, IdealGeometryRecord, edm::Transition::BeginRun>(edm::ESInputTag{"", name}));
     }
   }
   for (auto const& source : caloHitSources_) {
     tok_hits_.emplace_back(consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", source)));
   }
 }

References caloHitSources_, heRebuild_, TFileService::kSharedResource, Skims_PA_cff::name, nameDetectors_, source, tok_hgcGeom_, tok_hits_, and tok_hrndc_.

◆ ~HGCalSimHitStudy()

HGCalSimHitStudy::~HGCalSimHitStudy ( )

inlineoverride

Definition at line 60 of file HGCalSimHitStudy.cc.

60 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 149 of file HGCalSimHitStudy.cc.

                                                                                   {
   //Now the hits
   for (unsigned int k = 0; k < tok_hits_.size(); ++k) {
     edm::Handle<edm::PCaloHitContainer> theCaloHitContainers;
     iEvent.getByToken(tok_hits_[k], theCaloHitContainers);
     if (theCaloHitContainers.isValid()) {
       if (verbosity_ > 0)
         edm::LogVerbatim("HGCalValidation") << " PcalohitItr = " << theCaloHitContainers->size();
       std::vector<PCaloHit> caloHits;
       if (heRebuild_[k]) {
         for (auto const& hit : *(theCaloHitContainers.product())) {
           unsigned int id = hit.id();
           HcalDetId hid = HcalHitRelabeller::relabel(id, hcons_);
           if (hid.subdet() != static_cast<int>(HcalEndcap)) {
             caloHits.emplace_back(hit);
             caloHits.back().setID(hid.rawId());
             if (verbosity_ > 0)
               edm::LogVerbatim("HGCalValidation") << "Hit[" << caloHits.size() << "] " << hid;
           }
         }
       } else {
         caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), theCaloHitContainers->end());
       }
       analyzeHits(k, nameDetectors_[k], caloHits);
     } else if (verbosity_ > 0) {
       edm::LogVerbatim("HGCalValidation") << "PCaloHitContainer does not "
                                           << "exist for " << nameDetectors_[k];
     }
   }
 }

References analyzeHits(), HcalEndcap, hcons_, heRebuild_, hit::id, iEvent, edm::HandleBase::isValid(), dqmdumpme::k, nameDetectors_, edm::Handle< T >::product(), DetId::rawId(), HcalHitRelabeller::relabel(), HcalDetId::subdet(), tok_hits_, and verbosity_.

◆ analyzeHits()

void HGCalSimHitStudy::analyzeHits	(	int	ih,
		const std::string &	name,
		const std::vector< PCaloHit > &	hits
	)

private

Definition at line 180 of file HGCalSimHitStudy.cc.

                                                                                                  {
   if (verbosity_ > 0)
     edm::LogVerbatim("HGCalValidation") << name << " with " << hits.size() << " PcaloHit elements";
  
   std::map<uint32_t, hitsinfo> map_hits;
   map_hits.clear();
  
   unsigned int nused(0);
   for (auto const& hit : hits) {
     double energy = hit.energy();
     double time = hit.time();
     uint32_t id = hit.id();
     int cell, sector, sector2(0), layer, zside;
     int subdet(0), cell2(0), type(0);
     HepGeom::Point3D<float> gcoord;
     if (heRebuild_[ih]) {
       HcalDetId detId = HcalDetId(id);
       subdet = detId.subdet();
       cell = detId.ietaAbs();
       sector = detId.iphi();
       layer = detId.depth();
       zside = detId.zside();
       std::pair<double, double> etaphi = hcons_->getEtaPhi(subdet, zside * cell, sector);
       double rz = hcons_->getRZ(subdet, zside * cell, layer);
       if (verbosity_ > 2)
         edm::LogVerbatim("HGCalValidation") << "i/p " << subdet << ":" << zside << ":" << cell << ":" << sector << ":"
                                             << layer << " o/p " << etaphi.first << ":" << etaphi.second << ":" << rz;
       gcoord = HepGeom::Point3D<float>(rz * cos(etaphi.second) / cosh(etaphi.first),
                                        rz * sin(etaphi.second) / cosh(etaphi.first),
                                        rz * tanh(etaphi.first));
     } else {
       std::pair<float, float> xy;
       if (ifNose_) {
         HFNoseDetId detId = HFNoseDetId(id);
         subdet = detId.subdetId();
         cell = detId.cellU();
         cell2 = detId.cellV();
         sector = detId.waferU();
         sector2 = detId.waferV();
         type = detId.type();
         layer = detId.layer();
         zside = detId.zside();
         xy = hgcons_[ih]->locateCell(layer, sector, sector2, cell, cell2, false, true);
         h_W2_[ih]->Fill(sector2);
         h_C2_[ih]->Fill(cell2);
  
       } else if (hgcons_[ih]->waferHexagon8()) {
         HGCSiliconDetId detId = HGCSiliconDetId(id);
         subdet = static_cast<int>(detId.det());
         cell = detId.cellU();
         cell2 = detId.cellV();
         sector = detId.waferU();
         sector2 = detId.waferV();
         type = detId.type();
         layer = detId.layer();
         zside = detId.zside();
         xy = hgcons_[ih]->locateCell(layer, sector, sector2, cell, cell2, false, true);
         h_W2_[ih]->Fill(sector2);
         h_C2_[ih]->Fill(cell2);
       } else if (hgcons_[ih]->tileTrapezoid()) {
         HGCScintillatorDetId detId = HGCScintillatorDetId(id);
         subdet = static_cast<int>(detId.det());
         sector = detId.ieta();
         cell = detId.iphi();
         type = detId.type();
         layer = detId.layer();
         zside = detId.zside();
         xy = hgcons_[ih]->locateCellTrap(layer, sector, cell, false);
       } else {
         HGCalTestNumbering::unpackHexagonIndex(id, subdet, zside, layer, sector, type, cell);
         xy = hgcons_[ih]->locateCell(cell, layer, sector, false);
       }
       double zp = hgcons_[ih]->waferZ(layer, false);
       if (zside < 0)
         zp = -zp;
       double xp = (zp < 0) ? -xy.first : xy.first;
       gcoord = HepGeom::Point3D<float>(xp, xy.second, zp);
       if (verbosity_ > 2)
         edm::LogVerbatim("HGCalValidation")
             << "i/p " << subdet << ":" << zside << ":" << layer << ":" << sector << ":" << sector2 << ":" << cell << ":"
             << cell2 << " o/p " << xy.first << ":" << xy.second << ":" << zp;
     }
     nused++;
     double tof = (gcoord.mag() * CLHEP::mm) / CLHEP::c_light;
     if (verbosity_ > 1)
       edm::LogVerbatim("HGCalValidation")
           << "Detector " << name << " zside = " << zside << " layer = " << layer << " type = " << type
           << " wafer = " << sector << ":" << sector2 << " cell = " << cell << ":" << cell2 << " positon = " << gcoord
           << " energy = " << energy << " time = " << time << ":" << tof;
     time -= tof;
     if (time < 0)
       time = 0;
     hitsinfo hinfo;
     if (map_hits.count(id) != 0) {
       hinfo = map_hits[id];
     } else {
       hinfo.layer = layer + layerFront_[ih];
       hinfo.phi = gcoord.getPhi();
       hinfo.eta = gcoord.getEta();
       hinfo.time = time;
     }
     hinfo.energy += energy;
     map_hits[id] = hinfo;
  
     //Fill in histograms
     h_RZ_[0]->Fill(std::abs(gcoord.z()), gcoord.rho());
     h_RZ_[ih + 1]->Fill(std::abs(gcoord.z()), gcoord.rho());
     if (ifLayer_) {
       if (hinfo.layer <= static_cast<int>(h_XY_.size()))
         h_XY_[hinfo.layer - 1]->Fill(gcoord.x(), gcoord.y());
     } else {
       h_EtaPhi_[0]->Fill(std::abs(hinfo.eta), hinfo.phi);
       h_EtaPhi_[ih + 1]->Fill(std::abs(hinfo.eta), hinfo.phi);
     }
     h_Ly_[ih]->Fill(layer);
     h_W1_[ih]->Fill(sector);
     h_C1_[ih]->Fill(cell);
   }
   if (verbosity_ > 0)
     edm::LogVerbatim("HGCalValidation") << name << " with " << map_hits.size() << ":" << nused << " detector elements"
                                         << " being hit";
  
   for (auto const& hit : map_hits) {
     hitsinfo hinfo = hit.second;
     if (verbosity_ > 1)
       edm::LogVerbatim("HGCalValidation")
           << " ----------------------   eta = " << hinfo.eta << " phi = " << hinfo.phi << " layer = " << hinfo.layer
           << " E = " << hinfo.energy << " T = " << hinfo.time;
     h_E_[0]->Fill(hinfo.energy);
     h_E_[ih + 1]->Fill(hinfo.energy);
     h_T_[0]->Fill(hinfo.time);
     h_T_[ih + 1]->Fill(hinfo.time);
     if (hinfo.eta > 0) {
       if (!ifLayer_) {
         h_EtFiZp_[0]->Fill(std::abs(hinfo.eta), hinfo.phi, hinfo.energy);
         h_EtFiZp_[ih + 1]->Fill(std::abs(hinfo.eta), hinfo.phi, hinfo.energy);
       }
       h_LayerZp_[0]->Fill(hinfo.layer, hinfo.energy);
       h_LayerZp_[ih + 1]->Fill(hinfo.layer, hinfo.energy);
     } else {
       if (!ifLayer_) {
         h_EtFiZm_[0]->Fill(std::abs(hinfo.eta), hinfo.phi, hinfo.energy);
         h_EtFiZm_[ih + 1]->Fill(std::abs(hinfo.eta), hinfo.phi, hinfo.energy);
       }
       h_LayerZm_[0]->Fill(hinfo.layer, hinfo.energy);
       h_LayerZm_[ih + 1]->Fill(hinfo.layer, hinfo.energy);
     }
   }
 }

References funct::abs(), hgcalTopologyTester_cfi::cell2, HFNoseDetId::cellU(), HGCSiliconDetId::cellU(), HFNoseDetId::cellV(), HGCSiliconDetId::cellV(), funct::cos(), HcalDetId::depth(), DetId::det(), HCALHighEnergyHPDFilter_cfi::energy, fastmath::etaphi(), HcalDDDRecConstants::getEtaPhi(), HcalDDDRecConstants::getRZ(), h_C1_, h_C2_, h_E_, h_EtaPhi_, h_EtFiZm_, h_EtFiZp_, h_LayerZm_, h_LayerZp_, h_Ly_, h_RZ_, h_T_, h_W1_, h_W2_, h_XY_, hcons_, heRebuild_, hgcons_, hfClusterShapes_cfi::hits, triggerObjects_cff::id, hit::id, HGCScintillatorDetId::ieta(), HcalDetId::ietaAbs(), ifLayer_, ifNose_, HGCScintillatorDetId::iphi(), HcalDetId::iphi(), HGCScintillatorDetId::layer(), HFNoseDetId::layer(), HGCSiliconDetId::layer(), layerFront_, Skims_PA_cff::name, hgcalTopologyTester_cfi::sector2, funct::sin(), HcalDetId::subdet(), DetId::subdetId(), ntuplemaker::time, HGCScintillatorDetId::type(), HFNoseDetId::type(), HGCSiliconDetId::type(), HGCalTestNumbering::unpackHexagonIndex(), verbosity_, HFNoseDetId::waferU(), HGCSiliconDetId::waferU(), HGCSiliconDetId::waferV(), HFNoseDetId::waferV(), geometryCSVtoXML::xy, HGCScintillatorDetId::zside(), HFNoseDetId::zside(), HGCSiliconDetId::zside(), HcalDetId::zside(), and ecaldqm::zside().

Referenced by analyze().

◆ beginJob()

void HGCalSimHitStudy::beginJob ( void )

overrideprotectedvirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 354 of file HGCalSimHitStudy.cc.

                                 {
   edm::Service<TFileService> fs;
  
   std::ostringstream name, title;
   for (unsigned int ih = 0; ih <= nameDetectors_.size(); ++ih) {
     name.str("");
     title.str("");
     if (ih == 0) {
       name << "RZ_AllDetectors";
       title << "R vs Z for All Detectors";
     } else {
       name << "RZ_" << nameDetectors_[ih - 1];
       title << "R vs Z for " << nameDetectors_[ih - 1];
     }
     h_RZ_.emplace_back(
         fs->make<TH2D>(name.str().c_str(), title.str().c_str(), nbinZ_, zmin_, zmax_, nbinR_, rmin_, rmax_));
     if (ifLayer_) {
       if (ih == 0) {
         for (int ly = 0; ly < nLayers_; ++ly) {
           name.str("");
           title.str("");
           name << "XY_L" << (ly + 1);
           title << "Y vs X at Layer " << (ly + 1);
           h_XY_.emplace_back(
               fs->make<TH2D>(name.str().c_str(), title.str().c_str(), nbinR_, -rmax_, rmax_, nbinR_, -rmax_, rmax_));
         }
       }
     } else {
       name.str("");
       title.str("");
       if (ih == 0) {
         name << "EtaPhi_AllDetectors";
         title << "#phi vs #eta for All Detectors";
       } else {
         name << "EtaPhi_" << nameDetectors_[ih - 1];
         title << "#phi vs #eta for " << nameDetectors_[ih - 1];
       }
       h_EtaPhi_.emplace_back(
           fs->make<TH2D>(name.str().c_str(), title.str().c_str(), nbinEta_, etamin_, etamax_, 200, -M_PI, M_PI));
       name.str("");
       title.str("");
       if (ih == 0) {
         name << "EtFiZp_AllDetectors";
         title << "#phi vs #eta (+z) for All Detectors";
       } else {
         name << "EtFiZp_" << nameDetectors_[ih - 1];
         title << "#phi vs #eta (+z) for " << nameDetectors_[ih - 1];
       }
       h_EtFiZp_.emplace_back(
           fs->make<TH2D>(name.str().c_str(), title.str().c_str(), nbinEta_, etamin_, etamax_, 200, -M_PI, M_PI));
       name.str("");
       title.str("");
       if (ih == 0) {
         name << "EtFiZm_AllDetectors";
         title << "#phi vs #eta (-z) for All Detectors";
       } else {
         name << "EtFiZm_" << nameDetectors_[ih - 1];
         title << "#phi vs #eta (-z) for " << nameDetectors_[ih - 1];
       }
       h_EtFiZm_.emplace_back(
           fs->make<TH2D>(name.str().c_str(), title.str().c_str(), nbinEta_, etamin_, etamax_, 200, -M_PI, M_PI));
     }
     name.str("");
     title.str("");
     if (ih == 0) {
       name << "LayerZp_AllDetectors";
       title << "Energy vs Layer (+z) for All Detectors";
     } else {
       name << "LayerZp_" << nameDetectors_[ih - 1];
       title << "Energy vs Layer (+z) for " << nameDetectors_[ih - 1];
     }
     h_LayerZp_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 60, 0.0, 60.0));
     name.str("");
     title.str("");
     if (ih == 0) {
       name << "LayerZm_AllDetectors";
       title << "Energy vs Layer (-z) for All Detectors";
     } else {
       name << "LayerZm_" << nameDetectors_[ih - 1];
       title << "Energy vs Layer (-z) for " << nameDetectors_[ih - 1];
     }
     h_LayerZm_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 60, 0.0, 60.0));
  
     name.str("");
     title.str("");
     if (ih == 0) {
       name << "E_AllDetectors";
       title << "Energy Deposit for All Detectors";
     } else {
       name << "E_" << nameDetectors_[ih - 1];
       title << "Energy Deposit for " << nameDetectors_[ih - 1];
     }
     h_E_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 1000, 0.0, 1.0));
  
     name.str("");
     title.str("");
     if (ih == 0) {
       name << "T_AllDetectors";
       title << "Time for All Detectors";
     } else {
       name << "T_" << nameDetectors_[ih - 1];
       title << "Time for " << nameDetectors_[ih - 1];
     }
     h_T_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 1000, 0.0, 200.0));
   }
  
   for (unsigned int ih = 0; ih < nameDetectors_.size(); ++ih) {
     name.str("");
     title.str("");
     name << "LY_" << nameDetectors_[ih];
     title << "Layer number for " << nameDetectors_[ih];
     h_Ly_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 200, 0, 100));
     if ((nameDetectors_[ih] == "HGCalHEScintillatorSensitive") || heRebuild_[ih]) {
       name.str("");
       title.str("");
       name << "IR_" << nameDetectors_[ih];
       title << "Radius index for " << nameDetectors_[ih];
       h_W1_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 200, -50, 50));
       name.str("");
       title.str("");
       name << "FI_" << nameDetectors_[ih];
       title << "#phi index for " << nameDetectors_[ih];
       h_C1_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 720, 0, 360));
     } else {
       name.str("");
       title.str("");
       name << "WU_" << nameDetectors_[ih];
       title << "u index of wafers for " << nameDetectors_[ih];
       h_W1_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 200, -50, 50));
       name.str("");
       title.str("");
       name << "WV_" << nameDetectors_[ih];
       title << "v index of wafers for " << nameDetectors_[ih];
       h_W2_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 100, -50, 50));
       name.str("");
       title.str("");
       name << "CU_" << nameDetectors_[ih];
       title << "u index of cells for " << nameDetectors_[ih];
       h_C1_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 100, 0, 50));
       name.str("");
       title.str("");
       name << "CV_" << nameDetectors_[ih];
       title << "v index of cells for " << nameDetectors_[ih];
       h_C2_.emplace_back(fs->make<TH1D>(name.str().c_str(), title.str().c_str(), 100, 0, 50));
     }
   }
 }

References etamax_, etamin_, h_C1_, h_C2_, h_E_, h_EtaPhi_, h_EtFiZm_, h_EtFiZp_, h_LayerZm_, h_LayerZp_, h_Ly_, h_RZ_, h_T_, h_W1_, h_W2_, h_XY_, heRebuild_, ifLayer_, M_PI, TFileService::make(), Skims_PA_cff::name, nameDetectors_, nbinEta_, nbinR_, nbinZ_, nLayers_, rmax_, rmin_, runGCPTkAlMap::title, zmax_, and zmin_.

◆ beginRun()

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

overrideprotected

Definition at line 331 of file HGCalSimHitStudy.cc.

                                                                           {
   for (unsigned int k = 0; k < nameDetectors_.size(); ++k) {
     if (heRebuild_[k]) {
       edm::ESHandle<HcalDDDRecConstants> pHRNDC = iSetup.getHandle(tok_hrndc_);
       hcons_ = pHRNDC.product();
       layers_.emplace_back(hcons_->getMaxDepth(1));
       hgcons_.emplace_back(nullptr);
       layerFront_.emplace_back(40);
     } else {
       edm::ESHandle<HGCalDDDConstants> pHGDC = iSetup.getHandle(tok_hgcGeom_[k]);
       hgcons_.emplace_back(pHGDC.product());
       layers_.emplace_back(hgcons_.back()->layers(false));
       if (k == 0)
         layerFront_.emplace_back(0);
       else
         layerFront_.emplace_back(28);
     }
     if (verbosity_ > 0)
       edm::LogVerbatim("HGCalValidation") << nameDetectors_[k] << " defined with " << layers_.back() << " Layers with "
                                           << layerFront_.back() << " in front";
   }
 }

References edm::EventSetup::getHandle(), HcalDDDRecConstants::getMaxDepth(), hcons_, heRebuild_, hgcons_, dqmdumpme::k, layerFront_, layers_, nameDetectors_, edm::ESHandle< T >::product(), tok_hgcGeom_, tok_hrndc_, and verbosity_.

◆ endRun()

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

inlineoverrideprotected

Definition at line 67 of file HGCalSimHitStudy.cc.

67 {}

◆ fillDescriptions()

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

static

Definition at line 127 of file HGCalSimHitStudy.cc.

                                                                                   {
   edm::ParameterSetDescription desc;
   std::vector<std::string> names = {"HGCalEESensitive", "HGCalHESiliconSensitive", "HGCalHEScintillatorSensitive"};
   std::vector<std::string> sources = {"HGCHitsEE", "HGCHitsHEfront", "HGCHitsHEback"};
   desc.add<std::vector<std::string> >("detectorNames", names);
   desc.add<std::vector<std::string> >("caloHitSources", sources);
   desc.addUntracked<double>("rMin", 0.0);
   desc.addUntracked<double>("rMax", 3000.0);
   desc.addUntracked<double>("zMin", 3000.0);
   desc.addUntracked<double>("zMax", 6000.0);
   desc.addUntracked<double>("etaMin", 1.0);
   desc.addUntracked<double>("etaMax", 3.0);
   desc.addUntracked<int>("nBinR", 300);
   desc.addUntracked<int>("nBinZ", 300);
   desc.addUntracked<int>("nBinEta", 200);
   desc.addUntracked<int>("layers", 50);
   desc.addUntracked<int>("verbosity", 0);
   desc.addUntracked<bool>("ifNose", false);
   desc.addUntracked<bool>("ifLayer", false);
   descriptions.add("hgcalSimHitStudy", desc);
 }

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, names, and CalibrationSummaryClient_cfi::sources.

Member Data Documentation

◆ caloHitSources_

const std::vector<std::string> HGCalSimHitStudy::caloHitSources_

private

Definition at line 73 of file HGCalSimHitStudy.cc.

Referenced by HGCalSimHitStudy().

◆ etamax_

const double HGCalSimHitStudy::etamax_

private

Definition at line 75 of file HGCalSimHitStudy.cc.

Referenced by beginJob().

◆ etamin_

const double HGCalSimHitStudy::etamin_

private

Definition at line 75 of file HGCalSimHitStudy.cc.

Referenced by beginJob().

◆ h_C1_

std::vector<TH1D*> HGCalSimHitStudy::h_C1_

private

Definition at line 89 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_C2_

std::vector<TH1D*> HGCalSimHitStudy::h_C2_

private

Definition at line 89 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_E_

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

private

Definition at line 88 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_EtaPhi_

std::vector<TH2D*> HGCalSimHitStudy::h_EtaPhi_

private

Definition at line 87 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_EtFiZm_

std::vector<TH2D*> HGCalSimHitStudy::h_EtFiZm_

private

Definition at line 87 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_EtFiZp_

std::vector<TH2D*> HGCalSimHitStudy::h_EtFiZp_

private

Definition at line 87 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_LayerZm_

std::vector<TH1D*> HGCalSimHitStudy::h_LayerZm_

private

Definition at line 88 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_LayerZp_

std::vector<TH1D*> HGCalSimHitStudy::h_LayerZp_

private

Definition at line 88 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_Ly_

std::vector<TH1D*> HGCalSimHitStudy::h_Ly_

private

Definition at line 89 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_RZ_

std::vector<TH2D*> HGCalSimHitStudy::h_RZ_

private

Definition at line 87 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_T_

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

private

Definition at line 88 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_W1_

std::vector<TH1D*> HGCalSimHitStudy::h_W1_

private

Definition at line 89 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_W2_

std::vector<TH1D*> HGCalSimHitStudy::h_W2_

private

Definition at line 89 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ h_XY_

std::vector<TH2D*> HGCalSimHitStudy::h_XY_

private

Definition at line 87 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ hcons_

const HcalDDDRecConstants* HGCalSimHitStudy::hcons_

private

Definition at line 81 of file HGCalSimHitStudy.cc.

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

◆ heRebuild_

std::vector<bool> HGCalSimHitStudy::heRebuild_

private

Definition at line 82 of file HGCalSimHitStudy.cc.

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

◆ hgcons_

std::vector<const HGCalDDDConstants*> HGCalSimHitStudy::hgcons_

private

Definition at line 80 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginRun().

◆ ifLayer_

const bool HGCalSimHitStudy::ifLayer_

private

Definition at line 77 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginJob().

◆ ifNose_

const bool HGCalSimHitStudy::ifNose_

private

Definition at line 77 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits().

◆ layerFront_

std::vector<int> HGCalSimHitStudy::layerFront_

private

Definition at line 84 of file HGCalSimHitStudy.cc.

Referenced by analyzeHits(), and beginRun().

◆ layers_

std::vector<int> HGCalSimHitStudy::layers_

private

Definition at line 84 of file HGCalSimHitStudy.cc.

Referenced by beginRun().

◆ nameDetectors_

const std::vector<std::string> HGCalSimHitStudy::nameDetectors_

private

Definition at line 73 of file HGCalSimHitStudy.cc.

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

◆ nbinEta_

const int HGCalSimHitStudy::nbinEta_

private

Definition at line 76 of file HGCalSimHitStudy.cc.

Referenced by beginJob().

◆ nbinR_

const int HGCalSimHitStudy::nbinR_

private

Definition at line 76 of file HGCalSimHitStudy.cc.

Referenced by beginJob().

◆ nbinZ_

const int HGCalSimHitStudy::nbinZ_

private

Definition at line 76 of file HGCalSimHitStudy.cc.

Referenced by beginJob().

◆ nLayers_

const int HGCalSimHitStudy::nLayers_

private

Definition at line 76 of file HGCalSimHitStudy.cc.

Referenced by beginJob().

◆ rmax_

const double HGCalSimHitStudy::rmax_

private

Definition at line 74 of file HGCalSimHitStudy.cc.

Referenced by beginJob().

◆ rmin_

const double HGCalSimHitStudy::rmin_

private

Definition at line 74 of file HGCalSimHitStudy.cc.

Referenced by beginJob().

◆ tok_hgcGeom_

std::vector<edm::ESGetToken<HGCalDDDConstants, IdealGeometryRecord> > HGCalSimHitStudy::tok_hgcGeom_

private

Definition at line 79 of file HGCalSimHitStudy.cc.

Referenced by beginRun(), and HGCalSimHitStudy().

◆ tok_hits_

std::vector<edm::EDGetTokenT<edm::PCaloHitContainer> > HGCalSimHitStudy::tok_hits_

private

Definition at line 83 of file HGCalSimHitStudy.cc.

Referenced by analyze(), and HGCalSimHitStudy().

◆ tok_hrndc_

edm::ESGetToken<HcalDDDRecConstants, HcalRecNumberingRecord> HGCalSimHitStudy::tok_hrndc_

private

Definition at line 78 of file HGCalSimHitStudy.cc.

Referenced by beginRun(), and HGCalSimHitStudy().

◆ verbosity_

const int HGCalSimHitStudy::verbosity_

private

Definition at line 76 of file HGCalSimHitStudy.cc.

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

◆ zmax_

const double HGCalSimHitStudy::zmax_

private

Definition at line 74 of file HGCalSimHitStudy.cc.

Referenced by beginJob().

◆ zmin_

const double HGCalSimHitStudy::zmin_

private

Definition at line 74 of file HGCalSimHitStudy.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

◆ HGCalSimHitStudy()

◆ ~HGCalSimHitStudy()

Member Function Documentation

◆ analyze()

◆ analyzeHits()

◆ beginJob()

◆ beginRun()

◆ endRun()

◆ fillDescriptions()

Member Data Documentation

◆ caloHitSources_

◆ etamax_

◆ etamin_

◆ h_C1_

◆ h_C2_

◆ h_E_

◆ h_EtaPhi_

◆ h_EtFiZm_

◆ h_EtFiZp_

◆ h_LayerZm_

◆ h_LayerZp_

◆ h_Ly_

◆ h_RZ_

◆ h_T_

◆ h_W1_

◆ h_W2_

◆ h_XY_

◆ hcons_

◆ heRebuild_

◆ hgcons_

◆ ifLayer_

◆ ifNose_

◆ layerFront_

◆ layers_

◆ nameDetectors_

◆ nbinEta_

◆ nbinR_

◆ nbinZ_

◆ nLayers_

◆ rmax_

◆ rmin_

◆ tok_hgcGeom_

◆ tok_hits_

◆ tok_hrndc_

◆ verbosity_

◆ zmax_

◆ zmin_