Inheritance diagram for HGCalSimHitValidation:

Classes
struct	energysum

struct	hitsinfo

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

	~HGCalSimHitValidation () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	analyzeHits (std::vector< PCaloHit > &hits)

TH1F *	createHisto (std::string histname, const int nbins, float minIndexX, float maxIndexX, bool isLogX=true)

bool	defineGeometry (const DDCompactView *ddViewH)

bool	defineGeometry (const cms::DDCompactView *ddViewH)

void	fillHitsInfo (std::pair< hitsinfo, energysum > hit_, unsigned int itimeslice, double esum)

void	fillMuonTomoHistos (int partialType, std::pair< hitsinfo, energysum > hit_)

void	fillOccupancyMap (std::map< int, int > &OccupancyMap, int layer)

void	histoSetting (TH1F &histo, const char xTitle, const char *yTitle="", Color_t lineColor=kBlack, Color_t markerColor=kBlack, int linewidth=1)

void	histoSetting (TH2F &histo, const char xTitle, const char *yTitle="", Color_t lineColor=kBlack, Color_t markerColor=kBlack, int linewidth=1)

Private Attributes
const std::string	caloHitSource_

std::vector< MonitorElement * >	energy_ [maxTime_]

std::vector< MonitorElement * >	energyFWCK_

std::vector< MonitorElement * >	energyFWCN_

std::vector< MonitorElement * >	energyFWF_

std::vector< MonitorElement * >	energyPWCK_

std::vector< MonitorElement * >	energyPWCN_

std::vector< MonitorElement * >	energyPWF_

std::vector< MonitorElement * >	EtaPhi_Minus_

std::vector< MonitorElement * >	EtaPhi_Plus_

int	firstLayer_

const bool	fromDDD_

const HGCalDDDConstants *	hgcons_

std::vector< MonitorElement * >	HitOccupancy_Minus_

std::vector< MonitorElement * >	HitOccupancy_Plus_

std::vector< MonitorElement * >	hitXYB_

std::vector< MonitorElement * >	hitXYFWCK_

std::vector< MonitorElement * >	hitXYFWCN_

std::vector< MonitorElement * >	hitXYFWF_

unsigned int	layers_

MonitorElement *	MeanHitOccupancy_Minus_

MonitorElement *	MeanHitOccupancy_Plus_

const std::string	nameDetector_

unsigned int	nTimes_

bool	symmDet_

const std::vector< double >	times_

const edm::ESGetToken< DDCompactView, IdealGeometryRecord >	tok_cpv_

const edm::ESGetToken< cms::DDCompactView, IdealGeometryRecord >	tok_cpvc_

edm::EDGetTokenT< edm::HepMCProduct >	tok_hepMC_

const edm::ESGetToken< HGCalDDDConstants, IdealGeometryRecord >	tok_hgcal_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hits_

std::map< uint32_t, HepGeom::Transform3D >	transMap_

const int	verbosity_

Static Private Attributes
static const unsigned int	maxTime_ = 6

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 49 of file HGCalSimHitValidation.cc.

Constructor & Destructor Documentation

◆ HGCalSimHitValidation()

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

explicit

Definition at line 128 of file HGCalSimHitValidation.cc.

References nameDetector_.

     : nameDetector_(iConfig.getParameter<std::string>("DetectorName")),
       caloHitSource_(iConfig.getParameter<std::string>("CaloHitSource")),
       times_(iConfig.getParameter<std::vector<double> >("TimeSlices")),
       verbosity_(iConfig.getUntrackedParameter<int>("Verbosity", 0)),
       fromDDD_(iConfig.getUntrackedParameter<bool>("fromDDD", true)),
       tok_hgcal_(esConsumes<HGCalDDDConstants, IdealGeometryRecord, edm::Transition::BeginRun>(
           edm::ESInputTag{"", nameDetector_})),
       tok_cpv_(esConsumes<DDCompactView, IdealGeometryRecord, edm::Transition::BeginRun>()),
       tok_cpvc_(esConsumes<cms::DDCompactView, IdealGeometryRecord, edm::Transition::BeginRun>()),
       symmDet_(true),
       firstLayer_(1) {
   tok_hepMC_ = consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"));
   tok_hits_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", caloHitSource_));
   nTimes_ = (times_.size() > maxTime_) ? maxTime_ : times_.size();
 }

◆ ~HGCalSimHitValidation()

HGCalSimHitValidation::~HGCalSimHitValidation ( )

inlineoverride

Definition at line 70 of file HGCalSimHitValidation.cc.

70 {}

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 157 of file HGCalSimHitValidation.cc.

References analyzeHits(), edm::HepMCProduct::GetEvent(), iEvent, edm::HandleBase::isValid(), dqmdumpme::k, AlCaHLTBitMon_ParallelJobs::p, tok_hepMC_, tok_hits_, and verbosity_.

                                                                                        {
   //Generator input
   if (verbosity_ > 0) {
     edm::Handle<edm::HepMCProduct> evtMC;
     iEvent.getByToken(tok_hepMC_, evtMC);
     if (!evtMC.isValid()) {
       edm::LogVerbatim("HGCalValidation") << "no HepMCProduct found";
     } else {
       const HepMC::GenEvent* myGenEvent = evtMC->GetEvent();
       unsigned int k(0);
       for (HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end();
            ++p, ++k) {
         edm::LogVerbatim("HGCalValidation") << "Particle[" << k << "] with pt " << (*p)->momentum().perp() << " eta "
                                             << (*p)->momentum().eta() << " phi " << (*p)->momentum().phi();
       }
     }
   }
 
   //Now the hits
   edm::Handle<edm::PCaloHitContainer> theCaloHitContainers;
   iEvent.getByToken(tok_hits_, theCaloHitContainers);
   if (theCaloHitContainers.isValid()) {
     if (verbosity_ > 0)
       edm::LogVerbatim("HGCalValidation") << " PcalohitItr = " << theCaloHitContainers->size();
     std::vector<PCaloHit> caloHits;
     caloHits.insert(caloHits.end(), theCaloHitContainers->begin(), theCaloHitContainers->end());
     analyzeHits(caloHits);
   } else if (verbosity_ > 0) {
     edm::LogVerbatim("HGCalValidation") << "PCaloHitContainer does not exist!";
   }
 }

