CaloTowersAnalyzer Class Reference

#include <CaloTowersAnalyzer.h>

Inheritance diagram for CaloTowersAnalyzer:

Public Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &c) override
void	beginJob () override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	CaloTowersAnalyzer (edm::ParameterSet const &conf)
void	dqmBeginRun (const edm::Run &run, const edm::EventSetup &c) override
void	endJob () override
	~CaloTowersAnalyzer () 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 Types
typedef math::RhoEtaPhiVector	Vector

Private Member Functions
double	dR (double eta1, double phi1, double eta2, double phi2)

Private Attributes
MonitorElement *	emean_vs_ieta_E
MonitorElement *	emean_vs_ieta_E1
MonitorElement *	emean_vs_ieta_EH
MonitorElement *	emean_vs_ieta_EH1
MonitorElement *	emean_vs_ieta_H
MonitorElement *	emean_vs_ieta_H1
MonitorElement *	emEnergyTiming_HB
MonitorElement *	emEnergyTiming_HE
MonitorElement *	emEnergyTiming_HF
MonitorElement *	emEnergyTiming_High_HB
MonitorElement *	emEnergyTiming_Low_HB
MonitorElement *	emEnergyTiming_Low_HE
MonitorElement *	emEnergyTiming_profile_HB
MonitorElement *	emEnergyTiming_profile_HE
MonitorElement *	emEnergyTiming_profile_HF
MonitorElement *	emEnergyTiming_profile_High_HB
MonitorElement *	emEnergyTiming_profile_Low_HB
MonitorElement *	emEnergyTiming_profile_Low_HE
MonitorElement *	emTiming_HB
MonitorElement *	emTiming_HE
MonitorElement *	emTiming_HF
double	etaMax [3]
double	etaMin [3]
MonitorElement *	hadEnergyTiming_HB
MonitorElement *	hadEnergyTiming_HE
MonitorElement *	hadEnergyTiming_HF
MonitorElement *	hadEnergyTiming_High_HB
MonitorElement *	hadEnergyTiming_Low_HB
MonitorElement *	hadEnergyTiming_Low_HE
MonitorElement *	hadEnergyTiming_Low_HF
MonitorElement *	hadEnergyTiming_profile_HB
MonitorElement *	hadEnergyTiming_profile_HE
MonitorElement *	hadEnergyTiming_profile_HF
MonitorElement *	hadEnergyTiming_profile_High_HB
MonitorElement *	hadEnergyTiming_profile_Low_HB
MonitorElement *	hadEnergyTiming_profile_Low_HE
MonitorElement *	hadEnergyTiming_profile_Low_HF
MonitorElement *	hadTiming_HB
MonitorElement *	hadTiming_HE
MonitorElement *	hadTiming_HF
std::string	hcalselector_
const HcalDDDRecConstants *	hcons
int	ieta_bins_
float	ieta_max_
float	ieta_min_
int	imc
int	iphi_bins_
float	iphi_max_
float	iphi_min_
int	isub
MonitorElement *	mapEnergy_E
MonitorElement *	mapEnergy_EH
MonitorElement *	mapEnergy_H
MonitorElement *	mapEnergy_HB
MonitorElement *	mapEnergy_HE
MonitorElement *	mapEnergy_HF
MonitorElement *	mapEnergy_N
MonitorElement *	mapEnergyEcal_HB
MonitorElement *	mapEnergyEcal_HE
MonitorElement *	mapEnergyEcal_HF
MonitorElement *	mapEnergyHcal_HB
MonitorElement *	mapEnergyHcal_HE
MonitorElement *	mapEnergyHcal_HF
std::string	mc_
MonitorElement *	meEnergyEcal_HB
MonitorElement *	meEnergyEcal_HE
MonitorElement *	meEnergyEcal_HF
MonitorElement *	meEnergyEcalTower_HB
MonitorElement *	meEnergyEcalTower_HE
MonitorElement *	meEnergyEcalTower_HF
MonitorElement *	meEnergyHcal_HB
MonitorElement *	meEnergyHcal_HE
MonitorElement *	meEnergyHcal_HF
MonitorElement *	meEnergyHcalTower_HB
MonitorElement *	meEnergyHcalTower_HE
MonitorElement *	meEnergyHcalTower_HF
MonitorElement *	meEnergyHcalvsEcal_HB
MonitorElement *	meEnergyHcalvsEcal_HE
MonitorElement *	meEnergyHcalvsEcal_HF
MonitorElement *	meEnergyHO_HB
MonitorElement *	meEnergyHO_HE
MonitorElement *	meEnergyHO_HF
MonitorElement *	meIphiCaloTower_HFM
MonitorElement *	meIphiCaloTower_HFP
MonitorElement *	meIphiHcalTower_HBM
MonitorElement *	meIphiHcalTower_HBP
MonitorElement *	meIphiHcalTower_HEM
MonitorElement *	meIphiHcalTower_HEP
MonitorElement *	meNumFiredTowers_HB
MonitorElement *	meNumFiredTowers_HE
MonitorElement *	meNumFiredTowers_HF
MonitorElement *	MET_HB
MonitorElement *	MET_HE
MonitorElement *	MET_HF
MonitorElement *	meTotEnergy_HB
MonitorElement *	meTotEnergy_HE
MonitorElement *	meTotEnergy_HF
int	nevent
MonitorElement *	Ntowers_vs_ieta
MonitorElement *	numBadCellsEcal_EB
MonitorElement *	numBadCellsEcal_EE
MonitorElement *	numBadCellsHcal_HB
MonitorElement *	numBadCellsHcal_HE
MonitorElement *	numBadCellsHcal_HF
MonitorElement *	numPrbCellsEcal_EB
MonitorElement *	numPrbCellsEcal_EE
MonitorElement *	numPrbCellsHcal_HB
MonitorElement *	numPrbCellsHcal_HE
MonitorElement *	numPrbCellsHcal_HF
MonitorElement *	numRcvCellsEcal_EB
MonitorElement *	numRcvCellsEcal_EE
MonitorElement *	numRcvCellsHcal_HB
MonitorElement *	numRcvCellsHcal_HE
MonitorElement *	numRcvCellsHcal_HF
MonitorElement *	occupancy_map
MonitorElement *	occupancy_vs_ieta
std::string	outputFile_
MonitorElement *	phiMET_HB
MonitorElement *	phiMET_HE
MonitorElement *	phiMET_HF
MonitorElement *	SET_HB
MonitorElement *	SET_HE
MonitorElement *	SET_HF
edm::EDGetTokenT< CaloTowerCollection >	tok_towers_
bool	useAllHistos_

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 41 of file CaloTowersAnalyzer.h.

Member Typedef Documentation

typedef math::RhoEtaPhiVector CaloTowersAnalyzer::Vector

private

Definition at line 60 of file CaloTowersAnalyzer.h.

Constructor & Destructor Documentation

CaloTowersAnalyzer::CaloTowersAnalyzer

(

edm::ParameterSet const &

conf

)

Definition at line 6 of file CaloTowersAnalyzer.cc.

