CMS 3D CMS Logo
EcalPFRecHitThresholdsMaker.cc
Go to the documentation of this file.
1 #include "FWCore/ServiceRegistry/interface/Service.h"
2 #include "CondCore/DBOutputService/interface/PoolDBOutputService.h"
3 
4 #include "FWCore/ParameterSet/interface/ParameterSet.h"
5 
6 
7 #include "FWCore/Framework/interface/ESHandle.h"
8 #include "FWCore/Framework/interface/EventSetup.h"
9 
10 #include "FWCore/Framework/interface/Event.h"
11 
12 #include "CalibCalorimetry/EcalLaserCorrection/interface/EcalLaserDbService.h"
13 #include "CalibCalorimetry/EcalLaserCorrection/interface/EcalLaserDbRecord.h"
14 
15 #include "CondTools/Ecal/interface/EcalPFRecHitThresholdsMaker.h"
16 #include "CondFormats/EcalObjects/interface/EcalPedestals.h"
17 #include "CondFormats/DataRecord/interface/EcalPedestalsRcd.h"
18 #include "CondFormats/EcalObjects/interface/EcalADCToGeVConstant.h"
19 #include "CondFormats/DataRecord/interface/EcalADCToGeVConstantRcd.h"
20 #include "CondFormats/EcalObjects/interface/EcalTimeCalibConstants.h"
21 #include "CondFormats/DataRecord/interface/EcalTimeCalibConstantsRcd.h"
22 #include "CondFormats/EcalObjects/interface/EcalChannelStatus.h"
23 #include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"
24 #include "CondFormats/EcalObjects/interface/EcalLinearCorrections.h"
25 #include "CondFormats/DataRecord/interface/EcalLinearCorrectionsRcd.h"
26 #include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
27 #include "CondFormats/DataRecord/interface/EcalIntercalibConstantsRcd.h"
28 #include "CondFormats/EcalObjects/interface/EcalIntercalibConstantsMC.h"
29 #include "CondFormats/DataRecord/interface/EcalIntercalibConstantsMCRcd.h"
30 #include "CondFormats/EcalObjects/interface/EcalIntercalibErrors.h"
31 #include "CondFormats/DataRecord/interface/EcalIntercalibErrorsRcd.h"
32 #include "CondFormats/EcalObjects/interface/EcalGainRatios.h"
33 #include "CondFormats/DataRecord/interface/EcalGainRatiosRcd.h"
34 #include "CondFormats/EcalObjects/interface/EcalWeightXtalGroups.h"
35 #include "CondFormats/DataRecord/interface/EcalWeightXtalGroupsRcd.h"
36 #include "CondFormats/EcalObjects/interface/EcalTBWeights.h"
37 #include "CondFormats/DataRecord/interface/EcalTBWeightsRcd.h"
38 #include "CondFormats/EcalObjects/interface/EcalLaserAlphas.h"
39 #include "CondFormats/DataRecord/interface/EcalLaserAlphasRcd.h"
40 #include "CondFormats/EcalObjects/interface/EcalLaserAPDPNRatios.h"
41 #include "CondFormats/DataRecord/interface/EcalLaserAPDPNRatiosRcd.h"
42 #include "CondFormats/EcalObjects/interface/EcalLaserAPDPNRatiosRef.h"
43 #include "CondFormats/DataRecord/interface/EcalLaserAPDPNRatiosRefRcd.h"
44 #include "CondFormats/EcalObjects/interface/EcalTPGCrystalStatus.h"
45 #include "CondFormats/DataRecord/interface/EcalTPGCrystalStatusRcd.h"
46 #include "CondFormats/EcalObjects/interface/EcalTPGTowerStatus.h"
47 #include "CondFormats/DataRecord/interface/EcalTPGTowerStatusRcd.h"
48 
49 #include "CondFormats/EcalObjects/interface/EcalDCSTowerStatus.h"
50 #include "CondFormats/DataRecord/interface/EcalDCSTowerStatusRcd.h"
51 #include "CondFormats/EcalObjects/interface/EcalDAQTowerStatus.h"
52 #include "CondFormats/DataRecord/interface/EcalDAQTowerStatusRcd.h"
53 #include "CondFormats/EcalObjects/interface/EcalDQMTowerStatus.h"
54 #include "CondFormats/DataRecord/interface/EcalDQMTowerStatusRcd.h"
55 #include "CondFormats/EcalObjects/interface/EcalDQMChannelStatus.h"
56 #include "CondFormats/DataRecord/interface/EcalDQMChannelStatusRcd.h"
57 
58 #include "CondFormats/EcalObjects/interface/EcalClusterCrackCorrParameters.h"
59 #include "CondFormats/DataRecord/interface/EcalClusterCrackCorrParametersRcd.h"
60 
61 #include "CondFormats/EcalObjects/interface/EcalPFRecHitThresholds.h"
62 #include "CondFormats/DataRecord/interface/EcalPFRecHitThresholdsRcd.h"
63 
64 #include "CondFormats/EcalObjects/interface/EcalClusterEnergyCorrectionParameters.h"
65 #include "CondFormats/DataRecord/interface/EcalClusterEnergyCorrectionParametersRcd.h"
66 #include "CondFormats/EcalObjects/interface/EcalClusterEnergyUncertaintyParameters.h"
67 #include "CondFormats/DataRecord/interface/EcalClusterEnergyUncertaintyParametersRcd.h"
68 #include "CondFormats/EcalObjects/interface/EcalClusterEnergyCorrectionObjectSpecificParameters.h"
69 #include "CondFormats/DataRecord/interface/EcalClusterEnergyCorrectionObjectSpecificParametersRcd.h"
70 #include "CondFormats/EcalObjects/interface/EcalClusterLocalContCorrParameters.h"
71 #include "CondFormats/DataRecord/interface/EcalClusterLocalContCorrParametersRcd.h"
72 
73 #include "CondFormats/Alignment/interface/Alignments.h"
74 #include "CondFormats/AlignmentRecord/interface/EBAlignmentRcd.h"
75 #include "CondFormats/AlignmentRecord/interface/EEAlignmentRcd.h"
76 #include "CondFormats/AlignmentRecord/interface/ESAlignmentRcd.h"
77 #include "CondFormats/EcalObjects/interface/EcalTimeOffsetConstant.h"
78 #include "CondFormats/DataRecord/interface/EcalTimeOffsetConstantRcd.h"
79 
80 #include "CondFormats/EcalObjects/interface/EcalSampleMask.h"
81 #include "CondFormats/DataRecord/interface/EcalSampleMaskRcd.h"
82 
83 #include "CondFormats/EcalObjects/interface/EcalTimeBiasCorrections.h"
84 #include "CondFormats/DataRecord/interface/EcalTimeBiasCorrectionsRcd.h"
85 
86 #include "CondFormats/EcalObjects/interface/EcalSamplesCorrelation.h"
87 #include "CondFormats/DataRecord/interface/EcalSamplesCorrelationRcd.h"
88 
89 #include "DataFormats/Provenance/interface/Timestamp.h"
90 
91 #include "CondCore/DBOutputService/interface/PoolDBOutputService.h"
92 
93 #include <vector>
94 
95 EcalPFRecHitThresholdsMaker::EcalPFRecHitThresholdsMaker(const edm::ParameterSet& iConfig) :
96  m_timetype(iConfig.getParameter<std::string>("timetype"))
97 {
98 
99  std::string container;
100  std::string tag;
101  std::string record;
102 
103  m_nsigma = iConfig.getParameter<double>("NSigma");
104 
105 }
106 
107 
108 EcalPFRecHitThresholdsMaker::~EcalPFRecHitThresholdsMaker()
109 {
110 
111 }
112 
113 void EcalPFRecHitThresholdsMaker::analyze( const edm::Event& evt, const edm::EventSetup& evtSetup){
114 
115  edm::Service<cond::service::PoolDBOutputService> dbOutput;
116  if ( !dbOutput.isAvailable() ) {
117  throw cms::Exception("PoolDBOutputService is not available");
118  }
119 
120 
121 
122  edm::ESHandle<EcalPedestals> handle1;
123  evtSetup.get<EcalPedestalsRcd>().get(handle1);
124  const EcalPedestals* ped_db = handle1.product();
125  std::cout << "ped pointer is: "<< ped_db<< std::endl;
126 
127 
128  edm::ESHandle<EcalADCToGeVConstant> handle2;
129  evtSetup.get<EcalADCToGeVConstantRcd>().get(handle2);
130  const EcalADCToGeVConstant* adc_db = handle2.product();
131  std::cout << "adc pointer is: "<< adc_db<< std::endl;
132 
133  edm::ESHandle<EcalIntercalibConstants> handle3;
134  evtSetup.get<EcalIntercalibConstantsRcd>().get(handle3);
135  const EcalIntercalibConstants* ical_db = handle3.product();
136  std::cout << "inter pointer is: "<< ical_db<< std::endl;
137 
138 
139  edm::ESHandle<EcalLaserDbService> laser;
140  evtSetup.get<EcalLaserDbRecord>().get(laser);
141 
142  EcalPFRecHitThresholds* pfthresh = new EcalPFRecHitThresholds();
143 
144  // const EcalIntercalibConstantMap& icalMap = ical_db->getMap();
145 
146 
147  float adc_EB= float(adc_db->getEEValue()) ;
148  float adc_EE=float(adc_db->getEBValue());
149 
150 
151  //edm::Timestamp tsince;
152 
153 
154 
155  for(int iEta=-EBDetId::MAX_IETA; iEta<=EBDetId::MAX_IETA ;++iEta) {
156  if(iEta==0) continue;
157  for(int iPhi=EBDetId::MIN_IPHI; iPhi<=EBDetId::MAX_IPHI; ++iPhi) {
158  // make an EBDetId since we need EBDetId::rawId() to be used as the key for the pedestals
159  if (EBDetId::validDetId(iEta,iPhi)) {
160  EBDetId ebdetid(iEta,iPhi,EBDetId::ETAPHIMODE);
161  EcalPedestals::const_iterator it =ped_db->find(ebdetid.rawId());
162  EcalPedestals::Item aped = (*it);
163 
164  EcalIntercalibConstants::const_iterator itc =ical_db->find(ebdetid.rawId());
165  float calib = (*itc);
166 
167  // get laser coefficient
168  float lasercalib = 1.;
169  lasercalib = laser->getLaserCorrection( ebdetid, evt.time()); // TODO correct time
170 
171  EcalPFRecHitThreshold thresh= aped.rms_x12 * calib * adc_EB * lasercalib * m_nsigma;
172 
173  if(iPhi==100) std::cout<<"Thresh(GeV)="<<thresh<<std::endl;
174 
175  pfthresh->insert(std::make_pair(ebdetid.rawId(),thresh));
176  }
177  }
178  }
179 
180  for(int iX=EEDetId::IX_MIN; iX<=EEDetId::IX_MAX ;++iX) {
181  for(int iY=EEDetId::IY_MIN; iY<=EEDetId::IY_MAX; ++iY) {
182  // make an EEDetId since we need EEDetId::rawId() to be used as the key for the pedestals
183  if (EEDetId::validDetId(iX,iY,1)) {
184  EEDetId eedetid(iX,iY,1);
185 
186  EcalPedestals::const_iterator it =ped_db->find(eedetid.rawId());
187  EcalPedestals::Item aped = (*it);
188 
189  EcalIntercalibConstants::const_iterator itc =ical_db->find(eedetid.rawId());
190  float calib = (*itc);
191 
192  // get laser coefficient
193  float lasercalib = 1.;
194  lasercalib = laser->getLaserCorrection( eedetid, evt.time()); // TODO correct time
195 
196  EcalPFRecHitThreshold thresh= aped.rms_x12 * calib * adc_EE * lasercalib * m_nsigma;
197  pfthresh->insert(std::make_pair(eedetid.rawId(),thresh));
198  }
199  if(EEDetId::validDetId(iX,iY,-1)) {
200 
201  EEDetId eedetid(iX,iY,-1);
202 
203  EcalPedestals::const_iterator it =ped_db->find(eedetid.rawId());
204  EcalPedestals::Item aped = (*it);
205 
206  EcalIntercalibConstants::const_iterator itc =ical_db->find(eedetid.rawId());
207  float calib = (*itc);
208 
209  // get laser coefficient
210  float lasercalib = 1.;
211  lasercalib = laser->getLaserCorrection( eedetid, evt.time()); // TODO correct time
212 
213  EcalPFRecHitThreshold thresh= aped.rms_x12 * calib * adc_EE * lasercalib * m_nsigma;
214  pfthresh->insert(std::make_pair(eedetid.rawId(),thresh));
215 
216  if(iX==50) std::cout<<"Thresh(GeV)="<<thresh<<std::endl;
217 
218 
219  }
220 
221  }
222  }
223 
224 
225 
226 
227  dbOutput->createNewIOV<const EcalPFRecHitThresholds>( pfthresh , dbOutput->beginOfTime(), dbOutput->endOfTime(),"EcalPFRecHitThresholdsRcd");
228 
229 
230 
231  std::cout<< "EcalPFRecHitThresholdsMaker wrote it " << std::endl;
232 }
233 
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
EBDetId::MIN_IPHI
static const int MIN_IPHI
Definition: EBDetId.h:142
EcalPFRecHitThresholds
EcalPFRecHitThresholdsMap EcalPFRecHitThresholds
Definition: EcalPFRecHitThresholds.h:12
GlobalPosition_Frontier_DevDB_cff.tag
tag
Definition: GlobalPosition_Frontier_DevDB_cff.py:11
EcalPFRecHitThresholdsMaker::EcalPFRecHitThresholdsMaker
EcalPFRecHitThresholdsMaker(const edm::ParameterSet &iConfig)
Definition: EcalPFRecHitThresholdsMaker.cc:95
record
JetCorrectorParameters::Record record
Definition: classes.h:7
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
EcalLaserDbService::getLaserCorrection
float getLaserCorrection(DetId const &xid, edm::Timestamp const &iTime) const
Definition: EcalLaserDbService.cc:44
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:47
scrapingFilter_cfi.thresh
thresh
Definition: scrapingFilter_cfi.py:6
EEDetId::IX_MIN
static const int IX_MIN
Definition: EEDetId.h:294
EEDetId::IY_MIN
static const int IY_MIN
Definition: EEDetId.h:298
EBDetId::validDetId
static bool validDetId(int i, int j)
check if a valid index combination
Definition: EBDetId.h:124
EcalPFRecHitThresholdsMaker::analyze
void analyze(const edm::Event &evt, const edm::EventSetup &evtSetup) override
Definition: EcalPFRecHitThresholdsMaker.cc:113
EcalPFRecHitThreshold
float EcalPFRecHitThreshold
Definition: EcalPFRecHitThresholds.h:10
edm::Service::isAvailable
bool isAvailable() const
Definition: Service.h:46
EBDetId::ETAPHIMODE
static const int ETAPHIMODE
Definition: EBDetId.h:166
cond::service::PoolDBOutputService::createNewIOV
void createNewIOV(T *firstPayloadObj, cond::Time_t firstSinceTime, cond::Time_t firstTillTime, const std::string &recordName, bool withlogging=false)
Definition: PoolDBOutputService.h:73
EcalCondObjectContainer::insert
void insert(std::pair< uint32_t, Item > const &a)
Definition: EcalCondObjectContainer.h:35
EEDetId::IX_MAX
static const int IX_MAX
Definition: EEDetId.h:302
EBDetId::MAX_IPHI
static const int MAX_IPHI
Definition: EBDetId.h:144
EEDetId::validDetId
static bool validDetId(int crystal_ix, int crystal_iy, int iz)
Definition: EEDetId.h:248
EcalCondObjectContainer< EcalPedestal >::const_iterator
std::vector< Item >::const_iterator const_iterator
Definition: EcalCondObjectContainer.h:17
EBDetId::MAX_IETA
static const int MAX_IETA
Definition: EBDetId.h:143
edm::EventSetup::get
T get() const
Definition: EventSetup.h:63
EcalCondObjectContainer::find
const_iterator find(uint32_t rawId) const
Definition: EcalCondObjectContainer.h:54
EEDetId::IY_MAX
static const int IY_MAX
Definition: EEDetId.h:306
gather_cfg.cout
cout
Definition: gather_cfg.py:143
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
edm::EventBase::time
edm::Timestamp time() const
Definition: EventBase.h:61
Generated for CMSSW Reference Manual by doxygen 1.8.11