#include <EBClusterTask.h>
Definition at line 20 of file EBClusterTask.h.
EBClusterTask::EBClusterTask | ( | const edm::ParameterSet & | ps | ) |
Constructor.
Definition at line 44 of file EBClusterTask.cc.
References BasicClusterCollection_, dqmStore_, EcalRawDataCollection_, EcalRecHitCollection_, enableCleanup_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), init_, meBCEne_, meBCEneMap_, meBCEneMapProjEta_, meBCEneMapProjPhi_, meBCETMap_, meBCETMapProjEta_, meBCETMapProjPhi_, meBCNum_, meBCNumMap_, meBCNumMapProjEta_, meBCNumMapProjPhi_, meBCSiz_, meBCSizMap_, meBCSizMapProjEta_, meBCSizMapProjPhi_, meInvMassHigh_, meInvMassHighSel_, meInvMassJPsi_, meInvMassJPsiSel_, meInvMassPi0_, meInvMassPi0Sel_, meInvMassZ0_, meInvMassZ0Sel_, mergeRuns_, mes1s9_, mes1s9thr_, mes9s25_, meSCCrystalSiz_, meSCEne2_, meSCEne_, meSCEneLowScale_, meSCEneVsEMax_, meSCMapSingleCrystal_, meSCNum_, meSCSeedEne_, meSCSeedMapOcc_, meSCSiz_, cmsCodeRules::cppFunctionSkipper::operator, prefixME_, SuperClusterCollection_, thrCandEt_, thrClusEt_, and thrS4S9_.
{ init_ = false; dqmStore_ = edm::Service<DQMStore>().operator->(); prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", ""); enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false); mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false); // parameters... EcalRawDataCollection_ = ps.getParameter<edm::InputTag>("EcalRawDataCollection"); BasicClusterCollection_ = ps.getParameter<edm::InputTag>("BasicClusterCollection"); SuperClusterCollection_ = ps.getParameter<edm::InputTag>("SuperClusterCollection"); EcalRecHitCollection_ = ps.getParameter<edm::InputTag>("EcalRecHitCollection"); // histograms... meBCEne_ = 0; meBCNum_ = 0; meBCSiz_ = 0; meBCEneMap_ = 0; meBCNumMap_ = 0; meBCETMap_ = 0; meBCSizMap_ = 0; meBCEneMapProjEta_ = 0; meBCNumMapProjEta_ = 0; meBCETMapProjEta_ = 0; meBCSizMapProjEta_ = 0; meBCEneMapProjPhi_ = 0; meBCNumMapProjPhi_ = 0; meBCETMapProjPhi_ = 0; meBCSizMapProjPhi_ = 0; meSCEne_ = 0; meSCNum_ = 0; meSCSiz_ = 0; meSCCrystalSiz_ = 0; meSCSeedEne_ = 0; meSCEne2_ = 0; meSCEneVsEMax_ = 0; meSCEneLowScale_ = 0; meSCSeedMapOcc_ = 0; meSCMapSingleCrystal_ = 0; mes1s9_ = 0; mes1s9thr_ = 0; mes9s25_ = 0; meInvMassPi0_ = 0; meInvMassJPsi_ = 0; meInvMassZ0_ = 0; meInvMassHigh_ = 0; meInvMassPi0Sel_ = 0; meInvMassJPsiSel_ = 0; meInvMassZ0Sel_ = 0; meInvMassHighSel_ = 0; thrS4S9_ = 0.85; thrClusEt_ = 0.200; thrCandEt_ = 0.650; }
EBClusterTask::~EBClusterTask | ( | ) | [virtual] |
void EBClusterTask::analyze | ( | const edm::Event & | e, |
const edm::EventSetup & | c | ||
) | [protected, virtual] |
Analyze.
Implements edm::EDAnalyzer.
Definition at line 518 of file EBClusterTask.cc.
References abs, BasicClusterCollection_, EcalDCCHeaderBlock::BEAMH2, EcalDCCHeaderBlock::BEAMH4, edm::SortedCollection< T, SORT >::begin(), EcalDCCHeaderBlock::COSMIC, EcalDCCHeaderBlock::COSMICS_GLOBAL, EcalDCCHeaderBlock::COSMICS_LOCAL, EcalClusterTools::e2x2(), EcalClusterTools::e3x3(), EcalClusterTools::e5x5(), DetId::Ecal, EcalBarrel, EcalRawDataCollection_, EcalRecHitCollection_, jptDQMConfig_cff::eMax, edm::SortedCollection< T, SORT >::end(), MonitorElement::Fill(), edm::SortedCollection< T, SORT >::find(), edm::EventSetup::get(), edm::Event::getByLabel(), EBDetId::ieta(), ievt_, init_, EBDetId::iphi(), edm::HandleBase::isValid(), edm::ESHandleBase::isValid(), M_PI, meBCEne_, meBCEneMap_, meBCEneMapProjEta_, meBCEneMapProjPhi_, meBCETMap_, meBCETMapProjEta_, meBCETMapProjPhi_, meBCNum_, meBCNumMap_, meBCNumMapProjEta_, meBCNumMapProjPhi_, meBCSiz_, meBCSizMap_, meBCSizMapProjEta_, meBCSizMapProjPhi_, meInvMassHigh_, meInvMassHighSel_, meInvMassJPsi_, meInvMassJPsiSel_, meInvMassPi0_, meInvMassPi0Sel_, meInvMassZ0_, meInvMassZ0Sel_, mes1s9_, mes1s9thr_, mes9s25_, meSCCrystalSiz_, meSCEne2_, meSCEne_, meSCEneLowScale_, meSCEneVsEMax_, meSCMapSingleCrystal_, meSCNum_, meSCSeedEne_, meSCSeedMapOcc_, meSCSiz_, EcalDCCHeaderBlock::MTCC, EcalDCCHeaderBlock::PHYSICS_GLOBAL, EcalDCCHeaderBlock::PHYSICS_LOCAL, edm::ESHandle< T >::product(), edm::Handle< T >::product(), setup(), funct::sin(), Numbers::subDet(), SuperClusterCollection_, swap(), thrCandEt_, thrClusEt_, thrS4S9_, and EBDetId::zside().
{ bool enable = false; edm::Handle<EcalRawDataCollection> dcchs; if ( e.getByLabel(EcalRawDataCollection_, dcchs) ) { for ( EcalRawDataCollection::const_iterator dcchItr = dcchs->begin(); dcchItr != dcchs->end(); ++dcchItr ) { if ( Numbers::subDet( *dcchItr ) != EcalBarrel ) continue; if ( dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH4 || dcchItr->getRunType() == EcalDCCHeaderBlock::BEAMH2 ) enable = true; 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; break; } } else { enable = true; edm::LogWarning("EBClusterTask") << EcalRawDataCollection_ << " not available"; } if ( ! enable ) return; if ( ! init_ ) this->setup(); ievt_++; // ECAL topology edm::ESHandle<CaloTopology> pTopology; c.get<CaloTopologyRecord>().get(pTopology); if ( !pTopology.isValid() ) { edm::LogWarning("EBClusterTask") << "Topology not valid"; return; } const CaloTopology* topology = pTopology.product(); // recHits edm::Handle< EcalRecHitCollection > pEBRecHits; e.getByLabel( EcalRecHitCollection_, pEBRecHits ); if ( !pEBRecHits.isValid() ) { edm::LogWarning("EBClusterTask") << "RecHit collection " << EcalRecHitCollection_ << " not available."; return; } const EcalRecHitCollection* ebRecHits = pEBRecHits.product(); reco::BasicClusterCollection bcSel; // --- Barrel Basic Clusters --- edm::Handle<reco::BasicClusterCollection> pBasicClusters; if ( e.getByLabel(BasicClusterCollection_, pBasicClusters) ) { int nbcc = pBasicClusters->size(); if ( nbcc > 0 ) meBCNum_->Fill(float(nbcc)); for ( reco::BasicClusterCollection::const_iterator bCluster = pBasicClusters->begin(); bCluster != pBasicClusters->end(); ++bCluster ) { meBCEne_->Fill(bCluster->energy()); meBCSiz_->Fill(float(bCluster->size())); float xphi = bCluster->phi(); if ( xphi > M_PI*(9-1.5)/9 ) xphi = xphi - M_PI*2; meBCEneMap_->Fill(xphi, bCluster->eta(), bCluster->energy()); meBCEneMapProjEta_->Fill(bCluster->eta(), bCluster->energy()); meBCEneMapProjPhi_->Fill(xphi, bCluster->energy()); meBCNumMap_->Fill(xphi, bCluster->eta()); meBCNumMapProjEta_->Fill(bCluster->eta()); meBCNumMapProjPhi_->Fill(xphi); meBCSizMap_->Fill(xphi, bCluster->eta(), float(bCluster->size())); meBCSizMapProjEta_->Fill(bCluster->eta(), float(bCluster->size())); meBCSizMapProjPhi_->Fill(xphi, float(bCluster->size())); meBCETMap_->Fill(xphi, bCluster->eta(), float(bCluster->energy()) * sin(bCluster->position().theta())); meBCETMapProjEta_->Fill(bCluster->eta(), float(bCluster->energy()) * sin(bCluster->position().theta())); meBCETMapProjPhi_->Fill(xphi, float(bCluster->energy()) * sin(bCluster->position().theta())); float e2x2 = EcalClusterTools::e2x2( *bCluster, ebRecHits, topology ); float e3x3 = EcalClusterTools::e3x3( *bCluster, ebRecHits, topology ); // fill the selected cluster collection float pt = std::abs( bCluster->energy()*sin(bCluster->position().theta()) ); if ( pt > thrClusEt_ && e2x2/e3x3 > thrS4S9_ ) bcSel.push_back(*bCluster); } } else { edm::LogWarning("EBClusterTask") << BasicClusterCollection_ << " not available"; } for ( reco::BasicClusterCollection::const_iterator bc1 = bcSel.begin(); bc1 != bcSel.end(); ++bc1 ) { TLorentzVector bc1P; bc1P.SetPtEtaPhiE(std::abs(bc1->energy()*sin(bc1->position().theta())), bc1->eta(), bc1->phi(), bc1->energy()); for ( reco::BasicClusterCollection::const_iterator bc2 = bc1+1; bc2 != bcSel.end(); ++bc2 ) { TLorentzVector bc2P; bc2P.SetPtEtaPhiE(std::abs(bc2->energy()*sin(bc2->position().theta())), bc2->eta(), bc2->phi(), bc2->energy()); TLorentzVector candP = bc1P + bc2P; if ( candP.Pt() > thrCandEt_ ) { float mass = candP.M(); if ( mass < 0.500 ) { meInvMassPi0Sel_->Fill( mass ); } else if ( mass > 2.9 && mass < 3.3 ) { meInvMassJPsiSel_->Fill( mass ); } else if ( mass > 40 && mass < 110 ) { meInvMassZ0Sel_->Fill( mass ); } else if ( mass > 110 ) { meInvMassHighSel_->Fill( mass ); } } } } // --- Barrel Super Clusters --- edm::Handle<reco::SuperClusterCollection> pSuperClusters; if ( e.getByLabel(SuperClusterCollection_, pSuperClusters) ) { int nscc = pSuperClusters->size(); if ( nscc > 0 ) meSCNum_->Fill(float(nscc)); TLorentzVector sc1_p(0,0,0,0); TLorentzVector sc2_p(0,0,0,0); for ( reco::SuperClusterCollection::const_iterator sCluster = pSuperClusters->begin(); sCluster != pSuperClusters->end(); ++sCluster ) { // energy, size meSCEne_->Fill( sCluster->energy() ); meSCSiz_->Fill( float(sCluster->clustersSize()) ); reco::CaloClusterPtr theSeed = sCluster->seed(); // Find the seed rec hit std::vector< std::pair<DetId,float> > sIds = sCluster->hitsAndFractions(); float eMax, e2nd; EcalRecHitCollection::const_iterator seedItr = ebRecHits->begin(); EcalRecHitCollection::const_iterator secondItr = ebRecHits->begin(); for(std::vector< std::pair<DetId,float> >::const_iterator idItr = sIds.begin(); idItr != sIds.end(); ++idItr) { DetId id = idItr->first; if(id.det() != DetId::Ecal) { continue; } EcalRecHitCollection::const_iterator hitItr = ebRecHits->find(id); if(hitItr == ebRecHits->end()) { continue; } if(hitItr->energy() > secondItr->energy()) { secondItr = hitItr; } if(hitItr->energy() > seedItr->energy()) { std::swap(seedItr,secondItr); } } eMax = seedItr->energy(); e2nd = secondItr->energy(); EBDetId seedId = (EBDetId) seedItr->id(); float e3x3 = EcalClusterTools::e3x3( *theSeed, ebRecHits, topology ); float e5x5 = EcalClusterTools::e5x5( *theSeed, ebRecHits, topology ); meSCCrystalSiz_->Fill(sIds.size()); meSCSeedEne_->Fill(eMax); meSCEne2_->Fill(eMax+e2nd); meSCEneVsEMax_->Fill(eMax,sCluster->energy()); meSCEneLowScale_->Fill(sCluster->energy()); // Prepare to fill maps int ebeta = seedId.ieta(); int ebphi = seedId.iphi(); float xebeta = ebeta - 0.5 * seedId.zside(); float xebphi = ebphi - 0.5; meSCSeedMapOcc_->Fill(xebphi,xebeta); if(sIds.size() == 1) meSCMapSingleCrystal_->Fill(xebphi,xebeta); mes1s9_->Fill( eMax/e3x3 ); if ( eMax > 3.0 ) mes1s9thr_->Fill( eMax/e3x3 ); mes9s25_->Fill( e3x3/e5x5 ); if ( nscc >= 2 ) { // look for the two most energetic super clusters if ( sCluster->energy() > sc1_p.Energy() ) { sc2_p=sc1_p; sc1_p.SetPtEtaPhiE(std::abs(sCluster->energy()*sin(sCluster->position().theta())), sCluster->eta(), sCluster->phi(), sCluster->energy()); } else if ( sCluster->energy() > sc2_p.Energy() ) { sc2_p.SetPtEtaPhiE(std::abs(sCluster->energy()*sin(sCluster->position().theta())), sCluster->eta(), sCluster->phi(), sCluster->energy()); } } } // Get the invariant mass of the two most energetic super clusters if ( nscc >= 2 ) { TLorentzVector sum = sc1_p+sc2_p; float mass = sum.M(); if ( mass < 0.500 ) { meInvMassPi0_->Fill( mass ); } else if ( mass > 2.9 && mass < 3.3 ) { meInvMassJPsi_->Fill( mass ); } else if ( mass > 40 && mass < 110 ) { meInvMassZ0_->Fill( mass ); } else if ( mass > 110 ) { meInvMassHigh_->Fill( mass ); } } } else { edm::LogWarning("EBClusterTask") << SuperClusterCollection_ << " not available"; } }
void EBClusterTask::beginJob | ( | void | ) | [protected, virtual] |
BeginJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 118 of file EBClusterTask.cc.
References dqmStore_, ievt_, prefixME_, DQMStore::rmdir(), and DQMStore::setCurrentFolder().
void EBClusterTask::beginRun | ( | const edm::Run & | r, |
const edm::EventSetup & | c | ||
) | [protected, virtual] |
BeginRun.
Reimplemented from edm::EDAnalyzer.
Definition at line 129 of file EBClusterTask.cc.
References Numbers::initGeometry(), mergeRuns_, and reset().
{ Numbers::initGeometry(c, false); if ( ! mergeRuns_ ) this->reset(); }
void EBClusterTask::cleanup | ( | void | ) | [protected] |
Cleanup.
Definition at line 389 of file EBClusterTask.cc.
References dqmStore_, MonitorElement::getName(), init_, meBCEne_, meBCEneMap_, meBCEneMapProjEta_, meBCEneMapProjPhi_, meBCETMap_, meBCETMapProjEta_, meBCETMapProjPhi_, meBCNum_, meBCNumMap_, meBCNumMapProjEta_, meBCNumMapProjPhi_, meBCSiz_, meBCSizMap_, meBCSizMapProjEta_, meBCSizMapProjPhi_, meInvMassHigh_, meInvMassHighSel_, meInvMassJPsi_, meInvMassJPsiSel_, meInvMassPi0_, meInvMassPi0Sel_, meInvMassZ0_, meInvMassZ0Sel_, mes1s9_, mes1s9thr_, mes9s25_, meSCCrystalSiz_, meSCEne2_, meSCEne_, meSCEneLowScale_, meSCEneVsEMax_, meSCMapSingleCrystal_, meSCNum_, meSCSeedEne_, meSCSeedMapOcc_, meSCSiz_, prefixME_, DQMStore::removeElement(), and DQMStore::setCurrentFolder().
Referenced by endJob().
{ if ( ! init_ ) return; if ( dqmStore_ ) { dqmStore_->setCurrentFolder(prefixME_ + "/EBClusterTask"); if ( meBCEne_ ) dqmStore_->removeElement( meBCEne_->getName() ); meBCEne_ = 0; if ( meBCNum_ ) dqmStore_->removeElement( meBCNum_->getName() ); meBCNum_ = 0; if ( meBCSiz_ ) dqmStore_->removeElement( meBCSiz_->getName() ); meBCSiz_ = 0; if ( meBCEneMap_ ) dqmStore_->removeElement( meBCEneMap_->getName() ); meBCEneMap_ = 0; if ( meBCNumMap_ ) dqmStore_->removeElement( meBCNumMap_->getName() ); meBCNumMap_ = 0; if ( meBCETMap_ ) dqmStore_->removeElement( meBCETMap_->getName() ); meBCETMap_ = 0; if ( meBCSizMap_ ) dqmStore_->removeElement( meBCSizMap_->getName() ); meBCSizMap_ = 0; if ( meBCEneMapProjEta_ ) dqmStore_->removeElement( meBCEneMapProjEta_->getName() ); meBCEneMapProjEta_ = 0; if ( meBCEneMapProjPhi_ ) dqmStore_->removeElement( meBCEneMapProjPhi_->getName() ); meBCEneMapProjPhi_ = 0; if ( meBCNumMapProjEta_ ) dqmStore_->removeElement( meBCNumMapProjEta_->getName() ); meBCNumMapProjEta_ = 0; if ( meBCNumMapProjPhi_ ) dqmStore_->removeElement( meBCNumMapProjPhi_->getName() ); meBCNumMapProjPhi_ = 0; if ( meBCETMapProjEta_ ) dqmStore_->removeElement( meBCETMapProjEta_->getName() ); meBCETMapProjEta_ = 0; if ( meBCETMapProjPhi_ ) dqmStore_->removeElement( meBCETMapProjPhi_->getName() ); meBCETMapProjPhi_ = 0; if ( meBCSizMapProjEta_ ) dqmStore_->removeElement( meBCSizMapProjEta_->getName() ); meBCSizMapProjEta_ = 0; if ( meBCSizMapProjPhi_ ) dqmStore_->removeElement( meBCSizMapProjPhi_->getName() ); meBCSizMapProjPhi_ = 0; if ( meSCEne_ ) dqmStore_->removeElement( meSCEne_->getName() ); meSCEne_ = 0; if ( meSCNum_ ) dqmStore_->removeElement( meSCNum_->getName() ); meSCNum_ = 0; if ( meSCSiz_ ) dqmStore_->removeElement( meSCSiz_->getName() ); meSCSiz_ = 0; if ( meSCCrystalSiz_ ) dqmStore_->removeElement( meSCCrystalSiz_->getName() ); meSCCrystalSiz_ = 0; if ( meSCSeedEne_ ) dqmStore_->removeElement( meSCSeedEne_->getName() ); meSCSeedEne_ = 0; if ( meSCEne2_ ) dqmStore_->removeElement( meSCEne2_->getName() ); meSCEne2_ = 0; if ( meSCEneVsEMax_ ) dqmStore_->removeElement( meSCEneVsEMax_->getName() ); meSCEneVsEMax_ = 0; if ( meSCEneLowScale_ ) dqmStore_->removeElement( meSCEneLowScale_->getName() ); meSCEneLowScale_ = 0; if ( meSCSeedMapOcc_ ) dqmStore_->removeElement( meSCSeedMapOcc_->getName() ); meSCSeedMapOcc_ = 0; if ( meSCMapSingleCrystal_ ) dqmStore_->removeElement( meSCMapSingleCrystal_->getName() ); meSCMapSingleCrystal_ = 0; if ( mes1s9_ ) dqmStore_->removeElement( mes1s9_->getName() ); mes1s9_ = 0; if ( mes1s9thr_ ) dqmStore_->removeElement( mes1s9thr_->getName() ); mes1s9thr_ = 0; if ( mes9s25_ ) dqmStore_->removeElement( mes9s25_->getName() ); mes9s25_ = 0; if ( meInvMassPi0_ ) dqmStore_->removeElement( meInvMassPi0_->getName() ); meInvMassPi0_ = 0; if ( meInvMassJPsi_ ) dqmStore_->removeElement( meInvMassJPsi_->getName() ); meInvMassJPsi_ = 0; if ( meInvMassZ0_ ) dqmStore_->removeElement( meInvMassZ0_->getName() ); meInvMassZ0_ = 0; if ( meInvMassHigh_ ) dqmStore_->removeElement( meInvMassHigh_->getName() ); meInvMassHigh_ = 0; if ( meInvMassPi0Sel_ ) dqmStore_->removeElement( meInvMassPi0Sel_->getName() ); meInvMassPi0Sel_ = 0; if ( meInvMassJPsiSel_ ) dqmStore_->removeElement( meInvMassJPsiSel_->getName() ); meInvMassJPsiSel_ = 0; if ( meInvMassZ0Sel_ ) dqmStore_->removeElement( meInvMassZ0Sel_->getName() ); meInvMassZ0Sel_ = 0; if ( meInvMassHighSel_ ) dqmStore_->removeElement( meInvMassHighSel_->getName() ); meInvMassHighSel_ = 0; } init_ = false; }
void EBClusterTask::endJob | ( | void | ) | [protected, virtual] |
EndJob.
Reimplemented from edm::EDAnalyzer.
Definition at line 510 of file EBClusterTask.cc.
References cleanup(), enableCleanup_, and ievt_.
{ edm::LogInfo("EBClusterTask") << "analyzed " << ievt_ << " events"; if ( enableCleanup_ ) this->cleanup(); }
void EBClusterTask::endRun | ( | const edm::Run & | r, |
const edm::EventSetup & | c | ||
) | [protected, virtual] |
void EBClusterTask::reset | ( | void | ) | [protected] |
Reset.
Definition at line 141 of file EBClusterTask.cc.
References meBCEne_, meBCEneMap_, meBCEneMapProjEta_, meBCEneMapProjPhi_, meBCETMap_, meBCETMapProjEta_, meBCETMapProjPhi_, meBCNum_, meBCNumMap_, meBCNumMapProjEta_, meBCNumMapProjPhi_, meBCSiz_, meBCSizMap_, meBCSizMapProjEta_, meBCSizMapProjPhi_, meInvMassHigh_, meInvMassHighSel_, meInvMassJPsi_, meInvMassJPsiSel_, meInvMassPi0_, meInvMassPi0Sel_, meInvMassZ0_, meInvMassZ0Sel_, mes1s9_, mes1s9thr_, mes9s25_, meSCCrystalSiz_, meSCEne2_, meSCEne_, meSCEneLowScale_, meSCEneVsEMax_, meSCMapSingleCrystal_, meSCNum_, meSCSeedEne_, meSCSeedMapOcc_, meSCSiz_, and MonitorElement::Reset().
Referenced by beginRun().
{ if ( meBCEne_ ) meBCEne_->Reset(); if ( meBCNum_ ) meBCNum_->Reset(); if ( meBCSiz_ ) meBCSiz_->Reset(); if ( meBCEneMap_ ) meBCEneMap_->Reset(); if ( meBCNumMap_ ) meBCNumMap_->Reset(); if ( meBCETMap_ ) meBCETMap_->Reset(); if ( meBCSizMap_ ) meBCSizMap_->Reset(); if ( meBCEneMapProjEta_ ) meBCEneMapProjEta_->Reset(); if ( meBCEneMapProjPhi_ ) meBCEneMapProjPhi_->Reset(); if ( meBCNumMapProjEta_ ) meBCNumMapProjEta_->Reset(); if ( meBCNumMapProjPhi_ ) meBCNumMapProjPhi_->Reset(); if ( meBCETMapProjEta_ ) meBCETMapProjEta_->Reset(); if ( meBCETMapProjPhi_ ) meBCETMapProjPhi_->Reset(); if ( meBCSizMapProjEta_ ) meBCSizMapProjEta_->Reset(); if ( meBCSizMapProjPhi_ ) meBCSizMapProjPhi_->Reset(); if ( meSCEne_ ) meSCEne_->Reset(); if ( meSCNum_ ) meSCNum_->Reset(); if ( meSCSiz_ ) meSCSiz_->Reset(); if ( meSCCrystalSiz_ ) meSCCrystalSiz_->Reset(); if ( meSCSeedEne_ ) meSCSeedEne_->Reset(); if ( meSCEne2_ ) meSCEne2_->Reset(); if ( meSCEneVsEMax_ ) meSCEneVsEMax_->Reset(); if ( meSCEneLowScale_ ) meSCEneLowScale_->Reset(); if ( meSCSeedMapOcc_ ) meSCSeedMapOcc_->Reset(); if ( meSCMapSingleCrystal_ ) meSCMapSingleCrystal_->Reset(); if ( mes1s9_ ) mes1s9_->Reset(); if ( mes1s9thr_ ) mes1s9thr_->Reset(); if ( mes9s25_ ) mes9s25_->Reset(); if ( meInvMassPi0_ ) meInvMassPi0_->Reset(); if ( meInvMassJPsi_ ) meInvMassJPsi_->Reset(); if ( meInvMassZ0_ ) meInvMassZ0_->Reset(); if ( meInvMassHigh_ ) meInvMassHigh_->Reset(); if ( meInvMassPi0Sel_ ) meInvMassPi0Sel_->Reset(); if ( meInvMassJPsiSel_ ) meInvMassJPsiSel_->Reset(); if ( meInvMassZ0Sel_ ) meInvMassZ0Sel_->Reset(); if ( meInvMassHighSel_ ) meInvMassHighSel_->Reset(); }
void EBClusterTask::setup | ( | void | ) | [protected] |
Setup.
Definition at line 217 of file EBClusterTask.cc.
References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), DQMStore::bookProfile2D(), dqmStore_, trackerHits::histo, init_, M_PI, meBCEne_, meBCEneMap_, meBCEneMapProjEta_, meBCEneMapProjPhi_, meBCETMap_, meBCETMapProjEta_, meBCETMapProjPhi_, meBCNum_, meBCNumMap_, meBCNumMapProjEta_, meBCNumMapProjPhi_, meBCSiz_, meBCSizMap_, meBCSizMapProjEta_, meBCSizMapProjPhi_, meInvMassHigh_, meInvMassHighSel_, meInvMassJPsi_, meInvMassJPsiSel_, meInvMassPi0_, meInvMassPi0Sel_, meInvMassZ0_, meInvMassZ0Sel_, mes1s9_, mes1s9thr_, mes9s25_, meSCCrystalSiz_, meSCEne2_, meSCEne_, meSCEneLowScale_, meSCEneVsEMax_, meSCMapSingleCrystal_, meSCNum_, meSCSeedEne_, meSCSeedMapOcc_, meSCSiz_, prefixME_, MonitorElement::setAxisTitle(), and DQMStore::setCurrentFolder().
Referenced by analyze().
{ init_ = true; char histo[200]; if ( dqmStore_ ) { dqmStore_->setCurrentFolder(prefixME_ + "/EBClusterTask"); sprintf(histo, "EBCLT BC energy"); meBCEne_ = dqmStore_->book1D(histo, histo, 100, 0., 150.); meBCEne_->setAxisTitle("energy (GeV)", 1); sprintf(histo, "EBCLT BC number"); meBCNum_ = dqmStore_->book1D(histo, histo, 100, 0., 100.); meBCNum_->setAxisTitle("number of clusters", 1); sprintf(histo, "EBCLT BC size"); meBCSiz_ = dqmStore_->book1D(histo, histo, 100, 0., 100.); meBCSiz_->setAxisTitle("cluster size", 1); sprintf(histo, "EBCLT BC energy map"); meBCEneMap_ = dqmStore_->bookProfile2D(histo, histo, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 34, -1.479, 1.479, 100, 0., 500., "s"); meBCEneMap_->setAxisTitle("phi", 1); meBCEneMap_->setAxisTitle("eta", 2); sprintf(histo, "EBCLT BC number map"); meBCNumMap_ = dqmStore_->book2D(histo, histo, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 34, -1.479, 1.479); meBCNumMap_->setAxisTitle("phi", 1); meBCNumMap_->setAxisTitle("eta", 2); sprintf(histo, "EBCLT BC ET map"); meBCETMap_ = dqmStore_->bookProfile2D(histo, histo, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 34, -1.479, 1.479, 100, 0., 500., "s"); meBCETMap_->setAxisTitle("phi", 1); meBCETMap_->setAxisTitle("eta", 2); sprintf(histo, "EBCLT BC size map"); meBCSizMap_ = dqmStore_->bookProfile2D(histo, histo, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 34, -1.479, 1.479, 100, 0., 100., "s"); meBCSizMap_->setAxisTitle("phi", 1); meBCSizMap_->setAxisTitle("eta", 2); sprintf(histo, "EBCLT BC energy projection eta"); meBCEneMapProjEta_ = dqmStore_->bookProfile(histo, histo, 34, -1.479, 1.479, 100, 0., 500., "s"); meBCEneMapProjEta_->setAxisTitle("eta", 1); meBCEneMapProjEta_->setAxisTitle("energy (GeV)", 2); sprintf(histo, "EBCLT BC energy projection phi"); meBCEneMapProjPhi_ = dqmStore_->bookProfile(histo, histo, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 100, 0., 500., "s"); meBCEneMapProjPhi_->setAxisTitle("phi", 1); meBCEneMapProjPhi_->setAxisTitle("energy (GeV)", 2); sprintf(histo, "EBCLT BC number projection eta"); meBCNumMapProjEta_ = dqmStore_->book1D(histo, histo, 34, -1.479, 1.479); meBCNumMapProjEta_->setAxisTitle("eta", 1); meBCNumMapProjEta_->setAxisTitle("number of clusters", 2); sprintf(histo, "EBCLT BC number projection phi"); meBCNumMapProjPhi_ = dqmStore_->book1D(histo, histo, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9); meBCNumMapProjPhi_->setAxisTitle("phi", 1); meBCNumMapProjPhi_->setAxisTitle("number of clusters", 2); sprintf(histo, "EBCLT BC ET projection eta"); meBCETMapProjEta_ = dqmStore_->bookProfile(histo, histo, 34, -1.479, 1.479, 100, 0., 500., "s"); meBCETMapProjEta_->setAxisTitle("eta", 1); meBCETMapProjEta_->setAxisTitle("transverse energy (GeV)", 2); sprintf(histo, "EBCLT BC ET projection phi"); meBCETMapProjPhi_ = dqmStore_->bookProfile(histo, histo, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 100, 0., 500., "s"); meBCETMapProjPhi_->setAxisTitle("phi", 1); meBCETMapProjPhi_->setAxisTitle("transverse energy (GeV)", 2); sprintf(histo, "EBCLT BC size projection eta"); meBCSizMapProjEta_ = dqmStore_->bookProfile(histo, histo, 34, -1.479, 1.479, 100, 0., 100., "s"); meBCSizMapProjEta_->setAxisTitle("eta", 1); meBCSizMapProjEta_->setAxisTitle("cluster size", 2); sprintf(histo, "EBCLT BC size projection phi"); meBCSizMapProjPhi_ = dqmStore_->bookProfile(histo, histo, 72, -M_PI*(9+1.5)/9, M_PI*(9-1.5)/9, 100, 0., 100., "s"); meBCSizMapProjPhi_->setAxisTitle("phi", 1); meBCSizMapProjPhi_->setAxisTitle("cluster size", 2); sprintf(histo, "EBCLT SC energy"); meSCEne_ = dqmStore_->book1D(histo, histo, 100, 0., 150.); meSCEne_->setAxisTitle("energy (GeV)", 1); sprintf(histo, "EBCLT SC number"); meSCNum_ = dqmStore_->book1D(histo, histo, 50, 0., 50.); meSCNum_->setAxisTitle("number of clusters", 1); sprintf(histo, "EBCLT SC size"); meSCSiz_ = dqmStore_->book1D(histo, histo, 50, 0., 50.); meSCSiz_->setAxisTitle("cluster size", 1); sprintf(histo, "EBCLT SC size (crystal)"); meSCCrystalSiz_ = dqmStore_->book1D(histo, histo, 150, 0, 150); meSCCrystalSiz_->setAxisTitle("cluster size in crystals", 1); sprintf(histo, "EBCLT SC seed crystal energy"); meSCSeedEne_ = dqmStore_->book1D(histo, histo, 100, 0., 10.); meSCSeedEne_->setAxisTitle("seed crystal energy (GeV)", 1); sprintf(histo, "EBCLT SC e2"); meSCEne2_ = dqmStore_->book1D(histo, histo, 100, 0., 10.); meSCEne2_->setAxisTitle("seed + highest neighbor crystal energy (GeV)", 1); sprintf(histo, "EBCLT SC energy vs seed crystal energy"); meSCEneVsEMax_ = dqmStore_->book2D(histo, histo, 50, 0., 10., 50, 0., 10.); meSCEneVsEMax_->setAxisTitle("seed crystal energy (GeV)", 1); meSCEneVsEMax_->setAxisTitle("cluster energy (GeV)", 2); sprintf(histo, "EBCLT SC energy (low scale)"); meSCEneLowScale_ = dqmStore_->book1D(histo, histo, 200, 0., 10.); meSCEneLowScale_->setAxisTitle("cluster energy (GeV)", 1); sprintf(histo, "EBCLT SC seed occupancy map"); meSCSeedMapOcc_ = dqmStore_->book2D(histo, histo, 72, 0., 360., 34, -85, 85); meSCSeedMapOcc_->setAxisTitle("jphi", 1); meSCSeedMapOcc_->setAxisTitle("jeta", 2); sprintf(histo, "EBCLT SC single crystal cluster seed occupancy map"); meSCMapSingleCrystal_ = dqmStore_->book2D(histo, histo, 72, 0., 360., 34, -85, 85); meSCMapSingleCrystal_->setAxisTitle("jphi", 1); meSCMapSingleCrystal_->setAxisTitle("jeta", 2); sprintf(histo, "EBCLT s1s9"); mes1s9_ = dqmStore_->book1D(histo, histo, 50, 0., 1.5); mes1s9_->setAxisTitle("s1/s9", 1); sprintf(histo, "EBCLT s1s9 thr"); mes1s9thr_ = dqmStore_->book1D(histo, histo, 50, 0., 1.5); mes1s9thr_->setAxisTitle("s1/s9", 1); sprintf(histo, "EBCLT s9s25"); mes9s25_ = dqmStore_->book1D(histo, histo, 75, 0., 1.5); mes9s25_->setAxisTitle("s9/s25", 1); sprintf(histo, "EBCLT dicluster invariant mass Pi0"); meInvMassPi0_ = dqmStore_->book1D(histo, histo, 50, 0.0, 0.500); meInvMassPi0_->setAxisTitle("mass (GeV)", 1); sprintf(histo, "EBCLT dicluster invariant mass JPsi"); meInvMassJPsi_ = dqmStore_->book1D(histo, histo, 50, 2.9, 3.3); meInvMassJPsi_->setAxisTitle("mass (GeV)", 1); sprintf(histo, "EBCLT dicluster invariant mass Z0"); meInvMassZ0_ = dqmStore_->book1D(histo, histo, 50, 40, 110); meInvMassZ0_->setAxisTitle("mass (GeV)", 1); sprintf(histo, "EBCLT dicluster invariant mass high"); meInvMassHigh_ = dqmStore_->book1D(histo, histo, 500, 110, 3000); meInvMassHigh_->setAxisTitle("mass (GeV)", 1); sprintf(histo, "EBCLT dicluster invariant mass Pi0 sel"); meInvMassPi0Sel_ = dqmStore_->book1D(histo, histo, 50, 0.00, 0.500); meInvMassPi0Sel_->setAxisTitle("mass (GeV)", 1); sprintf(histo, "EBCLT dicluster invariant mass JPsi sel"); meInvMassJPsiSel_ = dqmStore_->book1D(histo, histo, 50, 2.9, 3.3); meInvMassJPsiSel_->setAxisTitle("mass (GeV)", 1); sprintf(histo, "EBCLT dicluster invariant mass Z0 sel"); meInvMassZ0Sel_ = dqmStore_->book1D(histo, histo, 50, 40, 110); meInvMassZ0Sel_->setAxisTitle("mass (GeV)", 1); sprintf(histo, "EBCLT dicluster invariant mass high sel"); meInvMassHighSel_ = dqmStore_->book1D(histo, histo, 500, 110, 3000); meInvMassHighSel_->setAxisTitle("mass (GeV)", 1); } }
Definition at line 69 of file EBClusterTask.h.
Referenced by analyze(), and EBClusterTask().
DQMStore* EBClusterTask::dqmStore_ [private] |
Definition at line 60 of file EBClusterTask.h.
Referenced by beginJob(), cleanup(), EBClusterTask(), and setup().
Definition at line 68 of file EBClusterTask.h.
Referenced by analyze(), and EBClusterTask().
Definition at line 71 of file EBClusterTask.h.
Referenced by analyze(), and EBClusterTask().
bool EBClusterTask::enableCleanup_ [private] |
Definition at line 64 of file EBClusterTask.h.
Referenced by EBClusterTask(), and endJob().
int EBClusterTask::ievt_ [private] |
Definition at line 58 of file EBClusterTask.h.
Referenced by analyze(), beginJob(), and endJob().
bool EBClusterTask::init_ [private] |
Definition at line 106 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), and setup().
MonitorElement* EBClusterTask::meBCEne_ [private] |
Definition at line 73 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meBCEneMap_ [private] |
Definition at line 77 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement * EBClusterTask::meBCEneMapProjEta_ [private] |
Definition at line 77 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement * EBClusterTask::meBCEneMapProjPhi_ [private] |
Definition at line 77 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meBCETMap_ [private] |
Definition at line 79 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement * EBClusterTask::meBCETMapProjEta_ [private] |
Definition at line 79 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement * EBClusterTask::meBCETMapProjPhi_ [private] |
Definition at line 79 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meBCNum_ [private] |
Definition at line 74 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meBCNumMap_ [private] |
Definition at line 78 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement * EBClusterTask::meBCNumMapProjEta_ [private] |
Definition at line 78 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement * EBClusterTask::meBCNumMapProjPhi_ [private] |
Definition at line 78 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meBCSiz_ [private] |
Definition at line 75 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meBCSizMap_ [private] |
Definition at line 80 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement * EBClusterTask::meBCSizMapProjEta_ [private] |
Definition at line 80 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement * EBClusterTask::meBCSizMapProjPhi_ [private] |
Definition at line 80 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meInvMassHigh_ [private] |
Definition at line 100 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meInvMassHighSel_ [private] |
Definition at line 104 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meInvMassJPsi_ [private] |
Definition at line 98 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meInvMassJPsiSel_ [private] |
Definition at line 102 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meInvMassPi0_ [private] |
Definition at line 97 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meInvMassPi0Sel_ [private] |
Definition at line 101 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meInvMassZ0_ [private] |
Definition at line 99 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meInvMassZ0Sel_ [private] |
Definition at line 103 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
bool EBClusterTask::mergeRuns_ [private] |
Definition at line 66 of file EBClusterTask.h.
Referenced by beginRun(), and EBClusterTask().
MonitorElement* EBClusterTask::mes1s9_ [private] |
Definition at line 94 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::mes1s9thr_ [private] |
Definition at line 95 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::mes9s25_ [private] |
Definition at line 96 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meSCCrystalSiz_ [private] |
Definition at line 86 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meSCEne2_ [private] |
Definition at line 88 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meSCEne_ [private] |
Definition at line 82 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meSCEneLowScale_ [private] |
Definition at line 90 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meSCEneVsEMax_ [private] |
Definition at line 89 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
Definition at line 92 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meSCNum_ [private] |
Definition at line 83 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meSCSeedEne_ [private] |
Definition at line 87 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meSCSeedMapOcc_ [private] |
Definition at line 91 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
MonitorElement* EBClusterTask::meSCSiz_ [private] |
Definition at line 84 of file EBClusterTask.h.
Referenced by analyze(), cleanup(), EBClusterTask(), reset(), and setup().
std::string EBClusterTask::prefixME_ [private] |
Definition at line 62 of file EBClusterTask.h.
Referenced by beginJob(), cleanup(), EBClusterTask(), and setup().
Definition at line 70 of file EBClusterTask.h.
Referenced by analyze(), and EBClusterTask().
float EBClusterTask::thrCandEt_ [private] |
Definition at line 108 of file EBClusterTask.h.
Referenced by analyze(), and EBClusterTask().
float EBClusterTask::thrClusEt_ [private] |
Definition at line 108 of file EBClusterTask.h.
Referenced by analyze(), and EBClusterTask().
float EBClusterTask::thrS4S9_ [private] |
Definition at line 108 of file EBClusterTask.h.
Referenced by analyze(), and EBClusterTask().