CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
EcalUncalibRecHitWorkerAnalFit Class Reference

#include <EcalUncalibRecHitWorkerAnalFit.h>

Inheritance diagram for EcalUncalibRecHitWorkerAnalFit:
EcalUncalibRecHitWorkerRunOneDigiBase EcalUncalibRecHitWorkerBaseClass

Public Member Functions

 EcalUncalibRecHitWorkerAnalFit (const edm::ParameterSet &ps, edm::ConsumesCollector &c)
 
 EcalUncalibRecHitWorkerAnalFit ()
 
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
 
 ~EcalUncalibRecHitWorkerAnalFit () 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 ()
 

Private Attributes

EcalUncalibRecHitRecAnalFitAlgo< EBDataFramealgoEB_
 
EcalUncalibRecHitRecAnalFitAlgo< EEDataFramealgoEE_
 
edm::ESHandle< EcalPedestalspedHandle
 
edm::ESHandle< EcalGainRatiospRatio
 

Detailed Description

produce ECAL uncalibrated rechits from dataframes with the analytical fit method

Author
Shahram Rahatlou, University of Rome & INFN, Sept 2005

Definition at line 27 of file EcalUncalibRecHitWorkerAnalFit.h.

Constructor & Destructor Documentation

EcalUncalibRecHitWorkerAnalFit::EcalUncalibRecHitWorkerAnalFit ( const edm::ParameterSet ps,
edm::ConsumesCollector c 
)
EcalUncalibRecHitWorkerAnalFit::EcalUncalibRecHitWorkerAnalFit ( )
inline

Definition at line 30 of file EcalUncalibRecHitWorkerAnalFit.h.

30 {};
EcalUncalibRecHitWorkerAnalFit::~EcalUncalibRecHitWorkerAnalFit ( )
inlineoverride

Definition at line 31 of file EcalUncalibRecHitWorkerAnalFit.h.

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

31 {};

Member Function Documentation

edm::ParameterSetDescription EcalUncalibRecHitWorkerAnalFit::getAlgoDescription ( )
overridevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 131 of file EcalUncalibRecHitWorkerAnalFit.cc.

References DEFINE_EDM_PLUGIN.

131  {
133  //psd.addNode(edm::ParameterSet<edm::InputTag>("EBdigiCollection", edm::InputTag("ecalDigis","ebDigis")) and
134  // edm::ParameterSet<std::string>("EEhitCollection", "EcalUncalibRecHitsEE") and
135  // edm::ParameterSet<edm::InputTag>("EEdigiCollection", edm::InputTag("ecalDigis","eeDigis")) and
136  // edm::ParameterSet<std::string>("algo", "EcalUncalibRecHitWorkerAnalFit") and
137  // edm::ParameterSet<std::string>("EBhitCollection", "EcalUncalibRecHitsEB"));
138  //
139  return psd;
140 }
bool EcalUncalibRecHitWorkerAnalFit::run ( const edm::Event evt,
const EcalDigiCollection::const_iterator digi,
EcalUncalibratedRecHitCollection result 
)
overridevirtual

Implements EcalUncalibRecHitWorkerRunOneDigiBase.

Definition at line 55 of file EcalUncalibRecHitWorkerAnalFit.cc.

References algoEB_, algoEE_, EcalUncalibratedRecHit::amplitude(), EcalBarrel, EcalCondObjectContainer< T >::end(), EcalCondObjectContainer< T >::find(), EcalMGPAGainRatio::gain12Over6(), EcalMGPAGainRatio::gain6Over1(), LogDebug, EcalUncalibRecHitRecAnalFitAlgo< C >::makeRecHit(), EcalPedestal::mean_x1, EcalPedestal::mean_x12, EcalPedestal::mean_x6, pedHandle, pRatio, and edm::SortedCollection< T, SORT >::push_back().

57  {
58  using namespace edm;
59 
60  const EcalGainRatioMap& gainMap = pRatio.product()->getMap(); // map of gain ratios
61  const EcalPedestalsMap& pedMap = pedHandle.product()->getMap(); // map of pedestals
62 
63  EcalPedestalsMapIterator pedIter; // pedestal iterator
64  EcalPedestals::Item aped; // pedestal object for a single xtal
65 
66  EcalGainRatioMap::const_iterator gainIter; // gain iterator
67  EcalMGPAGainRatio aGain; // gain object for a single xtal
68 
69  DetId detid(itdg->id());
70 
71  // find pedestals for this channel
72  //LogDebug("EcalUncalibRecHitDebug") << "looking up pedestal for crystal: " << itdg->id(); // FIXME
73  pedIter = pedMap.find(detid);
74  if (pedIter != pedMap.end()) {
75  aped = (*pedIter);
76  } else {
77  edm::LogError("EcalUncalibRecHitWorkerAnalFit") << "error!! could not find pedestals for channel: ";
78  if (detid.subdetId() == EcalBarrel) {
79  edm::LogError("EcalUncalibRecHitWorkerAnalFit") << EBDetId(detid);
80  } else {
81  edm::LogError("EcalUncalibRecHitWorkerAnalFit") << EEDetId(detid);
82  }
83  edm::LogError("EcalUncalibRecHitWorkerAnalFit") << "\n no uncalib rechit will be made for this digi!";
84  return false;
85  }
86  double pedVec[3];
87  pedVec[0] = aped.mean_x12;
88  pedVec[1] = aped.mean_x6;
89  pedVec[2] = aped.mean_x1;
90 
91  // find gain ratios
92  //LogDebug("EcalUncalibRecHitDebug") << "looking up gainRatios for crystal: " << itdg->id(); // FIXME
93  gainIter = gainMap.find(detid);
94  if (gainIter != gainMap.end()) {
95  aGain = (*gainIter);
96  } else {
97  edm::LogError("EcalUncalibRecHitWorkerAnalFit") << "error!! could not find gain ratios for channel: ";
98  if (detid.subdetId() == EcalBarrel) {
99  edm::LogError("EcalUncalibRecHitWorkerAnalFit") << EBDetId(detid);
100  } else {
101  edm::LogError("EcalUncalibRecHitWorkerAnalFit") << EEDetId(detid);
102  }
103  edm::LogError("EcalUncalibRecHitWorkerAnalFit") << "\n no uncalib rechit will be made for this digi!";
104  return false;
105  }
106  double gainRatios[3];
107  gainRatios[0] = 1.;
108  gainRatios[1] = aGain.gain12Over6();
109  gainRatios[2] = aGain.gain6Over1() * aGain.gain12Over6();
110 
111  if (detid.subdetId() == EcalBarrel) {
112  EcalUncalibratedRecHit aHit = algoEB_.makeRecHit(*itdg, pedVec, gainRatios, nullptr, nullptr);
113  result.push_back(aHit);
114  if (aHit.amplitude() > 0.) {
115  LogDebug("EcalUncalibRecHitInfo") << "EcalUncalibRecHitWorkerAnalFit: processed EBDataFrame with id: "
116  << EBDetId(detid) << "\n"
117  << "uncalib rechit amplitude: " << aHit.amplitude();
118  }
119  } else {
120  EcalUncalibratedRecHit aHit = algoEE_.makeRecHit(*itdg, pedVec, gainRatios, nullptr, nullptr);
121  result.push_back(aHit);
122  if (aHit.amplitude() > 0.) {
123  LogDebug("EcalUncalibRecHitInfo") << "EcalUncalibRecHitWorkerAnalFit: processed EEDataFrame with id: "
124  << EEDetId(detid) << "\n"
125  << "uncalib rechit amplitude: " << aHit.amplitude();
126  }
127  }
128  return true;
129 }
#define LogDebug(id)
void push_back(T const &t)
edm::ESHandle< EcalPedestals > pedHandle
EcalPedestalsMap::const_iterator EcalPedestalsMapIterator
Definition: EcalPedestals.h:49
EcalUncalibRecHitRecAnalFitAlgo< EBDataFrame > algoEB_
EcalUncalibRecHitRecAnalFitAlgo< EEDataFrame > algoEE_
float gain6Over1() const
Definition: DetId.h:17
std::vector< Item >::const_iterator const_iterator
EcalUncalibratedRecHit makeRecHit(const C &dataFrame, const double *pedestals, const double *gainRatios, const EcalWeightSet::EcalWeightMatrix **weights, const EcalWeightSet::EcalChi2WeightMatrix **chi2Matrix) override
Compute parameters.
float gain12Over6() const
edm::ESHandle< EcalGainRatios > pRatio
HLT enums.
const_iterator find(uint32_t rawId) const
const_iterator end() const
void EcalUncalibRecHitWorkerAnalFit::set ( const edm::EventSetup es)
overridevirtual

Implements EcalUncalibRecHitWorkerBaseClass.

Definition at line 43 of file EcalUncalibRecHitWorkerAnalFit.cc.

References edm::EventSetup::get(), LogDebug, pedHandle, and pRatio.

43  {
44  // Gain Ratios
45  LogDebug("EcalUncalibRecHitDebug") << "fetching gainRatios....";
46  es.get<EcalGainRatiosRcd>().get(pRatio);
47  LogDebug("EcalUncalibRecHitDebug") << "done.";
48 
49  // fetch the pedestals from the cond DB via EventSetup
50  LogDebug("EcalUncalibRecHitDebug") << "fetching pedestals....";
51  es.get<EcalPedestalsRcd>().get(pedHandle);
52  LogDebug("EcalUncalibRecHitDebug") << "done.";
53 }
#define LogDebug(id)
edm::ESHandle< EcalPedestals > pedHandle
edm::ESHandle< EcalGainRatios > pRatio
T get() const
Definition: EventSetup.h:73

Member Data Documentation

EcalUncalibRecHitRecAnalFitAlgo<EBDataFrame> EcalUncalibRecHitWorkerAnalFit::algoEB_
private

Definition at line 41 of file EcalUncalibRecHitWorkerAnalFit.h.

Referenced by run().

EcalUncalibRecHitRecAnalFitAlgo<EEDataFrame> EcalUncalibRecHitWorkerAnalFit::algoEE_
private

Definition at line 42 of file EcalUncalibRecHitWorkerAnalFit.h.

Referenced by run().

edm::ESHandle<EcalPedestals> EcalUncalibRecHitWorkerAnalFit::pedHandle
private

Definition at line 45 of file EcalUncalibRecHitWorkerAnalFit.h.

Referenced by run(), and set().

edm::ESHandle<EcalGainRatios> EcalUncalibRecHitWorkerAnalFit::pRatio
private

Definition at line 44 of file EcalUncalibRecHitWorkerAnalFit.h.

Referenced by run(), and set().