Inheritance diagram for HcalHitValidation:

Public Member Functions
	HcalHitValidation (const edm::ParameterSet &ps)

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

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

	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

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

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

void	analyzeJets (const edm::Handle< PHcalValidInfoJets > &)

void	analyzeLayer (const edm::Handle< PHcalValidInfoLayer > &)

void	analyzeNxN (const edm::Handle< PHcalValidInfoNxN > &)

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

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Attributes
const bool	checkHit_

const bool	checkJet_

const bool	checkLay_

const bool	checkNxN_

const std::string	g4Label

const std::string	hcalHits

const std::string	jetsInfo

const std::string	layerInfo

MonitorElement *	meAllNHit_

MonitorElement *	meBadDetHit_

MonitorElement *	meBadIdHit_

MonitorElement *	meBadSubHit_

MonitorElement *	meDepHlay_

MonitorElement *	meDepLay_ [nDepthsMAX]

MonitorElement *	meDepthHit_

MonitorElement *	meDetaJet_

MonitorElement *	meDetectHit_

MonitorElement *	meDphiJet_

MonitorElement *	meDrJet_

MonitorElement *	meEcalJet_

MonitorElement *	meEcalNxN_

MonitorElement *	meEcalRNxN_

MonitorElement *	meEcHcJet_

MonitorElement *	meEFibLLay_

MonitorElement *	meEFibSLay_

MonitorElement *	meEHBHELay_

MonitorElement *	meEHFEmLay_

MonitorElement *	meEHFHdLay_

MonitorElement *	meEHOLay_

MonitorElement *	meEiNxN_

MonitorElement *	meEJet_

MonitorElement *	meEneDLay_

MonitorElement *	meEneHLay_

MonitorElement *	meEneJet_

MonitorElement *	meEneLay_ [nLayersMAX]

MonitorElement *	meEneLLay_

MonitorElement *	meEnergyHit_

MonitorElement *	meEtaHit_

MonitorElement *	meEtaHLay_

MonitorElement *	meEtaJet_

MonitorElement *	meEtaPhi_

MonitorElement *	meEtotJet_

MonitorElement *	meEtotLay_

MonitorElement *	meEtotNxN_

MonitorElement *	meEtotRNxN_

MonitorElement *	meHBDepHit_

MonitorElement *	meHBEneHit_

MonitorElement *	meHBEtaHit_

MonitorElement *	meHBNHit_

MonitorElement *	meHBPhiHit_

MonitorElement *	meHBTimHit_

MonitorElement *	meHcalJet_

MonitorElement *	meHcalNxN_

MonitorElement *	meHcalRNxN_

MonitorElement *	meHEDepHit_

MonitorElement *	meHEEneHit_

MonitorElement *	meHEEtaHit_

MonitorElement *	meHENHit_

MonitorElement *	meHEPhiHit_

MonitorElement *	meHETimHit_

MonitorElement *	meHFDepHit_

MonitorElement *	meHFEneHit_

MonitorElement *	meHFEtaHit_

MonitorElement *	meHFNHit_

MonitorElement *	meHFPhiHit_

MonitorElement *	meHFTimHit_

MonitorElement *	meHitELay_

MonitorElement *	meHitHLay_

MonitorElement *	meHitTLay_

MonitorElement *	meHODepHit_

MonitorElement *	meHOEneHit_

MonitorElement *	meHOEtaHit_

MonitorElement *	meHoJet_

MonitorElement *	meHONHit_

MonitorElement *	meHoNxN_

MonitorElement *	meHOPhiHit_

MonitorElement *	meHoRNxN_

MonitorElement *	meHOTimHit_

MonitorElement *	meLayerLay_

MonitorElement *	meLngLay_

MonitorElement *	meMassJet_

MonitorElement *	mePhiHit_

MonitorElement *	mePhiHLay_

MonitorElement *	mePhiJet_

MonitorElement *	mePMTDepHit_

MonitorElement *	mePMTEn1Hit_

MonitorElement *	mePMTEn2Hit_

MonitorElement *	mePMTEtaHit_

MonitorElement *	mePMTHit_

MonitorElement *	mePMTPhiHit_

MonitorElement *	mePMTTimHit_

MonitorElement *	meRJet_

MonitorElement *	meSubdetHit_

MonitorElement *	meTimeHit_

MonitorElement *	meTimeWHit_

MonitorElement *	meTimHLay_

MonitorElement *	meTimWLay_

MonitorElement *	meTiNxN_

MonitorElement *	meTJet_

MonitorElement *	meTrNxN_

const std::string	nxNInfo

const std::string	outFile_

const bool	scheme_

const edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hh_

const edm::EDGetTokenT< PHcalValidInfoJets >	tok_iJ_

const edm::EDGetTokenT< PHcalValidInfoLayer >	tok_iL_