◆ analyzeHits()

void HGCalSimHitValidation::analyzeHits ( std::vector< PCaloHit > & hits )

private

Definition at line 189 of file HGCalSimHitValidation.cc.

References hgcalTopologyTester_cfi::cell2, HGCSiliconDetId::cellU(), HGCSiliconDetId::cellV(), TauDecayModes::dec, HCALHighEnergyHPDFilter_cfi::energy, PVValHelper::eta, HGCalSimHitValidation::energysum::eTime, HGCalSimHitValidation::energysum::etotal, fillHitsInfo(), fillMuonTomoHistos(), fillOccupancyMap(), firstLayer_, hgcons_, HitOccupancy_Minus_, HitOccupancy_Plus_, hfClusterShapes_cfi::hits, mps_fire::i, HGCScintillatorDetId::ietaAbs(), HGCScintillatorDetId::iphi(), dqmdumpme::k, HGCScintillatorDetId::layer(), HGCSiliconDetId::layer(), phase1PixelTopology::layer, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellTrap(), nameDetector_, nTimes_, or, hgcalPerformanceValidation::orient, HGCalDDDConstants::tileTrapezoid(), protons_cff::time, times_, HGCScintillatorDetId::type(), HGCSiliconDetId::type(), HGCalTestNumbering::unpackHexagonIndex(), verbosity_, HGCalDDDConstants::waferHexagon8(), HGCalDDDConstants::waferType(), HGCSiliconDetId::waferU(), HGCSiliconDetId::waferV(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xy, HGCScintillatorDetId::zside(), HGCSiliconDetId::zside(), and ecaldqm::zside().

Referenced by analyze().

                                                                  {
   std::map<int, int> OccupancyMap_plus, OccupancyMap_minus;
   OccupancyMap_plus.clear();
   OccupancyMap_minus.clear();
 
   std::map<uint32_t, std::pair<hitsinfo, energysum> > map_hits;
   map_hits.clear();
 
   if (verbosity_ > 0)
     edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " << hits.size() << " PcaloHit elements";
   unsigned int nused(0);
   for (unsigned int i = 0; i < hits.size(); i++) {
     double energy = hits[i].energy();
     double time = hits[i].time();
     uint32_t id_ = hits[i].id();
     int cell, sector, subsector(0), layer, zside;
     int cell2(0), type(0);
     if (hgcons_->waferHexagon8()) {
       HGCSiliconDetId detId = HGCSiliconDetId(id_);
       cell = detId.cellU();
       cell2 = detId.cellV();
       sector = detId.waferU();
       subsector = detId.waferV();
       type = detId.type();
       layer = detId.layer();
       zside = detId.zside();
     } else if (hgcons_->tileTrapezoid()) {
       HGCScintillatorDetId detId = HGCScintillatorDetId(id_);
       sector = detId.ietaAbs();
       cell = detId.iphi();
       subsector = 1;
       type = detId.type();
       layer = detId.layer();
       zside = detId.zside();
     } else {
       int subdet;
       HGCalTestNumbering::unpackHexagonIndex(id_, subdet, zside, layer, sector, type, cell);
     }
     nused++;
     if (verbosity_ > 1)
       edm::LogVerbatim("HGCalValidation")
           << "Detector " << nameDetector_ << " zside = " << zside << " sector|wafer = " << sector << ":" << subsector
           << " type = " << type << " layer = " << layer << " cell = " << cell << ":" << cell2 << " energy = " << energy
           << " energyem = " << hits[i].energyEM() << " energyhad = " << hits[i].energyHad() << " time = " << time;
 
     HepGeom::Point3D<float> gcoord;
     std::pair<float, float> xy;
     if (hgcons_->waferHexagon8()) {
       xy = hgcons_->locateCell(layer, sector, subsector, cell, cell2, false, true);
     } else if (hgcons_->tileTrapezoid()) {
       xy = hgcons_->locateCellTrap(layer, sector, cell, false);
     } else {
       xy = hgcons_->locateCell(cell, layer, sector, false);
     }
     double zp = hgcons_->waferZ(layer, false);
     if (zside < 0)
       zp = -zp;
     float xp = (zp < 0) ? -xy.first : xy.first;
     gcoord = HepGeom::Point3D<float>(xp, xy.second, zp);
     double tof = (gcoord.mag() * CLHEP::mm) / CLHEP::c_light;
     if (verbosity_ > 1)
       edm::LogVerbatim("HGCalValidation")
           << std::hex << id_ << std::dec << " global coordinate " << gcoord << " time " << time << ":" << tof;
     time -= tof;
 
     energysum esum;
     hitsinfo hinfo;
     if (map_hits.count(id_) != 0) {
       hinfo = map_hits[id_].first;
       esum = map_hits[id_].second;
     } else {
       hinfo.x = gcoord.x();
       hinfo.y = gcoord.y();
       hinfo.z = gcoord.z();
       hinfo.sector = sector;
       hinfo.sector2 = subsector;
       hinfo.cell = cell;
       hinfo.cell2 = cell2;
       hinfo.type = type;
       hinfo.layer = layer - firstLayer_;
       hinfo.phi = gcoord.getPhi();
       hinfo.eta = gcoord.getEta();
     }
     esum.etotal += energy;
     for (unsigned int k = 0; k < nTimes_; ++k) {
       if (time > 0 && time < times_[k])
         esum.eTime[k] += energy;
     }
 
     if (verbosity_ > 1)
       edm::LogVerbatim("HGCalValidation") << " -----------------------   gx = " << hinfo.x << " gy = " << hinfo.y
                                           << " gz = " << hinfo.z << " phi = " << hinfo.phi << " eta = " << hinfo.eta;
     map_hits[id_] = std::pair<hitsinfo, energysum>(hinfo, esum);
   }
   if (verbosity_ > 0)
     edm::LogVerbatim("HGCalValidation") << nameDetector_ << " with " << map_hits.size()
                                         << " detector elements being hit";
 
   std::map<uint32_t, std::pair<hitsinfo, energysum> >::iterator itr;
   for (itr = map_hits.begin(); itr != map_hits.end(); ++itr) {
     hitsinfo hinfo = (*itr).second.first;
     energysum esum = (*itr).second.second;
     int layer = hinfo.layer;
     double eta = hinfo.eta;
     int type, part, orient;
     int partialType = -1;
     if (nameDetector_ == "HGCalEESensitive" or nameDetector_ == "HGCalHESiliconSensitive") {
       HGCSiliconDetId detId = HGCSiliconDetId((*itr).first);
       std::tie(type, part, orient) = hgcons_->waferType(detId);
       partialType = part;
     }
 
     for (unsigned int itimeslice = 0; itimeslice < nTimes_; itimeslice++) {
       fillHitsInfo((*itr).second, itimeslice, esum.eTime[itimeslice]);
     }
 
     if (eta > 0.0)
       fillOccupancyMap(OccupancyMap_plus, layer);
     else
       fillOccupancyMap(OccupancyMap_minus, layer);
 
     fillMuonTomoHistos(partialType, (*itr).second);
   }
   if (verbosity_ > 0)
     edm::LogVerbatim("HGCalValidation") << "With map:used:total " << hits.size() << "|" << nused << "|"
                                         << map_hits.size() << " hits";
 
   for (auto const& itr : OccupancyMap_plus) {
     int layer = itr.first;
     int occupancy = itr.second;
     HitOccupancy_Plus_.at(layer)->Fill(occupancy);
   }
   for (auto const& itr : OccupancyMap_minus) {
     int layer = itr.first;
     int occupancy = itr.second;
     HitOccupancy_Minus_.at(layer)->Fill(occupancy);
   }
 }

◆ bookHistograms()

void HGCalSimHitValidation::bookHistograms	(	DQMStore::IBooker &	iB,
		edm::Run const &	,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 511 of file HGCalSimHitValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), createHisto(), energy_, energyFWCK_, energyFWCN_, energyFWF_, energyPWCK_, energyPWCN_, energyPWF_, EtaPhi_Minus_, EtaPhi_Plus_, firstLayer_, histoSetting(), HitOccupancy_Minus_, HitOccupancy_Plus_, hitXYB_, hitXYFWCK_, hitXYFWCN_, hitXYFWF_, layers_, MeanHitOccupancy_Minus_, MeanHitOccupancy_Plus_, nameDetector_, nTimes_, or, pi, dqm::implementation::NavigatorBase::setCurrentFolder(), and cond::impl::to_string().

                                                                                                    {
   iB.setCurrentFolder("HGCAL/HGCalSimHitsV/" + nameDetector_);
 
   std::ostringstream histoname;
   for (unsigned int il = 0; il < layers_; ++il) {
     int ilayer = firstLayer_ + static_cast<int>(il);
     auto istr1 = std::to_string(ilayer);
     while (istr1.size() < 2) {
       istr1.insert(0, "0");
     }
     histoname.str("");
     histoname << "HitOccupancy_Plus_layer_" << istr1;
     HitOccupancy_Plus_.push_back(iB.book1D(histoname.str().c_str(), "HitOccupancy_Plus", 501, -0.5, 500.5));
     histoname.str("");
     histoname << "HitOccupancy_Minus_layer_" << istr1;
     HitOccupancy_Minus_.push_back(iB.book1D(histoname.str().c_str(), "HitOccupancy_Minus", 501, -0.5, 500.5));
 
     histoname.str("");
     histoname << "EtaPhi_Plus_"
               << "layer_" << istr1;
     EtaPhi_Plus_.push_back(iB.book2D(histoname.str().c_str(), "Occupancy", 31, 1.45, 3.0, 72, -CLHEP::pi, CLHEP::pi));
     histoname.str("");
     histoname << "EtaPhi_Minus_"
               << "layer_" << istr1;
     EtaPhi_Minus_.push_back(
         iB.book2D(histoname.str().c_str(), "Occupancy", 31, -3.0, -1.45, 72, -CLHEP::pi, CLHEP::pi));
 
     for (unsigned int itimeslice = 0; itimeslice < nTimes_; itimeslice++) {
       histoname.str("");
       histoname << "energy_time_" << itimeslice << "_layer_" << istr1;
       energy_[itimeslice].push_back(iB.book1D(histoname.str().c_str(), "energy_", 100, 0, 0.1));
     }
 
     if (nameDetector_ == "HGCalEESensitive" or nameDetector_ == "HGCalHESiliconSensitive") {
       histoname.str("");
       histoname << "energy_FullWafer_Fine_layer_" << istr1;
       TH1F* hEdepFWF = createHisto(histoname.str(), 100, 0., 400., false);
       histoSetting(hEdepFWF, "Eloss (keV)", "", kRed, kRed, 2);
       energyFWF_.push_back(iB.book1D(histoname.str().c_str(), hEdepFWF));
       hEdepFWF->Delete();
 
       histoname.str("");
       histoname << "energy_FullWafer_CoarseThin_layer_" << istr1;
       TH1F* hEdepFWCN = createHisto(histoname.str(), 100, 0., 400., false);
       histoSetting(hEdepFWCN, "Eloss (keV)", "", kGreen + 1, kGreen + 1, 2);
       energyFWCN_.push_back(iB.book1D(histoname.str().c_str(), hEdepFWCN));
       hEdepFWCN->Delete();
 
       histoname.str("");
       histoname << "energy_FullWafer_CoarseThick_layer_" << istr1;
       TH1F* hEdepFWCK = createHisto(histoname.str(), 100, 0., 400., false);
       histoSetting(hEdepFWCK, "Eloss (keV)", "", kMagenta, kMagenta, 2);
       energyFWCK_.push_back(iB.book1D(histoname.str().c_str(), hEdepFWCK));
       hEdepFWCK->Delete();
     }
 
 
     if (nameDetector_ == "HGCalEESensitive" or nameDetector_ == "HGCalHESiliconSensitive") {
       histoname.str("");
       histoname << "energy_PartialWafer_Fine_layer_" << istr1;
       TH1F* hEdepPWF = createHisto(histoname.str(), 100, 0., 400., false);
       histoSetting(hEdepPWF, "Eloss (keV)", "", kRed, kRed, 2);
       energyPWF_.push_back(iB.book1D(histoname.str().c_str(), hEdepPWF));
       hEdepPWF->Delete();
 
       histoname.str("");
       histoname << "energy_PartialWafer_CoarseThin_layer_" << istr1;
       TH1F* hEdepPWCN = createHisto(histoname.str(), 100, 0., 400., false);
       histoSetting(hEdepPWCN, "Eloss (keV)", "", kGreen + 1, kGreen + 1, 2);
       energyPWCN_.push_back(iB.book1D(histoname.str().c_str(), hEdepPWCN));
       hEdepPWCN->Delete();
 
       histoname.str("");
       histoname << "energy_PartialWafer_CoarseThick_layer_" << istr1;
       TH1F* hEdepPWCK = createHisto(histoname.str(), 100, 0., 400., false);
       histoSetting(hEdepPWCK, "Eloss (keV)", "", kMagenta, kMagenta, 2);
       energyPWCK_.push_back(iB.book1D(histoname.str().c_str(), hEdepPWCK));
       hEdepPWCK->Delete();
     }
 
     // ///////////// Histograms for the XY distribution of fired cells/scintillator tiles ///////////////
     if (nameDetector_ == "HGCalEESensitive" or nameDetector_ == "HGCalHESiliconSensitive") {
       histoname.str("");
       histoname << "hitXY_FullWafer_Fine_layer_" << istr1;
       TH2F* hitXYFWF = new TH2F(
           Form("hitXYFWF_%s", histoname.str().c_str()), histoname.str().c_str(), 100, -300., 300., 100, -300., 300.);
       histoSetting(hitXYFWF, "x (cm)", "y (cm)", kRed, kRed);
       hitXYFWF_.push_back(iB.book2D(histoname.str().c_str(), hitXYFWF));
       hitXYFWF->Delete();
 
       histoname.str("");
       histoname << "hitXY_FullWafer_CoarseThin_layer_" << istr1;
       TH2F* hitXYFWCN = new TH2F(
           Form("hitXYFWCN_%s", histoname.str().c_str()), histoname.str().c_str(), 100, -300., 300., 100, -300., 300.);
       histoSetting(hitXYFWCN, "x (cm)", "y (cm)", kGreen + 1, kGreen + 1);
       hitXYFWCN_.push_back(iB.book2D(histoname.str().c_str(), hitXYFWCN));
       hitXYFWCN->Delete();
 
       histoname.str("");
       histoname << "hitXY_FullWafer_CoarseThick_layer_" << istr1;
       TH2F* hitXYFWCK = new TH2F(
           Form("hitXYFWCK_%s", histoname.str().c_str()), histoname.str().c_str(), 100, -300., 300., 100, -300., 300.);
       histoSetting(hitXYFWCK, "x (cm)", "y (cm)", kMagenta, kMagenta);
       hitXYFWCK_.push_back(iB.book2D(histoname.str().c_str(), hitXYFWCK));
       hitXYFWCK->Delete();
     }
 
     if (nameDetector_ == "HGCalHEScintillatorSensitive") {
       histoname.str("");
       histoname << "hitXY_Scintillator_layer_" << istr1;
       TH2F* hitXYB = new TH2F(
           Form("hitXYB_%s", histoname.str().c_str()), histoname.str().c_str(), 100, -300., 300., 100, -300., 300.);
       histoSetting(hitXYB, "x (cm)", "y (cm)", kBlue, kBlue);
       hitXYB_.push_back(iB.book2D(histoname.str().c_str(), hitXYB));
       hitXYB->Delete();
     }
   }
   MeanHitOccupancy_Plus_ = iB.book1D("MeanHitOccupancy_Plus", "MeanHitOccupancy_Plus", layers_, 0.5, layers_ + 0.5);
   MeanHitOccupancy_Minus_ = iB.book1D("MeanHitOccupancy_Minus", "MeanHitOccupancy_Minus", layers_, 0.5, layers_ + 0.5);
 }

◆ createHisto()

TH1F * HGCalSimHitValidation::createHisto	(	std::string	histname,
		const int	nbins,
		float	minIndexX,
		float	maxIndexX,
		bool	isLogX = `true`
	)

private

Definition at line 637 of file HGCalSimHitValidation.cc.

References PVValHelper::dx, gpuVertexFinder::hist, mps_fire::i, LaserClient_cfi::nbins, and fw3dlego::xbins.

Referenced by bookHistograms().

                                                                                         {
   TH1F* hist = nullptr;
   if (isLogX) {
     Double_t xbins[nbins + 1];
     double dx = (maxIndexX - minIndexX) / nbins;
     for (int i = 0; i <= nbins; i++) {
       xbins[i] = TMath::Power(10, (minIndexX + i * dx));
     }
     hist = new TH1F(Form("hEdep_%s", histname.c_str()), histname.c_str(), nbins, xbins);
   } else {
     hist = new TH1F(Form("hEdep_%s", histname.c_str()), histname.c_str(), nbins, minIndexX, maxIndexX);
   }
   return hist;
 }

◆ defineGeometry() [1/2]

bool HGCalSimHitValidation::defineGeometry ( const DDCompactView * ddViewH )

private

Definition at line 404 of file HGCalSimHitValidation.cc.

References DDTrap::alpha1(), filterCSVwithJSON::copy, DDFilteredView::copyNumbers(), ALCARECOTkAlBeamHalo_cff::filter, DDFilteredView::firstChild(), HGCalDDDConstants::geomMode(), hgcons_, DDFilteredView::logicalPart(), Skims_PA_cff::name, nameDetector_, DDFilteredView::next(), HGCalTestNumbering::packSquareIndex(), idealTransformation::rotation, DDFilteredView::rotation(), fileinputsource_cfi::sec, mkfit::Const::sol, DDLogicalPart::solid(), HGCalGeometryMode::Square, AlCaHLTBitMon_QueryRunRegistry::string, symmDet_, DDFilteredView::translation(), transMap_, funct::true, verbosity_, x, y, and z.

Referenced by dqmBeginRun().

                                                                        {
   if (verbosity_ > 0)
     edm::LogVerbatim("HGCalValidation") << "Initialize HGCalDDDConstants (DDD) for " << nameDetector_ << " : "
                                         << hgcons_;
 
   if (hgcons_->geomMode() == HGCalGeometryMode::Square) {
     const DDCompactView& cview = *ddViewH;
     std::string attribute = "Volume";
     std::string value = nameDetector_;
 
     DDSpecificsMatchesValueFilter filter{DDValue(attribute, value, 0)};
     DDFilteredView fv(cview, filter);
     bool dodet = fv.firstChild();
 
     while (dodet) {
       const DDSolid& sol = fv.logicalPart().solid();
       const std::string& name = sol.name().fullname();
       int isd = (name.find(nameDetector_) == std::string::npos) ? -1 : 1;
       if (isd > 0) {
         std::vector<int> copy = fv.copyNumbers();
         int nsiz = static_cast<int>(copy.size());
         int lay = (nsiz > 0) ? copy[nsiz - 1] : -1;
         int sec = (nsiz > 1) ? copy[nsiz - 2] : -1;
         int zp = (nsiz > 3) ? copy[nsiz - 4] : -1;
         if (zp != 1)
           zp = -1;
         const DDTrap& trp = static_cast<DDTrap>(sol);
         int subs = (trp.alpha1() > 0 ? 1 : 0);
         symmDet_ = (trp.alpha1() == 0 ? true : false);
         uint32_t id = HGCalTestNumbering::packSquareIndex(zp, lay, sec, subs, 0);
         DD3Vector x, y, z;
         fv.rotation().GetComponents(x, y, z);
         const CLHEP::HepRep3x3 rotation(x.X(), y.X(), z.X(), x.Y(), y.Y(), z.Y(), x.Z(), y.Z(), z.Z());
         const CLHEP::HepRotation hr(rotation);
         const CLHEP::Hep3Vector h3v(fv.translation().X(), fv.translation().Y(), fv.translation().Z());
         const HepGeom::Transform3D ht3d(hr, h3v);
         transMap_.insert(std::make_pair(id, ht3d));
         if (verbosity_ > 2)
           edm::LogVerbatim("HGCalValidation") << HGCalDetId(id) << " Transform using " << h3v << " and " << hr;
       }
       dodet = fv.next();
     }
     if (verbosity_ > 0)
       edm::LogVerbatim("HGCalValidation") << "Finds " << transMap_.size() << " elements and SymmDet_ = " << symmDet_;
   }
   return true;
 }

◆ defineGeometry() [2/2]

bool HGCalSimHitValidation::defineGeometry ( const cms::DDCompactView * ddViewH )

private

Definition at line 452 of file HGCalSimHitValidation.cc.

References filterCSVwithJSON::copy, cms::DDFilteredView::copyNos(), ALCARECOTkAlBeamHalo_cff::filter, cms::DDFilteredView::firstChild(), HGCalDDDConstants::geomMode(), hgcons_, Skims_PA_cff::name, cms::DDFilteredView::name(), nameDetector_, HGCalTestNumbering::packSquareIndex(), cms::DDFilteredView::parameters(), idealTransformation::rotation, cms::DDFilteredView::rotation(), fileinputsource_cfi::sec, HGCalGeometryMode::Square, symmDet_, cms::DDFilteredView::translation(), transMap_, funct::true, verbosity_, x, y, and z.

                                                                           {
   if (verbosity_ > 0)
     edm::LogVerbatim("HGCalValidation") << "Initialize HGCalDDDConstants (DD4hep) for " << nameDetector_ << " : "
                                         << hgcons_;
 
   if (hgcons_->geomMode() == HGCalGeometryMode::Square) {
     const cms::DDCompactView& cview = *ddViewH;
     const cms::DDFilter filter("Volume", nameDetector_);
     cms::DDFilteredView fv(cview, filter);
 
     while (fv.firstChild()) {
       const auto& name = fv.name();
       int isd = (name.find(nameDetector_) == std::string::npos) ? -1 : 1;
       if (isd > 0) {
         const auto& copy = fv.copyNos();
         int nsiz = static_cast<int>(copy.size());
         int lay = (nsiz > 0) ? copy[0] : -1;
         int sec = (nsiz > 1) ? copy[1] : -1;
         int zp = (nsiz > 3) ? copy[3] : -1;
         if (zp != 1)
           zp = -1;
         const auto& pars = fv.parameters();
         int subs = (pars[6] > 0 ? 1 : 0);
         symmDet_ = (pars[6] == 0 ? true : false);
         uint32_t id = HGCalTestNumbering::packSquareIndex(zp, lay, sec, subs, 0);
         DD3Vector x, y, z;
         fv.rotation().GetComponents(x, y, z);
         const CLHEP::HepRep3x3 rotation(x.X(), y.X(), z.X(), x.Y(), y.Y(), z.Y(), x.Z(), y.Z(), z.Z());
         const CLHEP::HepRotation hr(rotation);
         const CLHEP::Hep3Vector h3v(fv.translation().X(), fv.translation().Y(), fv.translation().Z());
         const HepGeom::Transform3D ht3d(hr, h3v);
         transMap_.insert(std::make_pair(id, ht3d));
         if (verbosity_ > 2)
           edm::LogVerbatim("HGCalValidation") << HGCalDetId(id) << " Transform using " << h3v << " and " << hr;
       }
     }
     if (verbosity_ > 0)
       edm::LogVerbatim("HGCalValidation") << "Finds " << transMap_.size() << " elements and SymmDet_ = " << symmDet_;
   }
   return true;
 }

◆ dqmBeginRun()

void HGCalSimHitValidation::dqmBeginRun	(	const edm::Run &	,
		const edm::EventSetup &	iSetup
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 495 of file HGCalSimHitValidation.cc.

References defineGeometry(), HGCalDDDConstants::firstLayer(), firstLayer_, fromDDD_, edm::EventSetup::getData(), hgcons_, HGCalDDDConstants::layers(), layers_, nameDetector_, tok_cpv_, tok_cpvc_, tok_hgcal_, and verbosity_.

                                                                                   {
   hgcons_ = &iSetup.getData(tok_hgcal_);
   layers_ = hgcons_->layers(false);
   firstLayer_ = hgcons_->firstLayer();
   if (fromDDD_) {
     const DDCompactView* pDD = &iSetup.getData(tok_cpv_);
     defineGeometry(pDD);
   } else {
     const cms::DDCompactView* pDD = &iSetup.getData(tok_cpvc_);
     defineGeometry(pDD);
   }
   if (verbosity_ > 0)
     edm::LogVerbatim("HGCalValidation") << nameDetector_ << " defined with " << layers_ << " Layers with first at "
                                         << firstLayer_;
 }

◆ fillDescriptions()

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

static

Definition at line 145 of file HGCalSimHitValidation.cc.

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

                                                                                        {
   edm::ParameterSetDescription desc;
   std::vector<double> times = {25.0, 1000.0};
   desc.add<std::string>("DetectorName", "HGCalEESensitive");
   desc.add<std::string>("CaloHitSource", "HGCHitsEE");
   desc.add<std::vector<double> >("TimeSlices", times);
   desc.addUntracked<int>("Verbosity", 0);
   desc.addUntracked<bool>("TestNumber", true);
   desc.addUntracked<bool>("fromDDD", true);
   descriptions.add("hgcalSimHitValidationEE", desc);
 }

◆ fillHitsInfo()

void HGCalSimHitValidation::fillHitsInfo	(	std::pair< hitsinfo, energysum >	hit_,
		unsigned int	itimeslice,
		double	esum
	)

private

Definition at line 336 of file HGCalSimHitValidation.cc.

References energy_, EtaPhi_Minus_, EtaPhi_Plus_, hfClusterShapes_cfi::hits, layers_, nameDetector_, and verbosity_.

Referenced by analyzeHits().

                                                                                                                 {
   unsigned int ilayer = hits.first.layer;
   if (ilayer < layers_) {
     energy_[itimeslice].at(ilayer)->Fill(esum);
     if (itimeslice == 0) {
       EtaPhi_Plus_.at(ilayer)->Fill(hits.first.eta, hits.first.phi);
       EtaPhi_Minus_.at(ilayer)->Fill(hits.first.eta, hits.first.phi);
     }
   } else {
     if (verbosity_ > 0)
       edm::LogVerbatim("HGCalValidation")
           << "Problematic Hit for " << nameDetector_ << " at sector " << hits.first.sector << ":" << hits.first.sector2
           << " layer " << hits.first.layer << " cell " << hits.first.cell << ":" << hits.first.cell2 << " energy "
           << hits.second.etotal;
   }
 }

◆ fillMuonTomoHistos()

void HGCalSimHitValidation::fillMuonTomoHistos	(	int	partialType,
		std::pair< hitsinfo, energysum >	hit_
	)

private

Definition at line 353 of file HGCalSimHitValidation.cc.

References MillePedeFileConverter_cfg::e, energyFWCK_, energyFWCN_, energyFWF_, energyPWCK_, energyPWCN_, energyPWF_, HGCalSimHitValidation::energysum::eTime, HGCSiliconDetId::HGCalCoarseThick, HGCSiliconDetId::HGCalCoarseThin, HGCSiliconDetId::HGCalFine, hfClusterShapes_cfi::hits, hitXYB_, hitXYFWCK_, hitXYFWCN_, hitXYFWF_, layers_, nameDetector_, or, x, and y.

Referenced by analyzeHits().

                                                                                                  {
   hitsinfo hinfo = hits.first;
   energysum esum = hits.second;
   double edep =
       esum.eTime[0] * CLHEP::GeV /
       CLHEP::keV;  //index 0 and 1 corresponds to 25 ns and 1000 ns, respectively. In addititon, chaging energy loss unit to keV.
 
   unsigned int ilayer = hinfo.layer;
   double x = hinfo.x * CLHEP::mm / CLHEP::cm;  // chaging length unit to cm.
   double y = hinfo.y * CLHEP::mm / CLHEP::cm;
   if (ilayer < layers_) {
     if (nameDetector_ == "HGCalEESensitive" or nameDetector_ == "HGCalHESiliconSensitive") {
       // Fill the energy loss histograms for MIP
       if (!TMath::AreEqualAbs(edep, 0.0, 1.e-5)) {  //to avoid peak at zero due Eloss less than 10 mili eV.
         if (hinfo.type == HGCSiliconDetId::HGCalFine) {
           if (partialType == 0)
             energyFWF_.at(ilayer)->Fill(edep);
           if (partialType > 0)
             energyPWF_.at(ilayer)->Fill(edep);
         }
         if (hinfo.type == HGCSiliconDetId::HGCalCoarseThin) {
           if (partialType == 0)
             energyFWCN_.at(ilayer)->Fill(edep);
           if (partialType > 0)
             energyPWCN_.at(ilayer)->Fill(edep);
         }
         if (hinfo.type == HGCSiliconDetId::HGCalCoarseThick) {
           if (partialType == 0)
             energyFWCK_.at(ilayer)->Fill(edep);
           if (partialType > 0)
             energyPWCK_.at(ilayer)->Fill(edep);
         }
       }
 
       // Fill the XY distribution of detector hits
       if (hinfo.type == HGCSiliconDetId::HGCalFine)
         hitXYFWF_.at(ilayer)->Fill(x, y);
 
       if (hinfo.type == HGCSiliconDetId::HGCalCoarseThin)
         hitXYFWCN_.at(ilayer)->Fill(x, y);
 
       if (hinfo.type == HGCSiliconDetId::HGCalCoarseThick)
         hitXYFWCK_.at(ilayer)->Fill(x, y);
 
     }  //is Silicon
     if (nameDetector_ == "HGCalHEScintillatorSensitive") {
       hitXYB_.at(ilayer)->Fill(x, y);
     }  //is Scintillator
   }    //layer condition
 }

◆ fillOccupancyMap()

void HGCalSimHitValidation::fillOccupancyMap	(	std::map< int, int > &	OccupancyMap,
		int	layer
	)

private

Definition at line 328 of file HGCalSimHitValidation.cc.

References phase1PixelTopology::layer.

Referenced by analyzeHits().

                                                                                       {
   if (OccupancyMap.find(layer) != OccupancyMap.end()) {
     ++OccupancyMap[layer];
   } else {
     OccupancyMap[layer] = 1;
   }
 }

◆ histoSetting() [1/2]

void HGCalSimHitValidation::histoSetting	(	TH1F *&	histo,
		const char *	xTitle,
		const char *	yTitle = `""`,
		Color_t	lineColor = `kBlack`,
		Color_t	markerColor = `kBlack`,
		int	linewidth = `1`
	)

private

Definition at line 653 of file HGCalSimHitValidation.cc.

References timingPdfMaker::histo, RecoTauValidation_cfi::lineColor, RecoTauValidation_cfi::lineWidth, and RecoTauValidation_cfi::markerColor.

Referenced by bookHistograms().

                                                                                                                  {
   histo->SetStats();
   histo->SetLineColor(lineColor);
   histo->SetLineWidth(lineWidth);
   histo->SetMarkerColor(markerColor);
   histo->GetXaxis()->SetTitle(xTitle);
   histo->GetYaxis()->SetTitle(yTitle);
 }

◆ histoSetting() [2/2]

void HGCalSimHitValidation::histoSetting	(	TH2F *&	histo,
		const char *	xTitle,
		const char *	yTitle = `""`,
		Color_t	lineColor = `kBlack`,
		Color_t	markerColor = `kBlack`,
		int	linewidth = `1`
	)

private

Definition at line 663 of file HGCalSimHitValidation.cc.

References timingPdfMaker::histo, RecoTauValidation_cfi::lineColor, RecoTauValidation_cfi::lineWidth, and RecoTauValidation_cfi::markerColor.

                                                                                                                  {
   histo->SetStats();
   histo->SetLineColor(lineColor);
   histo->SetLineWidth(lineWidth);
   histo->SetMarkerColor(markerColor);
   histo->SetMarkerStyle(kFullCircle);
   histo->SetMarkerSize(0.6);
   histo->GetXaxis()->SetTitle(xTitle);
   histo->GetYaxis()->SetTitle(yTitle);
 }

Member Data Documentation

◆ caloHitSource_

const std::string HGCalSimHitValidation::caloHitSource_

private

Definition at line 102 of file HGCalSimHitValidation.cc.

◆ energy_

std::vector<MonitorElement*> HGCalSimHitValidation::energy_[maxTime_]

private

Definition at line 121 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

◆ energyFWCK_

std::vector<MonitorElement*> HGCalSimHitValidation::energyFWCK_

private

Definition at line 122 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ energyFWCN_

std::vector<MonitorElement*> HGCalSimHitValidation::energyFWCN_

private

Definition at line 122 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ energyFWF_

std::vector<MonitorElement*> HGCalSimHitValidation::energyFWF_

private

Definition at line 122 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ energyPWCK_

std::vector<MonitorElement*> HGCalSimHitValidation::energyPWCK_

private

Definition at line 123 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ energyPWCN_

std::vector<MonitorElement*> HGCalSimHitValidation::energyPWCN_

private

Definition at line 123 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ energyPWF_

std::vector<MonitorElement*> HGCalSimHitValidation::energyPWF_

private

Definition at line 123 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ EtaPhi_Minus_

std::vector<MonitorElement*> HGCalSimHitValidation::EtaPhi_Minus_

private

Definition at line 118 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

◆ EtaPhi_Plus_

std::vector<MonitorElement*> HGCalSimHitValidation::EtaPhi_Plus_

private

Definition at line 118 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillHitsInfo().

◆ firstLayer_

int HGCalSimHitValidation::firstLayer_

private

Definition at line 114 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), bookHistograms(), and dqmBeginRun().

◆ fromDDD_

const bool HGCalSimHitValidation::fromDDD_

private

Definition at line 106 of file HGCalSimHitValidation.cc.

Referenced by dqmBeginRun().

◆ hgcons_

const HGCalDDDConstants* HGCalSimHitValidation::hgcons_

private

Definition at line 103 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), defineGeometry(), and dqmBeginRun().

