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 Member Functions | Private Types | Private Attributes
EcalEndcapDigisValidation Class Reference

#include <EcalEndcapDigisValidation.h>

Inheritance diagram for EcalEndcapDigisValidation:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 EcalEndcapDigisValidation (const edm::ParameterSet &ps)
 Constructor. More...
 
 ~EcalEndcapDigisValidation ()
 Destructor. More...
 
- Public Member Functions inherited from edm::EDAnalyzer
 EDAnalyzer ()
 
std::string workerType () const
 
virtual ~EDAnalyzer ()
 
- Public Member Functions inherited from edm::EDConsumerBase
 EDConsumerBase ()
 
ProductHolderIndex indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductHolderIndex > &) const
 
void itemsToGet (BranchType, std::vector< ProductHolderIndex > &) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
 
virtual ~EDConsumerBase ()
 

Protected Member Functions

void analyze (edm::Event const &e, edm::EventSetup const &c)
 Analyze. More...
 
void beginRun (edm::Run const &, edm::EventSetup const &c)
 
void checkCalibrations (edm::EventSetup const &c)
 
void endJob (void)
 
- Protected Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
CurrentProcessingContext const * currentContext () const
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Private Types

typedef std::map< uint32_t,
float, std::less< uint32_t > > 		MapType
 

Private Attributes

double barrelADCtoGeV_
 
DQMStore * dbe_
 
edm::InputTag EEdigiCollection_
 
double endcapADCtoGeV_
 
std::map< int, double,
std::less< int > > 		gainConv_
 
MonitorElement * meEEDigiADCAnalog_ [10]
 
MonitorElement * meEEDigiADCg12_ [10]
 
MonitorElement * meEEDigiADCg1_ [10]
 
MonitorElement * meEEDigiADCg6_ [10]
 
MonitorElement * meEEDigiADCGlobal_
 
MonitorElement * meEEDigiADCgS_ [10]
 
MonitorElement * meEEDigiGain_ [10]
 
MonitorElement * meEEDigiMultiplicityzm_
 
MonitorElement * meEEDigiMultiplicityzp_
 
MonitorElement * meEEDigiOccupancyzm_
 
MonitorElement * meEEDigiOccupancyzp_
 
MonitorElement * meEEMaximumgt100ADC_
 
MonitorElement * meEEMaximumgt20ADC_
 
MonitorElement * meEEnADCafterSwitch_
 
MonitorElement * meEEPedestal_
 
std::string outputFile_
 
bool verbose_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
typedef WorkerT< EDAnalyzer > WorkerType
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 

Detailed Description

Definition at line 36 of file EcalEndcapDigisValidation.h.

Member Typedef Documentation

typedef std::map<uint32_t,float,std::less<uint32_t> > EcalEndcapDigisValidation::MapType
private

Definition at line 38 of file EcalEndcapDigisValidation.h.

Constructor & Destructor Documentation

EcalEndcapDigisValidation::EcalEndcapDigisValidation ( const edm::ParameterSet &  ps)

Constructor.

Definition at line 18 of file EcalEndcapDigisValidation.cc.

References barrelADCtoGeV_, dbe_, endcapADCtoGeV_, gainConv_, edm::ParameterSet::getUntrackedParameter(), timingPdfMaker::histo, i, meEEDigiADCAnalog_, meEEDigiADCg12_, meEEDigiADCg1_, meEEDigiADCg6_, meEEDigiADCGlobal_, meEEDigiADCgS_, meEEDigiGain_, meEEDigiMultiplicityzm_, meEEDigiMultiplicityzp_, meEEDigiOccupancyzm_, meEEDigiOccupancyzp_, meEEMaximumgt100ADC_, meEEMaximumgt20ADC_, meEEnADCafterSwitch_, meEEPedestal_, cppFunctionSkipper::operator, and verbose_.

18  :
19  EEdigiCollection_(ps.getParameter<edm::InputTag>("EEdigiCollection"))
20  {
21 
22  // verbosity switch
23  verbose_ = ps.getUntrackedParameter<bool>("verbose", false);
24 
25  dbe_ = 0;
26 
27  // get hold of back-end interface
28  dbe_ = Service<DQMStore>().operator->();
29 
30  if ( dbe_ ) {
31  if ( verbose_ ) {
32  dbe_->setVerbose(1);
33  } else {
34  dbe_->setVerbose(0);
35  }
36  }
37 
38  if ( dbe_ ) {
39  if ( verbose_ ) dbe_->showDirStructure();
40  }
41 
42  gainConv_[1] = 1.;
43  gainConv_[2] = 2.;
44  gainConv_[3] = 12.;
45  gainConv_[0] = 12.; // saturated channels
46  barrelADCtoGeV_ = 0.035;
47  endcapADCtoGeV_ = 0.06;
48 
49  meEEDigiOccupancyzp_ = 0;
50  meEEDigiOccupancyzm_ = 0;
51 
52  meEEDigiMultiplicityzp_ = 0;
53  meEEDigiMultiplicityzm_ = 0;
54 
55  meEEDigiADCGlobal_ = 0;
56 
57  for (int i = 0; i < 10 ; i++ ) {
58  meEEDigiADCAnalog_[i] = 0;
59  meEEDigiADCgS_[i] = 0;
60  meEEDigiADCg1_[i] = 0;
61  meEEDigiADCg6_[i] = 0;
62  meEEDigiADCg12_[i] = 0;
63  meEEDigiGain_[i] = 0;
64  }
65 
66  meEEPedestal_ = 0;
67 
68  meEEMaximumgt100ADC_ = 0;
69 
70  meEEMaximumgt20ADC_ = 0;
71 
72  meEEnADCafterSwitch_ = 0;
73 
74  Char_t histo[200];
75 
76 
77  if ( dbe_ ) {
78  dbe_->setCurrentFolder("EcalDigisV/EcalDigiTask");
79 
80  sprintf (histo, "EcalDigiTask Endcap occupancy z+" ) ;
81  meEEDigiOccupancyzp_ = dbe_->book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
82 
83  sprintf (histo, "EcalDigiTask Endcap occupancy z-" ) ;
84  meEEDigiOccupancyzm_ = dbe_->book2D(histo, histo, 100, 0., 100., 100, 0., 100.);
85 
86  sprintf (histo, "EcalDigiTask Endcap multiplicity z+" ) ;
87  meEEDigiMultiplicityzp_ = dbe_->book1D(histo, histo, 100, 0., 7324.);
88 
89  sprintf (histo, "EcalDigiTask Endcap multiplicity z-" ) ;
90  meEEDigiMultiplicityzm_ = dbe_->book1D(histo, histo, 100, 0., 7324.);
91 
92  sprintf (histo, "EcalDigiTask Endcap global pulse shape" ) ;
93  meEEDigiADCGlobal_ = dbe_->bookProfile(histo, histo, 10, 0, 10, 10000, 0., 1000.) ;
94 
95  for (int i = 0; i < 10 ; i++ ) {
96 
97  sprintf (histo, "EcalDigiTask Endcap analog pulse %02d", i+1) ;
98  meEEDigiADCAnalog_[i] = dbe_->book1D(histo, histo, 4000, 0., 400.);
99 
100  sprintf (histo, "EcalDigiTask Endcap ADC pulse %02d Gain 0 - Saturated", i+1) ;
101  meEEDigiADCgS_[i] = dbe_->book1D(histo, histo, 4096, -0.5, 4095.5);
102 
103  sprintf (histo, "EcalDigiTask Endcap ADC pulse %02d Gain 1", i+1) ;
104  meEEDigiADCg1_[i] = dbe_->book1D(histo, histo, 4096, -0.5, 4095.5);
105 
106  sprintf (histo, "EcalDigiTask Endcap ADC pulse %02d Gain 6", i+1) ;
107  meEEDigiADCg6_[i] = dbe_->book1D(histo, histo, 4096, -0.5, 4095.5);
108 
109  sprintf (histo, "EcalDigiTask Endcap ADC pulse %02d Gain 12", i+1) ;
110  meEEDigiADCg12_[i] = dbe_->book1D(histo, histo, 4096, -0.5, 4095.5);
111 
112  sprintf (histo, "EcalDigiTask Endcap gain pulse %02d", i+1) ;
113  meEEDigiGain_[i] = dbe_->book1D(histo, histo, 4, 0, 4);
114  }
115 
116  sprintf (histo, "EcalDigiTask Endcap pedestal for pre-sample" ) ;
117  meEEPedestal_ = dbe_->book1D(histo, histo, 4096, -0.5, 4095.5) ;
118 
119  sprintf (histo, "EcalDigiTask Endcap maximum position gt 100 ADC" ) ;
120  meEEMaximumgt100ADC_ = dbe_->book1D(histo, histo, 10, 0., 10.) ;
121 
122  sprintf (histo, "EcalDigiTask Endcap maximum position gt 20 ADC" ) ;
123  meEEMaximumgt20ADC_ = dbe_->book1D(histo, histo, 10, 0., 10.) ;
124 
125  sprintf (histo, "EcalDigiTask Endcap ADC counts after gain switch" ) ;
126  meEEnADCafterSwitch_ = dbe_->book1D(histo, histo, 10, 0., 10.) ;
127 
128  }
129 
130 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:722
EcalEndcapDigisValidation::meEEDigiADCAnalog_
MonitorElement * meEEDigiADCAnalog_[10]
Definition: EcalEndcapDigisValidation.h:85
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:1036
DQMStore::setVerbose
void setVerbose(unsigned level)
Definition: DQMStore.cc:398
EcalEndcapDigisValidation::meEEDigiADCg12_
MonitorElement * meEEDigiADCg12_[10]
Definition: EcalEndcapDigisValidation.h:90
EcalEndcapDigisValidation::gainConv_
std::map< int, double, std::less< int > > gainConv_
Definition: EcalEndcapDigisValidation.h:72
DQMStore::showDirStructure
void showDirStructure(void) const
Definition: DQMStore.cc:2766
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:850
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:434
EcalEndcapDigisValidation::~EcalEndcapDigisValidation ( )

Destructor.

Definition at line 132 of file EcalEndcapDigisValidation.cc.

132  {
133 
134 }

Member Function Documentation

void EcalEndcapDigisValidation::analyze ( edm::Event const &  e,
edm::EventSetup const &  c 
)
protectedvirtual

Analyze.

Implements edm::EDAnalyzer.

Definition at line 146 of file EcalEndcapDigisValidation.cc.

References EcalMGPASample::adc(), EEdigiCollection_, endcapADCtoGeV_, MonitorElement::Fill(), gainConv_, EcalMGPASample::gainId(), edm::Event::getByLabel(), i, EEDataFrame::id(), edm::HandleBase::isValid(), EEDetId::ix(), EEDetId::iy(), LogDebug, EcalDataFrame::MAXSAMPLES, meEEDigiADCAnalog_, meEEDigiADCg12_, meEEDigiADCg1_, meEEDigiADCg6_, meEEDigiADCGlobal_, meEEDigiADCgS_, meEEDigiGain_, meEEDigiMultiplicityzm_, meEEDigiMultiplicityzp_, meEEDigiOccupancyzm_, meEEDigiOccupancyzp_, meEEMaximumgt100ADC_, meEEMaximumgt20ADC_, meEEnADCafterSwitch_, meEEPedestal_, ecalTB2006H4_GenSimDigiReco_cfg::mySample, edm::Handle< T >::product(), DetId::rawId(), compare_using_db::sample, EcalDataFrame::size(), and EEDetId::zside().

146  {
147 
148  //LogInfo("EventInfo") << " Run = " << e.id().run() << " Event = " << e.id().event();
149 
150  Handle<EEDigiCollection> EcalDigiEE;
151 
152  e.getByLabel( EEdigiCollection_ , EcalDigiEE );
153 
154  // Return if no Endcap data available
155  if( !EcalDigiEE.isValid() ) return;
156 
157  // ENDCAP
158 
159  // loop over Digis
160 
161  const EEDigiCollection * endcapDigi = EcalDigiEE.product () ;
162 
163  std::vector<double> eeAnalogSignal ;
164  std::vector<double> eeADCCounts ;
165  std::vector<double> eeADCGains ;
166  eeAnalogSignal.reserve(EEDataFrame::MAXSAMPLES);
167  eeADCCounts.reserve(EEDataFrame::MAXSAMPLES);
168  eeADCGains.reserve(EEDataFrame::MAXSAMPLES);
169 
170  int nDigiszp = 0;
171  int nDigiszm = 0;
172 
173  for (unsigned int digis=0; digis<EcalDigiEE->size(); ++digis) {
174 
175  EEDataFrame eedf=(*endcapDigi)[digis];
176  int nrSamples=eedf.size();
177 
178  EEDetId eeid = eedf.id () ;
179 
180  if (eeid.zside() > 0 ) {
181  if (meEEDigiOccupancyzp_) meEEDigiOccupancyzp_->Fill( eeid.ix(), eeid.iy() );
182  nDigiszp++;
183  }
184  else if (eeid.zside() < 0 ) {
185  if (meEEDigiOccupancyzm_) meEEDigiOccupancyzm_->Fill( eeid.ix(), eeid.iy() );
186  nDigiszm++;
187  }
188 
189  double Emax = 0. ;
190  int Pmax = 0 ;
191  double pedestalPreSample = 0.;
192  double pedestalPreSampleAnalog = 0.;
193  int countsAfterGainSwitch = -1;
194  double higherGain = 1.;
195  int higherGainSample = 0;
196 
197  for (int sample = 0 ; sample < nrSamples; ++sample) {
198  eeAnalogSignal[sample] = 0.;
199  eeADCCounts[sample] = 0.;
200  eeADCGains[sample] = 0.;
201  }
202 
203  for (int sample = 0 ; sample < nrSamples; ++sample) {
204 
205  EcalMGPASample mySample = eedf[sample];
206 
207  eeADCCounts[sample] = (mySample.adc());
208  eeADCGains[sample] = (mySample.gainId()) ;
209  eeAnalogSignal[sample] = (eeADCCounts[sample]*gainConv_[(int)eeADCGains[sample]]*endcapADCtoGeV_);
210 
211  if (Emax < eeAnalogSignal[sample] ) {
212  Emax = eeAnalogSignal[sample] ;
213  Pmax = sample ;
214  }
215 
216  if ( sample < 3 ) {
217  pedestalPreSample += eeADCCounts[sample] ;
218  pedestalPreSampleAnalog += eeADCCounts[sample]*gainConv_[(int)eeADCGains[sample]]*endcapADCtoGeV_ ;
219  }
220 
221  if (sample > 0 && ( ((eeADCGains[sample] > eeADCGains[sample-1]) && (eeADCGains[sample-1]!=0)) || (countsAfterGainSwitch<0 && eeADCGains[sample]==0)) ){
222  higherGain = eeADCGains[sample];
223  higherGainSample = sample;
224  countsAfterGainSwitch = 1;
225  }
226 
227  if ( (higherGain > 1 && (higherGainSample != sample) && (eeADCGains[sample] == higherGain)) || (higherGain==3 && (higherGainSample != sample) && (eeADCGains[sample]==0)) || (higherGain==0 && (higherGainSample != sample) && ((eeADCGains[sample]==0) || (eeADCGains[sample]==3))) ) countsAfterGainSwitch++ ;
228  }
229  pedestalPreSample /= 3. ;
230  pedestalPreSampleAnalog /= 3. ;
231 
232  LogDebug("DigiInfo") << "Endcap Digi for EEDetId = " << eeid.rawId() << " x,y " << eeid.ix() << " " << eeid.iy() ;
233  for ( int i = 0; i < 10 ; i++ ) {
234  LogDebug("DigiInfo") << "sample " << i << " ADC = " << eeADCCounts[i] << " gain = " << eeADCGains[i] << " Analog = " << eeAnalogSignal[i] ;
235  }
236  LogDebug("DigiInfo") << "Maximum energy = " << Emax << " in sample " << Pmax << " Pedestal from pre-sample = " << pedestalPreSampleAnalog;
237  if ( countsAfterGainSwitch > 0 ) LogDebug("DigiInfo") << "Counts after switch " << countsAfterGainSwitch;
238 
239  if ( countsAfterGainSwitch > 0 && countsAfterGainSwitch < 5 ) {
240  edm::LogWarning("DigiWarning") << "Wrong number of counts after gain switch before next switch! " << countsAfterGainSwitch ;
241  for ( int i = 0; i < 10 ; i++ ) {
242  edm::LogWarning("DigiWarning") << "sample " << i << " ADC = " << eeADCCounts[i] << " gain = " << eeADCGains[i] << " Analog = " << eeAnalogSignal[i];
243  }
244  }
245 
246  for ( int i = 0 ; i < 10 ; i++ ) {
247  if (meEEDigiADCGlobal_ && (Emax-pedestalPreSampleAnalog*gainConv_[(int)eeADCGains[Pmax]]) > 100.*endcapADCtoGeV_) meEEDigiADCGlobal_->Fill( i , eeAnalogSignal[i] ) ;
248  if (meEEDigiADCAnalog_[i]) meEEDigiADCAnalog_[i]->Fill( eeAnalogSignal[i] ) ;
249  if ( eeADCGains[i] == 0 ) {
250  if (meEEDigiADCgS_[i]) meEEDigiADCgS_[i]->Fill( eeADCCounts[i] ) ;
251  }
252  else if ( eeADCGains[i] == 3 ) {
253  if (meEEDigiADCg1_[i]) meEEDigiADCg1_[i]->Fill( eeADCCounts[i] ) ;
254  }
255  else if ( eeADCGains[i] == 2 ) {
256  if (meEEDigiADCg6_[i]) meEEDigiADCg6_[i]->Fill( eeADCCounts[i] ) ;
257  }
258  else if ( eeADCGains[i] == 1 ) {
259  if (meEEDigiADCg12_[i]) meEEDigiADCg12_[i]->Fill( eeADCCounts[i] ) ;
260  }
261  if (meEEDigiGain_[i]) meEEDigiGain_[i]->Fill( eeADCGains[i] ) ;
262  }
263 
264  if (meEEPedestal_) meEEPedestal_->Fill ( pedestalPreSample ) ;
265  if (meEEMaximumgt20ADC_ && (Emax-pedestalPreSampleAnalog*gainConv_[(int)eeADCGains[Pmax]]) > 20.*endcapADCtoGeV_) meEEMaximumgt20ADC_->Fill( Pmax ) ;
266  if (meEEMaximumgt100ADC_ && (Emax-pedestalPreSampleAnalog*gainConv_[(int)eeADCGains[Pmax]]) > 100.*endcapADCtoGeV_) meEEMaximumgt100ADC_->Fill( Pmax ) ;
267  if (meEEnADCafterSwitch_) meEEnADCafterSwitch_->Fill(countsAfterGainSwitch);
268 
269  }
270 
271  if ( meEEDigiMultiplicityzp_ ) meEEDigiMultiplicityzp_->Fill(nDigiszp);
272  if ( meEEDigiMultiplicityzm_ ) meEEDigiMultiplicityzm_->Fill(nDigiszm);
273 
274 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
i
int i
Definition: DBlmapReader.cc:9
EEDetId::ix
int ix() const
Definition: EEDetId.h:77
EcalMGPASample::gainId
int gainId() const
get the gainId (2 bits)
Definition: EcalMGPASample.h:37
EcalDataFrame::size
int size() const
Definition: EcalDataFrame.h:27
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:45
EcalEndcapDigisValidation::meEEDigiADCAnalog_
MonitorElement * meEEDigiADCAnalog_[10]
Definition: EcalEndcapDigisValidation.h:85
EEDetId::zside
int zside() const
Definition: EEDetId.h:71
EEDetId::iy
int iy() const
Definition: EEDetId.h:83
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
EcalEndcapDigisValidation::meEEDigiADCg12_
MonitorElement * meEEDigiADCg12_[10]
Definition: EcalEndcapDigisValidation.h:90
EEDataFrame::id
key_type id() const
Definition: EEDataFrame.h:29
edm::Handle::product
T const * product() const
Definition: Handle.h:74
ecalTB2006H4_GenSimDigiReco_cfg.mySample
string mySample
Definition: ecalTB2006H4_GenSimDigiReco_cfg.py:4
EcalEndcapDigisValidation::gainConv_
std::map< int, double, std::less< int > > gainConv_
Definition: EcalEndcapDigisValidation.h:72
EcalDataFrame::MAXSAMPLES
static const int MAXSAMPLES
Definition: EcalDataFrame.h:49
EcalMGPASample::adc
int adc() const
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:35
compare_using_db.sample
string sample
Definition: compare_using_db.py:29
void EcalEndcapDigisValidation::beginRun ( edm::Run const &  ,
edm::EventSetup const &  c 
)
protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 136 of file EcalEndcapDigisValidation.cc.

References checkCalibrations().

136  {
137 
138  checkCalibrations(c);
139 
140 }
EcalEndcapDigisValidation::checkCalibrations
void checkCalibrations(edm::EventSetup const &c)
Definition: EcalEndcapDigisValidation.cc:276
void EcalEndcapDigisValidation::checkCalibrations ( edm::EventSetup const &  c)
protected

Definition at line 276 of file EcalEndcapDigisValidation.cc.

References barrelADCtoGeV_, endcapADCtoGeV_, EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), gainConv_, edm::EventSetup::get(), EcalADCToGeVConstant::getEBValue(), EcalADCToGeVConstant::getEEValue(), LogDebug, and edm::ESHandle< class >::product().

Referenced by beginRun().

277 {
278 
279  // ADC -> GeV Scale
280  edm::ESHandle<EcalADCToGeVConstant> pAgc;
281  eventSetup.get<EcalADCToGeVConstantRcd>().get(pAgc);
282  const EcalADCToGeVConstant* agc = pAgc.product();
283 
284  EcalMGPAGainRatio * defaultRatios = new EcalMGPAGainRatio();
285 
286  gainConv_[1] = 1.;
287  gainConv_[2] = defaultRatios->gain12Over6() ;
288  gainConv_[3] = gainConv_[2]*(defaultRatios->gain6Over1()) ;
289  gainConv_[0] = gainConv_[2]*(defaultRatios->gain6Over1()) ;
290 
291  LogDebug("EcalDigi") << " Gains conversions: " << "\n" << " g0 = " << gainConv_[0] << "\n" << " g1 = " << gainConv_[1] << "\n" << " g2 = " << gainConv_[2] << "\n" << " g3 = " << gainConv_[3];
292 
293  delete defaultRatios;
294 
295  const double barrelADCtoGeV_ = agc->getEBValue();
296  LogDebug("EcalDigi") << " Barrel GeV/ADC = " << barrelADCtoGeV_;
297  const double endcapADCtoGeV_ = agc->getEEValue();
298  LogDebug("EcalDigi") << " Endcap GeV/ADC = " << endcapADCtoGeV_;
299 
300 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
EcalMGPAGainRatio::gain6Over1
float gain6Over1() const
Definition: EcalMGPAGainRatio.h:19
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
EcalMGPAGainRatio::gain12Over6
float gain12Over6() const
Definition: EcalMGPAGainRatio.h:18
EcalEndcapDigisValidation::gainConv_
std::map< int, double, std::less< int > > gainConv_
Definition: EcalEndcapDigisValidation.h:72
void EcalEndcapDigisValidation::endJob ( void  )
protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 142 of file EcalEndcapDigisValidation.cc.

142  {
143 
144 }

Member Data Documentation

double EcalEndcapDigisValidation::barrelADCtoGeV_
private

Definition at line 74 of file EcalEndcapDigisValidation.h.

Referenced by checkCalibrations(), and EcalEndcapDigisValidation().

DQMStore* EcalEndcapDigisValidation::dbe_
private

Definition at line 66 of file EcalEndcapDigisValidation.h.

Referenced by EcalEndcapDigisValidation().

edm::InputTag EcalEndcapDigisValidation::EEdigiCollection_
private

Definition at line 70 of file EcalEndcapDigisValidation.h.

Referenced by analyze().

double EcalEndcapDigisValidation::endcapADCtoGeV_
private

Definition at line 75 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), checkCalibrations(), and EcalEndcapDigisValidation().

std::map<int, double, std::less<int> > EcalEndcapDigisValidation::gainConv_
private

Definition at line 72 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), checkCalibrations(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiADCAnalog_[10]
private

Definition at line 85 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiADCg12_[10]
private

Definition at line 90 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiADCg1_[10]
private

Definition at line 88 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiADCg6_[10]
private

Definition at line 89 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiADCGlobal_
private

Definition at line 83 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiADCgS_[10]
private

Definition at line 87 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiGain_[10]
private

Definition at line 92 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiMultiplicityzm_
private

Definition at line 81 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiMultiplicityzp_
private

Definition at line 80 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiOccupancyzm_
private

Definition at line 78 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEDigiOccupancyzp_
private

Definition at line 77 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEMaximumgt100ADC_
private

Definition at line 96 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEMaximumgt20ADC_
private

Definition at line 98 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEnADCafterSwitch_
private

Definition at line 100 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

MonitorElement* EcalEndcapDigisValidation::meEEPedestal_
private

Definition at line 94 of file EcalEndcapDigisValidation.h.

Referenced by analyze(), and EcalEndcapDigisValidation().

std::string EcalEndcapDigisValidation::outputFile_
private

Definition at line 68 of file EcalEndcapDigisValidation.h.

bool EcalEndcapDigisValidation::verbose_
private

Definition at line 64 of file EcalEndcapDigisValidation.h.

Referenced by EcalEndcapDigisValidation().

Generated for CMSSW Reference Manual by doxygen 1.8.5