const edm::EDGetTokenT< PHcalValidInfoNxN >	tok_iN_

const bool	verbose_

Static Private Attributes
static const int	nDepthsMAX = 5

static const int	nLayersMAX = 20

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

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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 26 of file HcalHitValidation.cc.

Constructor & Destructor Documentation

◆ HcalHitValidation()

HcalHitValidation::HcalHitValidation ( const edm::ParameterSet & ps )

Definition at line 81 of file HcalHitValidation.cc.

References checkHit_, checkJet_, checkLay_, checkNxN_, g4Label, hcalHits, jetsInfo, layerInfo, nxNInfo, outFile_, and scheme_.

     : g4Label(ps.getUntrackedParameter<std::string>("moduleLabel", "g4SimHits")),
       hcalHits(ps.getUntrackedParameter<std::string>("HitCollection", "HcalHits")),
       layerInfo(ps.getUntrackedParameter<std::string>("LayerInfo", "PHcalValidInfoLayer")),
       nxNInfo(ps.getUntrackedParameter<std::string>("NxNInfo", "PHcalValidInfoNxN")),
       jetsInfo(ps.getUntrackedParameter<std::string>("JetsInfo", "PHcalValidInfoJets")),
       outFile_(ps.getUntrackedParameter<std::string>("outputFile", "hcValid.root")),
       verbose_(ps.getUntrackedParameter<bool>("Verbose", false)),
       scheme_(ps.getUntrackedParameter<bool>("TestNumbering", true)),
       checkHit_(ps.getUntrackedParameter<bool>("CheckHits", true)),
       checkLay_(ps.getUntrackedParameter<bool>("CheckLayer", true)),
       checkNxN_(ps.getUntrackedParameter<bool>("CheckNxN", true)),
       checkJet_(ps.getUntrackedParameter<bool>("CheckJets", true)),
       tok_hh_(consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label, hcalHits))),
       tok_iL_(consumes<PHcalValidInfoLayer>(edm::InputTag(g4Label, layerInfo))),
       tok_iN_(consumes<PHcalValidInfoNxN>(edm::InputTag(g4Label, nxNInfo))),
       tok_iJ_(consumes<PHcalValidInfoJets>(edm::InputTag(g4Label, jetsInfo))) {
   edm::LogVerbatim("HcalHitValid") << "Module Label: " << g4Label << "   Hits: " << hcalHits << " / " << checkHit_
                                    << "   LayerInfo: " << layerInfo << " / " << checkLay_ << "  NxNInfo: " << nxNInfo
                                    << " / " << checkNxN_ << "  jetsInfo: " << jetsInfo << " / " << checkJet_
                                    << "   Output: " << outFile_ << "   Usage of TestNumberingScheme " << scheme_;
 }

◆ ~HcalHitValidation()

HcalHitValidation::~HcalHitValidation ( )

overridedefault

Member Function Documentation

◆ analyze()

void HcalHitValidation::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 232 of file HcalHitValidation.cc.

References analyzeHits(), analyzeJets(), analyzeLayer(), analyzeNxN(), checkHit_, checkJet_, checkLay_, checkNxN_, MillePedeFileConverter_cfg::e, pixelTopology::getLayer(), edm::HandleBase::isValid(), tok_hh_, tok_iJ_, tok_iL_, and tok_iN_.

                                                                         {
   edm::LogVerbatim("HcalHitValid") << "Run = " << e.id().run() << " Event = " << e.id().event();
 
   bool getHits = false;
   if (checkHit_) {
     const edm::Handle<edm::PCaloHitContainer> &hitsHcal = e.getHandle(tok_hh_);
     if (hitsHcal.isValid())
       getHits = true;
     if (getHits) {
       std::vector<PCaloHit> caloHits;
       caloHits.insert(caloHits.end(), hitsHcal->begin(), hitsHcal->end());
       edm::LogVerbatim("HcalHitValid") << "HcalValidation: Hit buffer " << caloHits.size();
       analyzeHits(caloHits);
     }
   }
 
   bool getLayer = false;
   if (checkLay_) {
     const edm::Handle<PHcalValidInfoLayer> &infoLayer = e.getHandle(tok_iL_);
     if (infoLayer.isValid()) {
       getLayer = true;
       analyzeLayer(infoLayer);
     }
   }
 
   bool getNxN = false;
   if (checkNxN_) {
     const edm::Handle<PHcalValidInfoNxN> &infoNxN = e.getHandle(tok_iN_);
     if (infoNxN.isValid()) {
       getNxN = true;
       analyzeNxN(infoNxN);
     }
   }
 
   bool getJets = false;
   if (checkJet_) {
     const edm::Handle<PHcalValidInfoJets> &infoJets = e.getHandle(tok_iJ_);
     if (infoJets.isValid()) {
       getJets = true;
       analyzeJets(infoJets);
     }
   }
 
   edm::LogVerbatim("HcalHitValid") << "HcalValidation: Input flags Hits " << getHits << ", Layer " << getLayer
                                    << ", NxN " << getNxN << ", Jets " << getJets;
 }

