CastorRecHitMonitor Class Reference

#include <CastorRecHitMonitor.h>

Inheritance diagram for CastorRecHitMonitor:

Public Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
	CastorRecHitMonitor (const edm::ParameterSet &ps)
void	processEvent (const CastorRecHitCollection &castorHits)
void	processEventJets (const reco::BasicJetCollection &Jets)
void	processEventTowers (const reco::CastorTowerCollection &castorTowers)
void	setup (const edm::ParameterSet &ps)
	~CastorRecHitMonitor ()
Public Member Functions inherited from CastorBaseMonitor
	CastorBaseMonitor ()
virtual	~CastorBaseMonitor ()

Private Attributes
float	energyInEachChannel [14][16]
TH2F *	h2RecHitMap
MonitorElement *	h2RHchan
MonitorElement *	h2RHentriesMap
MonitorElement *	h2RHmap
MonitorElement *	h2RHoccmap
MonitorElement *	h2RHvsSec
MonitorElement *	h2TowerEMhad
MonitorElement *	hallchan
MonitorElement *	hJetEnergy
MonitorElement *	hJetEta
MonitorElement *	hJetPhi
MonitorElement *	hJetsMultipl
MonitorElement *	hRHtime
MonitorElement *	hTowerDepth
MonitorElement *	hTowerE
MonitorElement *	hTowerMultipl
int	ievt_
std::string	subsystemname

Additional Inherited Members
Protected Attributes inherited from CastorBaseMonitor
std::string	baseFolder_
edm::CPUTimer	cpu_timer
int	fVerbosity
std::string	rootFolder_
bool	showTiming

Detailed Description

Definition at line 21 of file CastorRecHitMonitor.h.

Constructor & Destructor Documentation

CastorRecHitMonitor::CastorRecHitMonitor

(

const edm::ParameterSet &

ps

)

Definition at line 15 of file CastorRecHitMonitor.cc.

References edm::ParameterSet::getUntrackedParameter(), ievt_, AlCaHLTBitMon_QueryRunRegistry::string, and subsystemname.

{
 subsystemname =
    ps.getUntrackedParameter<std::string>("subSystemFolder","Castor");
 ievt_=0; 
}

CastorRecHitMonitor::~CastorRecHitMonitor

(

)

Definition at line 23 of file CastorRecHitMonitor.cc.

{ }

Member Function Documentation

void CastorRecHitMonitor::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

Definition at line 35 of file CastorRecHitMonitor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), gather_cfg::cout, create_public_lumi_plots::exp, CastorBaseMonitor::fVerbosity, MonitorElement::getTH1F(), MonitorElement::getTH2F(), h2RHchan, h2RHentriesMap, h2RHmap, h2RHoccmap, h2RHvsSec, h2TowerEMhad, hallchan, hJetEnergy, hJetEta, hJetPhi, hJetsMultipl, hRHtime, hTowerDepth, hTowerE, hTowerMultipl, i, j, dqm-mbProfile::log, alignCSCRings::s, DQMStore::IBooker::setCurrentFolder(), and subsystemname.

Referenced by CastorMonitorModule::bookHistograms().

{
 char s[60];
 if(fVerbosity>0) 
  std::cout<<"CastorRecHitMonitor::bookHistograms"<<std::endl;
 ibooker.setCurrentFolder(subsystemname + "/CastorRecHitMonitor");

  const int N_Sec = 16;
  const int nySec = 20;
  float ySec[nySec+1];
  float xSec[N_Sec+1];
  double E0sec = 1./1024.;
  ySec[0] = 0.; ySec[1] = E0sec;
  double lnBsec = log(2.);
  for(int j=1; j<nySec; j++) ySec[j+1] = E0sec*exp(j*lnBsec);
  for(int i=0; i<=N_Sec; i++) xSec[i]=i;

  sprintf(s,"CastorRecHit by Sectors");
    h2RHvsSec = ibooker.book2D(s,s, N_Sec, xSec, nySec, ySec);
    h2RHvsSec->getTH2F()->GetXaxis()->SetTitle("sectorPhi");
    h2RHvsSec->getTH2F()->GetYaxis()->SetTitle("RecHit / GeV");
    h2RHvsSec->getTH2F()->SetOption("colz");


 const int nxCh = 224;
 const int nyE = 18;
 float xCh[nxCh+1];
 float yErh[nyE+1];
 for(int i=0; i<=nxCh; i++) xCh[i]=i;
 double E0 = 1./1024.;
 double lnA = log(2.);
 yErh[0] = 0.; yErh[1] = E0;
 for(int j=1; j<nyE; j++) yErh[j+1] = E0*exp(j*lnA);

  sprintf(s,"CastorTileRecHit");
    h2RHchan = ibooker.book2D(s,s, nxCh, xCh, nyE, yErh);
    h2RHchan->getTH2F()->GetXaxis()->SetTitle("sector*14+module");
    h2RHchan->getTH2F()->GetYaxis()->SetTitle("RecHit / GeV");
    h2RHchan->getTH2F()->SetOption("colz");  

  sprintf(s,"Reco all tiles");
   hallchan = ibooker.book1D(s,s,nyE,yErh);
   hallchan->getTH1F()->GetXaxis()->SetTitle("GeV");
 
  sprintf(s,"CastorRecHitMap(cumulative)");
    h2RHmap = ibooker.book2D(s,s,14, 0,14, 16, 0,16);
    h2RHmap->getTH2F()->GetXaxis()->SetTitle("moduleZ");  
    h2RHmap->getTH2F()->GetYaxis()->SetTitle("sectorPhi");
    h2RHmap->getTH2F()->SetOption("colz");

  sprintf(s,"CastorRecHitOccMap");
    h2RHoccmap = ibooker.book2D(s,s,14, 0,14, 16, 0,16);
    h2RHoccmap->getTH2F()->GetXaxis()->SetTitle("moduleZ");
    h2RHoccmap->getTH2F()->GetYaxis()->SetTitle("sectorPhi");
    h2RHoccmap->getTH2F()->SetOption("colz");

  sprintf(s,"CastorRecHitEntriesMap");
    h2RHentriesMap = ibooker.book2D(s,s,14, 0,14, 16, 0,16);
    h2RHentriesMap->getTH2F()->GetXaxis()->SetTitle("moduleZ");
    h2RHentriesMap->getTH2F()->GetYaxis()->SetTitle("sectorPhi");
    h2RHentriesMap->getTH2F()->SetOption("colz");

  sprintf(s,"CastorRecHitTime");
    hRHtime = ibooker.book1D(s,s,301, -101.,200.);

  sprintf(s,"CASTORTowerDepth");
   hTowerDepth = ibooker.book1D(s,s,130,-15500.,-14200.);
   hTowerDepth->getTH1F()->GetXaxis()->SetTitle("mm");

  sprintf(s,"CASTORTowerMultiplicity");
   hTowerMultipl = ibooker.book1D(s,s,20,0.,20.);

 const int NEtow = 20;
  float EhadTow[NEtow+1];
  float EMTow[NEtow+1];
  float ETower[NEtow+2];
  double E0tow = 1./1024.;
  EMTow[0] = 0.; EMTow[1] = E0tow;
  EhadTow[0] = 0.; EhadTow[1] = E0tow;
  ETower[0] = 0.; ETower[1] = E0tow;
  double lnBtow = log(2.);
  for(int j=1; j<NEtow; j++) EMTow[j+1] = E0tow*exp(j*lnBtow);
  for(int j=1; j<NEtow; j++) EhadTow[j+1] = E0tow*exp(j*lnBtow);
  for(int j=1; j<=NEtow; j++) ETower[j+1] = E0tow*exp(j*lnBtow);

  sprintf(s,"CASTORTowerEMvsEhad");
    h2TowerEMhad = ibooker.book2D(s,s, NEtow, EhadTow, NEtow, EMTow);
    h2TowerEMhad->getTH2F()->GetXaxis()->SetTitle("Ehad / GeV");
    h2TowerEMhad->getTH2F()->GetYaxis()->SetTitle("EM / GeV");
    h2TowerEMhad->getTH2F()->SetOption("colz");

  sprintf(s,"CASTORTowerTotalEnergy");
   hTowerE = ibooker.book1D(s,s,NEtow+1,ETower);
   hTowerE->getTH1F()->GetXaxis()->SetTitle("GeV");

  sprintf(s,"CASTORJetsMultiplicity");
   hJetsMultipl = ibooker.book1D(s,s,16, 0.,16.);

  sprintf(s,"CASTORJetEnergy");
   hJetEnergy = ibooker.book1D(s,s,5000, 0.,500.);

  sprintf(s,"CASTORJetEta");
   hJetEta = ibooker.book1D(s,s,126, -6.3, 6.3);

  sprintf(s,"CASTORJetPhi");
   hJetPhi = ibooker.book1D(s,s,63, -3.15,3.15);

  if(fVerbosity>0) 
    std::cout<<"CastorRecHitMonitor::bookHistograms(end)"<<std::endl;
  return;
}

void CastorRecHitMonitor::processEvent

(

const CastorRecHitCollection &

castorHits

)

Definition at line 165 of file CastorRecHitMonitor.cc.

References a, edm::SortedCollection< T, SORT >::begin(), gather_cfg::cout, CastorBaseMonitor::cpu_timer, edm::SortedCollection< T, SORT >::end(), relval_parameters_module::energy, energyInEachChannel, MonitorElement::Fill(), CastorBaseMonitor::fVerbosity, MonitorElement::getTH2F(), h2RHchan, h2RHentriesMap, h2RHmap, h2RHoccmap, h2RHvsSec, hallchan, hRHtime, ievt_, mod(), python.rootplot.argparse::module, phi, edm::CPUTimer::reset(), CastorBaseMonitor::showTiming, edm::SortedCollection< T, SORT >::size(), edm::CPUTimer::start(), and z.

Referenced by CastorMonitorModule::analyze().

{
 if(fVerbosity>0) std::cout << "CastorRecHitMonitor::processEvent (begin)"<< std::endl;
 ievt_++; 
 for (int z=0; z<14; z++) for (int phi=0; phi<16; phi++)
    energyInEachChannel[z][phi] = 0.;

 CastorRecHitCollection::const_iterator CASTORiter;
 if (showTiming)  { cpu_timer.reset(); cpu_timer.start(); } 

 if(castorHits.size() <= 0) return;

 //for(edm::TriggerResults::const_iterator iTrig= hltResults->begin();
//  iTrig!= hltResults->end(); iTrig++) {;}

 for(CASTORiter=castorHits.begin(); CASTORiter!=castorHits.end(); ++CASTORiter)
 { 
   float energy = CASTORiter->energy();    
   float time = CASTORiter->time();
   float time2 = time;
   if(time < -100.) time2 = -100.;
   hRHtime->Fill(time2);

   HcalCastorDetId id(CASTORiter->detid().rawId());
      //float zside  = id.zside(); 
   int module = (int)id.module(); //-- get module
   int sector = (int)id.sector(); //-- get sector 

   energyInEachChannel[module-1][sector-1] += energy;

   h2RHentriesMap->Fill(module-1,sector-1);
 } // end for(CASTORiter=castorHits.begin(); CASTORiter!= ...

  double etot = 0.;
  for(int phi=0; phi<16; phi++) {
    double es = 0.;
    for (int z=0; z<14; z++) {
      float rh = energyInEachChannel[z][phi]*0.001;
      int ind = phi*14 + z +1;
      h2RHchan->Fill(ind,rh);
      hallchan->Fill(rh);
      if(rh < 0.) continue;      
      h2RHmap->Fill(z,phi,rh); 
      es += rh;
    }
    h2RHvsSec->Fill(phi,es);
    etot += es;
  } // end for(int phi=0;

 if(ievt_ %100 == 0) 
  for(int mod=1; mod<=14; mod++) for(int sec=1; sec<=16;sec++) {
    double a= h2RHmap->getTH2F()->GetBinContent(mod,sec);
    h2RHoccmap->getTH2F()->SetBinContent(mod,sec,a/double(ievt_));
  }

  if(fVerbosity>0) std::cout << "CastorRecHitMonitor::processEvent (end)"<< std::endl;
  return;
}

void CastorRecHitMonitor::processEventJets

