#include <CastorDigiProducer.h>

Inheritance diagram for CastorDigiProducer:

Public Member Functions
virtual void	accumulate (edm::Event const &e, edm::EventSetup const &c) override

virtual void	accumulate (PileUpEventPrincipal const &e, edm::EventSetup const &c, edm::StreamID const &) override

	CastorDigiProducer (const edm::ParameterSet &ps, edm::one::EDProducerBase &mixMod, edm::ConsumesCollector &iC)

virtual void	finalizeEvent (edm::Event &e, edm::EventSetup const &c) override

virtual void	initializeEvent (edm::Event const &e, edm::EventSetup const &c) override

virtual	~CastorDigiProducer ()

Public Member Functions inherited from DigiAccumulatorMixMod
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup)

virtual void	beginRun (edm::Run const &run, edm::EventSetup const &setup)

	DigiAccumulatorMixMod ()

virtual void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup)

virtual void	endRun (edm::Run const &run, edm::EventSetup const &setup)

virtual void	finalizeBunchCrossing (edm::Event &event, edm::EventSetup const &setup, int bunchCrossing)

virtual PileupMixingContent *	getEventPileupInfo ()

virtual void	initializeBunchCrossing (edm::Event const &event, edm::EventSetup const &setup, int bunchCrossing)

virtual void	StorePileupInformation (std::vector< int > &numInteractionList, std::vector< int > &bunchCrossingList, std::vector< float > &TrueInteractionList, std::vector< edm::EventID > &eventList, int bunchSpace)

virtual	~DigiAccumulatorMixMod ()

Private Types
typedef CaloTDigitizer < CastorDigitizerTraits >	CastorDigitizer

Private Member Functions
void	accumulateCaloHits (std::vector< PCaloHit > const &, int bunchCrossing, CLHEP::HepRandomEngine *)

void	checkGeometry (const edm::EventSetup &eventSetup)

void	fillFakeHits ()
	some hits in each subdetector, just for testing purposes More...

CLHEP::HepRandomEngine *	randomEngine (edm::StreamID const &streamID)

void	sortHits (const edm::PCaloHitContainer &hits)
	fills the vectors for each subdetector More...

Private Attributes
std::vector < CLHEP::HepRandomEngine * >	randomEngines_

CastorAmplifier *	theAmplifier

CastorDigitizer *	theCastorDigitizer

CastorHitFilter	theCastorHitFilter

std::vector< PCaloHit >	theCastorHits

CaloVShape *	theCastorIntegratedShape

CaloHitResponse *	theCastorResponse

CaloVShape *	theCastorShape

CastorCoderFactory *	theCoderFactory

CastorElectronicsSim *	theElectronicsSim

CastorHitCorrection *	theHitCorrection

CastorSimParameterMap *	theParameterMap

Detailed Description

Definition at line 34 of file CastorDigiProducer.h.

Member Typedef Documentation

typedef CaloTDigitizer<CastorDigitizerTraits> CastorDigiProducer::CastorDigitizer

private

Reconstruction algorithm

Definition at line 59 of file CastorDigiProducer.h.

Constructor & Destructor Documentation

CastorDigiProducer::CastorDigiProducer	(	const edm::ParameterSet &	ps,
		edm::one::EDProducerBase &	mixMod,
		edm::ConsumesCollector &	iC
	)

explicit

Definition at line 22 of file CastorDigiProducer.cc.

References edm::ConsumesCollector::consumes(), CastorCoderFactory::DB, edm::hlt::Exception, edm::ParameterSet::getParameter(), edm::Service< T >::isAvailable(), CaloHitResponse::setHitCorrection(), CaloHitResponse::setHitFilter(), theAmplifier, theCastorDigitizer, theCastorHitFilter, theCastorResponse, theCoderFactory, theElectronicsSim, theHitCorrection, and theParameterMap.

 : theParameterMap(new CastorSimParameterMap(ps)),
   theCastorShape(new CastorShape()),
   theCastorIntegratedShape(new CaloShapeIntegrator(theCastorShape)),
   theCastorResponse(new CaloHitResponse(theParameterMap, theCastorIntegratedShape)),
   theAmplifier(0),
   theCoderFactory(0),
   theElectronicsSim(0),
   theHitCorrection(0),
   theCastorDigitizer(0),
   theCastorHits()
 {
   iC.consumes<std::vector<PCaloHit> >(edm::InputTag("g4SimHits", "CastorFI"));
   
   mixMod.produces<CastorDigiCollection>();
 
   theCastorResponse->setHitFilter(&theCastorHitFilter);
 
   bool doTimeSlew = ps.getParameter<bool>("doTimeSlew");
   if(doTimeSlew) {
     // no time slewing for HF
     theCastorResponse->setHitCorrection(theHitCorrection);
   }
 
   bool doNoise = ps.getParameter<bool>("doNoise");
   theAmplifier = new CastorAmplifier(theParameterMap, doNoise);
   theCoderFactory = new CastorCoderFactory(CastorCoderFactory::DB);
   theElectronicsSim = new CastorElectronicsSim(theAmplifier, theCoderFactory);
 
   theCastorDigitizer = new CastorDigitizer(theCastorResponse, theElectronicsSim, doNoise);
 
   edm::Service<edm::RandomNumberGenerator> rng;
   if ( ! rng.isAvailable()) {
     throw cms::Exception("Configuration")
       << "CastorDigiProducer requires the RandomNumberGeneratorService\n"
          "which is not present in the configuration file.  You must add the service\n"
          "in the configuration file or remove the modules that require it.";
   }
 }

CastorDigiProducer::~CastorDigiProducer ( )

virtual

Definition at line 63 of file CastorDigiProducer.cc.

References theAmplifier, theCastorDigitizer, theCastorIntegratedShape, theCastorResponse, theCastorShape, theCoderFactory, theElectronicsSim, theHitCorrection, and theParameterMap.

                                         {
   delete theCastorDigitizer;
   delete theParameterMap;
   delete theCastorShape;
   delete theCastorIntegratedShape;
   delete theCastorResponse;
   delete theElectronicsSim;
   delete theAmplifier;
   delete theCoderFactory;
   delete theHitCorrection;
 }

Member Function Documentation

void CastorDigiProducer::accumulate	(	edm::Event const &	e,
		edm::EventSetup const &	c
	)

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 102 of file CastorDigiProducer.cc.

References accumulateCaloHits(), edm::Event::getByLabel(), edm::Handle< T >::product(), randomEngine(), and edm::Event::streamID().

                                                                            {
   // Step A: Get and accumulate digitized hits 
   edm::Handle<std::vector<PCaloHit> > castorHandle;
   e.getByLabel(edm::InputTag("g4SimHits", "CastorFI"), castorHandle);
 
   accumulateCaloHits(*castorHandle.product(), 0, randomEngine(e.streamID()));
 }

void CastorDigiProducer::accumulate	(	PileUpEventPrincipal const &	e,
		edm::EventSetup const &	c,
		edm::StreamID const &	streamID
	)

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 110 of file CastorDigiProducer.cc.

References accumulateCaloHits(), PileUpEventPrincipal::bunchCrossing(), PileUpEventPrincipal::getByLabel(), edm::Handle< T >::product(), and randomEngine().

                                                                                                                     {
   // Step A: Get and accumulate digitized hits 
   edm::Handle<std::vector<PCaloHit> > castorHandle;
   e.getByLabel(edm::InputTag("g4SimHits", "CastorFI"), castorHandle);
 
   accumulateCaloHits(*castorHandle.product(), e.bunchCrossing(), randomEngine(streamID));
 }

void CastorDigiProducer::accumulateCaloHits	(	std::vector< PCaloHit > const &	hcalHits,
		int	bunchCrossing,
		CLHEP::HepRandomEngine *	engine
	)

private

Definition at line 93 of file CastorDigiProducer.cc.

References CaloTDigitizer< Traits >::add(), CastorHitCorrection::fillChargeSums(), theCastorDigitizer, and theHitCorrection.

Referenced by accumulate().

                                                                                                                                 {
   //fillFakeHits();
 
   if(theHitCorrection != 0) {
     theHitCorrection->fillChargeSums(hcalHits);
   }
   theCastorDigitizer->add(hcalHits, bunchCrossing, engine);
 }

void CastorDigiProducer::checkGeometry ( const edm::EventSetup & eventSetup )

private

make sure the digitizer has the correct list of all cells that exist in the geometry

Definition at line 153 of file CastorDigiProducer.cc.

References DetId::Calo, geometry, edm::EventSetup::get(), CaloTDigitizer< Traits >::setDetIds(), CaloHitResponse::setGeometry(), HcalCastorDetId::SubdetectorId, theCastorDigitizer, and theCastorResponse.

