#include <EcalUncalibRecHitWorkerWeights.h>

Inheritance diagram for EcalUncalibRecHitWorkerWeights:

Public Member Functions
	EcalUncalibRecHitWorkerWeights (const edm::ParameterSet &, edm::ConsumesCollector &c)

	EcalUncalibRecHitWorkerWeights ()=default

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

	~EcalUncalibRecHitWorkerWeights () 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 Attributes
const EcalWeightSet::EcalChi2WeightMatrix *	chi2mat [2]

double	gainRatios [3]

edm::ESHandle< EcalGainRatios >	gains_

edm::ESHandle < EcalWeightXtalGroups >	grps_

double	pedRMSVec [3]

edm::ESHandle< EcalPedestals >	peds_

double	pedVec [3]

EBShape	testbeamEBShape

EEShape	testbeamEEShape

edm::ESGetToken < EcalGainRatios, EcalGainRatiosRcd >	tokenGains_

edm::ESGetToken < EcalWeightXtalGroups, EcalWeightXtalGroupsRcd >	tokenGrps_

edm::ESGetToken< EcalPedestals, EcalPedestalsRcd >	tokenPeds_

edm::ESGetToken< EcalTBWeights, EcalTBWeightsRcd >	tokenWgts_

EcalUncalibRecHitRecWeightsAlgo < EBDataFrame >	uncalibMaker_barrel_

EcalUncalibRecHitRecWeightsAlgo < EEDataFrame >	uncalibMaker_endcap_

const EcalWeightSet::EcalWeightMatrix *	weights [2]

edm::ESHandle< EcalTBWeights >	wgts_

Detailed Description

Definition at line 34 of file EcalUncalibRecHitWorkerWeights.h.

Constructor & Destructor Documentation

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

Definition at line 12 of file EcalUncalibRecHitWorkerWeights.cc.

     : EcalUncalibRecHitWorkerRunOneDigiBase(ps, c),
       tokenPeds_(c.esConsumes<EcalPedestals, EcalPedestalsRcd>()),
       tokenGains_(c.esConsumes<EcalGainRatios, EcalGainRatiosRcd>()),
       tokenGrps_(c.esConsumes<EcalWeightXtalGroups, EcalWeightXtalGroupsRcd>()),
       tokenWgts_(c.esConsumes<EcalTBWeights, EcalTBWeightsRcd>()),
       testbeamEEShape(c),
       testbeamEBShape(c) {}

EcalUncalibRecHitWorkerWeights::EcalUncalibRecHitWorkerWeights ( )

default

EcalUncalibRecHitWorkerWeights::~EcalUncalibRecHitWorkerWeights ( )

inlineoverride

Definition at line 38 of file EcalUncalibRecHitWorkerWeights.h.

38 {};

Member Function Documentation

edm::ParameterSetDescription EcalUncalibRecHitWorkerWeights::getAlgoDescription ( )

overridevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 113 of file EcalUncalibRecHitWorkerWeights.cc.

                                                                               {
   edm::ParameterSetDescription psd;
   return psd;
 }

bool EcalUncalibRecHitWorkerWeights::run	(	const edm::Event &	evt,
		const EcalDigiCollection::const_iterator &	digi,
		EcalUncalibratedRecHitCollection &	result
	)

overridevirtual

Implements EcalUncalibRecHitWorkerRunOneDigiBase.

Definition at line 31 of file EcalUncalibRecHitWorkerWeights.cc.

References EcalEndcap, EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), gainRatios, gains_, EcalWeightSet::getWeightsAfterGainSwitch(), EcalWeightSet::getWeightsBeforeGainSwitch(), grps_, ecalpyutils::hashedIndex(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), EcalXtalGroupId::id(), EcalUncalibratedRecHit::kHasSwitchToGain1, EcalUncalibratedRecHit::kHasSwitchToGain6, EcalUncalibRecHitRecWeightsAlgo< C >::makeRecHit(), pedRMSVec, peds_, pedVec, EcalUncalibratedRecHit::setFlagBit(), testbeamEBShape, testbeamEEShape, uncalibMaker_endcap_, weights, and wgts_.

                                                                                    {
   DetId detid(itdg->id());
 
   const EcalPedestals::Item* aped = nullptr;
   const EcalMGPAGainRatio* aGain = nullptr;
   const EcalXtalGroupId* gid = nullptr;
   EcalTBWeights::EcalTDCId tdcid(1);
 
   if (detid.subdetId() == EcalEndcap) {
     unsigned int hashedIndex = EEDetId(detid).hashedIndex();
     aped = &peds_->endcap(hashedIndex);
     aGain = &gains_->endcap(hashedIndex);
     gid = &grps_->endcap(hashedIndex);
   } else {
     unsigned int hashedIndex = EBDetId(detid).hashedIndex();
     aped = &peds_->barrel(hashedIndex);
     aGain = &gains_->barrel(hashedIndex);
     gid = &grps_->barrel(hashedIndex);
   }
 
   pedVec[0] = aped->mean_x12;
   pedVec[1] = aped->mean_x6;
   pedVec[2] = aped->mean_x1;
   pedRMSVec[0] = aped->rms_x12;
   pedRMSVec[1] = aped->rms_x6;
   pedRMSVec[2] = aped->rms_x1;
   gainRatios[0] = 1.;
   gainRatios[1] = aGain->gain12Over6();
   gainRatios[2] = aGain->gain6Over1() * aGain->gain12Over6();
 
   // now lookup the correct weights in the map
   EcalTBWeights::EcalTBWeightMap const& wgtsMap = wgts_->getMap();
   EcalTBWeights::EcalTBWeightMap::const_iterator wit;
   wit = wgtsMap.find(std::make_pair(*gid, tdcid));
   if (wit == wgtsMap.end()) {
     edm::LogError("EcalUncalibRecHitError")
         << "No weights found for EcalGroupId: " << gid->id() << " and  EcalTDCId: " << tdcid
         << "\n  skipping digi with id: " << detid.rawId();
 
     return false;
   }
   const EcalWeightSet& wset = wit->second;  // this is the EcalWeightSet
 
   const EcalWeightSet::EcalWeightMatrix& mat1 = wset.getWeightsBeforeGainSwitch();
   const EcalWeightSet::EcalWeightMatrix& mat2 = wset.getWeightsAfterGainSwitch();
   //        const EcalWeightSet::EcalChi2WeightMatrix& mat3 = wset.getChi2WeightsBeforeGainSwitch();
   //        const EcalWeightSet::EcalChi2WeightMatrix& mat4 = wset.getChi2WeightsAfterGainSwitch();
 
   weights[0] = &mat1;
   weights[1] = &mat2;
 
   //        chi2mat[0] = &mat3;
   //        chi2mat[1] = &mat4;
   /*
         if (detid.subdetId()==EcalEndcap) {
                 result.push_back(uncalibMaker_endcap_.makeRecHit(*itdg, pedVec, gainRatios, weights, chi2mat));
         } else {
                 result.push_back(uncalibMaker_barrel_.makeRecHit(*itdg, pedVec, gainRatios, weights, chi2mat));
         }
 */
   if (detid.subdetId() == EcalEndcap) {
     EcalUncalibratedRecHit rhit =
         (uncalibMaker_endcap_.makeRecHit(*itdg, pedVec, pedRMSVec, gainRatios, weights, testbeamEEShape));
     if (((EcalDataFrame)(*itdg)).hasSwitchToGain6())
       rhit.setFlagBit(EcalUncalibratedRecHit::kHasSwitchToGain6);
     if (((EcalDataFrame)(*itdg)).hasSwitchToGain1())
       rhit.setFlagBit(EcalUncalibratedRecHit::kHasSwitchToGain1);
     result.emplace_back(rhit);
   } else {
     EcalUncalibratedRecHit rhit =
         (uncalibMaker_endcap_.makeRecHit(*itdg, pedVec, pedRMSVec, gainRatios, weights, testbeamEBShape));
     if (((EcalDataFrame)(*itdg)).hasSwitchToGain6())
       rhit.setFlagBit(EcalUncalibratedRecHit::kHasSwitchToGain6);
     if (((EcalDataFrame)(*itdg)).hasSwitchToGain1())
       rhit.setFlagBit(EcalUncalibratedRecHit::kHasSwitchToGain1);
     result.emplace_back(rhit);
   }
   return true;
 }

void EcalUncalibRecHitWorkerWeights::set ( const edm::EventSetup & es )

overridevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 21 of file EcalUncalibRecHitWorkerWeights.cc.

References gains_, edm::EventSetup::getHandle(), grps_, peds_, EcalShapeBase::setEventSetup(), testbeamEBShape, testbeamEEShape, tokenGains_, tokenGrps_, tokenPeds_, tokenWgts_, and wgts_.

                                                                 {
   gains_ = es.getHandle(tokenGains_);
   peds_ = es.getHandle(tokenPeds_);
   grps_ = es.getHandle(tokenGrps_);
   wgts_ = es.getHandle(tokenWgts_);
 
   testbeamEEShape.setEventSetup(es);
   testbeamEBShape.setEventSetup(es);
 }