#include <HcalCaloTowerMonitor.h>

Inheritance diagram for HcalCaloTowerMonitor:

Public Member Functions
int	getIeta (double eta)

	HcalCaloTowerMonitor ()

void	processEvent (const CaloTowerCollection &calotower)

void	reset ()

void	setup (const edm::ParameterSet &ps, DQMStore *dbe)

	~HcalCaloTowerMonitor ()

Public Member Functions inherited from HcalBaseMonitor
void	beginLuminosityBlock (int lb)

virtual void	beginRun ()

virtual void	clearME ()

virtual void	done ()

void	endLuminosityBlock ()

bool	getDiagnostics () const

int	getVerbosity () const

	HcalBaseMonitor ()

void	hideKnownBadCells ()

virtual void	periodicReset ()

void	processEvent ()

void	setDiagnostics (bool myval)

void	setMinMaxHists1D (std::vector< MonitorElement * > &hh, double min, double max)

void	setMinMaxHists2D (std::vector< MonitorElement * > &hh, double min, double max)

void	setupDepthHists1D (MonitorElement &h, std::vector< MonitorElement > &hh, std::string Name, std::string Units, int lowbound, int highbound, int Nbins)

void	setupDepthHists1D (std::vector< MonitorElement * > &hh, std::string Name, std::string Units, int lowbound, int highbound, int Nbins)

void	setupDepthHists2D (MonitorElement &h, std::vector< MonitorElement > &hh, std::string Name, std::string Units)

void	setupDepthHists2D (std::vector< MonitorElement * > &hh, std::string Name, std::string Units)

void	setupDepthHists2D (MonitorElement &h, std::vector< MonitorElement > &hh, std::string Name, std::string Units, int nbinsx, int lowboundx, int highboundx, int nbinsy, int lowboundy, int highboundy)

void	setupDepthHists2D (std::vector< MonitorElement * > &hh, std::string Name, std::string Units, int nbinsx, int lowboundx, int highboundx, int nbinsy, int lowboundy, int highboundy)

void	SetupEtaPhiHists (MonitorElement *&h, EtaPhiHists &hh, std::string Name, std::string Units)

void	SetupEtaPhiHists (EtaPhiHists &hh, std::string Name, std::string Units)

void	setVerbosity (int verb)

bool	vetoCell (HcalDetId &id)

virtual	~HcalBaseMonitor ()

Private Attributes
MonitorElement *	caloTowerEnergy

MonitorElement *	caloTowerEnergyMap

MonitorElement *	caloTowerMeanEnergyEta

MonitorElement *	caloTowerOccMap

MonitorElement *	caloTowerTime

bool	debug_
	Monitoring elements. More...

MonitorElement *	ecalEnergy

MonitorElement *	ecalEnergyMap

MonitorElement *	ecalMeanEnergyEta

MonitorElement *	ecalOccMap

MonitorElement *	ecalTime

MonitorElement *	energy_CaloTowervsEcal

MonitorElement *	energy_CaloTowervsHcal

MonitorElement *	energy_HcalvsEcal

int	etaBins_

double	etaMax_

double	etaMin_

MonitorElement *	hcalEnergy

MonitorElement *	hcalEnergyMap

MonitorElement *	hcalMeanEnergyEta

MonitorElement *	hcalOccMap

MonitorElement *	hcalTime

int	ievt_

MonitorElement *	meEVT_

int	phiBins_

double	phiMax_

double	phiMin_

MonitorElement *	time_CaloTowervsEcal

MonitorElement *	time_CaloTowervsHcal

MonitorElement *	time_HcalvsEcal

Additional Inherited Members
Protected Member Functions inherited from HcalBaseMonitor
void	LumiBlockUpdate (int lb)

Protected Attributes inherited from HcalBaseMonitor
std::vector< int >	AllowedCalibTypes_

std::vector< std::string >	badCells_

std::string	baseFolder_

bool	checkHB_

