Inheritance diagram for HGCGeometryValidation:

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

	~HGCGeometryValidation () override=default

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 (edm::Run const &, edm::EventSetup const &) override

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
const edm::EDGetTokenT< PHGCalValidInfo >	g4Token_

const std::vector< std::string >	geometrySource_

const std::vector< edm::ESGetToken< HGCalDDDConstants, IdealGeometryRecord > >	geomToken_

MonitorElement *	hebdX

MonitorElement *	hebdxVsLayer

MonitorElement *	hebdxVsX

MonitorElement *	hebdY

MonitorElement *	hebdyVsLayer

MonitorElement *	hebdyVsY

MonitorElement *	hebdZ

MonitorElement *	hebdzVsLayer

MonitorElement *	hebdzVsZ

MonitorElement *	hebLayerVsEnStep

MonitorElement *	hebTotEdepStep

MonitorElement *	hebXG4VsId

MonitorElement *	hebYG4VsId

MonitorElement *	hebZG4VsId

MonitorElement *	heedX

MonitorElement *	heedxVsLayer

MonitorElement *	heedxVsX

MonitorElement *	heedY

MonitorElement *	heedyVsLayer

MonitorElement *	heedyVsY

MonitorElement *	heedZ

MonitorElement *	heedzVsLayer

MonitorElement *	heedzVsZ

MonitorElement *	heeLayerVsEnStep

MonitorElement *	heeTotEdepStep

MonitorElement *	heeXG4VsId

MonitorElement *	heeYG4VsId

MonitorElement *	heeZG4VsId

MonitorElement *	hefdX

MonitorElement *	hefdxVsLayer

MonitorElement *	hefdxVsX

MonitorElement *	hefdY

MonitorElement *	hefdyVsLayer

MonitorElement *	hefdyVsY

MonitorElement *	hefdZ

MonitorElement *	hefdzVsLayer

MonitorElement *	hefdzVsZ

MonitorElement *	hefLayerVsEnStep

MonitorElement *	hefTotEdepStep

MonitorElement *	hefXG4VsId

MonitorElement *	hefYG4VsId

MonitorElement *	hefZG4VsId

std::vector< const HGCalDDDConstants * >	hgcGeometry_

const int	verbosity_

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 32 of file HGCGeometryValidation.cc.

Constructor & Destructor Documentation

◆ HGCGeometryValidation()

HGCGeometryValidation::HGCGeometryValidation ( const edm::ParameterSet & cfg )

explicit

Definition at line 70 of file HGCGeometryValidation.cc.

References geometrySource_, Skims_PA_cff::name, AlCaHLTBitMon_QueryRunRegistry::string, and edm::vector_transform().

     : g4Token_(consumes<PHGCalValidInfo>(cfg.getParameter<edm::InputTag>("g4Source"))),
       geometrySource_(cfg.getUntrackedParameter<std::vector<std::string>>("geometrySource")),
       verbosity_(cfg.getUntrackedParameter<int>("verbosity", 0)),
       geomToken_{edm::vector_transform(geometrySource_, [this](const std::string &name) {
         return esConsumes<HGCalDDDConstants, IdealGeometryRecord, edm::Transition::BeginRun>(edm::ESInputTag{"", name});
       })} {}

◆ ~HGCGeometryValidation()

HGCGeometryValidation::~HGCGeometryValidation ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 151 of file HGCGeometryValidation.cc.

