#include <Validation/HGCalValidation/plugins/HGCalHitValidation.cc>

Inheritance diagram for HGCalHitValidation:

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

	~HGCalHitValidation ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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 edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Private Types
typedef std::tuple< float, float, float, float >	HGCHitTuple

Private Attributes
edm::EDGetTokenT < HBHERecHitCollection >	bhRecHitToken_

edm::InputTag	bhSimHitSource

edm::EDGetTokenT< std::vector < PCaloHit > >	bhSimHitToken_

edm::EDGetTokenT < HGCeeRecHitCollection >	eeRecHitToken_

edm::InputTag	eeSimHitSource

edm::EDGetTokenT< std::vector < PCaloHit > >	eeSimHitToken_

edm::EDGetTokenT < HGChefRecHitCollection >	fhRecHitToken_

edm::InputTag	fhSimHitSource

edm::EDGetTokenT< std::vector < PCaloHit > >	fhSimHitToken_

std::vector< std::string >	geometrySource_

const HcalDDDRecConstants *	hcConr_

const HcalDDDSimConstants *	hcCons_

const CaloSubdetectorGeometry *	hcGeometry_

MonitorElement *	hebdEtaVsEta

MonitorElement *	hebdPhiVsPhi

MonitorElement *	hebdzVsZ

MonitorElement *	hebEnRec

MonitorElement *	hebEnSim

MonitorElement *	hebEnSimRec

MonitorElement *	hebRecVsSimX

MonitorElement *	hebRecVsSimY

MonitorElement *	hebRecVsSimZ

MonitorElement *	heedxVsX

MonitorElement *	heedyVsY

MonitorElement *	heedzVsZ

MonitorElement *	heeEnRec

MonitorElement *	heeEnSim

MonitorElement *	heeEnSimRec

MonitorElement *	heeRecVsSimX

MonitorElement *	heeRecVsSimY

MonitorElement *	heeRecVsSimZ

MonitorElement *	hefdxVsX

MonitorElement *	hefdyVsY

MonitorElement *	hefdzVsZ

MonitorElement *	hefEnRec

MonitorElement *	hefEnSim

MonitorElement *	hefEnSimRec

MonitorElement *	hefRecVsSimX

MonitorElement *	hefRecVsSimY

MonitorElement *	hefRecVsSimZ

std::vector< const HGCalDDDConstants * >	hgcCons_

std::vector< const HGCalGeometry * >	hgcGeometry_

std::vector< int >	ietaExcludeBH_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 62 of file HGCalHitValidation.cc.

Member Typedef Documentation

typedef std::tuple<float,float,float,float> HGCalHitValidation::HGCHitTuple

private

Definition at line 76 of file HGCalHitValidation.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 112 of file HGCalHitValidation.cc.

References bhRecHitToken_, bhSimHitToken_, eeRecHitToken_, eeSimHitToken_, fhRecHitToken_, fhSimHitToken_, geometrySource_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ietaExcludeBH_, and relval_2017::k.

                                                                  {
 
   geometrySource_ = cfg.getUntrackedParameter< std::vector<std::string> >("geometrySource");
   eeSimHitToken_  = consumes<std::vector<PCaloHit>>(cfg.getParameter<edm::InputTag>("eeSimHitSource"));
   fhSimHitToken_  = consumes<std::vector<PCaloHit>>(cfg.getParameter<edm::InputTag>("fhSimHitSource"));
   bhSimHitToken_  = consumes<std::vector<PCaloHit>>(cfg.getParameter<edm::InputTag>("bhSimHitSource"));
   eeRecHitToken_  = consumes<HGCeeRecHitCollection>(cfg.getParameter<edm::InputTag>("eeRecHitSource"));
   fhRecHitToken_  = consumes<HGChefRecHitCollection>(cfg.getParameter<edm::InputTag>("fhRecHitSource"));
   bhRecHitToken_  = consumes<HBHERecHitCollection>(cfg.getParameter<edm::InputTag>("bhRecHitSource"));
   ietaExcludeBH_  = cfg.getParameter<std::vector<int> >("ietaExcludeBH");
 
 #ifdef EDM_ML_DEBUG
   edm::LogInfo("HGCalValid") << "Exclude the following " 
                  << ietaExcludeBH_.size() 
                  << " ieta values from BH plots";
   for (unsigned int k=0; k<ietaExcludeBH_.size(); ++k) 
     edm::LogInfo("HGCalValid") << " [" << k << "] " << ietaExcludeBH_[k];
 #endif
 }

