#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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () 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));
 }