bool	checkHE_

bool	checkHF_

bool	checkHO_

int	checkNevents_

edm::CPUTimer	cpu_timer

bool	dump2database

int	etaBins_

double	etaMax_

double	etaMin_

int	fVerbosity

int	ievt_

bool	LBprocessed_

int	levt_

int	lumiblock

DQMStore *	m_dbe

bool	makeDiagnostics

MonitorElement *	meEVT_

MonitorElement *	meTOTALEVT_

double	minErrorFlag_

int	Nlumiblocks_

int	NumBadHB

int	NumBadHE

int	NumBadHF

int	NumBadHO

int	oldlumiblock

bool	Online_

int	phiBins_

double	phiMax_

double	phiMin_

MonitorElement *	ProblemCells

EtaPhiHists	ProblemCellsByDepth

MonitorElement *	ProblemsVsLB

MonitorElement *	ProblemsVsLB_HB

MonitorElement *	ProblemsVsLB_HBHEHF

MonitorElement *	ProblemsVsLB_HE

MonitorElement *	ProblemsVsLB_HF

MonitorElement *	ProblemsVsLB_HO

int	resetNevents_

std::string	rootFolder_

bool	showTiming

int	tevt_

Detailed Description

Date:: 2008/08/14 18:40:28

Revision:: 1.2

Author: J. Temple - Univ. of Maryland

Definition at line 19 of file HcalCaloTowerMonitor.h.

Constructor & Destructor Documentation

HcalCaloTowerMonitor::HcalCaloTowerMonitor ( )

Definition at line 4 of file HcalCaloTowerMonitor.cc.

4 {}

HcalCaloTowerMonitor::~HcalCaloTowerMonitor ( )

Definition at line 6 of file HcalCaloTowerMonitor.cc.

6 {}

Member Function Documentation

int HcalCaloTowerMonitor::getIeta ( double eta )

Definition at line 200 of file HcalCaloTowerMonitor.cc.

References neg.

Referenced by processEvent().

 {
   // return Ieta index given physical eta value
 
   double myeta=fabs(eta);
   int neg=int(eta/myeta); // eta can be negative or positive
 
   if (myeta<=1.740) // first 20 bins are spaced evenly in increments of .087 radians, starting at ieta=1
     return (1+int(myeta/.087))*neg;
 
   if (fabs(eta)<=1.83)
     return int(21*neg);
   if (fabs(eta)<=1.93)
     return int(22*neg);
   if (fabs(eta)<=2.043)
     return int(23*neg);
   if (fabs(eta)<=2.172)
     return int(24*neg);
   if (fabs(eta)<=2.322)
     return int(25*neg);
          
   if (fabs(eta)<=2.5)
     return int(26*neg);
   if (fabs(eta)<=2.65)
     return int(27*neg);
 
   // 28 and 29 aren't quite right -- eta values depend on depth
   if (fabs(eta)<=2.853)
     return int(28*neg);
   if (fabs(eta)<=2.964)
     return int(29*neg);
   if (fabs(eta)<=3.139)
     return int(30*neg);
          
   if (fabs(eta)<=3.314)
     return int(31*neg);
   if (fabs(eta)<=3.489)
     return int(32*neg);
   if (fabs(eta)<=3.664)
     return int(33*neg);
   if (fabs(eta)<=3.839)
     return int(34*neg);
   if (fabs(eta)<=4.013)
     return int(35*neg);
      
   if (fabs(eta)<=4.191)
     return int(36*neg);
   if (fabs(eta)<=4.363)
     return int(37*neg);
   if (fabs(eta)<=4.538)
     return int(38*neg);
   if (fabs(eta)<=4.716)
     return int(39*neg);
   if (fabs(eta)<=4.889)
     return int(40*neg);
   if (fabs(eta)<=5.191)
     return int(41*neg);
 
   // Anything with large eta is outside calorimeter; skip it?
   // Are there ZDC-based calotowers?
   return 42;
     
 } // int HcalCaloTowerMonitor::getIeta(double eta)