HGCalHitValidation::~HGCalHitValidation ( )

Definition at line 132 of file HGCalHitValidation.cc.

132 { }

Member Function Documentation

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

overrideprotectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 237 of file HGCalHitValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), bhRecHitToken_, bhSimHitToken_, HcalDDDSimConstants::cell(), funct::cos(), TauDecayModes::dec, HcalDDDRecConstants::HcalID::depth, HLT_25ns10e33_v2_cff::depth, HGCalTopology::detectorType(), HLT_25ns10e33_v2_cff::distance, eeRecHitToken_, eeSimHitToken_, edm::SortedCollection< T, SORT >::end(), relval_parameters_module::energy, CastorDataFrameFilter_impl::energySum(), HcalDDDRecConstants::HcalID::eta, eta, PV3DBase< T, PVType, FrameType >::eta(), fhRecHitToken_, fhSimHitToken_, MonitorElement::Fill(), spr::find(), edm::Event::getByToken(), CaloSubdetectorGeometry::getGeometry(), HcalDDDRecConstants::getHCID(), CaloCellGeometry::getPosition(), HcalEndcap, hcConr_, hcCons_, hcGeometry_, hebdEtaVsEta, hebdPhiVsPhi, hebdzVsZ, hebEnRec, hebEnSim, hebEnSimRec, hebRecVsSimX, hebRecVsSimY, hebRecVsSimZ, heedxVsX, heedyVsY, heedzVsZ, heeEnRec, heeEnSim, heeEnSimRec, heeRecVsSimX, heeRecVsSimY, heeRecVsSimZ, hefdxVsX, hefdyVsY, hefdzVsZ, hefEnRec, hefEnSim, hefEnSimRec, hefRecVsSimX, hefRecVsSimY, hefRecVsSimZ, hgcCons_, hgcGeometry_, ietaExcludeBH_, edm::HandleBase::isValid(), phi, HcalDDDRecConstants::HcalID::phi, PV3DBase< T, PVType, FrameType >::phi(), edm::Handle< T >::product(), HcalCellType::HcalCell::rz, jetcorrextractor::sign(), funct::sin(), HcalTestNumbering::unpackHcalIndex(), HGCalTestNumbering::unpackHexagonIndex(), Geom::Phi< T >::value(), PV3DBase< T, PVType, FrameType >::x(), create_public_lumi_plots::xy, PV3DBase< T, PVType, FrameType >::y(), z, PV3DBase< T, PVType, FrameType >::z(), and ecaldqm::zside().

                                                                                      {
   std::map<unsigned int, HGCHitTuple> eeHitRefs, fhHitRefs, bhHitRefs;
 
   //declare topology and DDD constants
   const HGCalTopology &heeTopo=hgcGeometry_[0]->topology();
   const HGCalTopology &hefTopo=hgcGeometry_[1]->topology();
 
   //Accesing ee simhits
   edm::Handle<std::vector<PCaloHit>> eeSimHits;
   iEvent.getByToken(eeSimHitToken_, eeSimHits);
 
   if (eeSimHits.isValid()) {
     for (std::vector<PCaloHit>::const_iterator simHit = eeSimHits->begin(); simHit != eeSimHits->end(); ++simHit) {
       int subdet, zside, layer, wafer, celltype, cell;
       HGCalTestNumbering::unpackHexagonIndex(simHit->id(), subdet, zside, layer, wafer, celltype, cell);
       std::pair<float, float> xy = hgcCons_[0]->locateCell(cell,layer,wafer,false);
       float zp = hgcCons_[0]->waferZ(layer,false);
       if (zside < 0) zp = -zp;
       float xp = (zp<0) ? -xy.first/10 : xy.first/10;
       float yp = xy.second/10.0;
 
       //skip this hit if after ganging it is not valid
       std::pair<int,int> recoLayerCell=hgcCons_[0]->simToReco(cell,layer,wafer,heeTopo.detectorType());
       cell  = recoLayerCell.first;
       layer = recoLayerCell.second;
 
       //skip this hit if after ganging it is not valid
       if (layer<0 || cell<0) {
       } else {
     
     //assign the RECO DetId
     HGCalDetId id = HGCalDetId((ForwardSubdetector)(subdet),zside,layer,celltype,wafer,cell);
     float energy = simHit->energy();
 
     float energySum(0);
     if (eeHitRefs.count(id.rawId()) != 0) energySum = std::get<0>(eeHitRefs[id.rawId()]);
     energySum += energy;
     eeHitRefs[id.rawId()] = std::make_tuple(energySum,xp,yp,zp);
     heeEnSim->Fill(energy);
       }
     }
 #ifdef EDM_ML_DEBUG
     for (std::map<unsigned int,HGCHitTuple>::iterator itr=eeHitRefs.begin();
      itr != eeHitRefs.end(); ++itr) {
       int idx = std::distance(eeHitRefs.begin(),itr);
       edm::LogInfo("HGCalValid") << "EEHit[" << idx << "] " << std::hex 
                  << itr->first << std::dec << "; Energy " 
                  << std::get<0>(itr->second) 
                  << "; Position (" << std::get<1>(itr->second) 
                  << ", " << std::get<2>(itr->second) <<", " 
                  << std::get<3>(itr->second) << ")" <<std::endl;
     }
 #endif
   } else {   
     edm::LogWarning("HGCalValid") << "No EE SimHit Found " << std::endl;
   }
 
   //Accesing fh simhits
   edm::Handle<std::vector<PCaloHit>> fhSimHits;
   iEvent.getByToken(fhSimHitToken_, fhSimHits);
   if (fhSimHits.isValid()) {
     for (std::vector<PCaloHit>::const_iterator simHit = fhSimHits->begin(); 
      simHit != fhSimHits->end();++simHit) {
       int subdet, zside, layer, wafer, celltype, cell;
       HGCalTestNumbering::unpackHexagonIndex(simHit->id(), subdet, zside, layer, wafer, celltype, cell);
       std::pair<float, float> xy = hgcCons_[1]->locateCell(cell,layer,wafer,false);
       float zp = hgcCons_[1]->waferZ(layer,false);
       if (zside < 0) zp = -zp;
       float xp = (zp<0) ? -xy.first/10 : xy.first/10;
       float yp = xy.second/10.0;
 
       //skip this hit if after ganging it is not valid
       std::pair<int,int> recoLayerCell = hgcCons_[1]->simToReco(cell,layer,wafer,hefTopo.detectorType());
       cell  = recoLayerCell.first;
       layer = recoLayerCell.second;
       //skip this hit if after ganging it is not valid
       if(layer<0 || cell<0) {
       } else {
     //assign the RECO DetId
     HGCalDetId id = HGCalDetId((ForwardSubdetector)(subdet),zside,layer,celltype,wafer,cell);
 
     float energy = simHit->energy();
     float energySum(0);
     if (fhHitRefs.count(id.rawId()) != 0) energySum = std::get<0>(fhHitRefs[id.rawId()]);
     energySum += energy;
     fhHitRefs[id.rawId()] = std::make_tuple(energySum,xp,yp,zp);
     hefEnSim->Fill(energy);
       }
     }
 #ifdef EDM_ML_DEBUG
     for (std::map<unsigned int,HGCHitTuple>::iterator itr=fhHitRefs.begin();
      itr != fhHitRefs.end(); ++itr) {
       int idx = std::distance(fhHitRefs.begin(),itr);
       edm::LogInfo("HGCalValid") << "FHHit[" << idx << "] " << std::hex 
                  << itr->first << std::dec << "; Energy "
                  << std::get<0>(itr->second) << "; Position (" 
                  << std::get<1>(itr->second) << ", "
                  << std::get<2>(itr->second) <<", " 
                  << std::get<3>(itr->second) << ")" <<std::endl;
     }
 #endif
   } else {
     edm::LogWarning("HGCalValid") << "No FH SimHit Found " << std::endl;
   }
     
   //Accessing bh simhits
   edm::Handle<std::vector<PCaloHit>> bhSimHits;
   iEvent.getByToken(bhSimHitToken_, bhSimHits);
   if (bhSimHits.isValid()) {
     for (std::vector<PCaloHit>::const_iterator simHit = bhSimHits->begin();
      simHit != bhSimHits->end(); ++simHit) {
       int subdet, z, depth, eta, phi, lay;
       HcalTestNumbering::unpackHcalIndex(simHit->id(), subdet, z, depth, eta, phi, lay);
       HcalCellType::HcalCell cell = hcCons_->cell(subdet, z, lay, eta, phi);
 
       double zp  = cell.rz/10; 
 
       if (subdet == static_cast<int>(HcalEndcap)) {
     HcalDDDRecConstants::HcalID idx = hcConr_->getHCID(subdet,eta,phi,lay,depth);
     int sign = (z==0)?(-1):(1);
     zp      *= sign;
     HcalDetId id = HcalDetId(HcalEndcap,sign*idx.eta,idx.phi,idx.depth);  
 
     float energy = simHit->energy();
     float energySum(0);
     if (bhHitRefs.count(id.rawId()) != 0) energySum = std::get<0>(bhHitRefs[id.rawId()]);
     energySum += energy;
     hebEnSim->Fill(energy);
     if (std::find(ietaExcludeBH_.begin(),ietaExcludeBH_.end(),idx.eta) ==
         ietaExcludeBH_.end()) {
       bhHitRefs[id.rawId()] = std::make_tuple(energySum,cell.eta,cell.phi,zp);
 #ifdef EDM_ML_DEBUG
       edm::LogInfo("HGCalValid") << "Accept " << id << std::endl;
     } else {
       edm::LogInfo("HGCalValid") << "Reject " << id << std::endl;
 #endif
     }
       }
     }
 #ifdef EDM_ML_DEBUG
     for (std::map<unsigned int,HGCHitTuple>::iterator itr=bhHitRefs.begin();
      itr != bhHitRefs.end(); ++itr) {
       int idx = std::distance(bhHitRefs.begin(),itr);
       edm::LogInfo("HGCalValid") << "BHHit[" << idx << "] " << std::hex 
                  << itr->first << std::dec << "; Energy " 
                  << std::get<0>(itr->second) << "; Position (" 
                  << std::get<1>(itr->second) << ", "
                  << std::get<2>(itr->second) <<", " 
                  << std::get<3>(itr->second) << ")" <<std::endl;
     }
 #endif
   } else {
     edm::LogWarning("HGCalValid") << "No BH SimHit Found " << std::endl;
   }
 
   //accessing EE Rechit information
   edm::Handle<HGCeeRecHitCollection> eeRecHit;  
   iEvent.getByToken(eeRecHitToken_, eeRecHit);
   if (eeRecHit.isValid()) {
     const HGCeeRecHitCollection* theHits = (eeRecHit.product());
     for (auto it = theHits->begin(); it != theHits->end(); ++it) {
       double energy = it->energy(); 
       heeEnRec->Fill(energy);
       std::map<unsigned int, HGCHitTuple>::const_iterator itr = eeHitRefs.find(it->id().rawId());
       if (itr != eeHitRefs.end()) {
     GlobalPoint xyz = hgcGeometry_[0]->getPosition(it->id());
     heeRecVsSimX->Fill(std::get<1>(itr->second),xyz.x());
     heeRecVsSimY->Fill(std::get<2>(itr->second),xyz.y());
     heeRecVsSimZ->Fill(std::get<3>(itr->second),xyz.z());
     heedxVsX->Fill(std::get<1>(itr->second),(xyz.x()-std::get<1>(itr->second)));
     heedyVsY->Fill(std::get<2>(itr->second),(xyz.y()-std::get<2>(itr->second)));
     heedzVsZ->Fill(std::get<3>(itr->second),(xyz.z()-std::get<3>(itr->second)));
     heeEnSimRec->Fill(std::get<0>(itr->second),energy);
 #ifdef EDM_ML_DEBUG
     edm::LogInfo("HGCalValid") << "EEHit: " << std::hex << it->id().rawId() 
                    << std::dec << " Sim (" 
                    << std::get<0>(itr->second) << ", "
                    << std::get<1>(itr->second) << ", " 
                    << std::get<2>(itr->second) << ", " 
                    << std::get<3>(itr->second) << ") Rec (" 
                    << energy << ", "  << xyz.x() << ", " 
                    << xyz.y() << ", " << xyz.z() << ")\n";
 #endif
       }
     }
   } else {
     edm::LogWarning("HGCalValid") << "No EE RecHit Found " << std::endl;
   }
 
   //accessing FH Rechit information
   edm::Handle<HGChefRecHitCollection> fhRecHit;
   iEvent.getByToken(fhRecHitToken_, fhRecHit);
   if (fhRecHit.isValid()) {
     const HGChefRecHitCollection* theHits = (fhRecHit.product());           
     for (auto it = theHits->begin(); it!=theHits->end(); ++it) {
       double energy = it->energy(); 
       hefEnRec->Fill(energy);
       std::map<unsigned int, HGCHitTuple>::const_iterator itr = fhHitRefs.find(it->id().rawId());
       if (itr != fhHitRefs.end()) {
     GlobalPoint xyz = hgcGeometry_[1]->getPosition(it->id());
 
     hefRecVsSimX->Fill(std::get<1>(itr->second),xyz.x());
         hefRecVsSimY->Fill(std::get<2>(itr->second),xyz.y());
         hefRecVsSimZ->Fill(std::get<3>(itr->second),xyz.z());
         hefdxVsX->Fill(std::get<1>(itr->second),(xyz.x()-std::get<1>(itr->second)));
         hefdyVsY->Fill(std::get<2>(itr->second),(xyz.y()-std::get<2>(itr->second)));
         hefdzVsZ->Fill(std::get<3>(itr->second),(xyz.z()-std::get<3>(itr->second)));
     hefEnSimRec->Fill(std::get<0>(itr->second),energy);
 #ifdef EDM_ML_DEBUG
     edm::LogInfo("HGCalValid") << "FHHit: " << std::hex << it->id().rawId()
                    << std::dec << " Sim (" 
                    << std::get<0>(itr->second) << ", "
                    << std::get<1>(itr->second) << ", " 
                    << std::get<2>(itr->second) << ", " 
                    << std::get<3>(itr->second) << ") Rec (" 
                    << energy << "," << xyz.x() << ", " 
                    << xyz.y() << ", " << xyz.z() << ")\n";
 #endif
       }
     }
   } else {
     edm::LogWarning("HGCalValid") << "No FH RecHit Found " << std::endl;
   }
 
 
   //accessing BH Rechit information
   edm::Handle<HBHERecHitCollection> bhRecHit;
   iEvent.getByToken(bhRecHitToken_, bhRecHit);
   if (bhRecHit.isValid()) {
     const HBHERecHitCollection* theHits = (bhRecHit.product());
     
     for (auto it = theHits->begin(); it!=theHits->end(); ++it) {
       DetId id = it->id();
       if (id.subdetId() == (int)(HcalEndcap)) {
     double energy = it->energy();
     hebEnRec->Fill(energy);
     GlobalPoint xyz = hcGeometry_->getGeometry(id)->getPosition();
 
     std::map<unsigned int, HGCHitTuple>::const_iterator itr = bhHitRefs.find(id.rawId());
     if (itr != bhHitRefs.end()) {
       float ang3 = xyz.phi().value(); // returns the phi in radians
       double fac = sinh(std::get<1>(itr->second));
       double pT  = std::get<3>(itr->second) / fac;
       double xp = pT * cos(std::get<2>(itr->second));
       double yp = pT * sin(std::get<2>(itr->second));
       hebRecVsSimX->Fill(xp,xyz.x());
       hebRecVsSimY->Fill(yp,xyz.y());
       hebRecVsSimZ->Fill(std::get<3>(itr->second),xyz.z());
       hebdEtaVsEta->Fill(std::get<1>(itr->second),(xyz.eta()-std::get<1>(itr->second)));
       hebdPhiVsPhi->Fill(std::get<2>(itr->second),(ang3-std::get<2>(itr->second)));
       hebdzVsZ->Fill(std::get<3>(itr->second),(xyz.z()-std::get<3>(itr->second)));
       hebEnSimRec->Fill(std::get<0>(itr->second),energy);
 
 #ifdef EDM_ML_DEBUG
       edm::LogInfo("HGCalValid") << "BHHit: " << std::hex << id.rawId()
                      << std::dec << " Sim (" 
                      << std::get<0>(itr->second) << ", "
                      << std::get<1>(itr->second) << ", " 
                      << std::get<2>(itr->second) << ", " 
                      << std::get<3>(itr->second) << ") Rec ("
                      << energy << ", "  << xyz.x() << ", " 
                      << xyz.y() << ", " << xyz.z() << ")\n";
 #endif
     }
       }
     }
   } else {
     edm::LogWarning("HGCalValid") << "No BH RecHit Found " << std::endl;
   }
   
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 142 of file HGCalHitValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), hebdEtaVsEta, hebdPhiVsPhi, hebdzVsZ, hebEnRec, hebEnSim, hebEnSimRec, hebRecVsSimX, hebRecVsSimY, hebRecVsSimZ, heedxVsX, heedyVsY, heedzVsZ, heeEnRec, heeEnSim, heeEnSimRec, heeRecVsSimX, heeRecVsSimY, heeRecVsSimZ, hefdxVsX, hefdyVsY, hefdzVsZ, hefEnRec, hefEnSim, hefEnSimRec, hefRecVsSimX, hefRecVsSimY, hefRecVsSimZ, M_PI, and DQMStore::IBooker::setCurrentFolder().

                                           {
 
   iB.setCurrentFolder("HGCalSimHitsV/HitValidation");
 
   //initiating histograms
   heedzVsZ     = iB.book2D("heedzVsZ","",7200,-360,360,100,-0.1,0.1);
   heedyVsY     = iB.book2D("heedyVsY","",400,-200,200,100,-0.02,0.02);
   heedxVsX     = iB.book2D("heedxVsX","",400,-200,200,100,-0.02,0.02);
   heeRecVsSimZ = iB.book2D("heeRecVsSimZ","",7200,-360,360,7200,-360,360);
   heeRecVsSimY = iB.book2D("heeRecVsSimY","",400,-200,200,400,-200,200);
   heeRecVsSimX = iB.book2D("heeRecVsSimX","",400,-200,200,400,-200,200);
 
   hefdzVsZ     = iB.book2D("hefdzVsZ","",8200,-410,410,100,-0.1,0.1);
   hefdyVsY     = iB.book2D("hefdyVsY","",400,-200,200,100,-0.02,0.02);
   hefdxVsX     = iB.book2D("hefdxVsX","",400,-200,200,100,-0.02,0.02);
   hefRecVsSimZ = iB.book2D("hefRecVsSimZ","",8200,-410,410,8200,-410,410);
   hefRecVsSimY = iB.book2D("hefRecVsSimY","",400,-200,200,400,-200,200);
   hefRecVsSimX = iB.book2D("hefRecVsSimX","",400,-200,200,400,-200,200);
 
   hebdzVsZ     = iB.book2D("hebdzVsZ","",1080,-540,540,100,-1.0,1.0);
   hebdPhiVsPhi = iB.book2D("hebdPhiVsPhi","",M_PI*100,-0.5,M_PI+0.5,200,-0.2,0.2);
   hebdEtaVsEta = iB.book2D("hebdEtaVsEta","",1000,-5,5,200,-0.1,0.1);
   hebRecVsSimZ = iB.book2D("hebRecVsSimZ","",1080,-540,540,1080,-540,540);
   hebRecVsSimY = iB.book2D("hebRecVsSimY","",400,-200,200,400,-200,200);
   hebRecVsSimX = iB.book2D("hebRecVsSimX","",400,-200,200,400,-200,200);
 
   heeEnRec     = iB.book1D("heeEnRec","",1000,0,10);
   heeEnSim     = iB.book1D("heeEnSim","",1000,0,0.01);
   heeEnSimRec  = iB.book2D("heeEnSimRec","",200,0,0.002,200,0,0.2);
 
   hefEnRec     = iB.book1D("hefEnRec","",1000,0,10);
   hefEnSim     = iB.book1D("hefEnSim","",1000,0,0.01);
   hefEnSimRec  = iB.book2D("hefEnSimRec","",200,0,0.001,200,0,0.5);
 
   hebEnRec     = iB.book1D("hebEnRec","",1000,0,15);
   hebEnSim     = iB.book1D("hebEnSim","",1000,0,0.01);
   hebEnSimRec  = iB.book2D("hebEnSimRec","",200,0,0.02,200,0,4);
 
 }

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 184 of file HGCalHitValidation.cc.

