#include <EcalUncalibRecHitWorkerWeights.h>

Inheritance diagram for EcalUncalibRecHitWorkerWeights:

Public Member Functions
	EcalUncalibRecHitWorkerWeights ()

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

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 ()

	EcalUncalibRecHitWorkerRunOneDigiBase (const edm::ParameterSet &p)

	EcalUncalibRecHitWorkerRunOneDigiBase (const edm::ParameterSet &p, edm::ConsumesCollector &c)

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

	~EcalUncalibRecHitWorkerRunOneDigiBase () override

Public Member Functions inherited from EcalUncalibRecHitWorkerBaseClass
	EcalUncalibRecHitWorkerBaseClass ()

	EcalUncalibRecHitWorkerBaseClass (const edm::ParameterSet &)

	EcalUncalibRecHitWorkerBaseClass (const edm::ParameterSet &, edm::ConsumesCollector &c)

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

EcalUncalibRecHitRecWeightsAlgo< EBDataFrame >	uncalibMaker_barrel_

EcalUncalibRecHitRecWeightsAlgo< EEDataFrame >	uncalibMaker_endcap_

const EcalWeightSet::EcalWeightMatrix *	weights [2]

edm::ESHandle< EcalTBWeights >	wgts

Detailed Description

Definition at line 29 of file EcalUncalibRecHitWorkerWeights.h.

Constructor & Destructor Documentation

◆ EcalUncalibRecHitWorkerWeights() [1/2]

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

Definition at line 17 of file EcalUncalibRecHitWorkerWeights.cc.

18 : EcalUncalibRecHitWorkerRunOneDigiBase(ps, c), testbeamEEShape(EEShape(true)), testbeamEBShape(EBShape(true)) {}

◆ EcalUncalibRecHitWorkerWeights() [2/2]

EcalUncalibRecHitWorkerWeights::EcalUncalibRecHitWorkerWeights ( )

inline

Definition at line 32 of file EcalUncalibRecHitWorkerWeights.h.

32 : testbeamEEShape(EEShape(true)), testbeamEBShape(EBShape(true)) { ; }

◆ ~EcalUncalibRecHitWorkerWeights()

EcalUncalibRecHitWorkerWeights::~EcalUncalibRecHitWorkerWeights ( )

inlineoverride

Definition at line 33 of file EcalUncalibRecHitWorkerWeights.h.

33 {};

Member Function Documentation

◆ getAlgoDescription()

edm::ParameterSetDescription EcalUncalibRecHitWorkerWeights::getAlgoDescription ( )

overridevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 112 of file EcalUncalibRecHitWorkerWeights.cc.

                                                                               {
   edm::ParameterSetDescription psd;
   return psd;
 }

◆ run()

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

overridevirtual

Implements EcalUncalibRecHitWorkerRunOneDigiBase.

Definition at line 30 of file EcalUncalibRecHitWorkerWeights.cc.

                                                                                    {
   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;
 }

References EcalEndcap, EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), gainRatios, gains, EcalTBWeights::getMap(), EcalWeightSet::getWeightsAfterGainSwitch(), EcalWeightSet::getWeightsBeforeGainSwitch(), grps, ecalpyutils::hashedIndex(), EBDetId::hashedIndex(), EEDetId::hashedIndex(), EcalXtalGroupId::id(), EcalUncalibratedRecHit::kHasSwitchToGain1, EcalUncalibratedRecHit::kHasSwitchToGain6, EcalUncalibRecHitRecWeightsAlgo< C >::makeRecHit(), pedRMSVec, peds, pedVec, mps_fire::result, EcalUncalibratedRecHit::setFlagBit(), testbeamEBShape, testbeamEEShape, uncalibMaker_endcap_, weights, and wgts.

◆ set()

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

overridevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 20 of file EcalUncalibRecHitWorkerWeights.cc.

                                                                 {
   es.get<EcalGainRatiosRcd>().get(gains);
   es.get<EcalPedestalsRcd>().get(peds);
   es.get<EcalWeightXtalGroupsRcd>().get(grps);
   es.get<EcalTBWeightsRcd>().get(wgts);
  
   testbeamEEShape.setEventSetup(es);
   testbeamEBShape.setEventSetup(es);
 }