void HcalCaloTowerMonitor::processEvent ( const CaloTowerCollection & calotower )

Definition at line 110 of file HcalCaloTowerMonitor.cc.

References edm::SortedCollection< T, SORT >::begin(), caloTowerEnergy, caloTowerEnergyMap, caloTowerMeanEnergyEta, caloTowerOccMap, caloTowerTime, gather_cfg::cout, ecalEnergy, ecalEnergyMap, ecalMeanEnergyEta, ecalOccMap, ecalTime, edm::SortedCollection< T, SORT >::end(), energy_HcalvsEcal, eta(), MonitorElement::Fill(), HcalBaseMonitor::fVerbosity, getIeta(), hcalEnergy, hcalEnergyMap, hcalMeanEnergyEta, hcalOccMap, hcalTime, i, HcalBaseMonitor::m_dbe, phi, and time_HcalvsEcal.

 {
   // fill histograms for each event
   if(!m_dbe)  
     { 
       if(fVerbosity) std::cout <<"<HcalCaloTowerMonitor::processEvent>    DQMStore not instantiated!!!\n"; 
       return; 
     }  // if(!m_dbe)
   
   if (fVerbosity) std::cout <<"Processing calotower"<<std::endl;
 
 
   // Store values of hcal, energy for forming averages at each eta
   float hcalenergy[90]={0.};
   float ecalenergy[90]={0.};
   int etacounts[90]={0};
 
   for (CaloTowerCollection::const_iterator it=calotower.begin();
        it!=calotower.end();++it)
     {
       //caloTowers give eta and phi; need to map to ieta,iphi
       // phi is mapped into 72 equally spaced segments (2pi/72=.087 radians each), starting at iphi=1
       // at ieta=21 (eta=1.74), iphi spaced in 10 degree increments 
       // at ieta=40 (eta=4.716), iphi spaced in 20 degree increments
       // However, segmentation values are similar
       // (iphi=1,2,3,4,5 for 5 degree segments,
       //  iphi 1,3,5,... for 10 degree segments, and iphi=3,7,... for 20 deg.)
       int iphi;
       double eta=it->eta();
       double phi=it->phi();
       // if phi<0; shift it by 2pi (so that phi runs from 0->2pi, rather than -pi->pi)
       if (phi<0)
     phi+=2*3.14159265359;
       iphi = int(phi/.087)+1;
       if (fabs(eta)>1.74)
     {
       if (fabs(eta)<=4.716)
         iphi=(iphi-1)/2*2+1;
       else
         iphi=(iphi-1)/4*4+3;
       
     }
       //ieta is more complicated -- it runs in segments of .087 for the first 20 segments, then the segmentation becomes non-uniform
       int ieta=getIeta(eta);
 
       // Add 42 to the indices.  Only index[1...] will have values within [HB...HF].  Index[0] and [84] store values outside the expected calorimeter range.
       hcalenergy[ieta+42]+=it->hadEnergy();
       ecalenergy[ieta+42]+=it->emEnergy();
       etacounts[ieta+42]++;
 
       /*
       std::cout <<"CALOTOWER eta = \t"<<eta<<"\tphi = "<<it->phi()<<std::endl;
       std::cout <<"\t ieta = \t"<<ieta<<"\tiphi = "<<iphi<<std::endl;
       */
 
       // Fill histograms
 
       // calotowers
       caloTowerOccMap->Fill(ieta,iphi);
       caloTowerEnergyMap->Fill(ieta,iphi,it->energy());
       caloTowerTime->Fill((it->ecalTime()+it->hcalTime())/2.);
       caloTowerEnergy->Fill(it->energy());
       //std::cout <<"calotower times: "<<it->ecalTime()<<"  "<<it->hcalTime()<<std::endlx;
 
       //hcal
       hcalOccMap->Fill(ieta,iphi);
       hcalEnergyMap->Fill(ieta,iphi,it->hadEnergy());
       hcalTime->Fill(it->hcalTime());
       hcalEnergy->Fill(it->hadEnergy());
 
       //ecal
       ecalOccMap->Fill(ieta,iphi);
       ecalEnergyMap->Fill(ieta,iphi,it->emEnergy());
       ecalTime->Fill(it->ecalTime());
       ecalEnergy->Fill(it->emEnergy());
 
       time_HcalvsEcal->Fill(it->ecalTime(),it->hcalTime());
       energy_HcalvsEcal->Fill(it->emEnergy(),it->hadEnergy());
     } // for CaloTowerCollection::const_iterator it...
 
   // Now form average for each eta bin;
   for (int i=0;i<90;++i)
     {
       if (etacounts[i]==0) continue; 
       caloTowerMeanEnergyEta->Fill(i-42, 1.*(hcalenergy[i]+ecalenergy[i])/etacounts[i]);
       hcalMeanEnergyEta->Fill(i-42, 1.*hcalenergy[i]/etacounts[i]);
       ecalMeanEnergyEta->Fill(i-42,1.*ecalenergy[i]/etacounts[i]);
     }
 }

