#include <HcalRecHitsValidation.h>

Inheritance diagram for HcalRecHitsValidation:

Public Member Functions
void	analyze (edm::Event const &ev, edm::EventSetup const &c) override

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

	HcalRecHitsValidation (edm::ParameterSet const &conf)

	~HcalRecHitsValidation () override

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

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

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

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

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

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

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase 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 &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
double	dPhiWsign (double phi1, double phi2)

double	dR (double eta1, double phi1, double eta2, double phi2)

virtual void	fillRecHitsTmp (int subdet_, edm::Event const &ev)

double	phi12 (double phi1, double en1, double phi2, double en2)

Private Attributes
std::vector< uint32_t >	cauxstwd

std::vector< int >	cdepth

std::vector< double >	cen

std::vector< double >	ceta

std::vector< int >	cieta

std::vector< int >	ciphi

std::vector< double >	cphi

std::vector< uint32_t >	cstwd

std::vector< int >	csub

std::vector< double >	ctime

std::vector< double >	cz

std::string	ecalselector_

edm::ESHandle< CaloGeometry >	geometry

std::string	hcalselector_

int	imc

int	iz

std::string	mc_

MonitorElement *	meEnConeEtaProfile

MonitorElement *	meEnConeEtaProfile_E

MonitorElement *	meEnConeEtaProfile_EH

MonitorElement *	meEnergyHcalVsEcalHB

MonitorElement *	meEnergyHcalVsEcalHE

MonitorElement *	meRecHitsEnergyHB

MonitorElement *	meRecHitsEnergyHE

MonitorElement *	meRecHitsEnergyHF

MonitorElement *	meRecHitsEnergyHO

MonitorElement *	meRecHitSimHitHB

MonitorElement *	meRecHitSimHitHE

MonitorElement *	meRecHitSimHitHF

MonitorElement *	meRecHitSimHitHFL

MonitorElement *	meRecHitSimHitHFS

MonitorElement *	meRecHitSimHitHO

MonitorElement *	meRecHitSimHitProfileHB

MonitorElement *	meRecHitSimHitProfileHE

MonitorElement *	meRecHitSimHitProfileHF

MonitorElement *	meRecHitSimHitProfileHFL

MonitorElement *	meRecHitSimHitProfileHFS

MonitorElement *	meRecHitSimHitProfileHO

MonitorElement *	meTEprofileHB

MonitorElement *	meTEprofileHB_High

MonitorElement *	meTEprofileHB_Low

MonitorElement *	meTEprofileHE

MonitorElement *	meTEprofileHE_Low

MonitorElement *	meTEprofileHF

MonitorElement *	meTEprofileHF_Low

MonitorElement *	meTEprofileHO

MonitorElement *	meTEprofileHO_High

int	nevtot

std::string	outputFile_

std::string	sign_

int	subdet_

bool	testNumber_

edm::EDGetTokenT< EBRecHitCollection >	tok_EB_

edm::EDGetTokenT< EERecHitCollection >	tok_EE_

edm::EDGetTokenT< edm::HepMCProduct >	tok_evt_

edm::EDGetTokenT< HBHERecHitCollection >	tok_hbhe_

edm::EDGetTokenT< HFRecHitCollection >	tok_hf_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hh_

edm::EDGetTokenT< HORecHitCollection >	tok_ho_

std::string	topFolderName_

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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)

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 56 of file HcalRecHitsValidation.h.

Constructor & Destructor Documentation

HcalRecHitsValidation::HcalRecHitsValidation ( edm::ParameterSet const & conf )

Definition at line 8 of file HcalRecHitsValidation.cc.

