#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

	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

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 25 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));
 }