#include <EcalDigiProducer_Ph2.h>

Inheritance diagram for EcalDigiProducer_Ph2:

Public Types
typedef EcalTDigitizer < EBDigitizerTraits_Ph2 >	EBDigitizer_Ph2

typedef EBDigitizerTraits_Ph2::ElectronicsSim	EcalElectronicsSim_Ph2

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

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

	EcalDigiProducer_Ph2 (const edm::ParameterSet &params, edm::ProducesCollector producesCollector, edm::ConsumesCollector &iC)

	EcalDigiProducer_Ph2 (const edm::ParameterSet &params, edm::ConsumesCollector &iC)

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

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

void	setEBNoiseSignalGenerator (EcalBaseSignalGenerator *noiseGenerator)

	~EcalDigiProducer_Ph2 () override

Public Member Functions inherited from DigiAccumulatorMixMod
virtual void	beginRun (edm::Run const &run, edm::EventSetup const &setup)

	DigiAccumulatorMixMod ()

	DigiAccumulatorMixMod (DigiAccumulatorMixMod const &)=delete

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)

DigiAccumulatorMixMod const &	operator= (DigiAccumulatorMixMod const &)=delete

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

virtual	~DigiAccumulatorMixMod ()

Protected Attributes
std::unique_ptr < EBHitResponse_Ph2 >	m_EBResponse

std::unique_ptr< const EcalSimParameterMap >	m_ParameterMap

Private Types
typedef CaloTSamples< float, ecalPh2::sampleSize >	EcalSamples_Ph2

typedef edm::Handle < std::vector< PCaloHit > >	HitsHandle

Private Member Functions
void	accumulateCaloHits (HitsHandle const &ebHandle, int bunchCrossing)

virtual void	cacheEBDigis (const EBDigiCollectionPh2 *ebDigiPtr) const

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

void	checkGeometry (const edm::EventSetup &eventSetup)

void	updateGeometry ()

Private Attributes
edm::ESGetToken < EcalADCToGeVConstant, EcalADCToGeVConstantRcd >	agcToken_

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	geom_token_

edm::ESGetToken < EcalIntercalibConstants, EcalIntercalibConstantsRcd >	icalToken_

edm::ESGetToken < EcalLaserDbService, EcalLaserDbRecord >	laserToken_

std::unique_ptr< EcalLiteDTUCoder >	m_APDCoder

const std::string	m_apdDigiTag

std::unique_ptr< EBDigitizer_Ph2 >	m_APDDigitizer

std::unique_ptr < EcalElectronicsSim < EcalLiteDTUCoder, EcalSamples_Ph2, EcalDataFrame_Ph2 > >	m_APDElectronicsSim

std::unique_ptr< const APDSimParameters >	m_apdParameters

std::unique_ptr < EBHitResponse_Ph2 >	m_APDResponse

const bool	m_apdSeparateDigi

APDShape	m_APDShape

std::unique_ptr< EBDigitizer_Ph2 >	m_BarrelDigitizer

std::unique_ptr< EcalLiteDTUCoder >	m_Coder

std::array< std::unique_ptr < CorrelatedNoisifier < EcalCorrMatrix_Ph2 > >, 2 >	m_EBCorrNoise

const std::string	m_EBdigiCollection

const double	m_EBs25notCont

EBShape	m_EBShape

std::unique_ptr < EcalElectronicsSim_Ph2 >	m_ElectronicsSim

const CaloGeometry *	m_Geometry

const std::string	m_hitsProducerTag

const bool	m_PreMix1

const bool	m_PreMix2

const unsigned int	m_readoutFrameSize

bool	m_useLCcorrection

edm::ESGetToken < EcalLiteDTUPedestalsMap, EcalLiteDTUPedestalsRcd >	pedestalToken_

CLHEP::HepRandomEngine *	randomEngine_ = nullptr

Detailed Description

Definition at line 51 of file EcalDigiProducer_Ph2.h.

Member Typedef Documentation

typedef EcalTDigitizer<EBDigitizerTraits_Ph2> EcalDigiProducer_Ph2::EBDigitizer_Ph2

Definition at line 53 of file EcalDigiProducer_Ph2.h.

typedef EBDigitizerTraits_Ph2::ElectronicsSim EcalDigiProducer_Ph2::EcalElectronicsSim_Ph2

Definition at line 54 of file EcalDigiProducer_Ph2.h.

typedef CaloTSamples<float, ecalPh2::sampleSize> EcalDigiProducer_Ph2::EcalSamples_Ph2

private

Definition at line 124 of file EcalDigiProducer_Ph2.h.

