#include <HcalSimHitStudy.h>

Inheritance diagram for HcalSimHitStudy:

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

	HcalSimHitStudy (const edm::ParameterSet &ps)

	~HcalSimHitStudy () override

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

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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 () noexcept(false) 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)

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

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

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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

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

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

void	consumesMany (const TypeToGet &id)

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

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

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

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

Private Attributes
bool	checkHit_

std::string	g4Label

std::string	hcalHits

const HcalDDDRecConstants *	hcons

bool	hep17_

int	ieta_bins_HB

int	ieta_bins_HE

int	ieta_bins_HF

int	ieta_bins_HO

float	ieta_max_HB

float	ieta_max_HE

float	ieta_max_HF

float	ieta_max_HO

float	ieta_min_HB

float	ieta_min_HE

float	ieta_min_HF

float	ieta_min_HO

int	iphi_bins

float	iphi_max

float	iphi_min

int	maxDepth_

int	maxDepthHB_

int	maxDepthHE_

int	maxDepthHF_

int	maxDepthHO_

MonitorElement *	meAllNHit_

MonitorElement *	meBadDetHit_

MonitorElement *	meBadIdHit_

MonitorElement *	meBadSubHit_

MonitorElement *	meDepthHit_

MonitorElement *	meDetectHit_

MonitorElement *	meEnergyHit_

MonitorElement *	meEtaHit_

MonitorElement *	meEtaPhiHit_

std::vector< MonitorElement * >	meEtaPhiHitDepth_

MonitorElement *	meHBDepHit_

MonitorElement *	meHBEneHit2_

MonitorElement *	meHBEneHit_

MonitorElement *	meHBEneMap_

MonitorElement *	meHBEneSum_

MonitorElement *	meHBEneSum_vs_ieta_

MonitorElement *	meHBEtaHit_

MonitorElement *	meHBL10Ene_

MonitorElement *	meHBL10EneP_

MonitorElement *	meHBNHit_

MonitorElement *	meHBPhiHit_

MonitorElement *	meHBTimHit_

MonitorElement *	meHEDepHit_

MonitorElement *	meHEEneHit2_

MonitorElement *	meHEEneHit_

MonitorElement *	meHEEneMap_

MonitorElement *	meHEEneSum_

MonitorElement *	meHEEneSum_vs_ieta_

MonitorElement *	meHEEtaHit_

MonitorElement *	meHEL10Ene_

MonitorElement *	meHEL10EneP_

MonitorElement *	meHENHit_

MonitorElement *	meHEP17EneHit2_

MonitorElement *	meHEP17EneHit_

MonitorElement *	meHEPhiHit_

MonitorElement *	meHETimHit_

MonitorElement *	meHFDepHit_

MonitorElement *	meHFDepHitw_

MonitorElement *	meHFEneHit2_

MonitorElement *	meHFEneHit_

MonitorElement *	meHFEneMap_

MonitorElement *	meHFEneSum_

MonitorElement *	meHFEneSum_vs_ieta_

MonitorElement *	meHFEtaHit_

MonitorElement *	meHFL10Ene_

MonitorElement *	meHFL10EneP_

MonitorElement *	meHFNHit_

MonitorElement *	meHFPhiHit_

MonitorElement *	meHFTimHit_

MonitorElement *	meHODepHit_

MonitorElement *	meHOEneHit2_

MonitorElement *	meHOEneHit_

MonitorElement *	meHOEneMap_

MonitorElement *	meHOEneSum_

MonitorElement *	meHOEneSum_vs_ieta_

MonitorElement *	meHOEtaHit_

MonitorElement *	meHOL10Ene_

MonitorElement *	meHOL10EneP_

MonitorElement *	meHONHit_

MonitorElement *	meHOPhiHit_

MonitorElement *	meHOTimHit_

MonitorElement *	mePhiHit_

MonitorElement *	mePhiHitb_

MonitorElement *	meSubdetHit_

MonitorElement *	meTimeHit_

MonitorElement *	meTimeWHit_

std::string	outFile_

bool	testNumber_

edm::EDGetTokenT< edm::PCaloHitContainer >	tok_hits_

bool	verbose_

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 31 of file HcalSimHitStudy.h.

Constructor & Destructor Documentation

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

Definition at line 7 of file HcalSimHitStudy.cc.

