CMS 3D CMS Logo
SiStripGainFromData.cc
Go to the documentation of this file.
1 // Original Author: Loic QUERTENMONT
2 // Created: Wed Feb 6 08:55:18 CET 2008
3 
4 #include <memory>
5 
6 #include "FWCore/Framework/interface/Frameworkfwd.h"
7 #include "FWCore/Framework/interface/EDAnalyzer.h"
8 #include "FWCore/Framework/interface/Event.h"
9 #include "FWCore/Framework/interface/MakerMacros.h"
10 #include "FWCore/Framework/interface/ESHandle.h"
11 #include "FWCore/ParameterSet/interface/ParameterSet.h"
12 #include "FWCore/MessageLogger/interface/MessageLogger.h"
13 #include "FWCore/Utilities/interface/Exception.h"
14 
15 #include "Geometry/CommonDetUnit/interface/GeomDet.h"
16 #include "Geometry/CommonDetUnit/interface/GeomDetType.h"
17 #include "Geometry/CommonTopologies/interface/StripTopology.h"
18 #include "DataFormats/GeometrySurface/interface/TrapezoidalPlaneBounds.h"
19 #include "DataFormats/GeometrySurface/interface/RectangularPlaneBounds.h"
20 
21 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
22 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
23 #include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
24 #include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
25 #include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"
26 #include "Geometry/CommonDetUnit/interface/TrackingGeometry.h"
27 
28 #include "DataFormats/SiStripCluster/interface/SiStripClusterCollection.h"
29 
30 #include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
31 #include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
32 
33 #include "DataFormats/FEDRawData/interface/FEDNumbering.h"
34 #include "DataFormats/TrackReco/interface/Track.h"
35 #include "DataFormats/TrackReco/interface/TrackFwd.h"
36 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit1D.h"
37 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2D.h"
38 #include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2D.h"
39 #include "DataFormats/DetId/interface/DetId.h"
40 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
41 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
42 #include "Geometry/Records/interface/TrackerTopologyRcd.h"
43 #include "DataFormats/TrackReco/interface/DeDxHit.h"
44 #include "DataFormats/TrackReco/interface/TrackDeDxHits.h"
45 
46 #include "CommonTools/ConditionDBWriter/interface/ConditionDBWriter.h"
47 #include "CondFormats/SiStripObjects/interface/SiStripApvGain.h"
48 
49 #include "TrackingTools/PatternTools/interface/Trajectory.h"
50 #include "TrackingTools/PatternTools/interface/TrajTrackAssociation.h"
51 
52 #include "DQMServices/Core/interface/DQMStore.h"
53 #include "DQMServices/Core/interface/MonitorElement.h"
54 
55 #include "CalibFormats/SiStripObjects/interface/SiStripGain.h"
56 #include "CalibTracker/Records/interface/SiStripGainRcd.h"
57 
58 
59 #include "TFile.h"
60 #include "TObjString.h"
61 #include "TString.h"
62 #include "TH1F.h"
63 #include "TH2F.h"
64 #include "TProfile.h"
65 #include "TF1.h"
66 #include "TROOT.h"
67 
68 #include <unordered_map>
69 
70 
71 
72 using namespace edm;
73 using namespace reco;
74 using namespace std;
75 
76 
77 struct stAPVGain{unsigned int Index; int DetId; int APVId; int SubDet; float Eta; float R; float Phi; float Thickness; double MPV; double Gain; double PreviousGain; char Side;};
78 
79 class SiStripGainFromData : public ConditionDBWriter<SiStripApvGain> {
80  public:
81  explicit SiStripGainFromData(const edm::ParameterSet&);
82  ~SiStripGainFromData() override;
83 
84 
85  private:
86  void algoBeginJob(const edm::EventSetup&) override ;
87  void algoEndJob() override ;
88  void algoBeginRun(const edm::Run &, const edm::EventSetup &) override;
89 // virtual void algoBeginRun(const edm::Event& iEvent, const edm::EventSetup& iSetup);
90  void algoAnalyze(const edm::Event &, const edm::EventSetup &) override;
91 
92  SiStripApvGain* getNewObject() override;
93  DQMStore* dqmStore_;
94  DQMStore* dqmStore_infile;
95 
96  double ComputeChargeOverPath(const SiStripCluster* Cluster,TrajectoryStateOnSurface trajState, const edm::EventSetup* iSetup, const Track* track, double trajChi2OverN);
97  bool IsFarFromBorder(TrajectoryStateOnSurface trajState, const uint32_t detid, const edm::EventSetup* iSetup);
98 
99  void getPeakOfLandau(TH1* InputHisto, double* FitResults, double LowRange=0, double HighRange=5400);
100 
101  const edm::EventSetup* iSetup_;
102  const edm::Event* iEvent_;
103 
104  bool CheckLocalAngle;
105  unsigned int MinNrEntries;
106  double MaxMPVError;
107  double MaxChi2OverNDF;
108  double MinTrackMomentum;
109  double MaxTrackMomentum;
110  double MinTrackEta;
111  double MaxTrackEta;
112  unsigned int MaxNrStrips;
113  unsigned int MinTrackHits;
114  double MaxTrackChiOverNdf;
115  bool AllowSaturation;
116  bool FirstSetOfConstants;
117  bool Validation;
118  int CalibrationLevel;
119  bool CheckIfFileExist;
120 
121  std::string AlgoMode;
122  std::string OutputGains;
123  std::string OutputHistos;
124  std::string TrajToTrackProducer;
125  std::string TrajToTrackLabel;
126 
127  vector<string> VInputFiles;
128 
129  MonitorElement* tmp;
130 
131  TH2F* Tracks_P_Vs_Eta;
132  TH2F* Tracks_Pt_Vs_Eta;
133 
134  TH1F* NumberOfEntriesByAPV;
135  TH1F* HChi2OverNDF;
136  TH1F* HTrackChi2OverNDF;
137  TH1F* HTrackHits;
138 
139  TH2F* MPV_Vs_EtaTIB;
140  TH2F* MPV_Vs_EtaTID;
141  TH2F* MPV_Vs_EtaTOB;
142  TH2F* MPV_Vs_EtaTEC;
143  TH2F* MPV_Vs_EtaTEC1;
144  TH2F* MPV_Vs_EtaTEC2;
145 
146  TH2F* MPV_Vs_PhiTIB;
147  TH2F* MPV_Vs_PhiTID;
148  TH2F* MPV_Vs_PhiTOB;
149  TH2F* MPV_Vs_PhiTEC;
150  TH2F* MPV_Vs_PhiTEC1;
151  TH2F* MPV_Vs_PhiTEC2;
152 
153  TH2F* Charge_Vs_PathTIB;
154  TH2F* Charge_Vs_PathTID;
155  TH2F* Charge_Vs_PathTOB;
156  TH2F* Charge_Vs_PathTEC;
157  TH2F* Charge_Vs_PathTEC1;
158  TH2F* Charge_Vs_PathTEC2;
159 
160  TH1F* MPV_Vs_PathTIB;
161  TH1F* MPV_Vs_PathTID;
162  TH1F* MPV_Vs_PathTOB;
163  TH1F* MPV_Vs_PathTEC;
164  TH1F* MPV_Vs_PathTEC1;
165  TH1F* MPV_Vs_PathTEC2;
166 
167 
168  TH1F* Charge_TIB;
169  TH1F* Charge_TID;
170  TH1F* Charge_TIDP;
171  TH1F* Charge_TIDM;
172  TH1F* Charge_TOB;
173  TH1F* Charge_TEC;
174  TH1F* Charge_TEC1;
175  TH1F* Charge_TEC2;
176  TH1F* Charge_TECP;
177  TH1F* Charge_TECM;
178 
179  TH2F* MPV_Vs_Phi;
180  TH2F* MPV_Vs_Eta;
181  TH2F* MPV_Vs_R;
182 
183 // TH2F* PD_Vs_Eta;
184 // TH2F* PD_Vs_R;
185 
186  TH1F* APV_DetId;
187  TH1F* APV_Id;
188  TH1F* APV_Eta;
189  TH1F* APV_R;
190  TH1F* APV_SubDet;
191  TH2F* APV_Momentum;
192  TH2F* APV_Charge;
193  TH2F* APV_PathLength;
194  TH1F* APV_PathLengthM;
195  TH1F* APV_MPV;
196  TH1F* APV_Gain;
197  TH1F* APV_CumulGain;
198  TH1F* APV_PrevGain;
199  TH1F* APV_Thickness;
200 
201  TH1F* MPVs;
202  TH1F* MPVs320;
203  TH1F* MPVs500;
204 
205 // TH2F* MPV_vs_10RplusEta;
206 
207 
208  TH1F* NSatStripInCluster;
209  TH1F* NHighStripInCluster;
210 // TH2F* Charge_Vs_PathLength_CS1;
211 // TH2F* Charge_Vs_PathLength_CS2;
212 // TH2F* Charge_Vs_PathLength_CS3;
213 // TH2F* Charge_Vs_PathLength_CS4;
214 // TH2F* Charge_Vs_PathLength_CS5;
215 
216 // TH1F* MPV_Vs_PathLength_CS1;
217 // TH1F* MPV_Vs_PathLength_CS2;
218 // TH1F* MPV_Vs_PathLength_CS3;
219 // TH1F* MPV_Vs_PathLength_CS4;
220 // TH1F* MPV_Vs_PathLength_CS5;
221 
222 // TH1F* FWHM_Vs_PathLength_CS1;
223 // TH1F* FWHM_Vs_PathLength_CS2;
224 // TH1F* FWHM_Vs_PathLength_CS3;
225 // TH1F* FWHM_Vs_PathLength_CS4;
226 // TH1F* FWHM_Vs_PathLength_CS5;
227 
228 
229  TH2F* Charge_Vs_PathLength;
230  TH2F* Charge_Vs_PathLength320;
231  TH2F* Charge_Vs_PathLength500;
232 
233  TH1F* MPV_Vs_PathLength;
234  TH1F* MPV_Vs_PathLength320;
235  TH1F* MPV_Vs_PathLength500;
236 
237  TH1F* FWHM_Vs_PathLength;
238  TH1F* FWHM_Vs_PathLength320;
239  TH1F* FWHM_Vs_PathLength500;
240 
241 
242  TH2F* Charge_Vs_TransversAngle;
243  TH1F* MPV_Vs_TransversAngle;
244  TH2F* NStrips_Vs_TransversAngle;
245 
246  TH2F* Charge_Vs_Alpha;
247  TH1F* MPV_Vs_Alpha;
248  TH2F* NStrips_Vs_Alpha;
249 
250  TH2F* Charge_Vs_Beta;
251  TH1F* MPV_Vs_Beta;
252  TH2F* NStrips_Vs_Beta;
253 
254  TH2F* Error_Vs_MPV;
255  TH2F* Error_Vs_Entries;
256  TH2F* Error_Vs_Eta;
257  TH2F* Error_Vs_Phi;
258 
259  TH2F* NoMPV_Vs_EtaPhi;
260 
261  TH2F* HitLocalPosition;
262  TH2F* HitLocalPositionBefCut;
263 
264  TH1F* JobInfo;
265 
266  TH1F* HFirstStrip;
267 
268  unsigned int NEvent;
269  unsigned int SRun;
270  unsigned int SEvent;
271  TimeValue_t STimestamp;
272  unsigned int ERun;
273  unsigned int EEvent;
274  TimeValue_t ETimestamp;
275 
276  private :
277  class isEqual{
278  public:
279  template <class T> bool operator () (const T& PseudoDetId1, const T& PseudoDetId2) { return PseudoDetId1==PseudoDetId2; }
280  };
281 
282  std::vector<stAPVGain*> APVsCollOrdered;
283  std::unordered_map<unsigned int, stAPVGain*> APVsColl;
284 };
285 
286 SiStripGainFromData::SiStripGainFromData(const edm::ParameterSet& iConfig) : ConditionDBWriter<SiStripApvGain>(iConfig)
287 {
288  AlgoMode = iConfig.getParameter<std::string>("AlgoMode");
289 
290  OutputGains = iConfig.getParameter<std::string>("OutputGains");
291  OutputHistos = iConfig.getParameter<std::string>("OutputHistos");
292 
293  TrajToTrackProducer = iConfig.getParameter<std::string>("TrajToTrackProducer");
294  TrajToTrackLabel = iConfig.getParameter<std::string>("TrajToTrackLabel");
295 
296  CheckLocalAngle = iConfig.getUntrackedParameter<bool> ("checkLocalAngle" , false);
297  MinNrEntries = iConfig.getUntrackedParameter<unsigned>("minNrEntries" , 20);
298  MaxMPVError = iConfig.getUntrackedParameter<double> ("maxMPVError" , 500.0);
299  MaxChi2OverNDF = iConfig.getUntrackedParameter<double> ("maxChi2OverNDF" , 5.0);
300  MinTrackMomentum = iConfig.getUntrackedParameter<double> ("minTrackMomentum" , 3.0);
301  MaxTrackMomentum = iConfig.getUntrackedParameter<double> ("maxTrackMomentum" , 99999.0);
302  MinTrackEta = iConfig.getUntrackedParameter<double> ("minTrackEta" , -5.0);
303  MaxTrackEta = iConfig.getUntrackedParameter<double> ("maxTrackEta" , 5.0);
304  MaxNrStrips = iConfig.getUntrackedParameter<unsigned>("maxNrStrips" , 2);
305  MinTrackHits = iConfig.getUntrackedParameter<unsigned>("MinTrackHits" , 8);
306  MaxTrackChiOverNdf = iConfig.getUntrackedParameter<double> ("MaxTrackChiOverNdf" , 3);
307  AllowSaturation = iConfig.getUntrackedParameter<bool> ("AllowSaturation" , false);
308  FirstSetOfConstants = iConfig.getUntrackedParameter<bool> ("FirstSetOfConstants", true);
309  Validation = iConfig.getUntrackedParameter<bool> ("Validation" , false);
310  CheckIfFileExist = iConfig.getUntrackedParameter<bool> ("CheckIfFileExist" , false);
311 
312  CalibrationLevel = iConfig.getUntrackedParameter<int> ("CalibrationLevel" , 0);
313 
314 
315  if( strcmp(AlgoMode.c_str(),"WriteOnDB")==0 )
316  VInputFiles = iConfig.getParameter<vector<string> >("VInputFiles");
317 
318  dqmStore_ = edm::Service<DQMStore>().operator->();
319  dqmStore_infile = edm::Service<DQMStore>().operator->();
320 
321  //if( OutputHistos!="" )
322  // dqmStore_->open(OutputHistos.c_str(), true);
323 }
324 
325 
326 SiStripGainFromData::~SiStripGainFromData()
327 {
328 }
329 
330 
331 void
332 SiStripGainFromData::algoBeginJob(const edm::EventSetup& iSetup)
333 {
334  //Retrieve tracker topology from geometry
335  edm::ESHandle<TrackerTopology> tTopoHandle;
336  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
337  const TrackerTopology* const tTopo = tTopoHandle.product();
338 
339  iSetup_ = &iSetup;
340 
341 // TH1::AddDirectory(kTRUE);
342 
343  tmp = dqmStore_->book1D ("JobInfo" , "JobInfo", 20,0,20); JobInfo = tmp->getTH1F();
344 
345  tmp = dqmStore_->book1D ("APV_DetId" , "APV_DetId" , 72785,0,72784); APV_DetId = tmp->getTH1F();
346  tmp = dqmStore_->book1D ("APV_Id" , "APV_Id" , 72785,0,72784); APV_Id = tmp->getTH1F();
347  tmp = dqmStore_->book1D ("APV_Eta" , "APV_Eta" , 72785,0,72784); APV_Eta = tmp->getTH1F();
348  tmp = dqmStore_->book1D ("APV_R" , "APV_R" , 72785,0,72784); APV_R = tmp->getTH1F();
349  tmp = dqmStore_->book1D ("APV_SubDet" , "APV_SubDet" , 72785,0,72784); APV_SubDet = tmp->getTH1F();
350  tmp = dqmStore_->book2D ("APV_Momentum" , "APV_Momentum" , 72785,0,72784, 50,0,200); APV_Momentum = tmp->getTH2F();
351  tmp = dqmStore_->book2D ("APV_Charge" , "APV_Charge" , 72785,0,72784, 1000,0,2000); APV_Charge = tmp->getTH2F();
352  tmp = dqmStore_->book2D ("APV_PathLength" , "APV_PathLength" , 72785,0,72784, 100,0.2,1.4); APV_PathLength = tmp->getTH2F();
353  tmp = dqmStore_->book1D ("APV_PathLengthM", "APV_PathLengthM", 72785,0,72784); APV_PathLengthM = tmp->getTH1F();
354  tmp = dqmStore_->book1D ("APV_MPV" , "APV_MPV" , 72785,0,72784); APV_MPV = tmp->getTH1F();
355  tmp = dqmStore_->book1D ("APV_Gain" , "APV_Gain" , 72785,0,72784); APV_Gain = tmp->getTH1F();
356  tmp = dqmStore_->book1D ("APV_PrevGain" , "APV_PrevGain" , 72785,0,72784); APV_PrevGain = tmp->getTH1F();
357  tmp = dqmStore_->book1D ("APV_CumulGain" , "APV_CumulGain" , 72785,0,72784); APV_CumulGain = tmp->getTH1F();
358  tmp = dqmStore_->book1D ("APV_Thickness" , "APV_Thicknes" , 72785,0,72784); APV_Thickness = tmp->getTH1F();
359 
360 
361  tmp = dqmStore_->book2D ("Tracks_P_Vs_Eta" , "Tracks_P_Vs_Eta" , 30, 0,3,50,0,200); Tracks_P_Vs_Eta = tmp->getTH2F();
362  tmp = dqmStore_->book2D ("Tracks_Pt_Vs_Eta" , "Tracks_Pt_Vs_Eta", 30, 0,3,50,0,200); Tracks_Pt_Vs_Eta = tmp->getTH2F();
363 
364  tmp = dqmStore_->book2D ("Charge_Vs_PathTIB" , "Charge_Vs_PathTIB" ,100,0.2,1.4, 500,0,2000); Charge_Vs_PathTIB = tmp->getTH2F();
365  tmp = dqmStore_->book2D ("Charge_Vs_PathTID" , "Charge_Vs_PathTID" ,100,0.2,1.4, 500,0,2000); Charge_Vs_PathTID = tmp->getTH2F();
366  tmp = dqmStore_->book2D ("Charge_Vs_PathTOB" , "Charge_Vs_PathTOB" ,100,0.2,1.4, 500,0,2000); Charge_Vs_PathTOB = tmp->getTH2F();
367  tmp = dqmStore_->book2D ("Charge_Vs_PathTEC" , "Charge_Vs_PathTEC" ,100,0.2,1.4, 500,0,2000); Charge_Vs_PathTEC = tmp->getTH2F();
368  tmp = dqmStore_->book2D ("Charge_Vs_PathTEC1", "Charge_Vs_PathTEC1",100,0.2,1.4, 500,0,2000); Charge_Vs_PathTEC1 = tmp->getTH2F();
369  tmp = dqmStore_->book2D ("Charge_Vs_PathTEC2", "Charge_Vs_PathTEC2",100,0.2,1.4, 500,0,2000); Charge_Vs_PathTEC2 = tmp->getTH2F();
370 
371 
372  tmp = dqmStore_->book1D ("Charge_TIB" , "Charge_TIB" ,1000,0,2000); Charge_TIB = tmp->getTH1F();
373  tmp = dqmStore_->book1D ("Charge_TID" , "Charge_TID" ,1000,0,2000); Charge_TID = tmp->getTH1F();
374  tmp = dqmStore_->book1D ("Charge_TID+", "Charge_TID+",1000,0,2000); Charge_TIDP = tmp->getTH1F();
375  tmp = dqmStore_->book1D ("Charge_TID-", "Charge_TID-",1000,0,2000); Charge_TIDM = tmp->getTH1F();
376  tmp = dqmStore_->book1D ("Charge_TOB" , "Charge_TOB" ,1000,0,2000); Charge_TOB = tmp->getTH1F();
377  tmp = dqmStore_->book1D ("Charge_TEC" , "Charge_TEC" ,1000,0,2000); Charge_TEC = tmp->getTH1F();
378  tmp = dqmStore_->book1D ("Charge_TEC1", "Charge_TEC1",1000,0,2000); Charge_TEC1 = tmp->getTH1F();
379  tmp = dqmStore_->book1D ("Charge_TEC2", "Charge_TEC2",1000,0,2000); Charge_TEC2 = tmp->getTH1F();
380  tmp = dqmStore_->book1D ("Charge_TEC+", "Charge_TEC+",1000,0,2000); Charge_TECP = tmp->getTH1F();
381  tmp = dqmStore_->book1D ("Charge_TEC-", "Charge_TEC-",1000,0,2000); Charge_TECM = tmp->getTH1F();
382 
383 
384 /*
385  tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS1", "Charge_Vs_PathLength_CS1" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS1 = tmp->getTH2F();
386  tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS2", "Charge_Vs_PathLength_CS2" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS2 = tmp->getTH2F();
387  tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS3", "Charge_Vs_PathLength_CS3" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS3 = tmp->getTH2F();
388  tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS4", "Charge_Vs_PathLength_CS4" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS4 = tmp->getTH2F();
389  tmp = dqmStore_->book2D ("Charge_Vs_PathLength_CS5", "Charge_Vs_PathLength_CS5" , 250,0.2,1.4, 500,0,2000); Charge_Vs_PathLength_CS5 = tmp->getTH2F();
390 */
391  tmp = dqmStore_->book2D ("Charge_Vs_PathLength" , "Charge_Vs_PathLength" , 100,0.2,1.4, 500,0,2000); Charge_Vs_PathLength = tmp->getTH2F();
392  tmp = dqmStore_->book2D ("Charge_Vs_PathLength320" , "Charge_Vs_PathLength" , 100,0.2,1.4, 500,0,2000); Charge_Vs_PathLength320 = tmp->getTH2F();
393  tmp = dqmStore_->book2D ("Charge_Vs_PathLength500" , "Charge_Vs_PathLength" , 100,0.2,1.4, 500,0,2000); Charge_Vs_PathLength500 = tmp->getTH2F();
394 
395  tmp = dqmStore_->book2D ("Charge_Vs_TransversAngle" , "Charge_Vs_TransversAngle" , 220,-20,200, 500,0,2000); Charge_Vs_TransversAngle = tmp->getTH2F();
396  tmp = dqmStore_->book2D ("Charge_Vs_Alpha" , "Charge_Vs_Alpha" , 220,-20,200, 500,0,2000); Charge_Vs_Alpha = tmp->getTH2F();
397  tmp = dqmStore_->book2D ("Charge_Vs_Beta" , "Charge_Vs_Beta" , 220,-20,200, 500,0,2000); Charge_Vs_Beta = tmp->getTH2F();
398 
399  tmp = dqmStore_->book2D ("NStrips_Vs_TransversAngle", "NStrips_Vs_TransversAngle", 220,-20,200, 10,0,10); NStrips_Vs_TransversAngle = tmp->getTH2F();
400  tmp = dqmStore_->book2D ("NStrips_Vs_Alpha" , "NStrips_Vs_Alpha" , 220,-20,200, 10,0,10); NStrips_Vs_Alpha = tmp->getTH2F();
401  tmp = dqmStore_->book2D ("NStrips_Vs_Beta" , "NStrips_Vs_Beta" , 220,-20,200, 10,0,10); NStrips_Vs_Beta = tmp->getTH2F();
402  tmp = dqmStore_->book1D ("NHighStripInCluster" , "NHighStripInCluster" , 15,0,14); NHighStripInCluster = tmp->getTH1F();
403  tmp = dqmStore_->book1D ("NSatStripInCluster" , "NSatStripInCluster" , 50,0,50); NSatStripInCluster = tmp->getTH1F();
404 
405  tmp = dqmStore_->book1D ("TrackChi2OverNDF","TrackChi2OverNDF", 500, 0,10); HTrackChi2OverNDF = tmp->getTH1F();
406  tmp = dqmStore_->book1D ("TrackHits","TrackHits", 40, 0,40); HTrackHits = tmp->getTH1F();
407 
408  tmp = dqmStore_->book1D ("FirstStrip","FirstStrip", 800, 0,800); HFirstStrip = tmp->getTH1F();
409 
410  if( strcmp(AlgoMode.c_str(),"MultiJob")!=0 ){
411 
412  tmp = dqmStore_->book2D ("MPV_Vs_EtaTIB" , "MPV_Vs_EtaTIB" , 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTIB = tmp->getTH2F();
413  tmp = dqmStore_->book2D ("MPV_Vs_EtaTID" , "MPV_Vs_EtaTID" , 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTID = tmp->getTH2F();
414  tmp = dqmStore_->book2D ("MPV_Vs_EtaTOB" , "MPV_Vs_EtaTOB" , 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTOB = tmp->getTH2F();
415  tmp = dqmStore_->book2D ("MPV_Vs_EtaTEC" , "MPV_Vs_EtaTEC" , 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTEC = tmp->getTH2F();
416  tmp = dqmStore_->book2D ("MPV_Vs_EtaTEC1" , "MPV_Vs_EtaTEC1", 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTEC1 = tmp->getTH2F();
417  tmp = dqmStore_->book2D ("MPV_Vs_EtaTEC2" , "MPV_Vs_EtaTEC2", 50, -3.0, 3.0, 300, 0, 600); MPV_Vs_EtaTEC2 = tmp->getTH2F();
418 
419  tmp = dqmStore_->book2D ("MPV_Vs_PhiTIB" , "MPV_Vs_PhiTIB" , 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTIB = tmp->getTH2F();
420  tmp = dqmStore_->book2D ("MPV_Vs_PhiTID" , "MPV_Vs_PhiTID" , 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTID = tmp->getTH2F();
421  tmp = dqmStore_->book2D ("MPV_Vs_PhiTOB" , "MPV_Vs_PhiTOB" , 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTOB = tmp->getTH2F();
422  tmp = dqmStore_->book2D ("MPV_Vs_PhiTEC" , "MPV_Vs_PhiTEC" , 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTEC = tmp->getTH2F();
423  tmp = dqmStore_->book2D ("MPV_Vs_PhiTEC1" , "MPV_Vs_PhiTEC1", 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTEC1 = tmp->getTH2F();
424  tmp = dqmStore_->book2D ("MPV_Vs_PhiTEC2" , "MPV_Vs_PhiTEC2", 50, -3.2, 3.2, 300, 0, 600); MPV_Vs_PhiTEC2 = tmp->getTH2F();
425 
426 
427  tmp = dqmStore_->book1D ("MPV_Vs_PathTIB" , "MPV_Vs_PathTIB" ,100,0.2,1.4); MPV_Vs_PathTIB = tmp->getTH1F();
428  tmp = dqmStore_->book1D ("MPV_Vs_PathTID" , "MPV_Vs_PathTID" ,100,0.2,1.4); MPV_Vs_PathTID = tmp->getTH1F();
429  tmp = dqmStore_->book1D ("MPV_Vs_PathTOB" , "MPV_Vs_PathTOB" ,100,0.2,1.4); MPV_Vs_PathTOB = tmp->getTH1F();
430  tmp = dqmStore_->book1D ("MPV_Vs_PathTEC" , "MPV_Vs_PathTEC" ,100,0.2,1.4); MPV_Vs_PathTEC = tmp->getTH1F();
431  tmp = dqmStore_->book1D ("MPV_Vs_PathTEC1" , "MPV_Vs_PathTEC1" ,100,0.2,1.4); MPV_Vs_PathTEC1 = tmp->getTH1F();
432  tmp = dqmStore_->book1D ("MPV_Vs_PathTEC2" , "MPV_Vs_PathTEC2" ,100,0.2,1.4); MPV_Vs_PathTEC2 = tmp->getTH1F();
433 
434  tmp = dqmStore_->book2D ("MPV_Vs_Phi", "MPV_Vs_Phi", 50, -3.2, 3.2 , 300, 0, 600); MPV_Vs_Phi = tmp->getTH2F();
435  tmp = dqmStore_->book2D ("MPV_Vs_Eta", "MPV_Vs_Eta", 50, -3.0, 3.0 , 300, 0, 600); MPV_Vs_Eta = tmp->getTH2F();
436  tmp = dqmStore_->book2D ("MPV_Vs_R" , "MPV_Vs_R" , 150, 0.0, 150.0, 300, 0, 600); MPV_Vs_R = tmp->getTH2F();
437 /*
438  tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS1" , "MPV_Vs_PathLength_CS1" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS1 = tmp->getTH1F();
439  tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS2" , "MPV_Vs_PathLength_CS2" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS2 = tmp->getTH1F();
440  tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS3" , "MPV_Vs_PathLength_CS3" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS3 = tmp->getTH1F();
441  tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS4" , "MPV_Vs_PathLength_CS4" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS4 = tmp->getTH1F();
442  tmp = dqmStore_->book1D ("MPV_Vs_PathLength_CS5" , "MPV_Vs_PathLength_CS5" , 250, 0.2, 1.4); MPV_Vs_PathLength_CS5 = tmp->getTH1F();
443 
444  tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS1" , "FWHM_Vs_PathLength_CS1", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS1 = tmp->getTH1F();
445  tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS2" , "FWHM_Vs_PathLength_CS2", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS2 = tmp->getTH1F();
446  tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS3" , "FWHM_Vs_PathLength_CS3", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS3 = tmp->getTH1F();
447  tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS4" , "FWHM_Vs_PathLength_CS4", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS4 = tmp->getTH1F();
448  tmp = dqmStore_->book1D ("FWHM_Vs_PathLength_CS5" , "FWHM_Vs_PathLength_CS5", 250, 0.2, 1.4); FWHM_Vs_PathLength_CS5 = tmp->getTH1F();
449 */
450  tmp = dqmStore_->book1D ("MPV_Vs_PathLength" , "MPV_Vs_PathLength" , 100, 0.2, 1.4); MPV_Vs_PathLength = tmp->getTH1F();
451  tmp = dqmStore_->book1D ("MPV_Vs_PathLength320" , "MPV_Vs_PathLength" , 100, 0.2, 1.4); MPV_Vs_PathLength320 = tmp->getTH1F();
452  tmp = dqmStore_->book1D ("MPV_Vs_PathLength500" , "MPV_Vs_PathLength" , 100, 0.2, 1.4); MPV_Vs_PathLength500 = tmp->getTH1F();
453 
454  tmp = dqmStore_->book1D ("FWHM_Vs_PathLength" , "FWHM_Vs_PathLength" , 100, 0.2, 1.4); FWHM_Vs_PathLength = tmp->getTH1F();
455  tmp = dqmStore_->book1D ("FWHM_Vs_PathLength320" , "FWHM_Vs_PathLength" , 100, 0.2, 1.4); FWHM_Vs_PathLength320 = tmp->getTH1F();
456  tmp = dqmStore_->book1D ("FWHM_Vs_PathLength500" , "FWHM_Vs_PathLength" , 100, 0.2, 1.4); FWHM_Vs_PathLength500 = tmp->getTH1F();
457 
458  tmp = dqmStore_->book1D ("MPV_Vs_TransversAngle" , "MPV_Vs_TransversAngle" , 220, -20, 200); MPV_Vs_TransversAngle = tmp->getTH1F();
459  tmp = dqmStore_->book1D ("MPV_Vs_Alpha" , "MPV_Vs_Alpha" , 220, -20, 200); MPV_Vs_Alpha = tmp->getTH1F();
460  tmp = dqmStore_->book1D ("MPV_Vs_Beta" , "MPV_Vs_Beta" , 220, -20, 200); MPV_Vs_Beta = tmp->getTH1F();
461 
462  tmp = dqmStore_->book2D ("Error_Vs_MPV" , "Error_Vs_MPV" ,600,0,600 ,100 ,0 ,50); Error_Vs_MPV = tmp->getTH2F();
463  tmp = dqmStore_->book2D ("Error_Vs_Entries","Error_Vs_Entries",500,0,10000 ,100 ,0 ,50); Error_Vs_Entries = tmp->getTH2F();
464  tmp = dqmStore_->book2D ("Error_Vs_Eta" , "Error_Vs_Eta" ,50 ,-3.0,3.0 ,100 ,0 ,50 ); Error_Vs_Eta = tmp->getTH2F();
465  tmp = dqmStore_->book2D ("Error_Vs_Phi" , "Error_Vs_Phi" ,50 ,-3.2,3.2 ,100 ,0 ,50); Error_Vs_Phi = tmp->getTH2F();
466 
467 
468  tmp = dqmStore_->book2D ("NoMPV_Vs_EtaPhi" , "NoMPV_Vs_EtaPhi" ,50,-3.0,3.0 ,50 ,-3.2,3.2); NoMPV_Vs_EtaPhi = tmp->getTH2F();
469 
470 
471 
472  tmp = dqmStore_->book1D ("NumberOfEntriesByAPV" , "NumberOfEntriesByAPV" , 1000, 0,10000); NumberOfEntriesByAPV = tmp->getTH1F();
473  tmp = dqmStore_->book1D ("Chi2OverNDF","Chi2OverNDF", 500, 0,25); HChi2OverNDF = tmp->getTH1F();
474 
475  tmp = dqmStore_->book1D ("MPVs", "MPVs", 600,0,600); MPVs = tmp->getTH1F();
476  tmp = dqmStore_->book1D ("MPVs320", "MPVs320", 600,0,600); MPVs320 = tmp->getTH1F();
477  tmp = dqmStore_->book1D ("MPVs500", "MPVs500", 600,0,600); MPVs500 = tmp->getTH1F();
478 
479 // MPV_vs_10RplusEta tmp = dqmStore_->book2D ("MPV_vs_10RplusEta","MPV_vs_10RplusEta", 48000,0,2400, 800,100,500);
480  }
481 
482  gROOT->cd();
483 
484  edm::ESHandle<TrackerGeometry> tkGeom;
485  iSetup.get<TrackerDigiGeometryRecord>().get( tkGeom );
486  auto const & Det = tkGeom->dets();
487 
488 
489  edm::ESHandle<SiStripGain> gainHandle;
490 // if(strcmp(AlgoMode.c_str(),"MultiJob")!=0 && !FirstSetOfConstants){
491  iSetup.get<SiStripGainRcd>().get(gainHandle);
492  if(!gainHandle.isValid()){printf("\n#####################\n\nERROR --> gainHandle is not valid\n\n#####################\n\n");exit(0);}
493 // }
494 
495  unsigned int Id=0;
496  for(unsigned int i=0;i<Det.size();i++){
497  DetId Detid = Det[i]->geographicalId();
498  int SubDet = Detid.subdetId();
499 
500  if( SubDet == StripSubdetector::TIB || SubDet == StripSubdetector::TID ||
501  SubDet == StripSubdetector::TOB || SubDet == StripSubdetector::TEC ){
502 
503  auto DetUnit = dynamic_cast<StripGeomDetUnit const*> (Det[i]);
504  if(!DetUnit)continue;
505 
506  const StripTopology& Topo = DetUnit->specificTopology();
507  unsigned int NAPV = Topo.nstrips()/128;
508 
509  double Phi = DetUnit->position().basicVector().phi();
510  double Eta = DetUnit->position().basicVector().eta();
511  double R = DetUnit->position().basicVector().transverse();
512  double Thick = DetUnit->surface().bounds().thickness();
513 
514  for(unsigned int j=0;j<NAPV;j++){
515  stAPVGain* APV = new stAPVGain;
516  APV->Index = Id;
517  APV->DetId = Detid.rawId();
518  APV->APVId = j;
519  APV->SubDet = SubDet;
520  APV->MPV = -1;
521  APV->Gain = -1;
522  APV->PreviousGain = 1;
523  APV->Eta = Eta;
524  APV->Phi = Phi;
525  APV->R = R;
526  APV->Thickness = Thick;
527  APV->Side = 0;
528 
529  if(SubDet==StripSubdetector::TID){
530 
531  APV->Side = tTopo->tecSide(Detid);
532  }else if(SubDet==StripSubdetector::TEC){
533 
534  APV->Side = tTopo->tecSide(Detid);
535  }
536 
537  APVsCollOrdered.push_back(APV);
538  APVsColl[(APV->DetId<<3) | APV->APVId] = APV;
539  Id++;
540 
541  APV_DetId ->Fill(Id,APV->DetId);
542  APV_Id ->Fill(Id,APV->APVId);
543  APV_Eta ->Fill(Id,APV->Eta);
544  APV_R ->Fill(Id,APV->R);
545  APV_SubDet ->Fill(Id,APV->SubDet);
546  APV_Thickness->Fill(Id,APV->Thickness);
547  }
548  }
549  }
550 
551  NEvent = 0;
552  SRun = 0;
553  SEvent = 0;
554  STimestamp = 0;
555  ERun = 0;
556  EEvent = 0;
557  ETimestamp = 0;
558 }
559 
560 void
561 SiStripGainFromData::algoEndJob() {
562 
563 
564  unsigned int I=0;
565 
566  if( strcmp(AlgoMode.c_str(),"WriteOnDB")==0 || strcmp(AlgoMode.c_str(),"Merge")==0){
567  TH1::AddDirectory(kTRUE);
568 
569  TFile* file = nullptr;
570  for(unsigned int f=0;f<VInputFiles.size();f++){
571  printf("Loading New Input File : %s\n", VInputFiles[f].c_str());
572  if(CheckIfFileExist){
573  FILE* doesFileExist = fopen( VInputFiles[f].c_str(), "r" );
574  if(!doesFileExist){
575  printf("File %s doesn't exist\n",VInputFiles[f].c_str());
576  continue;
577  }else{
578  fclose(doesFileExist);
579  }
580  }
581  file = new TFile( VInputFiles[f].c_str() ); if(!file || file->IsZombie() ){printf("### Bug With File %s\n### File will be skipped \n",VInputFiles[f].c_str()); continue;}
582  APV_Charge ->Add( (TH1*) file->FindObjectAny("APV_Charge") , 1);
583  APV_Momentum ->Add( (TH1*) file->FindObjectAny("APV_Momentum") , 1);
584  APV_PathLength ->Add( (TH1*) file->FindObjectAny("APV_PathLength") , 1);
585 
586  Tracks_P_Vs_Eta ->Add( (TH1*) file->FindObjectAny("Tracks_P_Vs_Eta") , 1);
587  Tracks_Pt_Vs_Eta ->Add( (TH1*) file->FindObjectAny("Tracks_Pt_Vs_Eta") , 1);
588 
589  Charge_Vs_PathTIB ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTIB") , 1);
590  Charge_Vs_PathTID ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTID") , 1);
591  Charge_Vs_PathTOB ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTOB") , 1);
592  Charge_Vs_PathTEC ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTEC") , 1);
593  Charge_Vs_PathTEC1 ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTEC1") , 1);
594  Charge_Vs_PathTEC2 ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathTEC2") , 1);
595 
596  HTrackChi2OverNDF ->Add( (TH1*) file->FindObjectAny("TrackChi2OverNDF") , 1);
597  HTrackHits ->Add( (TH1*) file->FindObjectAny("TrackHits") , 1);
598 
599  NHighStripInCluster ->Add( (TH1*) file->FindObjectAny("NHighStripInCluster") , 1);
600  NSatStripInCluster ->Add( (TH1*) file->FindObjectAny("NSatStripInCluster") , 1);
601  Charge_Vs_PathLength ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathLength") , 1);
602  Charge_Vs_PathLength320 ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathLength320") , 1);
603  Charge_Vs_PathLength500 ->Add( (TH1*) file->FindObjectAny("Charge_Vs_PathLength500") , 1);
604  Charge_Vs_TransversAngle ->Add( (TH1*) file->FindObjectAny("Charge_Vs_TransversAngle") , 1);
605  NStrips_Vs_TransversAngle->Add( (TH1*) file->FindObjectAny("NStrips_Vs_TransversAngle"), 1);
606  Charge_Vs_Alpha ->Add( (TH1*) file->FindObjectAny("Charge_Vs_Alpha") , 1);
607  NStrips_Vs_Alpha ->Add( (TH1*) file->FindObjectAny("NStrips_Vs_Alpha") , 1);
608  HFirstStrip ->Add( (TH1*) file->FindObjectAny("FirstStrip") , 1);
609 
610  TH1F* JobInfo_tmp = (TH1F*) file->FindObjectAny("JobInfo");
611  NEvent += (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(1));
612  unsigned int tmp_SRun = (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(3));
613  unsigned int tmp_SEvent = (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(4));
614  unsigned int tmp_ERun = (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(6));
615  unsigned int tmp_EEvent = (unsigned int) JobInfo_tmp->GetBinContent(JobInfo_tmp->GetXaxis()->FindBin(7));
616 
617  if(SRun==0)SRun = tmp_SRun;
618 
619  if (tmp_SRun< SRun){SRun=tmp_SRun; SEvent=tmp_SEvent;}
620  else if(tmp_SRun==SRun && tmp_SEvent<SEvent){SEvent=tmp_SEvent;}
621 
622  if (tmp_ERun> ERun){ERun=tmp_ERun; EEvent=tmp_EEvent;}
623  else if(tmp_ERun==ERun && tmp_EEvent>EEvent){EEvent=tmp_EEvent;}
624 
625  printf("Deleting Current Input File\n");
626  file->Close();
627  delete file;
628  }
629  }
630 
631  JobInfo->Fill(1,NEvent);
632  JobInfo->Fill(3,SRun);
633  JobInfo->Fill(4,SEvent);
634  JobInfo->Fill(6,ERun);
635  JobInfo->Fill(7,EEvent);
636 
637 
638  if( strcmp(AlgoMode.c_str(),"MultiJob")!=0 ){
639  TH1D* Proj = nullptr;
640  double* FitResults = new double[5];
641  I=0;
642  for(auto it = APVsColl.begin();it!=APVsColl.end();it++){
643  if( I%3650==0 ) printf("Fitting Histograms \t %6.2f%%\n",(100.0*I)/APVsColl.size());
644  I++;
645  stAPVGain* APV = it->second;
646 
647  int bin = APV_Charge->GetXaxis()->FindBin(APV->Index);
648  Proj = APV_Charge->ProjectionY(" ",bin,bin,"e");
649  Proj = (TH1D*)Proj->Clone();
650  if(Proj==nullptr)continue;
651 
652  // ADD PROJECTTIONS COMMING FROM THE SECOND APV IN THE PAIR
653  if(CalibrationLevel==1){
654  int SecondAPVId = APV->APVId;
655  if(SecondAPVId%2==0){
656  SecondAPVId = SecondAPVId+1;
657  }else{
658  SecondAPVId = SecondAPVId-1;
659  }
660  stAPVGain* APV2 = APVsColl[(APV->DetId<<3) | SecondAPVId];
661 
662  int bin2 = APV_Charge->GetXaxis()->FindBin(APV2->Index);
663  TH1D* Proj2 = APV_Charge->ProjectionY(" ",bin2,bin2,"e");
664  if(Proj2!=nullptr){
665  Proj->Add(Proj2,1);
666  }
667  }else if(CalibrationLevel>1){
668 // printf("%8i %i--> %4.0f + %4.0f\n",APV->DetId, APV->APVId, 0.0, Proj->GetEntries());
669  for(auto it2 = APVsColl.begin();it2!=APVsColl.end();it2++){
670  stAPVGain* APV2 = it2->second;
671 
672  if(APV2->DetId != APV->DetId)continue;
673  if(APV2->APVId == APV->APVId)continue;
674 
675  int bin2 = APV_Charge->GetXaxis()->FindBin(APV2->Index);
676  TH1D* Proj2 = APV_Charge->ProjectionY(" ",bin2,bin2,"e");
677  if(Proj2!=nullptr){
678 // printf("%8i %i--> %4.0f + %4.0f\n",APV2->DetId, APV2->APVId, Proj->GetEntries(), Proj2->GetEntries());
679  Proj->Add(Proj2,1);
680  }
681  }
682 // printf("%8i %i--> %4.0f Full\n",APV->DetId, APV->APVId, Proj->GetEntries());
683  }
684 
685 
686  //std::cout << "Proj->GetEntries(): " << Proj->GetEntries() << ", Proj->GetMean(): " << Proj->GetMean() << std::endl;
687 
688  getPeakOfLandau(Proj,FitResults);
689  APV->MPV = FitResults[0];
690 // printf("MPV = %f - %f\n",FitResults[0], FitResults[1]);
691  if(FitResults[0]!=-0.5 && FitResults[1]<MaxMPVError){
692  APV_MPV->Fill(APV->Index,APV->MPV);
693  MPVs ->Fill(APV->MPV);
694  if(APV->Thickness<0.04) MPVs320->Fill(APV->MPV);
695  if(APV->Thickness>0.04) MPVs500->Fill(APV->MPV);
696 
697  MPV_Vs_R->Fill(APV->R,APV->MPV);
698  MPV_Vs_Eta->Fill(APV->Eta,APV->MPV);
699  if(APV->SubDet==StripSubdetector::TIB) MPV_Vs_EtaTIB ->Fill(APV->Eta,APV->MPV);
700  if(APV->SubDet==StripSubdetector::TID) MPV_Vs_EtaTID ->Fill(APV->Eta,APV->MPV);
701  if(APV->SubDet==StripSubdetector::TOB) MPV_Vs_EtaTOB ->Fill(APV->Eta,APV->MPV);
702  if(APV->SubDet==StripSubdetector::TEC){ MPV_Vs_EtaTEC ->Fill(APV->Eta,APV->MPV);
703  if(APV->Thickness<0.04) MPV_Vs_EtaTEC1->Fill(APV->Eta,APV->MPV);
704  if(APV->Thickness>0.04) MPV_Vs_EtaTEC2->Fill(APV->Eta,APV->MPV);
705  }
706  MPV_Vs_Phi->Fill(APV->Phi,APV->MPV);
707  if(APV->SubDet==StripSubdetector::TIB) MPV_Vs_PhiTIB ->Fill(APV->Phi,APV->MPV);
708  if(APV->SubDet==StripSubdetector::TID) MPV_Vs_PhiTID ->Fill(APV->Phi,APV->MPV);
709  if(APV->SubDet==StripSubdetector::TOB) MPV_Vs_PhiTOB ->Fill(APV->Phi,APV->MPV);
710  if(APV->SubDet==StripSubdetector::TEC){ MPV_Vs_PhiTEC ->Fill(APV->Phi,APV->MPV);
711  if(APV->Thickness<0.04) MPV_Vs_PhiTEC1->Fill(APV->Phi,APV->MPV);
712  if(APV->Thickness>0.04) MPV_Vs_PhiTEC2->Fill(APV->Phi,APV->MPV);
713  }
714 
715  if(APV->SubDet==StripSubdetector::TIB) Charge_TIB ->Add(Proj,1);
716  if(APV->SubDet==StripSubdetector::TID){ Charge_TID ->Add(Proj,1);
717  if(APV->Side==1) Charge_TIDM->Add(Proj,1);
718  if(APV->Side==2) Charge_TIDP->Add(Proj,1);
719  }
720  if(APV->SubDet==StripSubdetector::TOB) Charge_TOB ->Add(Proj,1);
721  if(APV->SubDet==StripSubdetector::TEC){ Charge_TEC ->Add(Proj,1);
722  if(APV->Thickness<0.04) Charge_TEC1->Add(Proj,1);
723  if(APV->Thickness>0.04) Charge_TEC2->Add(Proj,1);
724  if(APV->Side==1) Charge_TECM->Add(Proj,1);
725  if(APV->Side==2) Charge_TECP->Add(Proj,1);
726  }
727  }
728 
729  if(APV->SubDet==StripSubdetector::TIB) Charge_TIB ->Add(Proj,1);
730  if(APV->SubDet==StripSubdetector::TID){ Charge_TID ->Add(Proj,1);
731  if(APV->Side==1) Charge_TIDM->Add(Proj,1);
732  if(APV->Side==2) Charge_TIDP->Add(Proj,1);
733  }
734  if(APV->SubDet==StripSubdetector::TOB) Charge_TOB ->Add(Proj,1);
735  if(APV->SubDet==StripSubdetector::TEC){ Charge_TEC ->Add(Proj,1);
736  if(APV->Thickness<0.04) Charge_TEC1->Add(Proj,1);
737  if(APV->Thickness>0.04) Charge_TEC2->Add(Proj,1);
738  if(APV->Side==1) Charge_TECM->Add(Proj,1);
739  if(APV->Side==2) Charge_TECP->Add(Proj,1);
740  }
741 
742  if(FitResults[0]!=-0.5){
743  HChi2OverNDF->Fill(FitResults[4]);
744  Error_Vs_MPV->Fill(FitResults[0],FitResults[1]);
745  Error_Vs_Entries->Fill(Proj->GetEntries(),FitResults[1]);
746  Error_Vs_Eta->Fill(APV->Eta,FitResults[1]);
747  Error_Vs_Phi->Fill(APV->Phi,FitResults[1]);
748  }
749  NumberOfEntriesByAPV->Fill(Proj->GetEntries());
750  delete Proj;
751 
752 
753  Proj = APV_PathLength->ProjectionY(" ",bin,bin,"e");
754  if(Proj==nullptr)continue;
755 
756  APV_PathLengthM->SetBinContent(APV->Index, Proj->GetMean(1) );
757  APV_PathLengthM->SetBinError (APV->Index, Proj->GetMeanError(1) );
758 // delete Proj;
759  }
760 
761  unsigned int GOOD = 0;
762  unsigned int BAD = 0;
763  double MPVmean = MPVs->GetMean();
764  MPVmean = 300;
765  for(auto it = APVsColl.begin();it!=APVsColl.end();it++){
766 
767  stAPVGain* APV = it->second;
768  if(APV->MPV>0){
769  APV->Gain = APV->MPV / MPVmean; // APV->MPV;
770  GOOD++;
771  }else{
772  NoMPV_Vs_EtaPhi->Fill(APV->Eta, APV->Phi);
773  APV->Gain = 1;
774  BAD++;
775  }
776  if(APV->Gain<=0) APV->Gain = 1;
777  APV_Gain->Fill(APV->Index,APV->Gain);
778 
779  if(!FirstSetOfConstants) APV->Gain *= APV->PreviousGain;
780  APV_CumulGain->Fill(APV->Index,APV->Gain);
781  }
782 
783  for(int j=0;j<Charge_Vs_PathLength->GetXaxis()->GetNbins();j++){
784  Proj = Charge_Vs_PathLength->ProjectionY(" ",j,j,"e");
785  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
786  MPV_Vs_PathLength->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength->GetXaxis()->GetBinCenter(j));
787  MPV_Vs_PathLength->SetBinError (j, FitResults[1]/Charge_Vs_PathLength->GetXaxis()->GetBinCenter(j));
788  FWHM_Vs_PathLength->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength->GetXaxis()->GetBinCenter(j)) );
789  FWHM_Vs_PathLength->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength->GetXaxis()->GetBinCenter(j)) );
790  delete Proj;
791  }
792 
793  for(int j=0;j<Charge_Vs_PathLength320->GetXaxis()->GetNbins();j++){
794  Proj = Charge_Vs_PathLength320->ProjectionY(" ",j,j,"e");
795  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
796  MPV_Vs_PathLength320->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength320->GetXaxis()->GetBinCenter(j));
797  MPV_Vs_PathLength320->SetBinError (j, FitResults[1]/Charge_Vs_PathLength320->GetXaxis()->GetBinCenter(j));
798  FWHM_Vs_PathLength320->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength320->GetXaxis()->GetBinCenter(j)) );
799  FWHM_Vs_PathLength320->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength320->GetXaxis()->GetBinCenter(j)) );
800  delete Proj;
801  }
802 
803  for(int j=0;j<Charge_Vs_PathLength500->GetXaxis()->GetNbins();j++){
804  Proj = Charge_Vs_PathLength500->ProjectionY(" ",j,j,"e");
805  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
806  MPV_Vs_PathLength500->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength500->GetXaxis()->GetBinCenter(j));
807  MPV_Vs_PathLength500->SetBinError (j, FitResults[1]/Charge_Vs_PathLength500->GetXaxis()->GetBinCenter(j));
808  FWHM_Vs_PathLength500->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength500->GetXaxis()->GetBinCenter(j) ));
809  FWHM_Vs_PathLength500->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength500->GetXaxis()->GetBinCenter(j) ));
810  delete Proj;
811  }
812 /*
813  for(int j=0;j<Charge_Vs_PathLength_CS1->GetXaxis()->GetNbins();j++){
814  Proj = Charge_Vs_PathLength_CS1->ProjectionY(" ",j,j,"e");
815  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
816  MPV_Vs_PathLength_CS1->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS1->GetXaxis()->GetBinCenter(j));
817  MPV_Vs_PathLength_CS1->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS1->GetXaxis()->GetBinCenter(j));
818  FWHM_Vs_PathLength_CS1->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS1->GetXaxis()->GetBinCenter(j) ));
819  FWHM_Vs_PathLength_CS1->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS1->GetXaxis()->GetBinCenter(j) ));
820  delete Proj;
821  }
822 
823  for(int j=0;j<Charge_Vs_PathLength_CS2->GetXaxis()->GetNbins();j++){
824  Proj = Charge_Vs_PathLength_CS2->ProjectionY(" ",j,j,"e");
825  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
826  MPV_Vs_PathLength_CS2->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS2->GetXaxis()->GetBinCenter(j));
827  MPV_Vs_PathLength_CS2->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS2->GetXaxis()->GetBinCenter(j));
828  FWHM_Vs_PathLength_CS2->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS2->GetXaxis()->GetBinCenter(j) ));
829  FWHM_Vs_PathLength_CS2->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS2->GetXaxis()->GetBinCenter(j) ));
830  delete Proj;
831  }
832 
833  for(int j=0;j<Charge_Vs_PathLength_CS3->GetXaxis()->GetNbins();j++){
834  Proj = Charge_Vs_PathLength_CS3->ProjectionY(" ",j,j,"e");
835  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
836  MPV_Vs_PathLength_CS3->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS3->GetXaxis()->GetBinCenter(j));
837  MPV_Vs_PathLength_CS3->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS3->GetXaxis()->GetBinCenter(j));
838  FWHM_Vs_PathLength_CS3->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS3->GetXaxis()->GetBinCenter(j) ));
839  FWHM_Vs_PathLength_CS3->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS3->GetXaxis()->GetBinCenter(j) ));
840  delete Proj;
841  }
842 
843  for(int j=0;j<Charge_Vs_PathLength_CS4->GetXaxis()->GetNbins();j++){
844  Proj = Charge_Vs_PathLength_CS4->ProjectionY(" ",j,j,"e");
845  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
846  MPV_Vs_PathLength_CS4->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS4->GetXaxis()->GetBinCenter(j));
847  MPV_Vs_PathLength_CS4->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS4->GetXaxis()->GetBinCenter(j));
848  FWHM_Vs_PathLength_CS4->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS4->GetXaxis()->GetBinCenter(j) ));
849  FWHM_Vs_PathLength_CS4->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS4->GetXaxis()->GetBinCenter(j) ));
850  delete Proj;
851  }
852 
853  for(int j=0;j<Charge_Vs_PathLength_CS5->GetXaxis()->GetNbins();j++){
854  Proj = Charge_Vs_PathLength_CS5->ProjectionY(" ",j,j,"e");
855  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
856  MPV_Vs_PathLength_CS5->SetBinContent (j, FitResults[0]/Charge_Vs_PathLength_CS5->GetXaxis()->GetBinCenter(j));
857  MPV_Vs_PathLength_CS5->SetBinError (j, FitResults[1]/Charge_Vs_PathLength_CS5->GetXaxis()->GetBinCenter(j));
858  FWHM_Vs_PathLength_CS5->SetBinContent(j, FitResults[2]/(FitResults[0]/Charge_Vs_PathLength_CS5->GetXaxis()->GetBinCenter(j) ));
859  FWHM_Vs_PathLength_CS5->SetBinError (j, FitResults[3]/(FitResults[0]/Charge_Vs_PathLength_CS5->GetXaxis()->GetBinCenter(j) ));
860  delete Proj;
861  }
862 */
863 
864 
865  for(int j=0;j<Charge_Vs_PathTIB->GetXaxis()->GetNbins();j++){
866  Proj = Charge_Vs_PathTIB->ProjectionY(" ",j,j,"e");
867  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
868  MPV_Vs_PathTIB->SetBinContent(j, FitResults[0]/Charge_Vs_PathTIB->GetXaxis()->GetBinCenter(j));
869  MPV_Vs_PathTIB->SetBinError (j, FitResults[1]/Charge_Vs_PathTIB->GetXaxis()->GetBinCenter(j));
870  delete Proj;
871  }
872 
873  for(int j=0;j<Charge_Vs_PathTID->GetXaxis()->GetNbins();j++){
874  Proj = Charge_Vs_PathTID->ProjectionY(" ",j,j,"e");
875  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
876  MPV_Vs_PathTID->SetBinContent(j, FitResults[0]/Charge_Vs_PathTID->GetXaxis()->GetBinCenter(j));
877  MPV_Vs_PathTID->SetBinError (j, FitResults[1]/Charge_Vs_PathTID->GetXaxis()->GetBinCenter(j));
878  delete Proj;
879  }
880 
881  for(int j=0;j<Charge_Vs_PathTOB->GetXaxis()->GetNbins();j++){
882  Proj = Charge_Vs_PathTOB->ProjectionY(" ",j,j,"e");
883  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
884  MPV_Vs_PathTOB->SetBinContent(j, FitResults[0]/Charge_Vs_PathTOB->GetXaxis()->GetBinCenter(j));
885  MPV_Vs_PathTOB->SetBinError (j, FitResults[1]/Charge_Vs_PathTOB->GetXaxis()->GetBinCenter(j));
886  delete Proj;
887  }
888 
889  for(int j=0;j<Charge_Vs_PathTEC->GetXaxis()->GetNbins();j++){
890  Proj = Charge_Vs_PathTEC->ProjectionY(" ",j,j,"e");
891  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
892  MPV_Vs_PathTEC->SetBinContent(j, FitResults[0]/Charge_Vs_PathTEC->GetXaxis()->GetBinCenter(j));
893  MPV_Vs_PathTEC->SetBinError (j, FitResults[1]/Charge_Vs_PathTEC->GetXaxis()->GetBinCenter(j));
894  delete Proj;
895  }
896 
897  for(int j=0;j<Charge_Vs_PathTEC1->GetXaxis()->GetNbins();j++){
898  Proj = Charge_Vs_PathTEC1->ProjectionY(" ",j,j,"e");
899  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
900  MPV_Vs_PathTEC1->SetBinContent(j, FitResults[0]/Charge_Vs_PathTEC1->GetXaxis()->GetBinCenter(j));
901  MPV_Vs_PathTEC1->SetBinError (j, FitResults[1]/Charge_Vs_PathTEC1->GetXaxis()->GetBinCenter(j));
902  delete Proj;
903  }
904 
905  for(int j=0;j<Charge_Vs_PathTEC2->GetXaxis()->GetNbins();j++){
906  Proj = Charge_Vs_PathTEC2->ProjectionY(" ",j,j,"e");
907  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
908  MPV_Vs_PathTEC2->SetBinContent(j, FitResults[0]/Charge_Vs_PathTEC2->GetXaxis()->GetBinCenter(j));
909  MPV_Vs_PathTEC2->SetBinError (j, FitResults[1]/Charge_Vs_PathTEC2->GetXaxis()->GetBinCenter(j));
910  delete Proj;
911  }
912 
913 
914  for(int j=1;j<Charge_Vs_TransversAngle->GetXaxis()->GetNbins();j++){
915  Proj = Charge_Vs_TransversAngle->ProjectionY(" ",j,j,"e");
916  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
917  MPV_Vs_TransversAngle->SetBinContent(j, FitResults[0]);
918  MPV_Vs_TransversAngle->SetBinError (j, FitResults[1]);
919  delete Proj;
920  }
921 
922  for(int j=1;j<Charge_Vs_Alpha->GetXaxis()->GetNbins();j++){
923  Proj = Charge_Vs_Alpha->ProjectionY(" ",j,j,"e");
924  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
925  MPV_Vs_Alpha->SetBinContent(j, FitResults[0]);
926  MPV_Vs_Alpha->SetBinError (j, FitResults[1]);
927  delete Proj;
928  }
929 
930  for(int j=1;j<Charge_Vs_Beta->GetXaxis()->GetNbins();j++){
931  Proj = Charge_Vs_Beta->ProjectionY(" ",j,j,"e");
932  getPeakOfLandau(Proj,FitResults); if(FitResults[0] ==-0.5)continue;
933  MPV_Vs_Beta->SetBinContent(j, FitResults[0]);
934  MPV_Vs_Beta->SetBinError (j, FitResults[1]);
935  delete Proj;
936  }
937 
938  FILE* Gains = fopen(OutputGains.c_str(),"w");
939  fprintf(Gains,"NEvents = %i\n",NEvent);
940  fprintf(Gains,"Number of APVs = %lu\n",static_cast<unsigned long>(APVsColl.size()));
941  fprintf(Gains,"GoodFits = %i BadFits = %i ratio = %f\n",GOOD,BAD,(100.0*GOOD)/(GOOD+BAD));
942  for(std::vector<stAPVGain*>::iterator it = APVsCollOrdered.begin();it!=APVsCollOrdered.end();it++){
943  stAPVGain* APV = *it;
944  fprintf(Gains,"%i | %i | PreviousGain = %7.5f NewGain = %7.5f\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain);
945  }
946 
947  std::vector<int> DetIdOfBuggedAPV;
948  fprintf(Gains,"----------------------------------------------------------------------\n");
949  for(std::vector<stAPVGain*>::iterator it = APVsCollOrdered.begin();it!=APVsCollOrdered.end();it++){
950  stAPVGain* APV = *it;
951  if(APV->MPV>0 && APV->MPV<200){
952  bool tmpBug = false;
953  for(unsigned int b=0;b<DetIdOfBuggedAPV.size()&&!tmpBug;b++){if(DetIdOfBuggedAPV[b]==APV->DetId)tmpBug=true;}
954  if(!tmpBug){fprintf(Gains,"%i,\n",APV->DetId);DetIdOfBuggedAPV.push_back(APV->DetId);}
955  }
956  }
957 
958 
959  fclose(Gains);
960 
961 // delete [] FitResults;
962 // delete Proj;
963  }
964 
965  dqmStore_->cd();
966  dqmStore_->save(OutputHistos);
967 
968 }
969 
970 
971 void SiStripGainFromData::algoBeginRun(const edm::Run &, const edm::EventSetup &iSetup){
972 
973  edm::ESHandle<SiStripGain> gainHandle;
974  if((strcmp(AlgoMode.c_str(),"MultiJob")!=0 && !FirstSetOfConstants) || Validation){
975  iSetup.get<SiStripGainRcd>().get(gainHandle);
976  if(!gainHandle.isValid()){printf("\n#####################\n\nERROR --> gainHandle is not valid\n\n#####################\n\n");exit(0);}
977 
978 
979  for(std::vector<stAPVGain*>::iterator it = APVsCollOrdered.begin();it!=APVsCollOrdered.end();it++){
980  stAPVGain* APV = *it;
981 
982  if(gainHandle.isValid()){
983  SiStripApvGain::Range detGainRange = gainHandle->getRange(APV->DetId);
984  APV->PreviousGain = *(detGainRange.first + APV->APVId);
985 // APV_PrevGain->Fill(APV->Index,APV->PreviousGain);
986  APV_PrevGain->SetBinContent(APV_PrevGain->GetXaxis()->FindBin(APV->Index),APV->PreviousGain);
987  if(APV->PreviousGain<0)APV->PreviousGain = 1;
988  }else{
989  printf("GAIN HANDLE IS NOT VALID\n");
990  }
991  }
992  }
993 
994 
995 }
996 
997 
998 
999 void
1000 SiStripGainFromData::algoAnalyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
1001 {
1002 
1003  if( strcmp(AlgoMode.c_str(),"WriteOnDB")==0 ) return;
1004 
1005  if(NEvent==0){
1006  SRun = iEvent.id().run();
1007  SEvent = iEvent.id().event();
1008  STimestamp = iEvent.time().value();
1009  }
1010  ERun = iEvent.id().run();
1011  EEvent = iEvent.id().event();
1012  ETimestamp = iEvent.time().value();
1013  NEvent++;
1014 
1015  iEvent_ = &iEvent;
1016 
1017  Handle<TrajTrackAssociationCollection> trajTrackAssociationHandle;
1018  iEvent.getByLabel(TrajToTrackProducer, TrajToTrackLabel, trajTrackAssociationHandle);
1019  const TrajTrackAssociationCollection TrajToTrackMap = *trajTrackAssociationHandle.product();
1020 
1021 
1022  for(TrajTrackAssociationCollection::const_iterator it = TrajToTrackMap.begin(); it!=TrajToTrackMap.end(); ++it) {
1023  const Track track = *it->val;
1024  const Trajectory traj = *it->key;
1025 
1026 
1027  if(track.p()<MinTrackMomentum || track.p()>MaxTrackMomentum || track.eta()<MinTrackEta || track.eta()>MaxTrackEta) continue;
1028 
1029  Tracks_Pt_Vs_Eta->Fill(fabs(track.eta()),track.pt());
1030  Tracks_P_Vs_Eta ->Fill(fabs(track.eta()),track.p());
1031 
1032  //BEGIN TO COMPUTE NDOF FOR TRACKS NO IMPLEMENTED BEFORE 200pre3
1033  int ndof =0;
1034  const Trajectory::RecHitContainer transRecHits = traj.recHits();
1035 
1036  for(Trajectory::RecHitContainer::const_iterator rechit = transRecHits.begin(); rechit != transRecHits.end(); ++rechit)
1037  if ((*rechit)->isValid()) ndof += (*rechit)->dimension();
1038  ndof -= 5;
1039  //END TO COMPUTE NDOF FOR TRACKS NO IMPLEMENTED BEFORE 200pre3
1040 
1041  HTrackChi2OverNDF->Fill(traj.chiSquared()/ndof);
1042  if(traj.chiSquared()/ndof>MaxTrackChiOverNdf)continue;
1043 
1044  vector<TrajectoryMeasurement> measurements = traj.measurements();
1045  HTrackHits->Fill(traj.foundHits());
1046  if(traj.foundHits()<(int)MinTrackHits)continue;
1047 /*
1048  //BEGIN TO COMPUTE #MATCHEDRECHIT IN THE TRACK
1049  int NMatchedHit = 0;
1050  for(vector<TrajectoryMeasurement>::const_iterator measurement_it = measurements.begin(); measurement_it!=measurements.end(); measurement_it++){
1051  TrajectoryStateOnSurface trajState = measurement_it->updatedState();
1052  if( !trajState.isValid() ) continue;
1053  const TrackingRecHit* hit = (*measurement_it->recHit()).hit();
1054  const SiStripMatchedRecHit2D* sistripmatchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>(hit);
1055  if(sistripmatchedhit)NMatchedHit++;
1056 // NMatchedHit++;
1057 
1058  }
1059  //END TO COMPUTE #MATCHEDRECHIT IN THE TRACK
1060 
1061  if(NMatchedHit<2){
1062  printf("NOT ENOUGH MATCHED RECHITS : %i\n",NMatchedHit);
1063  continue;
1064  }
1065 */
1066 
1067  for(vector<TrajectoryMeasurement>::const_iterator measurement_it = measurements.begin(); measurement_it!=measurements.end(); measurement_it++){
1068 
1069  TrajectoryStateOnSurface trajState = measurement_it->updatedState();
1070  if( !trajState.isValid() ) continue;
1071 
1072  const TrackingRecHit* hit = (*measurement_it->recHit()).hit();
1073  const SiStripRecHit1D* sistripsimple1dhit = dynamic_cast<const SiStripRecHit1D*>(hit);
1074  const SiStripRecHit2D* sistripsimplehit = dynamic_cast<const SiStripRecHit2D*>(hit);
1075  const SiStripMatchedRecHit2D* sistripmatchedhit = dynamic_cast<const SiStripMatchedRecHit2D*>(hit);
1076 
1077  if(sistripsimplehit){
1078  ComputeChargeOverPath((sistripsimplehit->cluster()).get(), trajState, &iSetup, &track, traj.chiSquared()/ndof);
1079  }else if(sistripmatchedhit){
1080  ComputeChargeOverPath(&sistripmatchedhit->monoCluster(),trajState, &iSetup, &track, traj.chiSquared()/ndof);
1081  ComputeChargeOverPath(&sistripmatchedhit->stereoCluster(),trajState, &iSetup, &track, traj.chiSquared()/ndof);
1082  }else if(sistripsimple1dhit){
1083  ComputeChargeOverPath((sistripsimple1dhit->cluster()).get(), trajState, &iSetup, &track, traj.chiSquared()/ndof);
1084  }else{
1085  }
1086 
1087  }
1088 
1089  }
1090 }
1091 
1092 double
1093 SiStripGainFromData::ComputeChargeOverPath(const SiStripCluster* Cluster ,TrajectoryStateOnSurface trajState, const edm::EventSetup* iSetup, const Track* track, double trajChi2OverN)
1094 {
1095  LocalVector trackDirection = trajState.localDirection();
1096  double cosine = trackDirection.z()/trackDirection.mag();
1097  auto const & Ampls = Cluster->amplitudes();
1098  uint32_t DetId = 0; // is 0 since long time Cluster->geographicalId();
1099  int FirstStrip = Cluster->firstStrip();
1100  int APVId = FirstStrip/128;
1101  stAPVGain* APV = APVsColl[(DetId<<3) | APVId];
1102  int Saturation = 0;
1103  bool Overlaping = false;
1104  int Charge = 0;
1105  unsigned int NHighStrip = 0;
1106 
1107  if(!IsFarFromBorder(trajState,DetId, iSetup))return -1;
1108 
1109 
1110  if(FirstStrip==0 )Overlaping=true;
1111  if(FirstStrip==128 )Overlaping=true;
1112  if(FirstStrip==256 )Overlaping=true;
1113  if(FirstStrip==384 )Overlaping=true;
1114  if(FirstStrip==512 )Overlaping=true;
1115  if(FirstStrip==640 )Overlaping=true;
1116 
1117  if(FirstStrip<=127 && FirstStrip+Ampls.size()>127)Overlaping=true;
1118  if(FirstStrip<=255 && FirstStrip+Ampls.size()>255)Overlaping=true;
1119  if(FirstStrip<=383 && FirstStrip+Ampls.size()>383)Overlaping=true;
1120  if(FirstStrip<=511 && FirstStrip+Ampls.size()>511)Overlaping=true;
1121  if(FirstStrip<=639 && FirstStrip+Ampls.size()>639)Overlaping=true;
1122 
1123  if(FirstStrip+Ampls.size()==127 )Overlaping=true;
1124  if(FirstStrip+Ampls.size()==255 )Overlaping=true;
1125  if(FirstStrip+Ampls.size()==383 )Overlaping=true;
1126  if(FirstStrip+Ampls.size()==511 )Overlaping=true;
1127  if(FirstStrip+Ampls.size()==639 )Overlaping=true;
1128  if(FirstStrip+Ampls.size()==767 )Overlaping=true;
1129  if(Overlaping)return -1;
1130 
1131 
1132 /*
1133  if(FirstStrip==0 )Overlaping=true;
1134  if(FirstStrip<=255 && FirstStrip+Ampls.size()>255)Overlaping=true;
1135  if(FirstStrip<=511 && FirstStrip+Ampls.size()>511)Overlaping=true;
1136  if(FirstStrip+Ampls.size()==511 )Overlaping=true;
1137  if(FirstStrip+Ampls.size()==767 )Overlaping=true;
1138  if(Overlaping)return -1;
1139 */
1140 
1141  for(unsigned int a=0;a<Ampls.size();a++){Charge+=Ampls[a];if(Ampls[a]>=254)Saturation++;if(Ampls[a]>=20)NHighStrip++;}
1142  double path = (10.0*APV->Thickness)/fabs(cosine);
1143  double ClusterChargeOverPath = (double)Charge / path ;
1144 
1145  NSatStripInCluster->Fill(Saturation);
1146 
1147  if(Ampls.size()>MaxNrStrips) return -1;
1148  if(Saturation>0 && !AllowSaturation)return -1;
1149  Charge_Vs_PathLength ->Fill(path,Charge);
1150  if(APV->Thickness<0.04) Charge_Vs_PathLength320->Fill(path,Charge);
1151  if(APV->Thickness>0.04) Charge_Vs_PathLength500->Fill(path,Charge);
1152  if(APV->SubDet==StripSubdetector::TIB) Charge_Vs_PathTIB ->Fill(path,Charge);
1153  if(APV->SubDet==StripSubdetector::TID) Charge_Vs_PathTID ->Fill(path,Charge);
1154  if(APV->SubDet==StripSubdetector::TOB) Charge_Vs_PathTOB ->Fill(path,Charge);
1155  if(APV->SubDet==StripSubdetector::TEC){ Charge_Vs_PathTEC ->Fill(path,Charge);
1156  if(APV->Thickness<0.04) Charge_Vs_PathTEC1 ->Fill(path,Charge);
1157  if(APV->Thickness>0.04) Charge_Vs_PathTEC2 ->Fill(path,Charge); }
1158 
1159  double trans = atan2(trackDirection.y(),trackDirection.x())*(180/3.14159265);
1160  double alpha = acos(trackDirection.x() / sqrt( pow(trackDirection.x(),2) + pow(trackDirection.z(),2) ) ) * (180/3.14159265);
1161  double beta = acos(trackDirection.y() / sqrt( pow(trackDirection.x(),2) + pow(trackDirection.z(),2) ) ) * (180/3.14159265);
1162 
1163  if(path>0.4 && path<0.45){
1164  Charge_Vs_TransversAngle->Fill(trans,Charge/path);
1165  Charge_Vs_Alpha->Fill(alpha,Charge/path);
1166  Charge_Vs_Beta->Fill(beta,Charge/path);
1167  }
1168 
1169  NStrips_Vs_TransversAngle->Fill(trans,Ampls.size());
1170  NStrips_Vs_Alpha ->Fill(alpha,Ampls.size());
1171  NStrips_Vs_Beta ->Fill(beta ,Ampls.size());
1172 
1173  NHighStripInCluster->Fill(NHighStrip);
1174 // if(NHighStrip==1) Charge_Vs_PathLength_CS1->Fill(path, Charge );
1175 // if(NHighStrip==2) Charge_Vs_PathLength_CS2->Fill(path, Charge );
1176 // if(NHighStrip==3) Charge_Vs_PathLength_CS3->Fill(path, Charge );
1177 // if(NHighStrip==4) Charge_Vs_PathLength_CS4->Fill(path, Charge );
1178 // if(NHighStrip==5) Charge_Vs_PathLength_CS5->Fill(path, Charge );
1179 
1180  HFirstStrip ->Fill(FirstStrip);
1181 
1182 
1183  if(Validation){ClusterChargeOverPath=ClusterChargeOverPath/APV->PreviousGain;}
1184 
1185  APV_Charge ->Fill(APV->Index,ClusterChargeOverPath);
1186  APV_Momentum ->Fill(APV->Index,trajState.globalMomentum().mag());
1187  APV_PathLength->Fill(APV->Index,path);
1188 
1189  return ClusterChargeOverPath;
1190 }
1191 
1192 bool SiStripGainFromData::IsFarFromBorder(TrajectoryStateOnSurface trajState, const uint32_t detid, const edm::EventSetup* iSetup)
1193 {
1194  edm::ESHandle<TrackerGeometry> tkGeom; iSetup->get<TrackerDigiGeometryRecord>().get( tkGeom );
1195 
1196  LocalPoint HitLocalPos = trajState.localPosition();
1197  LocalError HitLocalError = trajState.localError().positionError() ;
1198 
1199  const GeomDetUnit* it = tkGeom->idToDetUnit(DetId(detid));
1200  if (dynamic_cast<const StripGeomDetUnit*>(it)==nullptr && dynamic_cast<const PixelGeomDetUnit*>(it)==nullptr) {
1201  std::cout << "this detID doesn't seem to belong to the Tracker" << std::endl;
1202  return false;
1203  }
1204 
1205  const BoundPlane plane = it->surface();
1206  const TrapezoidalPlaneBounds* trapezoidalBounds( dynamic_cast<const TrapezoidalPlaneBounds*>(&(plane.bounds())));
1207  const RectangularPlaneBounds* rectangularBounds( dynamic_cast<const RectangularPlaneBounds*>(&(plane.bounds())));
1208 
1209  double DistFromBorder = 1.0;
1210  //double HalfWidth = it->surface().bounds().width() /2.0;
1211  double HalfLength = it->surface().bounds().length() /2.0;
1212 
1213  if(trapezoidalBounds)
1214  {
1215  std::array<const float, 4> const & parameters = (*trapezoidalBounds).parameters();
1216  HalfLength = parameters[3];
1217  //double t = (HalfLength + HitLocalPos.y()) / (2*HalfLength) ;
1218  //HalfWidth = parameters[0] + (parameters[1]-parameters[0]) * t;
1219  }else if(rectangularBounds){
1220  //HalfWidth = it->surface().bounds().width() /2.0;
1221  HalfLength = it->surface().bounds().length() /2.0;
1222  }else{return false;}
1223 
1224 // if (fabs(HitLocalPos.x())+HitLocalError.xx() >= (HalfWidth - DistFromBorder) ) return false;//Don't think is really necessary
1225  if (fabs(HitLocalPos.y())+HitLocalError.yy() >= (HalfLength - DistFromBorder) ) return false;
1226 
1227  return true;
1228 }
1229 
1230 void SiStripGainFromData::getPeakOfLandau(TH1* InputHisto, double* FitResults, double LowRange, double HighRange)
1231 {
1232  double adcs = -0.5;
1233  double adcs_err = 0.;
1234  double width = -0.5;
1235  double width_err = 0;
1236  double chi2overndf = -0.5;
1237 
1238  double nr_of_entries = InputHisto->GetEntries();
1239 
1240  if( (unsigned int)nr_of_entries < MinNrEntries){
1241  FitResults[0] = adcs;
1242  FitResults[1] = adcs_err;
1243  FitResults[2] = width;
1244  FitResults[3] = width_err;
1245  FitResults[4] = chi2overndf;
1246  return;
1247  }
1248 
1249  // perform fit with standard landau
1250  TF1* MyLandau = new TF1("MyLandau","landau",LowRange, HighRange);
1251  MyLandau->SetParameter("MPV",300);
1252 
1253  InputHisto->Fit(MyLandau,"QR WW");
1254 
1255  // MPV is parameter 1 (0=constant, 1=MPV, 2=Sigma)
1256  adcs = MyLandau->GetParameter("MPV");
1257  adcs_err = MyLandau->GetParError(1);
1258  width = MyLandau->GetParameter(2);
1259  width_err = MyLandau->GetParError(2);
1260  chi2overndf = MyLandau->GetChisquare() / MyLandau->GetNDF();
1261 
1262  // if still wrong, give up
1263  if(adcs<2. || chi2overndf>MaxChi2OverNDF){
1264  adcs = -0.5; adcs_err = 0.;
1265  width = -0.5; width_err = 0;
1266  chi2overndf = -0.5;
1267  }
1268 
1269  FitResults[0] = adcs;
1270  FitResults[1] = adcs_err;
1271  FitResults[2] = width;
1272  FitResults[3] = width_err;
1273  FitResults[4] = chi2overndf;
1274 
1275  delete MyLandau;
1276 }
1277 
1278 
1279 
1280 
1281 SiStripApvGain* SiStripGainFromData::getNewObject()
1282 {
1283  cout << "START getNewObject\n";
1284 
1285 // if( !(strcmp(AlgoMode.c_str(),"WriteOnDB")==0 || strcmp(AlgoMode.c_str(),"SingleJob")==0) )return NULL;
1286  if( !(strcmp(AlgoMode.c_str(),"WriteOnDB")==0 || strcmp(AlgoMode.c_str(),"SingleJob")==0) )return new SiStripApvGain();
1287 
1288 
1289  SiStripApvGain * obj = new SiStripApvGain();
1290  std::vector<float>* theSiStripVector = nullptr;
1291  int PreviousDetId = -1;
1292  for(unsigned int a=0;a<APVsCollOrdered.size();a++)
1293  {
1294  stAPVGain* APV = APVsCollOrdered[a];
1295  if(APV==nullptr){ printf("Bug\n"); continue; }
1296  if(APV->DetId != PreviousDetId){
1297  if(theSiStripVector!=nullptr){
1298  SiStripApvGain::Range range(theSiStripVector->begin(),theSiStripVector->end());
1299  if ( !obj->put(PreviousDetId,range) ) printf("Bug to put detId = %i\n",PreviousDetId);
1300  }
1301 
1302  theSiStripVector = new std::vector<float>;
1303  PreviousDetId = APV->DetId;
1304  }
1305  printf("%i | %i | PreviousGain = %7.5f NewGain = %7.5f\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain);
1306  theSiStripVector->push_back(APV->Gain);
1307 // theSiStripVector->push_back(APV->Gain);
1308  }
1309 
1310  if(theSiStripVector!=nullptr){
1311  SiStripApvGain::Range range(theSiStripVector->begin(),theSiStripVector->end());
1312  if ( !obj->put(PreviousDetId,range) ) printf("Bug to put detId = %i\n",PreviousDetId);
1313  }
1314 
1315  cout << "END getNewObject\n";
1316  return obj;
1317 }
1318 
1319 
1320 
1321 DEFINE_FWK_MODULE(SiStripGainFromData);
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
SiStripRecHit1D::cluster
ClusterRef cluster() const
Definition: SiStripRecHit1D.h:26
reco::TrackBase::p
double p() const
momentum vector magnitude
Definition: TrackBase.h:615
DeDxTools::IsFarFromBorder
bool IsFarFromBorder(const TrajectoryStateOnSurface &trajState, const GeomDetUnit *it)
Definition: DeDxTools.cc:311
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
Trajectory::foundHits
int foundHits() const
Definition: Trajectory.h:225
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
alpha
float alpha
Definition: AMPTWrapper.h:95
Bounds::length
virtual float length() const =0
SiStripGainFromData::algoEndJob
void algoEndJob() override
Definition: SiStripGainFromData.cc:561
mps_fire.i
i
Definition: mps_fire.py:269
SiStripGainFromData::SiStripGainFromData
SiStripGainFromData(const edm::ParameterSet &)
Definition: SiStripGainFromData.cc:286
edm::AssociationMap::end
const_iterator end() const
last iterator over the map (read only)
Definition: AssociationMap.h:160
pfBoostedDoubleSVAK8TagInfos_cfi.beta
beta
Definition: pfBoostedDoubleSVAK8TagInfos_cfi.py:7
MonitorElement::getTH1F
TH1F * getTH1F() const
Definition: MonitorElement.cc:1483
SiStripGainFromData::APVsCollOrdered
std::vector< stAPVGain * > APVsCollOrdered
Definition: SiStripGainFromData.cc:282
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
SiStripMatchedRecHit2D::monoCluster
SiStripCluster const & monoCluster() const
Definition: SiStripMatchedRecHit2D.h:39
TrajectoryStateOnSurface::localDirection
LocalVector localDirection() const
Definition: TrajectoryStateOnSurface.h:116
stAPVGain::APVId
unsigned int APVId
Definition: APVGainStruct.h:11
stAPVGain::Thickness
float Thickness
Definition: APVGainStruct.h:19
SiStripGainFromData::IsFarFromBorder
bool IsFarFromBorder(TrajectoryStateOnSurface trajState, const uint32_t detid, const edm::EventSetup *iSetup)
Definition: SiStripGainFromData.cc:1192
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:47
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
SiStripGainFromData::algoBeginRun
void algoBeginRun(const edm::Run &, const edm::EventSetup &) override
Definition: SiStripGainFromData.cc:971
Trajectory::recHits
ConstRecHitContainer recHits() const
Definition: Trajectory.h:204
Surface::bounds
const Bounds & bounds() const
Definition: Surface.h:120
SiStripCluster::firstStrip
uint16_t firstStrip() const
Definition: SiStripCluster.h:41
stAPVGain::DetId
unsigned int DetId
Definition: APVGainStruct.h:10
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:42
SiStripGainFromData::iSetup_
const edm::EventSetup * iSetup_
Definition: SiStripGainFromData.cc:101
LocalTrajectoryError::positionError
LocalError positionError() const
Definition: LocalTrajectoryError.h:83
SiStripApvGain::put
bool put(const uint32_t &detID, Range input)
Definition: SiStripApvGain.cc:8
TrackerGeometry::dets
const DetContainer & dets() const override
Returm a vector of all GeomDet (including all GeomDetUnits)
Definition: TrackerGeometry.h:61
TrackerGeometry::idToDetUnit
const TrackerGeomDet * idToDetUnit(DetId) const override
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: TrackerGeometry.cc:202
Trajectory::measurements
DataContainer const & measurements() const
Definition: Trajectory.h:196
SiStripGainFromData::VInputFiles
vector< string > VInputFiles
Definition: SiStripGainFromData.cc:127
iEvent
int iEvent
Definition: GenABIO.cc:230
stAPVGain::Gain
double Gain
Definition: APVGainStruct.h:26
PV3DBase::mag
T mag() const
Definition: PV3DBase.h:67
reco::TrackBase::eta
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:651
fastjetJetProducer_validation_cfg.DQMStore
DQMStore
Definition: fastjetJetProducer_validation_cfg.py:29
SiStripGainFromData::APVsColl
std::unordered_map< unsigned int, stAPVGain * > APVsColl
Definition: SiStripGainFromData.cc:283
stAPVGain::PreviousGain
double PreviousGain
Definition: APVGainStruct.h:28
LocalError::yy
float yy() const
Definition: LocalError.h:26
stAPVGain::Index
unsigned int Index
Definition: APVGainStruct.h:8
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
reco::TrackBase::pt
double pt() const
track transverse momentum
Definition: TrackBase.h:621
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
FrontierConditions_GlobalTag_cff.file
file
Definition: FrontierConditions_GlobalTag_cff.py:13
DetId::subdetId
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:41
SiStripApvGain::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripApvGain.h:30
Exhume::I
const std::complex< double > I
Definition: I.h:8
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
SiStripRecHit2D::cluster
ClusterRef cluster() const
Definition: SiStripRecHit2D.h:29
hcaldqm::constants::GOOD
double const GOOD
Definition: Constants.h:15
TrajectoryStateOnSurface::localError
const LocalTrajectoryError & localError() const
Definition: TrajectoryStateOnSurface.h:119
SiStripGainFromData::algoBeginJob
void algoBeginJob(const edm::EventSetup &) override
Definition: SiStripGainFromData.cc:332
edm::TimeValue_t
unsigned long long TimeValue_t
Definition: Timestamp.h:28
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:535
MonitorElement::getTH2F
TH2F * getTH2F() const
Definition: MonitorElement.cc:1507
stringResolutionProvider_cfi.bin
bin
set the eta bin as selection string.
Definition: stringResolutionProvider_cfi.py:11
DetId
Definition: DetId.h:18
hcaldqm::constants::BAD
double const BAD
Definition: Constants.h:17
Trajectory::RecHitContainer
ConstRecHitContainer RecHitContainer
Definition: Trajectory.h:46
edm::Handle::product
T const * product() const
Definition: Handle.h:81
GetRecoTauVFromDQM_MC_cff.obj
obj
Definition: GetRecoTauVFromDQM_MC_cff.py:87
SiStripGainFromData::ComputeChargeOverPath
double ComputeChargeOverPath(const SiStripCluster *Cluster, TrajectoryStateOnSurface trajState, const edm::EventSetup *iSetup, const Track *track, double trajChi2OverN)
Definition: SiStripGainFromData.cc:1093
HiIsolationCommonParameters_cff.track
track
Definition: HiIsolationCommonParameters_cff.py:9
ApeEstimator_cff.width
width
Definition: ApeEstimator_cff.py:24
SiStripGainFromData::algoAnalyze
void algoAnalyze(const edm::Event &, const edm::EventSetup &) override
Definition: SiStripGainFromData.cc:1000
StripTopology::nstrips
virtual int nstrips() const =0
b
double b
Definition: hdecay.h:120
Trajectory::chiSquared
float chiSquared() const
Definition: Trajectory.h:262
SiStripGainFromData::getPeakOfLandau
void getPeakOfLandau(TH1 *InputHisto, double *FitResults, double LowRange=0, double HighRange=5400)
Definition: SiStripGainFromData.cc:1230
SiStripMatchedRecHit2D::stereoCluster
SiStripCluster const & stereoCluster() const
Definition: SiStripMatchedRecHit2D.h:36
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:60
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:43
edm
HLT enums.
Definition: AlignableModifier.h:17
TrajectoryStateOnSurface::globalMomentum
GlobalVector globalMomentum() const
Definition: TrajectoryStateOnSurface.h:86
a
double a
Definition: hdecay.h:121
edm::EventSetup::get
T get() const
Definition: EventSetup.h:63
parameters
parameters
Definition: BeamSpot_PayloadInspector.cc:14
gather_cfg.cout
cout
Definition: gather_cfg.py:143
edm::AssociationMap::begin
const_iterator begin() const
first iterator over the map (read only)
Definition: AssociationMap.h:158
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:47
g4SimHits_cfi.Gain
Gain
Definition: g4SimHits_cfi.py:33
T
long double T
Definition: Basic3DVectorLD.h:66
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
stAPVGain::SubDet
unsigned int SubDet
Definition: APVGainStruct.h:12
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
SiStripGainFromData::getNewObject
SiStripApvGain * getNewObject() override
Definition: SiStripGainFromData.cc:1281
edm::Timestamp::value
TimeValue_t value() const
Definition: Timestamp.h:56
SiStripCluster::amplitudes
const std::vector< uint8_t > & amplitudes() const
Definition: SiStripCluster.h:54
edm::EventBase::time
edm::Timestamp time() const
Definition: EventBase.h:61
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
SiStripGainFromData::iEvent_
const edm::Event * iEvent_
Definition: SiStripGainFromData.cc:102
edm::Run
Definition: Run.h:44
TrackerTopology::tecSide
unsigned int tecSide(const DetId &id) const
Definition: TrackerTopology.h:184
SiStripGain::getRange
const SiStripApvGain::Range getRange(uint32_t detID) const
Definition: SiStripGain.h:70
Generated for CMSSW Reference Manual by doxygen 1.8.11