TBHodoActiveVolumeRawInfoProducer Class Reference

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 26 of file TBHodoActiveVolumeRawInfoProducer.cc.

Constructor & Destructor Documentation

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

explicit

Constructor.

Definition at line 47 of file TBHodoActiveVolumeRawInfoProducer.cc.

References HLT_FULL_cff::InputTag.

                                                                                              {
  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 56 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 58 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.cc.

double TBHodoActiveVolumeRawInfoProducer::myThreshold

private

Definition at line 38 of file TBHodoActiveVolumeRawInfoProducer.cc.

EcalTBHodoscopeGeometry* TBHodoActiveVolumeRawInfoProducer::theTBHodoGeom_

private

Definition at line 41 of file TBHodoActiveVolumeRawInfoProducer.cc.