References isolatedTracks_cfi::EBRecHitCollectionLabel, ecalselector_, isolatedTracks_cfi::EERecHitCollectionLabel, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcalselector_, imc, iz, mc_, nevtot, outputFile_, sign_, AlCaHLTBitMon_QueryRunRegistry::string, subdet_, testNumber_, tok_EB_, tok_EE_, tok_evt_, tok_hbhe_, tok_hf_, tok_hh_, and tok_ho_.

   : topFolderName_ (conf.getParameter<std::string>("TopFolderName"))
 {
   // DQM ROOT output
   outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");
   
   if ( !outputFile_.empty() ) {
     edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will be saved to '" << outputFile_.c_str() << "'";
   } else {
     edm::LogInfo("OutputInfo") << " Hcal RecHit Task histograms will NOT be saved";
   }
   
   nevtot = 0;
   
   hcalselector_ = conf.getUntrackedParameter<std::string>("hcalselector", "all");
   ecalselector_ = conf.getUntrackedParameter<std::string>("ecalselector", "yes");
   sign_         = conf.getUntrackedParameter<std::string>("sign", "*");
   mc_           = conf.getUntrackedParameter<std::string>("mc", "yes");
   testNumber_   = conf.getParameter<bool>("TestNumber");
 
   //Collections
   tok_hbhe_ = consumes<HBHERecHitCollection>(conf.getUntrackedParameter<edm::InputTag>("HBHERecHitCollectionLabel"));
   tok_hf_   = consumes<HFRecHitCollection>(conf.getUntrackedParameter<edm::InputTag>("HFRecHitCollectionLabel"));
   tok_ho_   = consumes<HORecHitCollection>(conf.getUntrackedParameter<edm::InputTag>("HORecHitCollectionLabel"));
 
   // register for data access
   tok_evt_ = consumes<edm::HepMCProduct>(edm::InputTag("generatorSmeared"));
   edm::InputTag EBRecHitCollectionLabel = conf.getParameter<edm::InputTag>("EBRecHitCollectionLabel");
   tok_EB_ = consumes<EBRecHitCollection>(EBRecHitCollectionLabel);
   edm::InputTag EERecHitCollectionLabel = conf.getParameter<edm::InputTag>("EERecHitCollectionLabel");
   tok_EE_ = consumes<EERecHitCollection>(EERecHitCollectionLabel);
 
   tok_hh_ = consumes<edm::PCaloHitContainer>(conf.getUntrackedParameter<edm::InputTag>("SimHitCollectionLabel"));
 
   subdet_ = 5;
   if (hcalselector_ == "noise") subdet_ = 0;
   if (hcalselector_ == "HB"   ) subdet_ = 1;
   if (hcalselector_ == "HE"   ) subdet_ = 2;
   if (hcalselector_ == "HO"   ) subdet_ = 3;
   if (hcalselector_ == "HF"   ) subdet_ = 4;
   if (hcalselector_ == "all"  ) subdet_ = 5;
   if (hcalselector_ == "ZS"   ) subdet_ = 6;
 
   iz = 1;
   if(sign_ == "-") iz = -1;
   if(sign_ == "*") iz = 0;
 
   imc = 1;
   if(mc_ == "no") imc = 0;
 
 }

HcalRecHitsValidation::~HcalRecHitsValidation ( )

override

Definition at line 61 of file HcalRecHitsValidation.cc.

61 { }

Member Function Documentation

