#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 49 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

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

Definition at line 48 of file EcalSignalGenerator.h.

Member Enumeration Documentation

template<class ECALDIGITIZERTRAITS>

anonymous enum

private

Enumerator
NBITS
MAXADC
ADCGAINSWITCH
NGAINS

Definition at line 181 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 51 of file EcalSignalGenerator.h.

51 : EcalBaseSignalGenerator() {}

EcalBaseSignalGenerator::EcalBaseSignalGenerator

EcalBaseSignalGenerator()

Definition: EcalBaseSignalGenerator.h:10

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 53 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 69 of file EcalSignalGenerator.h.

69 {}

Member Function Documentation

template<class ECALDIGITIZERTRAITS>

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

inlinevirtual

Definition at line 133 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 176 of file EcalSignalGenerator.h.

176 {}

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 177 of file EcalSignalGenerator.h.

177 {}

template<class ECALDIGITIZERTRAITS>

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

inlineprivate

Definition at line 203 of file EcalSignalGenerator.h.

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

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

EcalSignalGenerator::m_maxEneEE

double m_maxEneEE

Definition: EcalSignalGenerator.h:238

DetId::subdetId

constexpr int subdetId() const

get the contents of the subdetector field (not cast into any detector's numbering enum) ...

Definition: DetId.h:41

EcalSignalGenerator::m_maxEneEB

double m_maxEneEB

Definition: EcalSignalGenerator.h:237

EcalBarrel

Definition: EcalSubdetector.h:10

template<class ECALDIGITIZERTRAITS>

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

inlineprivate

Definition at line 190 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 71 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 205 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 167 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 240 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

const ESGain* EcalSignalGenerator< ECALDIGITIZERTRAITS >::esgain

private

Definition at line 225 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

int EcalSignalGenerator< ECALDIGITIZERTRAITS >::ESgain

private

Definition at line 228 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

const ESIntercalibConstants* EcalSignalGenerator< ECALDIGITIZERTRAITS >::esmips

private

Definition at line 226 of file EcalSignalGenerator.h.

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

template<class ECALDIGITIZERTRAITS>

const ESMIPToGeVConstant* EcalSignalGenerator< ECALDIGITIZERTRAITS >::esMipToGeV

private

Definition at line 227 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::ESMIPToGeV

private

Definition at line 229 of file EcalSignalGenerator.h.

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

template<class ECALDIGITIZERTRAITS>

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

private

Definition at line 213 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

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

private

Definition at line 221 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

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

private

Definition at line 223 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

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

private

Definition at line 222 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

const EcalIntercalibConstantsMC* EcalSignalGenerator< ECALDIGITIZERTRAITS >::ical

private

Definition at line 241 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 231 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_EEs25notCont

private

Definition at line 232 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_maxEneEB

private

Definition at line 237 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_maxEneEE

private

Definition at line 238 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_peToABarrel

private

Definition at line 234 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::m_peToAEndcap

private

Definition at line 235 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

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

private

Definition at line 215 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

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

private

Definition at line 214 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

double EcalSignalGenerator< ECALDIGITIZERTRAITS >::theDefaultGains[NGAINS]

private

Definition at line 243 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

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

private

these fields are set in initializeEvent()

Definition at line 210 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

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

private

Definition at line 211 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

edm::InputTag EcalSignalGenerator< ECALDIGITIZERTRAITS >::theInputTag

private

these come from the ParameterSet

Definition at line 218 of file EcalSignalGenerator.h.

template<class ECALDIGITIZERTRAITS>

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

private

Definition at line 219 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 >