typedef edm::Handle<std::vector<PCaloHit> > EcalDigiProducer_Ph2::HitsHandle

private

Definition at line 73 of file EcalDigiProducer_Ph2.h.

Constructor & Destructor Documentation

EcalDigiProducer_Ph2::EcalDigiProducer_Ph2	(	const edm::ParameterSet &	params,
		edm::ProducesCollector	producesCollector,
		edm::ConsumesCollector &	iC
	)

Definition at line 44 of file EcalDigiProducer_Ph2.cc.

References m_apdDigiTag, m_apdSeparateDigi, m_EBdigiCollection, and edm::ProducesCollector::produces().

     : EcalDigiProducer_Ph2(params, iC) {
   if (m_apdSeparateDigi)
     producesCollector.produces<EBDigiCollectionPh2>(m_apdDigiTag);
 
   producesCollector.produces<EBDigiCollectionPh2>(m_EBdigiCollection);
 }

EcalDigiProducer_Ph2::EcalDigiProducer_Ph2	(	const edm::ParameterSet &	params,
		edm::ConsumesCollector &	iC
	)

Definition at line 55 of file EcalDigiProducer_Ph2.cc.

     : DigiAccumulatorMixMod(),
       m_APDShape(iC),
       m_EBShape(iC),
 
       m_EBdigiCollection(params.getParameter<std::string>("EBdigiCollectionPh2")),
 
       m_hitsProducerTag(params.getParameter<std::string>("hitsProducer")),
       m_useLCcorrection(params.getUntrackedParameter<bool>("UseLCcorrection")),
       m_apdSeparateDigi(params.getParameter<bool>("apdSeparateDigi")),
 
       m_EBs25notCont(params.getParameter<double>("EBs25notContainment")),
 
       m_readoutFrameSize(ecalPh2::sampleSize),
 
       m_ParameterMap(std::make_unique<EcalSimParameterMap>(params.getParameter<double>("simHitToPhotoelectronsBarrel"),
                                                            0,  // endcap parameters not needed
                                                            params.getParameter<double>("photoelectronsToAnalogBarrel"),
                                                            0,
                                                            params.getParameter<double>("samplingFactor"),
                                                            params.getParameter<double>("timePhase"),
                                                            m_readoutFrameSize,
                                                            params.getParameter<int>("binOfMaximum"),
                                                            params.getParameter<bool>("doPhotostatistics"),
                                                            params.getParameter<bool>("syncPhase"))),
 
       m_apdDigiTag(params.getParameter<std::string>("apdDigiTag")),
       m_apdParameters(std::make_unique<APDSimParameters>(params.getParameter<bool>("apdAddToBarrel"),
                                                          m_apdSeparateDigi,
                                                          params.getParameter<double>("apdSimToPELow"),
                                                          params.getParameter<double>("apdSimToPEHigh"),
                                                          params.getParameter<double>("apdTimeOffset"),
                                                          params.getParameter<double>("apdTimeOffWidth"),
                                                          params.getParameter<bool>("apdDoPEStats"),
                                                          m_apdDigiTag,
                                                          params.getParameter<std::vector<double>>("apdNonlParms"))),
 
       m_APDResponse(!m_apdSeparateDigi
                         ? nullptr
                         : std::make_unique<EBHitResponse_Ph2>(
                               m_ParameterMap.get(), &m_EBShape, true, m_apdParameters.get(), &m_APDShape)),
 
       m_EBResponse(std::make_unique<EBHitResponse_Ph2>(m_ParameterMap.get(),
                                                        &m_EBShape,
                                                        false,  // barrel
                                                        m_apdParameters.get(),
                                                        &m_APDShape)),
 
       m_PreMix1(params.getParameter<bool>("EcalPreMixStage1")),
       m_PreMix2(params.getParameter<bool>("EcalPreMixStage2")),
 
       m_APDDigitizer(nullptr),
       m_BarrelDigitizer(nullptr),
       m_ElectronicsSim(nullptr),
       m_Coder(nullptr),
       m_APDElectronicsSim(nullptr),
       m_APDCoder(nullptr),
       m_Geometry(nullptr),
       m_EBCorrNoise({{nullptr, nullptr}})

EcalDigiProducer_Ph2::~EcalDigiProducer_Ph2 ( )

override

Definition at line 174 of file EcalDigiProducer_Ph2.cc.

174 {}

Member Function Documentation

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 199 of file EcalDigiProducer_Ph2.cc.

References accumulateCaloHits(), edm::Event::getByLabel(), m_APDShape, m_EBShape, m_hitsProducerTag, and EcalShapeBase::setEventSetup().

                                                                                         {
   // Step A: Get Inputs
   edm::Handle<std::vector<PCaloHit>> ebHandle;
 
   m_EBShape.setEventSetup(eventSetup);
   m_APDShape.setEventSetup(eventSetup);
   edm::InputTag ebTag(m_hitsProducerTag, "EcalHitsEB");
   e.getByLabel(ebTag, ebHandle);
 
   accumulateCaloHits(ebHandle, 0);
 }

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 211 of file EcalDigiProducer_Ph2.cc.

References accumulateCaloHits(), PileUpEventPrincipal::bunchCrossing(), PileUpEventPrincipal::getByLabel(), and m_hitsProducerTag.

                                                                    {
   // Step A: Get Inputs
   edm::Handle<std::vector<PCaloHit>> ebHandle;
 
   edm::InputTag ebTag(m_hitsProducerTag, "EcalHitsEB");
   e.getByLabel(ebTag, ebHandle);
 
   accumulateCaloHits(ebHandle, e.bunchCrossing());
 }

void EcalDigiProducer_Ph2::accumulateCaloHits	(	HitsHandle const &	ebHandle,
		int	bunchCrossing
	)

private

Definition at line 189 of file EcalDigiProducer_Ph2.cc.

References edm::HandleBase::isValid(), m_APDDigitizer, m_apdSeparateDigi, m_BarrelDigitizer, edm::Handle< T >::product(), and randomEngine_.

Referenced by accumulate().

                                                                                            {
   if (ebHandle.isValid()) {
     m_BarrelDigitizer->add(*ebHandle.product(), bunchCrossing, randomEngine_);
 
     if (m_apdSeparateDigi) {
       m_APDDigitizer->add(*ebHandle.product(), bunchCrossing, randomEngine_);
     }
   }
 }

void EcalDigiProducer_Ph2::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumi,
		edm::EventSetup const &	setup
	)

overridevirtual

Reimplemented from DigiAccumulatorMixMod.

Definition at line 249 of file EcalDigiProducer_Ph2.cc.

References Exception, edm::RandomNumberGenerator::getEngine(), edm::LuminosityBlock::index(), edm::Service< T >::isAvailable(), m_APDResponse, and m_EBResponse.

                                                                                                           {
   edm::Service<edm::RandomNumberGenerator> rng;
   if (!rng.isAvailable()) {
     throw cms::Exception("Configuration") << "RandomNumberGenerator service is not available.\n"
                                              "You must add the service in the configuration file\n"
                                              "or remove the module that requires it.";
   }
   CLHEP::HepRandomEngine* engine = &rng->getEngine(lumi.index());
 
   if (nullptr != m_APDResponse)
     m_APDResponse->initialize(engine);
   m_EBResponse->initialize(engine);
 }

virtual void EcalDigiProducer_Ph2::cacheEBDigis ( const EBDigiCollectionPh2 * ebDigiPtr ) const

inlineprivatevirtual

Definition at line 71 of file EcalDigiProducer_Ph2.h.

Referenced by finalizeEvent().

71 {}

void EcalDigiProducer_Ph2::checkCalibrations	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

private

Definition at line 263 of file EcalDigiProducer_Ph2.cc.

References agcToken_, ecalPh2::gains, edm::EventSetup::getData(), icalToken_, EcalCondDBWriter_cfi::laser, laserToken_, LogDebug, m_APDCoder, m_APDResponse, m_Coder, m_EBResponse, m_EBs25notCont, m_useLCcorrection, ecalPh2::MAXADC, and pedestalToken_.

Referenced by initializeEvent().

                                                                                                    {
   // Pedestals from event setup
   auto pedestals = &eventSetup.getData(pedestalToken_);
 
   m_Coder->setPedestals(pedestals);
   if (nullptr != m_APDCoder)
     m_APDCoder->setPedestals(pedestals);
 
   // Ecal Intercalibration Constants
   auto ical = &eventSetup.getData(icalToken_);
 
   m_Coder->setIntercalibConstants(ical);
   if (nullptr != m_APDCoder)
     m_APDCoder->setIntercalibConstants(ical);
 
   m_EBResponse->setIntercal(ical);
   if (nullptr != m_APDResponse)
     m_APDResponse->setIntercal(ical);
 
   // Ecal LaserCorrection Constants
   auto laser = &eventSetup.getData(laserToken_);
 
   const edm::TimeValue_t eventTimeValue = event.time().value();
 
   m_EBResponse->setEventTime(eventTimeValue);
   m_EBResponse->setLaserConstants(laser, m_useLCcorrection);
 
   // ADC -> GeV Scale
   auto agc = &eventSetup.getData(agcToken_);
 
   m_Coder->setGainRatios(ecalPh2::gains[0] / ecalPh2::gains[1]);
   if (nullptr != m_APDCoder)
     m_APDCoder->setGainRatios(ecalPh2::gains[0] / ecalPh2::gains[1]);
 
   const double EBscale((agc->getEBValue()) * ecalPh2::gains[1] * (ecalPh2::MAXADC)*m_EBs25notCont);
 
   LogDebug("EcalDigi") << " GeV/ADC = " << agc->getEBValue() << "\n"
                        << " notCont = " << m_EBs25notCont << "\n"
                        << " saturation for EB = " << EBscale << ", " << m_EBs25notCont;
 
   m_Coder->setFullScaleEnergy(EBscale);
   if (nullptr != m_APDCoder)
     m_APDCoder->setFullScaleEnergy(EBscale);
 }

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

private

Definition at line 308 of file EcalDigiProducer_Ph2.cc.

References geom_token_, edm::EventSetup::getHandle(), m_Geometry, and updateGeometry().

Referenced by initializeEvent().

                                                                         {
   // TODO find a way to avoid doing this every event
   edm::ESHandle<CaloGeometry> hGeometry = eventSetup.getHandle(geom_token_);
   const CaloGeometry* pGeometry = &*hGeometry;
 
   if (pGeometry != m_Geometry) {
     m_Geometry = pGeometry;
     updateGeometry();
   }
 }

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 223 of file EcalDigiProducer_Ph2.cc.

References cacheEBDigis(), m_APDDigitizer, m_apdSeparateDigi, m_BarrelDigitizer, m_EBdigiCollection, eostools::move(), and randomEngine_.

                                                                                          {
   // Step B: Create empty output
   std::unique_ptr<EBDigiCollectionPh2> apdResult(nullptr);
   std::unique_ptr<EBDigiCollectionPh2> barrelResult = std::make_unique<EBDigiCollectionPh2>();
   if (m_apdSeparateDigi) {
     apdResult = std::make_unique<EBDigiCollectionPh2>();
   }
   // run the algorithm
 
   m_BarrelDigitizer->run(*barrelResult, randomEngine_);
   cacheEBDigis(&*barrelResult);
 
   edm::LogInfo("DigiInfo") << "EB Digis: " << barrelResult->size();
 
   if (m_apdSeparateDigi) {
     m_APDDigitizer->run(*apdResult, randomEngine_);
     edm::LogInfo("DigiInfo") << "APD Digis: " << apdResult->size();
   }
 
   // Step D: Put outputs into event
 
   event.put(std::move(barrelResult), m_EBdigiCollection);
 
   randomEngine_ = nullptr;  // to prevent access outside event
 }

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 176 of file EcalDigiProducer_Ph2.cc.

References checkCalibrations(), checkGeometry(), edm::RandomNumberGenerator::getEngine(), m_APDDigitizer, m_apdSeparateDigi, m_BarrelDigitizer, randomEngine_, and edm::Event::streamID().

                                                                                                  {
   edm::Service<edm::RandomNumberGenerator> rng;
   randomEngine_ = &rng->getEngine(event.streamID());
 
   checkGeometry(eventSetup);
   checkCalibrations(event, eventSetup);
 
   m_BarrelDigitizer->initializeHits();
   if (m_apdSeparateDigi) {
     m_APDDigitizer->initializeHits();
   }
 }

void EcalDigiProducer_Ph2::setEBNoiseSignalGenerator ( EcalBaseSignalGenerator * noiseGenerator )

Definition at line 325 of file EcalDigiProducer_Ph2.cc.

References m_BarrelDigitizer.

                                                                                             {
   m_BarrelDigitizer->setNoiseSignalGenerator(noiseGenerator);
 }

void EcalDigiProducer_Ph2::updateGeometry ( )

private

Definition at line 319 of file EcalDigiProducer_Ph2.cc.

References DetId::Ecal, EcalBarrel, CaloGeometry::getSubdetectorGeometry(), m_APDResponse, m_EBResponse, and m_Geometry.

Referenced by checkGeometry().

                                           {
   if (nullptr != m_APDResponse)
     m_APDResponse->setGeometry(m_Geometry->getSubdetectorGeometry(DetId::Ecal, EcalBarrel));
   m_EBResponse->setGeometry(m_Geometry->getSubdetectorGeometry(DetId::Ecal, EcalBarrel));
 }