(

const reco::BasicJetCollection &

Jets

)

Definition at line 224 of file CastorRecHitMonitor.cc.

References relval_parameters_module::energy, MonitorElement::Fill(), hJetEnergy, hJetEta, hJetPhi, and hJetsMultipl.

Referenced by CastorMonitorModule::analyze().

{
 int nJets=0;
 for(reco::BasicJetCollection::const_iterator ibegin = Jets.begin(),
  iend = Jets.end(), ijet = ibegin; ijet!= iend; ++ijet) {
  nJets++;
  float energy = ijet->energy()*0.001;
  hJetEnergy->Fill(energy);
  hJetEta->Fill(ijet->eta());
  hJetPhi->Fill(ijet->phi());
 }
 hJetsMultipl->Fill(nJets);
}

void CastorRecHitMonitor::processEventTowers

(

const reco::CastorTowerCollection &

castorTowers

)

Definition at line 148 of file CastorRecHitMonitor.cc.

References MonitorElement::Fill(), h2TowerEMhad, hTowerDepth, hTowerE, hTowerMultipl, and ecaldqm::nTowers.

Referenced by CastorMonitorModule::analyze().

{
 if(castorTowers.size() <= 0) return;
 int nTowers = 0;

 for(reco::CastorTowerCollection::const_iterator iTower= castorTowers.begin();
  iTower!= castorTowers.end(); iTower++) {

   hTowerE->Fill(iTower->energy()*0.001);
   h2TowerEMhad->Fill(iTower->hadEnergy()*0.001,iTower->emEnergy()*0.001);
   hTowerDepth->Fill(iTower->depth());
   nTowers++;
 }
 hTowerMultipl->Fill(nTowers);
}

void CastorRecHitMonitor::setup ( const edm::ParameterSet &  ps )

virtual

Reimplemented from CastorBaseMonitor.

Definition at line 27 of file CastorRecHitMonitor.cc.

Referenced by CastorMonitorModule::CastorMonitorModule().

{
//  CastorBaseMonitor::setup(ps);
  return;
}

Member Data Documentation

float CastorRecHitMonitor::energyInEachChannel[14][16]

private

Definition at line 33 of file CastorRecHitMonitor.h.

Referenced by processEvent().

TH2F* CastorRecHitMonitor::h2RecHitMap

private

Definition at line 45 of file CastorRecHitMonitor.h.

MonitorElement* CastorRecHitMonitor::h2RHchan

private

Definition at line 46 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorRecHitMonitor::h2RHentriesMap

private

Definition at line 50 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorRecHitMonitor::h2RHmap

private

Definition at line 48 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorRecHitMonitor::h2RHoccmap

private

Definition at line 49 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorRecHitMonitor::h2RHvsSec

private

Definition at line 47 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorRecHitMonitor::h2TowerEMhad

private

Definition at line 38 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEventTowers().

MonitorElement* CastorRecHitMonitor::hallchan

private

Definition at line 52 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorRecHitMonitor::hJetEnergy

private

Definition at line 41 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEventJets().

MonitorElement* CastorRecHitMonitor::hJetEta

private

Definition at line 42 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEventJets().

MonitorElement* CastorRecHitMonitor::hJetPhi

private

Definition at line 43 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEventJets().

MonitorElement* CastorRecHitMonitor::hJetsMultipl

private

Definition at line 40 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEventJets().

MonitorElement* CastorRecHitMonitor::hRHtime

private

Definition at line 51 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEvent().

MonitorElement* CastorRecHitMonitor::hTowerDepth

private

Definition at line 37 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEventTowers().

MonitorElement* CastorRecHitMonitor::hTowerE

private

Definition at line 36 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEventTowers().

MonitorElement* CastorRecHitMonitor::hTowerMultipl

private

Definition at line 39 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and processEventTowers().

int CastorRecHitMonitor::ievt_

private

Definition at line 32 of file CastorRecHitMonitor.h.

Referenced by CastorRecHitMonitor(), and processEvent().

std::string CastorRecHitMonitor::subsystemname

private

Definition at line 34 of file CastorRecHitMonitor.h.

Referenced by bookHistograms(), and CastorRecHitMonitor().