#include <TrigPrimTask.h>

Inheritance diagram for ecaldqm::TrigPrimTask:

Public Types
enum	Constants { nBXBins = 15 }

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

void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

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

template<typename DigiCollection >
void	runOnDigis (DigiCollection const &)

void	runOnEmulTPs (EcalTrigPrimDigiCollection const &)

void	runOnRealTPs (EcalTrigPrimDigiCollection const &)

	TrigPrimTask ()

	~TrigPrimTask ()

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

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

virtual bool	filterRunType (short const *)

virtual bool	filterTrigger (edm::TriggerResultsByName const &)

void	recoverStats ()

virtual void	setTokens (edm::ConsumesCollector &)

void	softReset ()

virtual	~DQWorkerTask ()

Public Member Functions inherited from ecaldqm::DQWorker
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 ()

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

Private Attributes
double	bxBin_

int	bxBinEdges_ [nBXBins+1]

EcalTrigPrimDigiCollection const *	realTps_

bool	runOnEmul_

edm::ESHandle< EcalTPGStripStatus >	StripStatusRcd

std::map< uint32_t, unsigned >	towerReadouts_

edm::ESHandle< EcalTPGTowerStatus >	TTStatusRcd

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 15 of file TrigPrimTask.h.

Member Enumeration Documentation

enum ecaldqm::TrigPrimTask::Constants

Enumerator
nBXBins

Definition at line 32 of file TrigPrimTask.h.

                    {
       nBXBins = 15
     };

Constructor & Destructor Documentation

ecaldqm::TrigPrimTask::TrigPrimTask ( )

Definition at line 17 of file TrigPrimTask.cc.

                              :
     DQWorkerTask(),
     realTps_(0),
     runOnEmul_(false),
     //     HLTCaloPath_(""),
     //     HLTMuonPath_(""),
     //     HLTCaloBit_(false),
     //     HLTMuonBit_(false),
     bxBinEdges_{1, 271, 541, 892, 1162, 1432, 1783, 2053, 2323, 2674, 2944, 3214, 3446, 3490, 3491, 3565},
     bxBin_(0.),
     towerReadouts_()
   {
   }

ecaldqm::TrigPrimTask::~TrigPrimTask ( )

inline

Definition at line 18 of file TrigPrimTask.h.

18 {}

Member Function Documentation

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 45 of file TrigPrimTask.cc.

References ecaldqm::kEBDigi, ecaldqm::kEEDigi, ecaldqm::kTrigPrimDigi, ecaldqm::kTrigPrimEmulDigi, ecaldqm::DependencySet::push_back(), and runOnEmul_.

   {
     if(runOnEmul_) _dependencies.push_back(Dependency(kTrigPrimEmulDigi, kEBDigi, kEEDigi, kTrigPrimDigi));
   }

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

inlineoverridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 58 of file TrigPrimTask.h.

References ecaldqm::kEBDigi, ecaldqm::kEEDigi, ecaldqm::kTrigPrimDigi, ecaldqm::kTrigPrimEmulDigi, runOnDigis(), runOnEmul_, runOnEmulTPs(), and runOnRealTPs().

                                                                           {
     switch(_collection){
     case kTrigPrimDigi:
       if(_p) runOnRealTPs(*static_cast<EcalTrigPrimDigiCollection const*>(_p));
       return true;
       break;
     case kTrigPrimEmulDigi:
       if(_p && runOnEmul_) runOnEmulTPs(*static_cast<EcalTrigPrimDigiCollection const*>(_p));
       return runOnEmul_;
       break;
     case kEBDigi:
       if(_p) runOnDigis(*static_cast<EBDigiCollection const*>(_p));
       return true;
       break;
     case kEEDigi:
       if(_p) runOnDigis(*static_cast<EEDigiCollection const*>(_p));
       return true;
       break;
     default:
       break;
     }
     return false;
   }

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 68 of file TrigPrimTask.cc.