◆ HitOccupancy_Minus_

std::vector<MonitorElement*> HGCalSimHitValidation::HitOccupancy_Minus_

private

Definition at line 117 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), and bookHistograms().

◆ HitOccupancy_Plus_

std::vector<MonitorElement*> HGCalSimHitValidation::HitOccupancy_Plus_

private

Definition at line 117 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), and bookHistograms().

◆ hitXYB_

std::vector<MonitorElement*> HGCalSimHitValidation::hitXYB_

private

Definition at line 124 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ hitXYFWCK_

std::vector<MonitorElement*> HGCalSimHitValidation::hitXYFWCK_

private

Definition at line 124 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ hitXYFWCN_

std::vector<MonitorElement*> HGCalSimHitValidation::hitXYFWCN_

private

Definition at line 124 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ hitXYFWF_

std::vector<MonitorElement*> HGCalSimHitValidation::hitXYFWF_

private

Definition at line 124 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), and fillMuonTomoHistos().

◆ layers_

unsigned int HGCalSimHitValidation::layers_

private

Definition at line 113 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms(), dqmBeginRun(), fillHitsInfo(), and fillMuonTomoHistos().

◆ maxTime_

const unsigned int HGCalSimHitValidation::maxTime_ = 6

staticprivate

Definition at line 120 of file HGCalSimHitValidation.cc.