void HcalCaloTowerMonitor::reset ( void )

Definition at line 265 of file HcalCaloTowerMonitor.cc.

265 {}

void HcalCaloTowerMonitor::setup	(	const edm::ParameterSet &	ps,
		DQMStore *	dbe
	)

virtual

Reimplemented from HcalBaseMonitor.

Definition at line 8 of file HcalCaloTowerMonitor.cc.

References HcalBaseMonitor::baseFolder_, DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookInt(), DQMStore::bookProfile(), caloTowerEnergy, caloTowerEnergyMap, caloTowerMeanEnergyEta, caloTowerOccMap, caloTowerTime, gather_cfg::cout, ecalEnergy, ecalEnergyMap, ecalMeanEnergyEta, ecalOccMap, ecalTime, energy_CaloTowervsEcal, energy_CaloTowervsHcal, energy_HcalvsEcal, etaBins_, etaMax_, etaMin_, MonitorElement::Fill(), edm::ParameterSet::getUntrackedParameter(), hcalEnergy, hcalEnergyMap, hcalMeanEnergyEta, hcalOccMap, hcalTime, ievt_, HcalBaseMonitor::m_dbe, meEVT_, phiBins_, phiMax_, phiMin_, HcalBaseMonitor::rootFolder_, DQMStore::setCurrentFolder(), HcalBaseMonitor::setup(), time_CaloTowervsEcal, time_CaloTowervsHcal, and time_HcalvsEcal.

 {
   HcalBaseMonitor::setup(ps,dbe);
 
   baseFolder_ = rootFolder_+"CaloTowerMonitor";
   etaMax_ = ps.getUntrackedParameter<double>("MaxEta", 41.5);
   etaMin_ = ps.getUntrackedParameter<double>("MinEta", -41.5);
   etaBins_ = (int)(etaMax_ - etaMin_);
   std::cout << "CaloTower Monitor eta min/max set to " << etaMin_ << "/" << etaMax_ << std::endl;
 
   phiMax_ = ps.getUntrackedParameter<double>("MaxPhi", 73);
   phiMin_ = ps.getUntrackedParameter<double>("MinPhi", 0);
   phiBins_ = (int)(phiMax_ - phiMin_);
   std::cout << "CaloTower Monitor phi min/max set to " << phiMin_ << "/" << phiMax_ << std::endl;
 
   ievt_=0;
   // book histograms
   if (m_dbe)
     {
       m_dbe->setCurrentFolder(baseFolder_); 
       meEVT_ = m_dbe->bookInt("CaloTower Event Number");     
       meEVT_->Fill(ievt_); 
 
       // Calo Tower occupancy
       caloTowerOccMap=m_dbe->book2D("CaloTowerOccupancyMap", "Calo Tower Occupancy Map",
                     etaBins_,etaMin_,etaMax_,
                     phiBins_,phiMin_,phiMax_);
       // Calo Tower energy
       caloTowerEnergyMap=m_dbe->book2D("CaloTowerEnergyMap", "Calo Tower Energy Map",
                        etaBins_,etaMin_,etaMax_,
                        phiBins_,phiMin_,phiMax_);
 
       caloTowerTime = m_dbe->book1D("CaloTowerTime", "Calo Tower Time",
                     10,0,10);
       caloTowerEnergy = m_dbe->book1D("CalotowerEnergy", "Calo Tower Energy",
                       100,0,100);
       caloTowerMeanEnergyEta = m_dbe->bookProfile("CaloTowerMeanEnergyEta","Calo tower Mean Energy vs. Eta",
                           etaBins_,etaMin_,etaMax_,100,0,100);
 
       // Make plots of HCAL contributions to tower
       m_dbe->setCurrentFolder(baseFolder_+"/"+"HCAL");
       hcalOccMap=m_dbe->book2D("HcalOccupancyMap", "Calo Tower Occupancy Map",
                    etaBins_,etaMin_,etaMax_,
                    phiBins_,phiMin_,phiMax_);
       // Calo Tower energy
       hcalEnergyMap=m_dbe->book2D("HcalEnergyMap", "Calo Tower Energy Map",
                   etaBins_,etaMin_,etaMax_,
                   phiBins_,phiMin_,phiMax_);
 
       hcalTime = m_dbe->book1D("HcalTime", "Calo Tower Time",
                     10,0,10);
       hcalEnergy = m_dbe->book1D("HcalEnergy", "Calo Tower Energy",
                       100,0,100);
       hcalMeanEnergyEta = m_dbe->bookProfile("HcalMeanEnergyEta",
                          "Calo tower Mean Energy vs. Eta",
                          etaBins_,etaMin_,etaMax_,
                          100,0,100);
       
 
       m_dbe->setCurrentFolder(baseFolder_+"/"+"ECAL");
       ecalOccMap=m_dbe->book2D("EcalOccupancyMap", "Calo Tower Occupancy Map",
                    etaBins_,etaMin_,etaMax_,
                    phiBins_,phiMin_,phiMax_);
       // Calo Tower energy
       ecalEnergyMap=m_dbe->book2D("EcalEnergyMap", "Calo Tower Energy Map",
                   etaBins_,etaMin_,etaMax_,
                   phiBins_,phiMin_,phiMax_);
 
       ecalTime = m_dbe->book1D("EcalTime", "Calo Tower Time",
                    10,0,10);
       ecalEnergy = m_dbe->book1D("EcalEnergy", "Calo Tower Energy",
                  100,0,100);
       ecalMeanEnergyEta = m_dbe->bookProfile("EcalMeanEnergyEta",
                          "Calo tower Mean Energy vs. Eta",
                          etaBins_,etaMin_,etaMax_,
                          100,0,100);
 
       // Comparison Plots
       m_dbe->setCurrentFolder(baseFolder_+"/"+"ComparisonPlots");
       time_HcalvsEcal=m_dbe->book2D("HcalvsEcalTime",
                     "Hcal Time vs. Ecal time",
                     10,0,10,10,0,10);
       time_CaloTowervsEcal=m_dbe->book2D("CaloTowervsEcalTime",
                      "Calotower Time vs. Ecal time",
                      10,0,10,10,0,10);
       time_CaloTowervsHcal=m_dbe->book2D("CaloTowervsHcalTime",
                      "CaloTower Time vs. Hcal time",
                      10,0,10,10,0,10);
       energy_HcalvsEcal=m_dbe->book2D("HcalvsEcalEnergy",
                     "Hcal Energy vs. Ecal energy",
                     100,0,100,100,0,100);
       energy_CaloTowervsEcal=m_dbe->book2D("CaloTowervsEcalEnergy",
                      "Calotower Energy vs. Ecal energy",
                      100,0,100,100,0,100);
       energy_CaloTowervsHcal=m_dbe->book2D("CaloTowervsHcalEnergy",
                      "CaloTower Energy vs. Hcal energy",
                      100,0,100,100,0,100);
     } // if (m_dbe)
   
 } //void HcalCaloTowerMonitor::setup(...)

