Inheritance diagram for TBHodoActiveVolumeRawInfoProducer:

Public Member Functions
void	produce (edm::Event &event, const edm::EventSetup &eventSetup) override
	Produce digis out of raw data. More...

	TBHodoActiveVolumeRawInfoProducer (const edm::ParameterSet &ps)
	Constructor. More...

	~TBHodoActiveVolumeRawInfoProducer () override=default
	Destructor. More...

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 Attributes
edm::EDGetTokenT < edm::PCaloHitContainer >	m_EcalToken

double	myThreshold

std::unique_ptr < EcalTBHodoscopeGeometry >	theTBHodoGeom_

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

Definition at line 27 of file TBHodoActiveVolumeRawInfoProducer.cc.

Constructor & Destructor Documentation

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

explicit

Constructor.

Definition at line 48 of file TBHodoActiveVolumeRawInfoProducer.cc.

References theTBHodoGeom_.

     : myThreshold(0.05E-3),
       m_EcalToken(consumes<edm::PCaloHitContainer>(edm::InputTag("g4SimHits", "EcalTBH4BeamHits"))) {
   produces<EcalTBHodoscopeRawInfo>();
 
   theTBHodoGeom_ = std::make_unique<EcalTBHodoscopeGeometry>();
 }

TBHodoActiveVolumeRawInfoProducer::~TBHodoActiveVolumeRawInfoProducer ( )

overridedefault

Destructor.

Member Function Documentation

void TBHodoActiveVolumeRawInfoProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

override

Produce digis out of raw data.

Definition at line 56 of file TBHodoActiveVolumeRawInfoProducer.cc.

References HodoscopeDetId::fibrId(), edm::HandleBase::isValid(), LogDebug, m_EcalToken, eostools::move(), myThreshold, HodoscopeDetId::planeId(), edm::Handle< T >::product(), EcalTBHodoscopePlaneRawHits::setHit(), and theTBHodoGeom_.

                                                                                                 {
   std::unique_ptr<EcalTBHodoscopeRawInfo> product(new EcalTBHodoscopeRawInfo());
 
   // caloHit container
   const edm::Handle<edm::PCaloHitContainer> &pCaloHit = event.getHandle(m_EcalToken);
   const edm::PCaloHitContainer *caloHits = nullptr;
   if (pCaloHit.isValid()) {
     caloHits = pCaloHit.product();
     LogDebug("EcalTBHodo") << "total # caloHits: " << caloHits->size();
   } else {
     edm::LogError("EcalTBHodo") << "Error! can't get the caloHitContainer ";
   }
   if (!caloHits) {
     return;
   }
 
   // detid - energy_sum map
   std::map<unsigned int, double> energyMap;
 
   for (auto &&aHit : *caloHits) {
     double thisHitEne = aHit.energy();
 
     std::map<unsigned int, double>::iterator itmap = energyMap.find(aHit.id());
     if (itmap == energyMap.end())
       energyMap.insert(pair<unsigned int, double>(aHit.id(), thisHitEne));
     else {
       (*itmap).second += thisHitEne;
     }
   }
 
   // planes and fibers
   int nPlanes = theTBHodoGeom_->getNPlanes();
   int nFibers = theTBHodoGeom_->getNFibres();
   product->setPlanes(nPlanes);
 
   bool firedChannels[4][64];
   for (int iPlane = 0; iPlane < nPlanes; ++iPlane) {
     for (int iFiber = 0; iFiber < nFibers; ++iFiber) {
       firedChannels[iPlane][iFiber] = 0.;
     }
   }
   for (std::map<unsigned int, double>::const_iterator itmap = energyMap.begin(); itmap != energyMap.end(); ++itmap) {
     if ((*itmap).second > myThreshold) {
       HodoscopeDetId myHodoDetId = HodoscopeDetId((*itmap).first);
       firedChannels[myHodoDetId.planeId()][myHodoDetId.fibrId()] = true;
     }
   }
   for (int iPlane = 0; iPlane < nPlanes; ++iPlane) {
     EcalTBHodoscopePlaneRawHits planeHit(nFibers);
 
     for (int iFiber = 0; iFiber < nFibers; ++iFiber) {
       planeHit.setHit(iFiber, firedChannels[iPlane][iFiber]);
     }
     product->setPlane((unsigned int)iPlane, planeHit);
   }
 
   LogDebug("EcalTBHodo") << (*product);
 
   event.put(std::move(product));
 }