EcalDigiProducer Class Reference

#include <EcalDigiProducer.h>

Inheritance diagram for EcalDigiProducer:

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
virtual void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override
	EcalDigiProducer (const edm::ParameterSet &params, edm::stream::EDProducerBase &mixMod, edm::ConsumesCollector &iC)
	EcalDigiProducer (const edm::ParameterSet &params, 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
void	setEBNoiseSignalGenerator (EcalBaseSignalGenerator *noiseGenerator)
void	setEENoiseSignalGenerator (EcalBaseSignalGenerator *noiseGenerator)
void	setESNoiseSignalGenerator (EcalBaseSignalGenerator *noiseGenerator)
virtual	~EcalDigiProducer ()
Public Member Functions inherited from DigiAccumulatorMixMod
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 ()

Protected Attributes
std::unique_ptr< EBHitResponse >	m_EBResponse
std::unique_ptr< EEHitResponse >	m_EEResponse
std::unique_ptr< const EcalSimParameterMap >	m_ParameterMap

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

Private Member Functions
void	accumulateCaloHits (HitsHandle const &ebHandle, HitsHandle const &eeHandle, HitsHandle const &esHandle, int bunchCrossing, CLHEP::HepRandomEngine *)
virtual void	cacheEBDigis (const EBDigiCollection *ebDigiPtr) const
virtual void	cacheEEDigis (const EEDigiCollection *eeDigiPtr) const
void	checkCalibrations (const edm::Event &event, const edm::EventSetup &eventSetup)
void	checkGeometry (const edm::EventSetup &eventSetup)
CLHEP::HepRandomEngine *	randomEngine (edm::StreamID const &streamID)
void	updateGeometry ()

Private Attributes
const bool	m_addESNoise
std::unique_ptr< EcalCoder >	m_APDCoder
const std::string	m_apdDigiTag
std::unique_ptr< EBDigitizer >	m_APDDigitizer
std::unique_ptr < EcalElectronicsSim >	m_APDElectronicsSim
std::unique_ptr< const APDSimParameters >	m_apdParameters
std::unique_ptr< EBHitResponse >	m_APDResponse
const bool	m_apdSeparateDigi
const APDShape	m_APDShape
std::unique_ptr< EBDigitizer >	m_BarrelDigitizer
std::unique_ptr< EcalCoder >	m_Coder
const bool	m_doEB
const bool	m_doEE
const bool	m_doES
const bool	m_doFastES
std::array< std::unique_ptr < CorrelatedNoisifier < EcalCorrMatrix > >, 3 >	m_EBCorrNoise
const std::string	m_EBdigiCollection
const double	m_EBs25notCont
const EBShape	m_EBShape
std::array< std::unique_ptr < CorrelatedNoisifier < EcalCorrMatrix > >, 3 >	m_EECorrNoise
const std::string	m_EEdigiCollection
const double	m_EEs25notCont
const EEShape	m_EEShape
std::unique_ptr < EcalElectronicsSim >	m_ElectronicsSim
std::unique_ptr< EEDigitizer >	m_EndcapDigitizer
const std::string	m_ESdigiCollection
std::unique_ptr< ESDigitizer >	m_ESDigitizer
std::unique_ptr< ESElectronicsSim >	m_ESElectronicsSim
std::unique_ptr < ESElectronicsSimFast >	m_ESElectronicsSimFast
std::unique_ptr< ESOldDigitizer >	m_ESOldDigitizer
std::unique_ptr< CaloHitResponse >	m_ESOldResponse
std::unique_ptr< ESHitResponse >	m_ESResponse
ESShape	m_ESShape
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
std::vector < CLHEP::HepRandomEngine * >	randomEngines_

Detailed Description

Definition at line 60 of file EcalDigiProducer.h.

Member Typedef Documentation

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

private

Definition at line 82 of file EcalDigiProducer.h.

Constructor & Destructor Documentation

EcalDigiProducer::EcalDigiProducer	(	const edm::ParameterSet &	params,
		edm::stream::EDProducerBase &	mixMod,
		edm::ConsumesCollector &	iC
	)

Definition at line 57 of file EcalDigiProducer.cc.

                                                                                                                            :
   DigiAccumulatorMixMod(),
   m_APDShape         ( params.getParameter<double>( "apdShapeTstart" ) ,
            params.getParameter<double>( "apdShapeTau"    )   )  ,
   m_EBShape          (   ) ,
   m_EEShape          (   ) ,
   m_ESShape          (   ) ,
   m_EBdigiCollection ( params.getParameter<std::string>("EBdigiCollection") ) ,
   m_EEdigiCollection ( params.getParameter<std::string>("EEdigiCollection") ) ,
   m_ESdigiCollection ( params.getParameter<std::string>("ESdigiCollection") ) ,
   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_EEs25notCont     ( params.getParameter<double>     ("EEs25notContainment") ) ,

   m_readoutFrameSize ( params.getParameter<int>       ("readoutFrameSize") ) ,
   m_ParameterMap     ( new EcalSimParameterMap(
               params.getParameter<double> ("simHitToPhotoelectronsBarrel") ,
               params.getParameter<double> ("simHitToPhotoelectronsEndcap") , 
               params.getParameter<double> ("photoelectronsToAnalogBarrel") ,
               params.getParameter<double> ("photoelectronsToAnalogEndcap") , 
               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 ( new 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 :
           new EBHitResponse( m_ParameterMap.get()  ,
                      &m_EBShape      ,
                      true            ,
                      m_apdParameters.get() ,
                      &m_APDShape       ) ) ,
   
   m_EBResponse ( new EBHitResponse( m_ParameterMap.get()  ,
                     &m_EBShape      ,
                     false           , // barrel
                     m_apdParameters.get() ,
                     &m_APDShape       ) ) ,

   m_EEResponse ( new EEHitResponse( m_ParameterMap.get(),
                     &m_EEShape       ) ) ,
   m_ESResponse ( new ESHitResponse( m_ParameterMap.get(), &m_ESShape ) ) ,
   m_ESOldResponse ( new CaloHitResponse( m_ParameterMap.get(), &m_ESShape ) ) ,

   m_addESNoise           ( params.getParameter<bool> ("doESNoise") ) ,
   m_PreMix1              ( params.getParameter<bool> ("EcalPreMixStage1") ) ,
   m_PreMix2              ( params.getParameter<bool> ("EcalPreMixStage2") ) ,

   m_doFastES             ( params.getParameter<bool> ("doFast"   ) ) ,

   m_doEB                 ( params.getParameter<bool> ("doEB"     ) ) ,
   m_doEE                 ( params.getParameter<bool> ("doEE"     ) ) ,
   m_doES                 ( params.getParameter<bool> ("doES"     ) ) ,

   m_ESElectronicsSim     ( m_doFastES ? nullptr :
                new ESElectronicsSim( m_addESNoise ) ) ,
     
   m_ESOldDigitizer       ( m_doFastES ? nullptr :
                new ESOldDigitizer( m_ESOldResponse.get()    , 
                        m_ESElectronicsSim.get() ,
                        m_addESNoise         ) ) ,
   
   m_ESElectronicsSimFast ( !m_doFastES  ? nullptr :
                new ESElectronicsSimFast( m_addESNoise ,
                                                      m_PreMix1) ) ,

   m_ESDigitizer          ( !m_doFastES ? nullptr :
                new ESDigitizer( m_ESResponse.get()           ,
                         m_ESElectronicsSimFast.get() ,
                         m_addESNoise            ) ) ,

   m_APDDigitizer      ( nullptr ) ,
   m_BarrelDigitizer   ( nullptr ) ,
   m_EndcapDigitizer   ( nullptr ) ,
   m_ElectronicsSim    ( nullptr ) ,
   m_Coder             ( nullptr ) ,
   m_APDElectronicsSim ( nullptr ) ,
   m_APDCoder          ( nullptr ) ,
   m_Geometry          ( nullptr ) ,
   m_EBCorrNoise       ( { {nullptr, nullptr, nullptr} } ) ,

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

Definition at line 281 of file EcalDigiProducer.cc.

                                                                                            :
   DigiAccumulatorMixMod(),
   m_APDShape         ( params.getParameter<double>( "apdShapeTstart" ) ,
            params.getParameter<double>( "apdShapeTau"    )   )  ,
   m_EBShape          (   ) ,
   m_EEShape          (   ) ,
   m_ESShape          (   ) ,
   m_EBdigiCollection ( params.getParameter<std::string>("EBdigiCollection") ) ,
   m_EEdigiCollection ( params.getParameter<std::string>("EEdigiCollection") ) ,
   m_ESdigiCollection ( params.getParameter<std::string>("ESdigiCollection") ) ,
   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_EEs25notCont     ( params.getParameter<double>     ("EEs25notContainment") ) ,

   m_readoutFrameSize ( params.getParameter<int>       ("readoutFrameSize") ) ,
   m_ParameterMap     ( new EcalSimParameterMap(
               params.getParameter<double> ("simHitToPhotoelectronsBarrel") ,
               params.getParameter<double> ("simHitToPhotoelectronsEndcap") , 
               params.getParameter<double> ("photoelectronsToAnalogBarrel") ,
               params.getParameter<double> ("photoelectronsToAnalogEndcap") , 
               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 ( new 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 :
           new EBHitResponse( m_ParameterMap.get()  ,
                      &m_EBShape      ,
                      true            ,
                      m_apdParameters.get() ,
                      &m_APDShape       ) ) ,
   
   m_EBResponse ( new EBHitResponse( m_ParameterMap.get()  ,
                     &m_EBShape      ,
                     false           , // barrel
                     m_apdParameters.get() ,
                     &m_APDShape       ) ) ,

   m_EEResponse ( new EEHitResponse( m_ParameterMap.get(),
                     &m_EEShape       ) ) ,
   m_ESResponse ( new ESHitResponse( m_ParameterMap.get(), &m_ESShape ) ) ,
   m_ESOldResponse ( new CaloHitResponse( m_ParameterMap.get(), &m_ESShape ) ) ,

   m_addESNoise           ( params.getParameter<bool> ("doESNoise") ) ,
   m_PreMix1              ( params.getParameter<bool> ("EcalPreMixStage1") ) ,
   m_PreMix2              ( params.getParameter<bool> ("EcalPreMixStage2") ) ,

   m_doFastES             ( params.getParameter<bool> ("doFast"   ) ) ,

   m_doEB                 ( params.getParameter<bool> ("doEB"     ) ) ,
   m_doEE                 ( params.getParameter<bool> ("doEE"     ) ) ,
   m_doES                 ( params.getParameter<bool> ("doES"     ) ) ,

   m_ESElectronicsSim     ( m_doFastES ? nullptr :
                new ESElectronicsSim( m_addESNoise ) ) ,
     
   m_ESOldDigitizer       ( m_doFastES ? nullptr :
                new ESOldDigitizer( m_ESOldResponse.get()    , 
                        m_ESElectronicsSim.get() ,
                        m_addESNoise         ) ) ,
   
   m_ESElectronicsSimFast ( !m_doFastES ? nullptr :
                new ESElectronicsSimFast( m_addESNoise, 
                              m_PreMix1      ) ) ,

   m_ESDigitizer          ( !m_doFastES ? nullptr :
                new ESDigitizer( m_ESResponse.get()           ,
                         m_ESElectronicsSimFast.get() ,
                         m_addESNoise            ) ) ,

   m_APDDigitizer      ( nullptr ) ,
   m_BarrelDigitizer   ( nullptr ) ,
   m_EndcapDigitizer   ( nullptr ) ,
   m_ElectronicsSim    ( nullptr ) ,
   m_Coder             ( nullptr ) ,
   m_APDElectronicsSim ( nullptr ) ,
   m_APDCoder          ( nullptr ) ,
   m_Geometry          ( nullptr ) ,
   m_EBCorrNoise       ( { {nullptr, nullptr, nullptr} } ) ,

EcalDigiProducer::~EcalDigiProducer ( )

virtual

Definition at line 503 of file EcalDigiProducer.cc.

{}

Member Function Documentation

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 552 of file EcalDigiProducer.cc.

References accumulateCaloHits(), edm::Event::getByLabel(), m_hitsProducerTag, randomEngine(), and edm::Event::streamID().

Referenced by edm::DataMixingEcalDigiWorkerProd::addEcalSignals().

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

  edm::InputTag eeTag(m_hitsProducerTag, "EcalHitsEE");
  edm::Handle<std::vector<PCaloHit> > eeHandle;
  e.getByLabel(eeTag, eeHandle);

  edm::InputTag esTag(m_hitsProducerTag, "EcalHitsES");
  edm::Handle<std::vector<PCaloHit> > esHandle;
  e.getByLabel(esTag, esHandle);

  accumulateCaloHits(ebHandle, eeHandle, esHandle, 0, randomEngine(e.streamID()));
}

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

overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 570 of file EcalDigiProducer.cc.

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

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

  edm::InputTag eeTag(m_hitsProducerTag, "EcalHitsEE");
  edm::Handle<std::vector<PCaloHit> > eeHandle;
  e.getByLabel(eeTag, eeHandle);

  edm::InputTag esTag(m_hitsProducerTag, "EcalHitsES");
  edm::Handle<std::vector<PCaloHit> > esHandle;
  e.getByLabel(esTag, esHandle);

  accumulateCaloHits(ebHandle, eeHandle, esHandle, e.bunchCrossing(), randomEngine(streamID));
}

void EcalDigiProducer::accumulateCaloHits ( HitsHandle const &  ebHandle, HitsHandle const &  eeHandle, HitsHandle const &  esHandle, int  bunchCrossing, CLHEP::HepRandomEngine *  engine  )

private

Definition at line 529 of file EcalDigiProducer.cc.

References edm::HandleBase::isValid(), m_APDDigitizer, m_apdSeparateDigi, m_BarrelDigitizer, m_doEB, m_doEE, m_doES, m_doFastES, m_EndcapDigitizer, m_ESDigitizer, m_ESOldDigitizer, and edm::Handle< T >::product().

Referenced by accumulate().

                                                                                                                                                                        {
  if(m_doEB && ebHandle.isValid()) {
    m_BarrelDigitizer->add(*ebHandle.product(), bunchCrossing, engine);

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

  if(m_doEE && eeHandle.isValid()) {
    m_EndcapDigitizer->add(*eeHandle.product(), bunchCrossing, engine);
  }

  if(m_doES && esHandle.isValid()) {
    if(m_doFastES) {
      m_ESDigitizer->add(*esHandle.product(), bunchCrossing, engine);
    } else {
      m_ESOldDigitizer->add(*esHandle.product(), bunchCrossing, engine);
    }
  }
}

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

overridevirtual

Reimplemented from DigiAccumulatorMixMod.

Definition at line 636 of file EcalDigiProducer.cc.

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

Referenced by edm::DataMixingEcalDigiWorkerProd::beginLuminosityBlock().

{
   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( m_doEB ) {
     if( nullptr != m_APDResponse ) m_APDResponse->initialize(engine);
     m_EBResponse->initialize(engine);
   }
}

virtual void EcalDigiProducer::cacheEBDigis ( const EBDigiCollection *  ebDigiPtr ) const

inlineprivatevirtual

Reimplemented in EcalTBDigiProducer.

Definition at line 79 of file EcalDigiProducer.h.

Referenced by finalizeEvent().

{ }

virtual void EcalDigiProducer::cacheEEDigis ( const EEDigiCollection *  eeDigiPtr ) const

inlineprivatevirtual

Reimplemented in EcalTBDigiProducer.

Definition at line 80 of file EcalDigiProducer.h.

Referenced by finalizeEvent().

{ }

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

private

Definition at line 654 of file EcalDigiProducer.cc.

References EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), edm::EventSetup::get(), EcalADCToGeVConstant::getEBValue(), EcalADCToGeVConstant::getEEValue(), EcalCondDBWriter_cfi::laser, LogDebug, m_APDCoder, m_APDResponse, m_Coder, m_doES, m_doFastES, m_EBResponse, m_EBs25notCont, m_EEResponse, m_EEs25notCont, m_ESDigitizer, m_ESElectronicsSim, m_ESElectronicsSimFast, m_ESOldDigitizer, m_ESShape, m_useLCcorrection, edm::ESHandle< class >::product(), and ESShape::setGain().

Referenced by initializeEvent().

{
   // Pedestals from event setup

   edm::ESHandle<EcalPedestals>            dbPed   ;
   eventSetup.get<EcalPedestalsRcd>().get( dbPed ) ;
   const EcalPedestals* pedestals        ( dbPed.product() ) ;
  
   m_Coder->setPedestals( pedestals ) ;
   if( nullptr != m_APDCoder ) m_APDCoder->setPedestals( pedestals ) ;

   // Ecal Intercalibration Constants
   edm::ESHandle<EcalIntercalibConstantsMC>            pIcal   ;
   eventSetup.get<EcalIntercalibConstantsMCRcd>().get( pIcal ) ;
   const EcalIntercalibConstantsMC* ical             ( pIcal.product() ) ;
  
   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                                
   edm::ESHandle<EcalLaserDbService> laser;
   eventSetup.get<EcalLaserDbRecord>().get(laser);
   const edm::TimeValue_t eventTimeValue = event.time().value();

   m_EBResponse->setEventTime(eventTimeValue);
   m_EBResponse->setLaserConstants(laser.product(), m_useLCcorrection);

   m_EEResponse->setEventTime(eventTimeValue);
   m_EEResponse->setLaserConstants(laser.product(), m_useLCcorrection);

   // ADC -> GeV Scale
   edm::ESHandle<EcalADCToGeVConstant> pAgc;
   eventSetup.get<EcalADCToGeVConstantRcd>().get(pAgc);
   const EcalADCToGeVConstant* agc = pAgc.product();
  
   // Gain Ratios
   edm::ESHandle<EcalGainRatios> pRatio;
   eventSetup.get<EcalGainRatiosRcd>().get(pRatio);
   const EcalGainRatios* gr = pRatio.product();

   m_Coder->setGainRatios( gr );
   if( nullptr != m_APDCoder) m_APDCoder->setGainRatios( gr );

   EcalMGPAGainRatio * defaultRatios = new EcalMGPAGainRatio();

   double theGains[m_Coder->NGAINS+1];
   theGains[0] = 0.;
   theGains[3] = 1.;
   theGains[2] = defaultRatios->gain6Over1() ;
   theGains[1] = theGains[2]*(defaultRatios->gain12Over6()) ;

   LogDebug("EcalDigi") << " Gains: " << "\n" << " g1 = " << theGains[1] 
            << "\n" << " g2 = " << theGains[2] 
            << "\n" << " g3 = " << theGains[3] ;

   delete defaultRatios;

   const double EBscale (
      ( agc->getEBValue())*theGains[1]*(m_Coder->MAXADC)*m_EBs25notCont ) ;

   LogDebug("EcalDigi") << " GeV/ADC = " << agc->getEBValue() 
            << "\n" << " notCont = " << m_EBs25notCont 
            << "\n" << " saturation for EB = " << EBscale 
            << ", " << m_EBs25notCont ;

   const double EEscale (
      (agc->getEEValue())*theGains[1]*(m_Coder->MAXADC)*m_EEs25notCont ) ;

   LogDebug("EcalDigi") << " GeV/ADC = " << agc->getEEValue() 
            << "\n" << " notCont = " << m_EEs25notCont 
            << "\n" << " saturation for EB = " << EEscale 
            << ", " << m_EEs25notCont ;

   m_Coder->setFullScaleEnergy( EBscale , 
                EEscale   ) ;
   if( nullptr != m_APDCoder ) m_APDCoder->setFullScaleEnergy( EBscale ,
                             EEscale   ) ;

   if( nullptr != m_ESOldDigitizer ||
       nullptr != m_ESDigitizer       )
   {
      // ES condition objects
      edm::ESHandle<ESGain>                hesgain      ;
      edm::ESHandle<ESMIPToGeVConstant>    hesMIPToGeV  ;
      edm::ESHandle<ESPedestals>           hesPedestals ;
      edm::ESHandle<ESIntercalibConstants> hesMIPs      ;

      eventSetup.get<ESGainRcd>().               get( hesgain      ) ;
      eventSetup.get<ESMIPToGeVConstantRcd>().   get( hesMIPToGeV  ) ;
      eventSetup.get<ESPedestalsRcd>().          get( hesPedestals ) ;
      eventSetup.get<ESIntercalibConstantsRcd>().get( hesMIPs      ) ;

      const ESGain*                esgain     ( hesgain.product()      ) ;
      const ESPedestals*           espeds     ( hesPedestals.product() ) ;
      const ESIntercalibConstants* esmips     ( hesMIPs.product()      ) ;
      const ESMIPToGeVConstant*    esMipToGeV ( hesMIPToGeV.product()  ) ;
      const int ESGain ( 1.1 > esgain->getESGain() ? 1 : 2 ) ;
      const double ESMIPToGeV ( ( 1 == ESGain ) ?
                esMipToGeV->getESValueLow()  :
                esMipToGeV->getESValueHigh()   ) ; 
   
      if( m_doES ) {
        m_ESShape.setGain( ESGain );      
        if( !m_doFastES )
          {
            m_ESElectronicsSim->setGain(      ESGain     ) ;
            m_ESElectronicsSim->setPedestals( espeds     ) ;
            m_ESElectronicsSim->setMIPs(      esmips     ) ;
            m_ESElectronicsSim->setMIPToGeV(  ESMIPToGeV ) ;
          }
        else
          {
            m_ESDigitizer->setGain(               ESGain     ) ;
            m_ESElectronicsSimFast->setPedestals( espeds     ) ;
            m_ESElectronicsSimFast->setMIPs(      esmips     ) ;
            m_ESElectronicsSimFast->setMIPToGeV(  ESMIPToGeV ) ;
          }
      }
   }
}

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

private

Definition at line 779 of file EcalDigiProducer.cc.

References edm::EventSetup::get(), m_Geometry, and updateGeometry().

Referenced by initializeEvent().

{
   // TODO find a way to avoid doing this every event
   edm::ESHandle<CaloGeometry>               hGeometry   ;
   eventSetup.get<CaloGeometryRecord>().get( hGeometry ) ;

   const CaloGeometry* pGeometry = &*hGeometry;

   if( pGeometry != m_Geometry )
   {
      m_Geometry = pGeometry;
      updateGeometry();
   }
}

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

overridevirtual

Implements DigiAccumulatorMixMod.

Reimplemented in EcalTBDigiProducer.

Definition at line 588 of file EcalDigiProducer.cc.

References cacheEBDigis(), cacheEEDigis(), m_APDDigitizer, m_apdSeparateDigi, m_BarrelDigitizer, m_doEB, m_doEE, m_doES, m_doFastES, m_EBdigiCollection, m_EEdigiCollection, m_EndcapDigitizer, m_ESdigiCollection, m_ESDigitizer, m_ESOldDigitizer, randomEngine(), and edm::Event::streamID().

Referenced by EcalTBDigiProducer::finalizeEvent(), and edm::DataMixingEcalDigiWorkerProd::putEcal().

                                                                                {
   // Step B: Create empty output
   std::auto_ptr<EBDigiCollection> apdResult      ( !m_apdSeparateDigi || !m_doEB ? nullptr :
                            new EBDigiCollection() ) ;
   std::auto_ptr<EBDigiCollection> barrelResult   ( new EBDigiCollection() ) ;
   std::auto_ptr<EEDigiCollection> endcapResult   ( new EEDigiCollection() ) ;
   std::auto_ptr<ESDigiCollection> preshowerResult( new ESDigiCollection() ) ;
   
   // run the algorithm

   if( m_doEB ) {
     m_BarrelDigitizer->run( *barrelResult, randomEngine(event.streamID()) ) ;
     cacheEBDigis( &*barrelResult ) ;
     
     edm::LogInfo("DigiInfo") << "EB Digis: " << barrelResult->size() ;
     
     if( m_apdSeparateDigi ) {
       m_APDDigitizer->run( *apdResult, randomEngine(event.streamID()) ) ;
       edm::LogInfo("DigiInfo") << "APD Digis: " << apdResult->size() ;
     }
   }

   if( m_doEE ) {
     m_EndcapDigitizer->run( *endcapResult, randomEngine(event.streamID()) ) ;
     edm::LogInfo("EcalDigi") << "EE Digis: " << endcapResult->size() ;
     cacheEEDigis( &*endcapResult ) ;
   }
   if( m_doES ) {
     if(m_doFastES) {
       m_ESDigitizer->run( *preshowerResult, randomEngine(event.streamID()) ) ;
     } else {
       m_ESOldDigitizer->run( *preshowerResult, randomEngine(event.streamID()) ) ;
     }
     edm::LogInfo("EcalDigi") << "ES Digis: " << preshowerResult->size();
   }


   // Step D: Put outputs into event
   if( m_apdSeparateDigi ) {
     //event.put( apdResult,    m_apdDigiTag         ) ;
   }

   event.put( barrelResult,    m_EBdigiCollection ) ;
   event.put( endcapResult,    m_EEdigiCollection ) ;
   event.put( preshowerResult, m_ESdigiCollection ) ;
}

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

overridevirtual

Implements DigiAccumulatorMixMod.

Reimplemented in EcalTBDigiProducer.

Definition at line 507 of file EcalDigiProducer.cc.

References checkCalibrations(), checkGeometry(), m_APDDigitizer, m_apdSeparateDigi, m_BarrelDigitizer, m_doEB, m_doEE, m_doES, m_doFastES, m_EndcapDigitizer, m_ESDigitizer, and m_ESOldDigitizer.

Referenced by EcalTBDigiProducer::initializeEvent(), and edm::DataMixingEcalDigiWorkerProd::initializeEvent().

                                                                                        {
   checkGeometry( eventSetup );
   checkCalibrations( event, eventSetup );
   if( m_doEB ) {
     m_BarrelDigitizer->initializeHits();
     if(m_apdSeparateDigi) {
       m_APDDigitizer->initializeHits();
     }
   }
   if( m_doEE ) {
     m_EndcapDigitizer->initializeHits();
   }
   if( m_doES ) {
     if(m_doFastES) {
       m_ESDigitizer->initializeHits();
     } else {
       m_ESOldDigitizer->initializeHits();
     }
   }
}

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

private

Definition at line 831 of file EcalDigiProducer.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 EcalDigiProducer::setEBNoiseSignalGenerator

(

EcalBaseSignalGenerator *

noiseGenerator

)

Definition at line 845 of file EcalDigiProducer.cc.

References m_BarrelDigitizer.

                                                                                         {
  //noiseGenerator->setParameterMap(theParameterMap);
  if(nullptr != m_BarrelDigitizer) m_BarrelDigitizer->setNoiseSignalGenerator(noiseGenerator);
}

void EcalDigiProducer::setEENoiseSignalGenerator

(

EcalBaseSignalGenerator *

noiseGenerator

)

Definition at line 850 of file EcalDigiProducer.cc.

References m_EndcapDigitizer.

                                                                                         {
  //noiseGenerator->setParameterMap(theParameterMap);
  if(nullptr != m_EndcapDigitizer) m_EndcapDigitizer->setNoiseSignalGenerator(noiseGenerator);
}

void EcalDigiProducer::setESNoiseSignalGenerator

(

EcalBaseSignalGenerator *

noiseGenerator

)

Definition at line 855 of file EcalDigiProducer.cc.

References m_ESDigitizer.

                                                                                         {
  //noiseGenerator->setParameterMap(theParameterMap);
  if(nullptr != m_ESDigitizer) m_ESDigitizer->setNoiseSignalGenerator(noiseGenerator);  
}

void EcalDigiProducer::updateGeometry ( )

private

Definition at line 795 of file EcalDigiProducer.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, CaloGeometry::getSubdetectorGeometry(), CaloSubdetectorGeometry::getValidDetIds(), m_APDResponse, m_doEB, m_doEE, m_doES, m_doFastES, m_EBResponse, m_EEResponse, m_ESDigitizer, m_ESOldDigitizer, m_ESOldResponse, m_ESResponse, m_Geometry, and nullptr.

Referenced by checkGeometry().

{
  if( m_doEB ) {
    if( nullptr != m_APDResponse ) m_APDResponse->setGeometry(
      m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalBarrel    ) ) ;
    m_EBResponse->setGeometry(
      m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalBarrel    ) ) ;
  }
  if( m_doEE ) {
    m_EEResponse->setGeometry(
      m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalEndcap    ) ) ;
  }
  if( m_doES ) {
    m_ESResponse->setGeometry(
      m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalPreshower    ) ) ;
    m_ESOldResponse->setGeometry( m_Geometry ) ;
   
   const std::vector<DetId>* theESDets ( 
      nullptr != m_Geometry->getSubdetectorGeometry(DetId::Ecal, EcalPreshower) ?
      &m_Geometry->getSubdetectorGeometry(DetId::Ecal, EcalPreshower)->getValidDetIds() : nullptr ) ;

   if( !m_doFastES ) 
   {
      if( nullptr != m_ESOldDigitizer &&
      nullptr != theESDets             )
     m_ESOldDigitizer->setDetIds( *theESDets ) ;
   }
   else
   {
      if( nullptr != m_ESDigitizer &&
      nullptr != theESDets         )
     m_ESDigitizer->setDetIds( *theESDets ) ; 
   }
  }
}

Member Data Documentation

const bool EcalDigiProducer::m_addESNoise

private

Definition at line 125 of file EcalDigiProducer.h.

std::unique_ptr<EcalCoder> EcalDigiProducer::m_APDCoder

private

Definition at line 146 of file EcalDigiProducer.h.

Referenced by checkCalibrations().

const std::string EcalDigiProducer::m_apdDigiTag

private

Definition at line 114 of file EcalDigiProducer.h.

std::unique_ptr<EBDigitizer> EcalDigiProducer::m_APDDigitizer

private

Definition at line 138 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), finalizeEvent(), and initializeEvent().

std::unique_ptr<EcalElectronicsSim> EcalDigiProducer::m_APDElectronicsSim

private

Definition at line 145 of file EcalDigiProducer.h.

std::unique_ptr<const APDSimParameters> EcalDigiProducer::m_apdParameters

private

Definition at line 115 of file EcalDigiProducer.h.

std::unique_ptr<EBHitResponse> EcalDigiProducer::m_APDResponse

private

Definition at line 117 of file EcalDigiProducer.h.

Referenced by beginLuminosityBlock(), checkCalibrations(), and updateGeometry().

const bool EcalDigiProducer::m_apdSeparateDigi

private

Definition at line 105 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), finalizeEvent(), and initializeEvent().

const APDShape EcalDigiProducer::m_APDShape

private

Definition at line 93 of file EcalDigiProducer.h.

std::unique_ptr<EBDigitizer> EcalDigiProducer::m_BarrelDigitizer

private

Definition at line 139 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), finalizeEvent(), initializeEvent(), and setEBNoiseSignalGenerator().

std::unique_ptr<EcalCoder> EcalDigiProducer::m_Coder

private

Definition at line 143 of file EcalDigiProducer.h.

Referenced by checkCalibrations().

const bool EcalDigiProducer::m_doEB

private

Definition at line 131 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), beginLuminosityBlock(), finalizeEvent(), initializeEvent(), and updateGeometry().

const bool EcalDigiProducer::m_doEE

private

Definition at line 131 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), finalizeEvent(), initializeEvent(), and updateGeometry().

const bool EcalDigiProducer::m_doES

private

Definition at line 131 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), checkCalibrations(), finalizeEvent(), initializeEvent(), and updateGeometry().

const bool EcalDigiProducer::m_doFastES

private

Definition at line 129 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), checkCalibrations(), finalizeEvent(), initializeEvent(), and updateGeometry().

std::array< std::unique_ptr<CorrelatedNoisifier<EcalCorrMatrix> >, 3 > EcalDigiProducer::m_EBCorrNoise

private

Definition at line 150 of file EcalDigiProducer.h.

const std::string EcalDigiProducer::m_EBdigiCollection

private

Definition at line 98 of file EcalDigiProducer.h.

Referenced by finalizeEvent().

std::unique_ptr<EBHitResponse> EcalDigiProducer::m_EBResponse

protected

Definition at line 119 of file EcalDigiProducer.h.

Referenced by beginLuminosityBlock(), checkCalibrations(), EcalTBDigiProducer::setPhaseShift(), and updateGeometry().

const double EcalDigiProducer::m_EBs25notCont

private

Definition at line 107 of file EcalDigiProducer.h.

Referenced by checkCalibrations().

const EBShape EcalDigiProducer::m_EBShape

private

Definition at line 94 of file EcalDigiProducer.h.

std::array< std::unique_ptr<CorrelatedNoisifier<EcalCorrMatrix> >, 3 > EcalDigiProducer::m_EECorrNoise

private

Definition at line 151 of file EcalDigiProducer.h.

const std::string EcalDigiProducer::m_EEdigiCollection

private

Definition at line 99 of file EcalDigiProducer.h.

Referenced by finalizeEvent().

std::unique_ptr<EEHitResponse> EcalDigiProducer::m_EEResponse

protected

Definition at line 120 of file EcalDigiProducer.h.

Referenced by checkCalibrations(), EcalTBDigiProducer::setPhaseShift(), and updateGeometry().

const double EcalDigiProducer::m_EEs25notCont

private

Definition at line 108 of file EcalDigiProducer.h.

Referenced by checkCalibrations().

const EEShape EcalDigiProducer::m_EEShape

private

Definition at line 95 of file EcalDigiProducer.h.

std::unique_ptr<EcalElectronicsSim> EcalDigiProducer::m_ElectronicsSim

private

Definition at line 142 of file EcalDigiProducer.h.

std::unique_ptr<EEDigitizer> EcalDigiProducer::m_EndcapDigitizer

private

Definition at line 140 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), finalizeEvent(), initializeEvent(), and setEENoiseSignalGenerator().

const std::string EcalDigiProducer::m_ESdigiCollection

private

Definition at line 100 of file EcalDigiProducer.h.

Referenced by finalizeEvent().

std::unique_ptr<ESDigitizer> EcalDigiProducer::m_ESDigitizer

private

Definition at line 136 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), checkCalibrations(), finalizeEvent(), initializeEvent(), setESNoiseSignalGenerator(), and updateGeometry().

std::unique_ptr<ESElectronicsSim> EcalDigiProducer::m_ESElectronicsSim

private

Definition at line 133 of file EcalDigiProducer.h.

Referenced by checkCalibrations().

std::unique_ptr<ESElectronicsSimFast> EcalDigiProducer::m_ESElectronicsSimFast

private

Definition at line 135 of file EcalDigiProducer.h.

Referenced by checkCalibrations().

std::unique_ptr<ESOldDigitizer> EcalDigiProducer::m_ESOldDigitizer

private

Definition at line 134 of file EcalDigiProducer.h.

Referenced by accumulateCaloHits(), checkCalibrations(), finalizeEvent(), initializeEvent(), and updateGeometry().

std::unique_ptr<CaloHitResponse> EcalDigiProducer::m_ESOldResponse

private

Definition at line 123 of file EcalDigiProducer.h.

Referenced by updateGeometry().

std::unique_ptr<ESHitResponse> EcalDigiProducer::m_ESResponse

private

Definition at line 122 of file EcalDigiProducer.h.

Referenced by updateGeometry().

ESShape EcalDigiProducer::m_ESShape

private

Definition at line 96 of file EcalDigiProducer.h.

Referenced by checkCalibrations().

const CaloGeometry* EcalDigiProducer::m_Geometry

private

Definition at line 148 of file EcalDigiProducer.h.

Referenced by checkGeometry(), and updateGeometry().

const std::string EcalDigiProducer::m_hitsProducerTag

private

Definition at line 101 of file EcalDigiProducer.h.

Referenced by accumulate().

std::unique_ptr<const EcalSimParameterMap> EcalDigiProducer::m_ParameterMap

protected

Definition at line 112 of file EcalDigiProducer.h.

Referenced by EcalTBDigiProducer::setPhaseShift().

const bool EcalDigiProducer::m_PreMix1

private

Definition at line 126 of file EcalDigiProducer.h.

const bool EcalDigiProducer::m_PreMix2

private

Definition at line 127 of file EcalDigiProducer.h.

const unsigned int EcalDigiProducer::m_readoutFrameSize

private

Definition at line 110 of file EcalDigiProducer.h.

bool EcalDigiProducer::m_useLCcorrection

private

Definition at line 103 of file EcalDigiProducer.h.

Referenced by checkCalibrations().

std::vector<CLHEP::HepRandomEngine*> EcalDigiProducer::randomEngines_

private

Definition at line 153 of file EcalDigiProducer.h.

Referenced by randomEngine().