◆ analyzeHits()

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

protected

Definition at line 279 of file HcalHitValidation.cc.

References TauDecayModes::dec, hcalRecHitTable_cff::depth, hcalRecHitTable_cff::energy, PVValHelper::eta, dqm::impl::MonitorElement::Fill(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hfClusterShapes_cfi::hits, mps_fire::i, meAllNHit_, meBadDetHit_, meBadIdHit_, meBadSubHit_, meDepthHit_, meDetectHit_, meEnergyHit_, meEtaHit_, meHBDepHit_, meHBEneHit_, meHBEtaHit_, meHBNHit_, meHBPhiHit_, meHBTimHit_, meHEDepHit_, meHEEneHit_, meHEEtaHit_, meHENHit_, meHEPhiHit_, meHETimHit_, meHFDepHit_, meHFEneHit_, meHFEtaHit_, meHFNHit_, meHFPhiHit_, meHFTimHit_, meHODepHit_, meHOEneHit_, meHOEtaHit_, meHONHit_, meHOPhiHit_, meHOTimHit_, mePhiHit_, mePMTDepHit_, mePMTEn1Hit_, mePMTEn2Hit_, mePMTEtaHit_, mePMTHit_, mePMTPhiHit_, mePMTTimHit_, meSubdetHit_, meTimeHit_, meTimeWHit_, phi, scheme_, hcalRecHitTable_cff::time, and ecaldqm::zside().

Referenced by analyze().

                                                              {
   int nHit = hits.size();
   int nHB = 0, nHE = 0, nHO = 0, nHF = 0, nPMT = 0, nBad1 = 0, nBad2 = 0, nBad = 0;
   for (int i = 0; i < nHit; i++) {
     double energy = hits[i].energy();
     double time = hits[i].time();
     unsigned int id_ = hits[i].id();
     int det, subdet, depth, eta, phi;
     if (scheme_) {
       det = 4;
       subdet = (id_ >> 28) & 15;
       depth = (id_ >> 26) & 3;
       depth++;
       int zside = ((id_ & 0x100000) ? (1) : (-1));
       eta = zside * ((id_ >> 10) & 1023);
       phi = (id_ & 1023);
     } else {
       HcalDetId id = HcalDetId(id_);
       det = id.det();
       subdet = id.subdet();
       depth = id.depth();
       eta = id.ieta();
       phi = id.iphi();
     }
     uint16_t depth_ = hits[i].depth();
     double energyEM = hits[i].energyEM();
     double energyHad = hits[i].energyHad();
     edm::LogVerbatim("HcalHitValid") << "Hit[" << i << "] ID " << std::hex << id_ << std::dec << " Det " << det
                                      << " Sub " << subdet << " depth " << depth << " " << depth_ << " Eta " << eta
                                      << " Phi " << phi << " E " << energy << "(EM " << energyEM << ", Had " << energyHad
                                      << ") time " << time;
     if (det == 4) {  // Check DetId.h
       if (subdet == static_cast<int>(HcalBarrel)) {
         nHB++;
       } else if (subdet == static_cast<int>(HcalEndcap)) {
         nHE++;
       } else if (subdet == static_cast<int>(HcalOuter)) {
         nHO++;
       } else if (subdet == static_cast<int>(HcalForward)) {
         if (depth_ == 0)
           nHF++;
         else
           nPMT++;
       } else {
         nBad++;
         nBad2++;
       }
     } else {
       nBad++;
       nBad1++;
     }
 
     meDetectHit_->Fill(double(det));
     if (det == 4 && depth_ == 0) {
       meSubdetHit_->Fill(double(subdet));
       meDepthHit_->Fill(double(depth));
       meEtaHit_->Fill(double(eta));
       mePhiHit_->Fill(double(phi));
       meEnergyHit_->Fill(energy);
       meTimeHit_->Fill(time);
       meTimeWHit_->Fill(double(time), energy);
       if (subdet == static_cast<int>(HcalBarrel)) {
         meHBDepHit_->Fill(double(depth));
         meHBEtaHit_->Fill(double(eta));
         meHBPhiHit_->Fill(double(phi));
         meHBEneHit_->Fill(energy);
         meHBTimHit_->Fill(time);
       } else if (subdet == static_cast<int>(HcalEndcap)) {
         meHEDepHit_->Fill(double(depth));
         meHEEtaHit_->Fill(double(eta));
         meHEPhiHit_->Fill(double(phi));
         meHEEneHit_->Fill(energy);
         meHETimHit_->Fill(time);
       } else if (subdet == static_cast<int>(HcalOuter)) {
         meHODepHit_->Fill(double(depth));
         meHOEtaHit_->Fill(double(eta));
         meHOPhiHit_->Fill(double(phi));
         meHOEneHit_->Fill(energy);
         meHOTimHit_->Fill(time);
       } else if (subdet == static_cast<int>(HcalForward)) {
         meHFDepHit_->Fill(double(depth));
         meHFEtaHit_->Fill(double(eta));
         meHFPhiHit_->Fill(double(phi));
         meHFEneHit_->Fill(energy);
         meHFTimHit_->Fill(time);
       }
     } else if (det == 0 && subdet == static_cast<int>(HcalForward)) {
       mePMTDepHit_->Fill(double(depth));
       mePMTEtaHit_->Fill(double(eta));
       mePMTPhiHit_->Fill(double(phi));
       mePMTEn1Hit_->Fill(energyEM);
       mePMTEn2Hit_->Fill(energyHad);
       mePMTTimHit_->Fill(time);
     }
   }
   meAllNHit_->Fill(double(nHit));
   meBadDetHit_->Fill(double(nBad1));
   meBadSubHit_->Fill(double(nBad2));
   meBadIdHit_->Fill(double(nBad));
   meHBNHit_->Fill(double(nHB));
   meHENHit_->Fill(double(nHE));
   meHONHit_->Fill(double(nHO));
   meHFNHit_->Fill(double(nHF));
   mePMTHit_->Fill(double(nPMT));
 
   edm::LogVerbatim("HcalHitValid") << "HcalHitValidation::analyzeHits: HB " << nHB << " HE " << nHE << " HO " << nHO
                                    << " HF " << nHF << " PMT " << nPMT << " Bad " << nBad << " All " << nHit;
 }

◆ analyzeJets()

void HcalHitValidation::analyzeJets ( const edm::Handle< PHcalValidInfoJets > & infoJets )

protected

Definition at line 503 of file HcalHitValidation.cc.

Referenced by analyze().

                                                                                  {
   // -- Leading Jet
   int nJetHits = infoJets->njethit();
 
   std::vector<float> rJetHits = infoJets->jethitr();
   std::vector<float> tJetHits = infoJets->jethitt();
   std::vector<float> eJetHits = infoJets->jethite();
 
   double ecalJet = infoJets->ecaljet();
   double hcalJet = infoJets->hcaljet();
   double hoJet = infoJets->hojet();
   double etotJet = infoJets->etotjet();
 
   double detaJet = infoJets->detajet();
   double dphiJet = infoJets->dphijet();
   double drJet = infoJets->drjet();
   double dijetM = infoJets->dijetm();
 
   for (int j = 0; j < nJetHits; j++) {
     meRJet_->Fill(rJetHits[j]);
     meTJet_->Fill(tJetHits[j]);
     meEJet_->Fill(eJetHits[j]);
   }
 
   meEcalJet_->Fill(ecalJet);
   meHcalJet_->Fill(hcalJet);
   meHoJet_->Fill(hoJet);
   meEtotJet_->Fill(etotJet);
   meEcHcJet_->Fill(ecalJet, hcalJet);
 
   meDetaJet_->Fill(detaJet);
   meDphiJet_->Fill(dphiJet);
   meDrJet_->Fill(drJet);
   meMassJet_->Fill(dijetM);
 
   // All Jets
   int nJets = infoJets->njet();
   std::vector<float> jetE = infoJets->jete();
   std::vector<float> jetEta = infoJets->jeteta();
   std::vector<float> jetPhi = infoJets->jetphi();
 
   for (int j = 0; j < nJets; j++) {
     meEneJet_->Fill(jetE[j]);
     meEtaJet_->Fill(jetEta[j]);
     mePhiJet_->Fill(jetPhi[j]);
   }
   edm::LogVerbatim("HcalHitValid") << "HcalHitValidation::analyzeJets: " << nJets << " jets with " << nJetHits
                                    << " hits in the leading jet\n"
                                    << "   d(Eta) = " << detaJet << "  d(Phi) = " << dphiJet << "  d(R) = " << drJet
                                    << "  diJet Mass = " << dijetM;
 }

◆ analyzeLayer()

void HcalHitValidation::analyzeLayer ( const edm::Handle< PHcalValidInfoLayer > & infoLayer )

protected

Definition at line 388 of file HcalHitValidation.cc.

Referenced by analyze().

                                                                                     {
   // CaloHits from PHcalValidInfoLayer
   int nHits = infoLayer->nHit();
   std::vector<float> idHits = infoLayer->idHit();
   std::vector<float> phiHits = infoLayer->phiHit();
   std::vector<float> etaHits = infoLayer->etaHit();
   std::vector<float> layerHits = infoLayer->layerHit();
   std::vector<float> eHits = infoLayer->eHit();
   std::vector<float> tHits = infoLayer->tHit();
 
   int ne = 0, nh = 0;
   for (int j = 0; j < nHits; j++) {
     int layer = (int)(layerHits[j]) - 1;
     int id = (int)(idHits[j]);
 
     if (id >= 10) {
       ne++;
     } else {
       nh++;
     }
 
     edm::LogVerbatim("HcalHitValid") << "HcalHitValidation::analyzeLayer:Hit subdet = " << id << "  lay = " << layer;
 
     meLayerLay_->Fill(double(layer));
     meEtaHLay_->Fill(etaHits[j]);
     mePhiHLay_->Fill(phiHits[j]);
     meEneHLay_->Fill(eHits[j]);
     meDepHlay_->Fill(idHits[j]);
     meTimHLay_->Fill(tHits[j]);
     meTimWLay_->Fill(tHits[j], eHits[j]);
     if (id < 6)  // HCAL only. Depth is needed, not layer !!!
     {
       meEtaPhi_->Fill(etaHits[j], phiHits[j]);
     }
   }
 
   meHitELay_->Fill(double(ne));
   meHitHLay_->Fill(double(nh));
   meHitTLay_->Fill(double(nHits));
 
   // Layers and depths PHcalValidInfoLayer
   std::vector<float> eLayer = infoLayer->elayer();
   std::vector<float> eDepth = infoLayer->edepth();
   float eTot = 0.;
 
   for (int j = 0; j < nLayersMAX; j++) {
     eTot += eLayer[j];
     meEneLLay_->Fill(eLayer[j]);
     meEneLay_[j]->Fill(eLayer[j]);
     meLngLay_->Fill((double)(j), eLayer[j]);  // HCAL SimHits only
   }
   for (int j = 0; j < nDepthsMAX; j++) {
     meEneDLay_->Fill(eDepth[j]);
     meDepLay_[j]->Fill(eDepth[j]);
   }
   meEtotLay_->Fill(eTot);
 
   // The rest  PHcalValidInfoLayer
   double eHO = infoLayer->eho();
   double eHBHE = infoLayer->ehbhe();
   double elongHF = infoLayer->elonghf();
   double eshortHF = infoLayer->eshorthf();
   double eEcalHF = infoLayer->eecalhf();
   double eHcalHF = infoLayer->ehcalhf();
 
   meEHOLay_->Fill(eHO);
   meEHBHELay_->Fill(eHBHE);
   meEFibLLay_->Fill(elongHF);
   meEFibSLay_->Fill(eshortHF);
   meEHFEmLay_->Fill(eEcalHF);
   meEHFHdLay_->Fill(eHcalHF);
 
   edm::LogVerbatim("HcalHitValid") << "HcalHitValidation::analyzeLayer: eHO " << eHO << "  eHBHE = " << eHBHE
                                    << " elongHF = " << elongHF << " eshortHF = " << eshortHF
                                    << "  eEcalHF = " << eEcalHF << "  eHcalHF = " << eHcalHF;
 }

◆ analyzeNxN()

void HcalHitValidation::analyzeNxN ( const edm::Handle< PHcalValidInfoNxN > & infoNxN )

protected

Definition at line 465 of file HcalHitValidation.cc.

Referenced by analyze().

                                                                               {
   // NxN quantities
   double ecalNxNr = infoNxN->ecalnxnr();
   double hcalNxNr = infoNxN->hcalnxnr();
   double hoNxNr = infoNxN->honxnr();
   double etotNxNr = infoNxN->etotnxnr();
 
   double ecalNxN = infoNxN->ecalnxn();
   double hcalNxN = infoNxN->hcalnxn();
   double hoNxN = infoNxN->honxn();
   double etotNxN = infoNxN->etotnxn();
 
   meEcalRNxN_->Fill(ecalNxNr);
   meHcalRNxN_->Fill(hcalNxNr);
   meHoRNxN_->Fill(hoNxNr);
   meEtotRNxN_->Fill(etotNxNr);
 
   meEcalNxN_->Fill(ecalNxN);
   meHcalNxN_->Fill(hcalNxN);
   meHoNxN_->Fill(hoNxN);
   meEtotNxN_->Fill(etotNxN);
 
   int nIxI = infoNxN->nnxn();
   std::vector<float> idIxI = infoNxN->idnxn();
   std::vector<float> eIxI = infoNxN->enxn();
   std::vector<float> tIxI = infoNxN->tnxn();
 
   for (int j = 0; j < nIxI; j++)  // NB !!! j < nIxI
   {
     meEiNxN_->Fill(eIxI[j]);
     meTiNxN_->Fill(tIxI[j]);
     meTrNxN_->Fill(idIxI[j], eIxI[j]);  // transverse profile
   }
 
   edm::LogVerbatim("HcalHitValid") << "HcalHitValidation::analyzeNxN: " << nIxI
                                    << " hits in NxN analysis; Total Energy " << etotNxN << "/" << etotNxNr;
 }

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 104 of file HcalHitValidation.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), checkHit_, checkJet_, checkLay_, checkNxN_, mps_fire::i, meAllNHit_, meBadDetHit_, meBadIdHit_, meBadSubHit_, meDepHlay_, meDepLay_, meDepthHit_, meDetaJet_, meDetectHit_, meDphiJet_, meDrJet_, meEcalJet_, meEcalNxN_, meEcalRNxN_, meEcHcJet_, meEFibLLay_, meEFibSLay_, meEHBHELay_, meEHFEmLay_, meEHFHdLay_, meEHOLay_, meEiNxN_, meEJet_, meEneDLay_, meEneHLay_, meEneJet_, meEneLay_, meEneLLay_, meEnergyHit_, meEtaHit_, meEtaHLay_, meEtaJet_, meEtaPhi_, meEtotJet_, meEtotLay_, meEtotNxN_, meEtotRNxN_, meHBDepHit_, meHBEneHit_, meHBEtaHit_, meHBNHit_, meHBPhiHit_, meHBTimHit_, meHcalJet_, meHcalNxN_, meHcalRNxN_, meHEDepHit_, meHEEneHit_, meHEEtaHit_, meHENHit_, meHEPhiHit_, meHETimHit_, meHFDepHit_, meHFEneHit_, meHFEtaHit_, meHFNHit_, meHFPhiHit_, meHFTimHit_, meHitELay_, meHitHLay_, meHitTLay_, meHODepHit_, meHOEneHit_, meHOEtaHit_, meHoJet_, meHONHit_, meHoNxN_, meHOPhiHit_, meHoRNxN_, meHOTimHit_, meLayerLay_, meLngLay_, meMassJet_, mePhiHit_, mePhiHLay_, mePhiJet_, mePMTDepHit_, mePMTEn1Hit_, mePMTEn2Hit_, mePMTEtaHit_, mePMTHit_, mePMTPhiHit_, mePMTTimHit_, meRJet_, meSubdetHit_, meTimeHit_, meTimeWHit_, meTimHLay_, meTimWLay_, meTiNxN_, meTJet_, meTrNxN_, Skims_PA_cff::name, nDepthsMAX, nLayersMAX, groupFilesInBlocks::nn, dqm::implementation::NavigatorBase::setCurrentFolder(), and runGCPTkAlMap::title.

                                                                                                       {
   ibooker.setCurrentFolder("HcalHitValidation");
 
   char title[60], name[20];
   double my_pi = 3.1415926;
   // Histograms for Hits
   if (checkHit_) {
     meAllNHit_ = ibooker.book1D("Hit01", "Number of Hits in HCal", 1000, 0., 5000.);
     meBadDetHit_ = ibooker.book1D("Hit02", "Hits with wrong Det", 100, 0., 100.);
     meBadSubHit_ = ibooker.book1D("Hit03", "Hits with wrong Subdet", 100, 0., 100.);
     meBadIdHit_ = ibooker.book1D("Hit04", "Hits with wrong ID", 100, 0., 100.);
     meHBNHit_ = ibooker.book1D("Hit05", "Number of Hits in HB", 1000, 0., 5000.);
     meHENHit_ = ibooker.book1D("Hit06", "Number of Hits in HE", 1000, 0., 5000.);
     meHONHit_ = ibooker.book1D("Hit07", "Number of Hits in HO", 1000, 0., 5000.);
     meHFNHit_ = ibooker.book1D("Hit08", "Number of Hits in HF", 1000, 0., 5000.);
     meDetectHit_ = ibooker.book1D("Hit09", "Detector ID", 50, 0., 50.);
     meSubdetHit_ = ibooker.book1D("Hit10", "Subdetectors in HCal", 50, 0., 50.);
     meDepthHit_ = ibooker.book1D("Hit11", "Depths in HCal", 20, 0., 20.);
     meEtaHit_ = ibooker.book1D("Hit12", "Eta in HCal", 100, -50., 50.);
     mePhiHit_ = ibooker.book1D("Hit13", "Phi in HCal", 100, 0., 100.);
     meEnergyHit_ = ibooker.book1D("Hit14", "Energy in HCal", 100, 0., 1.);
     meTimeHit_ = ibooker.book1D("Hit15", "Time in HCal", 100, 0., 400.);
     meTimeWHit_ = ibooker.book1D("Hit16", "Time in HCal (E wtd)", 100, 0., 400.);
     meHBDepHit_ = ibooker.book1D("Hit17", "Depths in HB", 20, 0., 20.);
     meHEDepHit_ = ibooker.book1D("Hit18", "Depths in HE", 20, 0., 20.);
     meHODepHit_ = ibooker.book1D("Hit19", "Depths in HO", 20, 0., 20.);
     meHFDepHit_ = ibooker.book1D("Hit20", "Depths in HF", 20, 0., 20.);
     meHBEtaHit_ = ibooker.book1D("Hit21", "Eta in HB", 100, -50., 50.);
     meHEEtaHit_ = ibooker.book1D("Hit22", "Eta in HE", 100, -50., 50.);
     meHOEtaHit_ = ibooker.book1D("Hit23", "Eta in HO", 100, -50., 50.);
     meHFEtaHit_ = ibooker.book1D("Hit24", "Eta in HF", 100, -50., 50.);
     meHBPhiHit_ = ibooker.book1D("Hit25", "Phi in HB", 100, 0., 100.);
     meHEPhiHit_ = ibooker.book1D("Hit26", "Phi in HE", 100, 0., 100.);
     meHOPhiHit_ = ibooker.book1D("Hit27", "Phi in HO", 100, 0., 100.);
     meHFPhiHit_ = ibooker.book1D("Hit28", "Phi in HF", 100, 0., 100.);
     meHBEneHit_ = ibooker.book1D("Hit29", "Energy in HB", 100, 0., 1.);
     meHEEneHit_ = ibooker.book1D("Hit30", "Energy in HE", 100, 0., 1.);
     meHOEneHit_ = ibooker.book1D("Hit31", "Energy in HO", 100, 0., 1.);
     meHFEneHit_ = ibooker.book1D("Hit32", "Energy in HF", 100, 0., 100.);
     meHBTimHit_ = ibooker.book1D("Hit33", "Time in HB", 100, 0., 400.);
     meHETimHit_ = ibooker.book1D("Hit34", "Time in HE", 100, 0., 400.);
     meHOTimHit_ = ibooker.book1D("Hit35", "Time in HO", 100, 0., 400.);
     meHFTimHit_ = ibooker.book1D("Hit36", "Time in HF", 100, 0., 400.);
     mePMTHit_ = ibooker.book1D("Hit37", "Number of Hit in PMT", 1000, 0., 1000.);
     mePMTDepHit_ = ibooker.book1D("Hit38", "Depths in HF PMT", 20, 0., 20.);
     mePMTEtaHit_ = ibooker.book1D("Hit39", "Eta in HF PMT", 100, -50., 50.);
     mePMTPhiHit_ = ibooker.book1D("Hit40", "Phi in HF PMT", 100, 0., 100.);
     mePMTEn1Hit_ = ibooker.book1D("Hit41", "Energy (Ceren) in PMT", 100, 0., 100.);
     mePMTEn2Hit_ = ibooker.book1D("Hit42", "Energy (dE/dx) in PMT", 100, 0., 100.);
     mePMTTimHit_ = ibooker.book1D("Hit43", "Time in HF PMT", 100, 0., 400.);
   }
 
   // Histograms for Layers
   if (checkLay_) {
     meLayerLay_ = ibooker.book1D("Lay01", "Layer # of the Hits", 20, 0., 20.);
     meEtaHLay_ = ibooker.book1D("Lay02", "Eta of the Hits", 100, -5., 5.);
     mePhiHLay_ = ibooker.book1D("Lay03", "Phi of the Hits", 100, -my_pi, my_pi);
     meEneHLay_ = ibooker.book1D("Lay04", "Energy of the Hits", 100, 0., 2.);
     meDepHlay_ = ibooker.book1D("Lay05", "Depth  of the Hits", 10, 0., 10.);
     meTimHLay_ = ibooker.book1D("Lay06", "Time of the Hits", 100, 0., 400.);
     meTimWLay_ = ibooker.book1D("Lay07", "Time (wtd) of Hits", 100, 0., 400.);
     meEtaPhi_ = ibooker.book2D("Lay08", "Phi%Eta of the Hits", 100, -5., 5., 100, -my_pi, my_pi);
 
     meHitELay_ = ibooker.book1D("Lay09", "Hit in Ecal", 1000, 0., 2000.);
     meHitHLay_ = ibooker.book1D("Lay10", "Hit in Hcal", 1000, 0., 2000.);
     meHitTLay_ = ibooker.book1D("Lay11", "Total Hits", 1000, 0., 2000.);
     meEneLLay_ = ibooker.book1D("Lay12", "Energy per layer", 100, 0., 1.);
     int nn = 0;
     for (int i = 0; i < nLayersMAX; i++) {
       sprintf(name, "Layl%d", nn);
       nn++;
       sprintf(title, "Energy deposit in Layer %d", i);
       meEneLay_[i] = ibooker.book1D(name, title, 100, 0., 0.4);
     }
     meLngLay_ = ibooker.book1D("Lay13", "Lonitudinal Shower Profile", 20, 0, 20.);
     meEneDLay_ = ibooker.book1D("Lay14", "Energy per depth", 100, 0., 1.);
     for (int i = 0; i < nDepthsMAX; i++) {
       sprintf(name, "Layl%d", nn);
       nn++;
       sprintf(title, "Energy deposit in Depth %d", i);
       meDepLay_[i] = ibooker.book1D(name, title, 100, 0., 2.);
     }
 
     meEtotLay_ = ibooker.book1D("Lay15", "Total Energy", 100, 0., 1.);
     meEHOLay_ = ibooker.book1D("Lay16", "Energy in HO", 100, 0., 2000.);
     meEHBHELay_ = ibooker.book1D("Lay17", "Energy in HB/HE", 100, 0., 2000.);
     meEFibLLay_ = ibooker.book1D("Lay18", "Energy in HF (Long)", 100, 0., 100.);
     meEFibSLay_ = ibooker.book1D("Lay19", "Energy in HF (Short)", 100, 0., 100.);
     meEHFEmLay_ = ibooker.book1D("Lay20", "EM   energy in HF", 100, 0., 200.);
     meEHFHdLay_ = ibooker.book1D("Lay21", "Had. energy in HF", 100, 0., 200.);
   }
 
   // Histograms for NxN analysis
   if (checkNxN_) {
     meEcalRNxN_ = ibooker.book1D("NxN01", "Energy in ECal (NxN)r", 100, 0., 200.);
     meHcalRNxN_ = ibooker.book1D("NxN02", "Energy in HCal (NxN)r", 100, 0., 200.);
     meHoRNxN_ = ibooker.book1D("NxN03", "Energy in HO (NxN)r", 100, 0., 200.);
     meEtotRNxN_ = ibooker.book1D("NxN04", "Energy Total (NxN)r", 100, 0., 200.);
     meEcalNxN_ = ibooker.book1D("NxN05", "Energy in ECal (NxN)", 100, 0., 200.);
     meHcalNxN_ = ibooker.book1D("NxN06", "Energy in HCal (NxN)", 100, 0., 200.);
     meHoNxN_ = ibooker.book1D("NxN07", "Energy in HO (NxN)", 100, 0., 200.);
     meEtotNxN_ = ibooker.book1D("NxN08", "Energy Total (NxN)", 100, 0., 200.);
     meEiNxN_ = ibooker.book1D("NxN09", "Energy of Hits in (NxN)", 100, 0., 1.);
     meTiNxN_ = ibooker.book1D("NxN10", "Time   of Hits in (NxN)", 100, 0., 400.);
     meTrNxN_ = ibooker.book1D("NxN11", "Dist.  of Hits in (NxN)", 100, 0., 1.);
   }
 
   // Histograms for Jets
   if (checkJet_) {
     meRJet_ = ibooker.book1D("Jet01", "R of Hits", 100, 0., 1.);
     meTJet_ = ibooker.book1D("Jet02", "T of Hits", 100, 0., 200.);
     meEJet_ = ibooker.book1D("Jet03", "E of Hits", 100, 0., 200.);
     meEcalJet_ = ibooker.book1D("Jet04", "Ecal Energy (First Jet)", 100, 0., 200.);
     meHcalJet_ = ibooker.book1D("Jet05", "Hcal Energy (First Jet)", 100, 0., 200.);
     meHoJet_ = ibooker.book1D("Jet06", "Ho   Energy (First Jet)", 100, 0., 200.);
     meEtotJet_ = ibooker.book1D("Jet07", "Total Energy(First Jet)", 100, 0., 200.);
     meEcHcJet_ = ibooker.book2D("Jet08", "Energy in Hcal% Ecal", 100, 0., 200., 100, 0., 200.);
 
     meDetaJet_ = ibooker.book1D("Jet09", "Delta Eta", 100, 0., 2.);
     meDphiJet_ = ibooker.book1D("Jet10", "Delta Phi", 100, 0., 1.);
     meDrJet_ = ibooker.book1D("Jet11", "Delta R", 100, 0., 2.);
     meMassJet_ = ibooker.book1D("Jet12", "Di-jet mass", 100, 0., 200.);
     meEneJet_ = ibooker.book1D("Jet13", "Jet Energy", 100, 0., 200.);
     meEtaJet_ = ibooker.book1D("Jet14", "Jet Eta", 100, -5., 5.);
     mePhiJet_ = ibooker.book1D("Jet15", "Jet Phi", 100, -my_pi, my_pi);
   }
 }