void HcalRecHitsValidation::analyze	(	edm::Event const &	ev,
		edm::EventSetup const &	c
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 181 of file HcalRecHitsValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), cdepth, cen, ceta, cieta, cphi, csub, ctime, egammaForCoreTracking_cff::depth, dR(), ecalselector_, edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), ALCARECOTkAlBeamHalo_cff::etaMax, MonitorElement::Fill(), fillRecHitsTmp(), edm::EventSetup::get(), edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), HcalGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, imc, edm::HandleBase::isValid(), MuonErrorMatrixAnalyzer_cfi::maxPt, meEnConeEtaProfile, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meRecHitsEnergyHB, meRecHitsEnergyHE, meRecHitsEnergyHF, meRecHitsEnergyHO, meRecHitSimHitHB, meRecHitSimHitHE, meRecHitSimHitHF, meRecHitSimHitHFL, meRecHitSimHitHFS, meRecHitSimHitHO, meRecHitSimHitProfileHB, meRecHitSimHitProfileHE, meRecHitSimHitProfileHF, meRecHitSimHitProfileHFL, meRecHitSimHitProfileHFS, meRecHitSimHitProfileHO, meTEprofileHB, meTEprofileHB_High, meTEprofileHB_Low, meTEprofileHE, meTEprofileHE_Low, meTEprofileHF, meTEprofileHF_Low, meTEprofileHO, meTEprofileHO_High, nevtot, npart, AlCaHLTBitMon_ParallelJobs::p, phi, PV3DBase< T, PVType, FrameType >::phi(), edm::Handle< T >::product(), EnergyCorrector::pt, alignCSCRings::r, HcalHitRelabeller::relabel(), HcalDetId::subdet(), subdet_, lumiQTWidget::t, testNumber_, tok_EB_, tok_EE_, tok_evt_, and tok_hh_.

                                                                               {
 
   using namespace edm;
 
   edm::ESHandle<HcalDDDRecConstants> pHRNDC;
   c.get<HcalRecNumberingRecord>().get( pHRNDC );
   const HcalDDDRecConstants* hcons = &(*pHRNDC);
 
   // cuts for each subdet_ector mimiking  "Scheme B"
   //  double cutHB = 0.9, cutHE = 1.4, cutHO = 1.1, cutHFL = 1.2, cutHFS = 1.8; 
 
   // energy in HCAL
   double eHcal        = 0.;
   double eHcalCone    = 0.;  
   double eHcalConeHB  = 0.;  
   double eHcalConeHE  = 0.;  
   double eHcalConeHO  = 0.;  
   double eHcalConeHF  = 0.;  
   double eHcalConeHFL = 0.;  
   double eHcalConeHFS = 0.;  
 
   // Total numbet of RecHits in HCAL, in the cone, above 1 GeV theshold
   int nrechits       = 0;
   int nrechitsCone   = 0;
   int nrechitsThresh = 0;
 
   // energy in ECAL
   double eEcal       = 0.;
   double eEcalB      = 0.;
   double eEcalE      = 0.;
   double eEcalCone   = 0.;
   int numrechitsEcal = 0;
 
   // MC info 
   double phi_MC = -999999.;  // phi of initial particle from HepMC
   double eta_MC = -999999.;  // eta of initial particle from HepMC
 
   // HCAL energy around MC eta-phi at all depths;
   double partR = 0.3;
 
   if(imc != 0) { 
 
      edm::Handle<edm::HepMCProduct> evtMC;
      ev.getByToken(tok_evt_,evtMC);  // generator in late 310_preX
      if (!evtMC.isValid()) {
         edm::LogInfo("HcalRecHitsValidation") << "no HepMCProduct found";    
      } else {
         //    std::cout << "*** source HepMCProduct found"<< std::endl;
      }  
 
      // MC particle with highest pt is taken as a direction reference  
      double maxPt = -99999.;
      int npart    = 0;
      const HepMC::GenEvent * myGenEvent = evtMC->GetEvent();
      for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end(); ++p ) {
        double phip = (*p)->momentum().phi();
        double etap = (*p)->momentum().eta();
        //    phi_MC = phip;
        //    eta_MC = etap;
        double pt  = (*p)->momentum().perp();
        if(pt > maxPt) { npart++; maxPt = pt; phi_MC = phip; eta_MC = etap; }
      }
      //  std::cout << "*** Max pT = " << maxPt <<  std::endl;  
 
   }
 
   //   std::cout << "*** 2" << std::endl; 
 
   c.get<CaloGeometryRecord>().get (geometry);
 
   // Fill working vectors of HCAL RecHits quantities (all of these are drawn)
   fillRecHitsTmp(subdet_, ev); 
 
   //  std::cout << "*** 3" << std::endl; 
 
   //===========================================================================
   // IN ALL other CASES : ieta-iphi maps 
   //===========================================================================
 
   // ECAL 
   if(ecalselector_ == "yes" && (subdet_ == 1 || subdet_ == 2 || subdet_ == 5)) {
     Handle<EBRecHitCollection> rhitEB;
 
     EcalRecHitCollection::const_iterator RecHit;
     EcalRecHitCollection::const_iterator RecHitEnd;
     
 
     if ( ev.getByToken(tok_EB_, rhitEB) ) {
 
       RecHit = rhitEB.product()->begin();  
       RecHitEnd = rhitEB.product()->end();  
     
       for (; RecHit != RecHitEnd ; ++RecHit) {
     EBDetId EBid = EBDetId(RecHit->id());
        
     auto cellGeometry = geometry->getSubdetectorGeometry(EBid)->getGeometry (EBid) ;
     double eta = cellGeometry->getPosition ().eta () ;
     double phi = cellGeometry->getPosition ().phi () ;
     double en  = RecHit->energy();
     eEcal  += en;
     eEcalB += en;
     
     double r   = dR(eta_MC, phi_MC, eta, phi);
     if( r < partR)  {
       eEcalCone += en;
       numrechitsEcal++; 
     }
       }
       
     }
 
     Handle<EERecHitCollection> rhitEE;
  
     if ( ev.getByToken(tok_EE_, rhitEE) ) {
       
       RecHit = rhitEE.product()->begin();  
       RecHitEnd = rhitEE.product()->end();  
       
       for (; RecHit != RecHitEnd ; ++RecHit) {
     EEDetId EEid = EEDetId(RecHit->id());
     
     auto cellGeometry = geometry->getSubdetectorGeometry(EEid)->getGeometry (EEid) ;
     double eta = cellGeometry->getPosition ().eta () ;
     double phi = cellGeometry->getPosition ().phi () ;  
     double en   = RecHit->energy();
     eEcal  += en;
     eEcalE += en;
     
     double r   = dR(eta_MC, phi_MC, eta, phi);
     if( r < partR)  {
       eEcalCone += en;
       numrechitsEcal++; 
     }
       }
     }
   }     // end of ECAL selection 
 
 
   //     std::cout << "*** 4" << std::endl; 
 
 
   //===========================================================================
   // SUBSYSTEMS,  
   //===========================================================================
   
   if ((subdet_ != 6) && (subdet_ != 0)) {
 
     //       std::cout << "*** 6" << std::endl; 
     
     
     double HcalCone    = 0.;
 
     int ietaMax   =  9999;
     double etaMax =  9999.;
 
 
     //   CYCLE over cells ====================================================
 
     for (unsigned int i = 0; i < cen.size(); i++) {
       int sub    = csub[i];
       int depth  = cdepth[i];
       double eta = ceta[i]; 
       double phi = cphi[i]; 
       double en  = cen[i]; 
       double t   = ctime[i];
       int   ieta = cieta[i];
 
       nrechits++;       
       eHcal += en;
       if(en > 1. ) nrechitsThresh++;
 
       double r    = dR(eta_MC, phi_MC, eta, phi);
       if( r < partR ){
         if(sub == 1)   eHcalConeHB += en;
         if(sub == 2)   eHcalConeHE += en;
         if(sub == 3)   eHcalConeHO += en;
         if(sub == 4) {
       eHcalConeHF += en;
       if (depth == 1) eHcalConeHFL += en;
       else            eHcalConeHFS += en;
     }
     eHcalCone += en;
     nrechitsCone++;
 
         HcalCone += en;
 
         // alternative: ietamax -> closest to MC eta  !!!
     float eta_diff = fabs(eta_MC - eta);
     if(eta_diff < etaMax) {
       etaMax  = eta_diff; 
           ietaMax = ieta; 
     }
       }
       
       //The energy and overall timing histos are drawn while
       //the ones split by depth are not
       if(sub == 1 && (subdet_ == 1 || subdet_ == 5)) {  
     meRecHitsEnergyHB->Fill(en);
     
     meTEprofileHB_Low->Fill(en, t);
     meTEprofileHB->Fill(en, t);
     meTEprofileHB_High->Fill(en, t);
 
       }     
       if(sub == 2 && (subdet_ == 2 || subdet_ == 5)) {  
     meRecHitsEnergyHE->Fill(en);
 
     meTEprofileHE_Low->Fill(en, t);
     meTEprofileHE->Fill(en, t);
 
       }
       if(sub == 4 && (subdet_ == 4 || subdet_ == 5)) {  
     meRecHitsEnergyHF->Fill(en);      
 
     meTEprofileHF_Low->Fill(en, t);
     meTEprofileHF->Fill(en, t);
 
       }
       if(sub == 3 && (subdet_ == 3 || subdet_ == 5)) {  
     meRecHitsEnergyHO->Fill(en);
 
     meTEprofileHO->Fill(en, t);
     meTEprofileHO_High->Fill(en, t);
       }
     }
 
     if(imc != 0) {
       meEnConeEtaProfile       ->Fill(double(ietaMax),  HcalCone);   // 
       meEnConeEtaProfile_E     ->Fill(double(ietaMax), eEcalCone);   
       meEnConeEtaProfile_EH    ->Fill(double(ietaMax),  HcalCone+eEcalCone); 
     }
 
     //     std::cout << "*** 7" << std::endl; 
 
     
   }
 
   //SimHits vs. RecHits
   const CaloGeometry* geo = geometry.product();
   if(subdet_ > 0 && subdet_ < 6 && imc !=0) {  // not noise 
 
     edm::Handle<PCaloHitContainer> hcalHits;
     if ( ev.getByToken(tok_hh_,hcalHits) ) {
       const PCaloHitContainer * SimHitResult = hcalHits.product () ;
       
       double enSimHits    = 0.;
       double enSimHitsHB  = 0.;
       double enSimHitsHE  = 0.;
       double enSimHitsHO  = 0.;
       double enSimHitsHF  = 0.;
       double enSimHitsHFL = 0.;
       double enSimHitsHFS = 0.;
       // sum of SimHits in the cone 
       
       for (std::vector<PCaloHit>::const_iterator SimHits = SimHitResult->begin () ; SimHits != SimHitResult->end(); ++SimHits) {
     
     int sub, depth;
     HcalDetId cell; 
     
     if (testNumber_) cell = HcalHitRelabeller::relabel(SimHits->id(),hcons);
     else cell = HcalDetId(SimHits->id());
     
     sub          = cell.subdet();
     depth        = cell.depth();
     
     if(sub != subdet_ && subdet_ != 5) continue; //If we are not looking at all of the subdetectors and the simhit doesn't come from the specific subdetector of interest, then we won't do any thing with it 
     
     const HcalGeometry* cellGeometry = 
       dynamic_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,cell.subdet()));
     double etaS = cellGeometry->getPosition(cell).eta () ;
     double phiS = cellGeometry->getPosition(cell).phi () ;
     double en   = SimHits->energy();    
     
     double r  = dR(eta_MC, phi_MC, etaS, phiS);
     
     if ( r < partR ){ // just energy in the small cone
       
       enSimHits += en;
       if(sub == static_cast<int>(HcalBarrel)) enSimHitsHB += en; 
       if(sub == static_cast<int>(HcalEndcap)) enSimHitsHE += en; 
       if(sub == static_cast<int>(HcalOuter)) enSimHitsHO += en; 
       if(sub == static_cast<int>(HcalForward)) {
         enSimHitsHF += en;
         if(depth == 1) enSimHitsHFL += en;
         else           enSimHitsHFS += en;
       } 
     }
       }
       
       // Now some histos with SimHits
       
       if(subdet_ == 4 || subdet_ == 5) {
     meRecHitSimHitHF->Fill( enSimHitsHF, eHcalConeHF );
     meRecHitSimHitProfileHF->Fill( enSimHitsHF, eHcalConeHF);
     
     meRecHitSimHitHFL->Fill( enSimHitsHFL, eHcalConeHFL );
     meRecHitSimHitProfileHFL->Fill( enSimHitsHFL, eHcalConeHFL);
     meRecHitSimHitHFS->Fill( enSimHitsHFS, eHcalConeHFS );
     meRecHitSimHitProfileHFS->Fill( enSimHitsHFS, eHcalConeHFS);       
       }
       if(subdet_ == 1  || subdet_ == 5) { 
     meRecHitSimHitHB->Fill( enSimHitsHB,eHcalConeHB );
     meRecHitSimHitProfileHB->Fill( enSimHitsHB,eHcalConeHB);
       }
       if(subdet_ == 2  || subdet_ == 5) { 
     meRecHitSimHitHE->Fill( enSimHitsHE,eHcalConeHE );
     meRecHitSimHitProfileHE->Fill( enSimHitsHE,eHcalConeHE);
       }
       if(subdet_ == 3  || subdet_ == 5) { 
     meRecHitSimHitHO->Fill( enSimHitsHO,eHcalConeHO );
     meRecHitSimHitProfileHO->Fill( enSimHitsHO,eHcalConeHO);
       }
     
     }
   }
 
   nevtot++;
 }