References edm::EventBase::bunchCrossing(), bxBin_, bxBinEdges_, edm::EventSetup::get(), EcalTPGStripStatus::getMap(), EcalTPGTowerStatus::getMap(), ecaldqm::DQWorker::MEs_, nBXBins, realTps_, and towerReadouts_.

   {
     using namespace std;
 
     towerReadouts_.clear();
 
     realTps_ = 0;
 
     //     HLTCaloBit_ = false;
     //     HLTMuonBit_ = false;
 
     int* pBin(std::upper_bound(bxBinEdges_, bxBinEdges_ + nBXBins + 1, _evt.bunchCrossing()));
     bxBin_ = static_cast<int>(pBin - bxBinEdges_) - 0.5;
 
     edm::ESHandle<EcalTPGTowerStatus> TTStatusRcd_;
     _es.get<EcalTPGTowerStatusRcd>().get(TTStatusRcd_);
     const EcalTPGTowerStatus * TTStatus=TTStatusRcd_.product();
     const EcalTPGTowerStatusMap &towerMap=TTStatus->getMap();
 
     edm::ESHandle<EcalTPGStripStatus> StripStatusRcd_;
     _es.get<EcalTPGStripStatusRcd>().get(StripStatusRcd_);
     const EcalTPGStripStatus * StripStatus=StripStatusRcd_.product();
     const EcalTPGStripStatusMap &stripMap=StripStatus->getMap();
 
     MESet& meTTMaskMap(MEs_.at("TTMaskMap"));
 
     for(EcalTPGTowerStatusMap::const_iterator ttItr(towerMap.begin()); ttItr != towerMap.end(); ++ttItr){
 
        if ((*ttItr).second > 0)
        {
          const EcalTrigTowerDetId  ttid((*ttItr).first);
          //if(ttid.subDet() == EcalBarrel)
             meTTMaskMap.fill(ttid,1);
        }//masked   
     }//loop on towers
   
     for(EcalTPGStripStatusMap::const_iterator stItr(stripMap.begin()); stItr != stripMap.end(); ++stItr){
 
        if ((*stItr).second > 0)
        {
          const EcalElectronicsId stid((*stItr).first);
          //if(stid.subdet() == EcalEndcap);
             meTTMaskMap.fill(stid,1);
        }//masked   
     }//loop on pseudo-strips
   
     //     if(HLTCaloPath_.size() || HLTMuonPath_.size()){
     //       edm::TriggerResultsByName results(_evt.triggerResultsByName("HLT"));
     //       if(!results.isValid()) results = _evt.triggerResultsByName("RECO");
     //       if(results.isValid()){
     //  const vector<string>& pathNames(results.triggerNames());
 
     //  size_t caloStar(HLTCaloPath_.find('*'));
     //  if(caloStar != string::npos){
     //    string caloSub(HLTCaloPath_.substr(0, caloStar));
     //    bool found(false);
     //    for(unsigned iP(0); iP < pathNames.size(); ++iP){
     //      if(pathNames[iP].substr(0, caloStar) == caloSub){
     //        HLTCaloPath_ = pathNames[iP];
     //        found = true;
     //        break;
     //      }
     //    }
     //    if(!found) HLTCaloPath_ = "";
     //  }
 
     //  size_t muonStar(HLTMuonPath_.find('*'));
     //  if(muonStar != string::npos){
     //    string muonSub(HLTMuonPath_.substr(0, muonStar));
     //    bool found(false);
     //    for(unsigned iP(0); iP < pathNames.size(); ++iP){
     //      if(pathNames[iP].substr(0, muonStar) == muonSub){
     //        HLTMuonPath_ = pathNames[iP];
     //        found = true;
     //        break;
     //      }
     //    }
     //    if(!found) HLTMuonPath_ = "";
     //  }
 
     //  if(HLTCaloPath_.size()){
     //    try{
     //      HLTCaloBit_ = results.accept(HLTCaloPath_);
     //    }
     //    catch(cms::Exception e){
     //      if(e.category() != "LogicError") throw e;
     //      HLTCaloPath_ = "";
     //    }
     //  }
     //  if(HLTMuonPath_.size()){
     //    try{
     //      HLTMuonBit_ = results.accept(HLTMuonPath_);
     //    }
     //    catch(cms::Exception e){
     //      if(e.category() != "LogicError") throw e;
     //      HLTMuonPath_ = "";
     //    }
     //  }
     //       }
     //     }
   }

