#include <EcalTPSkimmer.h>

Inheritance diagram for EcalTPSkimmer:

Public Member Functions
	EcalTPSkimmer (const edm::ParameterSet &ps)

void	produce (edm::Event &evt, const edm::EventSetup &es) override

	~EcalTPSkimmer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
bool	alreadyInserted (EcalTrigTowerDetId ttId)

void	insertTP (EcalTrigTowerDetId ttId, edm::Handle< EcalTrigPrimDigiCollection > &in, EcalTrigPrimDigiCollection &out)

Private Attributes
edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd >	chStatusToken_

std::vector< uint32_t >	chStatusToSelectTP_

bool	doBarrel_

bool	doEndcap_

std::set< EcalTrigTowerDetId >	insertedTP_

bool	skipModule_

std::string	tpCollection_

edm::EDGetTokenT< EcalTrigPrimDigiCollection >	tpInputToken_

std::string	tpOutputCollection_

edm::ESHandle< EcalTrigTowerConstituentsMap >	ttMap_

edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	ttMapToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

produce a subset of TP information

Author: Federico Ferri, CEA/Saclay Irfu/SPP

Definition at line 26 of file EcalTPSkimmer.h.

Constructor & Destructor Documentation

◆ EcalTPSkimmer()

EcalTPSkimmer::EcalTPSkimmer ( const edm::ParameterSet & ps )

explicit

Definition at line 19 of file EcalTPSkimmer.cc.

                                                       {
   skipModule_ = ps.getParameter<bool>("skipModule");
  
   doBarrel_ = ps.getParameter<bool>("doBarrel");
   doEndcap_ = ps.getParameter<bool>("doEndcap");
  
   chStatusToSelectTP_ = ps.getParameter<std::vector<uint32_t> >("chStatusToSelectTP");
  
   tpOutputCollection_ = ps.getParameter<std::string>("tpOutputCollection");
   tpInputToken_ = consumes<EcalTrigPrimDigiCollection>(ps.getParameter<edm::InputTag>("tpInputCollection"));
   ttMapToken_ = esConsumes<EcalTrigTowerConstituentsMap, IdealGeometryRecord>();
   if (not skipModule_) {
     chStatusToken_ = esConsumes<EcalChannelStatus, EcalChannelStatusRcd>();
   }
   produces<EcalTrigPrimDigiCollection>(tpOutputCollection_);
 }

References chStatusToken_, chStatusToSelectTP_, doBarrel_, doEndcap_, edm::ParameterSet::getParameter(), skipModule_, AlCaHLTBitMon_QueryRunRegistry::string, tpInputToken_, tpOutputCollection_, and ttMapToken_.

◆ ~EcalTPSkimmer()

EcalTPSkimmer::~EcalTPSkimmer ( )

override

Definition at line 36 of file EcalTPSkimmer.cc.

36 {}

Member Function Documentation

◆ alreadyInserted()

bool EcalTPSkimmer::alreadyInserted ( EcalTrigTowerDetId ttId )

private

Definition at line 114 of file EcalTPSkimmer.cc.

114 { return (insertedTP_.find(ttId) != insertedTP_.end()); }

References insertedTP_, and ecaldqm::ttId().

Referenced by produce().

◆ insertTP()

void EcalTPSkimmer::insertTP	(	EcalTrigTowerDetId	ttId,
		edm::Handle< EcalTrigPrimDigiCollection > &	in,
		EcalTrigPrimDigiCollection &	out
	)

private

Definition at line 116 of file EcalTPSkimmer.cc.

                                                                 {
   EcalTrigPrimDigiCollection::const_iterator tpIt = tpIn->find(ttId);
   if (tpIt != tpIn->end()) {
     tpOut.push_back(*tpIt);
     insertedTP_.insert(ttId);
   }
 }

References edm::SortedCollection< T, SORT >::end(), edm::SortedCollection< T, SORT >::find(), insertedTP_, edm::SortedCollection< T, SORT >::push_back(), and ecaldqm::ttId().

Referenced by produce().

◆ produce()

void EcalTPSkimmer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 38 of file EcalTPSkimmer.cc.

                                                                   {
   insertedTP_.clear();
  
   using namespace edm;
  
   ttMap_ = es.getHandle(ttMapToken_);
  
   // collection of rechits to put in the event
   auto tpOut = std::make_unique<EcalTrigPrimDigiCollection>();
  
   if (skipModule_) {
     evt.put(std::move(tpOut), tpOutputCollection_);
     return;
   }
  
   edm::ESHandle<EcalChannelStatus> chStatus = es.getHandle(chStatusToken_);
  
   edm::Handle<EcalTrigPrimDigiCollection> tpIn;
   evt.getByToken(tpInputToken_, tpIn);
  
   if (doBarrel_) {
     EcalChannelStatusMap::const_iterator chit;
     uint16_t code = 0;
     for (int i = 0; i < EBDetId::kSizeForDenseIndexing; ++i) {
       if (!EBDetId::validDenseIndex(i))
         continue;
       EBDetId id = EBDetId::detIdFromDenseIndex(i);
       chit = chStatus->find(id);
       // check if the channel status means TP to be kept
       if (chit != chStatus->end()) {
         code = (*chit).getStatusCode();
         if (std::find(chStatusToSelectTP_.begin(), chStatusToSelectTP_.end(), code) != chStatusToSelectTP_.end()) {
           // retrieve the TP DetId
           EcalTrigTowerDetId ttDetId(((EBDetId)id).tower());
           // insert the TP if not done already
           if (!alreadyInserted(ttDetId))
             insertTP(ttDetId, tpIn, *tpOut);
         }
       } else {
         edm::LogError("EcalDetIdToBeRecoveredProducer") << "No channel status found for xtal " << id.rawId()
                                                         << "! something wrong with EcalChannelStatus in your DB? ";
       }
     }
   }
  
   if (doEndcap_) {
     EcalChannelStatusMap::const_iterator chit;
     uint16_t code = 0;
     for (int i = 0; i < EEDetId::kSizeForDenseIndexing; ++i) {
       if (!EEDetId::validDenseIndex(i))
         continue;
       EEDetId id = EEDetId::detIdFromDenseIndex(i);
       chit = chStatus->find(id);
       // check if the channel status means TP to be kept
       if (chit != chStatus->end()) {
         code = (*chit).getStatusCode();
         if (std::find(chStatusToSelectTP_.begin(), chStatusToSelectTP_.end(), code) != chStatusToSelectTP_.end()) {
           // retrieve the TP DetId
           EcalTrigTowerDetId ttDetId = ttMap_->towerOf(id);
           // insert the TP if not done already
           if (!alreadyInserted(ttDetId))
             insertTP(ttDetId, tpIn, *tpOut);
         }
       } else {
         edm::LogError("EcalDetIdToBeRecoveredProducer") << "No channel status found for xtal " << id.rawId()
                                                         << "! something wrong with EcalChannelStatus in your DB? ";
       }
     }
   }
  
   // put the collection of reconstructed hits in the event
   LogInfo("EcalTPSkimmer") << "total # of TP inserted: " << tpOut->size();
  
   evt.put(std::move(tpOut), tpOutputCollection_);
 }