◆ MeanHitOccupancy_Minus_

MonitorElement * HGCalSimHitValidation::MeanHitOccupancy_Minus_

private

Definition at line 119 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms().

◆ MeanHitOccupancy_Plus_

MonitorElement* HGCalSimHitValidation::MeanHitOccupancy_Plus_

private

Definition at line 119 of file HGCalSimHitValidation.cc.

Referenced by bookHistograms().

◆ nameDetector_

const std::string HGCalSimHitValidation::nameDetector_

private

Definition at line 102 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), bookHistograms(), defineGeometry(), dqmBeginRun(), fillHitsInfo(), fillMuonTomoHistos(), and HGCalSimHitValidation().

◆ nTimes_

unsigned int HGCalSimHitValidation::nTimes_

private

Definition at line 125 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits(), and bookHistograms().

◆ symmDet_

bool HGCalSimHitValidation::symmDet_

private

Definition at line 112 of file HGCalSimHitValidation.cc.

Referenced by defineGeometry().

◆ times_

const std::vector<double> HGCalSimHitValidation::times_

private

Definition at line 104 of file HGCalSimHitValidation.cc.

Referenced by analyzeHits().

◆ tok_cpv_

const edm::ESGetToken<DDCompactView, IdealGeometryRecord> HGCalSimHitValidation::tok_cpv_