void ecaldqm::TrigPrimTask::beginLuminosityBlock	(	edm::LuminosityBlock const &	,
		edm::EventSetup const &
	)

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 61 of file TrigPrimTask.cc.

References ecaldqm::DQWorker::MEs_.

   {
     // Reset by LS plots at beginning of every LS
     MEs_.at("EtSummaryByLumi").reset();
   }

void ecaldqm::TrigPrimTask::beginRun	(	edm::Run const &	,
		edm::EventSetup const &	_es
	)

overridevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 51 of file TrigPrimTask.cc.

References edm::EventSetup::get(), StripStatusRcd, and TTStatusRcd.

   {
     // Read-in Status records:
     // Status records stay constant over run so they are read-in only once here
     // but filled by LS in runOnRealTPs() because MEs are not yet booked at beginRun()
     _es.get<EcalTPGTowerStatusRcd>().get( TTStatusRcd );
     _es.get<EcalTPGStripStatusRcd>().get( StripStatusRcd );
   }

template<typename DigiCollection >

void ecaldqm::TrigPrimTask::runOnDigis ( DigiCollection const & _digis )

Definition at line 172 of file TrigPrimTask.cc.

References ecaldqm::getTrigTowerMap(), DetId::rawId(), and towerReadouts_.

Referenced by analyze().

   {
     for(typename DigiCollection::const_iterator digiItr(_digis.begin()); digiItr != _digis.end(); ++digiItr){
       EcalTrigTowerDetId ttid(getTrigTowerMap()->towerOf(digiItr->id()));
       towerReadouts_[ttid.rawId()]++;
     }
   }

void ecaldqm::TrigPrimTask::runOnEmulTPs ( EcalTrigPrimDigiCollection const & _tps )

Definition at line 282 of file TrigPrimTask.cc.

References edm::SortedCollection< T, SORT >::begin(), EcalTrigTowerConstituentsMap::constituentsOf(), edm::SortedCollection< T, SORT >::end(), edm::SortedCollection< T, SORT >::find(), ecaldqm::getTrigTowerMap(), match(), ecaldqm::DQWorker::MEs_, realTps_, and towerReadouts_.

Referenced by analyze().

   {
     MESet& meEtMaxEmul(MEs_.at("EtMaxEmul"));
     MESet& meEmulMaxIndex(MEs_.at("EmulMaxIndex"));
     MESet& meMatchedIndex(MEs_.at("MatchedIndex"));
     MESet& meEtEmulError(MEs_.at("EtEmulError"));
     MESet& meFGEmulError(MEs_.at("FGEmulError"));
     MESet& meRealvEmulEt(MEs_.at("RealvEmulEt"));
 
     for(EcalTrigPrimDigiCollection::const_iterator tpItr(_tps.begin()); tpItr != _tps.end(); ++tpItr){
       EcalTrigTowerDetId ttid(tpItr->id());
 
       int et(tpItr->compressedEt());
 
       float maxEt(0.);
       int iMax(0);
       for(int iDigi(0); iDigi < 5; iDigi++){
         float sampleEt((*tpItr)[iDigi].compressedEt());
 
         if(sampleEt > maxEt){
           maxEt = sampleEt;
           iMax = iDigi + 1;
         }
       }
 
       meEtMaxEmul.fill(ttid, maxEt);
       if(maxEt > 0.)
         meEmulMaxIndex.fill(ttid, iMax);
 
       bool match(true);
       bool matchFG(true);
 
       // Loop over real TPs and look for an emulated TP index with matching Et:
       // If an Et match is found, return TP index correpsonding to BX of emulated TP where match was found
       // Standard TPG comparison: { TP index:matched BX } = { no emul:No Et match, 0:BX-2, 1:BX-1, 2:in-time, 3:BX+1, 4:BX+2 }
       EcalTrigPrimDigiCollection::const_iterator realItr(realTps_->find(ttid));
       if(realItr != realTps_->end()){
 
         int realEt(realItr->compressedEt());
 
         if(realEt > 0){
 
           int interest(realItr->ttFlag() & 0x3);
           if((interest == 1 || interest == 3) && towerReadouts_[ttid.rawId()] == getTrigTowerMap()->constituentsOf(ttid).size()){
 
             if(et != realEt) match = false;
             if(tpItr->fineGrain() != realItr->fineGrain()) matchFG = false;
 
         // NOTE: matchedIndex comparison differs from Standard TPG comparison:
         // { matchedIndex:TP index } = { 0:no emul, 1:BX-2, 2:BX-1, 3:in-time, 4:BX+1, 5:BX+2 }
         std::vector<int> matchedIndex(0);
         // iDigi only loops over explicit Et matches:
         // { iDigi:TP index } = { 0:BX-2, 1:BX-1, 2:in-time, 3:BX+1, 4:BX+2 }
         for(int iDigi(0); iDigi < 5; iDigi++){
           if((*tpItr)[iDigi].compressedEt() == realEt) {
         // matchedIndex = iDigi + 1
         if (iDigi != 2) {
           matchedIndex.push_back(iDigi + 1);
         }
         // If an in-time match is found, exit loop and clear out any other matches:
         // Ensures multiple matches are not returned (e.g. during saturation)
         else {
           matchedIndex.clear();
           matchedIndex.push_back(3); // Et match is to in-time emulated TP
           break;
         }
           } // Et match found
         } // iDigi
         if(!matchedIndex.size()) matchedIndex.push_back(0); // no Et match found => no emul
 
         // Fill Real vs Emulated TP Et
         meRealvEmulEt.fill( ttid,realEt,(*tpItr)[2].compressedEt() ); // iDigi=2:in-time BX
 
         // Fill matchedIndex ME
         for(std::vector<int>::iterator matchItr(matchedIndex.begin()); matchItr != matchedIndex.end(); ++matchItr){
               meMatchedIndex.fill(ttid, *matchItr + 0.5);
 
               // timing information is only within emulated TPs (real TPs have one time sample)
               //      if(HLTCaloBit_) MEs_[kTimingCalo].fill(ttid, float(*matchItr));
               //      if(HLTMuonBit_) MEs_[kTimingMuon].fill(ttid, float(*matchItr));
             }
           }
 
         }
       }
       else{
         match = false;
         matchFG = false;
       }
 
       if(!match) meEtEmulError.fill(ttid);
       if(!matchFG) meFGEmulError.fill(ttid);
     }
   }

void ecaldqm::TrigPrimTask::runOnRealTPs ( EcalTrigPrimDigiCollection const & _tps )

Definition at line 181 of file TrigPrimTask.cc.

References edm::SortedCollection< T, SORT >::begin(), bxBin_, EcalTrigTowerConstituentsMap::constituentsOf(), EcalBarrel, EcalEndcap, edm::SortedCollection< T, SORT >::end(), ecaldqm::getElectronicsMap(), ecaldqm::getTrigTowerMap(), ecaldqm::DQWorker::MEs_, edm::ESHandle< class >::product(), realTps_, StripStatusRcd, EcalTriggerElectronicsId::tccId(), towerReadouts_, EcalTriggerElectronicsId::ttId(), and TTStatusRcd.

