#include <HcalRecHitsValidation.h>

Inheritance diagram for HcalRecHitsValidation:

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

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

	HcalRecHitsValidation (edm::ParameterSet const &conf)

	~HcalRecHitsValidation ()

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

virtual void	beginStream (edm::StreamID id) final

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

	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)

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_

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

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 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)

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 ecalselector_, 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_.

                                                                         {
   // DQM ROOT output
   outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");
   
   if ( outputFile_.size() != 0 ) {
     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"));
   tok_EB_ = consumes<EBRecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEB"));
   tok_EE_ = consumes<EERecHitCollection>(edm::InputTag("ecalRecHit","EcalRecHitsEE"));
   tok_hh_ = consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits","HcalHits"));
 
   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 ( )

Definition at line 56 of file HcalRecHitsValidation.cc.

56 { }

Member Function Documentation

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

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 176 of file HcalRecHitsValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), cdepth, cen, ceta, cieta, cphi, csub, ctime, particleFlowClusterECALTimeSelected_cfi::depth, dR(), ecalselector_, edm::SortedCollection< T, SORT >::end(), stringResolutionProvider_cfi::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, 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, 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;
 
 
       ev.getByToken(tok_EB_, rhitEB);
 
     EcalRecHitCollection::const_iterator RecHit = rhitEB.product()->begin();  
     EcalRecHitCollection::const_iterator RecHitEnd = rhitEB.product()->end();  
     
     for (; RecHit != RecHitEnd ; ++RecHit) {
       EBDetId EBid = EBDetId(RecHit->id());
        
       const CaloCellGeometry* 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;
  
       ev.getByToken(tok_EE_, rhitEE);
 
     RecHit = rhitEE.product()->begin();  
     RecHitEnd = rhitEE.product()->end();  
     
     for (; RecHit != RecHitEnd ; ++RecHit) {
       EEDetId EEid = EEDetId(RecHit->id());
       
       const CaloCellGeometry* 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;
     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 = 
     (HcalGeometry*)(geo->getSubdetectorGeometry(DetId::Hcal,cell.subdet()));
 //      const CaloCellGeometry* cellGeometry =
 //  geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
       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
	)

virtual

Implements DQMEDAnalyzer.

Definition at line 58 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(), and subdet_.

 {
 
   Char_t histo[200];
 
     ib.setCurrentFolder("HcalRecHitsV/HcalRecHitTask");
 
     //======================= 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, 82, -41., 41., -100., 2000., " ");  
       
       sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_E");
       meEnConeEtaProfile_E = ib.bookProfile(histo, histo, 82, -41., 41., -100., 2000., " ");  
       
       sprintf (histo, "HcalRecHitTask_En_rechits_cone_profile_vs_ieta_all_depths_EH");
       meEnConeEtaProfile_EH = ib.bookProfile(histo, histo, 82, -41., 41., -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 649 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 623 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 495 of file HcalRecHitsValidation.cc.

References edm::SortedCollection< T, SORT >::begin(), cdepth, cen, ceta, cieta, ciphi, cphi, csub, ctime, cz, particleFlowClusterECALTimeSelected_cfi::depth, edm::SortedCollection< T, SORT >::end(), stringResolutionProvider_cfi::eta, edm::Event::getByToken(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, iz, j, 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;
     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 = 
     (HcalGeometry*)(geo->getSubdetectorGeometry(DetId::Hcal,cell.subdet()));
 //      const CaloCellGeometry* cellGeometry =
 //  geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
       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();
       if(inteta > 0) inteta -= 1;
       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;
     ev.getByToken(tok_hf_, hfcoll);
 
     for (HFRecHitCollection::const_iterator j = hfcoll->begin(); j != hfcoll->end(); j++) {
       HcalDetId cell(j->id());
       const CaloCellGeometry* 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();
       if(inteta > 0) inteta -= 1;
       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;
     ev.getByToken(tok_ho_, hocoll);
     
     for (HORecHitCollection::const_iterator j = hocoll->begin(); j != hocoll->end(); j++) {
       HcalDetId cell(j->id());
       const CaloCellGeometry* 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();
       if(inteta > 0) inteta -= 1;
       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 633 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;
 
 }