References spr::find(), geometrySource_, edm::EventSetup::get(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, hcConr_, hcCons_, hcGeometry_, hgcCons_, hgcGeometry_, i, edm::ESHandleBase::isValid(), and edm::ESHandle< class >::product().

                                                   {
   //initiating hgc Geometry
   for (size_t i=0; i<geometrySource_.size(); i++) {
     if (geometrySource_[i].find("Hcal") != std::string::npos) {
       edm::ESHandle<HcalDDDSimConstants> pHSNDC;
       iSetup.get<HcalSimNumberingRecord>().get(pHSNDC);
       if (pHSNDC.isValid()) {
         hcCons_ = pHSNDC.product();
         hgcCons_.push_back(0);
       } else {
         edm::LogWarning("HGCalValid") << "Cannot initiate HcalDDDSimConstants: "
                                       << geometrySource_[i] << std::endl;
       }
       edm::ESHandle<HcalDDDRecConstants> pHRNDC;
       iSetup.get<HcalRecNumberingRecord>().get(pHRNDC);
       if (pHRNDC.isValid()) {
         hcConr_ = pHRNDC.product();
       } else {
         edm::LogWarning("HGCalValid") << "Cannot initiate HcalDDDRecConstants: "
                                       << geometrySource_[i] << std::endl;
       }
       edm::ESHandle<CaloGeometry> caloG;
       iSetup.get<CaloGeometryRecord>().get(caloG);
       if (caloG.isValid()) {
     const CaloGeometry* geo = caloG.product();
     hcGeometry_ = geo->getSubdetectorGeometry(DetId::Hcal,HcalBarrel);
     hgcGeometry_.push_back(0);
       } else {
         edm::LogWarning("HGCalValid") << "Cannot initiate HcalGeometry for "
                                       << geometrySource_[i] << std::endl;
       }
    } else {
       edm::ESHandle<HGCalDDDConstants> hgcCons;
       iSetup.get<IdealGeometryRecord>().get(geometrySource_[i],hgcCons);
       if (hgcCons.isValid()) {
         hgcCons_.push_back(hgcCons.product());
       } else {
         edm::LogWarning("HGCalValid") << "Cannot initiate HGCalDDDConstants for "
                                       << geometrySource_[i] << std::endl;
       }
       edm::ESHandle<HGCalGeometry> hgcGeom;
       iSetup.get<IdealGeometryRecord>().get(geometrySource_[i],hgcGeom);    
       if(hgcGeom.isValid()) {
     hgcGeometry_.push_back(hgcGeom.product());  
       } else {
     edm::LogWarning("HGCalValid") << "Cannot initiate HGCalGeometry for "
                       << geometrySource_[i] << std::endl;
       }
     }
   }
 }

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

static

Definition at line 134 of file HGCalHitValidation.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                                     {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }