#include <TBHodoActiveVolumeRawInfoProducer.h>

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
	Destructor. More...

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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Private Attributes
edm::EDGetTokenT< edm::PCaloHitContainer >	m_EcalToken

double	myThreshold

EcalTBHodoscopeGeometry *	theTBHodoGeom_

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

Definition at line 26 of file TBHodoActiveVolumeRawInfoProducer.h.

Constructor & Destructor Documentation

◆ TBHodoActiveVolumeRawInfoProducer()

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

explicit

Constructor.

Definition at line 14 of file TBHodoActiveVolumeRawInfoProducer.cc.

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

◆ ~TBHodoActiveVolumeRawInfoProducer()

TBHodoActiveVolumeRawInfoProducer::~TBHodoActiveVolumeRawInfoProducer ( )

override

Destructor.

Definition at line 23 of file TBHodoActiveVolumeRawInfoProducer.cc.

23 { delete theTBHodoGeom_; }

Member Function Documentation

◆ produce()

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

override

Produce digis out of raw data.

Definition at line 25 of file TBHodoActiveVolumeRawInfoProducer.cc.

                                                                                                 {
   unique_ptr<EcalTBHodoscopeRawInfo> product(new EcalTBHodoscopeRawInfo());
  
   // caloHit container
   edm::Handle<edm::PCaloHitContainer> pCaloHit;
   const edm::PCaloHitContainer *caloHits = nullptr;
   event.getByToken(m_EcalToken, pCaloHit);
   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));
 }