|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
|
#ifndef CalibCalorimetry_EcalTPGTools_EcalTPGParamBuilder_h
#define CalibCalorimetry_EcalTPGTools_EcalTPGParamBuilder_h
//Author: Pascal Paganini - LLR
//Date: 2006/07/10 15:58:06 $
// system include files
#include <memory>
// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/one/EDAnalyzer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Utilities/interface/ESGetToken.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/CaloTopology/interface/EcalTrigTowerConstituentsMap.h"
#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"
#include "Geometry/EcalMapping/interface/EcalMappingRcd.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"
#include "CondFormats/DataRecord/interface/EcalADCToGeVConstantRcd.h"
#include "CondFormats/DataRecord/interface/EcalIntercalibConstantsRcd.h"
#include "CondFormats/DataRecord/interface/EcalGainRatiosRcd.h"
#include "CondFormats/DataRecord/interface/EcalLaserAlphasRcd.h"
#include "CondFormats/DataRecord/interface/EcalLaserAPDPNRatiosRcd.h"
#include "CondFormats/DataRecord/interface/EcalPedestalsRcd.h"
#include "CondFormats/EcalObjects/interface/EcalADCToGeVConstant.h"
#include "CondFormats/EcalObjects/interface/EcalIntercalibConstants.h"
#include "CondFormats/EcalObjects/interface/EcalGainRatios.h"
#include "CondFormats/EcalObjects/interface/EcalLaserAPDPNRatios.h"
#include "CondFormats/EcalObjects/interface/EcalPedestals.h"
#include "OnlineDB/EcalCondDB/interface/all_monitoring_types.h"
#include "OnlineDB/EcalCondDB/interface/all_fe_config_types.h"
//modif-alex-27-july-2015
#include "CondFormats/EcalObjects/interface/EcalLaserAlphas.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EcalShapeBase.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EBShape.h"
#include "SimCalorimetry/EcalSimAlgos/interface/EEShape.h"
#include <TH1F.h>
#include <vector>
#include <string>
#include <map>
#include <iostream>
class EcalTPGDBApp;
class coeffStruc {
public:
coeffStruc() {}
double calibCoeff_;
double gainRatio_[3];
int pedestals_[3];
};
class linStruc {
public:
linStruc() {}
int pedestal_[3];
int mult_[3];
int shift_[3];
};
class EcalTPGParamBuilder : public edm::one::EDAnalyzer<> {
public:
explicit EcalTPGParamBuilder(edm::ParameterSet const& pSet);
~EcalTPGParamBuilder() override;
void analyze(const edm::Event& evt, const edm::EventSetup& evtSetup) override;
void beginJob() override;
private:
bool checkIfOK(EcalPedestals::Item item);
bool computeLinearizerParam(
double theta, double gainRatio, double calibCoeff, std::string subdet, int& mult, int& shift);
void create_header();
int uncodeWeight(double weight, int complement2 = 7);
double uncodeWeight(int iweight, int complement2 = 7);
std::vector<unsigned int> computeWeights(EcalShapeBase& shape, TH1F* histo);
void computeLUT(int* lut, std::string det = "EB");
//void getCoeff(coeffStruc & coeff, const EcalIntercalibConstantMap & calibMap, uint rawId) ; //modif-alex-27-july-2015 uncomment to go back
void getCoeff(coeffStruc& coeff, const EcalGainRatioMap& gainMap, uint rawId);
void getCoeff(coeffStruc& coeff, const EcalPedestalsMap& pedMap, uint rawId);
void getCoeff(coeffStruc& coeff, const std::map<EcalLogicID, MonPedestalsDat>& pedMap, const EcalLogicID& logicId);
//modif-alex-27-july-2015
void getCoeff(coeffStruc& coeff,
const EcalIntercalibConstantMap& calibMap,
const EcalLaserAlphaMap& laserAlphaMap,
uint rawId,
std::string& ss);
void computeFineGrainEBParameters(uint& lowRatio, uint& highRatio, uint& lowThreshold, uint& highThreshold, uint& lut);
void computeFineGrainEEParameters(uint& threshold, uint& lut_strip, uint& lut_tower);
int getEtaSlice(int tccId, int towerInTCC);
bool realignBaseline(linStruc& lin, float forceBase12, int cmsswid);
int getGCTRegionPhi(int ttphi);
int getGCTRegionEta(int tteta);
std::string getDet(int tcc);
std::pair<std::string, int> getCrate(int tcc);
edm::ESGetToken<CaloSubdetectorGeometry, EcalEndcapGeometryRecord> theEndcapGeometryToken_;
edm::ESGetToken<CaloSubdetectorGeometry, EcalBarrelGeometryRecord> theBarrelGeometryToken_;
const CaloSubdetectorGeometry* theEndcapGeometry_;
const CaloSubdetectorGeometry* theBarrelGeometry_;
edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> eTTmapToken_;
const EcalTrigTowerConstituentsMap* eTTmap_;
edm::ESGetToken<EcalElectronicsMapping, EcalMappingRcd> ecalmappingToken_;
const EcalElectronicsMapping* theMapping_;
edm::ESGetToken<EcalLaserAlphas, EcalLaserAlphasRcd> ecalLaserAlphasToken_;
edm::ESGetToken<EcalLaserAPDPNRatios, EcalLaserAPDPNRatiosRcd> ecalLaserAPDPNRatiosToken_;
edm::ESGetToken<EcalPedestals, EcalPedestalsRcd> ecalPedestalsToken_;
edm::ESGetToken<EcalIntercalibConstants, EcalIntercalibConstantsRcd> ecalIntercalibConstantsToken_;
edm::ESGetToken<EcalGainRatios, EcalGainRatiosRcd> ecalGainRatiosToken_;
edm::ESGetToken<EcalADCToGeVConstant, EcalADCToGeVConstantRcd> ecalADCToGeVConstantToken_;
bool useTransverseEnergy_;
double xtal_LSB_EB_, xtal_LSB_EE_;
double Et_sat_EB_, Et_sat_EE_;
unsigned int sliding_;
unsigned int sampleMax_;
double weight_timeShift_;
bool weight_unbias_recovery_;
bool weight_even_computeFromShape_;
std::string weight_even_idMapFile_;
std::string weight_even_weightGroupFile_;
std::string weight_odd_idMapFile_;
std::string weight_odd_weightGroupFile_;
static constexpr uint NWEIGHTS = 5; // It can be up to 6, but in Run 1,2,3 only 5 weights have been used
bool weight_useDoubleWeights_;
bool TPmode_EnableEBOddFilter_;
bool TPmode_EnableEEOddFilter_;
bool TPmode_EnableEBOddPeakFinder_;
bool TPmode_EnableEEOddPeakFinder_;
bool TPmode_DisableEBEvenPeakFinder_;
bool TPmode_DisableEEEvenPeakFinder_;
unsigned int TPmode_FenixEBStripOutput_;
unsigned int TPmode_FenixEEStripOutput_;
unsigned int TPmode_FenixEBStripInfobit2_;
unsigned int TPmode_FenixEEStripInfobit2_;
unsigned int TPmode_FenixEBTcpOutput_;
unsigned int TPmode_FenixEBTcpInfobit1_;
unsigned int TPmode_FenixEETcpOutput_;
unsigned int TPmode_FenixEETcpInfobit1_;
unsigned int nSample_;
unsigned int complement2_;
std::string LUT_option_;
double LUT_threshold_EB_, LUT_threshold_EE_;
double LUT_stochastic_EB_, LUT_noise_EB_, LUT_constant_EB_;
double LUT_stochastic_EE_, LUT_noise_EE_, LUT_constant_EE_;
double TTF_lowThreshold_EB_, TTF_highThreshold_EB_;
double TTF_lowThreshold_EE_, TTF_highThreshold_EE_;
double FG_lowThreshold_EB_, FG_highThreshold_EB_, FG_lowRatio_EB_, FG_highRatio_EB_;
unsigned int FG_lut_EB_;
double FG_Threshold_EE_;
unsigned int FG_lut_strip_EE_, FG_lut_tower_EE_;
int forcedPedestalValue_;
bool forceEtaSlice_;
unsigned int SFGVB_Threshold_, SFGVB_lut_, pedestal_offset_;
int SFGVB_SpikeKillingThreshold_; //modif-alex 01/21/11
bool useInterCalibration_, H2_;
//modif-alex-30/01/2012
std::string Transparency_Corr_;
bool useTransparencyCorr_;
//modif-alex-02/02/11
std::string TimingDelays_EB_;
std::string TimingDelays_EE_;
std::string TimingPhases_EB_;
std::string TimingPhases_EE_;
std::map<int, std::vector<int> > delays_EB_;
std::map<int, std::vector<int> > phases_EB_;
std::map<int, std::vector<int> > delays_EE_;
std::map<int, std::vector<int> > phases_EE_;
//modif-alex 30/01/2012
std::map<int, double> Transparency_Correction_;
EBShape shapeEB_;
EEShape shapeEE_;
std::ofstream* out_file_;
std::ofstream* geomFile_;
EcalTPGDBApp* db_;
bool writeToDB_;
bool writeToFiles_;
unsigned int DBrunNb_;
bool DBEE_;
int ped_conf_id_;
int lin_conf_id_;
int lut_conf_id_;
int fgr_conf_id_;
int sli_conf_id_;
int wei_conf_id_;
int wei2_conf_id_;
int spi_conf_id_; //modif-alex 21/01.11
int del_conf_id_; //modif-alex 21/01.11
int bxt_conf_id_;
int btt_conf_id_;
int bst_conf_id_;
int coke_conf_id_;
std::string tag_;
int version_;
int m_write_ped;
int m_write_lin;
int m_write_lut;
int m_write_wei;
int m_write_wei2;
int m_write_fgr;
int m_write_sli;
int m_write_spi; //modif-alex 21/01/11
int m_write_del; //modif-alex 21/01/11
int m_write_bxt;
int m_write_btt;
int m_write_bst;
Int_t* ntupleInts_;
Char_t ntupleDet_[10];
Char_t ntupleCrate_[10];
};
#endif
|