private

Definition at line 108 of file HGCalSimHitValidation.cc.

Referenced by dqmBeginRun().

◆ tok_cpvc_

const edm::ESGetToken<cms::DDCompactView, IdealGeometryRecord> HGCalSimHitValidation::tok_cpvc_

private

Definition at line 109 of file HGCalSimHitValidation.cc.

Referenced by dqmBeginRun().

◆ tok_hepMC_

edm::EDGetTokenT<edm::HepMCProduct> HGCalSimHitValidation::tok_hepMC_

private

Definition at line 111 of file HGCalSimHitValidation.cc.

Referenced by analyze().

◆ tok_hgcal_

const edm::ESGetToken<HGCalDDDConstants, IdealGeometryRecord> HGCalSimHitValidation::tok_hgcal_

private

Definition at line 107 of file HGCalSimHitValidation.cc.

Referenced by dqmBeginRun().

◆ tok_hits_

edm::EDGetTokenT<edm::PCaloHitContainer> HGCalSimHitValidation::tok_hits_

private

Definition at line 110 of file HGCalSimHitValidation.cc.

Referenced by analyze().

◆ transMap_

std::map<uint32_t, HepGeom::Transform3D> HGCalSimHitValidation::transMap_

private

Definition at line 115 of file HGCalSimHitValidation.cc.

