#include <CastorRecHitMonitor.h>

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)

	~CastorRecHitMonitor ()

Private Attributes
float	energyInEachChannel [14][16]

int	fVerbosity = 0

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

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 gather_cfg::cout, fVerbosity, edm::ParameterSet::getUntrackedParameter(), ievt_, AlCaHLTBitMon_QueryRunRegistry::string, and subsystemname.

 {
  fVerbosity = ps.getUntrackedParameter<int>("debug",0);
  if(fVerbosity>0)
    std::cout<<"CastorRecHitMonitor Constructor: "<<this<<std::endl;
  subsystemname =
     ps.getUntrackedParameter<std::string>("subSystemFolder","Castor");
  ievt_=0; 
 }

CastorRecHitMonitor::~CastorRecHitMonitor ( )

Definition at line 25 of file CastorRecHitMonitor.cc.

25 { }

Member Function Documentation

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

Definition at line 27 of file CastorRecHitMonitor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), gather_cfg::cout, JetChargeProducer_cfi::exp, fVerbosity, MonitorElement::getTH1F(), MonitorElement::getTH2F(), h2RHchan, h2RHentriesMap, h2RHmap, h2RHoccmap, h2RHvsSec, h2TowerEMhad, hallchan, hJetEnergy, hJetEta, hJetPhi, hJetsMultipl, hRHtime, hTowerDepth, hTowerE, hTowerMultipl, mps_fire::i, cmsBatch::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 157 of file CastorRecHitMonitor.cc.

References a, edm::SortedCollection< T, SORT >::begin(), gather_cfg::cout, edm::SortedCollection< T, SORT >::empty(), edm::SortedCollection< T, SORT >::end(), energyInEachChannel, MonitorElement::Fill(), fVerbosity, MonitorElement::getTH2F(), h2RHchan, h2RHentriesMap, h2RHmap, h2RHoccmap, h2RHvsSec, hallchan, hRHtime, triggerObjects_cff::id, ievt_, createfilelist::int, mod(), python.rootplot.argparse::module, phi, ntuplemaker::time, 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.empty()) 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 217 of file CastorRecHitMonitor.cc.

References 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 140 of file CastorRecHitMonitor.cc.

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

Referenced by CastorMonitorModule::analyze().

 {
  if(castorTowers.empty()) 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);
 }