Referenced by analyze().

   {
     MESet& meEtVsBx(MEs_.at("EtVsBx"));
     MESet& meEtReal(MEs_.at("EtReal"));
     MESet& meEtRealMap(MEs_.at("EtRealMap"));
     MESet& meEtSummary(MEs_.at("EtSummary"));
     MESet& meEtSummaryByLumi(MEs_.at("EtSummaryByLumi"));
     MESet& meLowIntMap(MEs_.at("LowIntMap"));
     MESet& meMedIntMap(MEs_.at("MedIntMap"));
     MESet& meHighIntMap(MEs_.at("HighIntMap"));
     MESet& meTTFlags(MEs_.at("TTFlags"));
     MESet& meTTFlagsVsEt(MEs_.at("TTFlagsVsEt"));
     MESet& meTTFlags4( MEs_.at("TTFlags4") );
     MESet& meTTFMismatch(MEs_.at("TTFMismatch"));
     MESet& meOccVsBx(MEs_.at("OccVsBx"));
 
     realTps_ = &_tps;
 
     double nTP[] = {0., 0., 0.};
 
     for(EcalTrigPrimDigiCollection::const_iterator tpItr(_tps.begin()); tpItr != _tps.end(); ++tpItr){
       EcalTrigTowerDetId ttid(tpItr->id());
       float et(tpItr->compressedEt());
 
       if(et > 0.){
         if(ttid.subDet() == EcalBarrel)
           nTP[0] += 1.;
         else if(ttid.zside() < 0)
           nTP[1] += 1.;
         else
           nTP[2] += 2.;
         meEtVsBx.fill(ttid, bxBin_, et);
       }
 
       meEtReal.fill(ttid, et);
       meEtRealMap.fill(ttid, et);
       meEtSummary.fill(ttid, et);
       meEtSummaryByLumi.fill(ttid, et);
 
       int interest(tpItr->ttFlag() & 0x3);
 
       switch(interest){
       case 0:
         meLowIntMap.fill(ttid);
         break;
       case 1:
         meMedIntMap.fill(ttid);
         break;
       case 3:
         meHighIntMap.fill(ttid);
         break;
       default:
         break;
       }
 
       // Fill TT Flag MEs
       float ttF( tpItr->ttFlag() );
       meTTFlags.fill( ttid, ttF );
       meTTFlagsVsEt.fill(ttid, et, ttF);
       // Monitor occupancy of TTF=4
       // which contains info about TT auto-masking
       if ( ttF == 4. )
         meTTFlags4.fill( ttid );
 
       if((interest == 1 || interest == 3) && towerReadouts_[ttid.rawId()] != getTrigTowerMap()->constituentsOf(ttid).size())
         meTTFMismatch.fill(ttid);
     }
 
     meOccVsBx.fill( EcalBarrel, bxBin_, nTP[0]);
     meOccVsBx.fill(-EcalEndcap, bxBin_, nTP[1]);
     meOccVsBx.fill( EcalEndcap, bxBin_, nTP[2]);
 
     // Set TT/Strip Masking status in Ecal3P view
     // Status Records are read-in at beginRun() but filled here
     // Requestied by ECAL Trigger in addition to TTMaskMap plots in SM view
     MESet& meTTMaskMapAll(MEs_.at("TTMaskMapAll"));
 
     // Fill from TT Status Rcd
     const EcalTPGTowerStatus *TTStatus( TTStatusRcd.product() );
     const EcalTPGTowerStatusMap &TTStatusMap( TTStatus->getMap() );
     for( EcalTPGTowerStatusMap::const_iterator ttItr(TTStatusMap.begin()); ttItr != TTStatusMap.end(); ++ttItr ){
       const EcalTrigTowerDetId ttid( ttItr->first );
       if ( ttItr->second > 0 )
         meTTMaskMapAll.setBinContent( ttid,1 ); // TT is masked
     } // TTs
 
     // Fill from Strip Status Rcd
     const EcalTPGStripStatus *StripStatus( StripStatusRcd.product() );
     const EcalTPGStripStatusMap &StripStatusMap( StripStatus->getMap() );
     for( EcalTPGStripStatusMap::const_iterator stItr(StripStatusMap.begin()); stItr != StripStatusMap.end(); ++stItr ){
       const EcalTriggerElectronicsId stid( stItr->first );
       // Since ME has kTriggerTower binning, convert to EcalTrigTowerDetId first
       // In principle, setBinContent() could be implemented for EcalTriggerElectronicsId class as well
       const EcalTrigTowerDetId ttid( getElectronicsMap()->getTrigTowerDetId(stid.tccId(), stid.ttId()) );
       if ( stItr->second > 0 )
         meTTMaskMapAll.setBinContent( ttid,1 ); // PseudoStrip is masked
     } // PseudoStrips
 
   } // TrigPrimTask::runOnRealTPs()

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 32 of file TrigPrimTask.cc.

References edm::ParameterSet::getUntrackedParameter(), ecaldqm::DQWorker::MEs_, runOnEmul_, and AlCaHLTBitMon_QueryRunRegistry::string.

   {
     runOnEmul_ = _params.getUntrackedParameter<bool>("runOnEmul");
     if(!runOnEmul_){
       MEs_.erase(std::string("EtMaxEmul"));
       MEs_.erase(std::string("EmulMaxIndex"));
       MEs_.erase(std::string("MatchedIndex"));
       MEs_.erase(std::string("EtEmulError"));
       MEs_.erase(std::string("FGEmulError"));
     }
   }