#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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Member Functions
bool	alreadyInserted (EcalTrigTowerDetId ttId)

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

Private Attributes
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_

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

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

produce a subset of TP information

Author: Federico Ferri, CEA/Saclay Irfu/SPP

Definition at line 22 of file EcalTPSkimmer.h.

Constructor & Destructor Documentation

◆ EcalTPSkimmer()

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

explicit

Definition at line 24 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"));
  
   produces<EcalTrigPrimDigiCollection>(tpOutputCollection_);
 }

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

◆ ~EcalTPSkimmer()

EcalTPSkimmer::~EcalTPSkimmer ( )

override

Definition at line 38 of file EcalTPSkimmer.cc.

38 {}

Member Function Documentation

◆ alreadyInserted()

bool EcalTPSkimmer::alreadyInserted ( EcalTrigTowerDetId ttId )

private

Definition at line 117 of file EcalTPSkimmer.cc.

117 { 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 119 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 40 of file EcalTPSkimmer.cc.

                                                                   {
   insertedTP_.clear();
  
   using namespace edm;
  
   es.get<IdealGeometryRecord>().get(ttMap_);
  
   // 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.get<EcalChannelStatusRcd>().get(chStatus);
  
   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_);
 }