#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

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	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 ()

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

Private Attributes
double	bxBin_

int	bxBinEdges_ [nBXBins+1]

EcalTrigPrimDigiCollection const *	realTps_

bool	runOnEmul_

std::map< uint32_t, unsigned >	towerReadouts_

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

Member Enumeration Documentation

enum ecaldqm::TrigPrimTask::Constants

Enumerator
nBXBins

Definition at line 25 of file TrigPrimTask.h.

                    {
       nBXBins = 15
     };

Constructor & Destructor Documentation

ecaldqm::TrigPrimTask::TrigPrimTask ( )

Definition at line 21 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 13 of file TrigPrimTask.h.

13 {}

Member Function Documentation

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

overridevirtual

Reimplemented from ecaldqm::DQWorkerTask.

Definition at line 49 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 47 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 55 of file TrigPrimTask.cc.

References edm::EventBase::bunchCrossing(), bxBin_, bxBinEdges_, EcalBarrel, EcalEndcap, edm::EventSetup::get(), EcalTPGTowerStatus::getMap(), EcalTPGStripStatus::getMap(), ecaldqm::DQWorker::MEs_, nBXBins, realTps_, EcalTrigTowerDetId::subDet(), EcalElectronicsId::subdet(), 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_ = "";
     //    }
     //  }
     //       }
     //     }
   }

template<typename DigiCollection >

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

Definition at line 159 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 231 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"));
 
     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);
 
       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;
 
             std::vector<int> matchedIndex(0);
             for(int iDigi(0); iDigi < 5; iDigi++){
               if((*tpItr)[iDigi].compressedEt() == realEt)
                 matchedIndex.push_back(iDigi + 1);
             }
 
             if(!matchedIndex.size()) matchedIndex.push_back(0);
             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 168 of file TrigPrimTask.cc.

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

Referenced by analyze().

   {
     MESet& meEtVsBx(MEs_.at("EtVsBx"));
     MESet& meEtReal(MEs_.at("EtReal"));
     MESet& meEtRealMap(MEs_.at("EtRealMap"));
     MESet& meEtSummary(MEs_.at("EtSummary"));
     MESet& meLowIntMap(MEs_.at("LowIntMap"));
     MESet& meMedIntMap(MEs_.at("MedIntMap"));
     MESet& meHighIntMap(MEs_.at("HighIntMap"));
     MESet& meTTFlags(MEs_.at("TTFlags"));
     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);
 
       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;
       }
 
       meTTFlags.fill(ttid, float(tpItr->ttFlag()));
 
       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]);
   }

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

overrideprivatevirtual

Reimplemented from ecaldqm::DQWorker.

Definition at line 36 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"));
     }
   }