References edm::ParameterSet::getUntrackedParameter(), hcalselector_, outputFile_, AlCaHLTBitMon_QueryRunRegistry::string, tok_towers_, and useAllHistos_.

                                                                 {

  tok_towers_ = consumes<CaloTowerCollection>(conf.getUntrackedParameter<edm::InputTag>("CaloTowerCollectionLabel"));
  // DQM ROOT output
  
  outputFile_ = conf.getUntrackedParameter<std::string>("outputFile", "myfile.root");

  hcalselector_ = conf.getUntrackedParameter<std::string>("hcalselector", "all");

  useAllHistos_ = conf.getUntrackedParameter<bool>("useAllHistos", false);
}

CaloTowersAnalyzer::~CaloTowersAnalyzer ( )

override

Definition at line 446 of file CaloTowersAnalyzer.cc.

                                        {
  
}

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 455 of file CaloTowersAnalyzer.cc.

References edm::SortedCollection< T, SORT >::begin(), DetId::det(), DetId::Ecal, EcalBarrel, EcalEndcap, emEnergyTiming_HB, emEnergyTiming_HE, emEnergyTiming_HF, emEnergyTiming_High_HB, emEnergyTiming_Low_HB, emEnergyTiming_Low_HE, emEnergyTiming_profile_HB, emEnergyTiming_profile_HE, emEnergyTiming_profile_HF, emEnergyTiming_profile_High_HB, emEnergyTiming_profile_Low_HB, emEnergyTiming_profile_Low_HE, emTiming_HB, emTiming_HE, emTiming_HF, edm::SortedCollection< T, SORT >::end(), etaMax, etaMin, MonitorElement::Fill(), HcalObjRepresent::Fill(), hadEnergyTiming_HB, hadEnergyTiming_HE, hadEnergyTiming_HF, hadEnergyTiming_High_HB, hadEnergyTiming_Low_HB, hadEnergyTiming_Low_HE, hadEnergyTiming_Low_HF, hadEnergyTiming_profile_HB, hadEnergyTiming_profile_HE, hadEnergyTiming_profile_HF, hadEnergyTiming_profile_High_HB, hadEnergyTiming_profile_Low_HB, hadEnergyTiming_profile_Low_HE, hadEnergyTiming_profile_Low_HF, hadTiming_HB, hadTiming_HE, hadTiming_HF, DetId::Hcal, HcalBarrel, HcalEndcap, HcalForward, CaloTowerDetId::ieta(), isub, mapEnergy_HB, mapEnergy_HE, mapEnergy_HF, mapEnergy_N, mapEnergyEcal_HB, mapEnergyEcal_HE, mapEnergyEcal_HF, mapEnergyHcal_HB, mapEnergyHcal_HE, mapEnergyHcal_HF, meEnergyEcal_HB, meEnergyEcal_HE, meEnergyEcal_HF, meEnergyEcalTower_HB, meEnergyEcalTower_HE, meEnergyEcalTower_HF, meEnergyHcal_HB, meEnergyHcal_HE, meEnergyHcal_HF, meEnergyHcalTower_HB, meEnergyHcalTower_HE, meEnergyHcalTower_HF, meEnergyHcalvsEcal_HB, meEnergyHcalvsEcal_HE, meEnergyHcalvsEcal_HF, meEnergyHO_HB, meEnergyHO_HE, meEnergyHO_HF, meIphiCaloTower_HFM, meIphiCaloTower_HFP, meIphiHcalTower_HBM, meIphiHcalTower_HBP, meIphiHcalTower_HEM, meIphiHcalTower_HEP, meNumFiredTowers_HB, meNumFiredTowers_HE, meNumFiredTowers_HF, RazorAnalyzer::met, MET_HB, MET_HE, MET_HF, meTotEnergy_HB, meTotEnergy_HE, meTotEnergy_HF, nevent, Ntowers_vs_ieta, numBadCellsEcal_EB, numBadCellsEcal_EE, numBadCellsHcal_HB, numBadCellsHcal_HE, numBadCellsHcal_HF, numPrbCellsEcal_EB, numPrbCellsEcal_EE, numPrbCellsHcal_HB, numPrbCellsHcal_HE, numPrbCellsHcal_HF, numRcvCellsEcal_EB, numRcvCellsEcal_EE, numRcvCellsHcal_HB, numRcvCellsHcal_HE, numRcvCellsHcal_HF, occupancy_map, phiMET_HB, phiMET_HE, phiMET_HF, SET_HB, SET_HE, SET_HF, mathSSE::sqrt(), HcalDetId::subdet(), DetId::subdetId(), tok_towers_, and useAllHistos_.

                                                                              {

  nevent++;

  edm::Handle<CaloTowerCollection> towers;
  event.getByToken(tok_towers_,towers);
  CaloTowerCollection::const_iterator cal;

  double met;
  double phimet;

  // ieta scan 
  //double partR  = 0.3;
  //double Rmin   = 9999.;
  //double Econe  = 0.;
  //double Hcone  = 0.;
  //double Ee1    = 0.;
  //double Eh1    = 0.;
  double ieta_MC = 9999;
  double iphi_MC = 9999;
 
  // HB   
  double sumEnergyHcal_HB = 0.;
  double sumEnergyEcal_HB = 0.;
  double sumEnergyHO_HB   = 0.;
  Int_t numFiredTowers_HB = 0;
  double metx_HB   =  0.;
  double mety_HB   =  0.;
  double metz_HB   =  0.;
  double sEt_HB    =  0.;
  // HE   
  double sumEnergyHcal_HE = 0.;
  double sumEnergyEcal_HE = 0.;
  double sumEnergyHO_HE   = 0.;
  Int_t numFiredTowers_HE = 0;
  double metx_HE   =  0.;
  double mety_HE   =  0.;
  double metz_HE   =  0.;
  double sEt_HE    =  0.;
  // HF   
  double sumEnergyHcal_HF = 0.;
  double sumEnergyEcal_HF = 0.;
  double sumEnergyHO_HF   = 0.;
  Int_t numFiredTowers_HF = 0;
  double metx_HF   =  0.;
  double mety_HF   =  0.;
  double metz_HF   =  0.;
  double sEt_HF    =  0.;

  for ( cal = towers->begin(); cal != towers->end(); ++cal ) {
    
    double eE     = cal->emEnergy();
    double eH     = cal->hadEnergy();
    double eHO    = cal->outerEnergy();
    double etaT   = cal->eta();
    //    double phiT   = cal->phi();
    double en     = cal->energy();
    double etT    = cal->et();
    double had_tm = cal->hcalTime();
    double em_tm  = cal->ecalTime();

    int numBadEcalCells = cal->numBadEcalCells();
    int numRcvEcalCells = cal->numRecoveredEcalCells();
    int numPrbEcalCells = cal->numProblematicEcalCells();

    int numBadHcalCells = cal->numBadHcalCells();
    int numRcvHcalCells = cal->numRecoveredHcalCells();
    int numPrbHcalCells = cal->numProblematicHcalCells(); 

    math::RhoEtaPhiVector mom(cal->et(), cal->eta(), cal->phi());
                  //  Vector mom  = cal->momentum(); 
  
    // cell properties    
    CaloTowerDetId idT = cal->id();
    int ieta = idT.ieta();
    int iphi = idT.iphi();

    // ecal:  0 EcalBarrel  1 EcalEndcap
    // hcal:  0 hcalBarrel  1 HcalEndcap  2 HcalForward
    std::vector<int> inEcals(2), inHcals(3);
    unsigned int constitSize = cal->constituentsSize();
    for(unsigned int ic=0; ic<constitSize; ic++){
       DetId detId = cal->constituent(ic);
       if(detId.det() == DetId::Ecal){
          if( detId.subdetId()==EcalBarrel ) inEcals[0] =1;
          else if( detId.subdetId()==EcalEndcap ) inEcals[1] =1;
       }
       if(detId.det() == DetId::Hcal){
          if( HcalDetId(detId).subdet()==HcalBarrel ) inHcals[0] =1;
          else if( HcalDetId(detId).subdet()==HcalEndcap ) inHcals[1] =1;
          else if( HcalDetId(detId).subdet()==HcalForward ) inHcals[2] =1;
       }
    }
    //All cell histos are used
    if( inEcals[0] ) {
       numBadCellsEcal_EB->Fill(numBadEcalCells);
       numRcvCellsEcal_EB->Fill(numRcvEcalCells);
       numPrbCellsEcal_EB->Fill(numPrbEcalCells);
    }
    if( inEcals[1] ) {
       numBadCellsEcal_EE->Fill(numBadEcalCells);
       numRcvCellsEcal_EE->Fill(numRcvEcalCells);
       numPrbCellsEcal_EE->Fill(numPrbEcalCells);
    } 
    
    //Ntowers is used in EndJob, occupancy_map is used for occupancy vs ieta
    Ntowers_vs_ieta -> Fill(double(ieta),1.);
    occupancy_map -> Fill(double(ieta),double(iphi)); 

    if((isub == 0 || isub == 1) 
       && (fabs(etaT) <  etaMax[0] && fabs(etaT) >= etaMin[0] )) {
      
      //All cell histos are used
      numBadCellsHcal_HB->Fill(numBadHcalCells);
      numRcvCellsHcal_HB->Fill(numRcvHcalCells);
      numPrbCellsHcal_HB->Fill(numPrbHcalCells);
      
      //Map energy histos are not used
      if (useAllHistos_){
    mapEnergy_HB     -> Fill(double(ieta), double(iphi), en); 
    mapEnergyHcal_HB -> Fill(double(ieta), double(iphi), eH); 
    mapEnergyEcal_HB -> Fill(double(ieta), double(iphi), eE); 
      }
      //      std::cout << " e_ecal = " << eE << std::endl;
      
      if (eH>0.){  // iphi of HCAL component of calotower
      if (ieta>0) meIphiHcalTower_HBP -> Fill(double(iphi)); 
      else        meIphiHcalTower_HBM -> Fill(double(iphi)); 
      }

      //  simple sums
      sumEnergyHcal_HB += eH;
      sumEnergyEcal_HB += eE;
      sumEnergyHO_HB   += eHO;
      
      numFiredTowers_HB++;

      //Not used
      if (useAllHistos_){
    meEnergyEcalTower_HB->Fill(eE);
    meEnergyHcalTower_HB->Fill(eH);    
      }
      
      // MET, SET & phimet
      //  double  etT = cal->et();
      metx_HB += mom.x();   
      mety_HB += mom.y();  //etT * sin(phiT);          
      sEt_HB  += etT;

      //Timing (all histos are used)
      emTiming_HB->Fill(em_tm);
      hadTiming_HB->Fill(had_tm);

      emEnergyTiming_Low_HB->Fill(eE, em_tm);
      emEnergyTiming_HB->Fill(eE, em_tm);
      emEnergyTiming_High_HB->Fill(eE, em_tm);
      emEnergyTiming_profile_Low_HB->Fill(eE, em_tm);
      emEnergyTiming_profile_HB->Fill(eE, em_tm);
      emEnergyTiming_profile_High_HB->Fill(eE, em_tm);

      hadEnergyTiming_Low_HB->Fill(eH, had_tm); 
      hadEnergyTiming_HB->Fill(eH, had_tm); 
      hadEnergyTiming_High_HB->Fill(eH, had_tm); 
      hadEnergyTiming_profile_Low_HB->Fill(eH, had_tm);     
      hadEnergyTiming_profile_HB->Fill(eH, had_tm);     
      hadEnergyTiming_profile_High_HB->Fill(eH, had_tm);     
    }
   
    if((isub == 0 || isub == 2) 
       && (fabs(etaT) <  etaMax[1] && fabs(etaT) >= etaMin[1] )) {
     
      //All cell histos are used
      numBadCellsHcal_HE->Fill(numBadHcalCells);
      numRcvCellsHcal_HE->Fill(numRcvHcalCells);
      numPrbCellsHcal_HE->Fill(numPrbHcalCells);
    
      //Map energy histos are not used
      if (useAllHistos_){
    mapEnergy_HE     -> Fill(double(ieta), double(iphi), en); 
    mapEnergyHcal_HE -> Fill(double(ieta), double(iphi), eH); 
    mapEnergyEcal_HE -> Fill(double(ieta), double(iphi), eE); 
      }
      //      std::cout << " e_ecal = " << eE << std::endl;
      
      if (eH>0.){  // iphi of HCAL component of calotower
      if (ieta>0) meIphiHcalTower_HEP -> Fill(double(iphi)); 
      else        meIphiHcalTower_HEM -> Fill(double(iphi)); 
      }

      //  simple sums
      sumEnergyHcal_HE += eH;
      sumEnergyEcal_HE += eE;
      sumEnergyHO_HE   += eHO;
      
      numFiredTowers_HE++;
      
      //Not used
      if (useAllHistos_){
    meEnergyEcalTower_HE->Fill(eE);
    meEnergyHcalTower_HE->Fill(eH);    
      }
      // MET, SET & phimet
      //  double  etT = cal->et();
      metx_HE += mom.x();   
      mety_HE += mom.y();  //etT * sin(phiT);          
      sEt_HE  += etT;    

      //Timing (all histos are used)
      emTiming_HE->Fill(em_tm);
      hadTiming_HE->Fill(had_tm);

      emEnergyTiming_Low_HE->Fill(eE, em_tm);
      emEnergyTiming_HE->Fill(eE, em_tm);
      emEnergyTiming_profile_Low_HE->Fill(eE, em_tm);
      emEnergyTiming_profile_HE->Fill(eE, em_tm);

      hadEnergyTiming_Low_HE->Fill(eH, had_tm);     
      hadEnergyTiming_HE->Fill(eH, had_tm);   
      hadEnergyTiming_profile_Low_HE->Fill(eH, had_tm);     
      hadEnergyTiming_profile_HE->Fill(eH, had_tm);     
    }

    if((isub == 0 || isub == 3) 
       && (fabs(etaT) <  etaMax[2] && fabs(etaT) >= etaMin[2] )) {
      
      //All cell histos are used
      numBadCellsHcal_HF->Fill(numBadHcalCells);
      numRcvCellsHcal_HF->Fill(numRcvHcalCells);
      numPrbCellsHcal_HF->Fill(numPrbHcalCells);
     
      //Map energy histos are not used
      if (useAllHistos_){
    mapEnergy_HF     -> Fill(double(ieta), double(iphi), en); 
    mapEnergyHcal_HF -> Fill(double(ieta), double(iphi), eH); 
    mapEnergyEcal_HF -> Fill(double(ieta), double(iphi), eE); 
      }
      //      std::cout << " e_ecal = " << eE << std::endl;

      if (ieta>0) meIphiCaloTower_HFP -> Fill(double(iphi)); 
      else        meIphiCaloTower_HFM -> Fill(double(iphi)); 

      //  simple sums
      sumEnergyHcal_HF += eH;
      sumEnergyEcal_HF += eE;
      sumEnergyHO_HF   += eHO;
      
      numFiredTowers_HF++;
     
      //Not used
      if (useAllHistos_){
    meEnergyEcalTower_HF->Fill(eE);
    meEnergyHcalTower_HF->Fill(eH);    
      }
      // MET, SET & phimet
      //  double  etT = cal->et();
      metx_HF += mom.x();   
      mety_HF += mom.y();  //etT * sin(phiT);          
      sEt_HF  += etT;    

      //Timing (all histos are used)
      emTiming_HF->Fill(em_tm);
      hadTiming_HF->Fill(had_tm);
      emEnergyTiming_HF->Fill(eE, em_tm);
      emEnergyTiming_profile_HF->Fill(eE, em_tm);

      hadEnergyTiming_Low_HF->Fill(eH, had_tm);     
      hadEnergyTiming_HF->Fill(eH, had_tm);     
      hadEnergyTiming_profile_Low_HF->Fill(eH, had_tm);     
      hadEnergyTiming_profile_HF->Fill(eH, had_tm);     
    }



  } // end of Towers cycle 

  //These are the six single pion histos; only the second set is used

  mapEnergy_N  -> Fill(double(ieta_MC), double(iphi_MC), 1.); 

  if(isub == 0 || isub == 1) {
    met    = sqrt(metx_HB*metx_HB + mety_HB*mety_HB);
    Vector metv(metx_HB,mety_HB,metz_HB);
    phimet = metv.phi();

    //Five oldest drawn histos first; the rest are not used
    meEnergyHcal_HB->      Fill(sumEnergyHcal_HB);
    meEnergyEcal_HB->      Fill(sumEnergyEcal_HB);
    meNumFiredTowers_HB->  Fill(numFiredTowers_HB);
    MET_HB    -> Fill (met); 
    SET_HB    -> Fill (sEt_HB); 
    
    if (useAllHistos_){
      meEnergyHcalvsEcal_HB->Fill(sumEnergyEcal_HB, sumEnergyHcal_HB);
      meEnergyHO_HB->        Fill(sumEnergyHO_HB);
      meTotEnergy_HB->       Fill(sumEnergyHcal_HB+sumEnergyEcal_HB
                  +sumEnergyHO_HB);    
      phiMET_HB -> Fill (phimet);
    } 
  }    


  if(isub == 0 || isub == 2) {
    met    = sqrt(metx_HE*metx_HE + mety_HE*mety_HE);
    Vector metv(metx_HE,mety_HE,metz_HE);
    phimet = metv.phi();

    //Five oldest drawn histos first; the rest are not used
    meEnergyHcal_HE->      Fill(sumEnergyHcal_HE);
    meEnergyEcal_HE->      Fill(sumEnergyEcal_HE);
    meNumFiredTowers_HE->  Fill(numFiredTowers_HE);
    MET_HE    -> Fill (met); 
    SET_HE    -> Fill (sEt_HE); 
    
    if (useAllHistos_){
      meEnergyHcalvsEcal_HE->Fill(sumEnergyEcal_HE, sumEnergyHcal_HE);
      meEnergyHO_HE->        Fill(sumEnergyHO_HE);
      meTotEnergy_HE->       Fill(sumEnergyHcal_HE+sumEnergyEcal_HE
                  +sumEnergyHO_HE);    
      phiMET_HE -> Fill (phimet); 
    }
  }

  if(isub == 0 || isub == 3) {
    met    = sqrt(metx_HF*metx_HF + mety_HF*mety_HF);
    Vector metv(metx_HF,mety_HF,metz_HF);
    phimet = metv.phi();
    
    //Five oldest drawn histos first; the rest are not used
    meEnergyHcal_HF->      Fill(sumEnergyHcal_HF);
    meEnergyEcal_HF->      Fill(sumEnergyEcal_HF);
    meNumFiredTowers_HF->  Fill(numFiredTowers_HF);
    MET_HF    -> Fill (met); 
    SET_HF    -> Fill (sEt_HF); 

    if (useAllHistos_){
      meEnergyHcalvsEcal_HF->Fill(sumEnergyEcal_HF, sumEnergyHcal_HF);
      meEnergyHO_HF->        Fill(sumEnergyHO_HF);
      meTotEnergy_HF->       Fill(sumEnergyHcal_HF+sumEnergyEcal_HF
                  +sumEnergyHO_HF);    
      phiMET_HF -> Fill (phimet);
    }
  }

}

void CaloTowersAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 453 of file CaloTowersAnalyzer.cc.

References nevent.

{ nevent = 0; }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 50 of file CaloTowersAnalyzer.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::bookProfile(), emean_vs_ieta_E, emean_vs_ieta_E1, emean_vs_ieta_EH, emean_vs_ieta_EH1, emean_vs_ieta_H, emean_vs_ieta_H1, emEnergyTiming_HB, emEnergyTiming_HE, emEnergyTiming_HF, emEnergyTiming_High_HB, emEnergyTiming_Low_HB, emEnergyTiming_Low_HE, emEnergyTiming_profile_HB, emEnergyTiming_profile_HE, emEnergyTiming_profile_HF, emEnergyTiming_profile_High_HB, emEnergyTiming_profile_Low_HB, emEnergyTiming_profile_Low_HE, emTiming_HB, emTiming_HE, emTiming_HF, etaMax, etaMin, hadEnergyTiming_HB, hadEnergyTiming_HE, hadEnergyTiming_HF, hadEnergyTiming_High_HB, hadEnergyTiming_Low_HB, hadEnergyTiming_Low_HE, hadEnergyTiming_Low_HF, hadEnergyTiming_profile_HB, hadEnergyTiming_profile_HE, hadEnergyTiming_profile_HF, hadEnergyTiming_profile_High_HB, hadEnergyTiming_profile_Low_HB, hadEnergyTiming_profile_Low_HE, hadEnergyTiming_profile_Low_HF, hadTiming_HB, hadTiming_HE, hadTiming_HF, hcalselector_, trackerHits::histo, ieta_bins_, ieta_max_, ieta_min_, iphi_bins_, iphi_max_, iphi_min_, isub, mapEnergy_E, mapEnergy_EH, mapEnergy_H, mapEnergy_HB, mapEnergy_HE, mapEnergy_HF, mapEnergy_N, mapEnergyEcal_HB, mapEnergyEcal_HE, mapEnergyEcal_HF, mapEnergyHcal_HB, mapEnergyHcal_HE, mapEnergyHcal_HF, meEnergyEcal_HB, meEnergyEcal_HE, meEnergyEcal_HF, meEnergyEcalTower_HB, meEnergyEcalTower_HE, meEnergyEcalTower_HF, meEnergyHcal_HB, meEnergyHcal_HE, meEnergyHcal_HF, meEnergyHcalTower_HB, meEnergyHcalTower_HE, meEnergyHcalTower_HF, meEnergyHcalvsEcal_HB, meEnergyHcalvsEcal_HE, meEnergyHcalvsEcal_HF, meEnergyHO_HB, meEnergyHO_HE, meEnergyHO_HF, meIphiCaloTower_HFM, meIphiCaloTower_HFP, meIphiHcalTower_HBM, meIphiHcalTower_HBP, meIphiHcalTower_HEM, meIphiHcalTower_HEP, meNumFiredTowers_HB, meNumFiredTowers_HE, meNumFiredTowers_HF, MET_HB, MET_HE, MET_HF, meTotEnergy_HB, meTotEnergy_HE, meTotEnergy_HF, Ntowers_vs_ieta, numBadCellsEcal_EB, numBadCellsEcal_EE, numBadCellsHcal_HB, numBadCellsHcal_HE, numBadCellsHcal_HF, numPrbCellsEcal_EB, numPrbCellsEcal_EE, numPrbCellsHcal_HB, numPrbCellsHcal_HE, numPrbCellsHcal_HF, numRcvCellsEcal_EB, numRcvCellsEcal_EE, numRcvCellsHcal_HB, numRcvCellsHcal_HE, numRcvCellsHcal_HF, occupancy_map, occupancy_vs_ieta, outputFile_, phiMET_HB, phiMET_HE, phiMET_HF, SET_HB, SET_HE, SET_HF, DQMStore::IBooker::setCurrentFolder(), and useAllHistos_.

{  

  etaMin[0] = 0.;
  etaMax[0] = 1.4;
  etaMin[1] = 1.4;
  etaMax[1] = 2.9;
  etaMin[2] = 2.9;
  etaMax[2] = 5.2;

  isub = 0;
  if(hcalselector_ == "HB") isub = 1;
  if(hcalselector_ == "HE") isub = 2;
  if(hcalselector_ == "HF") isub = 3;

  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";
  }
  
  Char_t histo[100];

  
  ibooker.setCurrentFolder("CaloTowersD/CaloTowersTask") ;
  
  

  //These two histos are not drawn by our macros, but they are used
  //in the EndJob for norms and such so I am leaving them alone for now
  //-------------------------------------------------------------------------------------------
  sprintf  (histo, "Ntowers_per_event_vs_ieta" );
  Ntowers_vs_ieta = ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_);
  
  sprintf  (histo, "CaloTowersTask_map_Nentries" );
  mapEnergy_N = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
  //-------------------------------------------------------------------------------------------

  //These the single pion scan histos
  //-------------------------------------------------------------------------------------------
  //The first three are not used
  if (useAllHistos_){
    sprintf  (histo, "emean_vs_ieta_E" );
    emean_vs_ieta_E = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
    sprintf  (histo, "emean_vs_ieta_H" );
    emean_vs_ieta_H = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
    sprintf  (histo, "emean_vs_ieta_EH" );
    emean_vs_ieta_EH = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
  }
  //These are drawn
  sprintf  (histo, "emean_vs_ieta_E1" );
  emean_vs_ieta_E1 = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
  sprintf  (histo, "emean_vs_ieta_H1" );
  emean_vs_ieta_H1 = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
  sprintf  (histo, "emean_vs_ieta_EH1" );
  emean_vs_ieta_EH1 = ibooker.bookProfile(histo, histo, ieta_bins_, ieta_min_, ieta_max_, 2100, -100., 2000.);
  //-------------------------------------------------------------------------------------------

  //Map energy histos are not drawn
  if (useAllHistos_){
    sprintf  (histo, "CaloTowersTask_map_energy_E" );
    mapEnergy_E = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
    sprintf  (histo, "CaloTowersTask_map_energy_H");
    mapEnergy_H = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
    sprintf  (histo, "CaloTowersTask_map_energy_EH" );
    mapEnergy_EH = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
  }

  //All ECAL cell histos are used
  // XXX: ECAL 0-25 [0-26, 26 bins]   HCAL 0-4 [0-5, 5 bins]
  sprintf  (histo, "number_of_bad_cells_Ecal_EB");
  numBadCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf  (histo, "number_of_bad_cells_Ecal_EE");
  numBadCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf  (histo, "number_of_recovered_cells_Ecal_EB");
  numRcvCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf  (histo, "number_of_recovered_cells_Ecal_EE");
  numRcvCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf  (histo, "number_of_problematic_cells_Ecal_EB");
  numPrbCellsEcal_EB = ibooker.book1D(histo, histo, 26, 0, 26);
  sprintf  (histo, "number_of_problematic_cells_Ecal_EE");
  numPrbCellsEcal_EE = ibooker.book1D(histo, histo, 26, 0, 26); 

  //Occupancy vs. ieta is drawn, occupancy map is needed to draw it 
  sprintf  (histo, "CaloTowersTask_map_occupancy" );
  occupancy_map = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);

  sprintf  (histo, "CaloTowersTask_occupancy_vs_ieta" );
  occupancy_vs_ieta = ibooker.book1D(histo, histo, ieta_bins_, ieta_min_, ieta_max_);
  
  if( isub == 1 || isub == 0) {
    //All cell histos are used
    sprintf  (histo, "number_of_bad_cells_Hcal_HB");
    numBadCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf  (histo, "number_of_recovered_cells_Hcal_HB");
    numRcvCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf  (histo, "number_of_problematic_cells_Hcal_HB");
    numPrbCellsHcal_HB = ibooker.book1D(histo, histo, 5, 0, 5); 
    
    //These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL, MET and SET
    //-------------------------------------------------------------------------------------------
    sprintf (histo, "CaloTowersTask_energy_HCAL_HB" ) ;
    meEnergyHcal_HB    = ibooker.book1D(histo, histo, 4100, -200, 8000);  
    
    sprintf (histo, "CaloTowersTask_energy_ECAL_HB" ) ;
    meEnergyEcal_HB    = ibooker.book1D(histo, histo, 3100, -200, 6000); 
    
    sprintf (histo, "CaloTowersTask_number_of_fired_towers_HB" ) ;
    meNumFiredTowers_HB = ibooker.book1D(histo, histo, 1000, 0, 2000); 

    sprintf  (histo, "CaloTowersTask_MET_HB" ) ;
    MET_HB = ibooker.book1D(histo, histo, 3000, 0. , 3000. ) ;
    
    sprintf  (histo, "CaloTowersTask_SET_HB" ) ;
    SET_HB = ibooker.book1D(histo, histo, 8000, 0. , 8000. ) ;
    //-------------------------------------------------------------------------------------------

    //Timing histos and profiles -- all six are necessary
    //-------------------------------------------------------------------------------------------
    sprintf  (histo, "CaloTowersTask_EM_Timing_HB" ) ;
    emTiming_HB = ibooker.book1D(histo, histo, 110, -120., 100. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Timing_HB" ) ;
    hadTiming_HB = ibooker.book1D(histo, histo, 70, -48., 92. ) ;

    //Energy-Timing histos are divided into low, medium and high to reduce memory usage
    //EM
    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_Low_HB" ) ;
    emEnergyTiming_Low_HB = ibooker.book2D(histo, histo, 40, 0. , 40., 110, -120., 100.  ) ;

    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_HB" ) ;
    emEnergyTiming_HB = ibooker.book2D(histo, histo, 200, 0. , 400., 110, -120., 100.  ) ;

    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_High_HB" ) ;
    emEnergyTiming_High_HB = ibooker.book2D(histo, histo, 200, 0. , 3000., 110, -120., 100.  ) ;

    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_Low_HB" ) ;
    emEnergyTiming_profile_Low_HB = ibooker.bookProfile(histo, histo, 40, 0. , 40., 110, -120., 100.  ) ;

    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_HB" ) ;
    emEnergyTiming_profile_HB = ibooker.bookProfile(histo, histo, 200, 0. , 400., 110, -120., 100.  ) ;

    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_High_HB" ) ;
    emEnergyTiming_profile_High_HB = ibooker.bookProfile(histo, histo, 200, 0. , 3000., 110, -120., 100.  ) ;

    //HAD
    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_Low_HB" ) ;
    hadEnergyTiming_Low_HB = ibooker.book2D(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_HB" ) ;
    hadEnergyTiming_HB = ibooker.book2D(histo, histo, 100, 0. , 200., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_High_HB" ) ;
    hadEnergyTiming_High_HB = ibooker.book2D(histo, histo, 300, 0. , 3000., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HB" ) ;
    hadEnergyTiming_profile_Low_HB = ibooker.bookProfile(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_HB" ) ;
    hadEnergyTiming_profile_HB = ibooker.bookProfile(histo, histo, 100, 0. , 200., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_High_HB" ) ;
    hadEnergyTiming_profile_High_HB = ibooker.bookProfile(histo, histo, 300, 0. , 3000., 70, -48., 92. ) ;
    //-------------------------------------------------------------------------------------------

    sprintf (histo, "CaloTowersTask_Iphi_HCAL_component_of_tower_HBP" ) ;
    meIphiHcalTower_HBP = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_); 
    sprintf (histo, "CaloTowersTask_Iphi_HCAL_component_of_tower_HBM" ) ;
    meIphiHcalTower_HBM = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_); 
        
    //-------------------------------------------------------------------------------------------

    //Everything else is not drawn
    if (useAllHistos_){
      sprintf (histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HB") ;
      meEnergyHcalvsEcal_HB    = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);
      
      sprintf (histo, "CaloTowersTask_energy_OUTER_HB" ) ;
      meEnergyHO_HB    = ibooker.book1D(histo, histo, 1640, -200, 8000);   
      
      sprintf (histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HB" ) ;
      meEnergyEcalTower_HB = ibooker.book1D(histo, histo, 440, -200, 2000); 
      
      sprintf (histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HB" ) ;
      meEnergyHcalTower_HB = ibooker.book1D(histo, histo, 440 , -200 , 2000); 
      
      sprintf  (histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HB" ) ;
      meTotEnergy_HB = ibooker.book1D(histo, histo,400, 0., 2000.) ;
      
      sprintf  (histo, "CaloTowersTask_map_energy_HB" );
      mapEnergy_HB = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf  (histo, "CaloTowersTask_map_energy_HCAL_HB");
      mapEnergyHcal_HB = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf  (histo, "CaloTowersTask_map_energy_ECAL_HB" );
      mapEnergyEcal_HB = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      
      sprintf  (histo, "CaloTowersTask_phi_MET_HB" ) ;
      phiMET_HB = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898 ) ;
    }
  } 
  
  if( isub == 2 || isub == 0) {
    //All cell histos are used
    sprintf  (histo, "number_of_bad_cells_Hcal_HE");
    numBadCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf  (histo, "number_of_recovered_cells_Hcal_HE");
    numRcvCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf  (histo, "number_of_problematic_cells_Hcal_HE");
    numPrbCellsHcal_HE = ibooker.book1D(histo, histo, 5, 0, 5); 

    //These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL, MET and SET
    //-------------------------------------------------------------------------------------------
    sprintf (histo, "CaloTowersTask_energy_HCAL_HE" ) ;
    meEnergyHcal_HE    = ibooker.book1D(histo, histo, 1240, -200, 6000);
    
    sprintf (histo, "CaloTowersTask_energy_ECAL_HE" ) ;
    meEnergyEcal_HE    = ibooker.book1D(histo, histo, 1240, -200, 6000);
    
    sprintf (histo, "CaloTowersTask_number_of_fired_towers_HE" ) ;
    meNumFiredTowers_HE = ibooker.book1D(histo, histo, 1000, 0, 2000); 

    sprintf  (histo, "CaloTowersTask_MET_HE" ) ;
    MET_HE = ibooker.book1D(histo, histo, 1000, 0. , 1000. ) ;
    
    sprintf  (histo, "CaloTowersTask_SET_HE" ) ;
    SET_HE = ibooker.book1D(histo, histo, 2000, 0. , 2000. ) ;
    //-------------------------------------------------------------------------------------------

    //Timing histos and profiles -- all six are necessary
    //-------------------------------------------------------------------------------------------
    sprintf  (histo, "CaloTowersTask_EM_Timing_HE" ) ;
    emTiming_HE = ibooker.book1D(histo, histo, 110, -120., 100. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Timing_HE" ) ;
    hadTiming_HE = ibooker.book1D(histo, histo, 70, -48., 92. ) ;

    //Energy-Timing histos are divided into low and normal to reduce memory usage
    //EM
    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_Low_HE" ) ;
    emEnergyTiming_Low_HE = ibooker.book2D(histo, histo, 160, 0. , 160., 110, -120., 100.  ) ;

    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_HE" ) ;
    emEnergyTiming_HE = ibooker.book2D(histo, histo, 200, 0. , 800., 110, -120., 100.  ) ;

    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_Low_HE" ) ;
    emEnergyTiming_profile_Low_HE = ibooker.bookProfile(histo, histo, 160, 0. , 160., 110, -120., 100.  ) ;

    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_HE" ) ;
    emEnergyTiming_profile_HE = ibooker.bookProfile(histo, histo, 200, 0. , 800., 110, -120., 100.  ) ;

    //HAD
    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_Low_HE" ) ;
    hadEnergyTiming_Low_HE = ibooker.book2D(histo, histo, 160, 0. , 160., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_HE" ) ;
    hadEnergyTiming_HE = ibooker.book2D(histo, histo, 200, 0. , 800., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HE" ) ;
    hadEnergyTiming_profile_Low_HE = ibooker.bookProfile(histo, histo, 160, 0. , 160., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_HE" ) ;
    hadEnergyTiming_profile_HE = ibooker.bookProfile(histo, histo, 200, 0. , 800., 70, -48., 92. ) ;
    //-------------------------------------------------------------------------------------------

    sprintf (histo, "CaloTowersTask_Iphi_HCAL_component_of_tower_HEP" ) ;
    meIphiHcalTower_HEP = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_); 
    sprintf (histo, "CaloTowersTask_Iphi_HCAL_component_of_tower_HEM" ) ;
    meIphiHcalTower_HEM = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_); 

    //-------------------------------------------------------------------------------------------

    //Everything else is not drawn
    if (useAllHistos_){
      sprintf (histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HE") ;
      meEnergyHcalvsEcal_HE    = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);
      
      sprintf (histo, "CaloTowersTask_energy_OUTER_HE" ) ;
      meEnergyHO_HE    = ibooker.book1D(histo, histo, 440, -200, 2000);   
      
      sprintf (histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HE" ) ;
      meEnergyEcalTower_HE = ibooker.book1D(histo, histo, 1100, -200, 2000); 
      
      sprintf (histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HE" ) ;
      meEnergyHcalTower_HE = ibooker.book1D(histo, histo, 1100 , -200 , 2000); 
      
      sprintf  (histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HE" ) ;
      meTotEnergy_HE = ibooker.book1D(histo, histo,400, 0., 2000.) ;
      
      sprintf  (histo, "CaloTowersTask_map_energy_HE" );
      mapEnergy_HE = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf  (histo, "CaloTowersTask_map_energy_HCAL_HE");
      mapEnergyHcal_HE = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf  (histo, "CaloTowersTask_map_energy_ECAL_HE" );
      mapEnergyEcal_HE = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      
      sprintf  (histo, "CaloTowersTask_phi_MET_HE" ) ;
      phiMET_HE = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898 ) ;
    }
  }


  if( isub == 3 || isub == 0) {
    //All cell histos are used
    sprintf  (histo, "number_of_bad_cells_Hcal_HF");
    numBadCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf  (histo, "number_of_recovered_cells_Hcal_HF");
    numRcvCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5);
    sprintf  (histo, "number_of_problematic_cells_Hcal_HF");
    numPrbCellsHcal_HF = ibooker.book1D(histo, histo, 5, 0, 5); 

    //These are the five oldest CaloTower histos used: NTowers, E in HCAL/ECAL, MET and SET
    //-------------------------------------------------------------------------------------------
    sprintf (histo, "CaloTowersTask_energy_HCAL_HF" ) ;
    meEnergyHcal_HF    = ibooker.book1D(histo, histo, 4040, -200, 20000);
    
    sprintf (histo, "CaloTowersTask_energy_ECAL_HF" ) ;
    meEnergyEcal_HF    = ibooker.book1D(histo, histo, 2440, -200, 12000);

    sprintf (histo, "CaloTowersTask_number_of_fired_towers_HF" ) ;
    meNumFiredTowers_HF = ibooker.book1D(histo, histo, 1000, 0, 2000); 

    sprintf  (histo, "CaloTowersTask_MET_HF" ) ;
    MET_HF = ibooker.book1D(histo, histo, 500, 0. , 500. ) ;
    
    sprintf  (histo, "CaloTowersTask_SET_HF" ) ;
    SET_HF = ibooker.book1D(histo, histo, 2000, 0. , 2000. ) ;
    //-------------------------------------------------------------------------------------------

    //Timing histos and profiles -- all six are necessary
    //-------------------------------------------------------------------------------------------
    sprintf  (histo, "CaloTowersTask_EM_Timing_HF" ) ;
    emTiming_HF = ibooker.book1D(histo, histo, 110, -120., 100. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Timing_HF" ) ;
    hadTiming_HF = ibooker.book1D(histo, histo, 70, -48., 92. ) ;

    //EM
    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_HF" ) ;
    emEnergyTiming_HF = ibooker.book2D(histo, histo, 150, 0. , 300., 110, -120., 100.  ) ;

    sprintf  (histo, "CaloTowersTask_EM_Energy_Timing_profile_HF" ) ;
    emEnergyTiming_profile_HF = ibooker.bookProfile(histo, histo, 150, 0. , 300., 110, -120., 100.  ) ;

    //HAD (requires two different sets of histograms to lower RAM usage)
    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_Low_HF" ) ;
    hadEnergyTiming_Low_HF = ibooker.book2D(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_HF" ) ;
    hadEnergyTiming_HF = ibooker.book2D(histo, histo, 200, 0. , 600., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_Low_HF" ) ;
    hadEnergyTiming_profile_Low_HF = ibooker.bookProfile(histo, histo, 40, 0. , 40., 70, -48., 92. ) ;

    sprintf  (histo, "CaloTowersTask_HAD_Energy_Timing_profile_HF" ) ;
    hadEnergyTiming_profile_HF = ibooker.bookProfile(histo, histo, 200, 0. , 600., 70, -48., 92. ) ;
    //-------------------------------------------------------------------------------------------

    sprintf (histo, "CaloTowersTask_Iphi_tower_HFP" ) ;
    meIphiCaloTower_HFP = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_); 
    sprintf (histo, "CaloTowersTask_Iphi_tower_HFM" ) ;
    meIphiCaloTower_HFM = ibooker.book1D(histo, histo, iphi_bins_, iphi_min_, iphi_max_); 

    //-------------------------------------------------------------------------------------------

    //Everything else is not drawn
    if (useAllHistos_){
      sprintf (histo, "CaloTowersTask_sum_of_energy_HCAL_vs_ECAL_HF") ;
      meEnergyHcalvsEcal_HF    = ibooker.book2D(histo, histo, 500, 0., 500., 500, 0., 500.);
      
      sprintf (histo, "CaloTowersTask_energy_OUTER_HF" ) ;
      meEnergyHO_HF    = ibooker.book1D(histo, histo, 440, -200, 2000);   
      
      sprintf (histo, "CaloTowersTask_energy_of_ECAL_component_of_tower_HF" ) ;
      meEnergyEcalTower_HF = ibooker.book1D(histo, histo, 440, -200, 2000); 
      
      sprintf (histo, "CaloTowersTask_energy_of_HCAL_component_of_tower_HF" ) ;
      meEnergyHcalTower_HF = ibooker.book1D(histo, histo, 440 , -200 , 2000); 
      
      sprintf  (histo, "CaloTowersTask_energy_HcalPlusEcalPlusHO_HF" ) ;
      meTotEnergy_HF = ibooker.book1D(histo, histo, 400, 0., 2000.) ;
      
      sprintf  (histo, "CaloTowersTask_map_energy_HF" );
      mapEnergy_HF = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf  (histo, "CaloTowersTask_map_energy_HCAL_HF");
      mapEnergyHcal_HF = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      sprintf  (histo, "CaloTowersTask_map_energy_ECAL_HF" );
      mapEnergyEcal_HF = ibooker.book2D(histo, histo, ieta_bins_, ieta_min_, ieta_max_, iphi_bins_, iphi_min_, iphi_max_);
      
      sprintf  (histo, "CaloTowersTask_phi_MET_HF" ) ;
      phiMET_HF = ibooker.book1D(histo, histo, 72, -3.1415926535898, 3.1415926535898 ) ;
    }
  }

}

void CaloTowersAnalyzer::dqmBeginRun ( const edm::Run &  run, const edm::EventSetup &  c  )

override

Definition at line 18 of file CaloTowersAnalyzer.cc.

References edm::EventSetup::get(), HcalDDDRecConstants::getEtaRange(), HcalDDDRecConstants::getNPhi(), hcons, ieta_bins_, ieta_max_, ieta_min_, hcalTTPDigis_cfi::iEtaMax, createfilelist::int, iphi_bins_, iphi_max_, and iphi_min_.

                                                                              {
  
    edm::ESHandle<HcalDDDRecConstants> pHRNDC;
    es.get<HcalRecNumberingRecord>().get( pHRNDC );
    hcons = &(*pHRNDC);

    //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_);

    //Retain classic behavior, all plots have same ieta range.

    int iEtaMax = (hcons->getEtaRange(0).second > hcons->getEtaRange(1).second ? hcons->getEtaRange(0).second : hcons->getEtaRange(1).second);
    iEtaMax = (iEtaMax > hcons->getEtaRange(2).second ? iEtaMax : hcons->getEtaRange(2).second);
    iEtaMax = (iEtaMax > hcons->getEtaRange(3).second ? iEtaMax : hcons->getEtaRange(3).second);

    //Give an empty bin around the subdet ieta range to make it clear that all ieta rings have been included
    ieta_min_ = -iEtaMax - 1.5;
    ieta_max_ = iEtaMax + 1.5;
    ieta_bins_ = (int) (ieta_max_ - ieta_min_);

}

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

private

Definition at line 800 of file CaloTowersAnalyzer.cc.

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

                                                                                { 
  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 CaloTowersAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDProducerBase.

Definition at line 450 of file CaloTowersAnalyzer.cc.

                                {
}

Member Data Documentation

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_E

private

Definition at line 83 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_E1

private

Definition at line 87 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_EH

private

Definition at line 85 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_EH1

private

Definition at line 89 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_H

private

Definition at line 84 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emean_vs_ieta_H1

private

Definition at line 88 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HB

private

Definition at line 144 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HE

private

Definition at line 182 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_HF

private

Definition at line 214 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_High_HB

private

Definition at line 145 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_Low_HB

private

Definition at line 143 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_Low_HE

private

Definition at line 181 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HB

private

Definition at line 147 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HE

private

Definition at line 184 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_HF

private

Definition at line 215 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_High_HB

private

Definition at line 148 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_Low_HB

private

Definition at line 146 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emEnergyTiming_profile_Low_HE

private

Definition at line 183 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emTiming_HB

private

Definition at line 140 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emTiming_HE

private

Definition at line 178 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::emTiming_HF

private

Definition at line 212 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

double CaloTowersAnalyzer::etaMax[3]

private

Definition at line 70 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

double CaloTowersAnalyzer::etaMin[3]

private

Definition at line 71 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HB

private

Definition at line 151 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HE

private

Definition at line 187 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_HF

private

Definition at line 218 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_High_HB

private

Definition at line 152 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HB

private

Definition at line 150 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HE

private

Definition at line 186 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_Low_HF

private

Definition at line 217 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HB

private

Definition at line 154 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HE

private

Definition at line 189 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_HF

private

Definition at line 220 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_High_HB

private

Definition at line 155 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HB

private

Definition at line 153 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HE

private

Definition at line 188 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadEnergyTiming_profile_Low_HF

private

Definition at line 219 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadTiming_HB

private

Definition at line 141 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadTiming_HE

private

Definition at line 179 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::hadTiming_HF

private

Definition at line 213 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::string CaloTowersAnalyzer::hcalselector_

private

Definition at line 56 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and CaloTowersAnalyzer().

const HcalDDDRecConstants* CaloTowersAnalyzer::hcons

private

Definition at line 74 of file CaloTowersAnalyzer.h.

Referenced by dqmBeginRun().

int CaloTowersAnalyzer::ieta_bins_

private

Definition at line 79 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float CaloTowersAnalyzer::ieta_max_

private

Definition at line 80 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float CaloTowersAnalyzer::ieta_min_

private

Definition at line 80 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int CaloTowersAnalyzer::imc

private

Definition at line 67 of file CaloTowersAnalyzer.h.

int CaloTowersAnalyzer::iphi_bins_

private

Definition at line 76 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float CaloTowersAnalyzer::iphi_max_

private

Definition at line 77 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

float CaloTowersAnalyzer::iphi_min_

private

Definition at line 77 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and dqmBeginRun().

int CaloTowersAnalyzer::isub

private

Definition at line 64 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_E

private

Definition at line 96 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_EH

private

Definition at line 98 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_H

private

Definition at line 97 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HB

private

Definition at line 133 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HE

private

Definition at line 171 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_HF

private

Definition at line 205 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergy_N

private

Definition at line 99 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HB

private

Definition at line 134 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HE

private

Definition at line 172 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyEcal_HF

private

Definition at line 206 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HB

private

Definition at line 135 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HE

private

Definition at line 173 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::mapEnergyHcal_HF

private

Definition at line 207 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

std::string CaloTowersAnalyzer::mc_

private

Definition at line 57 of file CaloTowersAnalyzer.h.

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HB

private

Definition at line 122 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HE

private

Definition at line 160 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcal_HF

private

Definition at line 194 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HB

private

Definition at line 126 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HE

private

Definition at line 164 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyEcalTower_HF

private

Definition at line 198 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HB

private

Definition at line 123 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HE

private

Definition at line 161 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcal_HF

private

Definition at line 195 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HB

private

Definition at line 127 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HE

private

Definition at line 165 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalTower_HF

private

Definition at line 199 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HB

private

Definition at line 120 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HE

private

Definition at line 158 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHcalvsEcal_HF

private

Definition at line 192 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HB

private

Definition at line 121 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HE

private

Definition at line 159 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meEnergyHO_HF

private

Definition at line 193 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiCaloTower_HFM

private

Definition at line 203 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiCaloTower_HFP

private

Definition at line 202 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiHcalTower_HBM

private

Definition at line 131 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiHcalTower_HBP

private

Definition at line 130 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiHcalTower_HEM

private

Definition at line 169 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meIphiHcalTower_HEP

private

Definition at line 168 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HB

private

Definition at line 124 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HE

private

Definition at line 162 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meNumFiredTowers_HF

private

Definition at line 196 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::MET_HB

private

Definition at line 136 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::MET_HE

private

Definition at line 174 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::MET_HF

private

Definition at line 208 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HB

private

Definition at line 128 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HE

private

Definition at line 166 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::meTotEnergy_HF

private

Definition at line 200 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

int CaloTowersAnalyzer::nevent

private

Definition at line 65 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and beginJob().

MonitorElement* CaloTowersAnalyzer::Ntowers_vs_ieta

private

Definition at line 91 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsEcal_EB

private

Definition at line 102 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsEcal_EE

private

Definition at line 103 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HB

private

Definition at line 109 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HE

private

Definition at line 110 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numBadCellsHcal_HF

private

Definition at line 111 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsEcal_EB

private

Definition at line 106 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsEcal_EE

private

Definition at line 107 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HB

private

Definition at line 115 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HE

private

Definition at line 116 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numPrbCellsHcal_HF

private

Definition at line 117 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsEcal_EB

private

Definition at line 104 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsEcal_EE

private

Definition at line 105 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HB

private

Definition at line 112 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HE

private

Definition at line 113 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::numRcvCellsHcal_HF

private

Definition at line 114 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::occupancy_map

private

Definition at line 92 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::occupancy_vs_ieta

private

Definition at line 93 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms().

std::string CaloTowersAnalyzer::outputFile_

private

Definition at line 55 of file CaloTowersAnalyzer.h.

Referenced by bookHistograms(), and CaloTowersAnalyzer().

MonitorElement* CaloTowersAnalyzer::phiMET_HB

private

Definition at line 138 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::phiMET_HE

private

Definition at line 176 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::phiMET_HF

private

Definition at line 210 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::SET_HB

private

Definition at line 137 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::SET_HE

private

Definition at line 175 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* CaloTowersAnalyzer::SET_HF

private

Definition at line 209 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<CaloTowerCollection> CaloTowersAnalyzer::tok_towers_

private

Definition at line 62 of file CaloTowersAnalyzer.h.

Referenced by analyze(), and CaloTowersAnalyzer().

bool CaloTowersAnalyzer::useAllHistos_

private

Definition at line 58 of file CaloTowersAnalyzer.h.

Referenced by analyze(), bookHistograms(), and CaloTowersAnalyzer().