Referenced by defineGeometry().

◆ verbosity_

const int HGCalSimHitValidation::verbosity_

private

Definition at line 105 of file HGCalSimHitValidation.cc.

Referenced by analyze(), analyzeHits(), defineGeometry(), dqmBeginRun(), and fillHitsInfo().

Classes

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGCalSimHitValidation()

◆ ~HGCalSimHitValidation()

Member Function Documentation

◆ analyze()

◆ analyzeHits()

◆ bookHistograms()

◆ createHisto()

◆ defineGeometry() [1/2]

◆ defineGeometry() [2/2]

◆ dqmBeginRun()

◆ fillDescriptions()

◆ fillHitsInfo()

◆ fillMuonTomoHistos()

◆ fillOccupancyMap()

◆ histoSetting() [1/2]

◆ histoSetting() [2/2]

Member Data Documentation

◆ caloHitSource_

◆ energy_

◆ energyFWCK_

◆ energyFWCN_

◆ energyFWF_

◆ energyPWCK_

◆ energyPWCN_

◆ energyPWF_

◆ EtaPhi_Minus_

◆ EtaPhi_Plus_

◆ firstLayer_

◆ fromDDD_

◆ hgcons_

◆ HitOccupancy_Minus_

◆ HitOccupancy_Plus_

◆ hitXYB_

◆ hitXYFWCK_

◆ hitXYFWCN_

◆ hitXYFWF_

◆ layers_

◆ maxTime_

◆ MeanHitOccupancy_Minus_

◆ MeanHitOccupancy_Plus_

◆ nameDetector_

◆ nTimes_

◆ symmDet_

◆ times_

◆ tok_cpv_

◆ tok_cpvc_

◆ tok_hepMC_

◆ tok_hgcal_

◆ tok_hits_

◆ transMap_

◆ verbosity_