#include <EETimingTask.h>

Inheritance diagram for EETimingTask:

Public Member Functions
	EETimingTask (const edm::ParameterSet &ps)
	Constructor. More...

virtual	~EETimingTask ()
	Destructor. More...

Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()

std::string	workerType () const

virtual	~EDAnalyzer ()

Static Public Attributes
static const float	shiftProf2D = 50.

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

void	beginJob (void)
	BeginJob. More...

void	beginRun (const edm::Run &r, const edm::EventSetup &c)
	BeginRun. More...

void	cleanup (void)
	Cleanup. More...

void	endJob (void)
	EndJob. More...

void	endRun (const edm::Run &r, const edm::EventSetup &c)
	EndRun. More...

void	reset (void)
	Reset. More...

void	setup (void)
	Setup. More...

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Attributes
DQMStore *	dqmStore_

edm::InputTag	EcalRawDataCollection_

edm::InputTag	EcalRecHitCollection_

bool	enableCleanup_

float	energyThreshold_

int	ievt_

bool	init_

bool	initCaloGeometry_

bool	mergeRuns_

MonitorElement *	meTime_ [18]

MonitorElement *	meTimeAmpli_ [18]

MonitorElement *	meTimeAmpliSummary_ [2]

MonitorElement *	meTimeDelta2D_

MonitorElement *	meTimeDelta_

MonitorElement *	meTimeMap_ [18]

MonitorElement *	meTimeSummary1D_ [2]

MonitorElement *	meTimeSummaryMap_ [2]

edm::ESHandle< CaloGeometry >	pGeometry_

std::string	prefixME_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

typedef WorkerT< EDAnalyzer >	WorkerType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 26 of file EETimingTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 33 of file EETimingTask.cc.

References dqmStore_, EcalRawDataCollection_, EcalRecHitCollection_, enableCleanup_, energyThreshold_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), i, init_, initCaloGeometry_, mergeRuns_, meTime_, meTimeAmpli_, meTimeAmpliSummary_, meTimeDelta2D_, meTimeDelta_, meTimeMap_, meTimeSummary1D_, meTimeSummaryMap_, cmsCodeRules.cppFunctionSkipper::operator, and prefixME_.

                                                    {
 
   init_ = false;
 
   initCaloGeometry_ = false;
 
   dqmStore_ = edm::Service<DQMStore>().operator->();
 
   prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
 
   enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
 
   mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
 
   energyThreshold_ = ps.getUntrackedParameter<double>("energyThreshold", 3.);
 
   EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection");
   EcalRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalRecHitCollection");
 
   for (int i = 0; i < 18; i++) {
     meTime_[i] = 0;
     meTimeMap_[i] = 0;
     meTimeAmpli_[i] = 0;
   }
 
   for (int i = 0; i < 2; i++) {
     meTimeAmpliSummary_[i] = 0;
     meTimeSummary1D_[i] = 0;
     meTimeSummaryMap_[i] = 0;
   }
 
   meTimeDelta_ = 0;
   meTimeDelta2D_ = 0;
 
 }

EETimingTask::~EETimingTask ( )

virtual

Destructor.

Definition at line 69 of file EETimingTask.cc.

                            {
 
 }

Member Function Documentation

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

protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 258 of file EETimingTask.cc.

References EcalDCCHeaderBlock::COSMIC, EcalDCCHeaderBlock::COSMICS_GLOBAL, EcalDCCHeaderBlock::COSMICS_LOCAL, EcalEndcap, EcalRawDataCollection_, EcalRecHitCollection_, energyThreshold_, MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByLabel(), i, ievt_, init_, ecalpyutils::ism(), Numbers::iSM(), EEDetId::ix(), EcalRecHit::kGood, EcalRecHit::kOutOfTime, EcalSeverityLevelAlgo::kWeird, LogDebug, meTime_, meTimeAmpli_, meTimeAmpliSummary_, meTimeDelta2D_, meTimeDelta_, meTimeMap_, meTimeSummary1D_, meTimeSummaryMap_, EcalDCCHeaderBlock::MTCC, EcalDCCHeaderBlock::PHYSICS_GLOBAL, EcalDCCHeaderBlock::PHYSICS_LOCAL, setup(), shiftProf2D, and Numbers::subDet().

                                                                    {
 
   bool isData = true;
   bool enable = false;
   int runType[18];
   for (int i=0; i<18; i++) runType[i] = -1;
 
   edm::Handle<EcalRawDataCollection> dcchs;
 
   if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) {
 
     for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) {
 
       if ( Numbers::subDet( *dcchItr ) != EcalEndcap ) continue;
 
       int ism = Numbers::iSM( *dcchItr, EcalEndcap );
 
       runType[ism-1] = dcchItr->getRunType();
 
       if ( dcchItr->getRunType() == EcalDCCHeaderBlock::COSMIC ||
            dcchItr->getRunType() == EcalDCCHeaderBlock::MTCC ||
            dcchItr->getRunType() == EcalDCCHeaderBlock::COSMICS_GLOBAL ||
            dcchItr->getRunType() == EcalDCCHeaderBlock::PHYSICS_GLOBAL ||
            dcchItr->getRunType() == EcalDCCHeaderBlock::COSMICS_LOCAL ||
            dcchItr->getRunType() == EcalDCCHeaderBlock::PHYSICS_LOCAL ) enable = true;
 
     }
 
   } else {
 
     isData = false; enable = true;
     edm::LogWarning("EETimingTask") << EcalRawDataCollection_ << " not available";
 
   }
 
   if ( ! enable ) return;
 
   if ( ! init_ ) this->setup();
 
   ievt_++;
 
   float sumTime_hithr[2] = {0.,0.};
   int n_hithr[2] = {0,0};
 
   edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
   c.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
 
   edm::Handle<EcalRecHitCollection> hits;
 
   if ( e.getByLabel(EcalRecHitCollection_, hits) ) {
 
     int neh = hits->size();
     LogDebug("EETimingTask") << "event " << ievt_ << " hits collection size " << neh;
 
     for ( EcalRecHitCollection::const_iterator hitItr = hits->begin(); hitItr != hits->end(); ++hitItr ) {
 
       EEDetId id = hitItr->id();
 
       int ix = id.ix();
       int iy = id.iy();
       int iz = ( id.positiveZ() ) ? 1 : 0;
 
       int ism = Numbers::iSM( id );
 
       if ( ism >= 1 && ism <= 9 ) ix = 101 - ix;
 
       float xix = ix - 0.5;
       float xiy = iy - 0.5;
 
       if ( isData ) {
 
         if ( ! ( runType[ism-1] == EcalDCCHeaderBlock::COSMIC ||
                  runType[ism-1] == EcalDCCHeaderBlock::MTCC ||
                  runType[ism-1] == EcalDCCHeaderBlock::COSMICS_GLOBAL ||
                  runType[ism-1] == EcalDCCHeaderBlock::PHYSICS_GLOBAL ||
                  runType[ism-1] == EcalDCCHeaderBlock::COSMICS_LOCAL ||
                  runType[ism-1] == EcalDCCHeaderBlock::PHYSICS_LOCAL ) ) continue;
 
       }
 
       MonitorElement* meTime = 0;
       MonitorElement* meTimeMap = 0;
       MonitorElement* meTimeAmpli = 0;
 
       meTime = meTime_[ism-1];
       meTimeMap = meTimeMap_[ism-1];
       meTimeAmpli = meTimeAmpli_[ism-1];
 
       float xval = hitItr->energy();
       float yval = hitItr->time();
 
       uint32_t flag = hitItr->recoFlag();
 
       uint32_t sev = sevlv->severityLevel(id, *hits );
 
       if ( (flag == EcalRecHit::kGood || flag == EcalRecHit::kOutOfTime) && sev != EcalSeverityLevelAlgo::kWeird ) {
         if ( meTimeAmpli ) meTimeAmpli->Fill(xval, yval);
         if ( meTimeAmpliSummary_[iz] ) meTimeAmpliSummary_[iz]->Fill(xval, yval);
         if ( hitItr->energy() > energyThreshold_ ) {
           if ( meTimeMap ) meTimeMap->Fill(xix, xiy, yval+shiftProf2D);
           if ( meTime ) meTime->Fill(yval);
           if ( meTimeSummary1D_[iz] ) meTimeSummary1D_[iz]->Fill(yval);
 
           if ( meTimeSummaryMap_[iz] ) meTimeSummaryMap_[iz]->Fill(id.ix()-0.5, xiy, yval+shiftProf2D);
 
           sumTime_hithr[iz] += yval;
           n_hithr[iz]++;
         }
       } // good rh for timing
     } // loop over rh
 
     if (n_hithr[0] > 0 && n_hithr[1] > 0 ) {
       if ( meTimeDelta_ ) meTimeDelta_->Fill( sumTime_hithr[1]/n_hithr[1] - sumTime_hithr[0]/n_hithr[0] );
       if ( meTimeDelta2D_ ) meTimeDelta2D_->Fill( sumTime_hithr[1]/n_hithr[1], sumTime_hithr[0]/n_hithr[0] );
     }
 
   } else {
 
     edm::LogWarning("EETimingTask") << EcalRecHitCollection_ << " not available";
 
   }
 
 }

void EETimingTask::beginJob ( void )

protectedvirtual

BeginJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 73 of file EETimingTask.cc.

References dqmStore_, ievt_, prefixME_, DQMStore::rmdir(), and DQMStore::setCurrentFolder().

                                {
 
   ievt_ = 0;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EETimingTask");
     dqmStore_->rmdir(prefixME_ + "/EETimingTask");
   }
 
 }

void EETimingTask::beginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

BeginRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 84 of file EETimingTask.cc.

References edm::EventSetup::get(), initCaloGeometry_, Numbers::initGeometry(), mergeRuns_, pGeometry_, and reset().

                                                                    {
 
   Numbers::initGeometry(c, true);
 
   if( !initCaloGeometry_ ) {
     c.get<CaloGeometryRecord>().get(pGeometry_);
     initCaloGeometry_ = true;
   }
 
   if ( ! mergeRuns_ ) this->reset();
 
 }

void EETimingTask::cleanup ( void )

protected

Cleanup.

Definition at line 208 of file EETimingTask.cc.

References dqmStore_, edm::getName(), MonitorElement::getName(), i, init_, meTime_, meTimeAmpli_, meTimeAmpliSummary_, meTimeDelta2D_, meTimeDelta_, meTimeMap_, meTimeSummary1D_, meTimeSummaryMap_, prefixME_, DQMStore::removeElement(), and DQMStore::setCurrentFolder().

Referenced by endJob().

                               {
 
   if ( ! init_ ) return;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EETimingTask");
 
     for ( int i = 0; i < 18; i++ ) {
       if ( meTime_[i] ) dqmStore_->removeElement( meTime_[i]->getName() );
       meTime_[i] = 0;
 
       if ( meTimeMap_[i] ) dqmStore_->removeElement( meTimeMap_[i]->getName() );
       meTimeMap_[i] = 0;
 
       if ( meTimeAmpli_[i] ) dqmStore_->removeElement( meTimeAmpli_[i]->getName() );
       meTimeAmpli_[i] = 0;
     }
 
     for (int i = 0; i < 2; i++) {
       if ( meTimeAmpliSummary_[i] ) dqmStore_->removeElement( meTimeAmpliSummary_[i]->getName() );
       meTimeAmpliSummary_[i] = 0;
 
       if ( meTimeSummary1D_[i] ) dqmStore_->removeElement( meTimeSummary1D_[i]->getName() );
       meTimeSummary1D_[i] = 0;
 
       if ( meTimeSummaryMap_[i] ) dqmStore_->removeElement( meTimeSummaryMap_[i]->getName() );
       meTimeSummaryMap_[i] = 0;
 
     }
 
     if ( meTimeDelta_ ) dqmStore_->removeElement( meTimeDelta_->getName() );
     meTimeDelta_ = 0;
 
     if ( meTimeDelta2D_ ) dqmStore_->removeElement( meTimeDelta2D_->getName() );
     meTimeDelta2D_ = 0;
 
   }
 
   init_ = false;
 
 }

void EETimingTask::endJob ( void )

protectedvirtual

EndJob.

Reimplemented from edm::EDAnalyzer.

Definition at line 250 of file EETimingTask.cc.

References cleanup(), enableCleanup_, and ievt_.

                              {
 
   edm::LogInfo("EETimingTask") << "analyzed " << ievt_ << " events";
 
   if ( enableCleanup_ ) this->cleanup();
 
 }

void EETimingTask::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

protectedvirtual

EndRun.

Reimplemented from edm::EDAnalyzer.

Definition at line 97 of file EETimingTask.cc.

                                                                  {
 
 }

void EETimingTask::reset ( void )

protected

Reset.

Definition at line 101 of file EETimingTask.cc.

References i, meTime_, meTimeAmpli_, meTimeAmpliSummary_, meTimeDelta2D_, meTimeDelta_, meTimeMap_, meTimeSummary1D_, meTimeSummaryMap_, and MonitorElement::Reset().

Referenced by beginRun().

                              {
 
   for (int i = 0; i < 18; i++) {
     if ( meTime_[i] ) meTime_[i]->Reset();
     if ( meTimeMap_[i] ) meTimeMap_[i]->Reset();
     if ( meTimeAmpli_[i] ) meTimeAmpli_[i]->Reset();
   }
 
   for (int i = 0; i < 2; i++) {
     if ( meTimeAmpliSummary_[i] ) meTimeAmpliSummary_[i]->Reset();
     if ( meTimeSummary1D_[i] ) meTimeSummary1D_[i]->Reset();
     if ( meTimeSummaryMap_[i] ) meTimeSummaryMap_[i]->Reset();
   }
 
   if ( meTimeDelta_ ) meTimeDelta_->Reset();
   if ( meTimeDelta2D_ ) meTimeDelta2D_->Reset();
 
 }

void EETimingTask::setup ( void )

protected

Setup.

Definition at line 120 of file EETimingTask.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile2D(), dqmStore_, i, init_, Numbers::ix0EE(), Numbers::iy0EE(), meTime_, meTimeAmpli_, meTimeAmpliSummary_, meTimeDelta2D_, meTimeDelta_, meTimeMap_, meTimeSummary1D_, meTimeSummaryMap_, mergeVDriftHistosByStation::name, funct::pow(), prefixME_, Numbers::sEE(), MonitorElement::setAxisTitle(), DQMStore::setCurrentFolder(), shiftProf2D, and DQMStore::tag().

Referenced by analyze().

                             {
 
   init_ = true;
 
   std::string name;
 
   //for timing vs amplitude plots
   const int nbinsE = 25;
   const float minlogE = -0.5;
   const float maxlogE = 2.;
   float binEdgesE[nbinsE + 1];
   for(int i = 0; i <= nbinsE; i++)
     binEdgesE[i] = std::pow(10., minlogE + (maxlogE - minlogE) / nbinsE * i);
 
   const int nbinsT = 200;
   const float minT = -50.;
   const float maxT = 50.;
   float binEdgesT[nbinsT + 1];
   for(int i = 0; i <= nbinsT; i++)
     binEdgesT[i] = minT + (maxT - minT) / nbinsT * i;
 
   if ( dqmStore_ ) {
     dqmStore_->setCurrentFolder(prefixME_ + "/EETimingTask");
 
     for (int i = 0; i < 18; i++) {
       name = "EETMT timing 1D " + Numbers::sEE(i+1);
       meTime_[i] = dqmStore_->book1D(name, name, 50, -25., 25.);
       meTime_[i]->setAxisTitle("time (ns)", 1);
       dqmStore_->tag(meTime_[i], i+1);
 
       name = "EETMT timing " + Numbers::sEE(i+1);
       meTimeMap_[i] = dqmStore_->bookProfile2D(name, name, 50, Numbers::ix0EE(i+1)+0., Numbers::ix0EE(i+1)+50., 50, Numbers::iy0EE(i+1)+0., Numbers::iy0EE(i+1)+50., -20.+shiftProf2D, 20.+shiftProf2D, "s");
       meTimeMap_[i]->setAxisTitle("ix", 1);
       if ( i+1 >= 1 && i+1 <= 9 ) meTimeMap_[i]->setAxisTitle("101-ix", 1);
       meTimeMap_[i]->setAxisTitle("iy", 2);
       meTimeMap_[i]->setAxisTitle("time (ns)", 3);
       dqmStore_->tag(meTimeMap_[i], i+1);
 
       name = "EETMT timing vs amplitude " + Numbers::sEE(i+1);
       meTimeAmpli_[i] = dqmStore_->book2D(name, name, nbinsE, binEdgesE, nbinsT, binEdgesT);
       meTimeAmpli_[i]->setAxisTitle("energy (GeV)", 1);
       meTimeAmpli_[i]->setAxisTitle("time (ns)", 2);
       dqmStore_->tag(meTimeAmpli_[i], i+1);
     }
 
     name = "EETMT timing vs amplitude summary EE -";
     meTimeAmpliSummary_[0] = dqmStore_->book2D(name, name, nbinsE, binEdgesE, nbinsT, binEdgesT);
     meTimeAmpliSummary_[0]->setAxisTitle("energy (GeV)", 1);
     meTimeAmpliSummary_[0]->setAxisTitle("time (ns)", 2);
 
     name = "EETMT timing vs amplitude summary EE +";
     meTimeAmpliSummary_[1] = dqmStore_->book2D(name, name, nbinsE, binEdgesE, nbinsT, binEdgesT);
     meTimeAmpliSummary_[1]->setAxisTitle("energy (GeV)", 1);
     meTimeAmpliSummary_[1]->setAxisTitle("time (ns)", 2);
 
     name = "EETMT timing 1D summary EE -";
     meTimeSummary1D_[0] = dqmStore_->book1D(name, name, 50, -25., 25.);
     meTimeSummary1D_[0]->setAxisTitle("time (ns)", 1);
 
     name = "EETMT timing 1D summary EE +";
     meTimeSummary1D_[1] = dqmStore_->book1D(name, name, 50, -25., 25.);
     meTimeSummary1D_[1]->setAxisTitle("time (ns)", 1);
 
     name = "EETMT timing map EE -";
     meTimeSummaryMap_[0] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., -20.+shiftProf2D, 20.+shiftProf2D, "s");
     meTimeSummaryMap_[0]->setAxisTitle("ix'", 1);
     meTimeSummaryMap_[0]->setAxisTitle("101-iy'", 2);
     meTimeSummaryMap_[0]->setAxisTitle("time (ns)", 3);
 
     name = "EETMT timing map EE +";
     meTimeSummaryMap_[1] = dqmStore_->bookProfile2D(name, name, 20, 0., 100., 20, 0., 100., -20.+shiftProf2D, 20.+shiftProf2D, "s");
     meTimeSummaryMap_[1]->setAxisTitle("ix'", 1);
     meTimeSummaryMap_[1]->setAxisTitle("iy'", 2);
     meTimeSummaryMap_[1]->setAxisTitle("time (ns)", 3);
 
     name = "EETMT timing EE+ - EE-";
     meTimeDelta_ = dqmStore_->book1D(name, name, 100, -3., 3.);
     meTimeDelta_->setAxisTitle("time (ns)", 1);
 
     name = "EETMT timing EE+ vs EE-";
     meTimeDelta2D_ = dqmStore_->book2D(name, name, 50, -25., 25., 50, -25., 25.);
     meTimeDelta2D_->setAxisTitle("EE+ average time (ns)", 1);
     meTimeDelta2D_->setAxisTitle("EE- average time (ns)", 2);
 
   }
 
 }