#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 &)

	~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 &)

const std::string &	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 () noexcept(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 Types inherited from ecaldqm::DQWorker
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

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::DQWorker
bool	booked_

MESetCollection	MEs_

std::string	name_

bool	onlineMode_

Timestamp	timestamp_

int	verbosity_

bool	willConvertToEDM_

Detailed Description

Definition at line 20 of file ClusterTask.h.

Member Enumeration Documentation

◆ TriggerTypes

enum ecaldqm::ClusterTask::TriggerTypes

Enumerator
kEcalTrigger
kHcalTrigger
kCSCTrigger
kDTTrigger
kRPCTrigger
nTriggerTypes

Definition at line 40 of file ClusterTask.h.

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

Constructor & Destructor Documentation

◆ ClusterTask()

ecaldqm::ClusterTask::ClusterTask ( )

Definition at line 21 of file ClusterTask.cc.

       : 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_() {}

◆ ~ClusterTask()

ecaldqm::ClusterTask::~ClusterTask ( )

inlineoverride

Definition at line 23 of file ClusterTask.h.

23 {}

Member Function Documentation

◆ addDependencies()

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 72 of file ClusterTask.cc.

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

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

◆ analyze()

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

inlineoverridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 62 of file ClusterTask.h.

                                                                           {
     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;
   }

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

◆ beginEvent()

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 77 of file ClusterTask.cc.

                                                                              {
     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);
     }
   }

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

◆ endEvent()

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 182 of file ClusterTask.cc.

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

References ebHits_, and eeHits_.

◆ filterRunType()

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 189 of file ClusterTask.cc.

                                                        {
     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;
   }

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

◆ runOnBasicClusters()

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

Definition at line 214 of file ClusterTask.cc.

                                                                                                     {
     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);
     //       }
     //     }
   }

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

Referenced by analyze().

◆ runOnRecHits()

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

Definition at line 201 of file ClusterTask.cc.

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

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

Referenced by analyze().

◆ runOnSuperClusters()

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

Definition at line 363 of file ClusterTask.cc.

                                                                                                       {
     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);
   }

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

Referenced by analyze().

◆ setParams()

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 37 of file ClusterTask.cc.

                                                             {
     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);
     }
   }

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

◆ setTokens()

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 506 of file ClusterTask.cc.

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

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

Member Data Documentation

◆ doExtra_

bool ecaldqm::ClusterTask::doExtra_

private

Definition at line 49 of file ClusterTask.h.

Referenced by beginEvent(), runOnSuperClusters(), and setParams().

◆ ebHits_

EcalRecHitCollection const* ecaldqm::ClusterTask::ebHits_

private

Definition at line 45 of file ClusterTask.h.

Referenced by endEvent(), runOnRecHits(), and runOnSuperClusters().

◆ eeHits_

EcalRecHitCollection const* ecaldqm::ClusterTask::eeHits_

private

Definition at line 46 of file ClusterTask.h.

Referenced by endEvent(), runOnRecHits(), and runOnSuperClusters().

◆ egTriggerAlgos_

std::vector<std::string> ecaldqm::ClusterTask::egTriggerAlgos_

private

Definition at line 52 of file ClusterTask.h.

Referenced by beginEvent(), and setParams().

◆ energyThreshold_

float ecaldqm::ClusterTask::energyThreshold_

private

Definition at line 50 of file ClusterTask.h.

Referenced by runOnSuperClusters(), and setParams().

◆ L1GlobalTriggerReadoutRecordTag_

edm::InputTag ecaldqm::ClusterTask::L1GlobalTriggerReadoutRecordTag_

private

Definition at line 56 of file ClusterTask.h.

Referenced by setParams(), and setTokens().

◆ L1GlobalTriggerReadoutRecordToken_

edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> ecaldqm::ClusterTask::L1GlobalTriggerReadoutRecordToken_

private

Definition at line 58 of file ClusterTask.h.

Referenced by beginEvent(), and setTokens().

◆ L1MuGMTReadoutCollectionTag_

edm::InputTag ecaldqm::ClusterTask::L1MuGMTReadoutCollectionTag_

private

Definition at line 57 of file ClusterTask.h.

Referenced by setParams(), and setTokens().

◆ L1MuGMTReadoutCollectionToken_

edm::EDGetTokenT<L1MuGMTReadoutCollection> ecaldqm::ClusterTask::L1MuGMTReadoutCollectionToken_

private

Definition at line 59 of file ClusterTask.h.

Referenced by beginEvent(), and setTokens().

◆ swissCrossMaxThreshold_

float ecaldqm::ClusterTask::swissCrossMaxThreshold_

private

Definition at line 51 of file ClusterTask.h.

Referenced by runOnSuperClusters(), and setParams().

◆ triggered_

std::bitset<nTriggerTypes> ecaldqm::ClusterTask::triggered_

private

Definition at line 53 of file ClusterTask.h.

Referenced by beginEvent(), and runOnSuperClusters().

◆ trigTypeToME_

unsigned ecaldqm::ClusterTask::trigTypeToME_[nTriggerTypes]

private

Definition at line 54 of file ClusterTask.h.

Referenced by runOnSuperClusters(), and setParams().

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ TriggerTypes

Constructor & Destructor Documentation

◆ ClusterTask()

◆ ~ClusterTask()

Member Function Documentation

◆ addDependencies()

◆ analyze()

◆ beginEvent()

◆ endEvent()

◆ filterRunType()

◆ runOnBasicClusters()

◆ runOnRecHits()

◆ runOnSuperClusters()

◆ setParams()

◆ setTokens()

Member Data Documentation

◆ doExtra_

◆ ebHits_

◆ eeHits_

◆ egTriggerAlgos_

◆ energyThreshold_

◆ L1GlobalTriggerReadoutRecordTag_

◆ L1GlobalTriggerReadoutRecordToken_

◆ L1MuGMTReadoutCollectionTag_

◆ L1MuGMTReadoutCollectionToken_

◆ swissCrossMaxThreshold_

◆ triggered_

◆ trigTypeToME_