Referenced by initializeEvent().

                                                                        {
   // TODO find a way to avoid doing this every event
   edm::ESHandle<CaloGeometry> geometry;
   eventSetup.get<CaloGeometryRecord>().get(geometry);
   theCastorResponse->setGeometry(&*geometry);
 
   const std::vector<DetId>& castorCells = geometry->getValidDetIds(DetId::Calo, HcalCastorDetId::SubdetectorId);
 
   //std::cout<<"CastorDigiProducer::CheckGeometry number of cells: "<<castorCells.size()<<std::endl;
   theCastorDigitizer->setDetIds(castorCells);
 }

void CastorDigiProducer::fillFakeHits ( )

private

some hits in each subdetector, just for testing purposes

Definition at line 146 of file CastorDigiProducer.cc.

References DetId::rawId(), and theCastorHits.

                                       {
   HcalCastorDetId castorDetId(HcalCastorDetId::Section(2),true,1,1);
 
   theCastorHits.emplace_back(castorDetId.rawId(), 50.0, 0.);
 }

void CastorDigiProducer::finalizeEvent	(	edm::Event &	e,
		edm::EventSetup const &	c
	)

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 118 of file CastorDigiProducer.cc.

References edm::Event::put(), randomEngine(), CaloTDigitizer< Traits >::run(), edm::Event::streamID(), and theCastorDigitizer.

                                                                                    {
   // Step B: Create empty output
 
   std::auto_ptr<CastorDigiCollection> castorResult(new CastorDigiCollection());
 
   // Step C: Invoke the algorithm, getting back outputs.
   theCastorDigitizer->run(*castorResult, randomEngine(e.streamID()));
 
   edm::LogInfo("CastorDigiProducer") << "HCAL/Castor digis   : " << castorResult->size();
 
   // Step D: Put outputs into event
   e.put(castorResult);
 }

void CastorDigiProducer::initializeEvent	(	edm::Event const &	e,
		edm::EventSetup const &	c
	)

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 75 of file CastorDigiProducer.cc.

References checkGeometry(), edm::EventSetup::get(), CaloTDigitizer< Traits >::initializeHits(), edm::ESHandle< class >::product(), CastorCoderFactory::setDbService(), CastorAmplifier::setDbService(), CastorSimParameterMap::setDbService(), theAmplifier, theCastorDigitizer, theCastorHits, theCoderFactory, and theParameterMap.

                                                                                          {
   // get the appropriate gains, noises, & widths for this event
   edm::ESHandle<CastorDbService> conditions;
   eventSetup.get<CastorDbRecord>().get(conditions);
   theAmplifier->setDbService(conditions.product());
   theCoderFactory->setDbService(conditions.product());
   theParameterMap->setDbService(conditions.product());
 
   edm::LogInfo("CastorDigiProducer") << "checking the geometry...";
 
   // get the correct geometry
   checkGeometry(eventSetup);
   
   theCastorHits.clear();
 
   theCastorDigitizer->initializeHits();
 }

CLHEP::HepRandomEngine * CastorDigiProducer::randomEngine ( edm::StreamID const & streamID )

private

Definition at line 165 of file CastorDigiProducer.cc.

References edm::RandomNumberGenerator::getEngine(), cmsHarvester::index, randomEngines_, and edm::StreamID::value().

Referenced by accumulate(), and finalizeEvent().

                                                                                   {
   unsigned int index = streamID.value();
   if(index >= randomEngines_.size()) {
     randomEngines_.resize(index + 1, nullptr);
   }
   CLHEP::HepRandomEngine* ptr = randomEngines_[index];
   if(!ptr) {
     edm::Service<edm::RandomNumberGenerator> rng;
     ptr = &rng->getEngine(streamID);
     randomEngines_[index] = ptr;
   }
   return ptr;
 }

void CastorDigiProducer::sortHits ( const edm::PCaloHitContainer & hits )

private

fills the vectors for each subdetector

Definition at line 133 of file CastorDigiProducer.cc.

References DetId::Calo, DetId::det(), HcalCastorDetId::SubdetectorId, DetId::subdetId(), and theCastorHits.

                                                                   {
   for(edm::PCaloHitContainer::const_iterator hitItr = hits.begin();
       hitItr != hits.end(); ++hitItr){
     DetId detId = hitItr->id();
     if (detId.det()==DetId::Calo && detId.subdetId()==HcalCastorDetId::SubdetectorId){
       theCastorHits.push_back(*hitItr);
     }
       else {
     edm::LogError("CastorDigiProducer") << "Bad Hit subdetector " << detId.subdetId();
       }
   }
 }