TBHodoActiveVolumeRawInfoProducer Class Reference

#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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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 ( 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 ( )

override

Destructor.

Definition at line 23 of file TBHodoActiveVolumeRawInfoProducer.cc.

{ delete theTBHodoGeom_; }

Member Function Documentation

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.

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

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

Member Data Documentation

edm::EDGetTokenT<edm::PCaloHitContainer> TBHodoActiveVolumeRawInfoProducer::m_EcalToken

private

Definition at line 40 of file TBHodoActiveVolumeRawInfoProducer.h.

double TBHodoActiveVolumeRawInfoProducer::myThreshold

private

Definition at line 38 of file TBHodoActiveVolumeRawInfoProducer.h.

EcalTBHodoscopeGeometry* TBHodoActiveVolumeRawInfoProducer::theTBHodoGeom_

private

Definition at line 41 of file TBHodoActiveVolumeRawInfoProducer.h.