void HcalRecHitsValidation::bookHistograms	(	DQMStore::IBooker &	ib,
		edm::Run const &	run,
		edm::EventSetup const &	es
	)

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 63 of file HcalRecHitsValidation.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), trackerHits::histo, imc, meEnConeEtaProfile, meEnConeEtaProfile_E, meEnConeEtaProfile_EH, meRecHitsEnergyHB, meRecHitsEnergyHE, meRecHitsEnergyHF, meRecHitsEnergyHO, meRecHitSimHitHB, meRecHitSimHitHE, meRecHitSimHitHF, meRecHitSimHitHFL, meRecHitSimHitHFS, meRecHitSimHitHO, meRecHitSimHitProfileHB, meRecHitSimHitProfileHE, meRecHitSimHitProfileHF, meRecHitSimHitProfileHFL, meRecHitSimHitProfileHFS, meRecHitSimHitProfileHO, meTEprofileHB, meTEprofileHB_High, meTEprofileHB_Low, meTEprofileHE, meTEprofileHE_Low, meTEprofileHF, meTEprofileHF_Low, meTEprofileHO, meTEprofileHO_High, DQMStore::IBooker::setCurrentFolder(), subdet_, and topFolderName_.

 {
 
   Char_t histo[200];
 
   ib.setCurrentFolder(topFolderName_);
 
     //======================= Now various cases one by one ===================
 
     //Histograms drawn for single pion scan
     if(subdet_ != 0 && imc != 0) { // just not for noise  
       sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths");
       meEnConeEtaProfile = ib.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");  
       
       sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_E");
       meEnConeEtaProfile_E = ib.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");  
       
       sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_EH");
       meEnConeEtaProfile_EH = ib.bookProfile(histo, histo, 83, -41.5, 41.5, -100., 2000., " ");  
     }
 
     // ************** HB **********************************
     if (subdet_ == 1 || subdet_ == 5 ){
 
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_HB" ) ;
       meRecHitsEnergyHB = ib.book1D(histo, histo, 2010 , -10. , 2000.); 
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HB" ) ;
       meTEprofileHB = ib.bookProfile(histo, histo, 150, -5., 295., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HB" ) ;
       meTEprofileHB_Low = ib.bookProfile(histo, histo, 150, -5., 295., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_High_HB" ) ;
       meTEprofileHB_High = ib.bookProfile(histo, histo, 150, -5., 295., 48., 92., " "); 
 
       if(imc != 0) {
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HB");
         meRecHitSimHitHB = ib.book2D(histo, histo, 120, 0., 1.2,  300, 0., 150.);
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HB");
         meRecHitSimHitProfileHB = ib.bookProfile(histo, histo, 120, 0., 1.2, 0., 500., " ");  
       } 
 
     }
     
     // ********************** HE ************************************
     if ( subdet_ == 2 || subdet_ == 5 ){
 
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_HE" ) ;
       meRecHitsEnergyHE = ib.book1D(histo, histo, 2010, -10., 2000.);
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HE" ) ;
       meTEprofileHE_Low = ib.bookProfile(histo, histo, 80, -5., 75., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HE" ) ;
       meTEprofileHE = ib.bookProfile(histo, histo, 200, -5., 2995., -48., 92., " "); 
       
       if(imc != 0) {
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HE");
         meRecHitSimHitHE = ib.book2D(histo, histo, 120, 0., 0.6,  300, 0., 150.);
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HE");
         meRecHitSimHitProfileHE = ib.bookProfile(histo, histo, 120, 0., 0.6, 0., 500., " ");  
       }
 
     }
 
     // ************** HO ****************************************
     if ( subdet_ == 3 || subdet_ == 5  ){
 
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_HO" ) ;
       meRecHitsEnergyHO = ib.book1D(histo, histo, 2010 , -10. , 2000.);
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HO" ) ;
       meTEprofileHO = ib.bookProfile(histo, histo, 60, -5., 55., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_High_HO" ) ;
       meTEprofileHO_High = ib.bookProfile(histo, histo, 100, -5., 995., -48., 92., " "); 
       
       if(imc != 0) {
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HO");
         meRecHitSimHitHO = ib.book2D(histo, histo, 150, 0., 1.5,  350, 0., 350.);
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HO");
         meRecHitSimHitProfileHO = ib.bookProfile(histo, histo, 150, 0., 1.5, 0., 500., " ");  
       }
 
     }   
   
     // ********************** HF ************************************
     if ( subdet_ == 4 || subdet_ == 5 ){
 
       sprintf (histo, "HcalRecHitTask_energy_of_rechits_HF" ) ;
       meRecHitsEnergyHF = ib.book1D(histo, histo, 2010 , -10. , 2000.); 
       
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_Low_HF" ) ;
       meTEprofileHF_Low = ib.bookProfile(histo, histo, 100, -5., 195., -48., 92., " "); 
 
       sprintf (histo, "HcalRecHitTask_timing_vs_energy_profile_HF" ) ;
       meTEprofileHF = ib.bookProfile(histo, histo, 200, -5., 995., -48., 92., " "); 
 
       if(imc != 0) {
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HF");
         meRecHitSimHitHF  = ib.book2D(histo, histo, 50, 0., 50., 150, 0., 150.);      
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HFL");
         meRecHitSimHitHFL = ib.book2D(histo, histo, 50, 0., 50., 150, 0., 150.);      
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_HFS");
         meRecHitSimHitHFS = ib.book2D(histo, histo, 50, 0., 50., 150, 0., 150.);
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HF");
         meRecHitSimHitProfileHF  = ib.bookProfile(histo, histo, 50, 0., 50., 0., 500., " ");  
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HFL");
         meRecHitSimHitProfileHFL = ib.bookProfile(histo, histo, 50, 0., 50., 0., 500., " ");  
         sprintf (histo, "HcalRecHitTask_energy_rechits_vs_simhits_profile_HFS");
         meRecHitSimHitProfileHFS = ib.bookProfile(histo, histo, 50, 0., 50., 0., 500., " ");  
       }
 
     }
 
 }

double HcalRecHitsValidation::dPhiWsign	(	double	phi1,
		double	phi2
	)

private

Definition at line 656 of file HcalRecHitsValidation.cc.

References DEFINE_FWK_MODULE, PI, and tmp.

                                                                 {
   // clockwise      phi2 w.r.t phi1 means "+" phi distance
   // anti-clockwise phi2 w.r.t phi1 means "-" phi distance 
 
   double PI = 3.1415926535898;
   double a1 = phi1; double a2  = phi2;
   double tmp =  a2 - a1;
   if( a1*a2 < 0.) {
     if(a1 > 0.5 * PI)  tmp += 2.*PI;
     if(a2 > 0.5 * PI)  tmp -= 2.*PI;
   }
   return tmp;
 
 }

double HcalRecHitsValidation::dR	(	double	eta1,
		double	phi1,
		double	eta2,
		double	phi2
	)

private

Definition at line 630 of file HcalRecHitsValidation.cc.

References PI, mathSSE::sqrt(), and tmp.

Referenced by analyze().

                                                                                    { 
   double PI = 3.1415926535898;
   double deltaphi= phi1 - phi2;
   if( phi2 > phi1 ) { deltaphi= phi2 - phi1;}
   if(deltaphi > PI) { deltaphi = 2.*PI - deltaphi;}
   double deltaeta = eta2 - eta1;
   double tmp = sqrt(deltaeta* deltaeta + deltaphi*deltaphi);
   return tmp;
 }

void HcalRecHitsValidation::fillRecHitsTmp	(	int	subdet_,
		edm::Event const &	ev
	)

privatevirtual

Definition at line 502 of file HcalRecHitsValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), cdepth, cen, ceta, cieta, ciphi, cphi, csub, ctime, cz, egammaForCoreTracking_cff::depth, edm::SortedCollection< T, SORT >::end(), PVValHelper::eta, edm::Event::getByToken(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, iz, phi, lumiQTWidget::t, tok_hbhe_, tok_hf_, and tok_ho_.

Referenced by analyze().

                                                                          {
 
   using namespace edm;
   
   // initialize data vectors
   csub.clear();
   cen.clear();
   ceta.clear();
   cphi.clear();
   ctime.clear();
   cieta.clear();
   ciphi.clear();
   cdepth.clear();
   cz.clear();
 
   if( subdet_ == 1 || subdet_ == 2  || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
     
     //HBHE
     edm::Handle<HBHERecHitCollection> hbhecoll;
     if ( ev.getByToken(tok_hbhe_, hbhecoll) ) {
       const CaloGeometry* geo = geometry.product();
       
       for (HBHERecHitCollection::const_iterator j=hbhecoll->begin(); j != hbhecoll->end(); j++) {
     
     HcalDetId cell(j->id());
     const HcalGeometry* cellGeometry = 
       dynamic_cast<const HcalGeometry*>(geo->getSubdetectorGeometry(DetId::Hcal,cell.subdet()));
     double eta  = cellGeometry->getPosition(cell).eta () ;
     double phi  = cellGeometry->getPosition(cell).phi () ;
     double zc   = cellGeometry->getPosition(cell).z ();
     int sub     = cell.subdet();
     int depth   = cell.depth();
     int inteta  = cell.ieta();
     int intphi  = cell.iphi()-1;
     double en   = j->energy();
     double t    = j->time();
     
     if((iz > 0 && eta > 0.) || (iz < 0 && eta <0.) || iz == 0) { 
       
       csub.push_back(sub);
       cen.push_back(en);
       ceta.push_back(eta);
       cphi.push_back(phi);
       ctime.push_back(t);
       cieta.push_back(inteta);
       ciphi.push_back(intphi);
       cdepth.push_back(depth);
       cz.push_back(zc);
     }
       }
     }
     
   }
 
   if( subdet_ == 4 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
 
     //HF
     edm::Handle<HFRecHitCollection> hfcoll;
     if ( ev.getByToken(tok_hf_, hfcoll) ) {
       
       for (HFRecHitCollection::const_iterator j = hfcoll->begin(); j != hfcoll->end(); j++) {
     
     HcalDetId cell(j->id());
     auto cellGeometry = geometry->getSubdetectorGeometry(cell)->getGeometry (cell) ;
     
     double eta   = cellGeometry->getPosition().eta () ;
     double phi   = cellGeometry->getPosition().phi () ;
     double zc     = cellGeometry->getPosition().z ();
     int sub      = cell.subdet();
     int depth    = cell.depth();
     int inteta   = cell.ieta();
     int intphi   = cell.iphi()-1;
     double en    = j->energy();
     double t     = j->time();
     
     if((iz > 0 && eta > 0.) || (iz < 0 && eta <0.) || iz == 0) { 
       
       csub.push_back(sub);
       cen.push_back(en);
       ceta.push_back(eta);
       cphi.push_back(phi);
       ctime.push_back(t);
       cieta.push_back(inteta);
       ciphi.push_back(intphi);
       cdepth.push_back(depth);
       cz.push_back(zc);
     }
       }
     }
   }
 
   //HO
   if( subdet_ == 3 || subdet_ == 5 || subdet_ == 6 || subdet_ == 0) {
   
     edm::Handle<HORecHitCollection> hocoll;
     if ( ev.getByToken(tok_ho_, hocoll) ) {
     
       for (HORecHitCollection::const_iterator j = hocoll->begin(); j != hocoll->end(); j++) {
     
     HcalDetId cell(j->id());
     auto cellGeometry = geometry->getSubdetectorGeometry(cell)->getGeometry (cell) ;
     
     double eta   = cellGeometry->getPosition().eta () ;
     double phi   = cellGeometry->getPosition().phi () ;
     double zc    = cellGeometry->getPosition().z ();
     int sub      = cell.subdet();
     int depth    = cell.depth();
     int inteta   = cell.ieta();
     int intphi   = cell.iphi()-1;
     double t     = j->time();
     double en    = j->energy();
     
     if((iz > 0 && eta > 0.) || (iz < 0 && eta <0.) || iz == 0) { 
       csub.push_back(sub);
       cen.push_back(en);
       ceta.push_back(eta);
       cphi.push_back(phi);
       ctime.push_back(t);
       cieta.push_back(inteta);
       ciphi.push_back(intphi);
       cdepth.push_back(depth);
       cz.push_back(zc);
     }
       }
     }
   }
 }

double HcalRecHitsValidation::phi12	(	double	phi1,
		double	en1,
		double	phi2,
		double	en2
	)

private

Definition at line 640 of file HcalRecHitsValidation.cc.

References PI, and tmp.

                                                                                     {
   // weighted mean value of phi1 and phi2
   
   double tmp;
   double PI = 3.1415926535898;
   double a1 = phi1; double a2  = phi2;
 
   if( a1 > 0.5*PI  && a2 < 0.) a2 += 2*PI; 
   if( a2 > 0.5*PI  && a1 < 0.) a1 += 2*PI; 
   tmp = (a1 * en1 + a2 * en2)/(en1 + en2);
   if(tmp > PI) tmp -= 2.*PI; 
  
   return tmp;
 
 }