References runTauDisplay::dtype, PHGCalValidInfo::eehgcEdep(), PHGCalValidInfo::eeTotEdep(), dqm::impl::MonitorElement::Fill(), g4Token_, hebdX, hebdxVsLayer, hebdxVsX, hebdY, hebdyVsLayer, hebdyVsY, hebdZ, hebdzVsLayer, hebdzVsZ, PHGCalValidInfo::hebhgcEdep(), hebLayerVsEnStep, PHGCalValidInfo::hebTotEdep(), hebTotEdepStep, hebXG4VsId, hebYG4VsId, hebZG4VsId, heedX, heedxVsLayer, heedxVsX, heedY, heedyVsLayer, heedyVsY, heedZ, heedzVsLayer, heedzVsZ, heeLayerVsEnStep, heeTotEdepStep, heeXG4VsId, heeYG4VsId, heeZG4VsId, hefdX, hefdxVsLayer, hefdxVsX, hefdY, hefdyVsLayer, hefdyVsY, hefdZ, hefdzVsLayer, hefdzVsZ, PHGCalValidInfo::hefhgcEdep(), hefLayerVsEnStep, PHGCalValidInfo::hefTotEdep(), hefTotEdepStep, hefXG4VsId, hefYG4VsId, hefZG4VsId, DetId::HGCalEE, DetId::HGCalHSi, hgcGeometry_, PHGCalValidInfo::hitDets(), PHGCalValidInfo::hitIndex(), PHGCalValidInfo::hitvtxX(), PHGCalValidInfo::hitvtxY(), PHGCalValidInfo::hitvtxZ(), mps_fire::i, l1ctLayer2EG_cff::id, ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::ietaAbs(), iEvent, ALPAKA_ACCELERATOR_NAMESPACE::caPixelDoublets::if(), hcalRecHitTable_cff::iphi, edm::HandleBase::isValid(), nano_mu_digi_cff::layer, mmtocm, verbosity_, HGCalWaferIndex::waferU(), HGCalWaferIndex::waferV(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, geometryCSVtoXML::yy, ecaldqm::zside(), and geometryCSVtoXML::zz.

                                                                                        {
   //Accessing G4 information
   const edm::Handle<PHGCalValidInfo> &infoLayer = iEvent.getHandle(g4Token_);
 
   if (infoLayer.isValid()) {
     //step vertex information
     std::vector<float> hitVtxX = infoLayer->hitvtxX();
     std::vector<float> hitVtxY = infoLayer->hitvtxY();
     std::vector<float> hitVtxZ = infoLayer->hitvtxZ();
     const std::vector<unsigned int> &hitDet = infoLayer->hitDets();
     const std::vector<unsigned int> &hitIdx = infoLayer->hitIndex();
 
     //energy information
     const std::vector<float> &edepLayerEE = infoLayer->eehgcEdep();
     const std::vector<float> &edepLayerHE = infoLayer->hefhgcEdep();
     const std::vector<float> &edepLayerHB = infoLayer->hebhgcEdep();
 
     unsigned int i;
     for (i = 0; i < edepLayerEE.size(); i++) {
       heeLayerVsEnStep->Fill(i, edepLayerEE[i]);
     }
 
     for (i = 0; i < edepLayerHE.size(); i++) {
       hefLayerVsEnStep->Fill(i, edepLayerHE[i]);
     }
 
     for (i = 0; i < edepLayerHB.size(); i++) {
       hebLayerVsEnStep->Fill(i, edepLayerHB[i]);
     }
 
     //fill total energy deposited
     heeTotEdepStep->Fill((double)infoLayer->eeTotEdep());
     hefTotEdepStep->Fill((double)infoLayer->hefTotEdep());
     hebTotEdepStep->Fill((double)infoLayer->hebTotEdep());
 
     //loop over all hits
     for (unsigned int i = 0; i < hitVtxX.size(); i++) {
       hitVtxX[i] *= mmtocm;
       hitVtxY[i] *= mmtocm;
       hitVtxZ[i] *= mmtocm;
 
       double xx, yy;
       int dtype(0), layer(0), zside(1);
       std::pair<float, float> xy;
       if ((hitDet[i] == static_cast<unsigned int>(DetId::HGCalEE)) ||
           (hitDet[i] == static_cast<unsigned int>(DetId::HGCalHSi))) {
         HGCSiliconDetId id(hitIdx[i]);
         dtype = (id.det() == DetId::HGCalEE) ? 0 : 1;
         layer = id.layer();
         zside = id.zside();
         xy = hgcGeometry_[dtype]->locateCell(
             zside, layer, id.waferU(), id.waferV(), id.cellU(), id.cellV(), true, true, false, false);
       } else {
         HGCScintillatorDetId id(hitIdx[i]);
         dtype = 2;
         layer = id.layer();
         zside = id.zside();
         xy = hgcGeometry_[dtype]->locateCellTrap(zside, layer, id.ietaAbs(), id.iphi(), true, false);
       }
       double zz = hgcGeometry_[dtype]->waferZ(layer, true);  //cm
       if (zside < 0)
         zz = -zz;
       xx = (zside < 0) ? -xy.first : xy.first;
       yy = xy.second;
 
       if (dtype == 0) {
         heedzVsZ->Fill(zz, (hitVtxZ[i] - zz));
         heedyVsY->Fill(yy, (hitVtxY[i] - yy));
         heedxVsX->Fill(xx, (hitVtxX[i] - xx));
 
         heeXG4VsId->Fill(hitVtxX[i], xx);
         heeYG4VsId->Fill(hitVtxY[i], yy);
         heeZG4VsId->Fill(hitVtxZ[i], zz);
 
         heedzVsLayer->Fill(layer, (hitVtxZ[i] - zz));
         heedyVsLayer->Fill(layer, (hitVtxY[i] - yy));
         heedxVsLayer->Fill(layer, (hitVtxX[i] - xx));
 
         heedX->Fill((hitVtxX[i] - xx));
         heedZ->Fill((hitVtxZ[i] - zz));
         heedY->Fill((hitVtxY[i] - yy));
 
       } else if (dtype == 1) {
         hefdzVsZ->Fill(zz, (hitVtxZ[i] - zz));
         hefdyVsY->Fill(yy, (hitVtxY[i] - yy));
         hefdxVsX->Fill(xx, (hitVtxX[i] - xx));
 
         hefXG4VsId->Fill(hitVtxX[i], xx);
         hefYG4VsId->Fill(hitVtxY[i], yy);
         hefZG4VsId->Fill(hitVtxZ[i], zz);
 
         hefdzVsLayer->Fill(layer, (hitVtxZ[i] - zz));
         hefdyVsLayer->Fill(layer, (hitVtxY[i] - yy));
         hefdxVsLayer->Fill(layer, (hitVtxX[i] - xx));
 
         hefdX->Fill((hitVtxX[i] - xx));
         hefdZ->Fill((hitVtxZ[i] - zz));
         hefdY->Fill((hitVtxY[i] - yy));
 
       } else {
         hebdzVsZ->Fill(zz, (hitVtxZ[i] - zz));
         hebdyVsY->Fill(yy, (hitVtxY[i] - yy));
         hebdxVsX->Fill(xx, (hitVtxX[i] - xx));
 
         hebXG4VsId->Fill(hitVtxX[i], xx);
         hebYG4VsId->Fill(hitVtxY[i], yy);
         hebZG4VsId->Fill(hitVtxZ[i], zz);
 
         hebdzVsLayer->Fill(layer, (hitVtxZ[i] - zz));
         hebdyVsLayer->Fill(layer, (hitVtxY[i] - yy));
         hebdxVsLayer->Fill(layer, (hitVtxX[i] - xx));
 
         hebdX->Fill((hitVtxX[i] - xx));
         hebdZ->Fill((hitVtxZ[i] - zz));
         hebdY->Fill((hitVtxY[i] - yy));
       }
     }  //end G4 hits
 
   } else {
     if (verbosity_ > 0)
       edm::LogVerbatim("HGCalValid") << "HGCGeometryValidation::No PHGCalInfo";
   }
 }

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 90 of file HGCGeometryValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), hebdX, hebdxVsLayer, hebdxVsX, hebdY, hebdyVsLayer, hebdyVsY, hebdZ, hebdzVsLayer, hebdzVsZ, hebLayerVsEnStep, hebTotEdepStep, hebXG4VsId, hebYG4VsId, hebZG4VsId, heedX, heedxVsLayer, heedxVsX, heedY, heedyVsLayer, heedyVsY, heedZ, heedzVsLayer, heedzVsZ, heeLayerVsEnStep, heeTotEdepStep, heeXG4VsId, heeYG4VsId, heeZG4VsId, hefdX, hefdxVsLayer, hefdxVsX, hefdY, hefdyVsLayer, hefdyVsY, hefdZ, hefdzVsLayer, hefdzVsZ, hefLayerVsEnStep, hefTotEdepStep, hefXG4VsId, hefYG4VsId, hefZG4VsId, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                                      {
   iB.setCurrentFolder("HGCAL/HGCalSimHitsV/Geometry");
 
   //initiating histograms
   heeTotEdepStep = iB.book1D("heeTotEdepStep", "", 100, 0, 100);
   hefTotEdepStep = iB.book1D("hefTotEdepStep", "", 100, 0, 100);
   hebTotEdepStep = iB.book1D("hebTotEdepStep", "", 100, 0, 100);
 
   hebLayerVsEnStep = iB.book2D("hebLayerVsEnStep", "", 25, 0, 25, 100, 0, 0.01);
   hefLayerVsEnStep = iB.book2D("hefLayerVsEnStep", "", 36, 0, 36, 100, 0, 0.01);
   heeLayerVsEnStep = iB.book2D("heeLayerVsEnStep", "", 84, 0, 84, 100, 0, 0.01);
 
   heeXG4VsId = iB.book2D("heeXG4VsId", "", 600, -300, 300, 600, -300, 300);
   heeYG4VsId = iB.book2D("heeYG4VsId", "", 600, -300, 300, 600, -300, 300);
   heeZG4VsId = iB.book2D("heeZG4VsId", "", 3000, 320, 350, 3000, 320, 350);
 
   hefXG4VsId = iB.book2D("hefXG4VsId", "", 600, -300, 300, 600, -300, 300);
   hefYG4VsId = iB.book2D("hefYG4VsId", "", 600, -300, 300, 600, -300, 300);
   hefZG4VsId = iB.book2D("hefZG4VsId", "", 6000, 350, 410, 6000, 350, 410);
 
   hebXG4VsId = iB.book2D("hebXG4VsId", "", 600, -300, 300, 600, -300, 300);
   hebYG4VsId = iB.book2D("hebYG4VsId", "", 600, -300, 300, 600, -300, 300);
   hebZG4VsId = iB.book2D("hebZG4VsId", "", 220, 400, 620, 220, 400, 620);
 
   heedzVsZ = iB.book2D("heedzVsZ", "", 600, 320, 350, 100, -1, 1);
   heedyVsY = iB.book2D("heedyVsY", "", 400, -200, 200, 100, -1, 1);
   heedxVsX = iB.book2D("heedxVsX", "", 400, -200, 200, 100, -1, 1);
 
   hefdzVsZ = iB.book2D("hefdzVsZ", "", 1200, 350, 410, 100, -1, 1);
   hefdyVsY = iB.book2D("hefdyVsY", "", 400, -200, 200, 100, -1, 1);
   hefdxVsX = iB.book2D("hefdxVsX", "", 400, -200, 200, 100, -1, 1);
 
   hebdzVsZ = iB.book2D("hebdzVsZ", "", 220, 400, 620, 100, -5, 5);
   hebdyVsY = iB.book2D("hebdyVsY", "", 400, -200, 200, 100, -5, 5);
   hebdxVsX = iB.book2D("hebdxVsX", "", 400, -200, 200, 100, -5, 5);
 
   heedzVsLayer = iB.book2D("heedzVsLayer", "", 100, 0, 100, 100, -1, 1);
   hefdzVsLayer = iB.book2D("hefdzVsLayer", "", 40, 0, 40, 100, -1, 1);
   hebdzVsLayer = iB.book2D("hebdzVsLayer", "", 50, 0, 25, 100, -5, 5);
 
   heedyVsLayer = iB.book2D("heedyVsLayer", "", 100, 0, 100, 100, -1, 1);
   hefdyVsLayer = iB.book2D("hefdyVsLayer", "", 40, 0, 40, 100, -1, 1);
   hebdyVsLayer = iB.book2D("hebdyVsLayer", "", 50, 0, 25, 100, -5, 5);
 
   heedxVsLayer = iB.book2D("heedxVsLayer", "", 100, 0, 100, 100, -1, 1);
   hefdxVsLayer = iB.book2D("hefdxVsLayer", "", 40, 0, 40, 500, -1, 1);
   hebdxVsLayer = iB.book2D("hebdxVsLayer", "", 50, 0, 25, 500, -5, 5.0);
 
   heedX = iB.book1D("heedX", "", 100, -1, 1);
   heedY = iB.book1D("heedY", "", 100, -1, 1);
   heedZ = iB.book1D("heedZ", "", 100, -1, 1);
 
   hefdX = iB.book1D("hefdX", "", 100, -1, 1);
   hefdY = iB.book1D("hefdY", "", 100, -1, 1);
   hefdZ = iB.book1D("hefdZ", "", 100, -1, 1);
 
   hebdX = iB.book1D("hebdX", "", 100, -1, 1);
   hebdY = iB.book1D("hebdY", "", 100, -1, 1);
   hebdZ = iB.book1D("hebdZ", "", 100, -1, 1);
 }

◆ dqmBeginRun()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 84 of file HGCGeometryValidation.cc.

References geometrySource_, geomToken_, edm::EventSetup::getData(), hgcGeometry_, and mps_fire::i.

                                                                                        {
   for (size_t i = 0; i < geometrySource_.size(); i++) {
     hgcGeometry_.emplace_back(&iSetup.getData(geomToken_[i]));
   }
 }

◆ fillDescriptions()

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

static

Definition at line 78 of file HGCGeometryValidation.cc.

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

                                                                                        {
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }