#include <ClusterTask.h>

Inheritance diagram for ecaldqm::ClusterTask:

Public Types
enum	TriggerTypes { kEcalTrigger, kHcalTrigger, kCSCTrigger, kDTTrigger, kRPCTrigger, nTriggerTypes }

Public Types inherited from ecaldqm::DQWorkerTask
typedef EcalDCCHeaderBlock::EcalDCCEventSettings	EventSettings

Public Member Functions
void	addDependencies (DependencySet &) override

bool	analyze (void const *, Collections) override

void	beginEvent (edm::Event const &, edm::EventSetup const &) override

	ClusterTask ()

void	endEvent (edm::Event const &, edm::EventSetup const &) override

bool	filterRunType (short const *) override

void	runOnBasicClusters (edm::View< reco::CaloCluster > const &, Collections)

void	runOnRecHits (EcalRecHitCollection const &, Collections)

void	runOnSuperClusters (reco::SuperClusterCollection const &, Collections)

void	setTokens (edm::ConsumesCollector &) override

	~ClusterTask () override

Public Member Functions inherited from ecaldqm::DQWorkerTask
	DQWorkerTask ()

virtual bool	filterTrigger (edm::TriggerResultsByName const &)

void	recoverStats ()

void	softReset ()

	~DQWorkerTask () override

Public Member Functions inherited from ecaldqm::DQWorker
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	beginRun (edm::Run const &, edm::EventSetup const &)

virtual void	bookMEs (DQMStore::IBooker &)

	DQWorker ()

virtual void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	endRun (edm::Run const &, edm::EventSetup const &)

std::string const &	getName () const

bool	onlineMode () const

virtual void	releaseMEs ()

void	setEventNumber (edm::EventNumber_t _e)

void	setLumiNumber (edm::LuminosityBlockNumber_t _l)

void	setRunNumber (edm::RunNumber_t _r)

void	setTime (time_t _t)

virtual	~DQWorker ()(false)

Private Member Functions
void	setParams (edm::ParameterSet const &) override

Private Attributes
bool	doExtra_

EcalRecHitCollection const *	ebHits_

EcalRecHitCollection const *	eeHits_

std::vector< std::string >	egTriggerAlgos_

float	energyThreshold_

edm::InputTag	L1GlobalTriggerReadoutRecordTag_

edm::EDGetTokenT< L1GlobalTriggerReadoutRecord >	L1GlobalTriggerReadoutRecordToken_

edm::InputTag	L1MuGMTReadoutCollectionTag_

edm::EDGetTokenT< L1MuGMTReadoutCollection >	L1MuGMTReadoutCollectionToken_

float	swissCrossMaxThreshold_

std::bitset< nTriggerTypes >	triggered_

unsigned	trigTypeToME_ [nTriggerTypes]

Additional Inherited Members
Static Public Member Functions inherited from ecaldqm::DQWorkerTask
static void	fillDescriptions (edm::ParameterSetDescription &)

Static Public Member Functions inherited from ecaldqm::DQWorker
static void	fillDescriptions (edm::ParameterSetDescription &_desc)

Protected Member Functions inherited from ecaldqm::DQWorkerTask
void	setME (edm::ParameterSet const &) final

Protected Member Functions inherited from ecaldqm::DQWorker
void	initialize (std::string const &_name, edm::ParameterSet const &)

void	print_ (std::string const &, int=0) const

virtual void	setSource (edm::ParameterSet const &)

void	setVerbosity (int _verbosity)

Protected Attributes inherited from ecaldqm::DQWorkerTask
std::set< std::string >	resettable_

Protected Attributes inherited from ecaldqm::DQWorker
bool	booked_

MESetCollection	MEs_

std::string	name_

bool	onlineMode_

Timestamp	timestamp_

int	verbosity_

bool	willConvertToEDM_

Detailed Description

Definition at line 21 of file ClusterTask.h.

Member Enumeration Documentation

enum ecaldqm::ClusterTask::TriggerTypes

Enumerator
kEcalTrigger
kHcalTrigger
kCSCTrigger
kDTTrigger
kRPCTrigger
nTriggerTypes

Definition at line 41 of file ClusterTask.h.

                       {
       kEcalTrigger,
       kHcalTrigger,
       kCSCTrigger,
       kDTTrigger,
       kRPCTrigger,
       nTriggerTypes
     };

Constructor & Destructor Documentation

ecaldqm::ClusterTask::ClusterTask ( )

Definition at line 22 of file ClusterTask.cc.

References L1GlobalTriggerReadoutRecordTag_, L1GlobalTriggerReadoutRecordToken_, L1MuGMTReadoutCollectionTag_, and L1MuGMTReadoutCollectionToken_.

                            :
     DQWorkerTask(),
     ebHits_(nullptr),
     eeHits_(nullptr),
     //    ievt_(0),
     //    massCalcPrescale_(_workerParams.getUntrackedParameter<int>("massCalcPrescale")),
     doExtra_(true),
     energyThreshold_(0.),
     swissCrossMaxThreshold_(3.),
     egTriggerAlgos_(),
     trigTypeToME_{0, 1, 2, 3, 4},
     L1GlobalTriggerReadoutRecordTag_(),
     L1MuGMTReadoutCollectionTag_(),
     L1GlobalTriggerReadoutRecordToken_(),
     L1MuGMTReadoutCollectionToken_()
   {
   }

ecaldqm::ClusterTask::~ClusterTask ( )

inlineoverride

Definition at line 24 of file ClusterTask.h.

References addDependencies(), analyze(), beginEvent(), endEvent(), filterRunType(), runOnBasicClusters(), runOnRecHits(), runOnSuperClusters(), and setTokens().

24 {}

Member Function Documentation

void ecaldqm::ClusterTask::addDependencies ( DependencySet & _dependencies )

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 84 of file ClusterTask.cc.

References ecaldqm::kEBRecHit, ecaldqm::kEBSuperCluster, ecaldqm::kEERecHit, ecaldqm::kEESuperCluster, and ecaldqm::DependencySet::push_back().

Referenced by ~ClusterTask().

   {
     _dependencies.push_back(Dependency(kEBSuperCluster, kEBRecHit));
     _dependencies.push_back(Dependency(kEESuperCluster, kEERecHit));
   }

bool ecaldqm::ClusterTask::analyze	(	void const *	_p,
		Collections	_collection
	)

inlineoverridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 70 of file ClusterTask.h.

References ecaldqm::kEBBasicCluster, ecaldqm::kEBRecHit, ecaldqm::kEBSuperCluster, ecaldqm::kEEBasicCluster, ecaldqm::kEERecHit, ecaldqm::kEESuperCluster, runOnBasicClusters(), runOnRecHits(), and runOnSuperClusters().

Referenced by ~ClusterTask().

                                                                          {
     switch(_collection){
     case kEBRecHit:
     case kEERecHit:
       if(_p) runOnRecHits(*static_cast<EcalRecHitCollection const*>(_p), _collection);
       return true;
       break;
     case kEBBasicCluster:
     case kEEBasicCluster:
       if(_p) runOnBasicClusters(*static_cast<edm::View<reco::CaloCluster> const*>(_p), _collection);
       return true;
       break;
     case kEBSuperCluster:
     case kEESuperCluster:
       if(_p) runOnSuperClusters(*static_cast<reco::SuperClusterCollection const*>(_p), _collection);
       return true;
       break;
     default:
       break;
     }
 
     return false;
   }

void ecaldqm::ClusterTask::beginEvent	(	edm::Event const &	_evt,
		edm::EventSetup const &	_es
	)

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 91 of file ClusterTask.cc.

References L1GtPsbWord::aData(), L1GtPsbWord::bData(), L1GlobalTriggerReadoutRecord::decisionWord(), doExtra_, egTriggerAlgos_, relativeConstraints::empty, ecaldqm::MESet::fill(), edm::EventSetup::get(), edm::Event::getByToken(), L1MuGMTReadoutCollection::getRecords(), L1GlobalTriggerReadoutRecord::gtPsbWord(), kCSCTrigger, kDTTrigger, kEcalTrigger, kHcalTrigger, kRPCTrigger, L1GlobalTriggerReadoutRecordToken_, L1MuGMTReadoutCollectionToken_, ecaldqm::DQWorker::MEs_, nTriggerTypes, and triggered_.

Referenced by ~ClusterTask().

   {
     if(!doExtra_) return;
 
     triggered_.reset();
 
     // TODO IS THIS THE MOST UP-TO-DATE ACCESSOR TO L1 INFO?
 
     edm::Handle<L1GlobalTriggerReadoutRecord> l1GTHndl;
     _evt.getByToken(L1GlobalTriggerReadoutRecordToken_, l1GTHndl);
     DecisionWord const& dWord(l1GTHndl->decisionWord());
 
     //Ecal
     edm::ESHandle<L1GtTriggerMenu> menuRcd;
     _es.get<L1GtTriggerMenuRcd>().get(menuRcd) ;
     L1GtTriggerMenu const* menu(menuRcd.product());
 
     if ( !dWord.empty())  { //protect against no L1GT in run
       for(unsigned iT(0); iT != egTriggerAlgos_.size(); ++iT){
     if(menu->gtAlgorithmResult(egTriggerAlgos_[iT], dWord)){
       triggered_.set(kEcalTrigger);
       break;
     }
       }
     }
 
     //Hcal
     bool hcal_top = false;
     bool hcal_bot = false;
     const L1GtPsbWord psb = l1GTHndl->gtPsbWord(0xbb0d, 0);
     std::vector<int> valid_phi; 
     if((psb.aData(4)&0x3f) >= 1) {valid_phi.push_back( (psb.aData(4)>>10)&0x1f ); }
     if((psb.bData(4)&0x3f) >= 1) {valid_phi.push_back( (psb.bData(4)>>10)&0x1f ); }
     if((psb.aData(5)&0x3f) >= 1) {valid_phi.push_back( (psb.aData(5)>>10)&0x1f ); }
     if((psb.bData(5)&0x3f) >= 1) {valid_phi.push_back( (psb.bData(5)>>10)&0x1f ); }
     std::vector<int>::const_iterator iphi;
     for(iphi=valid_phi.begin(); iphi!=valid_phi.end(); iphi++) {
       if(*iphi<9) hcal_top=true;
       if(*iphi>8) hcal_bot=true;
     }
     if(hcal_top && hcal_bot) triggered_.set(kHcalTrigger);
 
     //Muons
     edm::Handle<L1MuGMTReadoutCollection> l1MuHndl;
     if(!_evt.getByToken(L1MuGMTReadoutCollectionToken_, l1MuHndl)) return;
     std::vector<L1MuGMTReadoutRecord> const& records(l1MuHndl->getRecords());
 
     for(unsigned iR(0); iR != records.size(); ++iR){
       if(records[iR].getBxInEvent() != 0) continue;
 
       unsigned iC(0);
 
       //DT triggers
       std::vector<L1MuRegionalCand> dtBXCands(records[iR].getDTBXCands());
       for(iC = 0; iC != dtBXCands.size(); ++iC)
         if(!dtBXCands[iC].empty()) break;
       if(iC != dtBXCands.size()) triggered_.set(kDTTrigger);
 
       //RPC triggers
       std::vector<L1MuRegionalCand> brlRPCCands(records[iR].getBrlRPCCands());
       for(iC = 0; iC != brlRPCCands.size(); ++iC)
         if(!brlRPCCands[iC].empty()) break;
       if(iC != brlRPCCands.size()) triggered_.set(kRPCTrigger);
 
       //CSC Triggers
       std::vector<L1MuRegionalCand> cscCands(records[iR].getCSCCands());
       for(iC = 0; iC != cscCands.size(); ++iC)
         if(!cscCands[iC].empty()) break;
       if(iC != cscCands.size()) triggered_.set(kCSCTrigger);
     }
 
     if(triggered_.none()) return;
 
     MESet& meTriggers(MEs_.at("Triggers"));
     MESet& meExclusiveTriggers(MEs_.at("ExclusiveTriggers"));
 
     for(unsigned iT(0); iT != nTriggerTypes; ++iT){
       if(!triggered_[iT]) continue;
       meTriggers.fill(iT + 0.5);
       if(triggered_.count() == 1)
         meExclusiveTriggers.fill(iT + 0.5);
     }
   }

void ecaldqm::ClusterTask::endEvent	(	edm::Event const &	,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 176 of file ClusterTask.cc.

References ebHits_, and eeHits_.

Referenced by ~ClusterTask().

   {
     //    ++ievt_;
 
     ebHits_ = nullptr;
     eeHits_ = nullptr;
   }

bool ecaldqm::ClusterTask::filterRunType ( short const * _runType )

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 185 of file ClusterTask.cc.

References EcalDCCHeaderBlock::COSMIC, EcalDCCHeaderBlock::COSMICS_GLOBAL, EcalDCCHeaderBlock::COSMICS_LOCAL, EcalDCCHeaderBlock::MTCC, ecaldqm::nDCC, EcalDCCHeaderBlock::PHYSICS_GLOBAL, and EcalDCCHeaderBlock::PHYSICS_LOCAL.

Referenced by ~ClusterTask().

   {
     for(unsigned iFED(0); iFED != nDCC; iFED++){
       if(_runType[iFED] == EcalDCCHeaderBlock::COSMIC ||
          _runType[iFED] == EcalDCCHeaderBlock::MTCC ||
          _runType[iFED] == EcalDCCHeaderBlock::COSMICS_GLOBAL ||
          _runType[iFED] == EcalDCCHeaderBlock::PHYSICS_GLOBAL ||
          _runType[iFED] == EcalDCCHeaderBlock::COSMICS_LOCAL ||
          _runType[iFED] == EcalDCCHeaderBlock::PHYSICS_LOCAL) return true;
     }
 
     return false;
   }

void ecaldqm::ClusterTask::runOnBasicClusters	(	edm::View< reco::CaloCluster > const &	_bcs,
		Collections	_collection
	)

Definition at line 215 of file ClusterTask.cc.

References edm::View< T >::begin(), reco::CaloID::DET_ECAL_BARREL, reco::CaloID::DET_ECAL_ENDCAP, DetId::Ecal, EcalBarrel, EcalEndcap, edm::View< T >::end(), stringResolutionProvider_cfi::et, ecaldqm::getGeometry(), runTauDisplay::gp, triggerObjects_cff::id, ecaldqm::DQWorker::Timestamp::iLumi, gedGsfElectrons_cfi::isBarrel, ecaldqm::kEBBasicCluster, ecaldqm::DQWorker::MEs_, ecaldqm::phi(), position, findQualityFiles::size, ecaldqm::DQWorker::timestamp_, and ecaldqm::zside().

Referenced by analyze(), and ~ClusterTask().

   {
     MESet& meBCE(MEs_.at("BCE"));
     MESet& meBCEMap(MEs_.at("BCEMap"));
     MESet& meBCEMapProjEta(MEs_.at("BCEMapProjEta"));
     MESet& meBCEMapProjPhi(MEs_.at("BCEMapProjPhi"));
     MESet& meBCEtMapProjEta(MEs_.at("BCEtMapProjEta"));
     MESet& meBCEtMapProjPhi(MEs_.at("BCEtMapProjPhi"));
     MESet& meBCOccupancy(MEs_.at("BCOccupancy"));
     MESet& meBCOccupancyProjEta(MEs_.at("BCOccupancyProjEta"));
     MESet& meBCOccupancyProjPhi(MEs_.at("BCOccupancyProjPhi"));
     MESet& meBCSize(MEs_.at("BCSize"));
     MESet& meTrendBCSize(MEs_.at("TrendBCSize"));
     MESet& meBCSizeMap(MEs_.at("BCSizeMap"));
     MESet& meBCSizeMapProjEta(MEs_.at("BCSizeMapProjEta"));
     MESet& meBCSizeMapProjPhi(MEs_.at("BCSizeMapProjPhi"));
     MESet& meBCNum(MEs_.at("BCNum"));
     MESet& meTrendNBC(MEs_.at("TrendNBC"));
 
     int nBC[] = {0, 0};
     bool isBarrel(_collection == kEBBasicCluster);
 
     //    vector<reco::BasicCluster const*> lowMassCands;
 
     for(edm::View<reco::CaloCluster>::const_iterator bcItr(_bcs.begin()); bcItr != _bcs.end(); ++bcItr){
       if(bcItr->caloID().detectors() != 0){
         if(isBarrel && !bcItr->caloID().detector(reco::CaloID::DET_ECAL_BARREL)) continue;
         if(!isBarrel && !bcItr->caloID().detector(reco::CaloID::DET_ECAL_ENDCAP)) continue;
       }
 
       math::XYZPoint const& position(bcItr->position());
 
       DetId id(bcItr->seed());
       if(id.null()){
         GlobalPoint gp(position.x(), position.y(), position.z());
         CaloSubdetectorGeometry const* subgeom(getGeometry()->getSubdetectorGeometry(DetId::Ecal, isBarrel ? EcalBarrel : EcalEndcap));
 
         id = subgeom->getClosestCell(gp);
       }
 
       if(id.null() || (id.subdetId() == EcalBarrel && !isBarrel) || (id.subdetId() == EcalEndcap && isBarrel)) continue;
 
       float energy(bcItr->energy());
       float posEta(position.eta());
       float posPhi(phi(position.phi()));
       float et(energy / std::cosh(posEta));
       int subdet(isBarrel ? EcalBarrel : EcalEndcap);
       if(subdet == EcalEndcap && position.z() < 0.) subdet = -EcalEndcap;
 
       meBCE.fill(id, energy);
 
       meBCEMap.fill(id, energy);
       meBCEMapProjEta.fill(posEta, energy);
       meBCEMapProjPhi.fill(subdet, posPhi, energy);
       meBCEtMapProjEta.fill(posEta, et);
       meBCEtMapProjPhi.fill(subdet, posPhi, et);
 
       meBCOccupancy.fill(id);
       meBCOccupancyProjEta.fill(posEta);
       meBCOccupancyProjPhi.fill(subdet, posPhi);
 
       float size(bcItr->size());
 
       meBCSize.fill(id, size);
       meTrendBCSize.fill(id, double(timestamp_.iLumi), size);
 
       meBCSizeMap.fill(id, size);
       meBCSizeMapProjEta.fill(posEta, size);
       meBCSizeMapProjPhi.fill(subdet, posPhi, size);
 
       int zside(position.z() > 0 ? 1 : 0);
       nBC[zside]++;
 
       //       if(ievt_ % massCalcPrescale_ != 0) continue;
 
       //       if(energy > 10.) continue;
 
       //       EcalRecHitCollection::const_iterator hitItr(isBarrel ? ebHits_->find(id) : eeHits_->find(id));
       //       if(hitItr == (isBarrel ? ebHits_->end() : eeHits_->end())) continue;
 
       //       // cuts here must be parametrized
       //       if(hitItr->energy() < 0.5) continue;
 
       //       if(hitItr->energy() / energy > 0.95) continue;
 
       //       lowMassCands.push_back(&(*bcItr));
     }
 
     if(isBarrel){
       meBCNum.fill(EcalBarrel, nBC[0] + nBC[1]);
       meTrendNBC.fill(EcalBarrel, double(timestamp_.iLumi), nBC[0] + nBC[1]);
     }else{
       meBCNum.fill(-EcalEndcap, nBC[0]);
       meBCNum.fill(EcalEndcap, nBC[1]);
       meTrendNBC.fill(EcalEndcap, double(timestamp_.iLumi), nBC[0] + nBC[1]);
     }
 
     //     if(ievt_ % massCalcPrescale_ != 0) return;
 
     //     double const pi(3.14159265);
 
     //     for(vector<reco::BasicCluster const*>::iterator bcItr1(lowMassCands.begin()); bcItr1 != lowMassCands.end(); ++bcItr1){
     //       reco::BasicCluster const& bc1(**bcItr1);
     //       float energy1(bc1.energy());
     //       float px1(energy1 * sin(bc1.position().theta()) * cos(bc1.phi()));
     //       float py1(energy1 * sin(bc1.position().theta()) * sin(bc1.phi()));
     //       float pz1(energy1 * cos(bc1.position().theta()));
 
     //       for(vector<reco::BasicCluster const*>::iterator bcItr2(lowMassCands.begin()); bcItr2 != lowMassCands.end(); ++bcItr2){
     //  if(*bcItr1 == *bcItr2) continue;
     //  reco::BasicCluster const& bc2(**bcItr2);
     //  float energy2(bc2.energy());
     //  float px2(energy2 * sin(bc2.position().theta()) * cos(bc2.phi()));
     //  float py2(energy2 * sin(bc2.position().theta()) * sin(bc2.phi()));
     //  float pz2(energy2 * cos(bc2.position().theta()));
 
     //  float ptpair(sqrt((px1 + px2) * (px1 + px2) + (py1 + py2) * (py1 + py2)));
     //  if(ptpair < 2.5) continue;
 
     //  float epair(energy1 + energy2);
     //  float pzpair(abs(pz1 + pz2));
 
     //         float m2(epair * epair - pzpair * pzpair - ptpair * ptpair);
     //         if(m2 < 0.) continue;
  
     //  float eta(0.5 * log((epair + pzpair)/(epair - pzpair)));
     //  float phi(atan2(px1 + px2, py1 + py2));
 
     //  float iso(0.);
     //  for(reco::BasicClusterCollection::const_iterator bcItr(_bcs.begin()); bcItr != _bcs.end(); ++bcItr){
     //    float dEta(bcItr->eta() - eta);
     //    float dPhi(bcItr->phi() - phi);
     //           if(dPhi > 2. * pi) dPhi -= 2. * pi;
     //           else if(dPhi < -2. * pi) dPhi += 2. * pi;
     //    if(sqrt(dEta * dEta + dPhi * dPhi) < 0.2) iso += bcItr->energy() * sin(bcItr->position().theta());
     //  }
     //  if(iso > 0.5) continue;
 
     //  float mass(sqrt(m2));
     //  MEs_[kPi0]->fill(mass);
     //  MEs_[kJPsi]->fill(mass);
     //       }
     //     }
   }

void ecaldqm::ClusterTask::runOnRecHits	(	EcalRecHitCollection const &	_hits,
		Collections	_collection
	)

Definition at line 200 of file ClusterTask.cc.

References ebHits_, eeHits_, ecaldqm::kEBRecHit, and ecaldqm::kEERecHit.

Referenced by analyze(), and ~ClusterTask().

   {
     switch(_collection){
     case kEBRecHit:
       ebHits_ = &_hits;
       break;
     case kEERecHit:
       eeHits_ = &_hits;
       break;
     default:
       break;
     }
   }

void ecaldqm::ClusterTask::runOnSuperClusters	(	reco::SuperClusterCollection const &	_scs,
		Collections	_collection
	)

Definition at line 361 of file ClusterTask.cc.

References doExtra_, reco::e1, vertexPlots::e4, ebHits_, DetId::Ecal, EcalBarrel, EcalEndcap, eeHits_, photonPostprocessing_cfi::eMax, energyThreshold_, ecaldqm::MESetMulti::fill(), ecaldqm::getGeometry(), ecaldqm::getTopology(), runTauDisplay::gp, hfClusterShapes_cfi::hits, ecaldqm::DQWorker::Timestamp::iLumi, gedGsfElectrons_cfi::isBarrel, ecaldqm::kEBSuperCluster, ecaldqm::DQWorker::MEs_, nTriggerTypes, ecaldqm::phi(), position, findQualityFiles::size, swissCrossMaxThreshold_, ecaldqm::DQWorker::timestamp_, triggered_, and trigTypeToME_.

Referenced by analyze(), and ~ClusterTask().

   {
     bool isBarrel(_collection == kEBSuperCluster);
     EcalSubdetector subdet(isBarrel ? EcalBarrel : EcalEndcap);
 
     MESet& meSCE(MEs_.at("SCE"));
     MESet& meSCELow(MEs_.at("SCELow"));
     MESet& meSCNBCs(MEs_.at("SCNBCs"));
     MESet& meSCNcrystals(MEs_.at("SCNcrystals"));
     MESet& meTrendSCSize(MEs_.at("TrendSCSize"));
     MESet& meSCSeedEnergy(MEs_.at("SCSeedEnergy"));
     MESet& meSCClusterVsSeed(MEs_.at("SCClusterVsSeed"));
     MESet& meSCSeedOccupancy(MEs_.at("SCSeedOccupancy"));
     MESet& meSingleCrystalCluster(MEs_.at("SingleCrystalCluster"));
     MESet& meSCR9(MEs_.at("SCR9"));
 
     MESet* meSCSizeVsEnergy(doExtra_ ? &MEs_.at("SCSizeVsEnergy") : nullptr);
     MESet* meSCSeedOccupancyHighE(doExtra_ ? &MEs_.at("SCSeedOccupancyHighE") : nullptr);
     MESet* meSCSeedOccupancyTrig(doExtra_ ? &MEs_.at("SCSeedOccupancyTrig") : nullptr);
     MESet* meSCSeedTimeTrigEx(doExtra_ ? &MEs_.at("SCSeedTimeTrigEx") : nullptr);
     MESet* meSCSeedTimeMapTrigEx(doExtra_ ? &MEs_.at("SCSeedTimeMapTrigEx") : nullptr);
     MESet* meSCOccupancyProjEta(doExtra_ ? &MEs_.at("SCOccupancyProjEta") : nullptr);
     MESet* meSCOccupancyProjPhi(doExtra_ ? &MEs_.at("SCOccupancyProjPhi") : nullptr);
     MESet* meSCSwissCross(doExtra_ ? &MEs_.at("SCSwissCross") : nullptr);
 
     EcalRecHitCollection const* hits(isBarrel ? ebHits_ : eeHits_);
 
     //     reco::SuperCluster const* leading(0);
     //     reco::SuperCluster const* subLeading(0);
 
     int nSC(0);
 
     for(reco::SuperClusterCollection::const_iterator scItr(_scs.begin()); scItr != _scs.end(); ++scItr){
       DetId seedId(scItr->seed()->seed());
       if(seedId.null()){
         math::XYZPoint const& position(scItr->position());
 
     GlobalPoint gp(position.x(), position.y(), position.z());
 
     CaloSubdetectorGeometry const* subgeom(getGeometry()->getSubdetectorGeometry(DetId::Ecal, isBarrel ? EcalBarrel : EcalEndcap));
 
     seedId = subgeom->getClosestCell(gp);
       }
 
       if(seedId.null() || (seedId.subdetId() != subdet)) continue;
 
       EcalRecHitCollection::const_iterator seedItr(hits->find(seedId));
       if(seedItr == hits->end()) continue;
 
       ++nSC;
 
       float energy(scItr->energy());
       float size(scItr->size());
 
       meSCE.fill(seedId, energy);
       meSCELow.fill(seedId, energy);
 
       meSCNBCs.fill(seedId, scItr->clustersSize());
       meSCNcrystals.fill(seedId, size);
 
       if(doExtra_) meSCSizeVsEnergy->fill(subdet, energy, size);
 
       meTrendSCSize.fill(seedId, double(timestamp_.iLumi), size);
 
       meSCSeedEnergy.fill(seedId, seedItr->energy());
       meSCClusterVsSeed.fill(seedId, seedItr->energy(), energy);
 
       meSCSeedOccupancy.fill(seedId);
       if(doExtra_ && energy > energyThreshold_) meSCSeedOccupancyHighE->fill(seedId);
 
       if(scItr->size() == 1) meSingleCrystalCluster.fill(seedId);
 
       float e3x3(EcalClusterTools::e3x3(*scItr->seed(), hits, getTopology()));
       meSCR9.fill(seedId, e3x3 / energy);
 
       if(doExtra_){
         for(unsigned iT(0); iT != nTriggerTypes; ++iT){
           if(!triggered_[iT]) continue;
 
           static_cast<MESetMulti*>(meSCSeedOccupancyTrig)->use(trigTypeToME_[iT]);
           meSCSeedOccupancyTrig->fill(seedId);
 
           // exclusive
           if(triggered_.count() == 1){
             static_cast<MESetMulti*>(meSCSeedTimeTrigEx)->use(trigTypeToME_[iT]);
             static_cast<MESetMulti*>(meSCSeedTimeMapTrigEx)->use(trigTypeToME_[iT]);
             meSCSeedTimeTrigEx->fill(subdet, seedItr->time());
             meSCSeedTimeMapTrigEx->fill(seedId, seedItr->time());
           }
         }
 
         meSCOccupancyProjEta->fill(subdet, scItr->eta());
         meSCOccupancyProjPhi->fill(subdet, phi(scItr->phi()));
 
         if(isBarrel){
           float e1(EcalClusterTools::eMax(*scItr, ebHits_));
           if(e1 > swissCrossMaxThreshold_){
             float e4(EcalClusterTools::eTop(*scItr, ebHits_, getTopology()) +
                      EcalClusterTools::eRight(*scItr, ebHits_, getTopology()) +
                      EcalClusterTools::eBottom(*scItr, ebHits_, getTopology()) +
                      EcalClusterTools::eLeft(*scItr, ebHits_, getTopology()));
 
             meSCSwissCross->fill(1. - e4 / e1);
           }
         }
       }
 
       //       if(ievt_ % massCalcPrescale_ != 0) continue;
 
       //       float et(energy * sin(scItr->position().theta()));
       //       if(!leading || et > leading->energy() * sin(leading->position().theta())){
       //  subLeading = leading;
       //  leading = &(*scItr);
       //       }
       //       else if(!subLeading || et > subLeading->energy() * sin(subLeading->position().theta())){
       //  subLeading = &(*scItr);
       //       }
     }
 
     MEs_.at("SCNum").fill(subdet, nSC);
     MEs_.at("TrendNSC").fill(subdet, double(timestamp_.iLumi), nSC);
 
     //     if(ievt_ % massCalcPrescale_ != 0) return;
 
     //     // implement isolation & cuts
     //     if(!leading || !subLeading) return;
     //     float e(leading->energy() + subLeading->energy());
     //     float px(leading->energy() * sin(leading->position().theta()) * cos(leading->phi()) + subLeading->energy() * sin(subLeading->position().theta()) * cos(subLeading->phi()));
     //     float py(leading->energy() * sin(leading->position().theta()) * sin(leading->phi()) + subLeading->energy() * sin(subLeading->position().theta()) * sin(subLeading->phi()));
     //     float pz(leading->energy() * cos(leading->position().theta()) + subLeading->energy() * cos(subLeading->position().theta()));
     //     float m2(e * e - px * px - py * py - pz * pz);
     //     if(m2 < 0.) return;
     //     float mass(sqrt(m2));
     //     MEs_[kZ]->fill(mass);
     //     MEs_[kHighMass]->fill(mass);
 
   }

void ecaldqm::ClusterTask::setParams ( edm::ParameterSet const & _params )

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 41 of file ClusterTask.cc.

References doExtra_, egTriggerAlgos_, energyThreshold_, ecaldqm::MESetMulti::getIndex(), edm::ParameterSet::getUntrackedParameter(), L1GlobalTriggerReadoutRecordTag_, L1MuGMTReadoutCollectionTag_, ecaldqm::DQWorker::MEs_, nTriggerTypes, AlCaHLTBitMon_QueryRunRegistry::string, swissCrossMaxThreshold_, and trigTypeToME_.

   {
     doExtra_ = _params.getUntrackedParameter<bool>("doExtra");
 
     if(!doExtra_){
       MEs_.erase(std::string("SCSizeVsEnergy"));
       MEs_.erase(std::string("SCSeedOccupancyHighE"));
       MEs_.erase(std::string("SCSeedOccupancyTrig"));
       MEs_.erase(std::string("SCSeedTimeTrigEx"));
       MEs_.erase(std::string("SCSeedTimeMapTrigEx"));
       MEs_.erase(std::string("SCOccupancyProjEta"));
       MEs_.erase(std::string("SCOccupancyProjPhi"));
       MEs_.erase(std::string("SCSwissCross"));
       MEs_.erase(std::string("Triggers"));
       MEs_.erase(std::string("ExclusiveTriggers"));
 
       return;
     }
 
     energyThreshold_ = _params.getUntrackedParameter<double>("energyThreshold");
     swissCrossMaxThreshold_ = _params.getUntrackedParameter<double>("swissCrossMaxThreshold");
     egTriggerAlgos_ = _params.getUntrackedParameter<std::vector<std::string> >("egTriggerAlgos");
     L1GlobalTriggerReadoutRecordTag_ = _params.getUntrackedParameter<edm::InputTag>("L1GlobalTriggerReadoutRecordTag");
     L1MuGMTReadoutCollectionTag_ = _params.getUntrackedParameter<edm::InputTag>("L1MuGMTReadoutCollectionTag");
 
     MESet::PathReplacements repl;
 
     std::string triggerTypeNames[nTriggerTypes] = {
       "ECAL",
       "HCAL",
       "CSC",
       "DT",
       "RPC"
     };
 
     MESetMulti& occupancy(static_cast<MESetMulti&>(MEs_.at("SCSeedOccupancyTrig")));
     for(unsigned iT(0); iT != nTriggerTypes; ++iT){
       repl["trig"] = triggerTypeNames[iT];
       trigTypeToME_[iT] = occupancy.getIndex(repl);
     }
   }

void ecaldqm::ClusterTask::setTokens ( edm::ConsumesCollector & _collector )

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 500 of file ClusterTask.cc.

References edm::ConsumesCollector::consumes(), DEFINE_ECALDQM_WORKER, L1GlobalTriggerReadoutRecordTag_, L1GlobalTriggerReadoutRecordToken_, L1MuGMTReadoutCollectionTag_, and L1MuGMTReadoutCollectionToken_.

Referenced by ~ClusterTask().

   {
     L1GlobalTriggerReadoutRecordToken_ = _collector.consumes<L1GlobalTriggerReadoutRecord>(L1GlobalTriggerReadoutRecordTag_);
     L1MuGMTReadoutCollectionToken_ = _collector.consumes<L1MuGMTReadoutCollection>(L1MuGMTReadoutCollectionTag_);
   }