Member Data Documentation

MonitorElement* HcalCaloTowerMonitor::caloTowerEnergy

private

Definition at line 39 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::caloTowerEnergyMap

private

Definition at line 37 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::caloTowerMeanEnergyEta

private

Definition at line 40 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::caloTowerOccMap

private

Definition at line 36 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::caloTowerTime

private

Definition at line 38 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

bool HcalCaloTowerMonitor::debug_

private

Monitoring elements.

Definition at line 30 of file HcalCaloTowerMonitor.h.

MonitorElement* HcalCaloTowerMonitor::ecalEnergy

private

Definition at line 53 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::ecalEnergyMap

private

Definition at line 51 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::ecalMeanEnergyEta

private

Definition at line 54 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::ecalOccMap

private

Definition at line 50 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::ecalTime

private

Definition at line 52 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::energy_CaloTowervsEcal

private

Definition at line 61 of file HcalCaloTowerMonitor.h.

Referenced by setup().

MonitorElement* HcalCaloTowerMonitor::energy_CaloTowervsHcal

private

Definition at line 62 of file HcalCaloTowerMonitor.h.

Referenced by setup().

MonitorElement* HcalCaloTowerMonitor::energy_HcalvsEcal

private

Definition at line 60 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

int HcalCaloTowerMonitor::etaBins_

