TBHodoActiveVolumeRawInfoProducer Class Reference

#include <TBHodoActiveVolumeRawInfoProducer.h>

Inheritance diagram for TBHodoActiveVolumeRawInfoProducer:

Public Member Functions
void	produce (edm::Event &event, const edm::EventSetup &eventSetup)
	Produce digis out of raw data. More...
	TBHodoActiveVolumeRawInfoProducer (const edm::ParameterSet &ps)
	Constructor. More...
virtual	~TBHodoActiveVolumeRawInfoProducer ()
	Destructor. More...
Public Member Functions inherited from edm::EDProducer
	EDProducer ()
ModuleDescription const &	moduleDescription () const
virtual	~EDProducer ()
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const
std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Attributes
double	myThreshold
EcalTBHodoscopeGeometry *	theTBHodoGeom_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 27 of file TBHodoActiveVolumeRawInfoProducer.h.

Constructor & Destructor Documentation

TBHodoActiveVolumeRawInfoProducer::TBHodoActiveVolumeRawInfoProducer

(

const edm::ParameterSet &

ps

)

Constructor.

Definition at line 14 of file TBHodoActiveVolumeRawInfoProducer.cc.

                                                                                              {

  produces<EcalTBHodoscopeRawInfo>();

  theTBHodoGeom_ = new EcalTBHodoscopeGeometry();

  myThreshold = 0.05E-3;
}

TBHodoActiveVolumeRawInfoProducer::~TBHodoActiveVolumeRawInfoProducer ( )

virtual

Destructor.

Definition at line 23 of file TBHodoActiveVolumeRawInfoProducer.cc.

                                                                      {   
  delete theTBHodoGeom_; 
}

Member Function Documentation

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

virtual

Produce digis out of raw data.

Implements edm::EDProducer.

Definition at line 27 of file TBHodoActiveVolumeRawInfoProducer.cc.

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

{
  auto_ptr<EcalTBHodoscopeRawInfo> product(new EcalTBHodoscopeRawInfo());

  // caloHit container
  edm::Handle<edm::PCaloHitContainer> pCaloHit;
  const edm::PCaloHitContainer* caloHits =0;
  event.getByLabel( "g4SimHits", "EcalTBH4BeamHits", 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;  

  int myCount = 0;
  for(edm::PCaloHitContainer::const_iterator itch = caloHits->begin(); itch != caloHits->end(); ++itch) {
    
    double thisHitEne = itch->energy();

    std::map<unsigned int,double>::iterator itmap = energyMap.find(itch->id());
    if ( itmap == energyMap.end() )
      energyMap.insert(pair<unsigned int, double>( itch->id(), thisHitEne));  
    else{
      (*itmap).second+=thisHitEne;
    }

    myCount++;
  }

  
  // 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()] = 1;
    }

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

Member Data Documentation

double TBHodoActiveVolumeRawInfoProducer::myThreshold

private

Definition at line 43 of file TBHodoActiveVolumeRawInfoProducer.h.

EcalTBHodoscopeGeometry* TBHodoActiveVolumeRawInfoProducer::theTBHodoGeom_

private

Definition at line 45 of file TBHodoActiveVolumeRawInfoProducer.h.