References checkHit_, g4Label, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hcalHits, hep17_, outFile_, AlCaHLTBitMon_QueryRunRegistry::string, testNumber_, tok_hits_, and verbose_.

                                                           {
 
   g4Label  = ps.getUntrackedParameter<std::string>("moduleLabel","g4SimHits");
   hcalHits = ps.getUntrackedParameter<std::string>("HitCollection","HcalHits");
   outFile_ = ps.getUntrackedParameter<std::string>("outputFile", "hcHit.root");
   verbose_ = ps.getUntrackedParameter<bool>("Verbose", false);
   testNumber_= ps.getParameter<bool>("TestNumber");
   hep17_   = ps.getParameter<bool>("hep17");
   checkHit_= true;
 
   tok_hits_ = consumes<edm::PCaloHitContainer>(edm::InputTag(g4Label,hcalHits));
 
   edm::LogInfo("HcalSim") << "Module Label: " << g4Label << "   Hits: "
               << hcalHits << " / "<< checkHit_ 
               << "   Output: " << outFile_;
 
 }

HcalSimHitStudy::~HcalSimHitStudy ( )

override

Definition at line 25 of file HcalSimHitStudy.cc.

25 {}

Member Function Documentation

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

overrideprotected

Definition at line 183 of file HcalSimHitStudy.cc.

References analyzeHits(), checkHit_, edm::EventID::event(), edm::Event::getByToken(), edm::EventBase::id(), edm::HandleBase::isValid(), LogDebug, edm::EventID::run(), and tok_hits_.

                                                                       {
 
   LogDebug("HcalSim") << "Run = " << e.id().run() << " Event = " 
               << e.id().event();
 
   std::vector<PCaloHit>               caloHits;
   edm::Handle<edm::PCaloHitContainer> hitsHcal;
 
   bool getHits = false;
   if (checkHit_) {
     e.getByToken(tok_hits_,hitsHcal); 
     if (hitsHcal.isValid()) getHits = true;
   }
 
   LogDebug("HcalSim") << "HcalValidation: Input flags Hits " << getHits;
 
   if (getHits) {
     caloHits.insert(caloHits.end(),hitsHcal->begin(),hitsHcal->end());
     LogDebug("HcalSim") << "HcalValidation: Hit buffer " 
             << caloHits.size(); 
     analyzeHits (caloHits);
   }
 }

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

protected

Definition at line 207 of file HcalSimHitStudy.cc.

