CMS 3D CMS Logo
List of all members | Public Member Functions | Protected Member Functions | Protected Attributes
EcalUncalibRecHitWorkerGlobal Class Reference

#include <EcalUncalibRecHitWorkerGlobal.h>

Inheritance diagram for EcalUncalibRecHitWorkerGlobal:
EcalUncalibRecHitWorkerRunOneDigiBase EcalUncalibRecHitWorkerBaseClass

Public Member Functions

 EcalUncalibRecHitWorkerGlobal (const edm::ParameterSet &, edm::ConsumesCollector &c)
 
 EcalUncalibRecHitWorkerGlobal (const edm::ParameterSet &)
 
 EcalUncalibRecHitWorkerGlobal ()
 
edm::ParameterSetDescription getAlgoDescription () override
 
bool run (const edm::Event &evt, const EcalDigiCollection::const_iterator &digi, EcalUncalibratedRecHitCollection &result) override
 
void set (const edm::EventSetup &es) override
 
 ~EcalUncalibRecHitWorkerGlobal () override
 
- Public Member Functions inherited from EcalUncalibRecHitWorkerRunOneDigiBase
 EcalUncalibRecHitWorkerRunOneDigiBase (const edm::ParameterSet &p, edm::ConsumesCollector &c)
 
 EcalUncalibRecHitWorkerRunOneDigiBase (const edm::ParameterSet &p)
 
 EcalUncalibRecHitWorkerRunOneDigiBase ()
 
void run (const edm::Event &evt, const EcalDigiCollection &digis, EcalUncalibratedRecHitCollection &result) override
 
 ~EcalUncalibRecHitWorkerRunOneDigiBase () override
 
- Public Member Functions inherited from EcalUncalibRecHitWorkerBaseClass
 EcalUncalibRecHitWorkerBaseClass (const edm::ParameterSet &, edm::ConsumesCollector &c)
 
 EcalUncalibRecHitWorkerBaseClass (const edm::ParameterSet &)
 
 EcalUncalibRecHitWorkerBaseClass ()
 
virtual void set (const edm::Event &evt)
 
virtual ~EcalUncalibRecHitWorkerBaseClass ()
 

Protected Member Functions

template<class C >
int isSaturated (const C &digi)
 
double timeCorrection (float ampli, const std::vector< float > &amplitudeBins, const std::vector< float > &shiftBins)
 

Protected Attributes

double amplitudeThreshEB_
 
double amplitudeThreshEE_
 
const EcalWeightSet::EcalChi2WeightMatrix * chi2mat [2]
 
double chi2ThreshEB_
 
double chi2ThreshEE_
 
std::vector< double > EBamplitudeFitParameters_
 
std::vector< double > EBchi2Parameters_
 
std::vector< double > ebPulseShape_
 
double ebSpikeThresh_
 
double EBtimeConstantTerm_
 
std::pair< double, double > EBtimeFitLimits_
 
std::vector< double > EBtimeFitParameters_
 
double EBtimeNconst_
 
std::vector< double > EEamplitudeFitParameters_
 
std::vector< double > EEchi2Parameters_
 
std::vector< double > eePulseShape_
 
double EEtimeConstantTerm_
 
std::pair< double, double > EEtimeFitLimits_
 
std::vector< double > EEtimeFitParameters_
 
double EEtimeNconst_
 
double gainRatios [3]
 
edm::ESHandle< EcalGainRatios > gains
 
edm::ESHandle< EcalWeightXtalGroups > grps
 
edm::ESHandle< EcalTimeCalibConstants > itime
 
bool kPoorRecoFlagEB_
 
bool kPoorRecoFlagEE_
 
edm::ESHandle< EcalTimeOffsetConstant > offtime
 
double outOfTimeThreshG12mEB_
 
double outOfTimeThreshG12mEE_
 
double outOfTimeThreshG12pEB_
 
double outOfTimeThreshG12pEE_
 
double outOfTimeThreshG61mEB_
 
double outOfTimeThreshG61mEE_
 
double outOfTimeThreshG61pEB_
 
double outOfTimeThreshG61pEE_
 
double pedRMSVec [3]
 
edm::ESHandle< EcalPedestals > peds
 
double pedVec [3]
 
EcalUncalibRecHitRatioMethodAlgo< EBDataFrame > ratioMethod_barrel_
 
EcalUncalibRecHitRatioMethodAlgo< EEDataFrame > ratioMethod_endcap_
 
edm::ESHandle< EcalSampleMask > sampleMaskHand_
 
EBShape testbeamEBShape
 
EEShape testbeamEEShape
 
edm::ESHandle< EcalTimeBiasCorrections > timeCorrBias_
 
const EcalWeightSet::EcalWeightMatrix * weights [2]
 
EcalUncalibRecHitRecWeightsAlgo< EBDataFrame > weightsMethod_barrel_
 
EcalUncalibRecHitRecWeightsAlgo< EEDataFrame > weightsMethod_endcap_
 
edm::ESHandle< EcalTBWeights > wgts
 

Detailed Description

Definition at line 36 of file EcalUncalibRecHitWorkerGlobal.h.

Constructor & Destructor Documentation

EcalUncalibRecHitWorkerGlobal::EcalUncalibRecHitWorkerGlobal ( const edm::ParameterSet &  ps,
edm::ConsumesCollector &  c 
)

Definition at line 20 of file EcalUncalibRecHitWorkerGlobal.cc.

References amplitudeThreshEB_, amplitudeThreshEE_, chi2ThreshEB_, chi2ThreshEE_, EBamplitudeFitParameters_, EBchi2Parameters_, ebPulseShape_, ebSpikeThresh_, EBtimeConstantTerm_, EBtimeFitLimits_, EBtimeFitParameters_, EBtimeNconst_, EEamplitudeFitParameters_, EEchi2Parameters_, eePulseShape_, EEtimeConstantTerm_, EEtimeFitLimits_, EEtimeFitParameters_, EEtimeNconst_, edm::ParameterSet::getParameter(), kPoorRecoFlagEB_, kPoorRecoFlagEE_, outOfTimeThreshG12mEB_, outOfTimeThreshG12mEE_, outOfTimeThreshG12pEB_, outOfTimeThreshG12pEE_, outOfTimeThreshG61mEB_, outOfTimeThreshG61mEE_, outOfTimeThreshG61pEB_, and outOfTimeThreshG61pEE_.

20  :
21  EcalUncalibRecHitWorkerRunOneDigiBase(ps,c), testbeamEEShape(EEShape(true)), testbeamEBShape(EBShape(true))
22 {
23  // ratio method parameters
24  EBtimeFitParameters_ = ps.getParameter<std::vector<double> >("EBtimeFitParameters");
25  EEtimeFitParameters_ = ps.getParameter<std::vector<double> >("EEtimeFitParameters");
26  EBamplitudeFitParameters_ = ps.getParameter<std::vector<double> >("EBamplitudeFitParameters");
27  EEamplitudeFitParameters_ = ps.getParameter<std::vector<double> >("EEamplitudeFitParameters");
28  EBtimeFitLimits_.first = ps.getParameter<double>("EBtimeFitLimits_Lower");
29  EBtimeFitLimits_.second = ps.getParameter<double>("EBtimeFitLimits_Upper");
30  EEtimeFitLimits_.first = ps.getParameter<double>("EEtimeFitLimits_Lower");
31  EEtimeFitLimits_.second = ps.getParameter<double>("EEtimeFitLimits_Upper");
32  EBtimeConstantTerm_=ps.getParameter<double>("EBtimeConstantTerm");
33  EBtimeNconst_=ps.getParameter<double>("EBtimeNconst");
34  EEtimeConstantTerm_=ps.getParameter<double>("EEtimeConstantTerm");
35  EEtimeNconst_=ps.getParameter<double>("EEtimeNconst");
36  outOfTimeThreshG12pEB_ = ps.getParameter<double>("outOfTimeThresholdGain12pEB");
37  outOfTimeThreshG12mEB_ = ps.getParameter<double>("outOfTimeThresholdGain12mEB");
38  outOfTimeThreshG61pEB_ = ps.getParameter<double>("outOfTimeThresholdGain61pEB");
39  outOfTimeThreshG61mEB_ = ps.getParameter<double>("outOfTimeThresholdGain61mEB");
40  outOfTimeThreshG12pEE_ = ps.getParameter<double>("outOfTimeThresholdGain12pEE");
41  outOfTimeThreshG12mEE_ = ps.getParameter<double>("outOfTimeThresholdGain12mEE");
42  outOfTimeThreshG61pEE_ = ps.getParameter<double>("outOfTimeThresholdGain61pEE");
43  outOfTimeThreshG61mEE_ = ps.getParameter<double>("outOfTimeThresholdGain61mEE");
44  amplitudeThreshEB_ = ps.getParameter<double>("amplitudeThresholdEB");
45  amplitudeThreshEE_ = ps.getParameter<double>("amplitudeThresholdEE");
46 
47  // spike threshold
48  ebSpikeThresh_ = ps.getParameter<double>("ebSpikeThreshold");
49 
50  ebPulseShape_ = ps.getParameter<std::vector<double> >("ebPulseShape");
51  eePulseShape_ = ps.getParameter<std::vector<double> >("eePulseShape");
52 
53  // chi2 parameters
54  kPoorRecoFlagEB_ = ps.getParameter<bool>("kPoorRecoFlagEB");
55  kPoorRecoFlagEE_ = ps.getParameter<bool>("kPoorRecoFlagEE");;
56  chi2ThreshEB_=ps.getParameter<double>("chi2ThreshEB_");
57  chi2ThreshEE_=ps.getParameter<double>("chi2ThreshEE_");
58  EBchi2Parameters_ = ps.getParameter<std::vector<double> >("EBchi2Parameters");
59  EEchi2Parameters_ = ps.getParameter<std::vector<double> >("EEchi2Parameters");
60 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
EcalUncalibRecHitWorkerGlobal::EBtimeFitLimits_
std::pair< double, double > EBtimeFitLimits_
Definition: EcalUncalibRecHitWorkerGlobal.h:80
EcalUncalibRecHitWorkerGlobal::EEtimeFitLimits_
std::pair< double, double > EEtimeFitLimits_
Definition: EcalUncalibRecHitWorkerGlobal.h:81
EBShape
Definition: EBShape.h:6
EEShape
Definition: EEShape.h:6
EcalUncalibRecHitWorkerGlobal::EcalUncalibRecHitWorkerGlobal ( const edm::ParameterSet &  ps)

Definition at line 64 of file EcalUncalibRecHitWorkerGlobal.cc.

References amplitudeThreshEB_, amplitudeThreshEE_, chi2ThreshEB_, chi2ThreshEE_, EBamplitudeFitParameters_, EBchi2Parameters_, ebPulseShape_, ebSpikeThresh_, EBtimeConstantTerm_, EBtimeFitLimits_, EBtimeFitParameters_, EBtimeNconst_, EEamplitudeFitParameters_, EEchi2Parameters_, eePulseShape_, EEtimeConstantTerm_, EEtimeFitLimits_, EEtimeFitParameters_, EEtimeNconst_, edm::ParameterSet::getParameter(), kPoorRecoFlagEB_, kPoorRecoFlagEE_, outOfTimeThreshG12mEB_, outOfTimeThreshG12mEE_, outOfTimeThreshG12pEB_, outOfTimeThreshG12pEE_, outOfTimeThreshG61mEB_, outOfTimeThreshG61mEE_, outOfTimeThreshG61pEB_, and outOfTimeThreshG61pEE_.

64  :
65  EcalUncalibRecHitWorkerRunOneDigiBase(ps), testbeamEEShape(EEShape(true)), testbeamEBShape(EBShape(true))
66 {
67  // ratio method parameters
68  EBtimeFitParameters_ = ps.getParameter<std::vector<double> >("EBtimeFitParameters");
69  EEtimeFitParameters_ = ps.getParameter<std::vector<double> >("EEtimeFitParameters");
70  EBamplitudeFitParameters_ = ps.getParameter<std::vector<double> >("EBamplitudeFitParameters");
71  EEamplitudeFitParameters_ = ps.getParameter<std::vector<double> >("EEamplitudeFitParameters");
72  EBtimeFitLimits_.first = ps.getParameter<double>("EBtimeFitLimits_Lower");
73  EBtimeFitLimits_.second = ps.getParameter<double>("EBtimeFitLimits_Upper");
74  EEtimeFitLimits_.first = ps.getParameter<double>("EEtimeFitLimits_Lower");
75  EEtimeFitLimits_.second = ps.getParameter<double>("EEtimeFitLimits_Upper");
76  EBtimeConstantTerm_=ps.getParameter<double>("EBtimeConstantTerm");
77  EBtimeNconst_=ps.getParameter<double>("EBtimeNconst");
78  EEtimeConstantTerm_=ps.getParameter<double>("EEtimeConstantTerm");
79  EEtimeNconst_=ps.getParameter<double>("EEtimeNconst");
80  outOfTimeThreshG12pEB_ = ps.getParameter<double>("outOfTimeThresholdGain12pEB");
81  outOfTimeThreshG12mEB_ = ps.getParameter<double>("outOfTimeThresholdGain12mEB");
82  outOfTimeThreshG61pEB_ = ps.getParameter<double>("outOfTimeThresholdGain61pEB");
83  outOfTimeThreshG61mEB_ = ps.getParameter<double>("outOfTimeThresholdGain61mEB");
84  outOfTimeThreshG12pEE_ = ps.getParameter<double>("outOfTimeThresholdGain12pEE");
85  outOfTimeThreshG12mEE_ = ps.getParameter<double>("outOfTimeThresholdGain12mEE");
86  outOfTimeThreshG61pEE_ = ps.getParameter<double>("outOfTimeThresholdGain61pEE");
87  outOfTimeThreshG61mEE_ = ps.getParameter<double>("outOfTimeThresholdGain61mEE");
88  amplitudeThreshEB_ = ps.getParameter<double>("amplitudeThresholdEB");
89  amplitudeThreshEE_ = ps.getParameter<double>("amplitudeThresholdEE");
90 
91  // spike threshold
92  ebSpikeThresh_ = ps.getParameter<double>("ebSpikeThreshold");
93 
94  ebPulseShape_ = ps.getParameter<std::vector<double> >("ebPulseShape");
95  eePulseShape_ = ps.getParameter<std::vector<double> >("eePulseShape");
96  // chi2 parameters
97  kPoorRecoFlagEB_ = ps.getParameter<bool>("kPoorRecoFlagEB");
98  kPoorRecoFlagEE_ = ps.getParameter<bool>("kPoorRecoFlagEE");;
99  chi2ThreshEB_=ps.getParameter<double>("chi2ThreshEB_");
100  chi2ThreshEE_=ps.getParameter<double>("chi2ThreshEE_");
101  EBchi2Parameters_ = ps.getParameter<std::vector<double> >("EBchi2Parameters");
102  EEchi2Parameters_ = ps.getParameter<std::vector<double> >("EEchi2Parameters");
103 }
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
EcalUncalibRecHitWorkerGlobal::EBtimeFitLimits_
std::pair< double, double > EBtimeFitLimits_
Definition: EcalUncalibRecHitWorkerGlobal.h:80
EcalUncalibRecHitWorkerGlobal::EEtimeFitLimits_
std::pair< double, double > EEtimeFitLimits_
Definition: EcalUncalibRecHitWorkerGlobal.h:81
EBShape
Definition: EBShape.h:6
EEShape
Definition: EEShape.h:6
EcalUncalibRecHitWorkerGlobal::EcalUncalibRecHitWorkerGlobal ( )
inline

Definition at line 41 of file EcalUncalibRecHitWorkerGlobal.h.

41 :testbeamEEShape(EEShape(true)), testbeamEBShape(EBShape(true)){;}
EBShape
Definition: EBShape.h:6
EEShape
Definition: EEShape.h:6
EcalUncalibRecHitWorkerGlobal::~EcalUncalibRecHitWorkerGlobal ( )
inlineoverride

Definition at line 42 of file EcalUncalibRecHitWorkerGlobal.h.

References mps_fire::result, and writedatasetfile::run.

42 {};

Member Function Documentation

edm::ParameterSetDescription EcalUncalibRecHitWorkerGlobal::getAlgoDescription ( )
overridevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 506 of file EcalUncalibRecHitWorkerGlobal.cc.

References edm::ParameterSetDescription::addNode(), and DEFINE_EDM_PLUGIN.

506  {
507 
508  edm::ParameterSetDescription psd;
509  psd.addNode(edm::ParameterDescription<std::vector<double>>("eePulseShape", {5.2e-05,-5.26e-05,6.66e-05,0.1168,0.7575,1.0,0.8876,0.6732,0.4741,0.3194}, true) and
510  edm::ParameterDescription<std::vector<double>>("EBtimeFitParameters", {-2.015452,3.130702,-12.3473,41.88921,-82.83944,91.01147,-50.35761,11.05621}, true) and
511  edm::ParameterDescription<double>("outOfTimeThresholdGain61pEB", 5, true) and
512  edm::ParameterDescription<double>("amplitudeThresholdEE", 10, true) and
513  edm::ParameterDescription<double>("EBtimeConstantTerm", 0.6, true) and
514  edm::ParameterDescription<double>("outOfTimeThresholdGain61pEE", 1000, true) and
515  edm::ParameterDescription<double>("ebSpikeThreshold", 1.042, true) and
516  edm::ParameterDescription<double>("EBtimeNconst", 28.5, true) and
517  edm::ParameterDescription<bool>("kPoorRecoFlagEB", true, true) and
518  edm::ParameterDescription<std::vector<double>>("ebPulseShape", {5.2e-05,-5.26e-05,6.66e-05,0.1168,0.7575,1.0,0.8876,0.6732,0.4741,0.3194}, true) and
519  edm::ParameterDescription<double>("EBtimeFitLimits_Lower", 0.2, true) and
520  edm::ParameterDescription<bool>("kPoorRecoFlagEE", false, true) and
521  edm::ParameterDescription<double>("chi2ThreshEB_", 36.0, true) and
522  edm::ParameterDescription<std::vector<double>>("EEtimeFitParameters", {-2.390548,3.553628,-17.62341,67.67538,-133.213,140.7432,-75.41106,16.20277}, true) and
523  edm::ParameterDescription<double>("outOfTimeThresholdGain61mEE", 1000, true) and
524  edm::ParameterDescription<std::vector<double>>("EEchi2Parameters", {2.122,0.022,2.122,0.022}, true) and
525  edm::ParameterDescription<double>("outOfTimeThresholdGain12mEE", 1000, true) and
526  edm::ParameterDescription<double>("outOfTimeThresholdGain12mEB", 5, true) and
527  edm::ParameterDescription<double>("EEtimeFitLimits_Upper", 1.4, true) and
528  edm::ParameterDescription<double>("EEtimeFitLimits_Lower", 0.2, true) and
529  edm::ParameterDescription<std::vector<double>>("EEamplitudeFitParameters", {1.89,1.4}, true) and
530  edm::ParameterDescription<std::vector<double>>("EBamplitudeFitParameters", {1.138,1.652}, true) and
531  edm::ParameterDescription<double>("amplitudeThresholdEB", 10, true) and
532  edm::ParameterDescription<double>("outOfTimeThresholdGain12pEE", 1000, true) and
533  edm::ParameterDescription<double>("outOfTimeThresholdGain12pEB", 5, true) and
534  edm::ParameterDescription<double>("EEtimeNconst", 31.8, true) and
535  edm::ParameterDescription<double>("outOfTimeThresholdGain61mEB", 5, true) and
536  edm::ParameterDescription<std::vector<double>>("EBchi2Parameters", {2.122,0.022,2.122,0.022}, true) and
537  edm::ParameterDescription<double>("EEtimeConstantTerm", 1.0, true) and
538  edm::ParameterDescription<double>("chi2ThreshEE_", 95.0, true) and
539  edm::ParameterDescription<double>("EBtimeFitLimits_Upper", 1.4, true));
540 
541  return psd;
542 }
edm::ParameterSetDescription::addNode
ParameterDescriptionNode * addNode(ParameterDescriptionNode const &node)
Definition: ParameterSetDescription.cc:41
template<class C >
int EcalUncalibRecHitWorkerGlobal::isSaturated ( const C &  digi)
protected

Definition at line 135 of file EcalUncalibRecHitWorkerGlobal.cc.

References mps_fire::i.

136 {
137  //bool saturated_ = 0;
138  int cnt;
139  for (int j = 0; j < C::MAXSAMPLES - 5; ++j) {
140  cnt = 0;
141  for (int i = j; i < (j + 5) && i < C::MAXSAMPLES; ++i) {
142  if ( dataFrame.sample(i).gainId() == 0 ) ++cnt;
143  }
144  if ( cnt == 5 ) return j-1 ; // the last unsaturated sample
145  }
146  return -1; // no saturation found
147 }
mps_fire.i
i
Definition: mps_fire.py:330
bool EcalUncalibRecHitWorkerGlobal::run ( const edm::Event &  evt,
const EcalDigiCollection::const_iterator &  digi,
EcalUncalibratedRecHitCollection &  result 
)
overridevirtual

Implements EcalUncalibRecHitWorkerRunOneDigiBase.

Definition at line 215 of file EcalUncalibRecHitWorkerGlobal.cc.

References CustomPhysics_cfi::amplitude, EcalUncalibratedRecHit::amplitude(), amplitudeThreshEB_, amplitudeThreshEE_, vertices_cff::chi2, EcalUncalibRecHitRecChi2Algo< C >::chi2(), chi2ThreshEB_, chi2ThreshEE_, EcalUncalibRecHitRatioMethodAlgo< C >::computeAmplitude(), EcalUncalibRecHitRatioMethodAlgo< C >::computeTime(), EBamplitudeFitParameters_, EBchi2Parameters_, EBtimeConstantTerm_, EcalTimeBiasCorrections::EBTimeCorrAmplitudeBins, EcalTimeBiasCorrections::EBTimeCorrShiftBins, EBtimeFitLimits_, EBtimeFitParameters_, EBtimeNconst_, EcalEndcap, EEamplitudeFitParameters_, EEchi2Parameters_, EEtimeConstantTerm_, EcalTimeBiasCorrections::EETimeCorrAmplitudeBins, EcalTimeBiasCorrections::EETimeCorrShiftBins, EEtimeFitLimits_, EEtimeFitParameters_, EEtimeNconst_, EcalUncalibRecHitRatioMethodAlgo< C >::fixMGPAslew(), EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), ecalMGPA::gainId(), gainRatios, gains, EcalUncalibRecHitRatioMethodAlgo< C >::getCalculatedRecHit(), EcalTimeOffsetConstant::getEBValue(), EcalTimeOffsetConstant::getEEValue(), EcalTBWeights::getMap(), EcalWeightSet::getWeightsAfterGainSwitch(), EcalWeightSet::getWeightsBeforeGainSwitch(), grps, ecalpyutils::hashedIndex(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), EcalXtalGroupId::id(), EcalUncalibRecHitRatioMethodAlgo< C >::init(), itime, EcalUncalibratedRecHit::jitter(), EcalUncalibratedRecHit::kHasSwitchToGain1, EcalUncalibratedRecHit::kHasSwitchToGain6, EcalUncalibratedRecHit::kOutOfTime, EcalUncalibratedRecHit::kPoorReco, kPoorRecoFlagEB_, kPoorRecoFlagEE_, EcalUncalibratedRecHit::kSaturated, EcalUncalibRecHitRecWeightsAlgo< C >::makeRecHit(), EcalDataFrame::MAXSAMPLES, EcalPedestal::mean_x1, EcalPedestal::mean_x12, EcalPedestal::mean_x6, offtime, outOfTimeThreshG12mEB_, outOfTimeThreshG12mEE_, outOfTimeThreshG12pEB_, outOfTimeThreshG12pEE_, outOfTimeThreshG61mEB_, outOfTimeThreshG61mEE_, outOfTimeThreshG61pEB_, outOfTimeThreshG61pEE_, muonCSCDigis_cfi::pedestal, pedRMSVec, peds, pedVec, funct::pow(), edm::ESHandle< T >::product(), edm::SortedCollection< T, SORT >::push_back(), ratioMethod_barrel_, ratioMethod_endcap_, EcalPedestal::rms_x1, EcalPedestal::rms_x12, EcalPedestal::rms_x6, simplePhotonAnalyzer_cfi::sample, sampleMaskHand_, EcalUncalibratedRecHit::setChi2(), EcalUncalibratedRecHit::setFlagBit(), EcalUncalibratedRecHit::setJitter(), EcalUncalibratedRecHit::setJitterError(), mathSSE::sqrt(), testbeamEBShape, testbeamEEShape, timeCorrBias_, timeCorrection(), weights, weightsMethod_barrel_, weightsMethod_endcap_, and wgts.

218 {
219  DetId detid(itdg->id());
220 
221  const EcalSampleMask *sampleMask_ = sampleMaskHand_.product();
222 
223  // intelligence for recHit computation
224  EcalUncalibratedRecHit uncalibRecHit;
225 
226 
227  const EcalPedestals::Item * aped = nullptr;
228  const EcalMGPAGainRatio * aGain = nullptr;
229  const EcalXtalGroupId * gid = nullptr;
230  float offsetTime = 0;
231 
232  if (detid.subdetId()==EcalEndcap) {
233  unsigned int hashedIndex = EEDetId(detid).hashedIndex();
234  aped = &peds->endcap(hashedIndex);
235  aGain = &gains->endcap(hashedIndex);
236  gid = &grps->endcap(hashedIndex);
237  offsetTime = offtime->getEEValue();
238  } else {
239  unsigned int hashedIndex = EBDetId(detid).hashedIndex();
240  aped = &peds->barrel(hashedIndex);
241  aGain = &gains->barrel(hashedIndex);
242  gid = &grps->barrel(hashedIndex);
243  offsetTime = offtime->getEBValue();
244  }
245 
246  pedVec[0] = aped->mean_x12;
247  pedVec[1] = aped->mean_x6;
248  pedVec[2] = aped->mean_x1;
249  pedRMSVec[0] = aped->rms_x12;
250  pedRMSVec[1] = aped->rms_x6;
251  pedRMSVec[2] = aped->rms_x1;
252  gainRatios[0] = 1.;
253  gainRatios[1] = aGain->gain12Over6();
254  gainRatios[2] = aGain->gain6Over1()*aGain->gain12Over6();
255 
256  // compute the right bin of the pulse shape using time calibration constants
257  EcalTimeCalibConstantMap::const_iterator it = itime->find( detid );
258  EcalTimeCalibConstant itimeconst = 0;
259  if( it != itime->end() ) {
260  itimeconst = (*it);
261  } else {
262  edm::LogError("EcalRecHitError") << "No time intercalib const found for xtal "
263  << detid.rawId()
264  << "! something wrong with EcalTimeCalibConstants in your DB? ";
265  }
266 
267 
268  // === amplitude computation ===
269  int leadingSample = -1;
270  if (detid.subdetId()==EcalEndcap) {
271  leadingSample = ((EcalDataFrame)(*itdg)).lastUnsaturatedSample();
272  } else {
273  leadingSample = ((EcalDataFrame)(*itdg)).lastUnsaturatedSample();
274  }
275 
276  if ( leadingSample == 4 ) { // saturation on the expected max sample
277  uncalibRecHit = EcalUncalibratedRecHit( (*itdg).id(), 4095*12, 0, 0, 0);
278  uncalibRecHit.setFlagBit( EcalUncalibratedRecHit::kSaturated );
279  // do not propagate the default chi2 = -1 value to the calib rechit (mapped to 64), set it to 0 when saturation
280  uncalibRecHit.setChi2(0);
281  } else if ( leadingSample >= 0 ) { // saturation on other samples: cannot extrapolate from the fourth one
282  double pedestal = 0.;
283  double gainratio = 1.;
284  int gainId = ((EcalDataFrame)(*itdg)).sample(5).gainId();
285 
286  if (gainId==0 || gainId==3) {
287  pedestal = aped->mean_x1;
288  gainratio = aGain->gain6Over1()*aGain->gain12Over6();
289  }
290  else if (gainId==1) {
291  pedestal = aped->mean_x12;
292  gainratio = 1.;
293  }
294  else if (gainId==2) {
295  pedestal = aped->mean_x6;
296  gainratio = aGain->gain12Over6();
297  }
298  double amplitude = ((double)(((EcalDataFrame)(*itdg)).sample(5).adc()) - pedestal) * gainratio;
299  uncalibRecHit = EcalUncalibratedRecHit( (*itdg).id(), amplitude, 0, 0, 0);
300  uncalibRecHit.setFlagBit( EcalUncalibratedRecHit::kSaturated );
301  // do not propagate the default chi2 = -1 value to the calib rechit (mapped to 64), set it to 0 when saturation
302  uncalibRecHit.setChi2(0);
303  } else {
304  // weights method
305  EcalTBWeights::EcalTDCId tdcid(1);
306  EcalTBWeights::EcalTBWeightMap const & wgtsMap = wgts->getMap();
307  EcalTBWeights::EcalTBWeightMap::const_iterator wit;
308  wit = wgtsMap.find( std::make_pair(*gid,tdcid) );
309  if( wit == wgtsMap.end() ) {
310  edm::LogError("EcalUncalibRecHitError") << "No weights found for EcalGroupId: "
311  << gid->id() << " and EcalTDCId: " << tdcid
312  << "\n skipping digi with id: " << detid.rawId();
313 
314  return false;
315  }
316  const EcalWeightSet& wset = wit->second; // this is the EcalWeightSet
317 
318  const EcalWeightSet::EcalWeightMatrix& mat1 = wset.getWeightsBeforeGainSwitch();
319  const EcalWeightSet::EcalWeightMatrix& mat2 = wset.getWeightsAfterGainSwitch();
320 
321  weights[0] = &mat1;
322  weights[1] = &mat2;
323 
324  // get uncalibrated recHit from weights
325  if (detid.subdetId()==EcalEndcap) {
326  uncalibRecHit = weightsMethod_endcap_.makeRecHit(*itdg, pedVec, pedRMSVec, gainRatios, weights, testbeamEEShape);
327  } else {
328  uncalibRecHit = weightsMethod_barrel_.makeRecHit(*itdg, pedVec, pedRMSVec, gainRatios, weights, testbeamEBShape);
329  }
330 
331  // === time computation ===
332  // ratio method
333  float const clockToNsConstant = 25.;
334  if (detid.subdetId()==EcalEndcap) {
335  ratioMethod_endcap_.init( *itdg, *sampleMask_, pedVec, pedRMSVec, gainRatios );
336  ratioMethod_endcap_.computeTime( EEtimeFitParameters_, EEtimeFitLimits_, EEamplitudeFitParameters_ );
337  ratioMethod_endcap_.computeAmplitude( EEamplitudeFitParameters_);
338  EcalUncalibRecHitRatioMethodAlgo<EEDataFrame>::CalculatedRecHit crh = ratioMethod_endcap_.getCalculatedRecHit();
339  double theTimeCorrectionEE = timeCorrection(uncalibRecHit.amplitude(),
340  timeCorrBias_->EETimeCorrAmplitudeBins, timeCorrBias_->EETimeCorrShiftBins);
341 
342  uncalibRecHit.setJitter( crh.timeMax - 5 + theTimeCorrectionEE);
343  uncalibRecHit.setJitterError( std::sqrt(pow(crh.timeError,2) + std::pow(EEtimeConstantTerm_,2)/std::pow(clockToNsConstant,2)) );
344  // consider flagging as kOutOfTime only if above noise
345  if (uncalibRecHit.amplitude() > pedRMSVec[0] * amplitudeThreshEE_){
346  float outOfTimeThreshP = outOfTimeThreshG12pEE_;
347  float outOfTimeThreshM = outOfTimeThreshG12mEE_;
348  // determine if gain has switched away from gainId==1 (x12 gain)
349  // and determine cuts (number of 'sigmas') to ose for kOutOfTime
350  // >3k ADC is necessasry condition for gain switch to occur
351  if (uncalibRecHit.amplitude() > 3000.){
352  for (int iSample = 0; iSample < EEDataFrame::MAXSAMPLES; iSample++) {
353  int GainId = ((EcalDataFrame)(*itdg)).sample(iSample).gainId();
354  if (GainId!=1) {
355  outOfTimeThreshP = outOfTimeThreshG61pEE_;
356  outOfTimeThreshM = outOfTimeThreshG61mEE_;
357  break;
358  }
359  }}
360  float correctedTime = (crh.timeMax-5) * clockToNsConstant + itimeconst + offsetTime;
361  float cterm = EEtimeConstantTerm_;
362  float sigmaped = pedRMSVec[0]; // approx for lower gains
363  float nterm = EEtimeNconst_*sigmaped/uncalibRecHit.amplitude();
364  float sigmat = std::sqrt( nterm*nterm + cterm*cterm );
365  if ( ( correctedTime > sigmat*outOfTimeThreshP ) ||
366  ( correctedTime < (-1.*sigmat*outOfTimeThreshM) ))
367  { uncalibRecHit.setFlagBit( EcalUncalibratedRecHit::kOutOfTime ); }
368  }
369 
370  } else {
371  ratioMethod_barrel_.init( *itdg, *sampleMask_, pedVec, pedRMSVec, gainRatios );
372  ratioMethod_barrel_.fixMGPAslew(*itdg);
373  ratioMethod_barrel_.computeTime( EBtimeFitParameters_, EBtimeFitLimits_, EBamplitudeFitParameters_ );
374  ratioMethod_barrel_.computeAmplitude( EBamplitudeFitParameters_);
375  EcalUncalibRecHitRatioMethodAlgo<EBDataFrame>::CalculatedRecHit crh = ratioMethod_barrel_.getCalculatedRecHit();
376 
377  double theTimeCorrectionEB = timeCorrection(uncalibRecHit.amplitude(),
378  timeCorrBias_->EBTimeCorrAmplitudeBins, timeCorrBias_->EBTimeCorrShiftBins);
379 
380  uncalibRecHit.setJitter( crh.timeMax - 5 + theTimeCorrectionEB);
381 
382  uncalibRecHit.setJitterError( std::sqrt(std::pow(crh.timeError,2) + std::pow(EBtimeConstantTerm_,2)/std::pow(clockToNsConstant,2)) );
383  // consider flagging as kOutOfTime only if above noise
384  if (uncalibRecHit.amplitude() > pedRMSVec[0] * amplitudeThreshEB_){
385  float outOfTimeThreshP = outOfTimeThreshG12pEB_;
386  float outOfTimeThreshM = outOfTimeThreshG12mEB_;
387  // determine if gain has switched away from gainId==1 (x12 gain)
388  // and determine cuts (number of 'sigmas') to ose for kOutOfTime
389  // >3k ADC is necessasry condition for gain switch to occur
390  if (uncalibRecHit.amplitude() > 3000.){
391  for (int iSample = 0; iSample < EBDataFrame::MAXSAMPLES; iSample++) {
392  int GainId = ((EcalDataFrame)(*itdg)).sample(iSample).gainId();
393  if (GainId!=1) {
394  outOfTimeThreshP = outOfTimeThreshG61pEB_;
395  outOfTimeThreshM = outOfTimeThreshG61mEB_;
396  break;}
397  } }
398  float correctedTime = (crh.timeMax-5) * clockToNsConstant + itimeconst + offsetTime;
399  float cterm = EBtimeConstantTerm_;
400  float sigmaped = pedRMSVec[0]; // approx for lower gains
401  float nterm = EBtimeNconst_*sigmaped/uncalibRecHit.amplitude();
402  float sigmat = std::sqrt( nterm*nterm + cterm*cterm );
403  if ( ( correctedTime > sigmat*outOfTimeThreshP ) ||
404  ( correctedTime < (-1.*sigmat*outOfTimeThreshM) ))
405  { uncalibRecHit.setFlagBit( EcalUncalibratedRecHit::kOutOfTime ); }
406  }
407  }
408 
409  // === chi2express ===
410  if (detid.subdetId()==EcalEndcap) {
411 
412  double amplitude = uncalibRecHit.amplitude();
413  double amplitudeOutOfTime = 0.;
414  double jitter= uncalibRecHit.jitter();
415 
416 
417 
418  EcalUncalibRecHitRecChi2Algo<EEDataFrame>chi2expressEE_(
419  *itdg,
420  amplitude,
421  (itimeconst + offsetTime),
422  amplitudeOutOfTime,
423  jitter,
424  pedVec,
425  pedRMSVec,
426  gainRatios,
427  testbeamEEShape,
428  EEchi2Parameters_
429  );
430  double chi2 = chi2expressEE_.chi2();
431  uncalibRecHit.setChi2(chi2);
432 
433  if(kPoorRecoFlagEE_)
434  {
435 
436  if (chi2>chi2ThreshEE_) {
437 
438  // first check if all samples are ok, if not don't use chi2 to flag
439  bool samplesok = true;
440  for (int sample =0; sample < EcalDataFrame::MAXSAMPLES; ++sample) {
441  if (!sampleMask_->useSampleEE(sample)) {
442  samplesok = false;
443  break;
444  }
445  }
446  if (samplesok) uncalibRecHit.setFlagBit(EcalUncalibratedRecHit::kPoorReco);
447  }
448 
449 
450  }
451 
452  } else {
453  double amplitude = uncalibRecHit.amplitude();
454  double amplitudeOutOfTime = 0.;
455  double jitter= uncalibRecHit.jitter();
456 
457  EcalUncalibRecHitRecChi2Algo<EBDataFrame>chi2expressEB_(
458  *itdg,
459  amplitude,
460  (itimeconst + offsetTime),
461  amplitudeOutOfTime,
462  jitter,
463  pedVec,
464  pedRMSVec,
465  gainRatios,
466  testbeamEBShape,
467  EBchi2Parameters_
468  );
469  double chi2 = chi2expressEB_.chi2();
470  uncalibRecHit.setChi2(chi2);
471 
472  if(kPoorRecoFlagEB_)
473  {
474 
475  if(chi2>chi2ThreshEB_){
476  // first check if all samples are ok, if not don't use chi2 to flag
477  bool samplesok = true;
478  for (int sample =0; sample < EcalDataFrame::MAXSAMPLES; ++sample) {
479  if (!sampleMask_->useSampleEB(sample)) {
480  samplesok = false;
481  break;
482  }
483  }
484  if (samplesok) uncalibRecHit.setFlagBit(EcalUncalibratedRecHit::kPoorReco);
485  }
486 
487  }
488  }
489  }
490 
491  // set flags if gain switch has occurred
492  if( ((EcalDataFrame)(*itdg)).hasSwitchToGain6() ) uncalibRecHit.setFlagBit( EcalUncalibratedRecHit::kHasSwitchToGain6 );
493  if( ((EcalDataFrame)(*itdg)).hasSwitchToGain1() ) uncalibRecHit.setFlagBit( EcalUncalibratedRecHit::kHasSwitchToGain1 );
494 
495  // put the recHit in the collection
496  if (detid.subdetId()==EcalEndcap) {
497  result.push_back( uncalibRecHit );
498  } else {
499  result.push_back( uncalibRecHit );
500  }
501 
502  return true;
503 }
EcalUncalibRecHitWorkerGlobal::offtime
edm::ESHandle< EcalTimeOffsetConstant > offtime
Definition: EcalUncalibRecHitWorkerGlobal.h:105
EcalTBWeights::EcalTBWeightMap
std::map< std::pair< EcalXtalGroupId, EcalTDCId >, EcalWeightSet > EcalTBWeightMap
Definition: EcalTBWeights.h:21
EcalUncalibRecHitWorkerGlobal::timeCorrBias_
edm::ESHandle< EcalTimeBiasCorrections > timeCorrBias_
Definition: EcalUncalibRecHitWorkerGlobal.h:102
EBDetId::hashedIndex
int hashedIndex() const
get a compact index for arrays
Definition: EBDetId.h:82
EcalUncalibRecHitRatioMethodAlgo::computeAmplitude
void computeAmplitude(std::vector< double > &amplitudeFitParameters)
Definition: EcalUncalibRecHitRatioMethodAlgo.h:560
EcalTimeBiasCorrections::EBTimeCorrShiftBins
std::vector< float > EBTimeCorrShiftBins
Definition: EcalTimeBiasCorrections.h:24
EcalUncalibRecHitWorkerGlobal::timeCorrection
double timeCorrection(float ampli, const std::vector< float > &amplitudeBins, const std::vector< float > &shiftBins)
Definition: EcalUncalibRecHitWorkerGlobal.cc:158
EcalUncalibratedRecHit::setJitter
void setJitter(float jitter)
Definition: EcalUncalibratedRecHit.h:46
simplePhotonAnalyzer_cfi.sample
sample
Definition: simplePhotonAnalyzer_cfi.py:12
EcalUncalibRecHitWorkerGlobal::EBtimeFitLimits_
std::pair< double, double > EBtimeFitLimits_
Definition: EcalUncalibRecHitWorkerGlobal.h:80
EcalUncalibRecHitRatioMethodAlgo::computeTime
void computeTime(std::vector< double > &timeFitParameters, std::pair< double, double > &timeFitLimits, std::vector< double > &amplitudeFitParameters)
Definition: EcalUncalibRecHitRatioMethodAlgo.h:224
edm::SortedCollection::push_back
void push_back(T const &t)
Definition: SortedCollection.h:199
EcalUncalibRecHitWorkerGlobal::weightsMethod_barrel_
EcalUncalibRecHitRecWeightsAlgo< EBDataFrame > weightsMethod_barrel_
Definition: EcalUncalibRecHitWorkerGlobal.h:67
EcalUncalibratedRecHit::setJitterError
void setJitterError(float jitterErr)
Definition: EcalUncalibratedRecHit.cc:41
CustomPhysics_cfi.amplitude
amplitude
Definition: CustomPhysics_cfi.py:12
EcalUncalibRecHitWorkerGlobal::EEtimeFitLimits_
std::pair< double, double > EEtimeFitLimits_
Definition: EcalUncalibRecHitWorkerGlobal.h:81
EcalUncalibRecHitWorkerGlobal::ratioMethod_endcap_
EcalUncalibRecHitRatioMethodAlgo< EEDataFrame > ratioMethod_endcap_
Definition: EcalUncalibRecHitWorkerGlobal.h:84
EcalXtalGroupId::id
const unsigned int id() const
Definition: EcalXtalGroupId.h:22
EcalTimeBiasCorrections::EBTimeCorrAmplitudeBins
std::vector< float > EBTimeCorrAmplitudeBins
Definition: EcalTimeBiasCorrections.h:23
EcalTimeBiasCorrections::EETimeCorrShiftBins
std::vector< float > EETimeCorrShiftBins
Definition: EcalTimeBiasCorrections.h:27
ecalpyutils::hashedIndex
int hashedIndex(int ieta, int iphi)
Definition: EcalPyUtils.cc:42
EcalUncalibRecHitWorkerGlobal::ratioMethod_barrel_
EcalUncalibRecHitRatioMethodAlgo< EBDataFrame > ratioMethod_barrel_
Definition: EcalUncalibRecHitWorkerGlobal.h:83
EcalUncalibRecHitWorkerGlobal::gains
edm::ESHandle< EcalGainRatios > gains
Definition: EcalUncalibRecHitWorkerGlobal.h:55
EcalTimeBiasCorrections::EETimeCorrAmplitudeBins
std::vector< float > EETimeCorrAmplitudeBins
Definition: EcalTimeBiasCorrections.h:26
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
EcalUncalibRecHitWorkerGlobal::peds
edm::ESHandle< EcalPedestals > peds
Definition: EcalUncalibRecHitWorkerGlobal.h:54
EcalUncalibRecHitRecWeightsAlgo::makeRecHit
virtual EcalUncalibratedRecHit makeRecHit(const C &dataFrame, const double *pedestals, const double *pedestalsRMS, const double *gainRatios, const EcalWeightSet::EcalWeightMatrix **weights, const EcalShapeBase &testbeamPulseShape)
Compute parameters.
Definition: EcalUncalibRecHitRecWeightsAlgo.h:30
EcalUncalibRecHitWorkerGlobal::wgts
edm::ESHandle< EcalTBWeights > wgts
Definition: EcalUncalibRecHitWorkerGlobal.h:64
EcalMGPAGainRatio::gain6Over1
float gain6Over1() const
Definition: EcalMGPAGainRatio.h:21
EcalUncalibRecHitWorkerGlobal::grps
edm::ESHandle< EcalWeightXtalGroups > grps
Definition: EcalUncalibRecHitWorkerGlobal.h:63
EcalWeightSet::getWeightsAfterGainSwitch
EcalWeightMatrix & getWeightsAfterGainSwitch()
Definition: EcalWeightSet.h:30
DetId
Definition: DetId.h:18
EEDetId::hashedIndex
int hashedIndex() const
Definition: EEDetId.h:183
EcalWeightSet::getWeightsBeforeGainSwitch
EcalWeightMatrix & getWeightsBeforeGainSwitch()
Definition: EcalWeightSet.h:29
EcalCondObjectContainer< float >::const_iterator
std::vector< Item >::const_iterator const_iterator
Definition: EcalCondObjectContainer.h:17
EcalMGPAGainRatio::gain12Over6
float gain12Over6() const
Definition: EcalMGPAGainRatio.h:20
muonCSCDigis_cfi.pedestal
pedestal
Definition: muonCSCDigis_cfi.py:12
EcalUncalibRecHitWorkerGlobal::weights
const EcalWeightSet::EcalWeightMatrix * weights[2]
Definition: EcalUncalibRecHitWorkerGlobal.h:65
ecalMGPA::gainId
constexpr int gainId(sample_type sample)
get the gainId (2 bits)
Definition: EcalMGPASample.h:13
EcalTimeCalibConstant
float EcalTimeCalibConstant
Definition: EcalTimeCalibConstants.h:10
EcalUncalibRecHitWorkerGlobal::itime
edm::ESHandle< EcalTimeCalibConstants > itime
Definition: EcalUncalibRecHitWorkerGlobal.h:104
EcalUncalibRecHitWorkerGlobal::weightsMethod_endcap_
EcalUncalibRecHitRecWeightsAlgo< EEDataFrame > weightsMethod_endcap_
Definition: EcalUncalibRecHitWorkerGlobal.h:68
EcalTBWeights::getMap
const EcalTBWeightMap & getMap() const
Definition: EcalTBWeights.h:31
EcalWeightSet::EcalWeightMatrix
math::Matrix< 3, 10 >::type EcalWeightMatrix
Definition: EcalWeightSet.h:22
EcalDataFrame::MAXSAMPLES
static constexpr int MAXSAMPLES
Definition: EcalDataFrame.h:48
EcalUncalibRecHitWorkerGlobal::sampleMaskHand_
edm::ESHandle< EcalSampleMask > sampleMaskHand_
Definition: EcalUncalibRecHitWorkerGlobal.h:73
EcalUncalibRecHitRatioMethodAlgo::init
void init(const C &dataFrame, const EcalSampleMask &sampleMask, const double *pedestals, const double *pedestalRMSes, const double *gainRatios)
Definition: EcalUncalibRecHitRatioMethodAlgo.h:95
vertices_cff.chi2
chi2
Definition: vertices_cff.py:33
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
void EcalUncalibRecHitWorkerGlobal::set ( const edm::EventSetup &  es)
overridevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 106 of file EcalUncalibRecHitWorkerGlobal.cc.

References gains, edm::EventSetup::get(), grps, itime, offtime, peds, sampleMaskHand_, EcalShapeBase::setEventSetup(), testbeamEBShape, testbeamEEShape, timeCorrBias_, and wgts.

107 {
108  // common setup
109  es.get<EcalGainRatiosRcd>().get(gains);
110  es.get<EcalPedestalsRcd>().get(peds);
111 
112  // for the weights method
113  es.get<EcalWeightXtalGroupsRcd>().get(grps);
114  es.get<EcalTBWeightsRcd>().get(wgts);
115 
116  // which of the samples need be used
117  es.get<EcalSampleMaskRcd>().get(sampleMaskHand_);
118 
119  // for the ratio method
120 
121  es.get<EcalTimeCalibConstantsRcd>().get(itime);
122  es.get<EcalTimeOffsetConstantRcd>().get(offtime);
123 
124  // for the time correction methods
125  es.get<EcalTimeBiasCorrectionsRcd>().get(timeCorrBias_);
126 
127  // for the DB Ecal Pulse Sim Shape
128  testbeamEEShape.setEventSetup(es);
129  testbeamEBShape.setEventSetup(es);
130 }
EcalUncalibRecHitWorkerGlobal::offtime
edm::ESHandle< EcalTimeOffsetConstant > offtime
Definition: EcalUncalibRecHitWorkerGlobal.h:105
EcalUncalibRecHitWorkerGlobal::timeCorrBias_
edm::ESHandle< EcalTimeBiasCorrections > timeCorrBias_
Definition: EcalUncalibRecHitWorkerGlobal.h:102
EcalUncalibRecHitWorkerGlobal::gains
edm::ESHandle< EcalGainRatios > gains
Definition: EcalUncalibRecHitWorkerGlobal.h:55
EcalUncalibRecHitWorkerGlobal::peds
edm::ESHandle< EcalPedestals > peds
Definition: EcalUncalibRecHitWorkerGlobal.h:54
EcalUncalibRecHitWorkerGlobal::wgts
edm::ESHandle< EcalTBWeights > wgts
Definition: EcalUncalibRecHitWorkerGlobal.h:64
EcalUncalibRecHitWorkerGlobal::grps
edm::ESHandle< EcalWeightXtalGroups > grps
Definition: EcalUncalibRecHitWorkerGlobal.h:63
edm::EventSetup::get
T get() const
Definition: EventSetup.h:71
EcalUncalibRecHitWorkerGlobal::itime
edm::ESHandle< EcalTimeCalibConstants > itime
Definition: EcalUncalibRecHitWorkerGlobal.h:104
EcalShapeBase::setEventSetup
void setEventSetup(const edm::EventSetup &evtSetup)
Definition: EcalShapeBase.cc:19
EcalUncalibRecHitWorkerGlobal::sampleMaskHand_
edm::ESHandle< EcalSampleMask > sampleMaskHand_
Definition: EcalUncalibRecHitWorkerGlobal.h:73
double EcalUncalibRecHitWorkerGlobal::timeCorrection ( float  ampli,
const std::vector< float > &  amplitudeBins,
const std::vector< float > &  shiftBins 
)
protected

Amplitude-dependent time corrections; EE and EB have separate corrections: EXtimeCorrAmplitudes (ADC) and EXtimeCorrShifts (ns) need to have the same number of elements Bins must be ordered in amplitude. First-last bins take care of under-overflows.

The algorithm is the same for EE and EB, only the correction vectors are different.

Returns
Jitter (in clock cycles) which will be added to UncalibRechit.setJitter(), 0 if no correction is applied.

Definition at line 158 of file EcalUncalibRecHitWorkerGlobal.cc.

References stringResolutionProvider_cfi::bin, and createfilelist::int.

Referenced by run().

161  {
162 
163  // computed initially in ns. Than turned in the BX's, as
164  // EcalUncalibratedRecHit need be.
165  double theCorrection = 0;
166 
167  // sanity check for arrays
168  if (amplitudeBins.empty()) {
169  edm::LogError("EcalRecHitError")
170  << "timeCorrAmplitudeBins is empty, forcing no time bias corrections.";
171 
172  return 0;
173  }
174 
175  if (amplitudeBins.size() != shiftBins.size()) {
176  edm::LogError("EcalRecHitError")
177  << "Size of timeCorrAmplitudeBins different from "
178  "timeCorrShiftBins. Forcing no time bias corrections. ";
179 
180  return 0;
181  }
182 
183  int myBin = -1;
184  for (int bin = 0; bin < (int) amplitudeBins.size(); bin++) {
185  if (ampli > amplitudeBins.at(bin)) {
186  myBin = bin;
187  } else {
188  break;
189  }
190  }
191 
192  if (myBin == -1) {
193  theCorrection = shiftBins.at(0);
194  } else if (myBin == ((int)(amplitudeBins.size() - 1))) {
195  theCorrection = shiftBins.at(myBin);
196  } else if (-1 < myBin && myBin < ((int) amplitudeBins.size() - 1)) {
197  // interpolate linearly between two assingned points
198  theCorrection = (shiftBins.at(myBin + 1) - shiftBins.at(myBin));
199  theCorrection *= (((double) ampli) - amplitudeBins.at(myBin)) /
200  (amplitudeBins.at(myBin + 1) - amplitudeBins.at(myBin));
201  theCorrection += shiftBins.at(myBin);
202  } else {
203  edm::LogError("EcalRecHitError")
204  << "Assigning time correction impossible. Setting it to 0 ";
205  theCorrection = 0.;
206  }
207 
208  // convert ns into clocks
209  return theCorrection / 25.;
210 }
stringResolutionProvider_cfi.bin
bin
set the eta bin as selection string.
Definition: stringResolutionProvider_cfi.py:11

Member Data Documentation

double EcalUncalibRecHitWorkerGlobal::amplitudeThreshEB_
protected

Definition at line 98 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::amplitudeThreshEE_
protected

Definition at line 99 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

const EcalWeightSet::EcalChi2WeightMatrix* EcalUncalibRecHitWorkerGlobal::chi2mat[2]
protected

Definition at line 66 of file EcalUncalibRecHitWorkerGlobal.h.

double EcalUncalibRecHitWorkerGlobal::chi2ThreshEB_
protected

Definition at line 112 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::chi2ThreshEE_
protected

Definition at line 113 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::vector<double> EcalUncalibRecHitWorkerGlobal::EBamplitudeFitParameters_
protected

Definition at line 78 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::vector<double> EcalUncalibRecHitWorkerGlobal::EBchi2Parameters_
protected

Definition at line 114 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::vector<double> EcalUncalibRecHitWorkerGlobal::ebPulseShape_
protected

Definition at line 106 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal().

double EcalUncalibRecHitWorkerGlobal::ebSpikeThresh_
protected

Definition at line 100 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal().

double EcalUncalibRecHitWorkerGlobal::EBtimeConstantTerm_
protected

Definition at line 86 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::pair<double,double> EcalUncalibRecHitWorkerGlobal::EBtimeFitLimits_
protected

Definition at line 80 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::vector<double> EcalUncalibRecHitWorkerGlobal::EBtimeFitParameters_
protected

Definition at line 76 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::EBtimeNconst_
protected

Definition at line 87 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::vector<double> EcalUncalibRecHitWorkerGlobal::EEamplitudeFitParameters_
protected

Definition at line 79 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::vector<double> EcalUncalibRecHitWorkerGlobal::EEchi2Parameters_
protected

Definition at line 115 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::vector<double> EcalUncalibRecHitWorkerGlobal::eePulseShape_
protected

Definition at line 107 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal().

double EcalUncalibRecHitWorkerGlobal::EEtimeConstantTerm_
protected

Definition at line 88 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::pair<double,double> EcalUncalibRecHitWorkerGlobal::EEtimeFitLimits_
protected

Definition at line 81 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

std::vector<double> EcalUncalibRecHitWorkerGlobal::EEtimeFitParameters_
protected

Definition at line 77 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::EEtimeNconst_
protected

Definition at line 89 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::gainRatios[3]
protected

Definition at line 52 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run().

edm::ESHandle<EcalGainRatios> EcalUncalibRecHitWorkerGlobal::gains
protected

Definition at line 55 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

edm::ESHandle<EcalWeightXtalGroups> EcalUncalibRecHitWorkerGlobal::grps
protected

Definition at line 63 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

edm::ESHandle<EcalTimeCalibConstants> EcalUncalibRecHitWorkerGlobal::itime
protected

Definition at line 104 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

bool EcalUncalibRecHitWorkerGlobal::kPoorRecoFlagEB_
protected

Definition at line 110 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

bool EcalUncalibRecHitWorkerGlobal::kPoorRecoFlagEE_
protected

Definition at line 111 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

edm::ESHandle<EcalTimeOffsetConstant> EcalUncalibRecHitWorkerGlobal::offtime
protected

Definition at line 105 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

double EcalUncalibRecHitWorkerGlobal::outOfTimeThreshG12mEB_
protected

Definition at line 91 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::outOfTimeThreshG12mEE_
protected

Definition at line 95 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::outOfTimeThreshG12pEB_
protected

Definition at line 90 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::outOfTimeThreshG12pEE_
protected

Definition at line 94 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::outOfTimeThreshG61mEB_
protected

Definition at line 93 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::outOfTimeThreshG61mEE_
protected

Definition at line 97 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::outOfTimeThreshG61pEB_
protected

Definition at line 92 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::outOfTimeThreshG61pEE_
protected

Definition at line 96 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by EcalUncalibRecHitWorkerGlobal(), and run().

double EcalUncalibRecHitWorkerGlobal::pedRMSVec[3]
protected

Definition at line 51 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run().

edm::ESHandle<EcalPedestals> EcalUncalibRecHitWorkerGlobal::peds
protected

Definition at line 54 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

double EcalUncalibRecHitWorkerGlobal::pedVec[3]
protected

Definition at line 50 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run().

EcalUncalibRecHitRatioMethodAlgo<EBDataFrame> EcalUncalibRecHitWorkerGlobal::ratioMethod_barrel_
protected

Definition at line 83 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run().

EcalUncalibRecHitRatioMethodAlgo<EEDataFrame> EcalUncalibRecHitWorkerGlobal::ratioMethod_endcap_
protected

Definition at line 84 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run().

edm::ESHandle<EcalSampleMask> EcalUncalibRecHitWorkerGlobal::sampleMaskHand_
protected

Definition at line 73 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

EBShape EcalUncalibRecHitWorkerGlobal::testbeamEBShape
protected

Definition at line 70 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

EEShape EcalUncalibRecHitWorkerGlobal::testbeamEEShape
protected

Definition at line 69 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

edm::ESHandle<EcalTimeBiasCorrections> EcalUncalibRecHitWorkerGlobal::timeCorrBias_
protected

Definition at line 102 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

const EcalWeightSet::EcalWeightMatrix* EcalUncalibRecHitWorkerGlobal::weights[2]
protected

Definition at line 65 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run().

EcalUncalibRecHitRecWeightsAlgo<EBDataFrame> EcalUncalibRecHitWorkerGlobal::weightsMethod_barrel_
protected

Definition at line 67 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run().

EcalUncalibRecHitRecWeightsAlgo<EEDataFrame> EcalUncalibRecHitWorkerGlobal::weightsMethod_endcap_
protected

Definition at line 68 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run().

edm::ESHandle<EcalTBWeights> EcalUncalibRecHitWorkerGlobal::wgts
protected

Definition at line 64 of file EcalUncalibRecHitWorkerGlobal.h.

Referenced by run(), and set().

Generated for CMSSW Reference Manual by doxygen 1.8.11