CMS 3D CMS Logo
SiStripGainFromCalibTree.cc
Go to the documentation of this file.
1 // Original Author: Loic QUERTENMONT
2 // Created: Mon Nov 16 08:55:18 CET 2009
3 
4 #include <memory>
5 #include <iostream>
6 
7 #include "FWCore/Framework/interface/Frameworkfwd.h"
8 #include "FWCore/Framework/interface/EDAnalyzer.h"
9 #include "FWCore/Framework/interface/Event.h"
10 #include "FWCore/Framework/interface/MakerMacros.h"
11 #include "FWCore/Framework/interface/ESHandle.h"
12 #include "FWCore/ParameterSet/interface/ParameterSet.h"
13 #include "FWCore/MessageLogger/interface/MessageLogger.h"
14 #include "FWCore/Utilities/interface/Exception.h"
15 #include "FWCore/Utilities/interface/EDGetToken.h"
16 
17 #include "CondFormats/RunInfo/interface/RunInfo.h"
18 
19 #include "Geometry/CommonDetUnit/interface/GeomDet.h"
20 #include "Geometry/CommonDetUnit/interface/GeomDetType.h"
21 #include "Geometry/CommonTopologies/interface/StripTopology.h"
22 #include "DataFormats/GeometrySurface/interface/TrapezoidalPlaneBounds.h"
23 #include "DataFormats/GeometrySurface/interface/RectangularPlaneBounds.h"
24 
25 #include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
26 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
27 #include "Geometry/Records/interface/TrackerTopologyRcd.h"
28 #include "Geometry/TrackerGeometryBuilder/interface/StripGeomDetUnit.h"
29 #include "Geometry/TrackerGeometryBuilder/interface/PixelGeomDetUnit.h"
30 #include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"
31 #include "Geometry/CommonDetUnit/interface/TrackingGeometry.h"
32 
33 #include "DataFormats/SiStripCluster/interface/SiStripClusterCollection.h"
34 
35 #include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
36 #include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
37 
38 #include "DataFormats/FEDRawData/interface/FEDNumbering.h"
39 #include "DataFormats/TrackReco/interface/Track.h"
40 #include "DataFormats/TrackReco/interface/TrackFwd.h"
41 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHit.h"
42 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit1D.h"
43 #include "DataFormats/TrackerRecHit2D/interface/SiStripRecHit2D.h"
44 #include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2D.h"
45 #include "DataFormats/DetId/interface/DetId.h"
46 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
47 #include "DataFormats/TrackReco/interface/DeDxHit.h"
48 #include "DataFormats/TrackReco/interface/TrackDeDxHits.h"
49 #include "DataFormats/TrackerCommon/interface/TrackerTopology.h"
50 
51 #include "CommonTools/ConditionDBWriter/interface/ConditionDBWriter.h"
52 #include "CondFormats/SiStripObjects/interface/SiStripApvGain.h"
53 
54 #include "TrackingTools/PatternTools/interface/Trajectory.h"
55 #include "TrackingTools/PatternTools/interface/TrajTrackAssociation.h"
56 
57 #include "DQMServices/Core/interface/DQMStore.h"
58 #include "FWCore/ServiceRegistry/interface/Service.h"
59 #include "CommonTools/UtilAlgos/interface/TFileService.h"
60 
61 #include "CalibFormats/SiStripObjects/interface/SiStripGain.h"
62 #include "CalibTracker/Records/interface/SiStripGainRcd.h"
63 
64 #include "CalibFormats/SiStripObjects/interface/SiStripQuality.h"
65 #include "CalibTracker/Records/interface/SiStripQualityRcd.h"
66 
67 
68 #include "TFile.h"
69 #include "TObjString.h"
70 #include "TString.h"
71 #include "TH1F.h"
72 #include "TH2S.h"
73 #include "TProfile.h"
74 #include "TF1.h"
75 #include "TROOT.h"
76 #include "TTree.h"
77 #include "TChain.h"
78 
79 // user includes
80 #include "CalibTracker/SiStripChannelGain/interface/APVGainStruct.h"
81 #include "CalibTracker/SiStripChannelGain/interface/APVGainHelpers.h"
82 
83 #include <unordered_map>
84 #include <array>
85 
86 
87 
88 using namespace edm;
89 using namespace reco;
90 using namespace std;
91 using namespace APVGain;
92 
93 class SiStripGainFromCalibTree : public ConditionDBWriter<SiStripApvGain> {
94 public:
95  explicit SiStripGainFromCalibTree(const edm::ParameterSet&);
96  ~SiStripGainFromCalibTree() override;
97 
98 
99 private:
100 
101 
102  void algoBeginRun(const edm::Run& run, const edm::EventSetup& iSetup) override;
103  void algoEndRun (const edm::Run& run, const edm::EventSetup& iSetup) override;
104  void algoBeginJob (const edm::EventSetup& iSetup) override;
105  void algoEndJob () override;
106  void algoAnalyze (const edm::Event &, const edm::EventSetup &) override;
107 
108  int statCollectionFromMode(const char* tag) const;
109  void bookDQMHistos(const char* dqm_dir, const char* tag);
110 
111  bool isBFieldConsistentWithMode( const edm::EventSetup& iSetup) const;
112  void swapBFieldMode(void);
113 
114  void merge(TH2* A, TH2* B); //needed to add histograms with different number of bins
115  void algoAnalyzeTheTree();
116  void algoComputeMPVandGain();
117  void processEvent(); //what really does the job
118 
119  void getPeakOfLandau(TH1* InputHisto, double* FitResults, double LowRange=50, double HighRange=5400);
120  bool IsGoodLandauFit(double* FitResults);
121  void storeOnTree(TFileService* tfs);
122  void qualityMonitor();
123  void MakeCalibrationMap();
124  bool produceTagFilter();
125 
126  template<typename T>
127  inline edm::Handle<T> connect(const T* &ptr, edm::EDGetTokenT<T> token, const edm::Event &evt) {
128  edm::Handle<T> handle;
129  evt.getByToken(token, handle);
130  ptr = handle.product();
131  return handle; //return handle to keep alive pointer (safety first)
132  }
133 
134  SiStripApvGain* getNewObject() override;
135 
136  TFileService *tfs;
137  DQMStore* dbe;
138  double MagFieldCurrentTh;
139  double MinNrEntries;
140  double MaxMPVError;
141  double MaxChi2OverNDF;
142  double MinTrackMomentum;
143  double MaxTrackMomentum;
144  double MinTrackEta;
145  double MaxTrackEta;
146  unsigned int MaxNrStrips;
147  unsigned int MinTrackHits;
148  double MaxTrackChiOverNdf;
149  int MaxTrackingIteration;
150  bool AllowSaturation;
151  bool FirstSetOfConstants;
152  bool Validation;
153  bool OldGainRemoving;
154  int CalibrationLevel;
155 
156  bool saveSummary;
157  bool useCalibration;
158  bool m_harvestingMode;
159  bool m_splitDQMstat;
160  bool doChargeMonitorPerPlane;
161  std::string m_calibrationMode;
162  std::string m_calibrationPath;
163  std::string m_DQMdir;
164  std::vector<std::string> VChargeHisto;
166  double tagCondition_NClusters;
167  double tagCondition_GoodFrac;
168 
169  std::string AlgoMode;
170  std::string OutputGains;
171  vector<string> VInputFiles;
172 
173  //enum statistic_type {None=-1, StdBunch, StdBunch0T, FaABunch, FaABunch0T, IsoBunch, IsoBunch0T, Harvest};
174 
175  std::vector<string> dqm_tag_;
176  std::string booked_dir_;
177 
178  std::vector<MonitorElement*> Charge_Vs_Index;
179  std::array< std::vector<APVGain::APVmon>,7 > Charge_1;
180  std::array< std::vector<APVGain::APVmon>,7 > Charge_2;
181  std::array< std::vector<APVGain::APVmon>,7 > Charge_3;
182  std::array< std::vector<APVGain::APVmon>,7 > Charge_4;
184  std::vector<MonitorElement*> Charge_Vs_PathlengthTIB;
185  std::vector<MonitorElement*> Charge_Vs_PathlengthTOB;
186  std::vector<MonitorElement*> Charge_Vs_PathlengthTIDP;
187  std::vector<MonitorElement*> Charge_Vs_PathlengthTIDM;
188  std::vector<MonitorElement*> Charge_Vs_PathlengthTECP1;
189  std::vector<MonitorElement*> Charge_Vs_PathlengthTECP2;
190  std::vector<MonitorElement*> Charge_Vs_PathlengthTECM1;
191  std::vector<MonitorElement*> Charge_Vs_PathlengthTECM2;
194  //std::vector<MonitorElement*> Charge_Vs_Index_Absolute;
195 
196  //Validation histograms
197  MonitorElement* MPV_Vs_EtaTIB;
198  MonitorElement* MPV_Vs_EtaTID;
199  MonitorElement* MPV_Vs_EtaTOB;
200  MonitorElement* MPV_Vs_EtaTEC;
201  MonitorElement* MPV_Vs_EtaTECthin;
202  MonitorElement* MPV_Vs_EtaTECthick;
204  MonitorElement* MPV_Vs_PhiTIB;
205  MonitorElement* MPV_Vs_PhiTID;
206  MonitorElement* MPV_Vs_PhiTOB;
207  MonitorElement* MPV_Vs_PhiTEC;
208  MonitorElement* MPV_Vs_PhiTECthin;
209  MonitorElement* MPV_Vs_PhiTECthick;
211  MonitorElement* NoMPVmasked;
212  MonitorElement* NoMPVfit;
214  MonitorElement* Gains;
215  MonitorElement* MPVs;
216  MonitorElement* MPVs320;
217  MonitorElement* MPVs500;
218  MonitorElement* MPVsTIB;
219  MonitorElement* MPVsTID;
220  MonitorElement* MPVsTIDP;
221  MonitorElement* MPVsTIDM;
222  MonitorElement* MPVsTOB;
223  MonitorElement* MPVsTEC;
224  MonitorElement* MPVsTECP;
225  MonitorElement* MPVsTECM;
226  MonitorElement* MPVsTECthin;
227  MonitorElement* MPVsTECthick;
228  MonitorElement* MPVsTECP1;
229  MonitorElement* MPVsTECP2;
230  MonitorElement* MPVsTECM1;
231  MonitorElement* MPVsTECM2;
233  MonitorElement* MPVError;
234  MonitorElement* MPVErrorVsMPV;
235  MonitorElement* MPVErrorVsEta;
236  MonitorElement* MPVErrorVsPhi;
237  MonitorElement* MPVErrorVsN;
239  MonitorElement* DiffWRTPrevGainTIB;
240  MonitorElement* DiffWRTPrevGainTID;
241  MonitorElement* DiffWRTPrevGainTOB;
242  MonitorElement* DiffWRTPrevGainTEC;
244  MonitorElement* GainVsPrevGainTIB;
245  MonitorElement* GainVsPrevGainTID;
246  MonitorElement* GainVsPrevGainTOB;
247  MonitorElement* GainVsPrevGainTEC;
249  std::vector<APVGain::APVmon> newCharge;
250 
251 
252  unsigned int NEvent;
253  unsigned int NTrack;
254  unsigned int NClusterStrip;
255  unsigned int NClusterPixel;
256  int NStripAPVs;
257  int NPixelDets;
258  unsigned int SRun;
259  unsigned int ERun;
260  unsigned int GOOD;
261  unsigned int BAD;
262  unsigned int MASKED;
263 
264  //Data members for processing
265 
266  //Event data
267  unsigned int eventnumber =0;
268  unsigned int runnumber =0;
269  const std::vector<bool>* TrigTech =nullptr; edm::EDGetTokenT<std::vector<bool> > TrigTech_token_;
270 
271  // Track data
272  const std::vector<double>* trackchi2ndof =nullptr; edm::EDGetTokenT<std::vector<double> > trackchi2ndof_token_;
273  const std::vector<float>* trackp =nullptr; edm::EDGetTokenT<std::vector<float> > trackp_token_;
274  const std::vector<float>* trackpt =nullptr; edm::EDGetTokenT<std::vector<float> > trackpt_token_;
275  const std::vector<double>* tracketa =nullptr; edm::EDGetTokenT<std::vector<double> > tracketa_token_;
276  const std::vector<double>* trackphi =nullptr; edm::EDGetTokenT<std::vector<double> > trackphi_token_;
277  const std::vector<unsigned int>* trackhitsvalid =nullptr; edm::EDGetTokenT<std::vector<unsigned int> > trackhitsvalid_token_;
278  const std::vector<int>* trackalgo =nullptr; edm::EDGetTokenT<std::vector<int> > trackalgo_token_;
279 
280  // CalibTree data
281  const std::vector<int>* trackindex =nullptr; edm::EDGetTokenT<std::vector<int> > trackindex_token_;
282  const std::vector<unsigned int>* rawid =nullptr; edm::EDGetTokenT<std::vector<unsigned int> > rawid_token_;
283  const std::vector<double>* localdirx =nullptr; edm::EDGetTokenT<std::vector<double> > localdirx_token_;
284  const std::vector<double>* localdiry =nullptr; edm::EDGetTokenT<std::vector<double> > localdiry_token_;
285  const std::vector<double>* localdirz =nullptr; edm::EDGetTokenT<std::vector<double> > localdirz_token_;
286  const std::vector<unsigned short>* firststrip =nullptr; edm::EDGetTokenT<std::vector<unsigned short> > firststrip_token_;
287  const std::vector<unsigned short>* nstrips =nullptr; edm::EDGetTokenT<std::vector<unsigned short> > nstrips_token_;
288  const std::vector<bool>* saturation =nullptr; edm::EDGetTokenT<std::vector<bool> > saturation_token_;
289  const std::vector<bool>* overlapping =nullptr; edm::EDGetTokenT<std::vector<bool> > overlapping_token_;
290  const std::vector<bool>* farfromedge =nullptr; edm::EDGetTokenT<std::vector<bool> > farfromedge_token_;
291  const std::vector<unsigned int>* charge =nullptr; edm::EDGetTokenT<std::vector<unsigned int> > charge_token_;
292  const std::vector<double>* path =nullptr; edm::EDGetTokenT<std::vector<double> > path_token_;
293  const std::vector<double>* chargeoverpath =nullptr; edm::EDGetTokenT<std::vector<double> > chargeoverpath_token_;
294  const std::vector<unsigned char>* amplitude =nullptr; edm::EDGetTokenT<std::vector<unsigned char> > amplitude_token_;
295  const std::vector<double>* gainused =nullptr; edm::EDGetTokenT<std::vector<double> > gainused_token_;
296  const std::vector<double>* gainusedTick =nullptr; edm::EDGetTokenT<std::vector<double> > gainusedTick_token_;
297 
298 
299  string EventPrefix_; //("");
300  string EventSuffix_; //("");
301  string TrackPrefix_; //("track");
302  string TrackSuffix_; //("");
303  string CalibPrefix_; //("GainCalibration");
304  string CalibSuffix_; //("");
305 
306 private :
307 
308  std::vector<stAPVGain*> APVsCollOrdered;
309  std::unordered_map<unsigned int, stAPVGain*> APVsColl;
310  const TrackerTopology* tTopo_;
311 };
312 
313 inline int
314 SiStripGainFromCalibTree::statCollectionFromMode(const char* tag) const
315 {
316  std::vector<string>::const_iterator it=dqm_tag_.begin();
317  while(it!=dqm_tag_.end()) {
318  if(*it==std::string(tag)) return it-dqm_tag_.begin();
319  it++;
320  }
321 
322  if (std::string(tag)=="") return 0; // return StdBunch calibration mode for backward compatibility
323 
324  return None;
325 }
326 
327 void SiStripGainFromCalibTree::merge(TH2* A, TH2* B){
328  if(A->GetNbinsX() == B->GetNbinsX()){
329  A->Add(B);
330  }else{
331  for(int x=0;x<=B->GetNbinsX()+1; x++){
332  for(int y=0;y<=B->GetNbinsY()+1; y++){
333  A->SetBinContent(x,y,A->GetBinContent(x,y)+B->GetBinContent(x,y));
334  }}
335  }
336 }
337 
338 
339 
340 SiStripGainFromCalibTree::SiStripGainFromCalibTree(const edm::ParameterSet& iConfig) : ConditionDBWriter<SiStripApvGain>(iConfig)
341 {
342  OutputGains = iConfig.getParameter<std::string>("OutputGains");
343 
344  AlgoMode = iConfig.getUntrackedParameter<std::string>("AlgoMode", "CalibTree");
345  MagFieldCurrentTh = iConfig.getUntrackedParameter<double> ("MagFieldCurrentTh" , 2000.);
346  MinNrEntries = iConfig.getUntrackedParameter<double> ("minNrEntries" , 20);
347  MaxMPVError = iConfig.getUntrackedParameter<double> ("maxMPVError" , 500.0);
348  MaxChi2OverNDF = iConfig.getUntrackedParameter<double> ("maxChi2OverNDF" , 5.0);
349  MinTrackMomentum = iConfig.getUntrackedParameter<double> ("minTrackMomentum" , 3.0);
350  MaxTrackMomentum = iConfig.getUntrackedParameter<double> ("maxTrackMomentum" , 99999.0);
351  MinTrackEta = iConfig.getUntrackedParameter<double> ("minTrackEta" , -5.0);
352  MaxTrackEta = iConfig.getUntrackedParameter<double> ("maxTrackEta" , 5.0);
353  MaxNrStrips = iConfig.getUntrackedParameter<unsigned>("maxNrStrips" , 2);
354  MinTrackHits = iConfig.getUntrackedParameter<unsigned>("MinTrackHits" , 8);
355  MaxTrackChiOverNdf = iConfig.getUntrackedParameter<double> ("MaxTrackChiOverNdf" , 3);
356  MaxTrackingIteration = iConfig.getUntrackedParameter<int> ("MaxTrackingIteration" , 7);
357  AllowSaturation = iConfig.getUntrackedParameter<bool> ("AllowSaturation" , false);
358  FirstSetOfConstants = iConfig.getUntrackedParameter<bool> ("FirstSetOfConstants" , true);
359  Validation = iConfig.getUntrackedParameter<bool> ("Validation" , false);
360  OldGainRemoving = iConfig.getUntrackedParameter<bool> ("OldGainRemoving" , false);
361 
362  CalibrationLevel = iConfig.getUntrackedParameter<int> ("CalibrationLevel" , 0);
363  VInputFiles = iConfig.getUntrackedParameter<vector<string> > ("InputFiles");
364  VChargeHisto = iConfig.getUntrackedParameter<vector<string> > ("ChargeHisto");
365 
366 
367  useCalibration = iConfig.getUntrackedParameter<bool> ("UseCalibration" , false);
368  m_harvestingMode = iConfig.getUntrackedParameter<bool> ("harvestingMode" , false);
369  m_splitDQMstat = iConfig.getUntrackedParameter<bool> ("splitDQMstat" , false);
370  m_calibrationMode = iConfig.getUntrackedParameter<string> ("calibrationMode" , "StdBunch");
371  m_calibrationPath = iConfig.getUntrackedParameter<string> ("calibrationPath");
372  m_DQMdir = iConfig.getUntrackedParameter<string> ("DQMdir" , "AlCaReco/SiStripGains");
373 
374  tagCondition_NClusters = iConfig.getUntrackedParameter<double> ("NClustersForTagProd" , 2E8);
375  tagCondition_GoodFrac = iConfig.getUntrackedParameter<double> ("GoodFracForTagProd" , 0.95);
376 
377  saveSummary = iConfig.getUntrackedParameter<bool> ("saveSummary" , false);
378 
379  doChargeMonitorPerPlane = iConfig.getUntrackedParameter<bool> ("doChargeMonitorPerPlane", false);
380 
381 
382  // Gather DQM Service
383  dbe = edm::Service<DQMStore>().operator->();
384  dbe->setVerbose(10);
385 
386  //Set the monitoring element tag and store
387  dqm_tag_.reserve(7);
388  dqm_tag_.clear();
389  dqm_tag_.push_back( "StdBunch" ); // statistic collection from Standard Collision Bunch @ 3.8 T
390  dqm_tag_.push_back( "StdBunch0T" ); // statistic collection from Standard Collision Bunch @ 0 T
391  dqm_tag_.push_back( "AagBunch" ); // statistic collection from First Collision After Abort Gap @ 3.8 T
392  dqm_tag_.push_back( "AagBunch0T" ); // statistic collection from First Collision After Abort Gap @ 0 T
393  dqm_tag_.push_back( "IsoMuon" ); // statistic collection from Isolated Muon @ 3.8 T
394  dqm_tag_.push_back( "IsoMuon0T" ); // statistic collection from Isolated Muon @ 0 T
395  dqm_tag_.push_back( "Harvest" ); // statistic collection: Harvest
396 
397  Charge_Vs_Index.insert( Charge_Vs_Index.begin(), dqm_tag_.size(), nullptr);
398  //Charge_Vs_Index_Absolute.insert( Charge_Vs_Index_Absolute.begin(), dqm_tag_.size(), 0);
399  Charge_Vs_PathlengthTIB.insert( Charge_Vs_PathlengthTIB.begin(), dqm_tag_.size(), nullptr);
400  Charge_Vs_PathlengthTOB.insert( Charge_Vs_PathlengthTOB.begin(), dqm_tag_.size(), nullptr);
401  Charge_Vs_PathlengthTIDP.insert( Charge_Vs_PathlengthTIDP.begin(), dqm_tag_.size(), nullptr);
402  Charge_Vs_PathlengthTIDM.insert( Charge_Vs_PathlengthTIDM.begin(), dqm_tag_.size(), nullptr);
403  Charge_Vs_PathlengthTECP1.insert( Charge_Vs_PathlengthTECP1.begin(), dqm_tag_.size(), nullptr);
404  Charge_Vs_PathlengthTECP2.insert( Charge_Vs_PathlengthTECP2.begin(), dqm_tag_.size(), nullptr);
405  Charge_Vs_PathlengthTECM1.insert( Charge_Vs_PathlengthTECM1.begin(), dqm_tag_.size(), nullptr);
406  Charge_Vs_PathlengthTECM2.insert( Charge_Vs_PathlengthTECM2.begin(), dqm_tag_.size(), nullptr);
407 
408 
409 
410 
411  // configure token for gathering the ntuple variables
412  edm::ParameterSet swhallowgain_pset = iConfig.getUntrackedParameter<edm::ParameterSet>("gain");
413 
414  string label = swhallowgain_pset.getUntrackedParameter<string>("label");
415  CalibPrefix_ = swhallowgain_pset.getUntrackedParameter<string>("prefix");
416  CalibSuffix_ = swhallowgain_pset.getUntrackedParameter<string>("suffix");
417 
418  trackindex_token_ = consumes<std::vector<int> >(edm::InputTag(label, CalibPrefix_ + "trackindex" + CalibSuffix_));
419  rawid_token_ = consumes<std::vector<unsigned int> >(edm::InputTag(label, CalibPrefix_ + "rawid" + CalibSuffix_));
420  localdirx_token_ = consumes<std::vector<double> >(edm::InputTag(label, CalibPrefix_ + "localdirx" + CalibSuffix_));
421  localdiry_token_ = consumes<std::vector<double> >(edm::InputTag(label, CalibPrefix_ + "localdiry" + CalibSuffix_));
422  localdirz_token_ = consumes<std::vector<double> >(edm::InputTag(label, CalibPrefix_ + "localdirz" + CalibSuffix_));
423  firststrip_token_ = consumes<std::vector<unsigned short> >(edm::InputTag(label, CalibPrefix_ + "firststrip" + CalibSuffix_));
424  nstrips_token_ = consumes<std::vector<unsigned short> >(edm::InputTag(label, CalibPrefix_ + "nstrips" + CalibSuffix_));
425  saturation_token_ = consumes<std::vector<bool> >(edm::InputTag(label, CalibPrefix_ + "saturation" + CalibSuffix_));
426  overlapping_token_ = consumes<std::vector<bool> >(edm::InputTag(label, CalibPrefix_ + "overlapping" + CalibSuffix_));
427  farfromedge_token_ = consumes<std::vector<bool> >(edm::InputTag(label, CalibPrefix_ + "farfromedge" + CalibSuffix_));
428  charge_token_ = consumes<std::vector<unsigned int> >(edm::InputTag(label, CalibPrefix_ + "charge" + CalibSuffix_));
429  path_token_ = consumes<std::vector<double> >(edm::InputTag(label, CalibPrefix_ + "path" + CalibSuffix_));
430  chargeoverpath_token_ = consumes<std::vector<double> >(edm::InputTag(label, CalibPrefix_ + "chargeoverpath"+ CalibSuffix_));
431  amplitude_token_ = consumes<std::vector<unsigned char> >(edm::InputTag(label, CalibPrefix_ + "amplitude" + CalibSuffix_));
432  gainused_token_ = consumes<std::vector<double> >(edm::InputTag(label, CalibPrefix_ + "gainused" + CalibSuffix_));
433  gainusedTick_token_ = consumes<std::vector<double> >(edm::InputTag(label, CalibPrefix_ + "gainusedTick" + CalibSuffix_));
434 
435  edm::ParameterSet evtinfo_pset = iConfig.getUntrackedParameter<edm::ParameterSet>("evtinfo");
436  label = evtinfo_pset.getUntrackedParameter<string>("label");
437  EventPrefix_ = evtinfo_pset.getUntrackedParameter<string>("prefix");
438  EventSuffix_ = evtinfo_pset.getUntrackedParameter<string>("suffix");
439  TrigTech_token_ = consumes<std::vector<bool> >(edm::InputTag(label, EventPrefix_ + "TrigTech" + EventSuffix_));
440 
441  edm::ParameterSet track_pset = iConfig.getUntrackedParameter<edm::ParameterSet>("tracks");
442  label = track_pset.getUntrackedParameter<string>("label");
443  TrackPrefix_ = track_pset.getUntrackedParameter<string>("prefix");
444  TrackSuffix_ = track_pset.getUntrackedParameter<string>("suffix");
445 
446  trackchi2ndof_token_ = consumes<std::vector<double> >(edm::InputTag(label, TrackPrefix_ + "chi2ndof" + TrackSuffix_));
447  trackp_token_ = consumes<std::vector<float> >(edm::InputTag(label, TrackPrefix_ + "momentum" + TrackSuffix_));
448  trackpt_token_ = consumes<std::vector<float> >(edm::InputTag(label, TrackPrefix_ + "pt" + TrackSuffix_));
449  tracketa_token_ = consumes<std::vector<double> >(edm::InputTag(label, TrackPrefix_ + "eta" + TrackSuffix_));
450  trackphi_token_ = consumes<std::vector<double> >(edm::InputTag(label, TrackPrefix_ + "phi" + TrackSuffix_));
451  trackhitsvalid_token_ = consumes<std::vector<unsigned int> >(edm::InputTag(label, TrackPrefix_ + "hitsvalid" + TrackSuffix_));
452  trackalgo_token_ = consumes<std::vector<int> >(edm::InputTag(label, TrackPrefix_ + "algo" + TrackSuffix_));
453 
454  tTopo_ = nullptr;
455 }
456 
457 void SiStripGainFromCalibTree::bookDQMHistos(const char* dqm_dir, const char* tag)
458 {
459  edm::LogInfo("SiStripGainFromCalibTree") << "Setting " << dqm_dir << "in DQM and booking histograms for tag "
460  << tag << std::endl;
461 
462  if ( strcmp(booked_dir_.c_str(),dqm_dir)!=0 ) {
463  booked_dir_ = dqm_dir;
464  dbe->setCurrentFolder(dqm_dir);
465  }
466 
467  std::string stag(tag);
468  if(!stag.empty() && stag[0]!='_') stag.insert(0,1,'_');
469 
470  std::string cvi = std::string("Charge_Vs_Index") + stag;
471  //std::string cviA = std::string("Charge_Vs_Index_Absolute") + stag;
472  std::string cvpTIB = std::string("Charge_Vs_PathlengthTIB") + stag;
473  std::string cvpTOB = std::string("Charge_Vs_PathlengthTOB") + stag;
474  std::string cvpTIDP = std::string("Charge_Vs_PathlengthTIDP") + stag;
475  std::string cvpTIDM = std::string("Charge_Vs_PathlengthTIDM") + stag;
476  std::string cvpTECP1 = std::string("Charge_Vs_PathlengthTECP1") + stag;
477  std::string cvpTECP2 = std::string("Charge_Vs_PathlengthTECP2") + stag;
478  std::string cvpTECM1 = std::string("Charge_Vs_PathlengthTECM1") + stag;
479  std::string cvpTECM2 = std::string("Charge_Vs_PathlengthTECM2") + stag;
480 
481  int elepos = (m_harvestingMode && AlgoMode=="PCL")? Harvest : statCollectionFromMode(tag);
482 
483  // The cluster charge is stored by exploiting a non uniform binning in order
484  // reduce the histogram memory size. The bin width is relaxed with a falling
485  // exponential function and the bin boundaries are stored in the binYarray.
486  // The binXarray is used to provide as many bins as the APVs.
487  //
488  // More details about this implementations are here:
489  // https://indico.cern.ch/event/649344/contributions/2672267/attachments/1498323/2332518/OptimizeChHisto.pdf
490 
491  std::vector<float> binXarray;
492  binXarray.reserve( NStripAPVs+1 );
493  for(int a=0;a<=NStripAPVs;a++){
494  binXarray.push_back( (float)a );
495  }
496 
497  std::array<float,688> binYarray;
498  double p0 = 5.445;
499  double p1 = 0.002113;
500  double p2 = 69.01576;
501  double y = 0.;
502  for(int b=0;b<687;b++) {
503  binYarray[b] = y;
504  if(y<=902.) y = y + 2.;
505  else y = ( p0 - log(exp(p0-p1*y) - p2*p1)) / p1;
506  }
507  binYarray[687] = 4000.;
508 
509  Charge_Vs_Index[elepos] = dbe->book2S(cvi.c_str() , cvi.c_str() , NStripAPVs, &binXarray[0], 687, binYarray.data());
510  //Charge_Vs_Index_Absolute[elepos] = dbe->book2S(cviA.c_str() , cviA.c_str() , 88625, 0 , 88624,1000,0,4000);
511  Charge_Vs_PathlengthTIB[elepos] = dbe->book2S(cvpTIB.c_str() , cvpTIB.c_str() , 20 , 0.3 , 1.3 , 250,0,2000);
512  Charge_Vs_PathlengthTOB[elepos] = dbe->book2S(cvpTOB.c_str() , cvpTOB.c_str() , 20 , 0.3 , 1.3 , 250,0,2000);
513  Charge_Vs_PathlengthTIDP[elepos] = dbe->book2S(cvpTIDP.c_str() , cvpTIDP.c_str() , 20 , 0.3 , 1.3 , 250,0,2000);
514  Charge_Vs_PathlengthTIDM[elepos] = dbe->book2S(cvpTIDM.c_str() , cvpTIDM.c_str() , 20 , 0.3 , 1.3 , 250,0,2000);
515  Charge_Vs_PathlengthTECP1[elepos] = dbe->book2S(cvpTECP1.c_str(), cvpTECP1.c_str(), 20 , 0.3 , 1.3 , 250,0,2000);
516  Charge_Vs_PathlengthTECP2[elepos] = dbe->book2S(cvpTECP2.c_str(), cvpTECP2.c_str(), 20 , 0.3 , 1.3 , 250,0,2000);
517  Charge_Vs_PathlengthTECM1[elepos] = dbe->book2S(cvpTECM1.c_str(), cvpTECM1.c_str(), 20 , 0.3 , 1.3 , 250,0,2000);
518  Charge_Vs_PathlengthTECM2[elepos] = dbe->book2S(cvpTECM2.c_str(), cvpTECM2.c_str(), 20 , 0.3 , 1.3 , 250,0,2000);
519 
520  //Book Charge monitoring histograms
521  std::vector<std::pair<std::string,std::string>> hnames = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"");
522  for (unsigned int i=0;i<hnames.size();i++){
523  std::string htag = (hnames[i]).first + stag;
524  MonitorElement* monitor = dbe->book1DD( htag.c_str(), (hnames[i]).second.c_str(), 100 , 0. , 1000. );
525  int id = APVGain::subdetectorId((hnames[i]).first);
526  int side = APVGain::subdetectorSide((hnames[i]).first);
527  int plane = APVGain::subdetectorPlane((hnames[i]).first);
528  Charge_1[elepos].push_back( APVGain::APVmon(id,side,plane,monitor) );
529  }
530 
531  hnames = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"woG2");
532  for (unsigned int i=0;i<hnames.size();i++){
533  std::string htag = (hnames[i]).first + stag;
534  MonitorElement* monitor = dbe->book1DD( htag.c_str(), (hnames[i]).second.c_str(), 100 , 0. , 1000. );
535  int id = APVGain::subdetectorId((hnames[i]).first);
536  int side = APVGain::subdetectorSide((hnames[i]).first);
537  int plane = APVGain::subdetectorPlane((hnames[i]).first);
538  Charge_2[elepos].push_back( APVGain::APVmon(id,side,plane,monitor) );
539  }
540 
541  hnames = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"woG1");
542  for (unsigned int i=0;i<hnames.size();i++){
543  std::string htag = (hnames[i]).first + stag;
544  MonitorElement* monitor = dbe->book1DD( htag.c_str(), (hnames[i]).second.c_str(), 100 , 0. , 1000. );
545  int id = APVGain::subdetectorId((hnames[i]).first);
546  int side = APVGain::subdetectorSide((hnames[i]).first);
547  int plane = APVGain::subdetectorPlane((hnames[i]).first);
548  Charge_3[elepos].push_back( APVGain::APVmon(id,side,plane,monitor) );
549  }
550 
551  hnames = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"woG1G2");
552  for (unsigned int i=0;i<hnames.size();i++){
553  std::string htag = (hnames[i]).first + stag;
554  MonitorElement* monitor = dbe->book1DD( htag.c_str(), (hnames[i]).second.c_str(), 100 , 0. , 1000. );
555  int id = APVGain::subdetectorId((hnames[i]).first);
556  int side = APVGain::subdetectorSide((hnames[i]).first);
557  int plane = APVGain::subdetectorPlane((hnames[i]).first);
558  Charge_4[elepos].push_back( APVGain::APVmon(id,side,plane,monitor) );
559  }
560 
561  //Book validation histograms
562  if (m_harvestingMode) {
563  int MPVbin = 300;
564  float MPVmin = 0.;
565  float MPVmax = 600.;
566 
567  MPV_Vs_EtaTIB = dbe->book2DD("MPV_vs_EtaTIB" ,"MPV vs Eta TIB" , 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
568  MPV_Vs_EtaTID = dbe->book2DD("MPV_vs_EtaTID" ,"MPV vs Eta TID" , 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
569  MPV_Vs_EtaTOB = dbe->book2DD("MPV_vs_EtaTOB" ,"MPV vs Eta TOB" , 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
570  MPV_Vs_EtaTEC = dbe->book2DD("MPV_vs_EtaTEC" ,"MPV vs Eta TEC" , 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
571  MPV_Vs_EtaTECthin = dbe->book2DD("MPV_vs_EtaTEC1","MPV vs Eta TEC-thin" , 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
572  MPV_Vs_EtaTECthick = dbe->book2DD("MPV_vs_EtaTEC2","MPV vs Eta TEC-thick", 50, -3.0, 3.0, MPVbin, MPVmin, MPVmax);
573 
574  MPV_Vs_PhiTIB = dbe->book2DD("MPV_vs_PhiTIB" ,"MPV vs Phi TIB" , 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
575  MPV_Vs_PhiTID = dbe->book2DD("MPV_vs_PhiTID" ,"MPV vs Phi TID" , 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
576  MPV_Vs_PhiTOB = dbe->book2DD("MPV_vs_PhiTOB" ,"MPV vs Phi TOB" , 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
577  MPV_Vs_PhiTEC = dbe->book2DD("MPV_vs_PhiTEC" ,"MPV vs Phi TEC" , 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
578  MPV_Vs_PhiTECthin = dbe->book2DD("MPV_vs_PhiTEC1","MPV vs Phi TEC-thin" , 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
579  MPV_Vs_PhiTECthick = dbe->book2DD("MPV_vs_PhiTEC2","MPV vs Phi TEC-thick", 50, -3.4, 3.4, MPVbin, MPVmin, MPVmax);
580 
581  NoMPVfit = dbe->book2DD("NoMPVfit" ,"Modules with bad Landau Fit",350, -350, 350, 240, 0, 120);
582  NoMPVmasked = dbe->book2DD("NoMPVmasked" ,"Masked Modules" ,350, -350, 350, 240, 0, 120);
583 
584  Gains = dbe->book1DD("Gains" ,"Gains" , 300, 0, 2);
585  MPVs = dbe->book1DD("MPVs" ,"MPVs" , MPVbin, MPVmin, MPVmax);
586  MPVs320 = dbe->book1DD("MPV_320" ,"MPV 320 thickness" , MPVbin, MPVmin, MPVmax);
587  MPVs500 = dbe->book1DD("MPV_500" ,"MPV 500 thickness" , MPVbin, MPVmin, MPVmax);
588  MPVsTIB = dbe->book1DD("MPV_TIB" ,"MPV TIB" , MPVbin, MPVmin, MPVmax);
589  MPVsTID = dbe->book1DD("MPV_TID" ,"MPV TID" , MPVbin, MPVmin, MPVmax);
590  MPVsTIDP = dbe->book1DD("MPV_TIDP" ,"MPV TIDP" , MPVbin, MPVmin, MPVmax);
591  MPVsTIDM = dbe->book1DD("MPV_TIDM" ,"MPV TIDM" , MPVbin, MPVmin, MPVmax);
592  MPVsTOB = dbe->book1DD("MPV_TOB" ,"MPV TOB" , MPVbin, MPVmin, MPVmax);
593  MPVsTEC = dbe->book1DD("MPV_TEC" ,"MPV TEC" , MPVbin, MPVmin, MPVmax);
594  MPVsTECP = dbe->book1DD("MPV_TECP" ,"MPV TECP" , MPVbin, MPVmin, MPVmax);
595  MPVsTECM = dbe->book1DD("MPV_TECM" ,"MPV TECM" , MPVbin, MPVmin, MPVmax);
596  MPVsTECthin = dbe->book1DD("MPV_TEC1" ,"MPV TEC1" , MPVbin, MPVmin, MPVmax);
597  MPVsTECthick = dbe->book1DD("MPV_TEC2" ,"MPV TEC2" , MPVbin, MPVmin, MPVmax);
598  MPVsTECP1 = dbe->book1DD("MPV_TECP1" ,"MPV TECP1" , MPVbin, MPVmin, MPVmax);
599  MPVsTECP2 = dbe->book1DD("MPV_TECP2" ,"MPV TECP2" , MPVbin, MPVmin, MPVmax);
600  MPVsTECM1 = dbe->book1DD("MPV_TECM1" ,"MPV TECM1" , MPVbin, MPVmin, MPVmax);
601  MPVsTECM2 = dbe->book1DD("MPV_TECM2" ,"MPV TECM2" , MPVbin, MPVmin, MPVmax);
602 
603  MPVError = dbe->book1DD("MPVError" ,"MPV Error" , 150, 0, 150);
604  MPVErrorVsMPV = dbe->book2DD("MPVErrorVsMPV" ,"MPV Error vs MPV" ,300, 0, 600, 150, 0, 150);
605  MPVErrorVsEta = dbe->book2DD("MPVErrorVsEta" ,"MPV Error vs Eta" , 50, -3.0, 3.0, 150, 0, 150);
606  MPVErrorVsPhi = dbe->book2DD("MPVErrorVsPhi" ,"MPV Error vs Phi" , 50, -3.4, 3.4, 150, 0, 150);
607  MPVErrorVsN = dbe->book2DD("MPVErrorVsN" ,"MPV Error vs N" ,500, 0,1000, 150, 0, 150);
608 
609  DiffWRTPrevGainTIB = dbe->book1DD("DiffWRTPrevGainTIB" ,"Diff w.r.t. PrevGain TIB" , 250, 0,2);
610  DiffWRTPrevGainTID = dbe->book1DD("DiffWRTPrevGainTID" ,"Diff w.r.t. PrevGain TID" , 250, 0,2);
611  DiffWRTPrevGainTOB = dbe->book1DD("DiffWRTPrevGainTOB" ,"Diff w.r.t. PrevGain TOB" , 250, 0,2);
612  DiffWRTPrevGainTEC = dbe->book1DD("DiffWRTPrevGainTEC" ,"Diff w.r.t. PrevGain TEC" , 250, 0,2);
613 
614  GainVsPrevGainTIB = dbe->book2DD("GainVsPrevGainTIB" ,"Gain vs PrevGain TIB" , 100, 0,2, 100, 0,2);
615  GainVsPrevGainTID = dbe->book2DD("GainVsPrevGainTID" ,"Gain vs PrevGain TID" , 100, 0,2, 100, 0,2);
616  GainVsPrevGainTOB = dbe->book2DD("GainVsPrevGainTOB" ,"Gain vs PrevGain TOB" , 100, 0,2, 100, 0,2);
617  GainVsPrevGainTEC = dbe->book2DD("GainVsPrevGainTEC" ,"Gain vs PrevGain TEC" , 100, 0,2, 100, 0,2);
618 
619  std::vector<std::pair<std::string,std::string>> hnames=APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"newG2");
620  for (unsigned int i=0;i<hnames.size();i++){
621  MonitorElement* monitor = dbe->book1DD( (hnames[i]).first.c_str(), (hnames[i]).second.c_str(), 100 , 0. , 1000. );
622  int id = APVGain::subdetectorId((hnames[i]).first);
623  int side = APVGain::subdetectorSide((hnames[i]).first);
624  int plane = APVGain::subdetectorPlane((hnames[i]).first);
625  newCharge.push_back( APVGain::APVmon(id,side,plane,monitor) );
626  }
627 
628  }
629 }
630 
631 void SiStripGainFromCalibTree::algoBeginJob(const edm::EventSetup& iSetup)
632 {
633  edm::LogInfo("SiStripGainFromCalibTree") << "AlgoMode : " << AlgoMode << "\n"
634  << "CalibrationMode : " << m_calibrationMode << "\n"
635  << "HarvestingMode : " << m_harvestingMode << std::endl;
636  //Setup DQM histograms
637  if(AlgoMode != "PCL" or m_harvestingMode) {
638  const char * dqm_dir = "AlCaReco/SiStripGainsHarvesting/";
639  this->bookDQMHistos( dqm_dir, dqm_tag_[statCollectionFromMode(m_calibrationMode.c_str())].c_str() );
640  } else {
641  //Check consistency of calibration Mode and BField only for the ALCAPROMPT in the PCL workflow
642  if (!isBFieldConsistentWithMode(iSetup)) {
643  string prevMode = m_calibrationMode;
644  swapBFieldMode();
645  edm::LogInfo("SiStripGainFromCalibTree") << "Switching calibration mode for endorsing BField status: "
646  << prevMode << " ==> " << m_calibrationMode << std::endl;
647  }
648  std::string dqm_dir = m_DQMdir + ((m_splitDQMstat)? m_calibrationMode:"") + "/";
649  int elem = statCollectionFromMode(m_calibrationMode.c_str());
650  this->bookDQMHistos( dqm_dir.c_str(), dqm_tag_[elem].c_str() );
651  }
652 
653  edm::ESHandle<TrackerTopology> TopoHandle;
654  iSetup.get<TrackerTopologyRcd>().get( TopoHandle );
655  tTopo_ = TopoHandle.product();
656 
657  edm::ESHandle<TrackerGeometry> tkGeom;
658  iSetup.get<TrackerDigiGeometryRecord>().get( tkGeom );
659  auto const & Det = tkGeom->dets();
660 
661 
662  NPixelDets = 0;
663  NStripAPVs = 0;
664  unsigned int Index=0;
665  for(unsigned int i=0;i<Det.size();i++){
666  DetId Detid = Det[i]->geographicalId();
667  int SubDet = Detid.subdetId();
668 
669  if( SubDet == StripSubdetector::TIB || SubDet == StripSubdetector::TID ||
670  SubDet == StripSubdetector::TOB || SubDet == StripSubdetector::TEC ){
671 
672  auto DetUnit = dynamic_cast<const StripGeomDetUnit*> (Det[i]);
673  if(!DetUnit)continue;
674 
675  const StripTopology& Topo = DetUnit->specificTopology();
676  unsigned int NAPV = Topo.nstrips()/128;
677 
678  for(unsigned int j=0;j<NAPV;j++){
679  stAPVGain* APV = new stAPVGain;
680  APV->Index = Index;
681  APV->Bin = -1;
682  APV->DetId = Detid.rawId();
683  APV->APVId = j;
684  APV->SubDet = SubDet;
685  APV->FitMPV = -1;
686  APV->FitMPVErr = -1;
687  APV->FitWidth = -1;
688  APV->FitWidthErr = -1;
689  APV->FitChi2 = -1;
690  APV->FitNorm = -1;
691  APV->Gain = -1;
692  APV->PreviousGain = 1;
693  APV->PreviousGainTick = 1;
694  APV->x = DetUnit->position().basicVector().x();
695  APV->y = DetUnit->position().basicVector().y();
696  APV->z = DetUnit->position().basicVector().z();
697  APV->Eta = DetUnit->position().basicVector().eta();
698  APV->Phi = DetUnit->position().basicVector().phi();
699  APV->R = DetUnit->position().basicVector().transverse();
700  APV->Thickness = DetUnit->surface().bounds().thickness();
701  APV->NEntries = 0;
702  APV->isMasked = false;
703 
704  APVsCollOrdered.push_back(APV);
705  APVsColl[(APV->DetId<<4) | APV->APVId] = APV;
706  Index++;
707  NStripAPVs++;
708  }
709  }
710  }
711 
712  for(unsigned int i=0;i<Det.size();i++){ //Make two loop such that the Pixel information is added at the end --> make transition simpler
713  DetId Detid = Det[i]->geographicalId();
714  int SubDet = Detid.subdetId();
715  if( SubDet == PixelSubdetector::PixelBarrel || SubDet == PixelSubdetector::PixelEndcap ){
716  auto DetUnit = dynamic_cast<const PixelGeomDetUnit*> (Det[i]);
717  if(!DetUnit) continue;
718 
719  const PixelTopology& Topo = DetUnit->specificTopology();
720  unsigned int NROCRow = Topo.nrows()/(80.);
721  unsigned int NROCCol = Topo.ncolumns()/(52.);
722 
723  for(unsigned int j=0;j<NROCRow;j++){
724  for(unsigned int i=0;i<NROCCol;i++){
725 
726  stAPVGain* APV = new stAPVGain;
727  APV->Index = Index;
728  APV->Bin = -1;
729  APV->DetId = Detid.rawId();
730  APV->APVId = (j<<3 | i);
731  APV->SubDet = SubDet;
732  APV->FitMPV = -1;
733  APV->FitMPVErr = -1;
734  APV->FitWidth = -1;
735  APV->FitWidthErr = -1;
736  APV->FitChi2 = -1;
737  APV->Gain = -1;
738  APV->PreviousGain = 1;
739  APV->PreviousGainTick = 1;
740  APV->x = DetUnit->position().basicVector().x();
741  APV->y = DetUnit->position().basicVector().y();
742  APV->z = DetUnit->position().basicVector().z();
743  APV->Eta = DetUnit->position().basicVector().eta();
744  APV->Phi = DetUnit->position().basicVector().phi();
745  APV->R = DetUnit->position().basicVector().transverse();
746  APV->Thickness = DetUnit->surface().bounds().thickness();
747  APV->isMasked = false; //SiPixelQuality_->IsModuleBad(Detid.rawId());
748  APV->NEntries = 0;
749 
750  APVsCollOrdered.push_back(APV);
751  APVsColl[(APV->DetId<<4) | APV->APVId] = APV;
752  Index++;
753  NPixelDets++;
754  }}
755  }
756  }
757 
758 
759  MakeCalibrationMap();
760 
761  NEvent = 0;
762  NTrack = 0;
763  NClusterStrip = 0;
764  NClusterPixel = 0;
765  SRun = 1<<31;
766  ERun = 0;
767  GOOD = 0;
768  BAD = 0;
769  MASKED = 0;
770 }
771 
772 bool SiStripGainFromCalibTree::isBFieldConsistentWithMode( const edm::EventSetup& iSetup) const {
773  edm::ESHandle<RunInfo> runInfo;
774  iSetup.get<RunInfoRcd>().get(runInfo);
775 
776  double average_current = runInfo.product()->m_avg_current;
777  bool isOn = (average_current > MagFieldCurrentTh);
778  bool is0T = (m_calibrationMode.substr( m_calibrationMode.length()-2 )=="0T");
779 
780  return ( (isOn && !is0T) || (!isOn && is0T) );
781 }
782 
783 void SiStripGainFromCalibTree::swapBFieldMode() {
784  if ( m_calibrationMode.substr( m_calibrationMode.length()-2 )=="0T" ) {
785  m_calibrationMode.erase( m_calibrationMode.length()-2,2);
786  } else {
787  m_calibrationMode.append( "0T" );
788  }
789 }
790 
791 void SiStripGainFromCalibTree::algoBeginRun(const edm::Run& run, const edm::EventSetup& iSetup)
792 {
793  if( !m_harvestingMode && AlgoMode=="PCL") {
794  //Check consistency of calibration Mode and BField only for the ALCAPROMPT in the PCL workflow
795  if (!isBFieldConsistentWithMode(iSetup)) {
796  string prevMode = m_calibrationMode;
797  swapBFieldMode();
798  edm::LogInfo("SiStripGainFromCalibTree") << "Switching calibration mode for endorsing BField status: "
799  << prevMode << " ==> " << m_calibrationMode << std::endl;
800  }
801  }
802 
803  edm::ESHandle<SiStripGain> gainHandle;
804  iSetup.get<SiStripGainRcd>().get(gainHandle);
805  if(!gainHandle.isValid()){edm::LogError("SiStripGainFromCalibTree")<< "gainHandle is not valid\n"; exit(0);}
806 
807  edm::ESHandle<SiStripQuality> SiStripQuality_;
808  iSetup.get<SiStripQualityRcd>().get(SiStripQuality_);
809 
810  for(unsigned int a=0;a<APVsCollOrdered.size();a++){
811 
812  stAPVGain* APV = APVsCollOrdered[a];
813 
814  // MM. 03/03/2017 all of this makes sense only for SiStrip (i.e. get me out of here if I am on a pixel ROC)
815  if(APV->SubDet==PixelSubdetector::PixelBarrel || APV->SubDet==PixelSubdetector::PixelEndcap) continue;
816 
817  APV->isMasked = SiStripQuality_->IsApvBad(APV->DetId,APV->APVId);
818  // if(!FirstSetOfConstants){
819 
820  if(gainHandle->getNumberOfTags()!=2){edm::LogError("SiStripGainFromCalibTree")<< "NUMBER OF GAIN TAG IS EXPECTED TO BE 2\n";fflush(stdout);exit(0);};
821  float newPreviousGain = gainHandle->getApvGain(APV->APVId,gainHandle->getRange(APV->DetId, 1),1);
822  if(APV->PreviousGain!=1 and newPreviousGain!=APV->PreviousGain)edm::LogWarning("SiStripGainFromCalibTree")<< "WARNING: ParticleGain in the global tag changed\n";
823  APV->PreviousGain = newPreviousGain;
824 
825  float newPreviousGainTick = gainHandle->getApvGain(APV->APVId,gainHandle->getRange(APV->DetId, 0),0);
826  if(APV->PreviousGainTick!=1 and newPreviousGainTick!=APV->PreviousGainTick){
827  edm::LogWarning("SiStripGainFromCalibTree")<< "WARNING: TickMarkGain in the global tag changed\n"<< std::endl
828  <<" APV->SubDet: "<< APV->SubDet << " APV->APVId:" << APV->APVId << std::endl
829  <<" APV->PreviousGainTick: "<<APV->PreviousGainTick<<" newPreviousGainTick: "<<newPreviousGainTick<<std::endl;
830  }
831  APV->PreviousGainTick = newPreviousGainTick;
832 
833  //printf("DETID = %7i APVID=%1i Previous Gain=%8.4f (G1) x %8.4f (G2)\n",APV->DetId,APV->APVId,APV->PreviousGainTick, APV->PreviousGain);
834  // }
835 
836  }
837 }
838 
839 void SiStripGainFromCalibTree::algoEndRun(const edm::Run& run, const edm::EventSetup& iSetup){
840  if(AlgoMode == "PCL" && !m_harvestingMode) return;//nothing to do in that case
841 
842  if(AlgoMode == "PCL" and m_harvestingMode){
843  // Load the 2D histograms from the DQM objects
844  // When running in AlCaHarvesting mode the histos are already booked and should be just retrieved from
845  // DQMStore so that they can be used in the fit
846 
847  edm::LogInfo("SiStripGainFromCalibTree") << "Starting harvesting statistics" << std::endl;
848 
849  // check the required tag before adding histograms
850  int elepos = statCollectionFromMode(m_calibrationMode.c_str());
851  if (elepos!=Harvest) {
852 
853  //collect statistics from DQM into the related monitored elements
854  std::string stag = m_calibrationMode;
855  if(!stag.empty() && stag[0]!='_') stag.insert(0,1,'_');
856 
857 
858  if (elepos==-1) {
859  //implememt backward compatibility
860  elepos=0;
861  stag="";
862  }
863 
864  std::string DQM_dir = m_DQMdir + ((m_splitDQMstat)? m_calibrationMode:"");
865 
866  std::string cvi = DQM_dir + std::string("/Charge_Vs_Index") + stag;
867  //std::string cviA = DQM_dir + std::string("/Charge_Vs_Index_Absolute") + stag;
868  std::string cvpTIB = DQM_dir + std::string("/Charge_Vs_PathlengthTIB") + stag;
869  std::string cvpTOB = DQM_dir + std::string("/Charge_Vs_PathlengthTOB") + stag;
870  std::string cvpTIDP = DQM_dir + std::string("/Charge_Vs_PathlengthTIDP") + stag;
871  std::string cvpTIDM = DQM_dir + std::string("/Charge_Vs_PathlengthTIDM") + stag;
872  std::string cvpTECP1 = DQM_dir + std::string("/Charge_Vs_PathlengthTECP1") + stag;
873  std::string cvpTECP2 = DQM_dir + std::string("/Charge_Vs_PathlengthTECP2") + stag;
874  std::string cvpTECM1 = DQM_dir + std::string("/Charge_Vs_PathlengthTECM1") + stag;
875  std::string cvpTECM2 = DQM_dir + std::string("/Charge_Vs_PathlengthTECM2") + stag;
876 
877  Charge_Vs_Index[elepos] = dbe->get(cvi);
878  //Charge_Vs_Index_Absolute[elepos] = dbe->get(cviA.c_str());
879  Charge_Vs_PathlengthTIB[elepos] = dbe->get(cvpTIB);
880  Charge_Vs_PathlengthTOB[elepos] = dbe->get(cvpTOB);
881  Charge_Vs_PathlengthTIDP[elepos] = dbe->get(cvpTIDP);
882  Charge_Vs_PathlengthTIDM[elepos] = dbe->get(cvpTIDM);
883  Charge_Vs_PathlengthTECP1[elepos] = dbe->get(cvpTECP1);
884  Charge_Vs_PathlengthTECP2[elepos] = dbe->get(cvpTECP2);
885  Charge_Vs_PathlengthTECM1[elepos] = dbe->get(cvpTECM1);
886  Charge_Vs_PathlengthTECM2[elepos] = dbe->get(cvpTECM2);
887 
888  if (Charge_Vs_Index[elepos]==nullptr) {
889  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cvi.c_str()
890  << ", statistics will not be summed!" << std::endl;
891  } else {
892  merge( (Charge_Vs_Index[Harvest])->getTH2S(), (Charge_Vs_Index[elepos])->getTH2S() );
893  edm::LogInfo("SiStripGainFromCalibTree") << "Harvesting "
894  << (Charge_Vs_Index[elepos])->getTH2S()->GetEntries() << " more clusters" << std::endl;
895  }
896 
897  //if (Charge_Vs_Index_Absolute[elepos]==0) {
898  // edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cviA.c_str()
899  // << ", statistics will not be summed!" << std::endl;
900  //} else merge( (Charge_Vs_Index_Absolute[Harvest])->getTH2S(), (Charge_Vs_Index_Absolute[elepos])->getTH2S() );
901 
902  if (Charge_Vs_PathlengthTIB[elepos]==nullptr) {
903  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cvpTIB.c_str()
904  << ", statistics will not be summed!" << std::endl;
905  } else (Charge_Vs_PathlengthTIB[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTIB[elepos])->getTH2S());
906 
907  if (Charge_Vs_PathlengthTOB[elepos]==nullptr) {
908  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cvpTOB.c_str()
909  << ", statistics will not be summed!" << std::endl;
910  } else (Charge_Vs_PathlengthTOB[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTOB[elepos])->getTH2S());
911 
912  if (Charge_Vs_PathlengthTIDP[elepos]==nullptr) {
913  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cvpTIDP.c_str()
914  << ", statistics will not be summed!" << std::endl;
915  } else (Charge_Vs_PathlengthTIDP[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTIDP[elepos])->getTH2S());
916 
917  if (Charge_Vs_PathlengthTIDM[elepos]==nullptr) {
918  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cvpTIDM.c_str()
919  << ", statistics will not be summed!" << std::endl;
920  } else (Charge_Vs_PathlengthTIDM[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTIDM[elepos])->getTH2S());
921 
922  if (Charge_Vs_PathlengthTECP1[elepos]==nullptr) {
923  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cvpTECP1.c_str()
924  << ", statistics will not be summed!" << std::endl;
925  } else (Charge_Vs_PathlengthTECP1[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTECP1[elepos])->getTH2S());
926 
927  if (Charge_Vs_PathlengthTECP2[elepos]==nullptr) {
928  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cvpTECP2.c_str()
929  << ", statistics will not be summed!" << std::endl;
930  } else (Charge_Vs_PathlengthTECP2[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTECP2[elepos])->getTH2S());
931 
932  if (Charge_Vs_PathlengthTECM1[elepos]==nullptr) {
933  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cvpTECM1.c_str()
934  << ", statistics will not be summed!" << std::endl;
935  } else (Charge_Vs_PathlengthTECM1[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTECM1[elepos])->getTH2S());
936 
937  if (Charge_Vs_PathlengthTECM2[elepos]==nullptr) {
938  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << cvpTECM2.c_str()
939  << ", statistics will not be summed!" << std::endl;
940  } else (Charge_Vs_PathlengthTECM2[Harvest])->getTH2S()->Add((Charge_Vs_PathlengthTECM2[elepos])->getTH2S());
941 
942 
943  // Gather Charge monitoring histograms
944  std::vector<std::pair<std::string,std::string>> tags = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"");
945  for (unsigned int i=0;i<tags.size();i++){
946  std::string tag = DQM_dir + "/" + (tags[i]).first + stag;
947  Charge_1[elepos].push_back( APVGain::APVmon(0,0,0,dbe->get( tag) ) );
948  if ( (Charge_1[elepos][i]).monitor==nullptr ) {
949  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << tag.c_str()
950  << ", statistics will not be summed!" << std::endl;
951  } else (Charge_1[Harvest][i]).monitor->getTH1D()->Add((Charge_1[elepos][i]).monitor->getTH1D());
952  }
953 
954  tags = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"woG2");
955  for (unsigned int i=0;i<tags.size();i++){
956  std::string tag = DQM_dir + "/" + (tags[i]).first + stag;
957  Charge_2[elepos].push_back( APVGain::APVmon(0,0,0,dbe->get( tag) ) );
958  if ( (Charge_2[elepos][i]).monitor==nullptr ) {
959  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << tag.c_str()
960  << ", statistics will not be summed!" << std::endl;
961  } else (Charge_2[Harvest][i]).monitor->getTH1D()->Add((Charge_2[elepos][i]).monitor->getTH1D());
962 
963  }
964 
965  tags = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"woG1");
966  for (unsigned int i=0;i<tags.size();i++){
967  std::string tag = DQM_dir + "/" + (tags[i]).first + stag;
968  Charge_3[elepos].push_back( APVGain::APVmon(0,0,0,dbe->get( tag) ) );
969  if ( (Charge_3[elepos][i]).monitor==nullptr ) {
970  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << tag.c_str()
971  << ", statistics will not be summed!" << std::endl;
972  } else (Charge_3[Harvest][i]).monitor->getTH1D()->Add((Charge_3[elepos][i]).monitor->getTH1D());
973 
974  }
975 
976  tags = APVGain::monHnames(VChargeHisto,doChargeMonitorPerPlane,"woG1G2");
977  for (unsigned int i=0;i<tags.size();i++){
978  std::string tag = DQM_dir + "/" + (tags[i]).first + stag;
979  Charge_4[elepos].push_back( APVGain::APVmon(0,0,0,dbe->get( tag) ) );
980  if ( (Charge_4[elepos][i]).monitor==nullptr ) {
981  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not retrieve " << tag.c_str()
982  << ", statistics will not be summed!" << std::endl;
983  } else (Charge_4[Harvest][i]).monitor->getTH1D()->Add((Charge_4[elepos][i]).monitor->getTH1D());
984  }
985 
986  }
987  }
988 }
989 
990 void
991 SiStripGainFromCalibTree::algoEndJob() {
992  if(AlgoMode == "PCL" && !m_harvestingMode) return;//nothing to do in that case
993 
994  if(AlgoMode == "CalibTree"){
995  edm::LogInfo("SiStripGainFromCalibTree") << "Analyzing calibration tree" << std::endl;
996  // Loop on calibTrees to fill the 2D histograms
997  algoAnalyzeTheTree();
998  }else if(m_harvestingMode){
999  NClusterStrip = (Charge_Vs_Index[Harvest])->getTH2S()->Integral(0,NStripAPVs+1, 0, 99999 );
1000  //NClusterPixel = (Charge_Vs_Index[Harvest])->getTH2S()->Integral(NStripAPVs+2, NStripAPVs+NPixelDets+2, 0, 99999 );
1001  }
1002 
1003  // Now that we have the full statistics we can extract the information of the 2D histograms
1004  algoComputeMPVandGain();
1005 
1006  // Result monitoring
1007  qualityMonitor();
1008 
1009  // Force the DB object writing,
1010  // thus setting the IOV as the first processed run (if timeFromEndRun is set to false)
1011  storeOnDbNow();
1012 
1013  if(AlgoMode != "PCL" or saveSummary){
1014  edm::LogInfo("SiStripGainFromCalibTree") << "Saving summary into root file" << std::endl;
1015 
1016  //also save the 2D monitor elements to this file as TH2D tfs
1017  tfs = edm::Service<TFileService>().operator->();
1018 
1019  //save only the statistics for the calibrationTag
1020  int elepos = statCollectionFromMode(m_calibrationMode.c_str());
1021 
1022  if( Charge_Vs_Index[elepos]!=nullptr ) tfs->make<TH2S> ( *(Charge_Vs_Index[elepos])->getTH2S() );
1023  //if( Charge_Vs_Index_Absolute[elepos]!=0 ) tfs->make<TH2S> ( *(Charge_Vs_Index_Absolute[elepos])->getTH2S() );
1024  if( Charge_Vs_PathlengthTIB[elepos]!=nullptr ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTIB[elepos])->getTH2S() );
1025  if( Charge_Vs_PathlengthTOB[elepos]!=nullptr ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTOB[elepos])->getTH2S() );
1026  if( Charge_Vs_PathlengthTIDP[elepos]!=nullptr ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTIDP[elepos])->getTH2S() );
1027  if( Charge_Vs_PathlengthTIDM[elepos]!=nullptr ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTIDM[elepos])->getTH2S() );
1028  if( Charge_Vs_PathlengthTECP1[elepos]!=nullptr ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTECP1[elepos])->getTH2S() );
1029  if( Charge_Vs_PathlengthTECP2[elepos]!=nullptr ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTECP2[elepos])->getTH2S() );
1030  if( Charge_Vs_PathlengthTECM1[elepos]!=nullptr ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTECM1[elepos])->getTH2S() );
1031  if( Charge_Vs_PathlengthTECM2[elepos]!=nullptr ) tfs->make<TH2S> ( *(Charge_Vs_PathlengthTECM2[elepos])->getTH2S() );
1032 
1033  storeOnTree(tfs);
1034  }
1035 }
1036 
1037 
1038 void SiStripGainFromCalibTree::getPeakOfLandau(TH1* InputHisto, double* FitResults, double LowRange, double HighRange)
1039 {
1040  FitResults[0] = -0.5; //MPV
1041  FitResults[1] = 0; //MPV error
1042  FitResults[2] = -0.5; //Width
1043  FitResults[3] = 0; //Width error
1044  FitResults[4] = -0.5; //Fit Chi2/NDF
1045  FitResults[5] = 0; //Normalization
1046 
1047  if( InputHisto->GetEntries() < MinNrEntries)return;
1048 
1049  // perform fit with standard landau
1050  TF1* MyLandau = new TF1("MyLandau","landau",LowRange, HighRange);
1051  MyLandau->SetParameter(1,300);
1052  InputHisto->Fit(MyLandau,"0QR WW");
1053 
1054  // MPV is parameter 1 (0=constant, 1=MPV, 2=Sigma)
1055  FitResults[0] = MyLandau->GetParameter(1); //MPV
1056  FitResults[1] = MyLandau->GetParError(1); //MPV error
1057  FitResults[2] = MyLandau->GetParameter(2); //Width
1058  FitResults[3] = MyLandau->GetParError(2); //Width error
1059  FitResults[4] = MyLandau->GetChisquare() / MyLandau->GetNDF(); //Fit Chi2/NDF
1060  FitResults[5] = MyLandau->GetParameter(0);
1061 
1062  delete MyLandau;
1063 }
1064 
1065 bool SiStripGainFromCalibTree::IsGoodLandauFit(double* FitResults){
1066  if(FitResults[0] <= 0 )return false;
1067 // if(FitResults[1] > MaxMPVError )return false;
1068 // if(FitResults[4] > MaxChi2OverNDF)return false;
1069  return true;
1070 }
1071 
1072 void SiStripGainFromCalibTree::processEvent() {
1073 
1074  edm::LogInfo("SiStripGainFromCalibTree") << "Processing run " << runnumber
1075  << " and event " << eventnumber
1076  << " for " << m_calibrationMode << " calibration." << std::endl;
1077 
1078  if(runnumber<SRun)SRun=runnumber;
1079  if(runnumber>ERun)ERun=runnumber;
1080 
1081  NEvent++;
1082  NTrack+=(*trackp).size();
1083 
1084  int elepos = statCollectionFromMode(m_calibrationMode.c_str());
1085 
1086  unsigned int FirstAmplitude=0;
1087  for(unsigned int i=0;i<(*chargeoverpath).size();i++){
1088  FirstAmplitude+=(*nstrips)[i];
1089  int TI = (*trackindex)[i];
1090 
1091  //printf("%i - %i - %i %i %i\n", (int)(*rawid)[i], (int)(*firststrip)[i]/128, (int)(*farfromedge)[i], (int)(*overlapping)[i], (int)(*saturation )[i] );
1092  if((*tracketa )[TI] < MinTrackEta )continue;
1093  if((*tracketa )[TI] > MaxTrackEta )continue;
1094  if((*trackp )[TI] < MinTrackMomentum )continue;
1095  if((*trackp )[TI] > MaxTrackMomentum )continue;
1096  if((*trackhitsvalid)[TI] < MinTrackHits )continue;
1097  if((*trackchi2ndof )[TI] > MaxTrackChiOverNdf )continue;
1098  if((*trackalgo )[TI] > MaxTrackingIteration )continue;
1099 
1100  stAPVGain* APV = APVsColl[((*rawid)[i]<<4) | ((*firststrip)[i]/128)]; //works for both strip and pixel thanks to firstStrip encoding for pixel in the calibTree
1101 
1102  if(APV->SubDet>2 && (*farfromedge)[i] == false )continue;
1103  if(APV->SubDet>2 && (*overlapping)[i] == true )continue;
1104  if(APV->SubDet>2 && (*saturation )[i] && !AllowSaturation)continue;
1105  if(APV->SubDet>2 && (*nstrips )[i] > MaxNrStrips )continue;
1106 
1107 
1108  //printf("detId=%7i run=%7i event=%9i charge=%5i cs=%3i\n",(*rawid)[i],runnumber,eventnumber,(*charge)[i],(*nstrips)[i]);
1109 
1110  //double trans = atan2((*localdiry)[i],(*localdirx)[i])*(180/3.14159265);
1111  //double alpha = acos ((*localdirx)[i] / sqrt( pow((*localdirx)[i],2) + pow((*localdirz)[i],2) ) ) * (180/3.14159265);
1112  //double beta = acos ((*localdiry)[i] / sqrt( pow((*localdirx)[i],2) + pow((*localdirz)[i],2) ) ) * (180/3.14159265);
1113 
1114  //printf("NStrip = %i : Charge = %i --> Path = %f --> ChargeOverPath=%f\n",(*nstrips)[i],(*charge)[i],(*path)[i],(*chargeoverpath)[i]);
1115  //printf("Amplitudes: ");
1116  //for(unsigned int a=0;a<(*nstrips)[i];a++){printf("%i ",(*amplitude)[FirstAmplitude+a]);}
1117  //printf("\n");
1118 
1119  if(APV->SubDet>2){NClusterStrip++;}else{NClusterPixel++;}
1120 
1121  int Charge = 0;
1122  if(APV->SubDet>2 && (useCalibration || !FirstSetOfConstants)){
1123  bool Saturation = false;
1124  for(unsigned int s=0;s<(*nstrips)[i];s++){
1125  int StripCharge = (*amplitude)[FirstAmplitude-(*nstrips)[i]+s];
1126  if(useCalibration && !FirstSetOfConstants){ StripCharge=(int)(StripCharge*(APV->PreviousGain/APV->CalibGain));
1127  }else if(useCalibration){ StripCharge=(int)(StripCharge/APV->CalibGain);
1128  }else if(!FirstSetOfConstants){ StripCharge=(int)(StripCharge*APV->PreviousGain);}
1129  if(StripCharge>1024){
1130  StripCharge = 255;
1131  Saturation = true;
1132  }else if(StripCharge>254){
1133  StripCharge = 254;
1134  Saturation = true;
1135  }
1136  Charge += StripCharge;
1137  }
1138  if(Saturation && !AllowSaturation)continue;
1139  }else if(APV->SubDet>2){
1140  Charge = (*charge)[i];
1141  }else{
1142  Charge = (*charge)[i]/265.0; //expected scale factor between pixel and strip charge
1143  }
1144 
1145  //printf("ChargeDifference = %i Vs %i with Gain = %f\n",(*charge)[i],Charge,APV->CalibGain);
1146 
1147  double ClusterChargeOverPath = ( (double) Charge )/(*path)[i] ;
1148  if(APV->SubDet>2){
1149  if(Validation) {ClusterChargeOverPath/=(*gainused)[i];}
1150  if(OldGainRemoving){ClusterChargeOverPath*=(*gainused)[i];}
1151  }
1152 
1153  // keep pixel cluster charge processing until here
1154  if(APV->SubDet<=2) continue;
1155 
1156  (Charge_Vs_Index[elepos]) ->Fill(APV->Index,ClusterChargeOverPath);
1157 
1158 
1159  // Compute the charge for monitoring and fill the relative histograms
1160  int mCharge1 = 0;
1161  int mCharge2 = 0;
1162  int mCharge3 = 0;
1163  int mCharge4 = 0;
1164  if(APV->SubDet>2) {
1165  for(unsigned int s=0;s<(*nstrips)[i];s++){
1166  int StripCharge = (*amplitude)[FirstAmplitude-(*nstrips)[i]+s];
1167  if(StripCharge>1024) StripCharge = 255;
1168  else if(StripCharge>254) StripCharge = 254;
1169  mCharge1 += StripCharge;
1170  mCharge2 += StripCharge;
1171  mCharge3 += StripCharge;
1172  mCharge4 += StripCharge;
1173  }
1174  // Revome gains for monitoring
1175  mCharge2 *= (*gainused)[i]; // remove G2
1176  mCharge3 *= (*gainusedTick)[i]; // remove G1
1177  mCharge4 *= ( (*gainused)[i] * (*gainusedTick)[i]); // remove G1 and G2
1178  }
1179  std::vector<APVGain::APVmon>& v1 = Charge_1[elepos];
1180  std::vector<MonitorElement*> cmon1 = APVGain::FetchMonitor(v1, (*rawid)[i], tTopo_);
1181  for(unsigned int m=0; m<cmon1.size(); m++) cmon1[m]->Fill(( (double) mCharge1 )/(*path)[i]);
1182 
1183  std::vector<APVGain::APVmon>& v2 = Charge_2[elepos];
1184  std::vector<MonitorElement*> cmon2 = APVGain::FetchMonitor(v2, (*rawid)[i], tTopo_);
1185  for(unsigned int m=0; m<cmon2.size(); m++) cmon2[m]->Fill(( (double) mCharge2 )/(*path)[i]);
1186 
1187  std::vector<APVGain::APVmon>& v3 = Charge_3[elepos];
1188  std::vector<MonitorElement*> cmon3 = APVGain::FetchMonitor(v3, (*rawid)[i], tTopo_);
1189  for(unsigned int m=0; m<cmon3.size(); m++) cmon3[m]->Fill(( (double) mCharge3 )/(*path)[i]);
1190 
1191  std::vector<APVGain::APVmon>& v4 = Charge_4[elepos];
1192  std::vector<MonitorElement*> cmon4 = APVGain::FetchMonitor(v4, (*rawid)[i], tTopo_);
1193  for(unsigned int m=0; m<cmon4.size(); m++) cmon4[m]->Fill(( (double) mCharge4 )/(*path)[i]);
1194 
1195 
1196 
1197  // Fill Charge Vs pathLenght histograms
1198  if(APV->SubDet==StripSubdetector::TIB){
1199  (Charge_Vs_PathlengthTIB[elepos])->Fill((*path)[i],Charge); // TIB
1200 
1201  }else if(APV->SubDet==StripSubdetector::TOB){
1202  (Charge_Vs_PathlengthTOB[elepos])->Fill((*path)[i],Charge); // TOB
1203 
1204  }else if(APV->SubDet==StripSubdetector::TID){
1205  if(APV->Eta<0) { (Charge_Vs_PathlengthTIDM[elepos])->Fill((*path)[i],Charge); } // TID minus
1206  else if(APV->Eta>0){ (Charge_Vs_PathlengthTIDP[elepos])->Fill((*path)[i],Charge); } // TID plus
1207 
1208  }else if(APV->SubDet==StripSubdetector::TEC){
1209  if(APV->Eta<0){
1210  if(APV->Thickness<0.04) { (Charge_Vs_PathlengthTECM1[elepos])->Fill((*path)[i],Charge); } // TEC minus, type 1
1211  else if(APV->Thickness>0.04){ (Charge_Vs_PathlengthTECM2[elepos])->Fill((*path)[i],Charge); } // TEC minus, type 2
1212  } else if(APV->Eta>0){
1213  if(APV->Thickness<0.04) { (Charge_Vs_PathlengthTECP1[elepos])->Fill((*path)[i],Charge); } // TEC plus, type 1
1214  else if(APV->Thickness>0.04){ (Charge_Vs_PathlengthTECP2[elepos])->Fill((*path)[i],Charge); } // TEC plus, type 2
1215  }
1216  }
1217 
1218  }// END OF ON-CLUSTER LOOP
1219 }//END OF processEvent()
1220 
1221 void SiStripGainFromCalibTree::algoAnalyzeTheTree()
1222 {
1223  for(unsigned int i=0;i<VInputFiles.size();i++){
1224  printf("Openning file %3i/%3i --> %s\n",i+1, (int)VInputFiles.size(), (char*)(VInputFiles[i].c_str())); fflush(stdout);
1225  TFile *tfile = TFile::Open(VInputFiles[i].c_str());
1226  TString tree_path = TString::Format("gainCalibrationTree%s/tree",m_calibrationMode.c_str());
1227  TTree* tree = dynamic_cast<TTree*> (tfile->Get(tree_path.Data()));
1228 
1229  tree->SetBranchAddress((EventPrefix_ + "event" + EventSuffix_).c_str(), &eventnumber , nullptr);
1230  tree->SetBranchAddress((EventPrefix_ + "run" + EventSuffix_).c_str(), &runnumber , nullptr);
1231  tree->SetBranchAddress((EventPrefix_ + "TrigTech" + EventSuffix_).c_str(), &TrigTech , nullptr);
1232 
1233  tree->SetBranchAddress((TrackPrefix_ + "chi2ndof" + TrackSuffix_).c_str(), &trackchi2ndof , nullptr);
1234  tree->SetBranchAddress((TrackPrefix_ + "momentum" + TrackSuffix_).c_str(), &trackp , nullptr);
1235  tree->SetBranchAddress((TrackPrefix_ + "pt" + TrackSuffix_).c_str(), &trackpt , nullptr);
1236  tree->SetBranchAddress((TrackPrefix_ + "eta" + TrackSuffix_).c_str(), &tracketa , nullptr);
1237  tree->SetBranchAddress((TrackPrefix_ + "phi" + TrackSuffix_).c_str(), &trackphi , nullptr);
1238  tree->SetBranchAddress((TrackPrefix_ + "hitsvalid" + TrackSuffix_).c_str(), &trackhitsvalid, nullptr);
1239  tree->SetBranchAddress((TrackPrefix_ + "algo" + TrackSuffix_).c_str(), &trackalgo , nullptr);
1240 
1241  tree->SetBranchAddress((CalibPrefix_ + "trackindex" + CalibSuffix_).c_str(), &trackindex , nullptr);
1242  tree->SetBranchAddress((CalibPrefix_ + "rawid" + CalibSuffix_).c_str(), &rawid , nullptr);
1243  tree->SetBranchAddress((CalibPrefix_ + "localdirx" + CalibSuffix_).c_str(), &localdirx , nullptr);
1244  tree->SetBranchAddress((CalibPrefix_ + "localdiry" + CalibSuffix_).c_str(), &localdiry , nullptr);
1245  tree->SetBranchAddress((CalibPrefix_ + "localdirz" + CalibSuffix_).c_str(), &localdirz , nullptr);
1246  tree->SetBranchAddress((CalibPrefix_ + "firststrip" + CalibSuffix_).c_str(), &firststrip , nullptr);
1247  tree->SetBranchAddress((CalibPrefix_ + "nstrips" + CalibSuffix_).c_str(), &nstrips , nullptr);
1248  tree->SetBranchAddress((CalibPrefix_ + "saturation" + CalibSuffix_).c_str(), &saturation , nullptr);
1249  tree->SetBranchAddress((CalibPrefix_ + "overlapping" + CalibSuffix_).c_str(), &overlapping , nullptr);
1250  tree->SetBranchAddress((CalibPrefix_ + "farfromedge" + CalibSuffix_).c_str(), &farfromedge , nullptr);
1251  tree->SetBranchAddress((CalibPrefix_ + "charge" + CalibSuffix_).c_str(), &charge , nullptr);
1252  tree->SetBranchAddress((CalibPrefix_ + "path" + CalibSuffix_).c_str(), &path , nullptr);
1253  tree->SetBranchAddress((CalibPrefix_ + "chargeoverpath" + CalibSuffix_).c_str(), &chargeoverpath, nullptr);
1254  tree->SetBranchAddress((CalibPrefix_ + "amplitude" + CalibSuffix_).c_str(), &amplitude , nullptr);
1255  tree->SetBranchAddress((CalibPrefix_ + "gainused" + CalibSuffix_).c_str(), &gainused , nullptr);
1256  tree->SetBranchAddress((CalibPrefix_ + "gainusedTick" + CalibSuffix_).c_str(), &gainusedTick , nullptr);
1257 
1258 
1259  unsigned int nentries = tree->GetEntries();
1260  printf("Number of Events = %i + %i = %i\n",NEvent,nentries,(NEvent+nentries));
1261  printf("Progressing Bar :0%% 20%% 40%% 60%% 80%% 100%%\n");
1262  printf("Looping on the Tree :");
1263  int TreeStep = nentries/50;if(TreeStep<=1)TreeStep=1;
1264  for (unsigned int ientry = 0; ientry < tree->GetEntries(); ientry++) {
1265  if(ientry%TreeStep==0){printf(".");fflush(stdout);}
1266  tree->GetEntry(ientry);
1267  processEvent();
1268  }printf("\n");// END OF EVENT LOOP
1269  }
1270 }
1271 
1272 
1273 
1274 void SiStripGainFromCalibTree::algoComputeMPVandGain() {
1275  unsigned int I=0;
1276  TH1F* Proj = nullptr;
1277  double FitResults[6];
1278  double MPVmean = 300;
1279 
1280  int elepos = (AlgoMode == "PCL")? Harvest : statCollectionFromMode(m_calibrationMode.c_str());
1281 
1282  if ( Charge_Vs_Index[elepos]==nullptr ) {
1283  edm::LogError("SiStripGainFromCalibTree") << "Harvesting: could not execute algoComputeMPVandGain method because "
1284  << m_calibrationMode.c_str() << " statistics cannot be retrieved.\n"
1285  << "Please check if input contains "
1286  << m_calibrationMode.c_str() << " data." << std::endl;
1287  return;
1288  }
1289 
1290  TH2S *chvsidx = (Charge_Vs_Index[elepos])->getTH2S();
1291 
1292 
1293  printf("Progressing Bar :0%% 20%% 40%% 60%% 80%% 100%%\n");
1294  printf("Fitting Charge Distribution :");
1295  int TreeStep = APVsColl.size()/50;
1296  for(auto it = APVsColl.begin();it!=APVsColl.end();it++,I++){
1297  if(I%TreeStep==0){printf(".");fflush(stdout);}
1298  stAPVGain* APV = it->second;
1299  if(APV->Bin<0) APV->Bin = chvsidx->GetXaxis()->FindBin(APV->Index);
1300 
1301  if(APV->isMasked){APV->Gain=APV->PreviousGain; MASKED++; continue;}
1302 
1303  Proj = (TH1F*)(chvsidx->ProjectionY("",chvsidx->GetXaxis()->FindBin(APV->Index),chvsidx->GetXaxis()->FindBin(APV->Index),"e"));
1304  if(!Proj)continue;
1305 
1306  if(CalibrationLevel==0){
1307  }else if(CalibrationLevel==1){
1308  int SecondAPVId = APV->APVId;
1309  if(SecondAPVId%2==0){ SecondAPVId = SecondAPVId+1; }else{ SecondAPVId = SecondAPVId-1; }
1310  stAPVGain* APV2 = APVsColl[(APV->DetId<<4) | SecondAPVId];
1311  if(APV2->Bin<0) APV2->Bin = chvsidx->GetXaxis()->FindBin(APV2->Index);
1312  TH1F* Proj2 = (TH1F*)(chvsidx->ProjectionY("",APV2->Bin,APV2->Bin,"e"));
1313  if(Proj2){Proj->Add(Proj2,1);delete Proj2;}
1314  }else if(CalibrationLevel==2){
1315  for(unsigned int i=0;i<16;i++){ //loop up to 6APV for Strip and up to 16 for Pixels
1316  auto tmpit = APVsColl.find((APV->DetId<<4) | i);
1317  if(tmpit==APVsColl.end())continue;
1318  stAPVGain* APV2 = tmpit->second;
1319  if(APV2->DetId != APV->DetId || APV2->APVId == APV->APVId)continue;
1320  if(APV2->Bin<0) APV2->Bin = chvsidx->GetXaxis()->FindBin(APV2->Index);
1321  TH1F* Proj2 = (TH1F*)(chvsidx->ProjectionY("",APV2->Bin,APV2->Bin,"e"));
1322  if(Proj2){Proj->Add(Proj2,1);delete Proj2;}
1323  }
1324  }else{
1325  CalibrationLevel = 0;
1326  printf("Unknown Calibration Level, will assume %i\n",CalibrationLevel);
1327  }
1328 
1329  getPeakOfLandau(Proj,FitResults);
1330  APV->FitMPV = FitResults[0];
1331  APV->FitMPVErr = FitResults[1];
1332  APV->FitWidth = FitResults[2];
1333  APV->FitWidthErr = FitResults[3];
1334  APV->FitChi2 = FitResults[4];
1335  APV->FitNorm = FitResults[5];
1336  APV->NEntries = Proj->GetEntries();
1337 
1338  if(IsGoodLandauFit(FitResults)){
1339  APV->Gain = APV->FitMPV / MPVmean;
1340  if(APV->SubDet>2)GOOD++;
1341  }else{
1342  APV->Gain = APV->PreviousGain;
1343  if(APV->SubDet>2)BAD++;
1344  }
1345  if(APV->Gain<=0) APV->Gain = 1;
1346 
1347  //printf("%5i/%5i: %6i - %1i %5E Entries --> MPV = %f +- %f\n",I,APVsColl.size(),APV->DetId, APV->APVId, Proj->GetEntries(), FitResults[0], FitResults[1]);fflush(stdout);
1348  delete Proj;
1349  }printf("\n");
1350 }
1351 
1352 void SiStripGainFromCalibTree::qualityMonitor() {
1353 
1354  int elepos = (AlgoMode == "PCL")? Harvest : statCollectionFromMode(m_calibrationMode.c_str());
1355 
1356  for(unsigned int a=0;a<APVsCollOrdered.size();a++) {
1357  stAPVGain* APV = APVsCollOrdered[a];
1358  if(APV==nullptr)continue;
1359 
1360  unsigned int Index = APV->Index;
1361  unsigned int SubDet = APV->SubDet;
1362  unsigned int DetId = APV->DetId;
1363  float z = APV->z;
1364  float Eta = APV->Eta;
1365  float R = APV->R;
1366  float Phi = APV->Phi;
1367  float Thickness = APV->Thickness;
1368  double FitMPV = APV->FitMPV;
1369  double FitMPVErr = APV->FitMPVErr;
1370  double Gain = APV->Gain;
1371  double NEntries = APV->NEntries;
1372  double PreviousGain = APV->PreviousGain;
1373 
1374  if (SubDet<3) continue; // avoid to loop over Pixel det id
1375 
1376  if (Gain!=1.) {
1377  std::vector<MonitorElement*> charge_histos = APVGain::FetchMonitor(newCharge, DetId, tTopo_);
1378  TH2S *chvsidx = (Charge_Vs_Index[elepos])->getTH2S();
1379  int bin = chvsidx->GetXaxis()->FindBin(Index);
1380  TH1D* Proj = chvsidx->ProjectionY("proj",bin,bin);
1381  for (int binId=0; binId<Proj->GetXaxis()->GetNbins();binId++) {
1382  double new_charge = Proj->GetXaxis()->GetBinCenter(binId) / Gain;
1383  if (Proj->GetBinContent(binId)!=0.) {
1384  for (unsigned int h=0;h<charge_histos.size();h++) {
1385  TH1D* chisto = (charge_histos[h])->getTH1D();
1386  for (int e=0;e<Proj->GetBinContent(binId);e++) chisto->Fill(new_charge);
1387  }
1388  }
1389  }
1390  }
1391 
1392 
1393  if (FitMPV<=0.) { // No fit of MPV
1394  if (APV->isMasked) NoMPVmasked->Fill(z,R);
1395  else NoMPVfit->Fill(z,R);
1396 
1397  } else { // Fit of MPV
1398  if(FitMPV>0.) Gains->Fill(Gain);
1399 
1400  MPVs->Fill(FitMPV);
1401  if(Thickness<0.04) MPVs320->Fill(FitMPV);
1402  if(Thickness>0.04) MPVs500->Fill(FitMPV);
1403 
1404  MPVError->Fill(FitMPVErr);
1405  MPVErrorVsMPV->Fill(FitMPV,FitMPVErr);
1406  MPVErrorVsEta->Fill(Eta,FitMPVErr);
1407  MPVErrorVsPhi->Fill(Phi,FitMPVErr);
1408  MPVErrorVsN->Fill(NEntries,FitMPVErr);
1409 
1410  if(SubDet==3) {
1411  MPV_Vs_EtaTIB->Fill(Eta,FitMPV);
1412  MPV_Vs_PhiTIB->Fill(Phi,FitMPV);
1413  MPVsTIB->Fill(FitMPV);
1414 
1415  } else if(SubDet==4) {
1416  MPV_Vs_EtaTID->Fill(Eta,FitMPV);
1417  MPV_Vs_PhiTID->Fill(Phi,FitMPV);
1418  MPVsTID->Fill(FitMPV);
1419  if(Eta<0.) MPVsTIDM->Fill(FitMPV);
1420  if(Eta>0.) MPVsTIDP->Fill(FitMPV);
1421 
1422  } else if (SubDet==5) {
1423  MPV_Vs_EtaTOB->Fill(Eta,FitMPV);
1424  MPV_Vs_PhiTOB->Fill(Phi,FitMPV);
1425  MPVsTOB->Fill(FitMPV);
1426 
1427  } else if (SubDet==6) {
1428  MPV_Vs_EtaTEC->Fill(Eta,FitMPV);
1429  MPV_Vs_PhiTEC->Fill(Phi,FitMPV);
1430  MPVsTEC->Fill(FitMPV);
1431  if(Eta<0.) MPVsTECM->Fill(FitMPV);
1432  if(Eta>0.) MPVsTECP->Fill(FitMPV);
1433  if(Thickness<0.04) {
1434  MPV_Vs_EtaTECthin->Fill(Eta,FitMPV);
1435  MPV_Vs_PhiTECthin->Fill(Phi,FitMPV);
1436  MPVsTECthin->Fill(FitMPV);
1437  if(Eta>0.) MPVsTECP1->Fill(FitMPV);
1438  if(Eta<0.) MPVsTECM1->Fill(FitMPV);
1439  }
1440  if(Thickness>0.04) {
1441  MPV_Vs_EtaTECthick->Fill(Eta,FitMPV);
1442  MPV_Vs_PhiTECthick->Fill(Phi,FitMPV);
1443  MPVsTECthick->Fill(FitMPV);
1444  if(Eta>0.) MPVsTECP2->Fill(FitMPV);
1445  if(Eta<0.) MPVsTECM2->Fill(FitMPV);
1446  }
1447  }
1448  }
1449 
1450  if(SubDet==3 && PreviousGain!=0. ) DiffWRTPrevGainTIB->Fill(Gain/PreviousGain);
1451  else if(SubDet==4 && PreviousGain!=0. ) DiffWRTPrevGainTID->Fill(Gain/PreviousGain);
1452  else if(SubDet==5 && PreviousGain!=0. ) DiffWRTPrevGainTOB->Fill(Gain/PreviousGain);
1453  else if(SubDet==6 && PreviousGain!=0. ) DiffWRTPrevGainTEC->Fill(Gain/PreviousGain);
1454 
1455  if(SubDet==3 ) GainVsPrevGainTIB->Fill(PreviousGain,Gain);
1456  else if(SubDet==4 ) GainVsPrevGainTID->Fill(PreviousGain,Gain);
1457  else if(SubDet==5 ) GainVsPrevGainTOB->Fill(PreviousGain,Gain);
1458  else if(SubDet==6 ) GainVsPrevGainTEC->Fill(PreviousGain,Gain);
1459  }
1460 }
1461 
1462 void SiStripGainFromCalibTree::storeOnTree(TFileService* tfs)
1463 {
1464  unsigned int tree_Index;
1465  unsigned int tree_Bin;
1466  unsigned int tree_DetId;
1467  unsigned char tree_APVId;
1468  unsigned char tree_SubDet;
1469  float tree_x;
1470  float tree_y;
1471  float tree_z;
1472  float tree_Eta;
1473  float tree_R;
1474  float tree_Phi;
1475  float tree_Thickness;
1476  float tree_FitMPV;
1477  float tree_FitMPVErr;
1478  float tree_FitWidth;
1479  float tree_FitWidthErr;
1480  float tree_FitChi2NDF;
1481  float tree_FitNorm;
1482  double tree_Gain;
1483  double tree_PrevGain;
1484  double tree_PrevGainTick;
1485  double tree_NEntries;
1486  bool tree_isMasked;
1487 
1488  TTree* MyTree;
1489  MyTree = tfs->make<TTree> ("APVGain","APVGain");
1490  MyTree->Branch("Index" ,&tree_Index ,"Index/i");
1491  MyTree->Branch("Bin" ,&tree_Bin ,"Bin/i");
1492  MyTree->Branch("DetId" ,&tree_DetId ,"DetId/i");
1493  MyTree->Branch("APVId" ,&tree_APVId ,"APVId/b");
1494  MyTree->Branch("SubDet" ,&tree_SubDet ,"SubDet/b");
1495  MyTree->Branch("x" ,&tree_x ,"x/F");
1496  MyTree->Branch("y" ,&tree_y ,"y/F");
1497  MyTree->Branch("z" ,&tree_z ,"z/F");
1498  MyTree->Branch("Eta" ,&tree_Eta ,"Eta/F");
1499  MyTree->Branch("R" ,&tree_R ,"R/F");
1500  MyTree->Branch("Phi" ,&tree_Phi ,"Phi/F");
1501  MyTree->Branch("Thickness" ,&tree_Thickness ,"Thickness/F");
1502  MyTree->Branch("FitMPV" ,&tree_FitMPV ,"FitMPV/F");
1503  MyTree->Branch("FitMPVErr" ,&tree_FitMPVErr ,"FitMPVErr/F");
1504  MyTree->Branch("FitWidth" ,&tree_FitWidth ,"FitWidth/F");
1505  MyTree->Branch("FitWidthErr" ,&tree_FitWidthErr,"FitWidthErr/F");
1506  MyTree->Branch("FitChi2NDF" ,&tree_FitChi2NDF ,"FitChi2NDF/F");
1507  MyTree->Branch("FitNorm" ,&tree_FitNorm ,"FitNorm/F");
1508  MyTree->Branch("Gain" ,&tree_Gain ,"Gain/D");
1509  MyTree->Branch("PrevGain" ,&tree_PrevGain ,"PrevGain/D");
1510  MyTree->Branch("PrevGainTick" ,&tree_PrevGainTick,"PrevGainTick/D");
1511  MyTree->Branch("NEntries" ,&tree_NEntries ,"NEntries/D");
1512  MyTree->Branch("isMasked" ,&tree_isMasked ,"isMasked/O");
1513 
1514 
1515  FILE* Gains = stdout;
1516  fprintf(Gains,"NEvents = %i\n",NEvent);
1517  fprintf(Gains,"NTracks = %i\n",NTrack);
1518  //fprintf(Gains,"NClustersPixel = %i\n",NClusterPixel);
1519  fprintf(Gains,"NClustersStrip = %i\n",NClusterStrip);
1520  //fprintf(Gains,"Number of Pixel Dets = %lu\n",static_cast<unsigned long>(NPixelDets));
1521  fprintf(Gains,"Number of Strip APVs = %lu\n",static_cast<unsigned long>(NStripAPVs));
1522  fprintf(Gains,"GoodFits = %i BadFits = %i ratio = %f%% (MASKED=%i)\n",GOOD,BAD,(100.0*GOOD)/(GOOD+BAD), MASKED);
1523 
1524  Gains=fopen(OutputGains.c_str(),"w");
1525  fprintf(Gains,"NEvents = %i\n",NEvent);
1526  fprintf(Gains,"NTracks = %i\n",NTrack);
1527  //fprintf(Gains,"NClustersPixel = %i\n",NClusterPixel);
1528  fprintf(Gains,"NClustersStrip = %i\n",NClusterStrip);
1529  fprintf(Gains,"Number of Strip APVs = %lu\n",static_cast<unsigned long>(NStripAPVs));
1530  //fprintf(Gains,"Number of Pixel Dets = %lu\n",static_cast<unsigned long>(NPixelDets));
1531  fprintf(Gains,"GoodFits = %i BadFits = %i ratio = %f%% (MASKED=%i)\n",GOOD,BAD,(100.0*GOOD)/(GOOD+BAD), MASKED);
1532 
1533  int elepos = statCollectionFromMode(m_calibrationMode.c_str());
1534 
1535  for(unsigned int a=0;a<APVsCollOrdered.size();a++){
1536  stAPVGain* APV = APVsCollOrdered[a];
1537  if(APV==nullptr)continue;
1538 // printf( "%i | %i | PreviousGain = %7.5f NewGain = %7.5f (#clusters=%8.0f)\n", APV->DetId,APV->APVId,APV->PreviousGain,APV->Gain, APV->NEntries);
1539  fprintf(Gains,"%i | %i | PreviousGain = %7.5f(tick) x %7.5f(particle) NewGain (particle) = %7.5f (#clusters=%8.0f)\n", APV->DetId,APV->APVId,APV->PreviousGainTick, APV->PreviousGain,APV->Gain, APV->NEntries);
1540 
1541  tree_Index = APV->Index;
1542  tree_Bin = (Charge_Vs_Index[elepos])->getTH2S()->GetXaxis()->FindBin(APV->Index);
1543  tree_DetId = APV->DetId;
1544  tree_APVId = APV->APVId;
1545  tree_SubDet = APV->SubDet;
1546  tree_x = APV->x;
1547  tree_y = APV->y;
1548  tree_z = APV->z;
1549  tree_Eta = APV->Eta;
1550  tree_R = APV->R;
1551  tree_Phi = APV->Phi;
1552  tree_Thickness = APV->Thickness;
1553  tree_FitMPV = APV->FitMPV;
1554  tree_FitMPVErr = APV->FitMPVErr;
1555  tree_FitWidth = APV->FitWidth;
1556  tree_FitWidthErr= APV->FitWidthErr;
1557  tree_FitChi2NDF = APV->FitChi2;
1558  tree_FitNorm = APV->FitNorm;
1559  tree_Gain = APV->Gain;
1560  tree_PrevGain = APV->PreviousGain;
1561  tree_PrevGainTick = APV->PreviousGainTick;
1562  tree_NEntries = APV->NEntries;
1563  tree_isMasked = APV->isMasked;
1564 
1565 
1566  if(tree_DetId==402673324){
1567  printf("%i | %i : %f --> %f (%f)\n", tree_DetId, tree_APVId, tree_PrevGain, tree_Gain, tree_NEntries);
1568  }
1569 
1570 
1571  MyTree->Fill();
1572  }
1573  if(Gains)fclose(Gains);
1574 
1575 
1576 }
1577 
1578 bool SiStripGainFromCalibTree::produceTagFilter(){
1579 
1580  // The goal of this function is to check wether or not there is enough statistics
1581  // to produce a meaningful tag for the DB
1582  int elepos = (AlgoMode == "PCL")? Harvest : statCollectionFromMode(m_calibrationMode.c_str());
1583  if( Charge_Vs_Index[elepos]==nullptr ) {
1584  edm::LogError("SiStripGainFromCalibTree") << "produceTagFilter -> Return false: could not retrieve the "
1585  << m_calibrationMode.c_str() << " statistics.\n"
1586  << "Please check if input contains "
1587  << m_calibrationMode.c_str() << " data." << std::endl;
1588  return false;
1589  }
1590 
1591 
1592  float integral = (Charge_Vs_Index[elepos])->getTH2S()->Integral();
1593  if( (Charge_Vs_Index[elepos])->getTH2S()->Integral(0,NStripAPVs+1, 0, 99999 ) < tagCondition_NClusters) {
1594  edm::LogWarning("SiStripGainFromCalibTree")
1595  << "calibrationMode -> " << m_calibrationMode << "\n"
1596  << "produceTagFilter -> Return false: Statistics is too low : " << integral << endl;
1597  return false;
1598  }
1599  if((1.0 * GOOD) / (GOOD+BAD) < tagCondition_GoodFrac) {
1600  edm::LogWarning("SiStripGainFromCalibTree")
1601  << "calibrationMode -> " << m_calibrationMode << "\n"
1602  << "produceTagFilter -> Return false: ratio of GOOD/TOTAL is too low: " << (1.0 * GOOD) / (GOOD+BAD) << endl;
1603  return false;
1604  }
1605  return true;
1606 }
1607 
1608 SiStripApvGain* SiStripGainFromCalibTree::getNewObject()
1609 {
1610  SiStripApvGain* obj = new SiStripApvGain();
1611  if(!m_harvestingMode) return obj;
1612 
1613  if(!produceTagFilter()){
1614  edm::LogWarning("SiStripGainFromCalibTree")<< "getNewObject -> will not produce a paylaod because produceTagFilter returned false " << endl;
1615  setDoStore(false);
1616  return obj;
1617  }
1618 
1619 
1620  std::vector<float>* theSiStripVector = nullptr;
1621  unsigned int PreviousDetId = 0;
1622  for(unsigned int a=0;a<APVsCollOrdered.size();a++){
1623  stAPVGain* APV = APVsCollOrdered[a];
1624  if(APV==nullptr){ printf("Bug\n"); continue; }
1625  if(APV->SubDet<=2)continue;
1626  if(APV->DetId != PreviousDetId){
1627  if(theSiStripVector!=nullptr){
1628  SiStripApvGain::Range range(theSiStripVector->begin(),theSiStripVector->end());
1629  if ( !obj->put(PreviousDetId,range) ) printf("Bug to put detId = %i\n",PreviousDetId);
1630  }
1631  theSiStripVector = new std::vector<float>;
1632  PreviousDetId = APV->DetId;
1633  }
1634  theSiStripVector->push_back(APV->Gain);
1635  }
1636  if(theSiStripVector!=nullptr){
1637  SiStripApvGain::Range range(theSiStripVector->begin(),theSiStripVector->end());
1638  if ( !obj->put(PreviousDetId,range) ) printf("Bug to put detId = %i\n",PreviousDetId);
1639  }
1640 
1641  if (theSiStripVector!=nullptr) delete theSiStripVector;
1642 
1643  return obj;
1644 }
1645 
1646 
1647 SiStripGainFromCalibTree::~SiStripGainFromCalibTree()
1648 {
1649  APVsColl.clear();
1650  for(unsigned int a=0;a<APVsCollOrdered.size();a++) {
1651  stAPVGain* APV = APVsCollOrdered[a];
1652  if(APV!=nullptr) delete APV;
1653  }
1654  APVsCollOrdered.clear();
1655 }
1656 
1657 
1658 void SiStripGainFromCalibTree::MakeCalibrationMap(){
1659  if(!useCalibration)return;
1660 
1661  TChain* t1 = new TChain("SiStripCalib/APVGain");
1662  t1->Add(m_calibrationPath.c_str());
1663 
1664  unsigned int tree_DetId;
1665  unsigned char tree_APVId;
1666  double tree_Gain;
1667 
1668  t1->SetBranchAddress("DetId" ,&tree_DetId );
1669  t1->SetBranchAddress("APVId" ,&tree_APVId );
1670  t1->SetBranchAddress("Gain" ,&tree_Gain );
1671 
1672  for (unsigned int ientry = 0; ientry < t1->GetEntries(); ientry++) {
1673  t1->GetEntry(ientry);
1674  stAPVGain* APV = APVsColl[(tree_DetId<<4) | (unsigned int)tree_APVId];
1675  APV->CalibGain = tree_Gain;
1676  }
1677 
1678  delete t1;
1679 }
1680 
1681 void
1682 SiStripGainFromCalibTree::algoAnalyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
1683 {
1684  // in AlCaHarvesting mode we just need to run the logic in the endJob step
1685  if(m_harvestingMode) return;
1686 
1687  if(AlgoMode=="CalibTree")return;
1688 
1689  eventnumber = iEvent.id().event();
1690  runnumber = iEvent.id().run();
1691  auto handle01 = connect(TrigTech , TrigTech_token_ , iEvent);
1692  auto handle02 = connect(trackchi2ndof , trackchi2ndof_token_ , iEvent);
1693  auto handle03 = connect(trackp , trackp_token_ , iEvent);
1694  auto handle04 = connect(trackpt , trackpt_token_ , iEvent);
1695  auto handle05 = connect(tracketa , tracketa_token_ , iEvent);
1696  auto handle06 = connect(trackphi , trackphi_token_ , iEvent);
1697  auto handle07 = connect(trackhitsvalid, trackhitsvalid_token_, iEvent);
1698  auto handle08 = connect(trackindex , trackindex_token_ , iEvent);
1699  auto handle09 = connect(rawid , rawid_token_ , iEvent);
1700  auto handle11 = connect(localdirx , localdirx_token_ , iEvent);
1701  auto handle12 = connect(localdiry , localdiry_token_ , iEvent);
1702  auto handle13 = connect(localdirz , localdirz_token_ , iEvent);
1703  auto handle14 = connect(firststrip , firststrip_token_ , iEvent);
1704  auto handle15 = connect(nstrips , nstrips_token_ , iEvent);
1705  auto handle16 = connect(saturation , saturation_token_ , iEvent);
1706  auto handle17 = connect(overlapping , overlapping_token_ , iEvent);
1707  auto handle18 = connect(farfromedge , farfromedge_token_ , iEvent);
1708  auto handle19 = connect(charge , charge_token_ , iEvent);
1709  auto handle21 = connect(path , path_token_ , iEvent);
1710  auto handle22 = connect(chargeoverpath, chargeoverpath_token_, iEvent);
1711  auto handle23 = connect(amplitude , amplitude_token_ , iEvent);
1712  auto handle24 = connect(gainused , gainused_token_ , iEvent);
1713  auto handle25 = connect(gainusedTick , gainusedTick_token_ , iEvent);
1714  auto handle26 = connect(trackalgo , trackalgo_token_ , iEvent);
1715 
1716  processEvent();
1717 }
1718 
1719 DEFINE_FWK_MODULE(SiStripGainFromCalibTree);
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
SiStripQuality::IsApvBad
bool IsApvBad(const uint32_t &detid, const short &apvNb) const
Definition: SiStripQuality.cc:617
SiStripGainFromCalibTree::trackindex_token_
edm::EDGetTokenT< std::vector< int > > trackindex_token_
Definition: SiStripGainFromCalibTree.cc:281
SiStripGainFromCalibTree::Charge_Vs_Index
std::vector< MonitorElement * > Charge_Vs_Index
Definition: SiStripGainFromCalibTree.cc:178
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::EventID::event
EventNumber_t event() const
Definition: EventID.h:41
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
PixelTopology::nrows
virtual int nrows() const =0
SiStripGainFromCalibTree::trackp
const std::vector< float > * trackp
Definition: SiStripGainFromCalibTree.cc:273
SiStripGainFromCalibTree::tracketa_token_
edm::EDGetTokenT< std::vector< double > > tracketa_token_
Definition: SiStripGainFromCalibTree.cc:275
SiStripGainFromCalibTree::amplitude_token_
edm::EDGetTokenT< std::vector< unsigned char > > amplitude_token_
Definition: SiStripGainFromCalibTree.cc:294
APVGain::FetchMonitor
std::vector< MonitorElement * > FetchMonitor(std::vector< APVmon >, uint32_t, const TrackerTopology *topo=nullptr)
Definition: APVGainHelpers.cc:108
GlobalPosition_Frontier_DevDB_cff.tag
tag
Definition: GlobalPosition_Frontier_DevDB_cff.py:11
SiStripGainFromCalibTree::newCharge
std::vector< APVGain::APVmon > newCharge
Definition: SiStripGainFromCalibTree.cc:249
SiStripGainFromCalibTree::chargeoverpath
const std::vector< double > * chargeoverpath
Definition: SiStripGainFromCalibTree.cc:293
mps_fire.i
i
Definition: mps_fire.py:269
stAPVGain::FitMPV
double FitMPV
Definition: APVGainStruct.h:19
SiStripGainFromCalibTree::rawid_token_
edm::EDGetTokenT< std::vector< unsigned int > > rawid_token_
Definition: SiStripGainFromCalibTree.cc:282
SiStripGainFromCalibTree::amplitude
const std::vector< unsigned char > * amplitude
Definition: SiStripGainFromCalibTree.cc:294
SiStripGainFromCalibTree::isBFieldConsistentWithMode
bool isBFieldConsistentWithMode(const edm::EventSetup &iSetup) const
Definition: SiStripGainFromCalibTree.cc:772
h
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
Definition: L1TUtmAlgorithmRcd.h:4
stAPVGain::PreviousGainTick
double PreviousGainTick
Definition: APVGainStruct.h:28
SiStripGainFromCalibTree::getNewObject
SiStripApvGain * getNewObject() override
Definition: SiStripGainFromCalibTree.cc:1608
SiStripGainFromCalibTree::localdirz
const std::vector< double > * localdirz
Definition: SiStripGainFromCalibTree.cc:285
SiStripGain::getNumberOfTags
size_t getNumberOfTags() const
Definition: SiStripGain.h:91
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:519
SiStripGainFromCalibTree::trackchi2ndof
const std::vector< double > * trackchi2ndof
Definition: SiStripGainFromCalibTree.cc:272
SiStripGainFromCalibTree::algoBeginRun
void algoBeginRun(const edm::Run &run, const edm::EventSetup &iSetup) override
Definition: SiStripGainFromCalibTree.cc:791
SiStripGainFromCalibTree::localdirx
const std::vector< double > * localdirx
Definition: SiStripGainFromCalibTree.cc:283
SiStripGainFromCalibTree::gainusedTick
const std::vector< double > * gainusedTick
Definition: SiStripGainFromCalibTree.cc:296
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
processEvent
dispatcher processEvent(e, inputTag, standby)
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
stAPVGain::APVId
unsigned int APVId
Definition: APVGainStruct.h:10
stAPVGain::Thickness
float Thickness
Definition: APVGainStruct.h:18
SiStripGainFromCalibTree::algoEndRun
void algoEndRun(const edm::Run &run, const edm::EventSetup &iSetup) override
Definition: SiStripGainFromCalibTree.cc:839
SiStripGainFromCalibTree::rawid
const std::vector< unsigned int > * rawid
Definition: SiStripGainFromCalibTree.cc:282
SiStripGainFromCalibTree::overlapping_token_
edm::EDGetTokenT< std::vector< bool > > overlapping_token_
Definition: SiStripGainFromCalibTree.cc:289
SiStripGainFromCalibTree::SiStripGainFromCalibTree
SiStripGainFromCalibTree(const edm::ParameterSet &)
Definition: SiStripGainFromCalibTree.cc:340
SiStripGainFromCalibTree::tTopo_
const TrackerTopology * tTopo_
Definition: SiStripGainFromCalibTree.cc:310
SiStripGainFromCalibTree::connect
edm::Handle< T > connect(const T *&ptr, edm::EDGetTokenT< T > token, const edm::Event &evt)
Definition: SiStripGainFromCalibTree.cc:127
APVGain::subdetectorPlane
int subdetectorPlane(uint32_t, const TrackerTopology *)
Definition: APVGainHelpers.cc:79
SiStripGainFromCalibTree::trackalgo
const std::vector< int > * trackalgo
Definition: SiStripGainFromCalibTree.cc:278
SiStripGainFromCalibTree::Charge_Vs_PathlengthTOB
std::vector< MonitorElement * > Charge_Vs_PathlengthTOB
Definition: SiStripGainFromCalibTree.cc:185
SiStripGainFromCalibTree::storeOnTree
void storeOnTree(TFileService *tfs)
Definition: SiStripGainFromCalibTree.cc:1462
TFileService::make
T * make(const Args &...args) const
make new ROOT object
Definition: TFileService.h:64
SiStripGainFromCalibTree::trackpt
const std::vector< float > * trackpt
Definition: SiStripGainFromCalibTree.cc:274
SiStripGainFromCalibTree::Charge_Vs_PathlengthTECM1
std::vector< MonitorElement * > Charge_Vs_PathlengthTECM1
Definition: SiStripGainFromCalibTree.cc:190
stAPVGain::DetId
unsigned int DetId
Definition: APVGainStruct.h:9
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
nullptr
#define nullptr
CustomPhysics_cfi.amplitude
amplitude
Definition: CustomPhysics_cfi.py:12
SiStripGainFromCalibTree::trackhitsvalid
const std::vector< unsigned int > * trackhitsvalid
Definition: SiStripGainFromCalibTree.cc:277
SiStripGainFromCalibTree::Charge_Vs_PathlengthTECM2
std::vector< MonitorElement * > Charge_Vs_PathlengthTECM2
Definition: SiStripGainFromCalibTree.cc:191
stAPVGain::isMasked
bool isMasked
Definition: APVGainStruct.h:32
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:168
SiStripGain::getApvGain
static float getApvGain(const uint16_t &apv, const SiStripApvGain::Range &range)
Definition: SiStripGain.h:73
SiStripGainFromCalibTree::farfromedge_token_
edm::EDGetTokenT< std::vector< bool > > farfromedge_token_
Definition: SiStripGainFromCalibTree.cc:290
SiStripApvGain::put
bool put(const uint32_t &detID, Range input)
Definition: SiStripApvGain.cc:8
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:250
SiStripGainFromCalibTree::trackhitsvalid_token_
edm::EDGetTokenT< std::vector< unsigned int > > trackhitsvalid_token_
Definition: SiStripGainFromCalibTree.cc:277
TrackerGeometry::dets
const DetContainer & dets() const override
Returm a vector of all GeomDet (including all GeomDetUnits)
Definition: TrackerGeometry.h:61
SiStripGainFromCalibTree::APVsCollOrdered
std::vector< stAPVGain * > APVsCollOrdered
Definition: SiStripGainFromCalibTree.cc:308
SiStripGainFromCalibTree::path_token_
edm::EDGetTokenT< std::vector< double > > path_token_
Definition: SiStripGainFromCalibTree.cc:292
stAPVGain::FitWidth
double FitWidth
Definition: APVGainStruct.h:21
iEvent
int iEvent
Definition: GenABIO.cc:230
stAPVGain::Gain
double Gain
Definition: APVGainStruct.h:25
SiStripGainFromCalibTree::trackalgo_token_
edm::EDGetTokenT< std::vector< int > > trackalgo_token_
Definition: SiStripGainFromCalibTree.cc:278
fastjetJetProducer_validation_cfg.DQMStore
DQMStore
Definition: fastjetJetProducer_validation_cfg.py:29
stAPVGain::PreviousGain
double PreviousGain
Definition: APVGainStruct.h:27
SiStripGainFromCalibTree::trackp_token_
edm::EDGetTokenT< std::vector< float > > trackp_token_
Definition: SiStripGainFromCalibTree.cc:273
SiStripGainFromCalibTree::IsGoodLandauFit
bool IsGoodLandauFit(double *FitResults)
Definition: SiStripGainFromCalibTree.cc:1065
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
SiStripGainFromCalibTree::APVsColl
std::unordered_map< unsigned int, stAPVGain * > APVsColl
Definition: SiStripGainFromCalibTree.cc:309
stAPVGain::Index
unsigned int Index
Definition: APVGainStruct.h:7
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
tfile
Definition: tfile.py:1
stAPVGain::FitNorm
double FitNorm
Definition: APVGainStruct.h:24
stAPVGain::CalibGain
double CalibGain
Definition: APVGainStruct.h:26
SiStripGainFromCalibTree::Charge_3
std::array< std::vector< APVGain::APVmon >, 7 > Charge_3
Definition: SiStripGainFromCalibTree.cc:181
SiStripGainFromCalibTree::algoBeginJob
void algoBeginJob(const edm::EventSetup &iSetup) override
Definition: SiStripGainFromCalibTree.cc:631
SiStripApvGain::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripApvGain.h:30
SiStripGainFromCalibTree::TrigTech_token_
edm::EDGetTokenT< std::vector< bool > > TrigTech_token_
Definition: SiStripGainFromCalibTree.cc:269
SiStripGainFromCalibTree::chargeoverpath_token_
edm::EDGetTokenT< std::vector< double > > chargeoverpath_token_
Definition: SiStripGainFromCalibTree.cc:293
Exhume::I
const std::complex< double > I
Definition: I.h:8
SiStripGainFromCalibTree::gainusedTick_token_
edm::EDGetTokenT< std::vector< double > > gainusedTick_token_
Definition: SiStripGainFromCalibTree.cc:296
SiStripGainFromCalibTree::gainused
const std::vector< double > * gainused
Definition: SiStripGainFromCalibTree.cc:295
stAPVGain::FitMPVErr
double FitMPVErr
Definition: APVGainStruct.h:20
SiStripGainFromCalibTree::Charge_Vs_PathlengthTECP1
std::vector< MonitorElement * > Charge_Vs_PathlengthTECP1
Definition: SiStripGainFromCalibTree.cc:188
SiStripGainFromCalibTree::localdirz_token_
edm::EDGetTokenT< std::vector< double > > localdirz_token_
Definition: SiStripGainFromCalibTree.cc:285
TtFullHadDaughter::B
static const std::string B
Definition: TtFullHadronicEvent.h:9
SiStripGainFromCalibTree::firststrip
const std::vector< unsigned short > * firststrip
Definition: SiStripGainFromCalibTree.cc:286
ALCARECOTkAlJpsiMuMu_cff.charge
charge
Definition: ALCARECOTkAlJpsiMuMu_cff.py:47
SiStripGainFromCalibTree::Charge_Vs_PathlengthTIDM
std::vector< MonitorElement * > Charge_Vs_PathlengthTIDM
Definition: SiStripGainFromCalibTree.cc:187
SiStripGainFromCalibTree::localdiry
const std::vector< double > * localdiry
Definition: SiStripGainFromCalibTree.cc:284
SiStripGainFromCalibTree::Charge_2
std::array< std::vector< APVGain::APVmon >, 7 > Charge_2
Definition: SiStripGainFromCalibTree.cc:180
SiStripGainFromCalibTree::charge_token_
edm::EDGetTokenT< std::vector< unsigned int > > charge_token_
Definition: SiStripGainFromCalibTree.cc:291
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37
p2
double p2[4]
Definition: TauolaWrapper.h:90
mps_setup.stdout
stdout
Definition: mps_setup.py:242
funct::integral
Integral< F, X >::type integral(const F &f)
Definition: Integral.h:69
HTMLExport.elem
def elem(elemtype, innerHTML='', html_class='', kwargs)
Definition: HTMLExport.py:18
stAPVGain::FitChi2
double FitChi2
Definition: APVGainStruct.h:23
stringResolutionProvider_cfi.bin
bin
set the eta bin as selection string.
Definition: stringResolutionProvider_cfi.py:11
SiStripGainFromCalibTree::nstrips_token_
edm::EDGetTokenT< std::vector< unsigned short > > nstrips_token_
Definition: SiStripGainFromCalibTree.cc:287
SiStripGainFromCalibTree::Charge_Vs_PathlengthTECP2
std::vector< MonitorElement * > Charge_Vs_PathlengthTECP2
Definition: SiStripGainFromCalibTree.cc:189
SiStripGainFromCalibTree::trackindex
const std::vector< int > * trackindex
Definition: SiStripGainFromCalibTree.cc:281
DetId
Definition: DetId.h:18
SiStripGainFromCalibTree::statCollectionFromMode
int statCollectionFromMode(const char *tag) const
Definition: SiStripGainFromCalibTree.cc:314
SiStripGainFromCalibTree::Charge_Vs_PathlengthTIB
std::vector< MonitorElement * > Charge_Vs_PathlengthTIB
Definition: SiStripGainFromCalibTree.cc:184
SiStripGainFromCalibTree::firststrip_token_
edm::EDGetTokenT< std::vector< unsigned short > > firststrip_token_
Definition: SiStripGainFromCalibTree.cc:286
SiStripGainFromCalibTree::algoAnalyze
void algoAnalyze(const edm::Event &, const edm::EventSetup &) override
Definition: SiStripGainFromCalibTree.cc:1682
SiStripGainFromCalibTree::trackpt_token_
edm::EDGetTokenT< std::vector< float > > trackpt_token_
Definition: SiStripGainFromCalibTree.cc:274
edm::Handle::product
T const * product() const
Definition: Handle.h:81
GetRecoTauVFromDQM_MC_cff.obj
obj
Definition: GetRecoTauVFromDQM_MC_cff.py:87
SiStripGainFromCalibTree::overlapping
const std::vector< bool > * overlapping
Definition: SiStripGainFromCalibTree.cc:289
SiStripGainFromCalibTree::saturation
const std::vector< bool > * saturation
Definition: SiStripGainFromCalibTree.cc:288
StripTopology::nstrips
virtual int nstrips() const =0
edm::merge
bool merge(LuminosityBlockRange &lh, LuminosityBlockRange &rh)
Definition: LuminosityBlockRange.cc:78
SiStripGainFromCalibTree::getPeakOfLandau
void getPeakOfLandau(TH1 *InputHisto, double *FitResults, double LowRange=50, double HighRange=5400)
Definition: SiStripGainFromCalibTree.cc:1038
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:58
APVGain::subdetectorSide
int subdetectorSide(uint32_t, const TrackerTopology *)
Definition: APVGainHelpers.cc:49
b
double b
Definition: hdecay.h:120
stAPVGain::NEntries
double NEntries
Definition: APVGainStruct.h:29
SiStripGainFromCalibTree::farfromedge
const std::vector< bool > * farfromedge
Definition: SiStripGainFromCalibTree.cc:290
RunInfo::m_avg_current
float m_avg_current
Definition: RunInfo.h:29
SiStripGainFromCalibTree::TrigTech
const std::vector< bool > * TrigTech
Definition: SiStripGainFromCalibTree.cc:269
ConditionDBWriter< SiStripApvGain >::setDoStore
void setDoStore(const bool doStore)
When set to false the payload will not be written to the db.
Definition: ConditionDBWriter.h:379
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
PixelTopology::ncolumns
virtual int ncolumns() const =0
p1
double p1[4]
Definition: TauolaWrapper.h:89
a
double a
Definition: hdecay.h:121
SiStripGainFromCalibTree::VChargeHisto
std::vector< std::string > VChargeHisto
Definition: SiStripGainFromCalibTree.cc:164
SiStripGainFromCalibTree::localdirx_token_
edm::EDGetTokenT< std::vector< double > > localdirx_token_
Definition: SiStripGainFromCalibTree.cc:283
SiStripGainFromCalibTree::saturation_token_
edm::EDGetTokenT< std::vector< bool > > saturation_token_
Definition: SiStripGainFromCalibTree.cc:288
patCaloMETCorrections_cff.A
A
Definition: patCaloMETCorrections_cff.py:43
SiStripGainFromCalibTree::trackchi2ndof_token_
edm::EDGetTokenT< std::vector< double > > trackchi2ndof_token_
Definition: SiStripGainFromCalibTree.cc:272
SiStripGainFromCalibTree::Charge_Vs_PathlengthTIDP
std::vector< MonitorElement * > Charge_Vs_PathlengthTIDP
Definition: SiStripGainFromCalibTree.cc:186
SiStripGainFromCalibTree::localdiry_token_
edm::EDGetTokenT< std::vector< double > > localdiry_token_
Definition: SiStripGainFromCalibTree.cc:284
tree
Definition: tree.py:1
JetChargeProducer_cfi.exp
exp
Definition: JetChargeProducer_cfi.py:6
SiStripGainFromCalibTree::charge
const std::vector< unsigned int > * charge
Definition: SiStripGainFromCalibTree.cc:291
APVGain::monHnames
std::vector< std::pair< std::string, std::string > > monHnames(std::vector< std::string >, bool, const char *tag)
Definition: APVGainHelpers.cc:137
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:47
g4SimHits_cfi.Gain
Gain
Definition: g4SimHits_cfi.py:33
SiStripGainFromCalibTree::path
const std::vector< double > * path
Definition: SiStripGainFromCalibTree.cc:292
T
long double T
Definition: Basic3DVectorLD.h:66
SiStripGainFromCalibTree::trackphi
const std::vector< double > * trackphi
Definition: SiStripGainFromCalibTree.cc:276
SiStripGainFromCalibTree::Charge_4
std::array< std::vector< APVGain::APVmon >, 7 > Charge_4
Definition: SiStripGainFromCalibTree.cc:182
stAPVGain::SubDet
unsigned int SubDet
Definition: APVGainStruct.h:11
SiStripGainFromCalibTree::nstrips
const std::vector< unsigned short > * nstrips
Definition: SiStripGainFromCalibTree.cc:287
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
SiStripGainFromCalibTree::gainused_token_
edm::EDGetTokenT< std::vector< double > > gainused_token_
Definition: SiStripGainFromCalibTree.cc:295
APVGain::subdetectorId
int subdetectorId(uint32_t)
Definition: APVGainHelpers.cc:17
SiStripGainFromCalibTree::bookDQMHistos
void bookDQMHistos(const char *dqm_dir, const char *tag)
Definition: SiStripGainFromCalibTree.cc:457
SiStripGainFromCalibTree::trackphi_token_
edm::EDGetTokenT< std::vector< double > > trackphi_token_
Definition: SiStripGainFromCalibTree.cc:276
SiStripGainFromCalibTree::tracketa
const std::vector< double > * tracketa
Definition: SiStripGainFromCalibTree.cc:275
edm::Run
Definition: Run.h:43
stAPVGain::FitWidthErr
double FitWidthErr
Definition: APVGainStruct.h:22
SiStripGain::getRange
const SiStripApvGain::Range getRange(uint32_t detID) const
Definition: SiStripGain.h:70
SiStripGainFromCalibTree::Charge_1
std::array< std::vector< APVGain::APVmon >, 7 > Charge_1
Definition: SiStripGainFromCalibTree.cc:179
Generated for CMSSW Reference Manual by doxygen 1.8.11