References TauDecayModes::dec, egammaForCoreTracking_cff::depth, HcalDetId::depth(), DetId::det(), PVValHelper::eta, MonitorElement::Fill(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, hcons, hep17_, mps_fire::i, HcalDetId::ieta(), ieta_bins_HB, ieta_bins_HE, ieta_bins_HF, ieta_bins_HO, ieta_min_HB, ieta_min_HE, ieta_min_HF, ieta_min_HO, createfilelist::int, HcalDetId::iphi(), iphi_bins, LogDebug, meAllNHit_, meBadDetHit_, meBadIdHit_, meBadSubHit_, meDepthHit_, meDetectHit_, meEnergyHit_, meEtaHit_, meEtaPhiHit_, meEtaPhiHitDepth_, meHBDepHit_, meHBEneHit2_, meHBEneHit_, meHBEneMap_, meHBEneSum_, meHBEneSum_vs_ieta_, meHBEtaHit_, meHBL10Ene_, meHBL10EneP_, meHBNHit_, meHBPhiHit_, meHBTimHit_, meHEDepHit_, meHEEneHit2_, meHEEneHit_, meHEEneMap_, meHEEneSum_, meHEEneSum_vs_ieta_, meHEEtaHit_, meHEL10Ene_, meHEL10EneP_, meHENHit_, meHEP17EneHit2_, meHEP17EneHit_, meHEPhiHit_, meHETimHit_, meHFDepHit_, meHFDepHitw_, meHFEneHit2_, meHFEneHit_, meHFEneMap_, meHFEneSum_, meHFEneSum_vs_ieta_, meHFEtaHit_, meHFL10Ene_, meHFL10EneP_, meHFNHit_, meHFPhiHit_, meHFTimHit_, meHODepHit_, meHOEneHit2_, meHOEneHit_, meHOEneMap_, meHOEneSum_, meHOEneSum_vs_ieta_, meHOEtaHit_, meHOL10Ene_, meHOL10EneP_, meHONHit_, meHOPhiHit_, meHOTimHit_, mePhiHit_, mePhiHitb_, meSubdetHit_, meTimeHit_, meTimeWHit_, phi, HcalHitRelabeller::relabel(), HcalDetId::subdet(), testNumber_, and ntuplemaker::time.

Referenced by analyze().

                                                             {
 
   int nHit = hits.size();
   int nHB=0, nHE=0, nHO=0, nHF=0, nBad1=0, nBad2=0, nBad=0;
   std::vector<double> encontHB(140, 0.);
   std::vector<double> encontHE(140, 0.);
   std::vector<double> encontHF(140, 0.);
   std::vector<double> encontHO(140, 0.);
   double entotHB = 0, entotHE = 0, entotHF = 0, entotHO = 0; 
 
   double HBEneMap[ieta_bins_HB][iphi_bins];
   double HEEneMap[ieta_bins_HE][iphi_bins];
   double HOEneMap[ieta_bins_HO][iphi_bins];
   double HFEneMap[ieta_bins_HF][iphi_bins];
  
   //Works in ieta_min_Hx is < 0
   int eta_offset_HB = -(int)ieta_min_HB;
   int eta_offset_HE = -(int)ieta_min_HE;
   int eta_offset_HO = -(int)ieta_min_HO;
   int eta_offset_HF = -(int)ieta_min_HF;
 
   for(int i = 0; i < ieta_bins_HB; i++){
      for(int j = 0; j < iphi_bins; j++){
        HBEneMap[i][j] = 0.;
      }
   }
 
   for(int i = 0; i < ieta_bins_HE; i++){
      for(int j = 0; j < iphi_bins; j++){
        HEEneMap[i][j] = 0.;
      }
   }
 
   for(int i = 0; i < ieta_bins_HO; i++){
      for(int j = 0; j < iphi_bins; j++){
        HOEneMap[i][j] = 0.;
      }
   }
 
   for(int i = 0; i < ieta_bins_HF; i++){
      for(int j = 0; j < iphi_bins; j++){
        HFEneMap[i][j] = 0.;
      }
   }
 
   for (int i=0; i<nHit; i++) {
     double energy    = hits[i].energy();
     double log10en   = log10(energy);
     int log10i       = int( (log10en+10.)*10. );
     double time      = hits[i].time();
     unsigned int id_ = hits[i].id();
     int det, subdet, depth, eta, phi;
     HcalDetId hid;
     if (testNumber_) hid = HcalHitRelabeller::relabel(id_,hcons);
     else hid = HcalDetId(id_);
     det      = hid.det();
     subdet   = hid.subdet();
     depth    = hid.depth();
     eta      = hid.ieta();
     phi      = hid.iphi();
 
     LogDebug("HcalSim") << "Hit[" << i << "] ID " << std::hex << id_ 
             << std::dec << " Det " << det << " Sub " 
             << subdet << " depth " << depth << " Eta " << eta
             << " Phi " << phi << " E " << energy << " 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)) nHF++;
       else    { nBad++;  nBad2++;}
     } else    { nBad++;  nBad1++;}
 
       meDetectHit_->Fill(double(det));
       if (det ==  4) {
     meSubdetHit_->Fill(double(subdet));
     meDepthHit_->Fill(double(depth));
     meEtaHit_->Fill(double(eta));
     meEtaPhiHit_->Fill(double(eta),double(phi));
     meEtaPhiHitDepth_[depth-1]->Fill(double(eta),double(phi));
 
     //We will group the phi plots by HB,HO and HE,HF since these groups share similar segmentation schemes
     if      (subdet == static_cast<int>(HcalBarrel))  mePhiHit_->Fill(double(phi));
     else if (subdet == static_cast<int>(HcalEndcap))  mePhiHitb_->Fill(double(phi));
     else if (subdet == static_cast<int>(HcalOuter))   mePhiHit_->Fill(double(phi));
     else if (subdet == static_cast<int>(HcalForward)) mePhiHitb_->Fill(double(phi));    
     
 
     //KC: HF energy is in photoelectrons rather than eV, so it will not be included in total HCal energy
     if (subdet != static_cast<int>(HcalForward)){
         meEnergyHit_->Fill(energy);
 
         //Since the HF energy is a different scale it does not make sense to include it in the Energy Weighted Plot
         meTimeWHit_->Fill(double(time),energy);
     }
     meTimeHit_->Fill(time);
     
     if      (subdet == static_cast<int>(HcalBarrel)) {
       meHBDepHit_->Fill(double(depth));
       meHBEtaHit_->Fill(double(eta));
       meHBPhiHit_->Fill(double(phi));
       meHBEneHit_->Fill(energy);
       meHBEneHit2_->Fill(energy);
       meHBTimHit_->Fill(time);
       meHBL10Ene_->Fill(log10en);
       if( log10i >=0 && log10i < 140 ) encontHB[log10i] += energy;
       entotHB += energy;
     
           HBEneMap[eta + eta_offset_HB][phi-1] += energy;
 
     } else if (subdet == static_cast<int>(HcalEndcap)) {
       meHEDepHit_->Fill(double(depth));
       meHEEtaHit_->Fill(double(eta));
       meHEPhiHit_->Fill(double(phi));
 
       bool isHEP17 = (phi>=63)&&(phi<=66)&&(eta>0); 
       if(hep17_){
         if(!isHEP17){
           meHEEneHit_->Fill(energy);
           meHEEneHit2_->Fill(energy);
         } else {
           meHEP17EneHit_->Fill(energy);
           meHEP17EneHit2_->Fill(energy);
         }
       } else {
         meHEEneHit_->Fill(energy);
         meHEEneHit2_->Fill(energy);
       }
 
       meHETimHit_->Fill(time);
       meHEL10Ene_->Fill(log10en);
       if( log10i >=0 && log10i < 140 ) encontHE[log10i] += energy;
       entotHE += energy;
     
           HEEneMap[eta + eta_offset_HE][phi-1] += energy;
 
     } else if (subdet == static_cast<int>(HcalOuter)) {
       meHODepHit_->Fill(double(depth));
       meHOEtaHit_->Fill(double(eta));
       meHOPhiHit_->Fill(double(phi));
       meHOEneHit_->Fill(energy);
       meHOEneHit2_->Fill(energy);
       meHOTimHit_->Fill(time);
       meHOL10Ene_->Fill(log10en);
       if( log10i >=0 && log10i < 140 ) encontHO[log10i] += energy;
       entotHO += energy;
     
           HOEneMap[eta + eta_offset_HO][phi-1] += energy;
 
     } else if (subdet == static_cast<int>(HcalForward)) {
       meHFDepHit_->Fill(double(depth));
       meHFDepHitw_->Fill(double(depth),energy);
       meHFEtaHit_->Fill(double(eta));
       meHFPhiHit_->Fill(double(phi));
       meHFEneHit_->Fill(energy);
       meHFEneHit2_->Fill(energy);
       meHFTimHit_->Fill(time);
       meHFL10Ene_->Fill(log10en);
       if( log10i >=0 && log10i < 140 ) encontHF[log10i] += energy;
       entotHF += energy;
     
           HFEneMap[eta + eta_offset_HF][phi-1] += energy;
 
     }
       }
   }
   if( entotHB != 0 ) for( int i=0; i<140; i++ ) meHBL10EneP_->Fill( -10.+(float(i)+0.5)/10., encontHB[i]/entotHB );
   if( entotHE != 0 ) for( int i=0; i<140; i++ ) meHEL10EneP_->Fill( -10.+(float(i)+0.5)/10., encontHE[i]/entotHE );
   if( entotHF != 0 ) for( int i=0; i<140; i++ ) meHFL10EneP_->Fill( -10.+(float(i)+0.5)/10., encontHF[i]/entotHF );
   if( entotHO != 0 ) for( int i=0; i<140; i++ ) meHOL10EneP_->Fill( -10.+(float(i)+0.5)/10., encontHO[i]/entotHO );
 
     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));
 
   for(int i = 0; i < ieta_bins_HB; i++){
      for(int j = 0; j < iphi_bins; j++){
 
         if(HBEneMap[i][j] != 0){
            meHBEneSum_->Fill(HBEneMap[i][j]);
            meHBEneSum_vs_ieta_->Fill((i - eta_offset_HB), HBEneMap[i][j]);
            meHBEneMap_->Fill((i - eta_offset_HB), j+1, HBEneMap[i][j]);
         }
 
      }
   }
   
   for(int i = 0; i < ieta_bins_HE; i++){
      for(int j = 0; j < iphi_bins; j++){
 
         if(HEEneMap[i][j] != 0){
            meHEEneSum_->Fill(HEEneMap[i][j]);
            meHEEneSum_vs_ieta_->Fill((i - eta_offset_HE), HEEneMap[i][j]);
            meHEEneMap_->Fill((i - eta_offset_HE), j+1, HEEneMap[i][j]);
         }
 
      }
   }
   
   for(int i = 0; i < ieta_bins_HO; i++){
      for(int j = 0; j < iphi_bins; j++){
 
         if(HOEneMap[i][j] != 0){
            meHOEneSum_->Fill(HOEneMap[i][j]);
            meHOEneSum_vs_ieta_->Fill((i - eta_offset_HO), HOEneMap[i][j]);
            meHOEneMap_->Fill((i - eta_offset_HO), j+1, HOEneMap[i][j]);
         }
 
      }
   }
   
   for(int i = 0; i < ieta_bins_HF; i++){
      for(int j = 0; j < iphi_bins; j++){
 
         if(HFEneMap[i][j] != 0){
            meHFEneSum_->Fill(HFEneMap[i][j]);
            meHFEneSum_vs_ieta_->Fill((i - eta_offset_HF), HFEneMap[i][j]);
            meHFEneMap_->Fill((i - eta_offset_HF), j+1, HFEneMap[i][j]);
         }
 
      }
   }
   
   LogDebug("HcalSim") << "HcalSimHitStudy::analyzeHits: HB " << nHB 
               << " HE " << nHE << " HO " << nHO << " HF " << nHF 
               << " Bad " << nBad << " All " << nHit;
 
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 27 of file HcalSimHitStudy.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), checkHit_, egammaForCoreTracking_cff::depth, edm::EventSetup::get(), HcalDDDRecConstants::getEtaRange(), HcalDDDRecConstants::getMaxDepth(), HcalDDDRecConstants::getNPhi(), hcons, ieta_bins_HB, ieta_bins_HE, ieta_bins_HF, ieta_bins_HO, ieta_max_HB, ieta_max_HE, ieta_max_HF, ieta_max_HO, ieta_min_HB, ieta_min_HE, ieta_min_HF, ieta_min_HO, hcalTTPDigis_cfi::iEtaMax, createfilelist::int, iphi_bins, iphi_max, iphi_min, maxDepth_, maxDepthHB_, maxDepthHE_, maxDepthHF_, maxDepthHO_, meAllNHit_, meBadDetHit_, meBadIdHit_, meBadSubHit_, meDepthHit_, meDetectHit_, meEnergyHit_, meEtaHit_, meEtaPhiHit_, meEtaPhiHitDepth_, meHBDepHit_, meHBEneHit2_, meHBEneHit_, meHBEneMap_, meHBEneSum_, meHBEneSum_vs_ieta_, meHBEtaHit_, meHBL10Ene_, meHBL10EneP_, meHBNHit_, meHBPhiHit_, meHBTimHit_, meHEDepHit_, meHEEneHit2_, meHEEneHit_, meHEEneMap_, meHEEneSum_, meHEEneSum_vs_ieta_, meHEEtaHit_, meHEL10Ene_, meHEL10EneP_, meHENHit_, meHEP17EneHit2_, meHEP17EneHit_, meHEPhiHit_, meHETimHit_, meHFDepHit_, meHFDepHitw_, meHFEneHit2_, meHFEneHit_, meHFEneMap_, meHFEneSum_, meHFEneSum_vs_ieta_, meHFEtaHit_, meHFL10Ene_, meHFL10EneP_, meHFNHit_, meHFPhiHit_, meHFTimHit_, meHODepHit_, meHOEneHit2_, meHOEneHit_, meHOEneMap_, meHOEneSum_, meHOEneSum_vs_ieta_, meHOEtaHit_, meHOL10Ene_, meHOL10EneP_, meHONHit_, meHOPhiHit_, meHOTimHit_, mePhiHit_, mePhiHitb_, meSubdetHit_, meTimeHit_, meTimeWHit_, and DQMStore::IBooker::setCurrentFolder().

 {
   edm::ESHandle<HcalDDDRecConstants> pHRNDC;
   es.get<HcalRecNumberingRecord>().get( pHRNDC );
   hcons = &(*pHRNDC);
   maxDepthHB_ = hcons->getMaxDepth(0);
   maxDepthHE_ = hcons->getMaxDepth(1);
   maxDepthHF_ = hcons->getMaxDepth(2);
   maxDepthHO_ = hcons->getMaxDepth(3);
   maxDepth_ = (maxDepthHB_ > maxDepthHE_ ? maxDepthHB_ : maxDepthHE_);
   maxDepth_ = (maxDepth_   > maxDepthHF_ ? maxDepth_   : maxDepthHF_);
   maxDepth_ = (maxDepth_   > maxDepthHO_ ? maxDepth_   : maxDepthHO_);
 
   //Get Phi segmentation from geometry, use the max phi number so that all iphi values are included.
 
   int NphiMax = hcons->getNPhi(0);
 
   NphiMax = (hcons->getNPhi(1) > NphiMax ? hcons->getNPhi(1) : NphiMax);
   NphiMax = (hcons->getNPhi(2) > NphiMax ? hcons->getNPhi(2) : NphiMax);
   NphiMax = (hcons->getNPhi(3) > NphiMax ? hcons->getNPhi(3) : NphiMax);
 
   //Center the iphi bins on the integers
   iphi_min = 0.5;
   iphi_max = NphiMax + 0.5;
   iphi_bins = (int) (iphi_max - iphi_min);
 
   int iEtaHBMax = hcons->getEtaRange(0).second;
   int iEtaHEMax = hcons->getEtaRange(1).second;
   int iEtaHFMax = hcons->getEtaRange(2).second;
   int iEtaHOMax = hcons->getEtaRange(3).second;
 
   //Retain classic behavior, all plots have same ieta range.
   //Comment out code to allow each subdetector to have its on range
 
   int iEtaMax = (iEtaHBMax > iEtaHEMax ? iEtaHBMax : iEtaHEMax);
   iEtaMax = (iEtaMax > iEtaHFMax ? iEtaMax : iEtaHFMax);
   iEtaMax = (iEtaMax > iEtaHOMax ? iEtaMax : iEtaHOMax);
 
   iEtaHBMax = iEtaMax;
   iEtaHEMax = iEtaMax;
   iEtaHFMax = iEtaMax;
   iEtaHOMax = iEtaMax;
 
   //Give an empty bin around the subdet ieta range to make it clear that all ieta rings have been included
   ieta_min_HB = -iEtaHBMax - 1.5;
   ieta_max_HB = iEtaHBMax + 1.5;
   ieta_bins_HB = (int) (ieta_max_HB - ieta_min_HB);
 
   ieta_min_HE = -iEtaHEMax - 1.5;
   ieta_max_HE = iEtaHEMax + 1.5;
   ieta_bins_HE = (int) (ieta_max_HE - ieta_min_HE);
 
   ieta_min_HF = -iEtaHFMax - 1.5;
   ieta_max_HF = iEtaHFMax + 1.5;
   ieta_bins_HF = (int) (ieta_max_HF - ieta_min_HF);
 
   ieta_min_HO = -iEtaHOMax - 1.5;
   ieta_max_HO = iEtaHOMax + 1.5;
   ieta_bins_HO = (int) (ieta_max_HO - ieta_min_HO);
 
   Char_t hname[100];
   Char_t htitle[100];
 
     ib.setCurrentFolder("HcalHitsV/HcalSimHitsTask");
 
     //Histograms for Hits
     if (checkHit_) {
       meAllNHit_  = ib.book1D("Hit01","Number of Hits in HCal",20000,0.,20000.);
       meBadDetHit_= ib.book1D("Hit02","Hits with wrong Det",   100,0.,100.);
       meBadSubHit_= ib.book1D("Hit03","Hits with wrong Subdet",100,0.,100.);
       meBadIdHit_ = ib.book1D("Hit04","Hits with wrong ID",    100,0.,100.);
       meHBNHit_   = ib.book1D("Hit05","Number of Hits in HB",20000,0.,20000.);
       meHENHit_   = ib.book1D("Hit06","Number of Hits in HE",10000,0.,10000.);
       meHONHit_   = ib.book1D("Hit07","Number of Hits in HO",10000,0.,10000.);
       meHFNHit_   = ib.book1D("Hit08","Number of Hits in HF",10000,0.,10000.);
       meDetectHit_= ib.book1D("Hit09","Detector ID",           50,0.,50.);
       meSubdetHit_= ib.book1D("Hit10","Subdetectors in HCal",  50,0.,50.);
       meDepthHit_ = ib.book1D("Hit11","Depths in HCal",        20,0.,20.);
       meEtaHit_   = ib.book1D("Hit12","Eta in HCal",          ieta_bins_HF,ieta_min_HF,ieta_max_HF);
       meEtaPhiHit_   = ib.book2D("Hit12b","Eta-phi in HCal",          ieta_bins_HF,ieta_min_HF,ieta_max_HF,iphi_bins,iphi_min,iphi_max);
       for (int depth = 1; depth <= maxDepth_; depth++) {
     sprintf(hname,"Hit12bd%d", depth);
     sprintf(htitle,"Eta-phi in HCal d%d", depth);
     meEtaPhiHitDepth_.push_back( ib.book2D(hname, htitle, 
                            ieta_bins_HF,ieta_min_HF,ieta_max_HF,
                            iphi_bins,iphi_min,iphi_max) );
       }
       //KC: There are different phi segmentation schemes, this plot uses wider bins to represent the most sparse segmentation
       mePhiHit_   = ib.book1D("Hit13","Phi in HCal (HB,HO)",  iphi_bins,iphi_min,iphi_max);
       mePhiHitb_  = ib.book1D("Hit13b","Phi in HCal (HE,HF)", iphi_bins,iphi_min,iphi_max);
       meEnergyHit_= ib.book1D("Hit14","Energy in HCal",       2000,0.,20.);
       meTimeHit_  = ib.book1D("Hit15","Time in HCal",         528,0.,528.);
       meTimeWHit_ = ib.book1D("Hit16","Time in HCal (E wtd)", 528,0.,528.);
       meHBDepHit_ = ib.book1D("Hit17","Depths in HB",          20,0.,20.);
       meHEDepHit_ = ib.book1D("Hit18","Depths in HE",          20,0.,20.);
       meHODepHit_ = ib.book1D("Hit19","Depths in HO",          20,0.,20.);
       meHFDepHit_ = ib.book1D("Hit20","Depths in HF",          20,0.,20.);
       meHFDepHitw_ = ib.book1D("Hit20b","Depths in HF (p.e. weighted)",          20,0.,20.);
       meHBEtaHit_ = ib.book1D("Hit21","Eta in HB",            ieta_bins_HB,ieta_min_HB,ieta_max_HB);
       meHEEtaHit_ = ib.book1D("Hit22","Eta in HE",            ieta_bins_HE,ieta_min_HE,ieta_max_HE);
       meHOEtaHit_ = ib.book1D("Hit23","Eta in HO",            ieta_bins_HO,ieta_min_HO,ieta_max_HO);
       meHFEtaHit_ = ib.book1D("Hit24","Eta in HF",            ieta_bins_HF,ieta_min_HF,ieta_max_HF);
       meHBPhiHit_ = ib.book1D("Hit25","Phi in HB",            iphi_bins,iphi_min,iphi_max);
       meHEPhiHit_ = ib.book1D("Hit26","Phi in HE",            iphi_bins,iphi_min,iphi_max); 
       meHOPhiHit_ = ib.book1D("Hit27","Phi in HO",            iphi_bins,iphi_min,iphi_max); 
       meHFPhiHit_ = ib.book1D("Hit28","Phi in HF",            iphi_bins,iphi_min,iphi_max); 
       meHBEneHit_ = ib.book1D("Hit29","Energy in HB",         2000,0.,20.);
       meHEEneHit_ = ib.book1D("Hit30","Energy in HE",         500,0.,5.);
       meHEP17EneHit_ = ib.book1D("Hit30b","Energy in HEP17",     500,0.,5.);
       meHOEneHit_ = ib.book1D("Hit31","Energy in HO",         500,0.,5.);
       meHFEneHit_ = ib.book1D("Hit32","Energy in HF",         1001,-0.5,1000.5);
 
       //HxEneMap, HxEneSum, HxEneSum_vs_ieta plot the sum of the simhits energy within a single ieta-iphi tower.
 
       meHBEneMap_ = ib.book2D("HBEneMap","HBEneMap", ieta_bins_HB, ieta_min_HB, ieta_max_HB, iphi_bins, iphi_min, iphi_max);
       meHEEneMap_ = ib.book2D("HEEneMap","HEEneMap", ieta_bins_HE, ieta_min_HE, ieta_max_HE, iphi_bins, iphi_min, iphi_max);
       meHOEneMap_ = ib.book2D("HOEneMap","HOEneMap", ieta_bins_HO, ieta_min_HO, ieta_max_HO, iphi_bins, iphi_min, iphi_max);
       meHFEneMap_ = ib.book2D("HFEneMap","HFEneMap", ieta_bins_HF, ieta_min_HF, ieta_max_HF, iphi_bins, iphi_min, iphi_max);
 
       meHBEneSum_ = ib.book1D("HBEneSum","HBEneSum", 2000, 0., 20.);
       meHEEneSum_ = ib.book1D("HEEneSum","HEEneSum", 500, 0., 5.);
       meHOEneSum_ = ib.book1D("HOEneSum","HOEneSum", 500, 0., 5.);
       meHFEneSum_ = ib.book1D("HFEneSum","HFEneSum", 1001, -0.5, 1000.5);
 
       meHBEneSum_vs_ieta_ = ib.bookProfile("HBEneSum_vs_ieta","HBEneSum_vs_ieta", ieta_bins_HB, ieta_min_HB, ieta_max_HB, 2011, -10.5, 2000.5, " ");
       meHEEneSum_vs_ieta_ = ib.bookProfile("HEEneSum_vs_ieta","HEEneSum_vs_ieta", ieta_bins_HE, ieta_min_HE, ieta_max_HE, 2011, -10.5, 2000.5, " ");
       meHOEneSum_vs_ieta_ = ib.bookProfile("HOEneSum_vs_ieta","HOEneSum_vs_ieta", ieta_bins_HO, ieta_min_HO, ieta_max_HO, 2011, -10.5, 2000.5, " ");
       meHFEneSum_vs_ieta_ = ib.bookProfile("HFEneSum_vs_ieta","HFEneSum_vs_ieta", ieta_bins_HF, ieta_min_HF, ieta_max_HF, 2011, -10.5, 2000.5, " ");
 
       meHBTimHit_ = ib.book1D("Hit33","Time in HB",           528,0.,528.);
       meHETimHit_ = ib.book1D("Hit34","Time in HE",           528,0.,528.);
       meHOTimHit_ = ib.book1D("Hit35","Time in HO",           528,0.,528.);
       meHFTimHit_ = ib.book1D("Hit36","Time in HF",           528,0.,528.);
       //These are the zoomed in energy ranges
       meHBEneHit2_ = ib.book1D("Hit37","Energy in HB 2",         100,0.,0.0001);
       meHEEneHit2_ = ib.book1D("Hit38","Energy in HE 2",         100,0.,0.0001);
       meHEP17EneHit2_ = ib.book1D("Hit38b","Energy in HEP17 2",     100,0.,0.0001);
       meHOEneHit2_ = ib.book1D("Hit39","Energy in HO 2",         100,0.,0.0001);
       meHFEneHit2_ = ib.book1D("Hit40","Energy in HF 2",         100,0.5,100.5);
       meHBL10Ene_ = ib.book1D("Hit41","Log10Energy in HB", 140, -10., 4. );
       meHEL10Ene_ = ib.book1D("Hit42","Log10Energy in HE", 140, -10., 4. );
       meHFL10Ene_ = ib.book1D("Hit43","Log10Energy in HF", 50, -1., 4. );
       meHOL10Ene_ = ib.book1D("Hit44","Log10Energy in HO", 140, -10., 4. );
       meHBL10EneP_ = ib.bookProfile("Hit45","Log10Energy in HB vs Hit contribution", 140, -10., 4., 100, 0., 1. );
       meHEL10EneP_ = ib.bookProfile("Hit46","Log10Energy in HE vs Hit contribution", 140, -10., 4., 100, 0., 1. );
       meHFL10EneP_ = ib.bookProfile("Hit47","Log10Energy in HF vs Hit contribution", 140, -10., 4., 100, 0., 1. );
       meHOL10EneP_ = ib.bookProfile("Hit48","Log10Energy in HO vs Hit contribution", 140, -10., 4., 100, 0., 1. );
     }
 
 }