#include <EcalSignalGenerator.h>

Inheritance diagram for EcalSignalGenerator< ECALDIGITIZERTRAITS >:

Public Types
typedef ECALDIGITIZERTRAITS::DigiCollection	COLLECTION

typedef ECALDIGITIZERTRAITS::Digi	DIGI

Public Member Functions
	EcalSignalGenerator ()

	EcalSignalGenerator (const edm::InputTag &inputTag, const edm::EDGetTokenT< COLLECTION > &t, const double EBs25notCont, const double EEs25notCont, const double peToABarrel, const double peToAEndcap)

virtual void	fill (edm::ModuleCallingContext const *mcc)

void	initializeEvent (const edm::Event event, const edm::EventSetup eventSetup)

void	initializeEvent (const edm::EventPrincipal eventPrincipal, const edm::EventSetup eventSetup)
	some users use EventPrincipals, not Events. We support both More...

	~EcalSignalGenerator () override

Public Member Functions inherited from EcalBaseSignalGenerator
	EcalBaseSignalGenerator ()

	~EcalBaseSignalGenerator () override

Public Member Functions inherited from CaloVNoiseSignalGenerator
	CaloVNoiseSignalGenerator ()

bool	contains (const DetId &detId) const

void	fillEvent (CLHEP::HepRandomEngine *)
	fill theNoiseSignals with one event's worth of noise, in units of pe More...

void	fillEvent ()

void	getNoiseSignals (std::vector< CaloSamples > &noiseSignals)

void	setNoiseSignals (const std::vector< CaloSamples > &noiseSignals)

virtual	~CaloVNoiseSignalGenerator ()

Private Types
enum	{ NBITS = 12, MAXADC = 4095, ADCGAINSWITCH = 4079, NGAINS = 3 }

Private Member Functions
void	fillNoiseSignals () override

void	fillNoiseSignals (CLHEP::HepRandomEngine *) override

double	fullScaleEnergy (const DetId &detId) const

const std::vector< float >	GetGainRatios (const DetId &detid)

double	peToAConversion (const DetId &detId) const

CaloSamples	samplesInPE (const DIGI &digi)

bool	validDigi (const DIGI &digi)

Private Attributes
const EcalADCToGeVConstant *	agc

const ESGain *	esgain

int	ESgain

const ESIntercalibConstants *	esmips

const ESMIPToGeVConstant *	esMipToGeV

double	ESMIPToGeV

edm::ESHandle< EcalGainRatios >	grHandle

edm::ESHandle< ESGain >	hesgain

edm::ESHandle< ESIntercalibConstants >	hesMIPs

edm::ESHandle< ESMIPToGeVConstant >	hesMIPToGeV

const EcalIntercalibConstantsMC *	ical

double	m_EBs25notCont

double	m_EEs25notCont

double	m_maxEneEB

double	m_maxEneEE

double	m_peToABarrel

double	m_peToAEndcap

edm::ESHandle< EcalADCToGeVConstant >	pAgc

edm::ESHandle< EcalIntercalibConstantsMC >	pIcal

double	theDefaultGains [NGAINS]

const edm::Event *	theEvent
	these fields are set in initializeEvent() More...

const edm::EventPrincipal *	theEventPrincipal

edm::InputTag	theInputTag
	these come from the ParameterSet More...

edm::EDGetTokenT< COLLECTION >	tok_

Additional Inherited Members
Protected Attributes inherited from CaloVNoiseSignalGenerator
std::vector< CaloSamples >	theNoiseSignals

Detailed Description

template<class ECALDIGITIZERTRAITS>
class EcalSignalGenerator< ECALDIGITIZERTRAITS >

Definition at line 46 of file EcalSignalGenerator.h.

Member Typedef Documentation

template<class ECALDIGITIZERTRAITS>

typedef ECALDIGITIZERTRAITS::DigiCollection EcalSignalGenerator< ECALDIGITIZERTRAITS >::COLLECTION

Definition at line 50 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

typedef ECALDIGITIZERTRAITS::Digi EcalSignalGenerator< ECALDIGITIZERTRAITS >::DIGI

Definition at line 49 of file EcalSignalGenerator.h.

Member Enumeration Documentation

template<class ECALDIGITIZERTRAITS>

anonymous enum

private

Enumerator
NBITS
MAXADC
ADCGAINSWITCH
NGAINS

Definition at line 191 of file EcalSignalGenerator.h.

        { NBITS         =   12 , // number of available bits
      MAXADC        = 4095 , // 2^12 -1,  adc max range
      ADCGAINSWITCH = 4079 , // adc gain switch
      NGAINS        =    3 };  // number of electronic gains

Constructor & Destructor Documentation

template<class ECALDIGITIZERTRAITS>

EcalSignalGenerator< ECALDIGITIZERTRAITS >::EcalSignalGenerator ( )

inline

Definition at line 52 of file EcalSignalGenerator.h.

52 :EcalBaseSignalGenerator() { }

EcalBaseSignalGenerator::EcalBaseSignalGenerator

EcalBaseSignalGenerator()

Definition: EcalBaseSignalGenerator.h:11

template<class ECALDIGITIZERTRAITS>

EcalSignalGenerator< ECALDIGITIZERTRAITS >::EcalSignalGenerator	(	const edm::InputTag &	inputTag,
		const edm::EDGetTokenT< COLLECTION > &	t,
		const double	EBs25notCont,
		const double	EEs25notCont,
		const double	peToABarrel,
		const double	peToAEndcap
	)

inline

Definition at line 54 of file EcalSignalGenerator.h.

     : EcalBaseSignalGenerator(), 
       theEvent(nullptr), 
       theEventPrincipal(nullptr), 
       theInputTag(inputTag), 
       tok_(t)
       { 
     EcalMGPAGainRatio * defaultRatios = new EcalMGPAGainRatio();
     theDefaultGains[2] = defaultRatios->gain6Over1() ;
     theDefaultGains[1] = theDefaultGains[2]*(defaultRatios->gain12Over6()) ;
     m_EBs25notCont = EBs25notCont; 
     m_EEs25notCont = EEs25notCont; 
     m_peToABarrel = peToABarrel;
     m_peToAEndcap = peToAEndcap;
 
       }

template<class ECALDIGITIZERTRAITS>

EcalSignalGenerator< ECALDIGITIZERTRAITS >::~EcalSignalGenerator ( )

inlineoverride

Definition at line 71 of file EcalSignalGenerator.h.

71 {}

Member Function Documentation

template<class ECALDIGITIZERTRAITS>

virtual void EcalSignalGenerator< ECALDIGITIZERTRAITS >::fill ( edm::ModuleCallingContext const * mcc )

inlinevirtual

Definition at line 130 of file EcalSignalGenerator.h.

Referenced by PreMixingEcalWorker::addPileups().

   {
 
     theNoiseSignals.clear();
     edm::Handle<COLLECTION> pDigis;
     const COLLECTION *  digis = nullptr;
     // try accessing by whatever is set, Event or EventPrincipal
     if(theEvent) 
      {
       if( theEvent->getByToken(tok_, pDigis) ) {
         digis = pDigis.product(); // get a ptr to the product
       }
       else
       {
         throw cms::Exception("EcalSignalGenerator") << "Cannot find input data " << theInputTag;
       }
     }
     else if(theEventPrincipal)
     {
        std::shared_ptr<edm::Wrapper<COLLECTION>  const> digisPTR =
           edm::getProductByTag<COLLECTION>(*theEventPrincipal, theInputTag, mcc );
        if(digisPTR) {
           digis = digisPTR->product();
        }
     }
     else
     {
       throw cms::Exception("EcalSignalGenerator") << "No Event or EventPrincipal was set";
     }
 
     if (digis)
     {
       // loop over digis, adding these to the existing maps
       for(typename COLLECTION::const_iterator it  = digis->begin();
           it != digis->end(); ++it) 
       {
     // need to convert to something useful
     if(validDigi(*it)){
       theNoiseSignals.push_back(samplesInPE(*it));
     }
       }
     }
     //else { std::cout << " NO digis for this input: " << theInputTag << std::endl;}
   }

template<class ECALDIGITIZERTRAITS>

void EcalSignalGenerator< ECALDIGITIZERTRAITS >::fillNoiseSignals ( )

inlineoverrideprivatevirtual

Implements CaloVNoiseSignalGenerator.

Definition at line 186 of file EcalSignalGenerator.h.

186 {}

template<class ECALDIGITIZERTRAITS>

void EcalSignalGenerator< ECALDIGITIZERTRAITS >::fillNoiseSignals ( CLHEP::HepRandomEngine * )

inlineoverrideprivatevirtual

if you want to fill signals on demand, override this subclass is responsible for clearing theNoiseSignals before adding

Implements CaloVNoiseSignalGenerator.

Definition at line 187 of file EcalSignalGenerator.h.

187 {}

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::fullScaleEnergy ( const DetId & detId ) const

inlineprivate

Definition at line 213 of file EcalSignalGenerator.h.

Referenced by EcalSignalGenerator< EBDigitizerTraits >::samplesInPE(), and EcalSignalGenerator< EEDigitizerTraits >::samplesInPE().

   {
      return detId.subdetId() == EcalBarrel ? m_maxEneEB : m_maxEneEE ;
   }

template<class ECALDIGITIZERTRAITS>

const std::vector<float> EcalSignalGenerator< ECALDIGITIZERTRAITS >::GetGainRatios ( const DetId & detid )

inlineprivate

Definition at line 198 of file EcalSignalGenerator.h.

Referenced by EcalSignalGenerator< EBDigitizerTraits >::samplesInPE(), and EcalSignalGenerator< EEDigitizerTraits >::samplesInPE().

                                                             {
 
     std::vector<float> gainRatios(4);
     // get gain ratios  
     EcalMGPAGainRatio theRatio= (*grHandle)[detid];
     
     gainRatios[0] = 0.;
     gainRatios[3] = 1.;
     gainRatios[2] = theRatio.gain6Over1();
     gainRatios[1] = theRatio.gain6Over1()  * theRatio.gain12Over6();
 
     return gainRatios;
   }

template<class ECALDIGITIZERTRAITS>

void EcalSignalGenerator< ECALDIGITIZERTRAITS >::initializeEvent	(	const edm::Event *	event,
		const edm::EventSetup *	eventSetup
	)

inline

Definition at line 74 of file EcalSignalGenerator.h.

Referenced by PreMixingEcalWorker::addPileups().

   {
     theEvent = event;
     eventSetup->get<EcalGainRatiosRcd>().get(grHandle); // find the gains
     // Ecal Intercalibration Constants
     eventSetup->get<EcalIntercalibConstantsMCRcd>().get( pIcal ) ;
     ical = pIcal.product();
     // adc to GeV
     eventSetup->get<EcalADCToGeVConstantRcd>().get(pAgc);
     agc = pAgc.product();
 
     m_maxEneEB = (agc->getEBValue())*theDefaultGains[1]*MAXADC*m_EBs25notCont  ;
     m_maxEneEE = (agc->getEEValue())*theDefaultGains[1]*MAXADC*m_EEs25notCont  ;
 
     //ES
     eventSetup->get<ESGainRcd>().               get( hesgain      ) ;
     eventSetup->get<ESMIPToGeVConstantRcd>().   get( hesMIPToGeV  ) ;
     eventSetup->get<ESIntercalibConstantsRcd>().get( hesMIPs      ) ;
 
     esgain     = hesgain.product()      ;
     esmips     = hesMIPs.product()      ;
     esMipToGeV = hesMIPToGeV.product()  ;
     if( 1.1 > esgain->getESGain() ) ESgain = 1;
     else ESgain = 2;
     if( ESgain ==1 ) ESMIPToGeV = esMipToGeV->getESValueLow();
     else ESMIPToGeV = esMipToGeV->getESValueHigh();
   }

template<class ECALDIGITIZERTRAITS>

void EcalSignalGenerator< ECALDIGITIZERTRAITS >::initializeEvent	(	const edm::EventPrincipal *	eventPrincipal,
		const edm::EventSetup *	eventSetup
	)

inline

some users use EventPrincipals, not Events. We support both

Definition at line 103 of file EcalSignalGenerator.h.

   {
     theEventPrincipal = eventPrincipal;
     eventSetup->get<EcalGainRatiosRcd>().get(grHandle);  // find the gains
     // Ecal Intercalibration Constants
     eventSetup->get<EcalIntercalibConstantsMCRcd>().get( pIcal ) ;
     ical = pIcal.product();
     // adc to GeV
     eventSetup->get<EcalADCToGeVConstantRcd>().get(pAgc);
     agc = pAgc.product();
     m_maxEneEB = (agc->getEBValue())*theDefaultGains[1]*MAXADC*m_EBs25notCont  ;
     m_maxEneEE = (agc->getEEValue())*theDefaultGains[1]*MAXADC*m_EEs25notCont  ;
 
     //ES
     eventSetup->get<ESGainRcd>().               get( hesgain      ) ;
     eventSetup->get<ESMIPToGeVConstantRcd>().   get( hesMIPToGeV  ) ;
     eventSetup->get<ESIntercalibConstantsRcd>().get( hesMIPs      ) ;
 
     esgain     = hesgain.product()      ;
     esmips     = hesMIPs.product()      ;
     esMipToGeV = hesMIPToGeV.product()  ;
     if( 1.1 > esgain->getESGain() ) ESgain = 1;
     else ESgain = 2;
     if( ESgain ==1 ) ESMIPToGeV = esMipToGeV->getESValueLow();
     else ESMIPToGeV = esMipToGeV->getESValueHigh();
   }

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::peToAConversion ( const DetId & detId ) const

inlineprivate

Definition at line 219 of file EcalSignalGenerator.h.

Referenced by EcalSignalGenerator< EBDigitizerTraits >::samplesInPE(), and EcalSignalGenerator< EEDigitizerTraits >::samplesInPE().

   {
      return detId.subdetId() == EcalBarrel ? m_peToABarrel : m_peToAEndcap ;
   }

template<class ECALDIGITIZERTRAITS>

CaloSamples EcalSignalGenerator< ECALDIGITIZERTRAITS >::samplesInPE ( const DIGI & digi )

private

template<class ECALDIGITIZERTRAITS>

bool EcalSignalGenerator< ECALDIGITIZERTRAITS >::validDigi ( const DIGI & digi )

inlineprivate

Definition at line 177 of file EcalSignalGenerator.h.

   {
     int DigiSum = 0;
     for(int id = 0; id<digi.size(); id++) {
       if(digi[id].adc() > 0) ++DigiSum;
     }
     return(DigiSum>0);
   }

Member Data Documentation

template<class ECALDIGITIZERTRAITS>

const EcalADCToGeVConstant* EcalSignalGenerator< ECALDIGITIZERTRAITS >::agc

private

Definition at line 256 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

const ESGain* EcalSignalGenerator< ECALDIGITIZERTRAITS >::esgain

private

Definition at line 241 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

int EcalSignalGenerator< ECALDIGITIZERTRAITS >::ESgain

private

Definition at line 244 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

const ESIntercalibConstants* EcalSignalGenerator< ECALDIGITIZERTRAITS >::esmips

private

Definition at line 242 of file EcalSignalGenerator.h.

Referenced by EcalSignalGenerator< ESDigitizerTraits >::samplesInPE().

template<class ECALDIGITIZERTRAITS>

const ESMIPToGeVConstant* EcalSignalGenerator< ECALDIGITIZERTRAITS >::esMipToGeV

private

Definition at line 243 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::ESMIPToGeV

private

Definition at line 245 of file EcalSignalGenerator.h.

Referenced by EcalSignalGenerator< ESDigitizerTraits >::samplesInPE().

template<class ECALDIGITIZERTRAITS>

edm::ESHandle<EcalGainRatios> EcalSignalGenerator< ECALDIGITIZERTRAITS >::grHandle

private

Definition at line 229 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

edm::ESHandle<ESGain> EcalSignalGenerator< ECALDIGITIZERTRAITS >::hesgain

private

Definition at line 237 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

edm::ESHandle<ESIntercalibConstants> EcalSignalGenerator< ECALDIGITIZERTRAITS >::hesMIPs

private

Definition at line 239 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

edm::ESHandle<ESMIPToGeVConstant> EcalSignalGenerator< ECALDIGITIZERTRAITS >::hesMIPToGeV

private

Definition at line 238 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

const EcalIntercalibConstantsMC* EcalSignalGenerator< ECALDIGITIZERTRAITS >::ical

private

Definition at line 257 of file EcalSignalGenerator.h.

Referenced by EcalSignalGenerator< EBDigitizerTraits >::samplesInPE(), and EcalSignalGenerator< EEDigitizerTraits >::samplesInPE().

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_EBs25notCont

private

Definition at line 247 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_EEs25notCont

private

Definition at line 248 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_maxEneEB

private

Definition at line 253 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_maxEneEE

private

Definition at line 254 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_peToABarrel

private

Definition at line 250 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_peToAEndcap

private

Definition at line 251 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

edm::ESHandle<EcalADCToGeVConstant> EcalSignalGenerator< ECALDIGITIZERTRAITS >::pAgc

private

Definition at line 231 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

edm::ESHandle<EcalIntercalibConstantsMC> EcalSignalGenerator< ECALDIGITIZERTRAITS >::pIcal

private

Definition at line 230 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::theDefaultGains[NGAINS]

private

Definition at line 259 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

const edm::Event* EcalSignalGenerator< ECALDIGITIZERTRAITS >::theEvent

private

these fields are set in initializeEvent()

Definition at line 226 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

const edm::EventPrincipal* EcalSignalGenerator< ECALDIGITIZERTRAITS >::theEventPrincipal

private

Definition at line 227 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

edm::InputTag EcalSignalGenerator< ECALDIGITIZERTRAITS >::theInputTag

private

these come from the ParameterSet

Definition at line 234 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

edm::EDGetTokenT<COLLECTION> EcalSignalGenerator< ECALDIGITIZERTRAITS >::tok_

private

Definition at line 235 of file EcalSignalGenerator.h.

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

template<class ECALDIGITIZERTRAITS> class EcalSignalGenerator< ECALDIGITIZERTRAITS >

Member Typedef Documentation

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

template<class ECALDIGITIZERTRAITS>
class EcalSignalGenerator< ECALDIGITIZERTRAITS >