CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Public Member Functions | Protected Attributes | Private Types | Private Member Functions | Private Attributes
EcalDigiProducer Class Reference

#include <EcalDigiProducer.h>

Inheritance diagram for EcalDigiProducer:
DigiAccumulatorMixMod EcalTBDigiProducer

Public Member Functions

virtual void accumulate (edm::Event const &e, edm::EventSetup const &c)
 
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::one::EDProducerBase &mixMod, edm::ConsumesCollector &iC)
 
 EcalDigiProducer (const edm::ParameterSet &params, edm::ConsumesCollector &iC)
 
virtual void finalizeEvent (edm::Event &e, edm::EventSetup const &c)
 
virtual void initializeEvent (edm::Event const &e, edm::EventSetup const &c)
 
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)
 
virtual ~DigiAccumulatorMixMod ()
 

Protected Attributes

EBHitResponse * m_EBResponse
 
EEHitResponse * m_EEResponse
 
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
 
EcalCoder * m_APDCoder
 
const std::string m_apdDigiTag
 
EBDigitizer * m_APDDigitizer
 
EcalElectronicsSim * m_APDElectronicsSim
 
const APDSimParameters * m_apdParameters
 
EBHitResponse * m_APDResponse
 
const bool m_apdSeparateDigi
 
const APDShape m_APDShape
 
EBDigitizer * m_BarrelDigitizer
 
EcalCoder * m_Coder
 
const bool m_doFastES
 
CorrelatedNoisifier
< EcalCorrMatrix > * 		m_EBCorrNoise [3]
 
const std::string m_EBdigiCollection
 
const double m_EBs25notCont
 
const EBShape m_EBShape
 
CorrelatedNoisifier
< EcalCorrMatrix > * 		m_EECorrNoise [3]
 
const std::string m_EEdigiCollection
 
const double m_EEs25notCont
 
const EEShape m_EEShape
 
EcalElectronicsSim * m_ElectronicsSim
 
EEDigitizer * m_EndcapDigitizer
 
const std::string m_ESdigiCollection
 
ESDigitizer * m_ESDigitizer
 
ESElectronicsSim * m_ESElectronicsSim
 
ESElectronicsSimFast * m_ESElectronicsSimFast
 
ESOldDigitizer * m_ESOldDigitizer
 
CaloHitResponse * m_ESOldResponse
 
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::one::EDProducerBase &  mixMod,
edm::ConsumesCollector &  iC 
)

Definition at line 57 of file EcalDigiProducer.cc.

References column, edm::ConsumesCollector::consumes(), alignCSCRings::e, edm::ParameterSet::getParameter(), cmsHarvester::index, m_APDCoder, m_apdDigiTag, m_APDDigitizer, m_APDElectronicsSim, m_APDResponse, m_apdSeparateDigi, m_BarrelDigitizer, m_Coder, m_EBCorrNoise, m_EBdigiCollection, m_EBResponse, m_EECorrNoise, m_EEdigiCollection, m_EEResponse, m_ElectronicsSim, m_EndcapDigitizer, m_ESdigiCollection, m_hitsProducerTag, m_ParameterMap, m_PreMix1, m_readoutFrameSize, and EcalHitResponse::setPhaseShift().

57  :
58  DigiAccumulatorMixMod(),
59  m_APDShape ( params.getParameter<double>( "apdShapeTstart" ) ,
60  params.getParameter<double>( "apdShapeTau" ) ) ,
61  m_EBShape ( ) ,
62  m_EEShape ( ) ,
63  m_ESShape ( ) ,
64  m_EBdigiCollection ( params.getParameter<std::string>("EBdigiCollection") ) ,
65  m_EEdigiCollection ( params.getParameter<std::string>("EEdigiCollection") ) ,
66  m_ESdigiCollection ( params.getParameter<std::string>("ESdigiCollection") ) ,
67  m_hitsProducerTag ( params.getParameter<std::string>("hitsProducer" ) ) ,
68  m_useLCcorrection ( params.getUntrackedParameter<bool>("UseLCcorrection") ) ,
69  m_apdSeparateDigi ( params.getParameter<bool> ("apdSeparateDigi") ) ,
70 
71  m_EBs25notCont ( params.getParameter<double> ("EBs25notContainment") ) ,
72  m_EEs25notCont ( params.getParameter<double> ("EEs25notContainment") ) ,
73 
74  m_readoutFrameSize ( params.getParameter<int> ("readoutFrameSize") ) ,
75  m_ParameterMap ( new EcalSimParameterMap(
76  params.getParameter<double> ("simHitToPhotoelectronsBarrel") ,
77  params.getParameter<double> ("simHitToPhotoelectronsEndcap") ,
78  params.getParameter<double> ("photoelectronsToAnalogBarrel") ,
79  params.getParameter<double> ("photoelectronsToAnalogEndcap") ,
80  params.getParameter<double> ("samplingFactor") ,
81  params.getParameter<double> ("timePhase") ,
82  m_readoutFrameSize ,
83  params.getParameter<int> ("binOfMaximum") ,
84  params.getParameter<bool> ("doPhotostatistics") ,
85  params.getParameter<bool> ("syncPhase") ) ) ,
86 
87  m_apdDigiTag ( params.getParameter<std::string> ("apdDigiTag" ) ) ,
88  m_apdParameters ( new APDSimParameters(
89  params.getParameter<bool> ("apdAddToBarrel" ) ,
90  m_apdSeparateDigi ,
91  params.getParameter<double> ("apdSimToPELow" ) ,
92  params.getParameter<double> ("apdSimToPEHigh" ) ,
93  params.getParameter<double> ("apdTimeOffset" ) ,
94  params.getParameter<double> ("apdTimeOffWidth" ) ,
95  params.getParameter<bool> ("apdDoPEStats" ) ,
96  m_apdDigiTag ,
97  params.getParameter<std::vector<double> > ( "apdNonlParms" ) ) ) ,
98 
99  m_APDResponse ( !m_apdSeparateDigi ? 0 :
100  new EBHitResponse( m_ParameterMap ,
101  &m_EBShape ,
102  true ,
103  m_apdParameters ,
104  &m_APDShape ) ) ,
105 
106  m_EBResponse ( new EBHitResponse( m_ParameterMap ,
107  &m_EBShape ,
108  false , // barrel
109  m_apdParameters ,
110  &m_APDShape ) ) ,
111 
112  m_EEResponse ( new EEHitResponse( m_ParameterMap,
113  &m_EEShape ) ) ,
114  m_ESResponse ( new ESHitResponse( m_ParameterMap, &m_ESShape ) ) ,
115  m_ESOldResponse ( new CaloHitResponse( m_ParameterMap, &m_ESShape ) ) ,
116 
117  m_addESNoise ( params.getParameter<bool> ("doESNoise") ) ,
118  m_PreMix1 ( params.getParameter<bool> ("EcalPreMixStage1") ) ,
119  m_PreMix2 ( params.getParameter<bool> ("EcalPreMixStage2") ) ,
120 
121  m_doFastES ( params.getParameter<bool> ("doFast" ) ) ,
122 
123  m_ESElectronicsSim ( m_doFastES ? 0 :
124  new ESElectronicsSim( m_addESNoise ) ) ,
125 
126  m_ESOldDigitizer ( m_doFastES ? 0 :
127  new ESOldDigitizer( m_ESOldResponse ,
128  m_ESElectronicsSim ,
129  m_addESNoise ) ) ,
130 
131  m_ESElectronicsSimFast ( !m_doFastES ? 0 :
132  new ESElectronicsSimFast( m_addESNoise ,
133  m_PreMix1) ) ,
134 
135  m_ESDigitizer ( !m_doFastES ? 0 :
136  new ESDigitizer( m_ESResponse ,
137  m_ESElectronicsSimFast ,
138  m_addESNoise ) ) ,
139 
140  m_APDDigitizer ( 0 ) ,
141  m_BarrelDigitizer ( 0 ) ,
142  m_EndcapDigitizer ( 0 ) ,
143  m_ElectronicsSim ( 0 ) ,
144  m_Coder ( 0 ) ,
145  m_APDElectronicsSim ( 0 ) ,
146  m_APDCoder ( 0 ) ,
147  m_Geometry ( 0 ) ,
148  m_EBCorrNoise ( ) ,
149  m_EECorrNoise ( )
150 {
151  if(m_apdSeparateDigi) mixMod.produces<EBDigiCollection>(m_apdDigiTag);
152  mixMod.produces<EBDigiCollection>(m_EBdigiCollection);
153  mixMod.produces<EEDigiCollection>(m_EEdigiCollection);
154  mixMod.produces<ESDigiCollection>(m_ESdigiCollection);
155  iC.consumes<std::vector<PCaloHit> >(edm::InputTag(m_hitsProducerTag, "EcalHitsEB"));
156  iC.consumes<std::vector<PCaloHit> >(edm::InputTag(m_hitsProducerTag, "EcalHitsEE"));
157  iC.consumes<std::vector<PCaloHit> >(edm::InputTag(m_hitsProducerTag, "EcalHitsES"));
158 
159  const std::vector<double> ebCorMatG12 = params.getParameter< std::vector<double> >("EBCorrNoiseMatrixG12");
160  const std::vector<double> eeCorMatG12 = params.getParameter< std::vector<double> >("EECorrNoiseMatrixG12");
161  const std::vector<double> ebCorMatG06 = params.getParameter< std::vector<double> >("EBCorrNoiseMatrixG06");
162  const std::vector<double> eeCorMatG06 = params.getParameter< std::vector<double> >("EECorrNoiseMatrixG06");
163  const std::vector<double> ebCorMatG01 = params.getParameter< std::vector<double> >("EBCorrNoiseMatrixG01");
164  const std::vector<double> eeCorMatG01 = params.getParameter< std::vector<double> >("EECorrNoiseMatrixG01");
165 
166  const bool applyConstantTerm = params.getParameter<bool> ("applyConstantTerm");
167  const double rmsConstantTerm = params.getParameter<double> ("ConstantTerm");
168 
169  const bool addNoise = params.getParameter<bool> ("doENoise");
170  const bool cosmicsPhase = params.getParameter<bool> ("cosmicsPhase");
171  const double cosmicsShift = params.getParameter<double> ("cosmicsShift");
172 
173 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
174 
175  // further phase for cosmics studies
176  if( cosmicsPhase )
177  {
178  m_EBResponse->setPhaseShift( 1. + cosmicsShift ) ;
179  m_EEResponse->setPhaseShift( 1. + cosmicsShift ) ;
180  }
181 
182  EcalCorrMatrix ebMatrix[ 3 ] ;
183  EcalCorrMatrix eeMatrix[ 3 ] ;
184 
185  assert( ebCorMatG12.size() == m_readoutFrameSize ) ;
186  assert( eeCorMatG12.size() == m_readoutFrameSize ) ;
187  assert( ebCorMatG06.size() == m_readoutFrameSize ) ;
188  assert( eeCorMatG06.size() == m_readoutFrameSize ) ;
189  assert( ebCorMatG01.size() == m_readoutFrameSize ) ;
190  assert( eeCorMatG01.size() == m_readoutFrameSize ) ;
191 
192  assert( 1.e-7 > fabs( ebCorMatG12[0] - 1.0 ) ) ;
193  assert( 1.e-7 > fabs( ebCorMatG06[0] - 1.0 ) ) ;
194  assert( 1.e-7 > fabs( ebCorMatG01[0] - 1.0 ) ) ;
195  assert( 1.e-7 > fabs( eeCorMatG12[0] - 1.0 ) ) ;
196  assert( 1.e-7 > fabs( eeCorMatG06[0] - 1.0 ) ) ;
197  assert( 1.e-7 > fabs( eeCorMatG01[0] - 1.0 ) ) ;
198 
199  for ( unsigned int row ( 0 ) ; row != m_readoutFrameSize ; ++row )
200  {
201  assert( 0 == row || 1. >= ebCorMatG12[row] ) ;
202  assert( 0 == row || 1. >= ebCorMatG06[row] ) ;
203  assert( 0 == row || 1. >= ebCorMatG01[row] ) ;
204  assert( 0 == row || 1. >= eeCorMatG12[row] ) ;
205  assert( 0 == row || 1. >= eeCorMatG06[row] ) ;
206  assert( 0 == row || 1. >= eeCorMatG01[row] ) ;
207  for ( unsigned int column ( 0 ) ; column <= row ; ++column )
208  {
209  const unsigned int index ( row - column ) ;
210  ebMatrix[0]( row, column ) = ebCorMatG12[ index ] ;
211  eeMatrix[0]( row, column ) = eeCorMatG12[ index ] ;
212  ebMatrix[1]( row, column ) = ebCorMatG06[ index ] ;
213  eeMatrix[1]( row, column ) = eeCorMatG06[ index ] ;
214  ebMatrix[2]( row, column ) = ebCorMatG01[ index ] ;
215  eeMatrix[2]( row, column ) = eeCorMatG01[ index ] ;
216  }
217  }
218 
219  m_EBCorrNoise[0] = new CorrelatedNoisifier<EcalCorrMatrix>( ebMatrix[0] ) ;
220  m_EECorrNoise[0] = new CorrelatedNoisifier<EcalCorrMatrix>( eeMatrix[0] ) ;
221  m_EBCorrNoise[1] = new CorrelatedNoisifier<EcalCorrMatrix>( ebMatrix[1] ) ;
222  m_EECorrNoise[1] = new CorrelatedNoisifier<EcalCorrMatrix>( eeMatrix[1] ) ;
223  m_EBCorrNoise[2] = new CorrelatedNoisifier<EcalCorrMatrix>( ebMatrix[2] ) ;
224  m_EECorrNoise[2] = new CorrelatedNoisifier<EcalCorrMatrix>( eeMatrix[2] ) ;
225 
226  m_Coder = new EcalCoder( addNoise ,
227  m_PreMix1 ,
228  m_EBCorrNoise[0] ,
229  m_EECorrNoise[0] ,
230  m_EBCorrNoise[1] ,
231  m_EECorrNoise[1] ,
232  m_EBCorrNoise[2] ,
233  m_EECorrNoise[2] ) ;
234 
235  m_ElectronicsSim = new EcalElectronicsSim( m_ParameterMap ,
236  m_Coder ,
237  applyConstantTerm ,
238  rmsConstantTerm ) ;
239 
240  if( m_apdSeparateDigi )
241  {
242  m_APDCoder = new EcalCoder( false ,
243  m_EBCorrNoise[0] ,
244  m_EECorrNoise[0] ,
245  m_EBCorrNoise[1] ,
246  m_EECorrNoise[1] ,
247  m_EBCorrNoise[2] ,
248  m_EECorrNoise[2] ) ;
249 
250  m_APDElectronicsSim = new EcalElectronicsSim( m_ParameterMap ,
251  m_APDCoder ,
252  applyConstantTerm ,
253  rmsConstantTerm ) ;
254 
255  m_APDDigitizer = new EBDigitizer( m_APDResponse ,
256  m_APDElectronicsSim ,
257  false ) ;
258  }
259 
260  m_BarrelDigitizer = new EBDigitizer( m_EBResponse ,
261  m_ElectronicsSim ,
262  addNoise ) ;
263 
264  m_EndcapDigitizer = new EEDigitizer( m_EEResponse ,
265  m_ElectronicsSim ,
266  addNoise ) ;
267 }
EcalDigiProducer::m_EndcapDigitizer
EEDigitizer * m_EndcapDigitizer
Definition: EcalDigiProducer.h:138
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:41
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
EcalDigiProducer::m_ESOldResponse
CaloHitResponse * m_ESOldResponse
Definition: EcalDigiProducer.h:123
ESOldDigitizer
CaloTDigitizer< ESOldDigitizerTraits > ESOldDigitizer
Definition: EcalDigiProducer.h:23
EEDigitizer
EcalTDigitizer< EEDigitizerTraits > EEDigitizer
Definition: EcalDigiProducer.h:22
EcalDigiProducer::m_EBCorrNoise
CorrelatedNoisifier< EcalCorrMatrix > * m_EBCorrNoise[3]
Definition: EcalDigiProducer.h:148
EcalDigiProducer::m_APDShape
const APDShape m_APDShape
Definition: EcalDigiProducer.h:93
EcalDigiProducer::m_ESdigiCollection
const std::string m_ESdigiCollection
Definition: EcalDigiProducer.h:100
EcalDigiProducer::m_ParameterMap
const EcalSimParameterMap * m_ParameterMap
Definition: EcalDigiProducer.h:112
EcalDigiProducer::m_hitsProducerTag
const std::string m_hitsProducerTag
Definition: EcalDigiProducer.h:101
EcalDigiProducer::m_ESOldDigitizer
ESOldDigitizer * m_ESOldDigitizer
Definition: EcalDigiProducer.h:132
EcalDigiProducer::m_EEdigiCollection
const std::string m_EEdigiCollection
Definition: EcalDigiProducer.h:99
EcalDigiProducer::m_ESElectronicsSimFast
ESElectronicsSimFast * m_ESElectronicsSimFast
Definition: EcalDigiProducer.h:133
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
EcalDigiProducer::m_ESDigitizer
ESDigitizer * m_ESDigitizer
Definition: EcalDigiProducer.h:134
EcalDigiProducer::m_apdParameters
const APDSimParameters * m_apdParameters
Definition: EcalDigiProducer.h:115
EcalDigiProducer::m_addESNoise
const bool m_addESNoise
Definition: EcalDigiProducer.h:125
EcalDigiProducer::m_EECorrNoise
CorrelatedNoisifier< EcalCorrMatrix > * m_EECorrNoise[3]
Definition: EcalDigiProducer.h:149
EcalDigiProducer::m_readoutFrameSize
const unsigned int m_readoutFrameSize
Definition: EcalDigiProducer.h:110
EcalDigiProducer::m_ElectronicsSim
EcalElectronicsSim * m_ElectronicsSim
Definition: EcalDigiProducer.h:140
EcalDigiProducer::m_APDResponse
EBHitResponse * m_APDResponse
Definition: EcalDigiProducer.h:117
CaloHitResponse
Creates electronics signals from hits.
Definition: CaloHitResponse.h:34
EcalDigiProducer::m_EBs25notCont
const double m_EBs25notCont
Definition: EcalDigiProducer.h:107
EcalDigiProducer::m_EBResponse
EBHitResponse * m_EBResponse
Definition: EcalDigiProducer.h:119
EcalHitResponse::setPhaseShift
void setPhaseShift(double phaseShift)
Definition: EcalHitResponse.cc:79
EcalDigiProducer::m_ESElectronicsSim
ESElectronicsSim * m_ESElectronicsSim
Definition: EcalDigiProducer.h:131
EcalDigiProducer::m_BarrelDigitizer
EBDigitizer * m_BarrelDigitizer
Definition: EcalDigiProducer.h:137
EcalCorrMatrix
math::ErrorD< CaloSamples::MAXSAMPLES >::type EcalCorrMatrix
Definition: EcalCorrelatedNoiseMatrix.h:7
EcalDigiProducer::m_Geometry
const CaloGeometry * m_Geometry
Definition: EcalDigiProducer.h:146
EcalDigiProducer::m_APDCoder
EcalCoder * m_APDCoder
Definition: EcalDigiProducer.h:144
EcalDigiProducer::m_ESResponse
ESHitResponse * m_ESResponse
Definition: EcalDigiProducer.h:122
column
#define column(...)
Definition: DbCore.h:74
EBDigitizer
EcalTDigitizer< EBDigitizerTraits > EBDigitizer
Definition: EcalDigiProducer.h:21
EcalDigiProducer::m_APDDigitizer
EBDigitizer * m_APDDigitizer
Definition: EcalDigiProducer.h:136
EcalDigiProducer::m_apdSeparateDigi
const bool m_apdSeparateDigi
Definition: EcalDigiProducer.h:105
EcalDigiProducer::m_apdDigiTag
const std::string m_apdDigiTag
Definition: EcalDigiProducer.h:114
EcalDigiProducer::m_EEShape
const EEShape m_EEShape
Definition: EcalDigiProducer.h:95
EcalDigiProducer::m_EBShape
const EBShape m_EBShape
Definition: EcalDigiProducer.h:94
EcalDigiProducer::m_EBdigiCollection
const std::string m_EBdigiCollection
Definition: EcalDigiProducer.h:98
EcalDigiProducer::m_APDElectronicsSim
EcalElectronicsSim * m_APDElectronicsSim
Definition: EcalDigiProducer.h:143
EcalDigiProducer::m_EEs25notCont
const double m_EEs25notCont
Definition: EcalDigiProducer.h:108
EcalDigiProducer::m_EEResponse
EEHitResponse * m_EEResponse
Definition: EcalDigiProducer.h:120
EcalDigiProducer::EcalDigiProducer ( const edm::ParameterSet &  params,
edm::ConsumesCollector &  iC 
)

Definition at line 270 of file EcalDigiProducer.cc.

References column, edm::ConsumesCollector::consumes(), alignCSCRings::e, edm::ParameterSet::getParameter(), cmsHarvester::index, m_APDCoder, m_APDDigitizer, m_APDElectronicsSim, m_APDResponse, m_apdSeparateDigi, m_BarrelDigitizer, m_Coder, m_EBCorrNoise, m_EBResponse, m_EECorrNoise, m_EEResponse, m_ElectronicsSim, m_EndcapDigitizer, m_hitsProducerTag, m_ParameterMap, m_PreMix1, m_readoutFrameSize, and EcalHitResponse::setPhaseShift().

270  :
271  DigiAccumulatorMixMod(),
272  m_APDShape ( params.getParameter<double>( "apdShapeTstart" ) ,
273  params.getParameter<double>( "apdShapeTau" ) ) ,
274  m_EBShape ( ) ,
275  m_EEShape ( ) ,
276  m_ESShape ( ) ,
277  m_EBdigiCollection ( params.getParameter<std::string>("EBdigiCollection") ) ,
278  m_EEdigiCollection ( params.getParameter<std::string>("EEdigiCollection") ) ,
279  m_ESdigiCollection ( params.getParameter<std::string>("ESdigiCollection") ) ,
280  m_hitsProducerTag ( params.getParameter<std::string>("hitsProducer" ) ) ,
281  m_useLCcorrection ( params.getUntrackedParameter<bool>("UseLCcorrection") ) ,
282  m_apdSeparateDigi ( params.getParameter<bool> ("apdSeparateDigi") ) ,
283 
284  m_EBs25notCont ( params.getParameter<double> ("EBs25notContainment") ) ,
285  m_EEs25notCont ( params.getParameter<double> ("EEs25notContainment") ) ,
286 
287  m_readoutFrameSize ( params.getParameter<int> ("readoutFrameSize") ) ,
288  m_ParameterMap ( new EcalSimParameterMap(
289  params.getParameter<double> ("simHitToPhotoelectronsBarrel") ,
290  params.getParameter<double> ("simHitToPhotoelectronsEndcap") ,
291  params.getParameter<double> ("photoelectronsToAnalogBarrel") ,
292  params.getParameter<double> ("photoelectronsToAnalogEndcap") ,
293  params.getParameter<double> ("samplingFactor") ,
294  params.getParameter<double> ("timePhase") ,
295  m_readoutFrameSize ,
296  params.getParameter<int> ("binOfMaximum") ,
297  params.getParameter<bool> ("doPhotostatistics") ,
298  params.getParameter<bool> ("syncPhase") ) ) ,
299 
300  m_apdDigiTag ( params.getParameter<std::string> ("apdDigiTag" ) ) ,
301  m_apdParameters ( new APDSimParameters(
302  params.getParameter<bool> ("apdAddToBarrel" ) ,
303  m_apdSeparateDigi ,
304  params.getParameter<double> ("apdSimToPELow" ) ,
305  params.getParameter<double> ("apdSimToPEHigh" ) ,
306  params.getParameter<double> ("apdTimeOffset" ) ,
307  params.getParameter<double> ("apdTimeOffWidth" ) ,
308  params.getParameter<bool> ("apdDoPEStats" ) ,
309  m_apdDigiTag ,
310  params.getParameter<std::vector<double> > ( "apdNonlParms" ) ) ) ,
311 
312  m_APDResponse ( !m_apdSeparateDigi ? 0 :
313  new EBHitResponse( m_ParameterMap ,
314  &m_EBShape ,
315  true ,
316  m_apdParameters ,
317  &m_APDShape ) ) ,
318 
319  m_EBResponse ( new EBHitResponse( m_ParameterMap ,
320  &m_EBShape ,
321  false , // barrel
322  m_apdParameters ,
323  &m_APDShape ) ) ,
324 
325  m_EEResponse ( new EEHitResponse( m_ParameterMap,
326  &m_EEShape ) ) ,
327  m_ESResponse ( new ESHitResponse( m_ParameterMap, &m_ESShape ) ) ,
328  m_ESOldResponse ( new CaloHitResponse( m_ParameterMap, &m_ESShape ) ) ,
329 
330  m_addESNoise ( params.getParameter<bool> ("doESNoise") ) ,
331  m_PreMix1 ( params.getParameter<bool> ("EcalPreMixStage1") ) ,
332  m_PreMix2 ( params.getParameter<bool> ("EcalPreMixStage2") ) ,
333 
334  m_doFastES ( params.getParameter<bool> ("doFast" ) ) ,
335 
336  m_ESElectronicsSim ( m_doFastES ? 0 :
337  new ESElectronicsSim( m_addESNoise ) ) ,
338 
339  m_ESOldDigitizer ( m_doFastES ? 0 :
340  new ESOldDigitizer( m_ESOldResponse ,
341  m_ESElectronicsSim ,
342  m_addESNoise ) ) ,
343 
344  m_ESElectronicsSimFast ( !m_doFastES ? 0 :
345  new ESElectronicsSimFast( m_addESNoise,
346  m_PreMix1 ) ) ,
347 
348  m_ESDigitizer ( !m_doFastES ? 0 :
349  new ESDigitizer( m_ESResponse ,
350  m_ESElectronicsSimFast ,
351  m_addESNoise ) ) ,
352 
353  m_APDDigitizer ( 0 ) ,
354  m_BarrelDigitizer ( 0 ) ,
355  m_EndcapDigitizer ( 0 ) ,
356  m_ElectronicsSim ( 0 ) ,
357  m_Coder ( 0 ) ,
358  m_APDElectronicsSim ( 0 ) ,
359  m_APDCoder ( 0 ) ,
360  m_Geometry ( 0 ) ,
361  m_EBCorrNoise ( ) ,
362  m_EECorrNoise ( )
363 {
364  // "produces" statements taken care of elsewhere.
365  // if(m_apdSeparateDigi) mixMod.produces<EBDigiCollection>(m_apdDigiTag);
366  // mixMod.produces<EBDigiCollection>(m_EBdigiCollection);
367  // mixMod.produces<EEDigiCollection>(m_EEdigiCollection);
368  // mixMod.produces<ESDigiCollection>(m_ESdigiCollection);
369  iC.consumes<std::vector<PCaloHit> >(edm::InputTag(m_hitsProducerTag, "EcalHitsEB"));
370  iC.consumes<std::vector<PCaloHit> >(edm::InputTag(m_hitsProducerTag, "EcalHitsEE"));
371  iC.consumes<std::vector<PCaloHit> >(edm::InputTag(m_hitsProducerTag, "EcalHitsES"));
372 
373  const std::vector<double> ebCorMatG12 = params.getParameter< std::vector<double> >("EBCorrNoiseMatrixG12");
374  const std::vector<double> eeCorMatG12 = params.getParameter< std::vector<double> >("EECorrNoiseMatrixG12");
375  const std::vector<double> ebCorMatG06 = params.getParameter< std::vector<double> >("EBCorrNoiseMatrixG06");
376  const std::vector<double> eeCorMatG06 = params.getParameter< std::vector<double> >("EECorrNoiseMatrixG06");
377  const std::vector<double> ebCorMatG01 = params.getParameter< std::vector<double> >("EBCorrNoiseMatrixG01");
378  const std::vector<double> eeCorMatG01 = params.getParameter< std::vector<double> >("EECorrNoiseMatrixG01");
379 
380  const bool applyConstantTerm = params.getParameter<bool> ("applyConstantTerm");
381  const double rmsConstantTerm = params.getParameter<double> ("ConstantTerm");
382 
383  const bool addNoise = params.getParameter<bool> ("doNoise");
384  const bool cosmicsPhase = params.getParameter<bool> ("cosmicsPhase");
385  const double cosmicsShift = params.getParameter<double> ("cosmicsShift");
386 
387 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
388 
389  // further phase for cosmics studies
390  if( cosmicsPhase )
391  {
392  m_EBResponse->setPhaseShift( 1. + cosmicsShift ) ;
393  m_EEResponse->setPhaseShift( 1. + cosmicsShift ) ;
394  }
395 
396  EcalCorrMatrix ebMatrix[ 3 ] ;
397  EcalCorrMatrix eeMatrix[ 3 ] ;
398 
399  assert( ebCorMatG12.size() == m_readoutFrameSize ) ;
400  assert( eeCorMatG12.size() == m_readoutFrameSize ) ;
401  assert( ebCorMatG06.size() == m_readoutFrameSize ) ;
402  assert( eeCorMatG06.size() == m_readoutFrameSize ) ;
403  assert( ebCorMatG01.size() == m_readoutFrameSize ) ;
404  assert( eeCorMatG01.size() == m_readoutFrameSize ) ;
405 
406  assert( 1.e-7 > fabs( ebCorMatG12[0] - 1.0 ) ) ;
407  assert( 1.e-7 > fabs( ebCorMatG06[0] - 1.0 ) ) ;
408  assert( 1.e-7 > fabs( ebCorMatG01[0] - 1.0 ) ) ;
409  assert( 1.e-7 > fabs( eeCorMatG12[0] - 1.0 ) ) ;
410  assert( 1.e-7 > fabs( eeCorMatG06[0] - 1.0 ) ) ;
411  assert( 1.e-7 > fabs( eeCorMatG01[0] - 1.0 ) ) ;
412 
413  for ( unsigned int row ( 0 ) ; row != m_readoutFrameSize ; ++row )
414  {
415  assert( 0 == row || 1. >= ebCorMatG12[row] ) ;
416  assert( 0 == row || 1. >= ebCorMatG06[row] ) ;
417  assert( 0 == row || 1. >= ebCorMatG01[row] ) ;
418  assert( 0 == row || 1. >= eeCorMatG12[row] ) ;
419  assert( 0 == row || 1. >= eeCorMatG06[row] ) ;
420  assert( 0 == row || 1. >= eeCorMatG01[row] ) ;
421  for ( unsigned int column ( 0 ) ; column <= row ; ++column )
422  {
423  const unsigned int index ( row - column ) ;
424  ebMatrix[0]( row, column ) = ebCorMatG12[ index ] ;
425  eeMatrix[0]( row, column ) = eeCorMatG12[ index ] ;
426  ebMatrix[1]( row, column ) = ebCorMatG06[ index ] ;
427  eeMatrix[1]( row, column ) = eeCorMatG06[ index ] ;
428  ebMatrix[2]( row, column ) = ebCorMatG01[ index ] ;
429  eeMatrix[2]( row, column ) = eeCorMatG01[ index ] ;
430  }
431  }
432 
433  m_EBCorrNoise[0] = new CorrelatedNoisifier<EcalCorrMatrix>( ebMatrix[0] ) ;
434  m_EECorrNoise[0] = new CorrelatedNoisifier<EcalCorrMatrix>( eeMatrix[0] ) ;
435  m_EBCorrNoise[1] = new CorrelatedNoisifier<EcalCorrMatrix>( ebMatrix[1] ) ;
436  m_EECorrNoise[1] = new CorrelatedNoisifier<EcalCorrMatrix>( eeMatrix[1] ) ;
437  m_EBCorrNoise[2] = new CorrelatedNoisifier<EcalCorrMatrix>( ebMatrix[2] ) ;
438  m_EECorrNoise[2] = new CorrelatedNoisifier<EcalCorrMatrix>( eeMatrix[2] ) ;
439 
440  m_Coder = new EcalCoder( addNoise ,
441  m_PreMix1 ,
442  m_EBCorrNoise[0] ,
443  m_EECorrNoise[0] ,
444  m_EBCorrNoise[1] ,
445  m_EECorrNoise[1] ,
446  m_EBCorrNoise[2] ,
447  m_EECorrNoise[2] ) ;
448 
449  m_ElectronicsSim = new EcalElectronicsSim( m_ParameterMap ,
450  m_Coder ,
451  applyConstantTerm ,
452  rmsConstantTerm ) ;
453 
454  if( m_apdSeparateDigi )
455  {
456  m_APDCoder = new EcalCoder( false ,
457  m_EBCorrNoise[0] ,
458  m_EECorrNoise[0] ,
459  m_EBCorrNoise[1] ,
460  m_EECorrNoise[1] ,
461  m_EBCorrNoise[2] ,
462  m_EECorrNoise[2] ) ;
463 
464  m_APDElectronicsSim = new EcalElectronicsSim( m_ParameterMap ,
465  m_APDCoder ,
466  applyConstantTerm ,
467  rmsConstantTerm ) ;
468 
469  m_APDDigitizer = new EBDigitizer( m_APDResponse ,
470  m_APDElectronicsSim ,
471  false ) ;
472  }
473 
474  m_BarrelDigitizer = new EBDigitizer( m_EBResponse ,
475  m_ElectronicsSim ,
476  addNoise ) ;
477 
478  m_EndcapDigitizer = new EEDigitizer( m_EEResponse ,
479  m_ElectronicsSim ,
480  addNoise ) ;
481 }
EcalDigiProducer::m_EndcapDigitizer
EEDigitizer * m_EndcapDigitizer
Definition: EcalDigiProducer.h:138
edm::ConsumesCollector::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: ConsumesCollector.h:41
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
EcalDigiProducer::m_ESOldResponse
CaloHitResponse * m_ESOldResponse
Definition: EcalDigiProducer.h:123
ESOldDigitizer
CaloTDigitizer< ESOldDigitizerTraits > ESOldDigitizer
Definition: EcalDigiProducer.h:23
EEDigitizer
EcalTDigitizer< EEDigitizerTraits > EEDigitizer
Definition: EcalDigiProducer.h:22
EcalDigiProducer::m_EBCorrNoise
CorrelatedNoisifier< EcalCorrMatrix > * m_EBCorrNoise[3]
Definition: EcalDigiProducer.h:148
EcalDigiProducer::m_APDShape
const APDShape m_APDShape
Definition: EcalDigiProducer.h:93
EcalDigiProducer::m_ESdigiCollection
const std::string m_ESdigiCollection
Definition: EcalDigiProducer.h:100
EcalDigiProducer::m_ParameterMap
const EcalSimParameterMap * m_ParameterMap
Definition: EcalDigiProducer.h:112
EcalDigiProducer::m_hitsProducerTag
const std::string m_hitsProducerTag
Definition: EcalDigiProducer.h:101
EcalDigiProducer::m_ESOldDigitizer
ESOldDigitizer * m_ESOldDigitizer
Definition: EcalDigiProducer.h:132
EcalDigiProducer::m_EEdigiCollection
const std::string m_EEdigiCollection
Definition: EcalDigiProducer.h:99
EcalDigiProducer::m_ESElectronicsSimFast
ESElectronicsSimFast * m_ESElectronicsSimFast
Definition: EcalDigiProducer.h:133
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
EcalDigiProducer::m_ESDigitizer
ESDigitizer * m_ESDigitizer
Definition: EcalDigiProducer.h:134
EcalDigiProducer::m_apdParameters
const APDSimParameters * m_apdParameters
Definition: EcalDigiProducer.h:115
EcalDigiProducer::m_addESNoise
const bool m_addESNoise
Definition: EcalDigiProducer.h:125
EcalDigiProducer::m_EECorrNoise
CorrelatedNoisifier< EcalCorrMatrix > * m_EECorrNoise[3]
Definition: EcalDigiProducer.h:149
EcalDigiProducer::m_readoutFrameSize
const unsigned int m_readoutFrameSize
Definition: EcalDigiProducer.h:110
EcalDigiProducer::m_ElectronicsSim
EcalElectronicsSim * m_ElectronicsSim
Definition: EcalDigiProducer.h:140
EcalDigiProducer::m_APDResponse
EBHitResponse * m_APDResponse
Definition: EcalDigiProducer.h:117
CaloHitResponse
Creates electronics signals from hits.
Definition: CaloHitResponse.h:34
EcalDigiProducer::m_EBs25notCont
const double m_EBs25notCont
Definition: EcalDigiProducer.h:107
EcalDigiProducer::m_EBResponse
EBHitResponse * m_EBResponse
Definition: EcalDigiProducer.h:119
EcalHitResponse::setPhaseShift
void setPhaseShift(double phaseShift)
Definition: EcalHitResponse.cc:79
EcalDigiProducer::m_ESElectronicsSim
ESElectronicsSim * m_ESElectronicsSim
Definition: EcalDigiProducer.h:131
EcalDigiProducer::m_BarrelDigitizer
EBDigitizer * m_BarrelDigitizer
Definition: EcalDigiProducer.h:137
EcalCorrMatrix
math::ErrorD< CaloSamples::MAXSAMPLES >::type EcalCorrMatrix
Definition: EcalCorrelatedNoiseMatrix.h:7
EcalDigiProducer::m_Geometry
const CaloGeometry * m_Geometry
Definition: EcalDigiProducer.h:146
EcalDigiProducer::m_APDCoder
EcalCoder * m_APDCoder
Definition: EcalDigiProducer.h:144
EcalDigiProducer::m_ESResponse
ESHitResponse * m_ESResponse
Definition: EcalDigiProducer.h:122
column
#define column(...)
Definition: DbCore.h:74
EBDigitizer
EcalTDigitizer< EBDigitizerTraits > EBDigitizer
Definition: EcalDigiProducer.h:21
EcalDigiProducer::m_APDDigitizer
EBDigitizer * m_APDDigitizer
Definition: EcalDigiProducer.h:136
EcalDigiProducer::m_apdSeparateDigi
const bool m_apdSeparateDigi
Definition: EcalDigiProducer.h:105
EcalDigiProducer::m_apdDigiTag
const std::string m_apdDigiTag
Definition: EcalDigiProducer.h:114
EcalDigiProducer::m_EEShape
const EEShape m_EEShape
Definition: EcalDigiProducer.h:95
EcalDigiProducer::m_EBShape
const EBShape m_EBShape
Definition: EcalDigiProducer.h:94
EcalDigiProducer::m_EBdigiCollection
const std::string m_EBdigiCollection
Definition: EcalDigiProducer.h:98
EcalDigiProducer::m_APDElectronicsSim
EcalElectronicsSim * m_APDElectronicsSim
Definition: EcalDigiProducer.h:143
EcalDigiProducer::m_EEs25notCont
const double m_EEs25notCont
Definition: EcalDigiProducer.h:108
EcalDigiProducer::m_EEResponse
EEHitResponse * m_EEResponse
Definition: EcalDigiProducer.h:120
EcalDigiProducer::~EcalDigiProducer ( )
virtual

Definition at line 484 of file EcalDigiProducer.cc.

References m_APDCoder, m_APDDigitizer, m_APDElectronicsSim, m_apdParameters, m_APDResponse, m_BarrelDigitizer, m_Coder, m_EBCorrNoise, m_EBResponse, m_EECorrNoise, m_EEResponse, m_ElectronicsSim, m_EndcapDigitizer, m_ESDigitizer, m_ESElectronicsSim, m_ESElectronicsSimFast, m_ESOldDigitizer, m_ESOldResponse, m_ESResponse, and m_ParameterMap.

485 {
486  delete m_EndcapDigitizer ;
487  delete m_BarrelDigitizer ;
488  delete m_APDDigitizer ;
489  delete m_APDElectronicsSim ;
490  delete m_APDCoder ;
491  delete m_ElectronicsSim ;
492  delete m_Coder ;
493  delete m_EBCorrNoise[0] ;
494  delete m_EECorrNoise[0] ;
495  delete m_EBCorrNoise[1] ;
496  delete m_EECorrNoise[1] ;
497  delete m_EBCorrNoise[2] ;
498  delete m_EECorrNoise[2] ;
499 
500  delete m_ESDigitizer ;
501  delete m_ESElectronicsSimFast ;
502  delete m_ESOldDigitizer ;
503  delete m_ESElectronicsSim ;
504 
505  delete m_ESOldResponse ;
506  delete m_ESResponse ;
507  delete m_EEResponse ;
508  delete m_EBResponse ;
509  delete m_APDResponse ;
510 
511  delete m_apdParameters ;
512  delete m_ParameterMap ;
513 }
EcalDigiProducer::m_EndcapDigitizer
EEDigitizer * m_EndcapDigitizer
Definition: EcalDigiProducer.h:138
EcalDigiProducer::m_ESOldResponse
CaloHitResponse * m_ESOldResponse
Definition: EcalDigiProducer.h:123
EcalDigiProducer::m_EBCorrNoise
CorrelatedNoisifier< EcalCorrMatrix > * m_EBCorrNoise[3]
Definition: EcalDigiProducer.h:148
EcalDigiProducer::m_ParameterMap
const EcalSimParameterMap * m_ParameterMap
Definition: EcalDigiProducer.h:112
EcalDigiProducer::m_ESOldDigitizer
ESOldDigitizer * m_ESOldDigitizer
Definition: EcalDigiProducer.h:132
EcalDigiProducer::m_ESElectronicsSimFast
ESElectronicsSimFast * m_ESElectronicsSimFast
Definition: EcalDigiProducer.h:133
EcalDigiProducer::m_ESDigitizer
ESDigitizer * m_ESDigitizer
Definition: EcalDigiProducer.h:134
EcalDigiProducer::m_apdParameters
const APDSimParameters * m_apdParameters
Definition: EcalDigiProducer.h:115
EcalDigiProducer::m_EECorrNoise
CorrelatedNoisifier< EcalCorrMatrix > * m_EECorrNoise[3]
Definition: EcalDigiProducer.h:149
EcalDigiProducer::m_ElectronicsSim
EcalElectronicsSim * m_ElectronicsSim
Definition: EcalDigiProducer.h:140
EcalDigiProducer::m_APDResponse
EBHitResponse * m_APDResponse
Definition: EcalDigiProducer.h:117
EcalDigiProducer::m_EBResponse
EBHitResponse * m_EBResponse
Definition: EcalDigiProducer.h:119
EcalDigiProducer::m_ESElectronicsSim
ESElectronicsSim * m_ESElectronicsSim
Definition: EcalDigiProducer.h:131
EcalDigiProducer::m_BarrelDigitizer
EBDigitizer * m_BarrelDigitizer
Definition: EcalDigiProducer.h:137
EcalDigiProducer::m_APDCoder
EcalCoder * m_APDCoder
Definition: EcalDigiProducer.h:144
EcalDigiProducer::m_ESResponse
ESHitResponse * m_ESResponse
Definition: EcalDigiProducer.h:122
EcalDigiProducer::m_APDDigitizer
EBDigitizer * m_APDDigitizer
Definition: EcalDigiProducer.h:136
EcalDigiProducer::m_APDElectronicsSim
EcalElectronicsSim * m_APDElectronicsSim
Definition: EcalDigiProducer.h:143
EcalDigiProducer::m_EEResponse
EEHitResponse * m_EEResponse
Definition: EcalDigiProducer.h:120

Member Function Documentation

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

Implements DigiAccumulatorMixMod.

Definition at line 555 of file EcalDigiProducer.cc.

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

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

555  {
556  // Step A: Get Inputs
557  edm::InputTag ebTag(m_hitsProducerTag, "EcalHitsEB");
558  edm::Handle<std::vector<PCaloHit> > ebHandle;
559  e.getByLabel(ebTag, ebHandle);
560 
561  edm::InputTag eeTag(m_hitsProducerTag, "EcalHitsEE");
562  edm::Handle<std::vector<PCaloHit> > eeHandle;
563  e.getByLabel(eeTag, eeHandle);
564 
565  edm::InputTag esTag(m_hitsProducerTag, "EcalHitsES");
566  edm::Handle<std::vector<PCaloHit> > esHandle;
567  e.getByLabel(esTag, esHandle);
568 
569  accumulateCaloHits(ebHandle, eeHandle, esHandle, 0, randomEngine(e.streamID()));
570 }
EcalDigiProducer::m_hitsProducerTag
const std::string m_hitsProducerTag
Definition: EcalDigiProducer.h:101
EcalDigiProducer::randomEngine
CLHEP::HepRandomEngine * randomEngine(edm::StreamID const &streamID)
Definition: EcalDigiProducer.cc:821
EcalDigiProducer::accumulateCaloHits
void accumulateCaloHits(HitsHandle const &ebHandle, HitsHandle const &eeHandle, HitsHandle const &esHandle, int bunchCrossing, CLHEP::HepRandomEngine *)
Definition: EcalDigiProducer.cc:532
void EcalDigiProducer::accumulate ( PileUpEventPrincipal const &  e,
edm::EventSetup const &  c,
edm::StreamID const &  streamID 
)
overridevirtual

Implements DigiAccumulatorMixMod.

Definition at line 573 of file EcalDigiProducer.cc.

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

573  {
574  // Step A: Get Inputs
575  edm::InputTag ebTag(m_hitsProducerTag, "EcalHitsEB");
576  edm::Handle<std::vector<PCaloHit> > ebHandle;
577  e.getByLabel(ebTag, ebHandle);
578 
579  edm::InputTag eeTag(m_hitsProducerTag, "EcalHitsEE");
580  edm::Handle<std::vector<PCaloHit> > eeHandle;
581  e.getByLabel(eeTag, eeHandle);
582 
583  edm::InputTag esTag(m_hitsProducerTag, "EcalHitsES");
584  edm::Handle<std::vector<PCaloHit> > esHandle;
585  e.getByLabel(esTag, esHandle);
586 
587  accumulateCaloHits(ebHandle, eeHandle, esHandle, e.bunchCrossing(), randomEngine(streamID));
588 }
EcalDigiProducer::m_hitsProducerTag
const std::string m_hitsProducerTag
Definition: EcalDigiProducer.h:101
EcalDigiProducer::randomEngine
CLHEP::HepRandomEngine * randomEngine(edm::StreamID const &streamID)
Definition: EcalDigiProducer.cc:821
EcalDigiProducer::accumulateCaloHits
void accumulateCaloHits(HitsHandle const &ebHandle, HitsHandle const &eeHandle, HitsHandle const &esHandle, int bunchCrossing, CLHEP::HepRandomEngine *)
Definition: EcalDigiProducer.cc:532
void EcalDigiProducer::accumulateCaloHits ( HitsHandle const &  ebHandle,
HitsHandle const &  eeHandle,
HitsHandle const &  esHandle,
int  bunchCrossing,
CLHEP::HepRandomEngine *  engine 
)
private

Definition at line 532 of file EcalDigiProducer.cc.

References EcalTDigitizer< Traits >::add(), CaloTDigitizer< Traits >::add(), edm::HandleBase::isValid(), m_APDDigitizer, m_apdSeparateDigi, m_BarrelDigitizer, m_doFastES, m_EndcapDigitizer, m_ESDigitizer, m_ESOldDigitizer, and edm::Handle< T >::product().

Referenced by accumulate().

532  {
533  if(ebHandle.isValid()) {
534  m_BarrelDigitizer->add(*ebHandle.product(), bunchCrossing, engine);
535 
536  if(m_apdSeparateDigi) {
537  m_APDDigitizer->add(*ebHandle.product(), bunchCrossing, engine);
538  }
539  }
540 
541  if(eeHandle.isValid()) {
542  m_EndcapDigitizer->add(*eeHandle.product(), bunchCrossing, engine);
543  }
544 
545  if(esHandle.isValid()) {
546  if(m_doFastES) {
547  m_ESDigitizer->add(*esHandle.product(), bunchCrossing, engine);
548  } else {
549  m_ESOldDigitizer->add(*esHandle.product(), bunchCrossing, engine);
550  }
551  }
552 }
EcalDigiProducer::m_EndcapDigitizer
EEDigitizer * m_EndcapDigitizer
Definition: EcalDigiProducer.h:138
EcalTDigitizer::add
void add(const std::vector< PCaloHit > &hits, int bunchCrossing, CLHEP::HepRandomEngine *)
CaloTDigitizer::add
void add(const std::vector< PCaloHit > &hits, int bunchCrossing, CLHEP::HepRandomEngine *engine)
Definition: CaloTDigitizer.h:60
EcalDigiProducer::m_ESOldDigitizer
ESOldDigitizer * m_ESOldDigitizer
Definition: EcalDigiProducer.h:132
EcalDigiProducer::m_ESDigitizer
ESDigitizer * m_ESDigitizer
Definition: EcalDigiProducer.h:134
EcalDigiProducer::m_BarrelDigitizer
EBDigitizer * m_BarrelDigitizer
Definition: EcalDigiProducer.h:137
EcalDigiProducer::m_APDDigitizer
EBDigitizer * m_APDDigitizer
Definition: EcalDigiProducer.h:136
EcalDigiProducer::m_apdSeparateDigi
const bool m_apdSeparateDigi
Definition: EcalDigiProducer.h:105
void EcalDigiProducer::beginLuminosityBlock ( edm::LuminosityBlock const &  lumi,
edm::EventSetup const &  setup 
)
overridevirtual

Reimplemented from DigiAccumulatorMixMod.

Definition at line 634 of file EcalDigiProducer.cc.

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

635 {
636  edm::Service<edm::RandomNumberGenerator> rng;
637  if ( ! rng.isAvailable() ) {
638  throw cms::Exception("Configuration") <<
639  "RandomNumberGenerator service is not available.\n"
640  "You must add the service in the configuration file\n"
641  "or remove the module that requires it.";
642  }
643  CLHEP::HepRandomEngine* engine = &rng->getEngine(lumi.index());
644 
645  if( 0 != m_APDResponse ) m_APDResponse->initialize(engine);
646  m_EBResponse->initialize(engine);
647 }
fjr2json.lumi
tuple lumi
Definition: fjr2json.py:35
EcalDigiProducer::m_APDResponse
EBHitResponse * m_APDResponse
Definition: EcalDigiProducer.h:117
EcalDigiProducer::m_EBResponse
EBHitResponse * m_EBResponse
Definition: EcalDigiProducer.h:119
edm::Service::isAvailable
bool isAvailable() const
Definition: Service.h:46
edm::RandomNumberGenerator::getEngine
virtual CLHEP::HepRandomEngine & getEngine(StreamID const &) const =0
Use this engine in event methods.
EBHitResponse::initialize
void initialize(CLHEP::HepRandomEngine *)
Definition: EBHitResponse.cc:62
virtual void EcalDigiProducer::cacheEBDigis ( const EBDigiCollection *  ebDigiPtr) const
inlineprivatevirtual

Reimplemented in EcalTBDigiProducer.

Definition at line 79 of file EcalDigiProducer.h.

Referenced by finalizeEvent().

79 { }
virtual void EcalDigiProducer::cacheEEDigis ( const EEDigiCollection *  eeDigiPtr) const
inlineprivatevirtual

Reimplemented in EcalTBDigiProducer.

Definition at line 80 of file EcalDigiProducer.h.

Referenced by finalizeEvent().

80 { }
void EcalDigiProducer::checkCalibrations ( const edm::Event &  event,
const edm::EventSetup &  eventSetup 
)
private

Definition at line 650 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_doFastES, m_EBResponse, m_EBs25notCont, m_EEResponse, m_EEs25notCont, m_ESDigitizer, m_ESElectronicsSim, m_ESElectronicsSimFast, m_ESOldDigitizer, m_ESShape, m_useLCcorrection, EcalCoder::MAXADC, EcalCoder::NGAINS, edm::ESHandle< class >::product(), EcalHitResponse::setEventTime(), EcalCoder::setFullScaleEnergy(), ESShape::setGain(), ESElectronicsSim::setGain(), ESDigitizer::setGain(), EcalCoder::setGainRatios(), EBHitResponse::setIntercal(), EcalCoder::setIntercalibConstants(), EcalHitResponse::setLaserConstants(), ESElectronicsSim::setMIPs(), ESElectronicsSimFast::setMIPs(), ESElectronicsSim::setMIPToGeV(), ESElectronicsSimFast::setMIPToGeV(), ESElectronicsSim::setPedestals(), ESElectronicsSimFast::setPedestals(), and EcalCoder::setPedestals().

Referenced by initializeEvent().

651 {
652  // Pedestals from event setup
653 
654  edm::ESHandle<EcalPedestals> dbPed ;
655  eventSetup.get<EcalPedestalsRcd>().get( dbPed ) ;
656  const EcalPedestals* pedestals ( dbPed.product() ) ;
657 
658  m_Coder->setPedestals( pedestals ) ;
659  if( 0 != m_APDCoder ) m_APDCoder->setPedestals( pedestals ) ;
660 
661  // Ecal Intercalibration Constants
662  edm::ESHandle<EcalIntercalibConstantsMC> pIcal ;
663  eventSetup.get<EcalIntercalibConstantsMCRcd>().get( pIcal ) ;
664  const EcalIntercalibConstantsMC* ical ( pIcal.product() ) ;
665 
666  m_Coder->setIntercalibConstants( ical ) ;
667  if( 0 != m_APDCoder) m_APDCoder->setIntercalibConstants( ical ) ;
668 
669  m_EBResponse->setIntercal( ical ) ;
670  if( 0 != m_APDResponse ) m_APDResponse->setIntercal( ical ) ;
671 
672  // Ecal LaserCorrection Constants
673  edm::ESHandle<EcalLaserDbService> laser;
674  eventSetup.get<EcalLaserDbRecord>().get(laser);
675  const edm::TimeValue_t eventTimeValue = event.time().value();
676 
677  m_EBResponse->setEventTime(eventTimeValue);
678  m_EBResponse->setLaserConstants(laser.product(), m_useLCcorrection);
679 
680  m_EEResponse->setEventTime(eventTimeValue);
681  m_EEResponse->setLaserConstants(laser.product(), m_useLCcorrection);
682 
683  // ADC -> GeV Scale
684  edm::ESHandle<EcalADCToGeVConstant> pAgc;
685  eventSetup.get<EcalADCToGeVConstantRcd>().get(pAgc);
686  const EcalADCToGeVConstant* agc = pAgc.product();
687 
688  // Gain Ratios
689  edm::ESHandle<EcalGainRatios> pRatio;
690  eventSetup.get<EcalGainRatiosRcd>().get(pRatio);
691  const EcalGainRatios* gr = pRatio.product();
692 
693  m_Coder->setGainRatios( gr );
694  if( 0 != m_APDCoder) m_APDCoder->setGainRatios( gr );
695 
696  EcalMGPAGainRatio * defaultRatios = new EcalMGPAGainRatio();
697 
698  double theGains[m_Coder->NGAINS+1];
699  theGains[0] = 0.;
700  theGains[3] = 1.;
701  theGains[2] = defaultRatios->gain6Over1() ;
702  theGains[1] = theGains[2]*(defaultRatios->gain12Over6()) ;
703 
704  LogDebug("EcalDigi") << " Gains: " << "\n" << " g1 = " << theGains[1]
705  << "\n" << " g2 = " << theGains[2]
706  << "\n" << " g3 = " << theGains[3] ;
707 
708  delete defaultRatios;
709 
710  const double EBscale (
711  ( agc->getEBValue())*theGains[1]*(m_Coder->MAXADC)*m_EBs25notCont ) ;
712 
713  LogDebug("EcalDigi") << " GeV/ADC = " << agc->getEBValue()
714  << "\n" << " notCont = " << m_EBs25notCont
715  << "\n" << " saturation for EB = " << EBscale
716  << ", " << m_EBs25notCont ;
717 
718  const double EEscale (
719  (agc->getEEValue())*theGains[1]*(m_Coder->MAXADC)*m_EEs25notCont ) ;
720 
721  LogDebug("EcalDigi") << " GeV/ADC = " << agc->getEEValue()
722  << "\n" << " notCont = " << m_EEs25notCont
723  << "\n" << " saturation for EB = " << EEscale
724  << ", " << m_EEs25notCont ;
725 
726  m_Coder->setFullScaleEnergy( EBscale ,
727  EEscale ) ;
728  if( 0 != m_APDCoder ) m_APDCoder->setFullScaleEnergy( EBscale ,
729  EEscale ) ;
730 
731  if( 0 != m_ESOldDigitizer ||
732  0 != m_ESDigitizer )
733  {
734  // ES condition objects
735  edm::ESHandle<ESGain> hesgain ;
736  edm::ESHandle<ESMIPToGeVConstant> hesMIPToGeV ;
737  edm::ESHandle<ESPedestals> hesPedestals ;
738  edm::ESHandle<ESIntercalibConstants> hesMIPs ;
739 
740  eventSetup.get<ESGainRcd>(). get( hesgain ) ;
741  eventSetup.get<ESMIPToGeVConstantRcd>(). get( hesMIPToGeV ) ;
742  eventSetup.get<ESPedestalsRcd>(). get( hesPedestals ) ;
743  eventSetup.get<ESIntercalibConstantsRcd>().get( hesMIPs ) ;
744 
745  const ESGain* esgain ( hesgain.product() ) ;
746  const ESPedestals* espeds ( hesPedestals.product() ) ;
747  const ESIntercalibConstants* esmips ( hesMIPs.product() ) ;
748  const ESMIPToGeVConstant* esMipToGeV ( hesMIPToGeV.product() ) ;
749  const int ESGain ( 1.1 > esgain->getESGain() ? 1 : 2 ) ;
750  const double ESMIPToGeV ( ( 1 == ESGain ) ?
751  esMipToGeV->getESValueLow() :
752  esMipToGeV->getESValueHigh() ) ;
753 
754  m_ESShape.setGain( ESGain );
755 
756  if( !m_doFastES )
757  {
758  m_ESElectronicsSim->setGain( ESGain ) ;
759  m_ESElectronicsSim->setPedestals( espeds ) ;
760  m_ESElectronicsSim->setMIPs( esmips ) ;
761  m_ESElectronicsSim->setMIPToGeV( ESMIPToGeV ) ;
762  }
763  else
764  {
765  m_ESDigitizer->setGain( ESGain ) ;
766  m_ESElectronicsSimFast->setPedestals( espeds ) ;
767  m_ESElectronicsSimFast->setMIPs( esmips ) ;
768  m_ESElectronicsSimFast->setMIPToGeV( ESMIPToGeV ) ;
769  }
770  }
771 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
EcalCondDBWriter_cfi.laser
string laser
Definition: EcalCondDBWriter_cfi.py:45
ESElectronicsSimFast::setMIPToGeV
void setMIPToGeV(double MIPToGeV)
Definition: ESElectronicsSimFast.cc:38
ESElectronicsSim::setGain
void setGain(const int gain)
Definition: ESElectronicsSim.h:26
ESGain
Definition: ESGain.h:7
ESShape::setGain
void setGain(const int gain)
Definition: ESShape.h:26
EcalHitResponse::setEventTime
void setEventTime(const edm::TimeValue_t &iTime)
Definition: EcalHitResponse.cc:109
ESElectronicsSim::setMIPs
void setMIPs(const ESIntercalibConstants *mips)
Definition: ESElectronicsSim.h:28
ESElectronicsSim::setMIPToGeV
void setMIPToGeV(const double MIPToGeV)
Definition: ESElectronicsSim.h:29
EcalCoder::setIntercalibConstants
void setIntercalibConstants(const EcalIntercalibConstantsMC *ical)
Definition: EcalCoder.cc:63
EcalDigiProducer::m_ESOldDigitizer
ESOldDigitizer * m_ESOldDigitizer
Definition: EcalDigiProducer.h:132
EcalDigiProducer::m_ESElectronicsSimFast
ESElectronicsSimFast * m_ESElectronicsSimFast
Definition: EcalDigiProducer.h:133
EcalDigiProducer::m_ESDigitizer
ESDigitizer * m_ESDigitizer
Definition: EcalDigiProducer.h:134
ESElectronicsSimFast::setPedestals
void setPedestals(const ESPedestals *peds)
Definition: ESElectronicsSimFast.cc:26
ESElectronicsSim::setPedestals
void setPedestals(const ESPedestals *peds)
Definition: ESElectronicsSim.h:27
EcalHitResponse::setLaserConstants
void setLaserConstants(const EcalLaserDbService *laser, bool &useLCcorrection)
Definition: EcalHitResponse.cc:115
EBHitResponse::setIntercal
void setIntercal(const EcalIntercalibConstantsMC *ical)
Definition: EBHitResponse.cc:158
EcalDigiProducer::m_APDResponse
EBHitResponse * m_APDResponse
Definition: EcalDigiProducer.h:117
EcalDigiProducer::m_EBs25notCont
const double m_EBs25notCont
Definition: EcalDigiProducer.h:107
EcalCoder::setPedestals
void setPedestals(const EcalPedestals *pedestals)
can be fetched every event from the EventSetup
Definition: EcalCoder.cc:51
EcalDigiProducer::m_EBResponse
EBHitResponse * m_EBResponse
Definition: EcalDigiProducer.h:119
ESDigitizer::setGain
void setGain(const int gain)
Definition: ESDigitizer.cc:44
EcalDigiProducer::m_ESElectronicsSim
ESElectronicsSim * m_ESElectronicsSim
Definition: EcalDigiProducer.h:131
EcalMGPAGainRatio::gain6Over1
float gain6Over1() const
Definition: EcalMGPAGainRatio.h:21
edm::TimeValue_t
unsigned long long TimeValue_t
Definition: Timestamp.h:28
EcalDigiProducer::m_APDCoder
EcalCoder * m_APDCoder
Definition: EcalDigiProducer.h:144
EcalCoder::setGainRatios
void setGainRatios(const EcalGainRatios *gainRatios)
Definition: EcalCoder.cc:57
EcalCoder::setFullScaleEnergy
void setFullScaleEnergy(double EBscale, double EEscale)
Definition: EcalCoder.cc:42
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
EcalMGPAGainRatio::gain12Over6
float gain12Over6() const
Definition: EcalMGPAGainRatio.h:20
ESElectronicsSimFast::setMIPs
void setMIPs(const ESIntercalibConstants *mips)
Definition: ESElectronicsSimFast.cc:32
EcalDigiProducer::m_EEs25notCont
const double m_EEs25notCont
Definition: EcalDigiProducer.h:108
EcalDigiProducer::m_EEResponse
EEHitResponse * m_EEResponse
Definition: EcalDigiProducer.h:120
void EcalDigiProducer::checkGeometry ( const edm::EventSetup &  eventSetup)
private

Definition at line 774 of file EcalDigiProducer.cc.

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

Referenced by initializeEvent().

775 {
776  // TODO find a way to avoid doing this every event
777  edm::ESHandle<CaloGeometry> hGeometry ;
778  eventSetup.get<CaloGeometryRecord>().get( hGeometry ) ;
779 
780  const CaloGeometry* pGeometry = &*hGeometry;
781 
782  if( pGeometry != m_Geometry )
783  {
784  m_Geometry = pGeometry;
785  updateGeometry();
786  }
787 }
EcalDigiProducer::m_Geometry
const CaloGeometry * m_Geometry
Definition: EcalDigiProducer.h:146
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
void EcalDigiProducer::finalizeEvent ( edm::Event &  e,
edm::EventSetup const &  c 
)
virtual

Implements DigiAccumulatorMixMod.

Reimplemented in EcalTBDigiProducer.

Definition at line 591 of file EcalDigiProducer.cc.

References cacheEBDigis(), cacheEEDigis(), m_APDDigitizer, m_apdSeparateDigi, m_BarrelDigitizer, m_doFastES, m_EBdigiCollection, m_EEdigiCollection, m_EndcapDigitizer, m_ESdigiCollection, m_ESDigitizer, m_ESOldDigitizer, randomEngine(), ESDigitizer::run(), EcalTDigitizer< Traits >::run(), CaloTDigitizer< Traits >::run(), and edm::Event::streamID().

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

591  {
592  // Step B: Create empty output
593  std::auto_ptr<EBDigiCollection> apdResult ( !m_apdSeparateDigi ? 0 :
594  new EBDigiCollection() ) ;
595  std::auto_ptr<EBDigiCollection> barrelResult ( new EBDigiCollection() ) ;
596  std::auto_ptr<EEDigiCollection> endcapResult ( new EEDigiCollection() ) ;
597  std::auto_ptr<ESDigiCollection> preshowerResult( new ESDigiCollection() ) ;
598 
599  // run the algorithm
600 
601  m_BarrelDigitizer->run( *barrelResult, randomEngine(event.streamID()) ) ;
602  cacheEBDigis( &*barrelResult ) ;
603 
604  edm::LogInfo("DigiInfo") << "EB Digis: " << barrelResult->size() ;
605 
606  if( m_apdSeparateDigi ) {
607  m_APDDigitizer->run( *apdResult, randomEngine(event.streamID()) ) ;
608  edm::LogInfo("DigiInfo") << "APD Digis: " << apdResult->size() ;
609  }
610 
611  m_EndcapDigitizer->run( *endcapResult, randomEngine(event.streamID()) ) ;
612  edm::LogInfo("EcalDigi") << "EE Digis: " << endcapResult->size() ;
613  cacheEEDigis( &*endcapResult ) ;
614 
615  if(m_doFastES) {
616  m_ESDigitizer->run( *preshowerResult, randomEngine(event.streamID()) ) ;
617  } else {
618  m_ESOldDigitizer->run( *preshowerResult, randomEngine(event.streamID()) ) ;
619  }
620  edm::LogInfo("EcalDigi") << "ES Digis: " << preshowerResult->size();
621 
622 
623  // Step D: Put outputs into event
624  if( m_apdSeparateDigi ) {
625  //event.put( apdResult, m_apdDigiTag ) ;
626  }
627 
628  event.put( barrelResult, m_EBdigiCollection ) ;
629  event.put( endcapResult, m_EEdigiCollection ) ;
630  event.put( preshowerResult, m_ESdigiCollection ) ;
631 }
EcalDigiProducer::m_EndcapDigitizer
EEDigitizer * m_EndcapDigitizer
Definition: EcalDigiProducer.h:138
EcalDigiProducer::m_ESdigiCollection
const std::string m_ESdigiCollection
Definition: EcalDigiProducer.h:100
EcalDigiProducer::m_ESOldDigitizer
ESOldDigitizer * m_ESOldDigitizer
Definition: EcalDigiProducer.h:132
EcalDigiProducer::m_EEdigiCollection
const std::string m_EEdigiCollection
Definition: EcalDigiProducer.h:99
EcalDigiProducer::m_ESDigitizer
ESDigitizer * m_ESDigitizer
Definition: EcalDigiProducer.h:134
EcalDigiProducer::m_BarrelDigitizer
EBDigitizer * m_BarrelDigitizer
Definition: EcalDigiProducer.h:137
event
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
CaloTDigitizer::run
void run(MixCollection< PCaloHit > &, DigiCollection &)
turns hits into digis
Definition: CaloTDigitizer.h:73
EcalDigiProducer::cacheEBDigis
virtual void cacheEBDigis(const EBDigiCollection *ebDigiPtr) const
Definition: EcalDigiProducer.h:79
EcalDigiProducer::m_APDDigitizer
EBDigitizer * m_APDDigitizer
Definition: EcalDigiProducer.h:136
EcalDigiProducer::m_apdSeparateDigi
const bool m_apdSeparateDigi
Definition: EcalDigiProducer.h:105
EcalDigiProducer::randomEngine
CLHEP::HepRandomEngine * randomEngine(edm::StreamID const &streamID)
Definition: EcalDigiProducer.cc:821
EcalTDigitizer::run
virtual void run(DigiCollection &output, CLHEP::HepRandomEngine *)
ESDigitizer::run
virtual void run(ESDigiCollection &output, CLHEP::HepRandomEngine *) override
turns hits into digis
Definition: ESDigitizer.cc:167
EcalDigiProducer::cacheEEDigis
virtual void cacheEEDigis(const EEDigiCollection *eeDigiPtr) const
Definition: EcalDigiProducer.h:80
EcalDigiProducer::m_EBdigiCollection
const std::string m_EBdigiCollection
Definition: EcalDigiProducer.h:98
void EcalDigiProducer::initializeEvent ( edm::Event const &  e,
edm::EventSetup const &  c 
)
virtual

Implements DigiAccumulatorMixMod.

Reimplemented in EcalTBDigiProducer.

Definition at line 516 of file EcalDigiProducer.cc.

References checkCalibrations(), checkGeometry(), EcalTDigitizer< Traits >::initializeHits(), CaloTDigitizer< Traits >::initializeHits(), m_APDDigitizer, m_apdSeparateDigi, m_BarrelDigitizer, m_doFastES, m_EndcapDigitizer, m_ESDigitizer, and m_ESOldDigitizer.

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

516  {
517  checkGeometry( eventSetup );
518  checkCalibrations( event, eventSetup );
519  m_BarrelDigitizer->initializeHits();
520  if(m_apdSeparateDigi) {
521  m_APDDigitizer->initializeHits();
522  }
523  m_EndcapDigitizer->initializeHits();
524  if(m_doFastES) {
525  m_ESDigitizer->initializeHits();
526  } else {
527  m_ESOldDigitizer->initializeHits();
528  }
529 }
EcalDigiProducer::m_EndcapDigitizer
EEDigitizer * m_EndcapDigitizer
Definition: EcalDigiProducer.h:138
CaloTDigitizer::initializeHits
void initializeHits()
Definition: CaloTDigitizer.h:68
EcalDigiProducer::m_ESOldDigitizer
ESOldDigitizer * m_ESOldDigitizer
Definition: EcalDigiProducer.h:132
EcalTDigitizer::initializeHits
virtual void initializeHits()
EcalDigiProducer::m_ESDigitizer
ESDigitizer * m_ESDigitizer
Definition: EcalDigiProducer.h:134
EcalDigiProducer::m_BarrelDigitizer
EBDigitizer * m_BarrelDigitizer
Definition: EcalDigiProducer.h:137
event
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
EcalDigiProducer::checkCalibrations
void checkCalibrations(const edm::Event &event, const edm::EventSetup &eventSetup)
Definition: EcalDigiProducer.cc:650
EcalDigiProducer::checkGeometry
void checkGeometry(const edm::EventSetup &eventSetup)
Definition: EcalDigiProducer.cc:774
EcalDigiProducer::m_APDDigitizer
EBDigitizer * m_APDDigitizer
Definition: EcalDigiProducer.h:136
EcalDigiProducer::m_apdSeparateDigi
const bool m_apdSeparateDigi
Definition: EcalDigiProducer.h:105
CLHEP::HepRandomEngine * EcalDigiProducer::randomEngine ( edm::StreamID const &  streamID)
private

Definition at line 821 of file EcalDigiProducer.cc.

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

Referenced by accumulate(), and finalizeEvent().

821  {
822  unsigned int index = streamID.value();
823  if(index >= randomEngines_.size()) {
824  randomEngines_.resize(index + 1, nullptr);
825  }
826  CLHEP::HepRandomEngine* ptr = randomEngines_[index];
827  if(!ptr) {
828  edm::Service<edm::RandomNumberGenerator> rng;
829  ptr = &rng->getEngine(streamID);
830  randomEngines_[index] = ptr;
831  }
832  return ptr;
833 }
edm::RandomNumberGenerator::getEngine
virtual CLHEP::HepRandomEngine & getEngine(StreamID const &) const =0
Use this engine in event methods.
EcalDigiProducer::randomEngines_
std::vector< CLHEP::HepRandomEngine * > randomEngines_
Definition: EcalDigiProducer.h:151
void EcalDigiProducer::setEBNoiseSignalGenerator ( EcalBaseSignalGenerator *  noiseGenerator)

Definition at line 835 of file EcalDigiProducer.cc.

References m_BarrelDigitizer, and EcalTDigitizer< Traits >::setNoiseSignalGenerator().

835  {
836  //noiseGenerator->setParameterMap(theParameterMap);
837  if(m_BarrelDigitizer) m_BarrelDigitizer->setNoiseSignalGenerator(noiseGenerator);
838 }
EcalTDigitizer::setNoiseSignalGenerator
void setNoiseSignalGenerator(EcalBaseSignalGenerator *noiseSignalGenerator)
EcalDigiProducer::m_BarrelDigitizer
EBDigitizer * m_BarrelDigitizer
Definition: EcalDigiProducer.h:137
void EcalDigiProducer::setEENoiseSignalGenerator ( EcalBaseSignalGenerator *  noiseGenerator)

Definition at line 840 of file EcalDigiProducer.cc.

References m_EndcapDigitizer, and EcalTDigitizer< Traits >::setNoiseSignalGenerator().

840  {
841  //noiseGenerator->setParameterMap(theParameterMap);
842  if(m_EndcapDigitizer) m_EndcapDigitizer->setNoiseSignalGenerator(noiseGenerator);
843 }
EcalDigiProducer::m_EndcapDigitizer
EEDigitizer * m_EndcapDigitizer
Definition: EcalDigiProducer.h:138
EcalTDigitizer::setNoiseSignalGenerator
void setNoiseSignalGenerator(EcalBaseSignalGenerator *noiseSignalGenerator)
void EcalDigiProducer::setESNoiseSignalGenerator ( EcalBaseSignalGenerator *  noiseGenerator)

Definition at line 845 of file EcalDigiProducer.cc.

References m_ESDigitizer, and EcalTDigitizer< Traits >::setNoiseSignalGenerator().

845  {
846  //noiseGenerator->setParameterMap(theParameterMap);
847  if(m_ESDigitizer) m_ESDigitizer->setNoiseSignalGenerator(noiseGenerator);
848 }
EcalTDigitizer::setNoiseSignalGenerator
void setNoiseSignalGenerator(EcalBaseSignalGenerator *noiseSignalGenerator)
EcalDigiProducer::m_ESDigitizer
ESDigitizer * m_ESDigitizer
Definition: EcalDigiProducer.h:134
void EcalDigiProducer::updateGeometry ( )
private

Definition at line 790 of file EcalDigiProducer.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, EcalPreshower, CaloGeometry::getSubdetectorGeometry(), CaloSubdetectorGeometry::getValidDetIds(), m_APDResponse, m_doFastES, m_EBResponse, m_EEResponse, m_ESDigitizer, m_ESOldDigitizer, m_ESOldResponse, m_ESResponse, m_Geometry, ESDigitizer::setDetIds(), CaloTDigitizer< Traits >::setDetIds(), EcalHitResponse::setGeometry(), and CaloHitResponse::setGeometry().

Referenced by checkGeometry().

791 {
792  if( 0 != m_APDResponse ) m_APDResponse->setGeometry(
793  m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalBarrel ) ) ;
794  m_EBResponse->setGeometry(
795  m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalBarrel ) ) ;
796  m_EEResponse->setGeometry(
797  m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalEndcap ) ) ;
798  m_ESResponse->setGeometry(
799  m_Geometry->getSubdetectorGeometry( DetId::Ecal, EcalPreshower ) ) ;
800  m_ESOldResponse->setGeometry( m_Geometry ) ;
801 
802  const std::vector<DetId>* theESDets (
803  0 != m_Geometry->getSubdetectorGeometry(DetId::Ecal, EcalPreshower) ?
804  &m_Geometry->getSubdetectorGeometry(DetId::Ecal, EcalPreshower)->getValidDetIds() : 0 ) ;
805 
806  if( !m_doFastES )
807  {
808  if( 0 != m_ESOldDigitizer &&
809  0 != theESDets )
810  m_ESOldDigitizer->setDetIds( *theESDets ) ;
811  }
812  else
813  {
814  if( 0 != m_ESDigitizer &&
815  0 != theESDets )
816  m_ESDigitizer->setDetIds( *theESDets ) ;
817  }
818 }
CaloGeometry::getSubdetectorGeometry
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:43
CaloHitResponse::setGeometry
void setGeometry(const CaloGeometry *geometry)
geometry needed for time-of-flight
Definition: CaloHitResponse.h:55
EcalDigiProducer::m_ESOldResponse
CaloHitResponse * m_ESOldResponse
Definition: EcalDigiProducer.h:123
ESDigitizer::setDetIds
void setDetIds(const std::vector< DetId > &detIds)
tell the digitizer which cells exist; cannot change during a run
Definition: ESDigitizer.cc:36
EcalDigiProducer::m_ESOldDigitizer
ESOldDigitizer * m_ESOldDigitizer
Definition: EcalDigiProducer.h:132
EcalDigiProducer::m_ESDigitizer
ESDigitizer * m_ESDigitizer
Definition: EcalDigiProducer.h:134
CaloSubdetectorGeometry::getValidDetIds
virtual const std::vector< DetId > & getValidDetIds(DetId::Detector det=DetId::Detector(0), int subdet=0) const
Get a list of valid detector ids (for the given subdetector)
Definition: CaloSubdetectorGeometry.cc:37
EcalDigiProducer::m_APDResponse
EBHitResponse * m_APDResponse
Definition: EcalDigiProducer.h:117
EcalDigiProducer::m_EBResponse
EBHitResponse * m_EBResponse
Definition: EcalDigiProducer.h:119
EcalDigiProducer::m_Geometry
const CaloGeometry * m_Geometry
Definition: EcalDigiProducer.h:146
EcalDigiProducer::m_ESResponse
ESHitResponse * m_ESResponse
Definition: EcalDigiProducer.h:122
EcalHitResponse::setGeometry
void setGeometry(const CaloSubdetectorGeometry *geometry)
Definition: EcalHitResponse.cc:73
CaloTDigitizer::setDetIds
void setDetIds(const std::vector< DetId > &detIds)
Definition: CaloTDigitizer.h:48
EcalDigiProducer::m_EEResponse
EEHitResponse * m_EEResponse
Definition: EcalDigiProducer.h:120

Member Data Documentation

const bool EcalDigiProducer::m_addESNoise
private

Definition at line 125 of file EcalDigiProducer.h.

EcalCoder* EcalDigiProducer::m_APDCoder
private

Definition at line 144 of file EcalDigiProducer.h.

Referenced by checkCalibrations(), EcalDigiProducer(), and ~EcalDigiProducer().

const std::string EcalDigiProducer::m_apdDigiTag
private

Definition at line 114 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer().

EBDigitizer* EcalDigiProducer::m_APDDigitizer
private

Definition at line 136 of file EcalDigiProducer.h.

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

EcalElectronicsSim* EcalDigiProducer::m_APDElectronicsSim
private

Definition at line 143 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer(), and ~EcalDigiProducer().

const APDSimParameters* EcalDigiProducer::m_apdParameters
private

Definition at line 115 of file EcalDigiProducer.h.

Referenced by ~EcalDigiProducer().

EBHitResponse* EcalDigiProducer::m_APDResponse
private

Definition at line 117 of file EcalDigiProducer.h.

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

const bool EcalDigiProducer::m_apdSeparateDigi
private

Definition at line 105 of file EcalDigiProducer.h.

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

const APDShape EcalDigiProducer::m_APDShape
private

Definition at line 93 of file EcalDigiProducer.h.

EBDigitizer* EcalDigiProducer::m_BarrelDigitizer
private

Definition at line 137 of file EcalDigiProducer.h.

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

EcalCoder* EcalDigiProducer::m_Coder
private

Definition at line 141 of file EcalDigiProducer.h.

Referenced by checkCalibrations(), EcalDigiProducer(), and ~EcalDigiProducer().

const bool EcalDigiProducer::m_doFastES
private

Definition at line 129 of file EcalDigiProducer.h.

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

CorrelatedNoisifier<EcalCorrMatrix>* EcalDigiProducer::m_EBCorrNoise[3]
private

Definition at line 148 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer(), and ~EcalDigiProducer().

const std::string EcalDigiProducer::m_EBdigiCollection
private

Definition at line 98 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer(), and finalizeEvent().

EBHitResponse* EcalDigiProducer::m_EBResponse
protected

Definition at line 119 of file EcalDigiProducer.h.

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

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.

CorrelatedNoisifier<EcalCorrMatrix>* EcalDigiProducer::m_EECorrNoise[3]
private

Definition at line 149 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer(), and ~EcalDigiProducer().

const std::string EcalDigiProducer::m_EEdigiCollection
private

Definition at line 99 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer(), and finalizeEvent().

EEHitResponse* EcalDigiProducer::m_EEResponse
protected

Definition at line 120 of file EcalDigiProducer.h.

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

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.

EcalElectronicsSim* EcalDigiProducer::m_ElectronicsSim
private

Definition at line 140 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer(), and ~EcalDigiProducer().

EEDigitizer* EcalDigiProducer::m_EndcapDigitizer
private

Definition at line 138 of file EcalDigiProducer.h.

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

const std::string EcalDigiProducer::m_ESdigiCollection
private

Definition at line 100 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer(), and finalizeEvent().

ESDigitizer* EcalDigiProducer::m_ESDigitizer
private

Definition at line 134 of file EcalDigiProducer.h.

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

ESElectronicsSim* EcalDigiProducer::m_ESElectronicsSim
private

Definition at line 131 of file EcalDigiProducer.h.

Referenced by checkCalibrations(), and ~EcalDigiProducer().

ESElectronicsSimFast* EcalDigiProducer::m_ESElectronicsSimFast
private

Definition at line 133 of file EcalDigiProducer.h.

Referenced by checkCalibrations(), and ~EcalDigiProducer().

ESOldDigitizer* EcalDigiProducer::m_ESOldDigitizer
private

Definition at line 132 of file EcalDigiProducer.h.

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

CaloHitResponse* EcalDigiProducer::m_ESOldResponse
private

Definition at line 123 of file EcalDigiProducer.h.

Referenced by updateGeometry(), and ~EcalDigiProducer().

ESHitResponse* EcalDigiProducer::m_ESResponse
private

Definition at line 122 of file EcalDigiProducer.h.

Referenced by updateGeometry(), and ~EcalDigiProducer().

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 146 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(), and EcalDigiProducer().

const EcalSimParameterMap* EcalDigiProducer::m_ParameterMap
protected

Definition at line 112 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer(), EcalTBDigiProducer::setPhaseShift(), and ~EcalDigiProducer().

const bool EcalDigiProducer::m_PreMix1
private

Definition at line 126 of file EcalDigiProducer.h.

Referenced by EcalDigiProducer().

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.

Referenced by EcalDigiProducer().

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 151 of file EcalDigiProducer.h.

Referenced by randomEngine().

Generated for CMSSW Reference Manual by doxygen 1.8.5