private

Definition at line 32 of file HcalCaloTowerMonitor.h.

Referenced by setup().

double HcalCaloTowerMonitor::etaMax_

private

Definition at line 31 of file HcalCaloTowerMonitor.h.

Referenced by setup().

double HcalCaloTowerMonitor::etaMin_

private

Definition at line 31 of file HcalCaloTowerMonitor.h.

Referenced by setup().

MonitorElement* HcalCaloTowerMonitor::hcalEnergy

private

Definition at line 46 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::hcalEnergyMap

private

Definition at line 44 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::hcalMeanEnergyEta

private

Definition at line 47 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::hcalOccMap

private

Definition at line 43 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

MonitorElement* HcalCaloTowerMonitor::hcalTime

private

Definition at line 45 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

int HcalCaloTowerMonitor::ievt_

private

Definition at line 33 of file HcalCaloTowerMonitor.h.

Referenced by setup().

MonitorElement* HcalCaloTowerMonitor::meEVT_

private

Definition at line 64 of file HcalCaloTowerMonitor.h.

Referenced by setup().

int HcalCaloTowerMonitor::phiBins_

private

Definition at line 32 of file HcalCaloTowerMonitor.h.

Referenced by setup().

double HcalCaloTowerMonitor::phiMax_

private

Definition at line 31 of file HcalCaloTowerMonitor.h.

Referenced by setup().

double HcalCaloTowerMonitor::phiMin_

private

Definition at line 31 of file HcalCaloTowerMonitor.h.

Referenced by setup().

MonitorElement* HcalCaloTowerMonitor::time_CaloTowervsEcal

private

Definition at line 58 of file HcalCaloTowerMonitor.h.

Referenced by setup().

MonitorElement* HcalCaloTowerMonitor::time_CaloTowervsHcal

private

Definition at line 59 of file HcalCaloTowerMonitor.h.

Referenced by setup().

MonitorElement* HcalCaloTowerMonitor::time_HcalvsEcal

private

Definition at line 57 of file HcalCaloTowerMonitor.h.

Referenced by processEvent(), and setup().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation