CMS 3D CMS Logo
List of all members | Public Member Functions | Private Types | Private Member Functions | Private Attributes
EcalUncalibRecHitWorkerMultiFit Class Referencefinal
Inheritance diagram for EcalUncalibRecHitWorkerMultiFit:
EcalUncalibRecHitWorkerBaseClass

Public Member Functions

 EcalUncalibRecHitWorkerMultiFit (const edm::ParameterSet &, edm::ConsumesCollector &c)
 
 EcalUncalibRecHitWorkerMultiFit ()
 
edm::ParameterSetDescription getAlgoDescription () override
 
- Public Member Functions inherited from EcalUncalibRecHitWorkerBaseClass
 EcalUncalibRecHitWorkerBaseClass (const edm::ParameterSet &, edm::ConsumesCollector &c)
 
 EcalUncalibRecHitWorkerBaseClass (const edm::ParameterSet &)
 
 EcalUncalibRecHitWorkerBaseClass ()
 
virtual ~EcalUncalibRecHitWorkerBaseClass ()
 

Private Types

enum  TimeAlgo { noMethod, ratioMethod, weightsMethod, crossCorrelationMethod }
 

Private Member Functions

const SampleMatrix & noisecor (bool barrel, int gain) const
 
const SampleMatrixGainArray & noisecor (bool barrel) const
 
void run (const edm::Event &evt, const EcalDigiCollection &digis, EcalUncalibratedRecHitCollection &result) override
 
void set (const edm::EventSetup &es) override
 
void set (const edm::Event &evt) override
 
double timeCorrection (float ampli, const std::vector< float > &amplitudeBins, const std::vector< float > &shiftBins)
 

Private Attributes

BXVector activeBX
 
double addPedestalUncertaintyEB_
 
double addPedestalUncertaintyEE_
 
bool ampErrorCalculation_
 
double amplitudeThreshEB_
 
double amplitudeThreshEE_
 
edm::EDGetTokenT< unsigned int > bunchSpacing_
 
int bunchSpacingManual_
 
double CCminTimeToBeLateMax_
 
double CCminTimeToBeLateMin_
 
double CCtargetTimePrecision_
 
double CCtargetTimePrecisionForDelayedPulses_
 
double CCTimeShiftWrtRations_
 
double chi2ThreshEB_
 
double chi2ThreshEE_
 
std::unique_ptr< EcalUncalibRecHitTimingCCAlgo > computeCC_
 
bool doPrefitEB_
 
bool doPrefitEE_
 
bool dynamicPedestalsEB_
 
bool dynamicPedestalsEE_
 
std::vector< double > EBamplitudeFitParameters_
 
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 > eePulseShape_
 
double EEtimeConstantTerm_
 
std::pair< double, double > EEtimeFitLimits_
 
std::vector< double > EEtimeFitParameters_
 
double EEtimeNconst_
 
edm::ESHandle< EcalGainRatios > gains
 
edm::ESGetToken< EcalGainRatios, EcalGainRatiosRcd > gainsToken_
 
bool gainSwitchUseMaxSampleEB_
 
bool gainSwitchUseMaxSampleEE_
 
edm::ESHandle< EcalWeightXtalGroups > grps
 
edm::ESGetToken< EcalWeightXtalGroups, EcalWeightXtalGroupsRcd > grpsToken_
 
edm::ESHandle< EcalTimeCalibConstants > itime
 
edm::ESGetToken< EcalTimeCalibConstants, EcalTimeCalibConstantsRcd > itimeToken_
 
bool kPoorRecoFlagEB_
 
bool kPoorRecoFlagEE_
 
bool mitigateBadSamplesEB_
 
bool mitigateBadSamplesEE_
 
EcalUncalibRecHitMultiFitAlgo multiFitMethod_
 
edm::ESGetToken< EcalSamplesCorrelation, EcalSamplesCorrelationRcd > noiseConvariancesToken_
 
std::array< SampleMatrixGainArray, 2 > noisecors_
 
edm::ESHandle< EcalSamplesCorrelation > noisecovariances
 
edm::ESHandle< EcalTimeOffsetConstant > offtime
 
edm::ESGetToken< EcalTimeOffsetConstant, EcalTimeOffsetConstantRcd > offtimeToken_
 
double outOfTimeThreshG12mEB_
 
double outOfTimeThreshG12mEE_
 
double outOfTimeThreshG12pEB_
 
double outOfTimeThreshG12pEE_
 
double outOfTimeThreshG61mEB_
 
double outOfTimeThreshG61mEE_
 
double outOfTimeThreshG61pEB_
 
double outOfTimeThreshG61pEE_
 
edm::ESHandle< EcalPedestals > peds
 
edm::ESGetToken< EcalPedestals, EcalPedestalsRcd > pedsToken_
 
double prefitMaxChiSqEB_
 
double prefitMaxChiSqEE_
 
edm::ESGetToken< EcalPulseCovariances, EcalPulseCovariancesRcd > pulseConvariancesToken_
 
edm::ESHandle< EcalPulseCovariances > pulsecovariances
 
edm::ESHandle< EcalPulseShapes > pulseshapes
 
edm::ESGetToken< EcalPulseShapes, EcalPulseShapesRcd > pulseShapesToken_
 
EcalUncalibRecHitRatioMethodAlgo< EBDataFrame > ratioMethod_barrel_
 
EcalUncalibRecHitRatioMethodAlgo< EEDataFrame > ratioMethod_endcap_
 
edm::ESHandle< EcalSampleMask > sampleMaskHand_
 
edm::ESGetToken< EcalSampleMask, EcalSampleMaskRcd > sampleMaskToken_
 
bool selectiveBadSampleCriteriaEB_
 
bool selectiveBadSampleCriteriaEE_
 
bool simplifiedNoiseModelForGainSwitch_
 
TimeAlgo timealgo_ = noMethod
 
edm::ESHandle< EcalTimeBiasCorrections > timeCorrBias_
 
edm::ESGetToken< EcalTimeBiasCorrections, EcalTimeBiasCorrectionsRcd > timeCorrBiasToken_
 
bool useLumiInfoRunHeader_
 
const EcalWeightSet::EcalWeightMatrix * weights [2]
 
EcalUncalibRecHitTimeWeightsAlgo< EBDataFrame > weightsMethod_barrel_
 
EcalUncalibRecHitTimeWeightsAlgo< EEDataFrame > weightsMethod_endcap_
 
edm::ESHandle< EcalTBWeights > wgts
 
edm::ESGetToken< EcalTBWeights, EcalTBWeightsRcd > wgtsToken_
 

Detailed Description

Definition at line 50 of file EcalUncalibRecHitWorkerMultiFit.cc.

Member Enumeration Documentation

◆ TimeAlgo

enum EcalUncalibRecHitWorkerMultiFit::TimeAlgo
private
Enumerator
noMethod 
ratioMethod 
weightsMethod 
crossCorrelationMethod 

Definition at line 95 of file EcalUncalibRecHitWorkerMultiFit.cc.

95 { noMethod, ratioMethod, weightsMethod, crossCorrelationMethod };

Constructor & Destructor Documentation

◆ EcalUncalibRecHitWorkerMultiFit() [1/2]

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

Definition at line 174 of file EcalUncalibRecHitWorkerMultiFit.cc.

References activeBX, HLT_2023v12_cff::activeBXs, addPedestalUncertaintyEB_, addPedestalUncertaintyEE_, ampErrorCalculation_, amplitudeThreshEB_, amplitudeThreshEE_, bunchSpacing_, bunchSpacingManual_, HltBtagPostValidation_cff::c, CCminTimeToBeLateMax_, CCminTimeToBeLateMin_, CCtargetTimePrecision_, CCtargetTimePrecisionForDelayedPulses_, CCTimeShiftWrtRations_, ecalcctiming::clockToNS, computeCC_, crossCorrelationMethod, doPrefitEB_, doPrefitEE_, dynamicPedestalsEB_, dynamicPedestalsEE_, EBamplitudeFitParameters_, EBtimeConstantTerm_, EBtimeFitLimits_, EBtimeFitParameters_, EBtimeNconst_, EEamplitudeFitParameters_, EEtimeConstantTerm_, EEtimeFitLimits_, EEtimeFitParameters_, EEtimeNconst_, gainsToken_, gainSwitchUseMaxSampleEB_, gainSwitchUseMaxSampleEE_, edm::ParameterSet::getParameter(), grpsToken_, ProducerED_cfi::InputTag, itimeToken_, mitigateBadSamplesEB_, mitigateBadSamplesEE_, noiseConvariancesToken_, offtimeToken_, outOfTimeThreshG12mEB_, outOfTimeThreshG12mEE_, outOfTimeThreshG12pEB_, outOfTimeThreshG12pEE_, outOfTimeThreshG61mEB_, outOfTimeThreshG61mEE_, outOfTimeThreshG61pEB_, outOfTimeThreshG61pEE_, pedsToken_, prefitMaxChiSqEB_, prefitMaxChiSqEE_, pulseConvariancesToken_, pulseShapesToken_, ratioMethod, BXVector< T >::resize(), sampleMaskToken_, selectiveBadSampleCriteriaEB_, selectiveBadSampleCriteriaEE_, simplifiedNoiseModelForGainSwitch_, beamSpotPI::startTime, AlCaHLTBitMon_QueryRunRegistry::string, timealgo_, timeCorrBiasToken_, useLumiInfoRunHeader_, weightsMethod, and wgtsToken_.

175  : EcalUncalibRecHitWorkerBaseClass(ps, c) {
176  // get the BX for the pulses to be activated
177  std::vector<int32_t> activeBXs = ps.getParameter<std::vector<int32_t>>("activeBXs");
178  activeBX.resize(activeBXs.size());
179  for (unsigned int ibx = 0; ibx < activeBXs.size(); ++ibx) {
180  activeBX.coeffRef(ibx) = activeBXs[ibx];
181  }
182 
183  // uncertainty calculation (CPU intensive)
184  ampErrorCalculation_ = ps.getParameter<bool>("ampErrorCalculation");
185  useLumiInfoRunHeader_ = ps.getParameter<bool>("useLumiInfoRunHeader");
186 
187  if (useLumiInfoRunHeader_) {
188  bunchSpacing_ = c.consumes<unsigned int>(edm::InputTag("bunchSpacingProducer"));
189  bunchSpacingManual_ = 0;
190  } else {
191  bunchSpacingManual_ = ps.getParameter<int>("bunchSpacing");
192  }
193 
194  doPrefitEB_ = ps.getParameter<bool>("doPrefitEB");
195  doPrefitEE_ = ps.getParameter<bool>("doPrefitEE");
196 
197  prefitMaxChiSqEB_ = ps.getParameter<double>("prefitMaxChiSqEB");
198  prefitMaxChiSqEE_ = ps.getParameter<double>("prefitMaxChiSqEE");
199 
200  dynamicPedestalsEB_ = ps.getParameter<bool>("dynamicPedestalsEB");
201  dynamicPedestalsEE_ = ps.getParameter<bool>("dynamicPedestalsEE");
202  mitigateBadSamplesEB_ = ps.getParameter<bool>("mitigateBadSamplesEB");
203  mitigateBadSamplesEE_ = ps.getParameter<bool>("mitigateBadSamplesEE");
204  gainSwitchUseMaxSampleEB_ = ps.getParameter<bool>("gainSwitchUseMaxSampleEB");
205  gainSwitchUseMaxSampleEE_ = ps.getParameter<bool>("gainSwitchUseMaxSampleEE");
206  selectiveBadSampleCriteriaEB_ = ps.getParameter<bool>("selectiveBadSampleCriteriaEB");
207  selectiveBadSampleCriteriaEE_ = ps.getParameter<bool>("selectiveBadSampleCriteriaEE");
208  addPedestalUncertaintyEB_ = ps.getParameter<double>("addPedestalUncertaintyEB");
209  addPedestalUncertaintyEE_ = ps.getParameter<double>("addPedestalUncertaintyEE");
210  simplifiedNoiseModelForGainSwitch_ = ps.getParameter<bool>("simplifiedNoiseModelForGainSwitch");
211  pedsToken_ = c.esConsumes<EcalPedestals, EcalPedestalsRcd>();
212  gainsToken_ = c.esConsumes<EcalGainRatios, EcalGainRatiosRcd>();
213  noiseConvariancesToken_ = c.esConsumes<EcalSamplesCorrelation, EcalSamplesCorrelationRcd>();
214  pulseShapesToken_ = c.esConsumes<EcalPulseShapes, EcalPulseShapesRcd>();
215  pulseConvariancesToken_ = c.esConsumes<EcalPulseCovariances, EcalPulseCovariancesRcd>();
216  sampleMaskToken_ = c.esConsumes<EcalSampleMask, EcalSampleMaskRcd>();
217  grpsToken_ = c.esConsumes<EcalWeightXtalGroups, EcalWeightXtalGroupsRcd>();
218  wgtsToken_ = c.esConsumes<EcalTBWeights, EcalTBWeightsRcd>();
219  timeCorrBiasToken_ = c.esConsumes<EcalTimeBiasCorrections, EcalTimeBiasCorrectionsRcd>();
220  itimeToken_ =
221  c.esConsumes<EcalTimeCalibConstants, EcalTimeCalibConstantsRcd>(ps.getParameter<edm::ESInputTag>("timeCalibTag"));
222  offtimeToken_ = c.esConsumes<EcalTimeOffsetConstant, EcalTimeOffsetConstantRcd>(
223  ps.getParameter<edm::ESInputTag>("timeOffsetTag"));
224 
225  // algorithm to be used for timing
226  auto const& timeAlgoName = ps.getParameter<std::string>("timealgo");
227  if (timeAlgoName == "RatioMethod")
228  timealgo_ = ratioMethod;
229  else if (timeAlgoName == "WeightsMethod")
230  timealgo_ = weightsMethod;
231  else if (timeAlgoName == "crossCorrelationMethod") {
232  timealgo_ = crossCorrelationMethod;
233  double startTime = ps.getParameter<double>("crossCorrelationStartTime");
234  double stopTime = ps.getParameter<double>("crossCorrelationStopTime");
235  CCtargetTimePrecision_ = ps.getParameter<double>("crossCorrelationTargetTimePrecision");
236  CCtargetTimePrecisionForDelayedPulses_ =
237  ps.getParameter<double>("crossCorrelationTargetTimePrecisionForDelayedPulses");
238  CCminTimeToBeLateMin_ = ps.getParameter<double>("crossCorrelationMinTimeToBeLateMin") / ecalcctiming::clockToNS;
239  CCminTimeToBeLateMax_ = ps.getParameter<double>("crossCorrelationMinTimeToBeLateMax") / ecalcctiming::clockToNS;
240  CCTimeShiftWrtRations_ = ps.getParameter<double>("crossCorrelationTimeShiftWrtRations");
241  computeCC_ = std::make_unique<EcalUncalibRecHitTimingCCAlgo>(startTime, stopTime);
242  } else if (timeAlgoName != "None")
243  edm::LogError("EcalUncalibRecHitError") << "No time estimation algorithm defined";
244 
245  // ratio method parameters
246  EBtimeFitParameters_ = ps.getParameter<std::vector<double>>("EBtimeFitParameters");
247  EEtimeFitParameters_ = ps.getParameter<std::vector<double>>("EEtimeFitParameters");
248  EBamplitudeFitParameters_ = ps.getParameter<std::vector<double>>("EBamplitudeFitParameters");
249  EEamplitudeFitParameters_ = ps.getParameter<std::vector<double>>("EEamplitudeFitParameters");
250  EBtimeFitLimits_.first = ps.getParameter<double>("EBtimeFitLimits_Lower");
251  EBtimeFitLimits_.second = ps.getParameter<double>("EBtimeFitLimits_Upper");
252  EEtimeFitLimits_.first = ps.getParameter<double>("EEtimeFitLimits_Lower");
253  EEtimeFitLimits_.second = ps.getParameter<double>("EEtimeFitLimits_Upper");
254  EBtimeConstantTerm_ = ps.getParameter<double>("EBtimeConstantTerm");
255  EEtimeConstantTerm_ = ps.getParameter<double>("EEtimeConstantTerm");
256  EBtimeNconst_ = ps.getParameter<double>("EBtimeNconst");
257  EEtimeNconst_ = ps.getParameter<double>("EEtimeNconst");
258  outOfTimeThreshG12pEB_ = ps.getParameter<double>("outOfTimeThresholdGain12pEB");
259  outOfTimeThreshG12mEB_ = ps.getParameter<double>("outOfTimeThresholdGain12mEB");
260  outOfTimeThreshG61pEB_ = ps.getParameter<double>("outOfTimeThresholdGain61pEB");
261  outOfTimeThreshG61mEB_ = ps.getParameter<double>("outOfTimeThresholdGain61mEB");
262  outOfTimeThreshG12pEE_ = ps.getParameter<double>("outOfTimeThresholdGain12pEE");
263  outOfTimeThreshG12mEE_ = ps.getParameter<double>("outOfTimeThresholdGain12mEE");
264  outOfTimeThreshG61pEE_ = ps.getParameter<double>("outOfTimeThresholdGain61pEE");
265  outOfTimeThreshG61mEE_ = ps.getParameter<double>("outOfTimeThresholdGain61mEE");
266  amplitudeThreshEB_ = ps.getParameter<double>("amplitudeThresholdEB");
267  amplitudeThreshEE_ = ps.getParameter<double>("amplitudeThresholdEE");
268 }
EcalUncalibRecHitWorkerMultiFit::pulseConvariancesToken_
edm::ESGetToken< EcalPulseCovariances, EcalPulseCovariancesRcd > pulseConvariancesToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:73
EcalPulseShapes
EcalPulseShapesMap EcalPulseShapes
Definition: EcalPulseShapes.h:23
ProducerED_cfi.InputTag
InputTag
Definition: ProducerED_cfi.py:5
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:307
EcalUncalibRecHitWorkerMultiFit::pulseShapesToken_
edm::ESGetToken< EcalPulseShapes, EcalPulseShapesRcd > pulseShapesToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:71
EcalUncalibRecHitWorkerMultiFit::gainsToken_
edm::ESGetToken< EcalGainRatios, EcalGainRatiosRcd > gainsToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:67
EcalWeightXtalGroups
EcalCondObjectContainer< EcalXtalGroupId > EcalWeightXtalGroups
Definition: EcalWeightXtalGroups.h:12
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
ecalcctiming::clockToNS
constexpr const float clockToNS
Definition: EcalConstants.h:10
EcalUncalibRecHitWorkerMultiFit::noiseConvariancesToken_
edm::ESGetToken< EcalSamplesCorrelation, EcalSamplesCorrelationRcd > noiseConvariancesToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:69
EcalUncalibRecHitWorkerMultiFit::bunchSpacing_
edm::EDGetTokenT< unsigned int > bunchSpacing_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:88
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:35
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
EcalUncalibRecHitWorkerMultiFit::pedsToken_
edm::ESGetToken< EcalPedestals, EcalPedestalsRcd > pedsToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:65
EcalUncalibRecHitWorkerMultiFit::wgtsToken_
edm::ESGetToken< EcalTBWeights, EcalTBWeightsRcd > wgtsToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:102
EcalGainRatios
EcalGainRatioMap EcalGainRatios
Definition: EcalGainRatios.h:12
EcalPulseCovariances
EcalPulseCovariancesMap EcalPulseCovariances
Definition: EcalPulseCovariances.h:22
EcalUncalibRecHitWorkerMultiFit::computeCC_
std::unique_ptr< EcalUncalibRecHitTimingCCAlgo > computeCC_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:166
EcalPedestals
EcalPedestalsMap EcalPedestals
Definition: EcalPedestals.h:50
EcalUncalibRecHitWorkerMultiFit::offtimeToken_
edm::ESGetToken< EcalTimeOffsetConstant, EcalTimeOffsetConstantRcd > offtimeToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:155
EcalTimeCalibConstants
EcalTimeCalibConstantMap EcalTimeCalibConstants
Definition: EcalTimeCalibConstants.h:12
HLT_2023v12_cff.activeBXs
activeBXs
Definition: HLT_2023v12_cff.py:6378
BXVector::resize
void resize(int bx, unsigned size)
EcalUncalibRecHitWorkerMultiFit::timeCorrBiasToken_
edm::ESGetToken< EcalTimeBiasCorrections, EcalTimeBiasCorrectionsRcd > timeCorrBiasToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:150
EcalUncalibRecHitWorkerMultiFit::sampleMaskToken_
edm::ESGetToken< EcalSampleMask, EcalSampleMaskRcd > sampleMaskToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:92
EcalUncalibRecHitWorkerMultiFit::itimeToken_
edm::ESGetToken< EcalTimeCalibConstants, EcalTimeCalibConstantsRcd > itimeToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:153
EcalUncalibRecHitWorkerMultiFit::grpsToken_
edm::ESGetToken< EcalWeightXtalGroups, EcalWeightXtalGroupsRcd > grpsToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:100

◆ EcalUncalibRecHitWorkerMultiFit() [2/2]

EcalUncalibRecHitWorkerMultiFit::EcalUncalibRecHitWorkerMultiFit ( )
inline

Definition at line 53 of file EcalUncalibRecHitWorkerMultiFit.cc.

53 {};

Member Function Documentation

◆ getAlgoDescription()

edm::ParameterSetDescription EcalUncalibRecHitWorkerMultiFit::getAlgoDescription ( )
overridevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 722 of file EcalUncalibRecHitWorkerMultiFit.cc.

References edm::ParameterSetDescription::addNode().

722  {
723  edm::ParameterSetDescription psd;
724  psd.addNode(edm::ParameterDescription<std::vector<int>>("activeBXs", {-5, -4, -3, -2, -1, 0, 1, 2, 3, 4}, true) and
725  edm::ParameterDescription<bool>("ampErrorCalculation", true, true) and
726  edm::ParameterDescription<bool>("useLumiInfoRunHeader", true, true) and
727  edm::ParameterDescription<int>("bunchSpacing", 0, true) and
728  edm::ParameterDescription<bool>("doPrefitEB", false, true) and
729  edm::ParameterDescription<bool>("doPrefitEE", false, true) and
730  edm::ParameterDescription<double>("prefitMaxChiSqEB", 25., true) and
731  edm::ParameterDescription<double>("prefitMaxChiSqEE", 10., true) and
732  edm::ParameterDescription<bool>("dynamicPedestalsEB", false, true) and
733  edm::ParameterDescription<bool>("dynamicPedestalsEE", false, true) and
734  edm::ParameterDescription<bool>("mitigateBadSamplesEB", false, true) and
735  edm::ParameterDescription<bool>("mitigateBadSamplesEE", false, true) and
736  edm::ParameterDescription<bool>("gainSwitchUseMaxSampleEB", true, true) and
737  edm::ParameterDescription<bool>("gainSwitchUseMaxSampleEE", false, true) and
738  edm::ParameterDescription<bool>("selectiveBadSampleCriteriaEB", false, true) and
739  edm::ParameterDescription<bool>("selectiveBadSampleCriteriaEE", false, true) and
740  edm::ParameterDescription<double>("addPedestalUncertaintyEB", 0., true) and
741  edm::ParameterDescription<double>("addPedestalUncertaintyEE", 0., true) and
742  edm::ParameterDescription<bool>("simplifiedNoiseModelForGainSwitch", true, true) and
743  edm::ParameterDescription<std::string>("timealgo", "RatioMethod", true) and
744  edm::ParameterDescription<std::vector<double>>("EBtimeFitParameters",
745  {-2.015452e+00,
746  3.130702e+00,
747  -1.234730e+01,
748  4.188921e+01,
749  -8.283944e+01,
750  9.101147e+01,
751  -5.035761e+01,
752  1.105621e+01},
753  true) and
754  edm::ParameterDescription<std::vector<double>>("EEtimeFitParameters",
755  {-2.390548e+00,
756  3.553628e+00,
757  -1.762341e+01,
758  6.767538e+01,
759  -1.332130e+02,
760  1.407432e+02,
761  -7.541106e+01,
762  1.620277e+01},
763  true) and
764  edm::ParameterDescription<std::vector<double>>("EBamplitudeFitParameters", {1.138, 1.652}, true) and
765  edm::ParameterDescription<std::vector<double>>("EEamplitudeFitParameters", {1.890, 1.400}, true) and
766  edm::ParameterDescription<edm::ESInputTag>("timeCalibTag", edm::ESInputTag(), true) and
767  edm::ParameterDescription<edm::ESInputTag>("timeOffsetTag", edm::ESInputTag(), true) and
768  edm::ParameterDescription<double>("EBtimeFitLimits_Lower", 0.2, true) and
769  edm::ParameterDescription<double>("EBtimeFitLimits_Upper", 1.4, true) and
770  edm::ParameterDescription<double>("EEtimeFitLimits_Lower", 0.2, true) and
771  edm::ParameterDescription<double>("EEtimeFitLimits_Upper", 1.4, true) and
772  edm::ParameterDescription<double>("EBtimeConstantTerm", .6, true) and
773  edm::ParameterDescription<double>("EEtimeConstantTerm", 1.0, true) and
774  edm::ParameterDescription<double>("EBtimeNconst", 28.5, true) and
775  edm::ParameterDescription<double>("EEtimeNconst", 31.8, true) and
776  edm::ParameterDescription<double>("outOfTimeThresholdGain12pEB", 5., true) and
777  edm::ParameterDescription<double>("outOfTimeThresholdGain12mEB", 5., true) and
778  edm::ParameterDescription<double>("outOfTimeThresholdGain61pEB", 5., true) and
779  edm::ParameterDescription<double>("outOfTimeThresholdGain61mEB", 5., true) and
780  edm::ParameterDescription<double>("outOfTimeThresholdGain12pEE", 1000, true) and
781  edm::ParameterDescription<double>("outOfTimeThresholdGain12mEE", 1000, true) and
782  edm::ParameterDescription<double>("outOfTimeThresholdGain61pEE", 1000, true) and
783  edm::ParameterDescription<double>("outOfTimeThresholdGain61mEE", 1000, true) and
784  edm::ParameterDescription<double>("amplitudeThresholdEB", 10, true) and
785  edm::ParameterDescription<double>("amplitudeThresholdEE", 10, true) and
786  edm::ParameterDescription<double>("crossCorrelationStartTime", -25.0, true) and
787  edm::ParameterDescription<double>("crossCorrelationStopTime", 25.0, true) and
788  edm::ParameterDescription<double>("crossCorrelationTargetTimePrecision", 0.01, true) and
789  edm::ParameterDescription<double>("crossCorrelationTargetTimePrecisionForDelayedPulses", 0.05, true) and
790  edm::ParameterDescription<double>("crossCorrelationTimeShiftWrtRations", 0., true) and
791  edm::ParameterDescription<double>("crossCorrelationMinTimeToBeLateMin", 2., true) and
792  edm::ParameterDescription<double>("crossCorrelationMinTimeToBeLateMax", 5., true));
793 
794  return psd;
795 }
edm::ParameterSetDescription::addNode
ParameterDescriptionNode * addNode(ParameterDescriptionNode const &node)
Definition: ParameterSetDescription.cc:41

◆ noisecor() [1/2]

const SampleMatrix& EcalUncalibRecHitWorkerMultiFit::noisecor ( bool  barrel,
int  gain 
) const
inlineprivate

Definition at line 77 of file EcalUncalibRecHitWorkerMultiFit.cc.

References Reference_intrackfit_cff::barrel, PedestalClient_cfi::gain, and noisecors_.

Referenced by run().

77 { return noisecors_[barrel ? 1 : 0][gain]; }
EcalUncalibRecHitWorkerMultiFit::noisecors_
std::array< SampleMatrixGainArray, 2 > noisecors_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:81
PedestalClient_cfi.gain
gain
Definition: PedestalClient_cfi.py:37
Reference_intrackfit_cff.barrel
list barrel
Definition: Reference_intrackfit_cff.py:37

◆ noisecor() [2/2]

const SampleMatrixGainArray& EcalUncalibRecHitWorkerMultiFit::noisecor ( bool  barrel) const
inlineprivate

Definition at line 78 of file EcalUncalibRecHitWorkerMultiFit.cc.

References Reference_intrackfit_cff::barrel, and noisecors_.

78 { return noisecors_[barrel ? 1 : 0]; }
EcalUncalibRecHitWorkerMultiFit::noisecors_
std::array< SampleMatrixGainArray, 2 > noisecors_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:81
Reference_intrackfit_cff.barrel
list barrel
Definition: Reference_intrackfit_cff.py:37

◆ run()

void EcalUncalibRecHitWorkerMultiFit::run ( const edm::Event &  evt,
const EcalDigiCollection &  digis,
EcalUncalibratedRecHitCollection &  result 
)
overrideprivatevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 396 of file EcalUncalibRecHitWorkerMultiFit.cc.

References funct::abs(), activeBX, addPedestalUncertaintyEB_, addPedestalUncertaintyEE_, CustomPhysics_cfi::amplitude, amplitudeThreshEB_, amplitudeThreshEE_, Reference_intrackfit_cff::barrel, edm::DataFrameContainer::begin(), CCtargetTimePrecision_, CCtargetTimePrecisionForDelayedPulses_, CCTimeShiftWrtRations_, ecalcctiming::clockToNS, EcalUncalibRecHitRatioMethodAlgo< C >::computeAmplitude(), computeCC_, EcalUncalibRecHitRatioMethodAlgo< C >::computeTime(), crossCorrelationMethod, doPrefitEB_, doPrefitEE_, dynamicPedestalsEB_, dynamicPedestalsEE_, EBamplitudeFitParameters_, EBtimeConstantTerm_, EcalTimeBiasCorrections::EBTimeCorrAmplitudeBins, EcalTimeBiasCorrections::EBTimeCorrShiftBins, EBtimeFitLimits_, EBtimeFitParameters_, EBtimeNconst_, EcalBarrel, EcalEndcap, EEamplitudeFitParameters_, EEtimeConstantTerm_, EcalTimeBiasCorrections::EETimeCorrAmplitudeBins, EcalTimeBiasCorrections::EETimeCorrShiftBins, EEtimeFitLimits_, EEtimeFitParameters_, EEtimeNconst_, edm::DataFrameContainer::empty(), edm::DataFrameContainer::end(), EcalUncalibRecHitRatioMethodAlgo< C >::fixMGPAslew(), EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), ecalLiteDTU::gainId(), gains, gainSwitchUseMaxSampleEB_, gainSwitchUseMaxSampleEE_, EcalUncalibRecHitRatioMethodAlgo< C >::getCalculatedRecHit(), EcalTimeOffsetConstant::getEBValue(), EcalTimeOffsetConstant::getEEValue(), EcalTBWeights::getMap(), EcalWeightSet::getWeightsAfterGainSwitch(), EcalWeightSet::getWeightsBeforeGainSwitch(), grps, ecalpyutils::hashedIndex(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), Matriplex::hypot(), mps_fire::i, EcalXtalGroupId::id(), EcalUncalibRecHitRatioMethodAlgo< C >::init(), itime, dqmiolumiharvest::j, EcalUncalibratedRecHit::kHasSwitchToGain1, EcalUncalibratedRecHit::kHasSwitchToGain6, EcalUncalibratedRecHit::kOutOfTime, EcalUncalibratedRecHit::kSaturated, EcalUncalibRecHitMultiFitAlgo::makeRecHit(), EcalDataFrame::MAXSAMPLES, mitigateBadSamplesEB_, mitigateBadSamplesEE_, multiFitMethod_, noisecor(), offtime, outOfTimeThreshG12mEB_, outOfTimeThreshG12mEE_, outOfTimeThreshG12pEB_, outOfTimeThreshG12pEE_, outOfTimeThreshG61mEB_, outOfTimeThreshG61mEE_, outOfTimeThreshG61pEB_, outOfTimeThreshG61pEE_, EcalCondDBWriter_cfi::pedestal, peds, conifer::pow(), prefitMaxChiSqEB_, prefitMaxChiSqEE_, edm::ESHandle< T >::product(), pulsecovariances, pulseshapes, ratioMethod, ratioMethod_barrel_, ratioMethod_endcap_, mps_fire::result, ecalGpuTask_cfi::sample, sampleMaskHand_, selectiveBadSampleCriteriaEB_, selectiveBadSampleCriteriaEE_, EcalUncalibRecHitMultiFitAlgo::setAddPedestalUncertainty(), EcalUncalibRecHitMultiFitAlgo::setDoPrefit(), EcalUncalibRecHitMultiFitAlgo::setDynamicPedestals(), EcalUncalibRecHitMultiFitAlgo::setGainSwitchUseMaxSample(), EcalUncalibRecHitMultiFitAlgo::setMitigateBadSamples(), EcalUncalibRecHitMultiFitAlgo::setPrefitMaxChiSq(), EcalUncalibRecHitMultiFitAlgo::setSelectiveBadSampleCriteria(), EcalUncalibRecHitMultiFitAlgo::setSimplifiedNoiseModelForGainSwitch(), simplifiedNoiseModelForGainSwitch_, BXVector< T >::size(), edm::DataFrameContainer::size(), mathSSE::sqrt(), EcalPulseShape::TEMPLATESAMPLES, remoteMonitoring_LASER_era2018_cfg::threshold, EcalUncalibRecHitTimeWeightsAlgo< C >::time(), timealgo_, timeCorrBias_, timeCorrection(), weights, weightsMethod, weightsMethod_barrel_, weightsMethod_endcap_, and wgts.

398  {
399  if (digis.empty())
400  return;
401 
402  // assume all digis come from the same subdetector (either barrel or endcap)
403  DetId detid(digis.begin()->id());
404  bool barrel = (detid.subdetId() == EcalBarrel);
405 
406  multiFitMethod_.setSimplifiedNoiseModelForGainSwitch(simplifiedNoiseModelForGainSwitch_);
407  if (barrel) {
408  multiFitMethod_.setDoPrefit(doPrefitEB_);
409  multiFitMethod_.setPrefitMaxChiSq(prefitMaxChiSqEB_);
410  multiFitMethod_.setDynamicPedestals(dynamicPedestalsEB_);
411  multiFitMethod_.setMitigateBadSamples(mitigateBadSamplesEB_);
412  multiFitMethod_.setGainSwitchUseMaxSample(gainSwitchUseMaxSampleEB_);
413  multiFitMethod_.setSelectiveBadSampleCriteria(selectiveBadSampleCriteriaEB_);
414  multiFitMethod_.setAddPedestalUncertainty(addPedestalUncertaintyEB_);
415  } else {
416  multiFitMethod_.setDoPrefit(doPrefitEE_);
417  multiFitMethod_.setPrefitMaxChiSq(prefitMaxChiSqEE_);
418  multiFitMethod_.setDynamicPedestals(dynamicPedestalsEE_);
419  multiFitMethod_.setMitigateBadSamples(mitigateBadSamplesEE_);
420  multiFitMethod_.setGainSwitchUseMaxSample(gainSwitchUseMaxSampleEE_);
421  multiFitMethod_.setSelectiveBadSampleCriteria(selectiveBadSampleCriteriaEE_);
422  multiFitMethod_.setAddPedestalUncertainty(addPedestalUncertaintyEE_);
423  }
424 
425  FullSampleVector fullpulse(FullSampleVector::Zero());
426  FullSampleMatrix fullpulsecov(FullSampleMatrix::Zero());
427 
428  result.reserve(result.size() + digis.size());
429  for (auto itdg = digis.begin(); itdg != digis.end(); ++itdg) {
430  DetId detid(itdg->id());
431 
432  const EcalSampleMask* sampleMask_ = sampleMaskHand_.product();
433 
434  // intelligence for recHit computation
435  float offsetTime = 0;
436 
437  const EcalPedestals::Item* aped = nullptr;
438  const EcalMGPAGainRatio* aGain = nullptr;
439  const EcalXtalGroupId* gid = nullptr;
440  const EcalPulseShapes::Item* aPulse = nullptr;
441  const EcalPulseCovariances::Item* aPulseCov = nullptr;
442 
443  if (barrel) {
444  unsigned int hashedIndex = EBDetId(detid).hashedIndex();
445  aped = &peds->barrel(hashedIndex);
446  aGain = &gains->barrel(hashedIndex);
447  gid = &grps->barrel(hashedIndex);
448  aPulse = &pulseshapes->barrel(hashedIndex);
449  aPulseCov = &pulsecovariances->barrel(hashedIndex);
450  offsetTime = offtime->getEBValue();
451  } else {
452  unsigned int hashedIndex = EEDetId(detid).hashedIndex();
453  aped = &peds->endcap(hashedIndex);
454  aGain = &gains->endcap(hashedIndex);
455  gid = &grps->endcap(hashedIndex);
456  aPulse = &pulseshapes->endcap(hashedIndex);
457  aPulseCov = &pulsecovariances->endcap(hashedIndex);
458  offsetTime = offtime->getEEValue();
459  }
460 
461  double pedVec[3] = {aped->mean_x12, aped->mean_x6, aped->mean_x1};
462  double pedRMSVec[3] = {aped->rms_x12, aped->rms_x6, aped->rms_x1};
463  double gainRatios[3] = {1., aGain->gain12Over6(), aGain->gain6Over1() * aGain->gain12Over6()};
464 
465  for (int i = 0; i < EcalPulseShape::TEMPLATESAMPLES; ++i)
466  fullpulse(i + 7) = aPulse->pdfval[i];
467 
468  for (int i = 0; i < EcalPulseShape::TEMPLATESAMPLES; i++)
469  for (int j = 0; j < EcalPulseShape::TEMPLATESAMPLES; j++)
470  fullpulsecov(i + 7, j + 7) = aPulseCov->covval[i][j];
471 
472  // compute the right bin of the pulse shape using time calibration constants
473  EcalTimeCalibConstantMap::const_iterator it = itime->find(detid);
474  EcalTimeCalibConstant itimeconst = 0;
475  if (it != itime->end()) {
476  itimeconst = (*it);
477  } else {
478  edm::LogError("EcalRecHitError") << "No time intercalib const found for xtal " << detid.rawId()
479  << "! something wrong with EcalTimeCalibConstants in your DB? ";
480  }
481 
482  int lastSampleBeforeSaturation = -2;
483  for (unsigned int iSample = 0; iSample < EcalDataFrame::MAXSAMPLES; iSample++) {
484  if (((EcalDataFrame)(*itdg)).sample(iSample).gainId() == 0) {
485  lastSampleBeforeSaturation = iSample - 1;
486  break;
487  }
488  }
489 
490  // === amplitude computation ===
491 
492  if (lastSampleBeforeSaturation == 4) { // saturation on the expected max sample
493  result.emplace_back((*itdg).id(), 4095 * 12, 0, 0, 0);
494  auto& uncalibRecHit = result.back();
495  uncalibRecHit.setFlagBit(EcalUncalibratedRecHit::kSaturated);
496  // do not propagate the default chi2 = -1 value to the calib rechit (mapped to 64), set it to 0 when saturation
497  uncalibRecHit.setChi2(0);
498  } else if (lastSampleBeforeSaturation >=
499  -1) { // saturation on other samples: cannot extrapolate from the fourth one
500  int gainId = ((EcalDataFrame)(*itdg)).sample(5).gainId();
501  if (gainId == 0)
502  gainId = 3;
503  auto pedestal = pedVec[gainId - 1];
504  auto gainratio = gainRatios[gainId - 1];
505  double amplitude = ((double)(((EcalDataFrame)(*itdg)).sample(5).adc()) - pedestal) * gainratio;
506  result.emplace_back((*itdg).id(), amplitude, 0, 0, 0);
507  auto& uncalibRecHit = result.back();
508  uncalibRecHit.setFlagBit(EcalUncalibratedRecHit::kSaturated);
509  // do not propagate the default chi2 = -1 value to the calib rechit (mapped to 64), set it to 0 when saturation
510  uncalibRecHit.setChi2(0);
511  } else {
512  // multifit
513  const SampleMatrixGainArray& noisecors = noisecor(barrel);
514 
515  result.push_back(multiFitMethod_.makeRecHit(*itdg, aped, aGain, noisecors, fullpulse, fullpulsecov, activeBX));
516  auto& uncalibRecHit = result.back();
517 
518  // === time computation ===
519  if (timealgo_ == ratioMethod) {
520  // ratio method
521  constexpr float clockToNsConstant = 25.;
522  constexpr float invClockToNs = 1. / clockToNsConstant;
523  if (not barrel) {
524  ratioMethod_endcap_.init(*itdg, *sampleMask_, pedVec, pedRMSVec, gainRatios);
525  ratioMethod_endcap_.computeTime(EEtimeFitParameters_, EEtimeFitLimits_, EEamplitudeFitParameters_);
526  ratioMethod_endcap_.computeAmplitude(EEamplitudeFitParameters_);
527  EcalUncalibRecHitRatioMethodAlgo<EEDataFrame>::CalculatedRecHit crh =
528  ratioMethod_endcap_.getCalculatedRecHit();
529  double theTimeCorrectionEE = timeCorrection(
530  uncalibRecHit.amplitude(), timeCorrBias_->EETimeCorrAmplitudeBins, timeCorrBias_->EETimeCorrShiftBins);
531 
532  uncalibRecHit.setJitter(crh.timeMax - 5 + theTimeCorrectionEE);
533  uncalibRecHit.setJitterError(
534  std::sqrt(std::pow(crh.timeError, 2) + std::pow(EEtimeConstantTerm_ * invClockToNs, 2)));
535 
536  // consider flagging as kOutOfTime only if above noise
537  if (uncalibRecHit.amplitude() > pedRMSVec[0] * amplitudeThreshEE_) {
538  float outOfTimeThreshP = outOfTimeThreshG12pEE_;
539  float outOfTimeThreshM = outOfTimeThreshG12mEE_;
540  // determine if gain has switched away from gainId==1 (x12 gain)
541  // and determine cuts (number of 'sigmas') to ose for kOutOfTime
542  // >3k ADC is necessasry condition for gain switch to occur
543  if (uncalibRecHit.amplitude() > 3000.) {
544  for (int iSample = 0; iSample < EEDataFrame::MAXSAMPLES; iSample++) {
545  int GainId = ((EcalDataFrame)(*itdg)).sample(iSample).gainId();
546  if (GainId != 1) {
547  outOfTimeThreshP = outOfTimeThreshG61pEE_;
548  outOfTimeThreshM = outOfTimeThreshG61mEE_;
549  break;
550  }
551  }
552  }
553  float correctedTime = (crh.timeMax - 5) * clockToNsConstant + itimeconst + offsetTime;
554  float cterm = EEtimeConstantTerm_;
555  float sigmaped = pedRMSVec[0]; // approx for lower gains
556  float nterm = EEtimeNconst_ * sigmaped / uncalibRecHit.amplitude();
557  float sigmat = std::sqrt(nterm * nterm + cterm * cterm);
558  if ((correctedTime > sigmat * outOfTimeThreshP) || (correctedTime < -sigmat * outOfTimeThreshM)) {
559  uncalibRecHit.setFlagBit(EcalUncalibratedRecHit::kOutOfTime);
560  }
561  }
562 
563  } else {
564  ratioMethod_barrel_.init(*itdg, *sampleMask_, pedVec, pedRMSVec, gainRatios);
565  ratioMethod_barrel_.fixMGPAslew(*itdg);
566  ratioMethod_barrel_.computeTime(EBtimeFitParameters_, EBtimeFitLimits_, EBamplitudeFitParameters_);
567  ratioMethod_barrel_.computeAmplitude(EBamplitudeFitParameters_);
568  EcalUncalibRecHitRatioMethodAlgo<EBDataFrame>::CalculatedRecHit crh =
569  ratioMethod_barrel_.getCalculatedRecHit();
570 
571  double theTimeCorrectionEB = timeCorrection(
572  uncalibRecHit.amplitude(), timeCorrBias_->EBTimeCorrAmplitudeBins, timeCorrBias_->EBTimeCorrShiftBins);
573 
574  uncalibRecHit.setJitter(crh.timeMax - 5 + theTimeCorrectionEB);
575  uncalibRecHit.setJitterError(std::hypot(crh.timeError, EBtimeConstantTerm_ / clockToNsConstant));
576 
577  // consider flagging as kOutOfTime only if above noise
578  if (uncalibRecHit.amplitude() > pedRMSVec[0] * amplitudeThreshEB_) {
579  float outOfTimeThreshP = outOfTimeThreshG12pEB_;
580  float outOfTimeThreshM = outOfTimeThreshG12mEB_;
581  // determine if gain has switched away from gainId==1 (x12 gain)
582  // and determine cuts (number of 'sigmas') to ose for kOutOfTime
583  // >3k ADC is necessasry condition for gain switch to occur
584  if (uncalibRecHit.amplitude() > 3000.) {
585  for (int iSample = 0; iSample < EBDataFrame::MAXSAMPLES; iSample++) {
586  int GainId = ((EcalDataFrame)(*itdg)).sample(iSample).gainId();
587  if (GainId != 1) {
588  outOfTimeThreshP = outOfTimeThreshG61pEB_;
589  outOfTimeThreshM = outOfTimeThreshG61mEB_;
590  break;
591  }
592  }
593  }
594  float correctedTime = (crh.timeMax - 5) * clockToNsConstant + itimeconst + offsetTime;
595  float cterm = EBtimeConstantTerm_;
596  float sigmaped = pedRMSVec[0]; // approx for lower gains
597  float nterm = EBtimeNconst_ * sigmaped / uncalibRecHit.amplitude();
598  float sigmat = std::sqrt(nterm * nterm + cterm * cterm);
599  if ((correctedTime > sigmat * outOfTimeThreshP) || (correctedTime < -sigmat * outOfTimeThreshM)) {
600  uncalibRecHit.setFlagBit(EcalUncalibratedRecHit::kOutOfTime);
601  }
602  }
603  }
604  } else if (timealgo_ == weightsMethod) {
605  // weights method on the PU subtracted pulse shape
606  std::vector<double> amplitudes;
607  for (unsigned int ibx = 0; ibx < activeBX.size(); ++ibx)
608  amplitudes.push_back(uncalibRecHit.outOfTimeAmplitude(ibx));
609 
610  EcalTBWeights::EcalTDCId tdcid(1);
611  EcalTBWeights::EcalTBWeightMap const& wgtsMap = wgts->getMap();
612  EcalTBWeights::EcalTBWeightMap::const_iterator wit;
613  wit = wgtsMap.find(std::make_pair(*gid, tdcid));
614  if (wit == wgtsMap.end()) {
615  edm::LogError("EcalUncalibRecHitError")
616  << "No weights found for EcalGroupId: " << gid->id() << " and EcalTDCId: " << tdcid
617  << "\n skipping digi with id: " << detid.rawId();
618  result.pop_back();
619  continue;
620  }
621  const EcalWeightSet& wset = wit->second; // this is the EcalWeightSet
622 
623  const EcalWeightSet::EcalWeightMatrix& mat1 = wset.getWeightsBeforeGainSwitch();
624  const EcalWeightSet::EcalWeightMatrix& mat2 = wset.getWeightsAfterGainSwitch();
625 
626  weights[0] = &mat1;
627  weights[1] = &mat2;
628 
629  double timerh;
630  if (detid.subdetId() == EcalEndcap) {
631  timerh = weightsMethod_endcap_.time(*itdg, amplitudes, aped, aGain, fullpulse, weights);
632  } else {
633  timerh = weightsMethod_barrel_.time(*itdg, amplitudes, aped, aGain, fullpulse, weights);
634  }
635  uncalibRecHit.setJitter(timerh);
636  uncalibRecHit.setJitterError(0.); // not computed with weights
637 
638  } else if (timealgo_ == crossCorrelationMethod) {
639  std::vector<double> amplitudes(activeBX.size());
640  for (unsigned int ibx = 0; ibx < activeBX.size(); ++ibx)
641  amplitudes[ibx] = uncalibRecHit.outOfTimeAmplitude(ibx);
642 
643  float jitter =
644  computeCC_->computeTimeCC(*itdg, amplitudes, aped, aGain, fullpulse, CCtargetTimePrecision_, true) +
645  CCTimeShiftWrtRations_ / ecalcctiming::clockToNS;
646  float noCorrectedJitter =
647  computeCC_->computeTimeCC(
648  *itdg, amplitudes, aped, aGain, fullpulse, CCtargetTimePrecisionForDelayedPulses_, false) +
649  CCTimeShiftWrtRations_ / ecalcctiming::clockToNS;
650 
651  uncalibRecHit.setJitter(jitter);
652  uncalibRecHit.setNonCorrectedTime(jitter, noCorrectedJitter);
653 
654  float retreivedNonCorrectedTime = uncalibRecHit.nonCorrectedTime();
655  float noCorrectedTime = ecalcctiming::clockToNS * noCorrectedJitter;
656  if (retreivedNonCorrectedTime > -29.0 && std::abs(retreivedNonCorrectedTime - noCorrectedTime) > 0.05) {
657  edm::LogError("EcalUncalibRecHitError") << "Problem with noCorrectedJitter: true value:" << noCorrectedTime
658  << "\t received: " << retreivedNonCorrectedTime << std::endl;
659  } //<<>>if (abs(retreivedNonCorrectedTime - noCorrectedJitter)>1);
660 
661  // consider flagging as kOutOfTime only if above noise
662  float threshold, cterm, timeNconst;
663  float timeThrP = 0.;
664  float timeThrM = 0.;
665  if (barrel) {
666  threshold = pedRMSVec[0] * amplitudeThreshEB_;
667  cterm = EBtimeConstantTerm_;
668  timeNconst = EBtimeNconst_;
669  timeThrP = outOfTimeThreshG12pEB_;
670  timeThrM = outOfTimeThreshG12mEB_;
671  if (uncalibRecHit.amplitude() > 3000.) { // Gain switch
672  for (int iSample = 0; iSample < EBDataFrame::MAXSAMPLES; iSample++) {
673  int GainId = ((EcalDataFrame)(*itdg)).sample(iSample).gainId();
674  if (GainId != 1) {
675  timeThrP = outOfTimeThreshG61pEB_;
676  timeThrM = outOfTimeThreshG61mEB_;
677  break;
678  }
679  }
680  }
681  } else { //EndCap
682  threshold = pedRMSVec[0] * amplitudeThreshEE_;
683  cterm = EEtimeConstantTerm_;
684  timeNconst = EEtimeNconst_;
685  timeThrP = outOfTimeThreshG12pEE_;
686  timeThrM = outOfTimeThreshG12mEE_;
687  if (uncalibRecHit.amplitude() > 3000.) { // Gain switch
688  for (int iSample = 0; iSample < EEDataFrame::MAXSAMPLES; iSample++) {
689  int GainId = ((EcalDataFrame)(*itdg)).sample(iSample).gainId();
690  if (GainId != 1) {
691  timeThrP = outOfTimeThreshG61pEE_;
692  timeThrM = outOfTimeThreshG61mEE_;
693  break;
694  }
695  }
696  }
697  }
698  if (uncalibRecHit.amplitude() > threshold) {
699  float correctedTime = noCorrectedJitter * ecalcctiming::clockToNS + itimeconst + offsetTime;
700  float sigmaped = pedRMSVec[0]; // approx for lower gains
701  float nterm = timeNconst * sigmaped / uncalibRecHit.amplitude();
702  float sigmat = std::sqrt(nterm * nterm + cterm * cterm);
703  if ((correctedTime > sigmat * timeThrP) || (correctedTime < -sigmat * timeThrM))
704  uncalibRecHit.setFlagBit(EcalUncalibratedRecHit::kOutOfTime);
705  }
706 
707  } else { // no time method;
708  uncalibRecHit.setJitter(0.);
709  uncalibRecHit.setJitterError(0.);
710  }
711  }
712 
713  // set flags if gain switch has occurred
714  auto& uncalibRecHit = result.back();
715  if (((EcalDataFrame)(*itdg)).hasSwitchToGain6())
716  uncalibRecHit.setFlagBit(EcalUncalibratedRecHit::kHasSwitchToGain6);
717  if (((EcalDataFrame)(*itdg)).hasSwitchToGain1())
718  uncalibRecHit.setFlagBit(EcalUncalibratedRecHit::kHasSwitchToGain1);
719  }
720 }
Matriplex::hypot
MPlex< T, D1, D2, N > hypot(const MPlex< T, D1, D2, N > &a, const MPlex< T, D1, D2, N > &b)
Definition: Matriplex.h:417
EcalUncalibRecHitRatioMethodAlgo::computeAmplitude
void computeAmplitude(std::vector< double > &amplitudeFitParameters)
Definition: EcalUncalibRecHitRatioMethodAlgo.h:527
mps_fire.i
i
Definition: mps_fire.py:429
EcalTimeBiasCorrections::EBTimeCorrShiftBins
std::vector< float > EBTimeCorrShiftBins
Definition: EcalTimeBiasCorrections.h:24
EcalUncalibRecHitWorkerMultiFit::multiFitMethod_
EcalUncalibRecHitMultiFitAlgo multiFitMethod_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:85
EcalTBWeights::getMap
const EcalTBWeightMap & getMap() const
Definition: EcalTBWeights.h:28
mps_fire.result
result
Definition: mps_fire.py:311
remoteMonitoring_LASER_era2018_cfg.threshold
threshold
Definition: remoteMonitoring_LASER_era2018_cfg.py:581
SampleMatrixGainArray
std::array< SampleMatrix, NGains > SampleMatrixGainArray
Definition: EigenMatrixTypes.h:28
EcalUncalibRecHitWorkerMultiFit::weightsMethod_endcap_
EcalUncalibRecHitTimeWeightsAlgo< EEDataFrame > weightsMethod_endcap_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:105
FullSampleMatrix
Eigen::Matrix< double, FullSampleVectorSize, FullSampleVectorSize > FullSampleMatrix
Definition: EigenMatrixTypes.h:18
EcalUncalibRecHitTimeWeightsAlgo::time
double time(const C &dataFrame, const std::vector< double > &amplitudes, const EcalPedestals::Item *aped, const EcalMGPAGainRatio *aGain, const FullSampleVector &fullpulse, const EcalWeightSet::EcalWeightMatrix **weights)
Compute time.
Definition: EcalUncalibRecHitTimeWeightsAlgo.h:33
FullSampleVector
Eigen::Matrix< double, FullSampleVectorSize, 1 > FullSampleVector
Definition: EigenMatrixTypes.h:13
EcalTBWeights::EcalTBWeightMap
std::map< std::pair< EcalXtalGroupId, EcalTDCId >, EcalWeightSet > EcalTBWeightMap
Definition: EcalTBWeights.h:18
EcalUncalibRecHitRatioMethodAlgo::computeTime
void computeTime(std::vector< double > &timeFitParameters, std::pair< double, double > &timeFitLimits, std::vector< double > &amplitudeFitParameters)
Definition: EcalUncalibRecHitRatioMethodAlgo.h:211
EcalPulseShape::TEMPLATESAMPLES
static const int TEMPLATESAMPLES
Definition: EcalPulseShapes.h:10
EcalUncalibRecHitWorkerMultiFit::grps
edm::ESHandle< EcalWeightXtalGroups > grps
Definition: EcalUncalibRecHitWorkerMultiFit.cc:99
conifer::pow
constexpr int pow(int x)
Definition: conifer.h:24
EcalWeightSet::EcalWeightMatrix
math::Matrix< 3, 10 >::type EcalWeightMatrix
Definition: EcalWeightSet.h:19
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129
ecalcctiming::clockToNS
constexpr const float clockToNS
Definition: EcalConstants.h:10
CustomPhysics_cfi.amplitude
amplitude
Definition: CustomPhysics_cfi.py:12
EcalCondDBWriter_cfi.pedestal
pedestal
Definition: EcalCondDBWriter_cfi.py:49
EcalUncalibRecHitWorkerMultiFit::ratioMethod_barrel_
EcalUncalibRecHitRatioMethodAlgo< EBDataFrame > ratioMethod_barrel_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:130
ecalGpuTask_cfi.sample
sample
Definition: ecalGpuTask_cfi.py:48
BXVector::size
unsigned size(int bx) const
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:36
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
EcalUncalibRecHitWorkerMultiFit::ratioMethod_endcap_
EcalUncalibRecHitRatioMethodAlgo< EEDataFrame > ratioMethod_endcap_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:131
EcalTimeBiasCorrections::EETimeCorrAmplitudeBins
std::vector< float > EETimeCorrAmplitudeBins
Definition: EcalTimeBiasCorrections.h:26
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
EcalUncalibRecHitWorkerMultiFit::pulseshapes
edm::ESHandle< EcalPulseShapes > pulseshapes
Definition: EcalUncalibRecHitWorkerMultiFit.cc:70
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
EcalUncalibRecHitWorkerMultiFit::itime
edm::ESHandle< EcalTimeCalibConstants > itime
Definition: EcalUncalibRecHitWorkerMultiFit.cc:152
EcalUncalibRecHitWorkerMultiFit::timeCorrection
double timeCorrection(float ampli, const std::vector< float > &amplitudeBins, const std::vector< float > &shiftBins)
Definition: EcalUncalibRecHitWorkerMultiFit.cc:349
EcalXtalGroupId::id
unsigned int id() const
Definition: EcalXtalGroupId.h:22
EcalUncalibRecHitWorkerMultiFit::computeCC_
std::unique_ptr< EcalUncalibRecHitTimingCCAlgo > computeCC_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:166
EcalWeightSet::getWeightsAfterGainSwitch
EcalWeightMatrix & getWeightsAfterGainSwitch()
Definition: EcalWeightSet.h:27
EcalUncalibRecHitWorkerMultiFit::noisecor
const SampleMatrix & noisecor(bool barrel, int gain) const
Definition: EcalUncalibRecHitWorkerMultiFit.cc:77
EcalUncalibRecHitMultiFitAlgo::makeRecHit
EcalUncalibratedRecHit makeRecHit(const EcalDataFrame &dataFrame, const EcalPedestals::Item *aped, const EcalMGPAGainRatio *aGain, const SampleMatrixGainArray &noisecors, const FullSampleVector &fullpulse, const FullSampleMatrix &fullpulsecov, const BXVector &activeBX)
compute rechits
Definition: EcalUncalibRecHitMultiFitAlgo.cc:27
edm::DataFrameContainer::end
const_iterator end() const
Definition: DataFrameContainer.h:146
EcalUncalibRecHitWorkerMultiFit::timeCorrBias_
edm::ESHandle< EcalTimeBiasCorrections > timeCorrBias_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:149
EcalMGPAGainRatio::gain12Over6
float gain12Over6() const
Definition: EcalMGPAGainRatio.h:19
DetId
Definition: DetId.h:17
EcalUncalibRecHitWorkerMultiFit::pulsecovariances
edm::ESHandle< EcalPulseCovariances > pulsecovariances
Definition: EcalUncalibRecHitWorkerMultiFit.cc:72
ecalLiteDTU::gainId
constexpr int gainId(sample_type sample)
get the gainId (2 bits)
Definition: EcalLiteDTUSample.h:14
EcalWeightSet::getWeightsBeforeGainSwitch
EcalWeightMatrix & getWeightsBeforeGainSwitch()
Definition: EcalWeightSet.h:26
edm::DataFrameContainer::begin
const_iterator begin() const
The iterator returned can not safely be used across threads.
Definition: DataFrameContainer.h:143
EcalCondObjectContainer::const_iterator
std::vector< Item >::const_iterator const_iterator
Definition: EcalCondObjectContainer.h:19
EcalMGPAGainRatio::gain6Over1
float gain6Over1() const
Definition: EcalMGPAGainRatio.h:20
EcalUncalibRecHitWorkerMultiFit::weightsMethod_barrel_
EcalUncalibRecHitTimeWeightsAlgo< EBDataFrame > weightsMethod_barrel_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:104
EcalUncalibRecHitWorkerMultiFit::offtime
edm::ESHandle< EcalTimeOffsetConstant > offtime
Definition: EcalUncalibRecHitWorkerMultiFit.cc:154
EcalTimeCalibConstant
float EcalTimeCalibConstant
Definition: EcalTimeCalibConstants.h:10
EBDetId::hashedIndex
int hashedIndex() const
get a compact index for arrays
Definition: EBDetId.h:82
EcalDataFrame::MAXSAMPLES
static constexpr int MAXSAMPLES
Definition: EcalDataFrame.h:48
EcalUncalibRecHitWorkerMultiFit::sampleMaskHand_
edm::ESHandle< EcalSampleMask > sampleMaskHand_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:91
EcalUncalibRecHitWorkerMultiFit::weights
const EcalWeightSet::EcalWeightMatrix * weights[2]
Definition: EcalUncalibRecHitWorkerMultiFit.cc:103
EcalUncalibRecHitRatioMethodAlgo::init
void init(const C &dataFrame, const EcalSampleMask &sampleMask, const double *pedestals, const double *pedestalRMSes, const double *gainRatios)
Definition: EcalUncalibRecHitRatioMethodAlgo.h:96
EEDetId::hashedIndex
int hashedIndex() const
Definition: EEDetId.h:183
Reference_intrackfit_cff.barrel
list barrel
Definition: Reference_intrackfit_cff.py:37

◆ set() [1/2]

void EcalUncalibRecHitWorkerMultiFit::set ( const edm::EventSetup &  es)
overrideprivatevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 270 of file EcalUncalibRecHitWorkerMultiFit.cc.

References funct::abs(), ampErrorCalculation_, EcalUncalibRecHitMultiFitAlgo::disableErrorCalculation(), EcalSamplesCorrelation::EBG12SamplesCorrelation, EcalSamplesCorrelation::EBG1SamplesCorrelation, EcalSamplesCorrelation::EBG6SamplesCorrelation, EcalSamplesCorrelation::EEG12SamplesCorrelation, EcalSamplesCorrelation::EEG1SamplesCorrelation, EcalSamplesCorrelation::EEG6SamplesCorrelation, gains, gainsToken_, edm::EventSetup::getHandle(), grps, grpsToken_, mps_fire::i, itime, itimeToken_, dqmiolumiharvest::j, multiFitMethod_, noiseConvariancesToken_, noisecors_, noisecovariances, offtime, offtimeToken_, peds, pedsToken_, pulseConvariancesToken_, pulsecovariances, pulseshapes, pulseShapesToken_, sampleMaskHand_, sampleMaskToken_, timeCorrBias_, timeCorrBiasToken_, wgts, and wgtsToken_.

270  {
271  // common setup
272  gains = es.getHandle(gainsToken_);
273  peds = es.getHandle(pedsToken_);
274 
275  // for the multifit method
276  if (!ampErrorCalculation_)
277  multiFitMethod_.disableErrorCalculation();
278  noisecovariances = es.getHandle(noiseConvariancesToken_);
279  pulseshapes = es.getHandle(pulseShapesToken_);
280  pulsecovariances = es.getHandle(pulseConvariancesToken_);
281 
282  // weights parameters for the time
283  grps = es.getHandle(grpsToken_);
284  wgts = es.getHandle(wgtsToken_);
285 
286  // which of the samples need be used
287  sampleMaskHand_ = es.getHandle(sampleMaskToken_);
288 
289  // for the ratio method
290  itime = es.getHandle(itimeToken_);
291  offtime = es.getHandle(offtimeToken_);
292 
293  // for the time correction methods
294  timeCorrBias_ = es.getHandle(timeCorrBiasToken_);
295 
296  int nnoise = SampleVector::RowsAtCompileTime;
297  SampleMatrix& noisecorEBg12 = noisecors_[1][0];
298  SampleMatrix& noisecorEBg6 = noisecors_[1][1];
299  SampleMatrix& noisecorEBg1 = noisecors_[1][2];
300  SampleMatrix& noisecorEEg12 = noisecors_[0][0];
301  SampleMatrix& noisecorEEg6 = noisecors_[0][1];
302  SampleMatrix& noisecorEEg1 = noisecors_[0][2];
303 
304  for (int i = 0; i < nnoise; ++i) {
305  for (int j = 0; j < nnoise; ++j) {
306  int vidx = std::abs(j - i);
307  noisecorEBg12(i, j) = noisecovariances->EBG12SamplesCorrelation[vidx];
308  noisecorEEg12(i, j) = noisecovariances->EEG12SamplesCorrelation[vidx];
309  noisecorEBg6(i, j) = noisecovariances->EBG6SamplesCorrelation[vidx];
310  noisecorEEg6(i, j) = noisecovariances->EEG6SamplesCorrelation[vidx];
311  noisecorEBg1(i, j) = noisecovariances->EBG1SamplesCorrelation[vidx];
312  noisecorEEg1(i, j) = noisecovariances->EEG1SamplesCorrelation[vidx];
313  }
314  }
315 }
EcalUncalibRecHitWorkerMultiFit::pulseConvariancesToken_
edm::ESGetToken< EcalPulseCovariances, EcalPulseCovariancesRcd > pulseConvariancesToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:73
EcalUncalibRecHitWorkerMultiFit::noisecors_
std::array< SampleMatrixGainArray, 2 > noisecors_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:81
mps_fire.i
i
Definition: mps_fire.py:429
EcalUncalibRecHitWorkerMultiFit::multiFitMethod_
EcalUncalibRecHitMultiFitAlgo multiFitMethod_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:85
EcalUncalibRecHitWorkerMultiFit::pulseShapesToken_
edm::ESGetToken< EcalPulseShapes, EcalPulseShapesRcd > pulseShapesToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:71
EcalSamplesCorrelation::EBG12SamplesCorrelation
std::vector< double > EBG12SamplesCorrelation
Definition: EcalSamplesCorrelation.h:22
EcalUncalibRecHitWorkerMultiFit::gainsToken_
edm::ESGetToken< EcalGainRatios, EcalGainRatiosRcd > gainsToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:67
EcalUncalibRecHitWorkerMultiFit::grps
edm::ESHandle< EcalWeightXtalGroups > grps
Definition: EcalUncalibRecHitWorkerMultiFit.cc:99
EcalUncalibRecHitWorkerMultiFit::noiseConvariancesToken_
edm::ESGetToken< EcalSamplesCorrelation, EcalSamplesCorrelationRcd > noiseConvariancesToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:69
EcalSamplesCorrelation::EBG6SamplesCorrelation
std::vector< double > EBG6SamplesCorrelation
Definition: EcalSamplesCorrelation.h:23
EcalUncalibRecHitWorkerMultiFit::pedsToken_
edm::ESGetToken< EcalPedestals, EcalPedestalsRcd > pedsToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:65
EcalUncalibRecHitWorkerMultiFit::wgtsToken_
edm::ESGetToken< EcalTBWeights, EcalTBWeightsRcd > wgtsToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:102
EcalUncalibRecHitWorkerMultiFit::pulseshapes
edm::ESHandle< EcalPulseShapes > pulseshapes
Definition: EcalUncalibRecHitWorkerMultiFit.cc:70
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
EcalUncalibRecHitWorkerMultiFit::itime
edm::ESHandle< EcalTimeCalibConstants > itime
Definition: EcalUncalibRecHitWorkerMultiFit.cc:152
edm::EventSetup::getHandle
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
Definition: EventSetup.h:130
EcalUncalibRecHitWorkerMultiFit::timeCorrBias_
edm::ESHandle< EcalTimeBiasCorrections > timeCorrBias_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:149
EcalUncalibRecHitWorkerMultiFit::offtimeToken_
edm::ESGetToken< EcalTimeOffsetConstant, EcalTimeOffsetConstantRcd > offtimeToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:155
EcalUncalibRecHitWorkerMultiFit::pulsecovariances
edm::ESHandle< EcalPulseCovariances > pulsecovariances
Definition: EcalUncalibRecHitWorkerMultiFit.cc:72
EcalSamplesCorrelation::EBG1SamplesCorrelation
std::vector< double > EBG1SamplesCorrelation
Definition: EcalSamplesCorrelation.h:24
EcalUncalibRecHitWorkerMultiFit::timeCorrBiasToken_
edm::ESGetToken< EcalTimeBiasCorrections, EcalTimeBiasCorrectionsRcd > timeCorrBiasToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:150
SampleMatrix
Eigen::Matrix< double, SampleVectorSize, SampleVectorSize > SampleMatrix
Definition: EigenMatrixTypes.h:17
EcalUncalibRecHitWorkerMultiFit::offtime
edm::ESHandle< EcalTimeOffsetConstant > offtime
Definition: EcalUncalibRecHitWorkerMultiFit.cc:154
EcalSamplesCorrelation::EEG1SamplesCorrelation
std::vector< double > EEG1SamplesCorrelation
Definition: EcalSamplesCorrelation.h:28
EcalUncalibRecHitWorkerMultiFit::sampleMaskToken_
edm::ESGetToken< EcalSampleMask, EcalSampleMaskRcd > sampleMaskToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:92
EcalUncalibRecHitWorkerMultiFit::itimeToken_
edm::ESGetToken< EcalTimeCalibConstants, EcalTimeCalibConstantsRcd > itimeToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:153
EcalUncalibRecHitWorkerMultiFit::sampleMaskHand_
edm::ESHandle< EcalSampleMask > sampleMaskHand_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:91
EcalSamplesCorrelation::EEG12SamplesCorrelation
std::vector< double > EEG12SamplesCorrelation
Definition: EcalSamplesCorrelation.h:26
EcalUncalibRecHitWorkerMultiFit::noisecovariances
edm::ESHandle< EcalSamplesCorrelation > noisecovariances
Definition: EcalUncalibRecHitWorkerMultiFit.cc:68
EcalSamplesCorrelation::EEG6SamplesCorrelation
std::vector< double > EEG6SamplesCorrelation
Definition: EcalSamplesCorrelation.h:27
EcalUncalibRecHitWorkerMultiFit::grpsToken_
edm::ESGetToken< EcalWeightXtalGroups, EcalWeightXtalGroupsRcd > grpsToken_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:100

◆ set() [2/2]

void EcalUncalibRecHitWorkerMultiFit::set ( const edm::Event &  evt)
overrideprivatevirtual

Reimplemented from EcalUncalibRecHitWorkerBaseClass.

Definition at line 317 of file EcalUncalibRecHitWorkerMultiFit.cc.

References activeBX, bunchSpacing_, bunchSpacingManual_, edm::Event::getByToken(), BXVector< T >::resize(), and useLumiInfoRunHeader_.

317  {
318  unsigned int bunchspacing = 450;
319 
320  if (useLumiInfoRunHeader_) {
321  edm::Handle<unsigned int> bunchSpacingH;
322  evt.getByToken(bunchSpacing_, bunchSpacingH);
323  bunchspacing = *bunchSpacingH;
324  } else {
325  bunchspacing = bunchSpacingManual_;
326  }
327 
328  if (useLumiInfoRunHeader_ || bunchSpacingManual_ > 0) {
329  if (bunchspacing == 25) {
330  activeBX.resize(10);
331  activeBX << -5, -4, -3, -2, -1, 0, 1, 2, 3, 4;
332  } else {
333  //50ns configuration otherwise (also for no pileup)
334  activeBX.resize(5);
335  activeBX << -4, -2, 0, 2, 4;
336  }
337  }
338 }
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:528
EcalUncalibRecHitWorkerMultiFit::bunchSpacing_
edm::EDGetTokenT< unsigned int > bunchSpacing_
Definition: EcalUncalibRecHitWorkerMultiFit.cc:88
BXVector::resize
void resize(int bx, unsigned size)

◆ timeCorrection()

double EcalUncalibRecHitWorkerMultiFit::timeCorrection ( float  ampli,
const std::vector< float > &  amplitudeBins,
const std::vector< float > &  shiftBins 
)
private

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 349 of file EcalUncalibRecHitWorkerMultiFit.cc.

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

Referenced by run().

351  {
352  // computed initially in ns. Than turned in the BX's, as
353  // EcalUncalibratedRecHit need be.
354  double theCorrection = 0;
355 
356  // sanity check for arrays
357  if (amplitudeBins.empty()) {
358  edm::LogError("EcalRecHitError") << "timeCorrAmplitudeBins is empty, forcing no time bias corrections.";
359 
360  return 0;
361  }
362 
363  if (amplitudeBins.size() != shiftBins.size()) {
364  edm::LogError("EcalRecHitError") << "Size of timeCorrAmplitudeBins different from "
365  "timeCorrShiftBins. Forcing no time bias corrections. ";
366 
367  return 0;
368  }
369 
370  // FIXME? what about a binary search?
371  int myBin = -1;
372  for (int bin = 0; bin < (int)amplitudeBins.size(); bin++) {
373  if (ampli > amplitudeBins[bin]) {
374  myBin = bin;
375  } else {
376  break;
377  }
378  }
379 
380  if (myBin == -1) {
381  theCorrection = shiftBins[0];
382  } else if (myBin == ((int)(amplitudeBins.size() - 1))) {
383  theCorrection = shiftBins[myBin];
384  } else {
385  // interpolate linearly between two assingned points
386  theCorrection = (shiftBins[myBin + 1] - shiftBins[myBin]);
387  theCorrection *= (((double)ampli) - amplitudeBins[myBin]) / (amplitudeBins[myBin + 1] - amplitudeBins[myBin]);
388  theCorrection += shiftBins[myBin];
389  }
390 
391  // convert ns into clocks
392  constexpr double inv25 = 1. / 25.;
393  return theCorrection * inv25;
394 }
edm::LogError
Log< level::Error, false > LogError
Definition: MessageLogger.h:129

Member Data Documentation

◆ activeBX

BXVector EcalUncalibRecHitWorkerMultiFit::activeBX
private

Definition at line 82 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), run(), and set().

◆ addPedestalUncertaintyEB_

double EcalUncalibRecHitWorkerMultiFit::addPedestalUncertaintyEB_
private

Definition at line 118 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ addPedestalUncertaintyEE_

double EcalUncalibRecHitWorkerMultiFit::addPedestalUncertaintyEE_
private

Definition at line 119 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ ampErrorCalculation_

bool EcalUncalibRecHitWorkerMultiFit::ampErrorCalculation_
private

Definition at line 83 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ amplitudeThreshEB_

double EcalUncalibRecHitWorkerMultiFit::amplitudeThreshEB_
private

Definition at line 145 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ amplitudeThreshEE_

double EcalUncalibRecHitWorkerMultiFit::amplitudeThreshEE_
private

Definition at line 146 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ bunchSpacing_

edm::EDGetTokenT<unsigned int> EcalUncalibRecHitWorkerMultiFit::bunchSpacing_
private

Definition at line 88 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ bunchSpacingManual_

int EcalUncalibRecHitWorkerMultiFit::bunchSpacingManual_
private

Definition at line 87 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ CCminTimeToBeLateMax_

double EcalUncalibRecHitWorkerMultiFit::CCminTimeToBeLateMax_
private

Definition at line 168 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit().

◆ CCminTimeToBeLateMin_

double EcalUncalibRecHitWorkerMultiFit::CCminTimeToBeLateMin_
private

Definition at line 167 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit().

◆ CCtargetTimePrecision_

double EcalUncalibRecHitWorkerMultiFit::CCtargetTimePrecision_
private

Definition at line 170 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ CCtargetTimePrecisionForDelayedPulses_

double EcalUncalibRecHitWorkerMultiFit::CCtargetTimePrecisionForDelayedPulses_
private

Definition at line 171 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ CCTimeShiftWrtRations_

double EcalUncalibRecHitWorkerMultiFit::CCTimeShiftWrtRations_
private

Definition at line 169 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ chi2ThreshEB_

double EcalUncalibRecHitWorkerMultiFit::chi2ThreshEB_
private

Definition at line 162 of file EcalUncalibRecHitWorkerMultiFit.cc.

◆ chi2ThreshEE_

double EcalUncalibRecHitWorkerMultiFit::chi2ThreshEE_
private

Definition at line 163 of file EcalUncalibRecHitWorkerMultiFit.cc.

◆ computeCC_

std::unique_ptr<EcalUncalibRecHitTimingCCAlgo> EcalUncalibRecHitWorkerMultiFit::computeCC_
private

Definition at line 166 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ doPrefitEB_

bool EcalUncalibRecHitWorkerMultiFit::doPrefitEB_
private

Definition at line 106 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ doPrefitEE_

bool EcalUncalibRecHitWorkerMultiFit::doPrefitEE_
private

Definition at line 107 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ dynamicPedestalsEB_

bool EcalUncalibRecHitWorkerMultiFit::dynamicPedestalsEB_
private

Definition at line 110 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ dynamicPedestalsEE_

bool EcalUncalibRecHitWorkerMultiFit::dynamicPedestalsEE_
private

Definition at line 111 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ EBamplitudeFitParameters_

std::vector<double> EcalUncalibRecHitWorkerMultiFit::EBamplitudeFitParameters_
private

Definition at line 125 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ ebPulseShape_

std::vector<double> EcalUncalibRecHitWorkerMultiFit::ebPulseShape_
private

Definition at line 156 of file EcalUncalibRecHitWorkerMultiFit.cc.

◆ ebSpikeThresh_

double EcalUncalibRecHitWorkerMultiFit::ebSpikeThresh_
private

Definition at line 147 of file EcalUncalibRecHitWorkerMultiFit.cc.

◆ EBtimeConstantTerm_

double EcalUncalibRecHitWorkerMultiFit::EBtimeConstantTerm_
private

Definition at line 133 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ EBtimeFitLimits_

std::pair<double, double> EcalUncalibRecHitWorkerMultiFit::EBtimeFitLimits_
private

Definition at line 127 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ EBtimeFitParameters_

std::vector<double> EcalUncalibRecHitWorkerMultiFit::EBtimeFitParameters_
private

Definition at line 123 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ EBtimeNconst_

double EcalUncalibRecHitWorkerMultiFit::EBtimeNconst_
private

Definition at line 135 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ EEamplitudeFitParameters_

std::vector<double> EcalUncalibRecHitWorkerMultiFit::EEamplitudeFitParameters_
private

Definition at line 126 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ eePulseShape_

std::vector<double> EcalUncalibRecHitWorkerMultiFit::eePulseShape_
private

Definition at line 157 of file EcalUncalibRecHitWorkerMultiFit.cc.

◆ EEtimeConstantTerm_

double EcalUncalibRecHitWorkerMultiFit::EEtimeConstantTerm_
private

Definition at line 134 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ EEtimeFitLimits_

std::pair<double, double> EcalUncalibRecHitWorkerMultiFit::EEtimeFitLimits_
private

Definition at line 128 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ EEtimeFitParameters_

std::vector<double> EcalUncalibRecHitWorkerMultiFit::EEtimeFitParameters_
private

Definition at line 124 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ EEtimeNconst_

double EcalUncalibRecHitWorkerMultiFit::EEtimeNconst_
private

Definition at line 136 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ gains

edm::ESHandle<EcalGainRatios> EcalUncalibRecHitWorkerMultiFit::gains
private

Definition at line 66 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ gainsToken_

edm::ESGetToken<EcalGainRatios, EcalGainRatiosRcd> EcalUncalibRecHitWorkerMultiFit::gainsToken_
private

Definition at line 67 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ gainSwitchUseMaxSampleEB_

bool EcalUncalibRecHitWorkerMultiFit::gainSwitchUseMaxSampleEB_
private

Definition at line 114 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ gainSwitchUseMaxSampleEE_

bool EcalUncalibRecHitWorkerMultiFit::gainSwitchUseMaxSampleEE_
private

Definition at line 115 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ grps

edm::ESHandle<EcalWeightXtalGroups> EcalUncalibRecHitWorkerMultiFit::grps
private

Definition at line 99 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ grpsToken_

edm::ESGetToken<EcalWeightXtalGroups, EcalWeightXtalGroupsRcd> EcalUncalibRecHitWorkerMultiFit::grpsToken_
private

Definition at line 100 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ itime

edm::ESHandle<EcalTimeCalibConstants> EcalUncalibRecHitWorkerMultiFit::itime
private

Definition at line 152 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ itimeToken_

edm::ESGetToken<EcalTimeCalibConstants, EcalTimeCalibConstantsRcd> EcalUncalibRecHitWorkerMultiFit::itimeToken_
private

Definition at line 153 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ kPoorRecoFlagEB_

bool EcalUncalibRecHitWorkerMultiFit::kPoorRecoFlagEB_
private

Definition at line 160 of file EcalUncalibRecHitWorkerMultiFit.cc.

◆ kPoorRecoFlagEE_

bool EcalUncalibRecHitWorkerMultiFit::kPoorRecoFlagEE_
private

Definition at line 161 of file EcalUncalibRecHitWorkerMultiFit.cc.

◆ mitigateBadSamplesEB_

bool EcalUncalibRecHitWorkerMultiFit::mitigateBadSamplesEB_
private

Definition at line 112 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ mitigateBadSamplesEE_

bool EcalUncalibRecHitWorkerMultiFit::mitigateBadSamplesEE_
private

Definition at line 113 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ multiFitMethod_

EcalUncalibRecHitMultiFitAlgo EcalUncalibRecHitWorkerMultiFit::multiFitMethod_
private

Definition at line 85 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ noiseConvariancesToken_

edm::ESGetToken<EcalSamplesCorrelation, EcalSamplesCorrelationRcd> EcalUncalibRecHitWorkerMultiFit::noiseConvariancesToken_
private

Definition at line 69 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ noisecors_

std::array<SampleMatrixGainArray, 2> EcalUncalibRecHitWorkerMultiFit::noisecors_
private

Definition at line 81 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by noisecor(), and set().

◆ noisecovariances

edm::ESHandle<EcalSamplesCorrelation> EcalUncalibRecHitWorkerMultiFit::noisecovariances
private

Definition at line 68 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by set().

◆ offtime

edm::ESHandle<EcalTimeOffsetConstant> EcalUncalibRecHitWorkerMultiFit::offtime
private

Definition at line 154 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ offtimeToken_

edm::ESGetToken<EcalTimeOffsetConstant, EcalTimeOffsetConstantRcd> EcalUncalibRecHitWorkerMultiFit::offtimeToken_
private

Definition at line 155 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ outOfTimeThreshG12mEB_

double EcalUncalibRecHitWorkerMultiFit::outOfTimeThreshG12mEB_
private

Definition at line 138 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ outOfTimeThreshG12mEE_

double EcalUncalibRecHitWorkerMultiFit::outOfTimeThreshG12mEE_
private

Definition at line 142 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ outOfTimeThreshG12pEB_

double EcalUncalibRecHitWorkerMultiFit::outOfTimeThreshG12pEB_
private

Definition at line 137 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ outOfTimeThreshG12pEE_

double EcalUncalibRecHitWorkerMultiFit::outOfTimeThreshG12pEE_
private

Definition at line 141 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ outOfTimeThreshG61mEB_

double EcalUncalibRecHitWorkerMultiFit::outOfTimeThreshG61mEB_
private

Definition at line 140 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ outOfTimeThreshG61mEE_

double EcalUncalibRecHitWorkerMultiFit::outOfTimeThreshG61mEE_
private

Definition at line 144 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ outOfTimeThreshG61pEB_

double EcalUncalibRecHitWorkerMultiFit::outOfTimeThreshG61pEB_
private

Definition at line 139 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ outOfTimeThreshG61pEE_

double EcalUncalibRecHitWorkerMultiFit::outOfTimeThreshG61pEE_
private

Definition at line 143 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ peds

edm::ESHandle<EcalPedestals> EcalUncalibRecHitWorkerMultiFit::peds
private

Definition at line 64 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ pedsToken_

edm::ESGetToken<EcalPedestals, EcalPedestalsRcd> EcalUncalibRecHitWorkerMultiFit::pedsToken_
private

Definition at line 65 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ prefitMaxChiSqEB_

double EcalUncalibRecHitWorkerMultiFit::prefitMaxChiSqEB_
private

Definition at line 108 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ prefitMaxChiSqEE_

double EcalUncalibRecHitWorkerMultiFit::prefitMaxChiSqEE_
private

Definition at line 109 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ pulseConvariancesToken_

edm::ESGetToken<EcalPulseCovariances, EcalPulseCovariancesRcd> EcalUncalibRecHitWorkerMultiFit::pulseConvariancesToken_
private

Definition at line 73 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ pulsecovariances

edm::ESHandle<EcalPulseCovariances> EcalUncalibRecHitWorkerMultiFit::pulsecovariances
private

Definition at line 72 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ pulseshapes

edm::ESHandle<EcalPulseShapes> EcalUncalibRecHitWorkerMultiFit::pulseshapes
private

Definition at line 70 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ pulseShapesToken_

edm::ESGetToken<EcalPulseShapes, EcalPulseShapesRcd> EcalUncalibRecHitWorkerMultiFit::pulseShapesToken_
private

Definition at line 71 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ ratioMethod_barrel_

EcalUncalibRecHitRatioMethodAlgo<EBDataFrame> EcalUncalibRecHitWorkerMultiFit::ratioMethod_barrel_
private

Definition at line 130 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run().

◆ ratioMethod_endcap_

EcalUncalibRecHitRatioMethodAlgo<EEDataFrame> EcalUncalibRecHitWorkerMultiFit::ratioMethod_endcap_
private

Definition at line 131 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run().

◆ sampleMaskHand_

edm::ESHandle<EcalSampleMask> EcalUncalibRecHitWorkerMultiFit::sampleMaskHand_
private

Definition at line 91 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ sampleMaskToken_

edm::ESGetToken<EcalSampleMask, EcalSampleMaskRcd> EcalUncalibRecHitWorkerMultiFit::sampleMaskToken_
private

Definition at line 92 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ selectiveBadSampleCriteriaEB_

bool EcalUncalibRecHitWorkerMultiFit::selectiveBadSampleCriteriaEB_
private

Definition at line 116 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ selectiveBadSampleCriteriaEE_

bool EcalUncalibRecHitWorkerMultiFit::selectiveBadSampleCriteriaEE_
private

Definition at line 117 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ simplifiedNoiseModelForGainSwitch_

bool EcalUncalibRecHitWorkerMultiFit::simplifiedNoiseModelForGainSwitch_
private

Definition at line 120 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ timealgo_

TimeAlgo EcalUncalibRecHitWorkerMultiFit::timealgo_ = noMethod
private

Definition at line 96 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and run().

◆ timeCorrBias_

edm::ESHandle<EcalTimeBiasCorrections> EcalUncalibRecHitWorkerMultiFit::timeCorrBias_
private

Definition at line 149 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ timeCorrBiasToken_

edm::ESGetToken<EcalTimeBiasCorrections, EcalTimeBiasCorrectionsRcd> EcalUncalibRecHitWorkerMultiFit::timeCorrBiasToken_
private

Definition at line 150 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ useLumiInfoRunHeader_

bool EcalUncalibRecHitWorkerMultiFit::useLumiInfoRunHeader_
private

Definition at line 84 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

◆ weights

const EcalWeightSet::EcalWeightMatrix* EcalUncalibRecHitWorkerMultiFit::weights[2]
private

Definition at line 103 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run().

◆ weightsMethod_barrel_

EcalUncalibRecHitTimeWeightsAlgo<EBDataFrame> EcalUncalibRecHitWorkerMultiFit::weightsMethod_barrel_
private

Definition at line 104 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run().

◆ weightsMethod_endcap_

EcalUncalibRecHitTimeWeightsAlgo<EEDataFrame> EcalUncalibRecHitWorkerMultiFit::weightsMethod_endcap_
private

Definition at line 105 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run().

◆ wgts

edm::ESHandle<EcalTBWeights> EcalUncalibRecHitWorkerMultiFit::wgts
private

Definition at line 101 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by run(), and set().

◆ wgtsToken_

edm::ESGetToken<EcalTBWeights, EcalTBWeightsRcd> EcalUncalibRecHitWorkerMultiFit::wgtsToken_
private

Definition at line 102 of file EcalUncalibRecHitWorkerMultiFit.cc.

Referenced by EcalUncalibRecHitWorkerMultiFit(), and set().

Generated for CMSSW Reference Manual by doxygen 1.8.14