CMS 3D CMS Logo
CalibrationAnalysis.cc
Go to the documentation of this file.
1 #include "CondFormats/SiStripObjects/interface/CalibrationAnalysis.h"
2 #include "DataFormats/SiStripCommon/interface/SiStripHistoTitle.h"
3 #include "DataFormats/SiStripCommon/interface/SiStripEnumsAndStrings.h"
4 #include "FWCore/MessageLogger/interface/MessageLogger.h"
5 #include <iostream>
6 #include <sstream>
7 #include <iomanip>
8 #include <cmath>
9 
10 using namespace sistrip;
11 
12 const float CalibrationAnalysis::minAmplitudeThreshold_ = 50;
13 const float CalibrationAnalysis::minBaselineThreshold_ = -50;
14 const float CalibrationAnalysis::maxBaselineThreshold_ = 50;
15 const float CalibrationAnalysis::maxChi2Threshold_ = 3;
16 const float CalibrationAnalysis::minDecayTimeThreshold_ = 30;
17 const float CalibrationAnalysis::maxDecayTimeThreshold_ = 250;
18 const float CalibrationAnalysis::minPeakTimeThreshold_ = 40;
19 const float CalibrationAnalysis::maxPeakTimeThreshold_ = 130;
20 const float CalibrationAnalysis::minRiseTimeThreshold_ = 10;
21 const float CalibrationAnalysis::maxRiseTimeThreshold_ = 100;
22 const float CalibrationAnalysis::minTurnOnThreshold_ = 5;
23 const float CalibrationAnalysis::maxTurnOnThreshold_ = 40;
24 
25 const float CalibrationAnalysis::minDecayTimeThresholdDeco_ = 10;
26 const float CalibrationAnalysis::maxDecayTimeThresholdDeco_ = 100;
27 const float CalibrationAnalysis::minPeakTimeThresholdDeco_ = 40;
28 const float CalibrationAnalysis::maxPeakTimeThresholdDeco_ = 130;
29 const float CalibrationAnalysis::minRiseTimeThresholdDeco_ = 10;
30 const float CalibrationAnalysis::maxRiseTimeThresholdDeco_ = 100;
31 const float CalibrationAnalysis::minTurnOnThresholdDeco_ = 10;
32 const float CalibrationAnalysis::maxTurnOnThresholdDeco_ = 80;
33 
34 CalibrationAnalysis::CalibrationAnalysis(const uint32_t& key, const bool& deconv)
35  : CommissioningAnalysis(key, "CalibrationAnalysis"),
36  amplitude_(2, VFloat(128, sistrip::invalid_)),
37  tail_(2, VFloat(128, sistrip::invalid_)),
38  riseTime_(2, VFloat(128, sistrip::invalid_)),
39  decayTime_(2, VFloat(128, sistrip::invalid_)),
40  turnOn_(2, VFloat(128, sistrip::invalid_)),
41  peakTime_(2, VFloat(128, sistrip::invalid_)),
42  undershoot_(2, VFloat(128, sistrip::invalid_)),
43  baseline_(2, VFloat(128, sistrip::invalid_)),
44  smearing_(2, VFloat(128, sistrip::invalid_)),
45  chi2_(2, VFloat(128, sistrip::invalid_)),
46  isvalid_(2, VBool(128, sistrip::invalid_)),
47  mean_amplitude_(2, sistrip::invalid_),
48  mean_tail_(2, sistrip::invalid_),
49  mean_riseTime_(2, sistrip::invalid_),
50  mean_decayTime_(2, sistrip::invalid_),
51  mean_turnOn_(2, sistrip::invalid_),
52  mean_peakTime_(2, sistrip::invalid_),
53  mean_undershoot_(2, sistrip::invalid_),
54  mean_baseline_(2, sistrip::invalid_),
55  mean_smearing_(2, sistrip::invalid_),
56  mean_chi2_(2, sistrip::invalid_),
57  min_amplitude_(2, sistrip::invalid_),
58  min_tail_(2, sistrip::invalid_),
59  min_riseTime_(2, sistrip::invalid_),
60  min_decayTime_(2, sistrip::invalid_),
61  min_turnOn_(2, sistrip::invalid_),
62  min_peakTime_(2, sistrip::invalid_),
63  min_undershoot_(2, sistrip::invalid_),
64  min_baseline_(2, sistrip::invalid_),
65  min_smearing_(2, sistrip::invalid_),
66  min_chi2_(2, sistrip::invalid_),
67  max_amplitude_(2, sistrip::invalid_),
68  max_tail_(2, sistrip::invalid_),
69  max_riseTime_(2, sistrip::invalid_),
70  max_decayTime_(2, sistrip::invalid_),
71  max_turnOn_(2, sistrip::invalid_),
72  max_peakTime_(2, sistrip::invalid_),
73  max_undershoot_(2, sistrip::invalid_),
74  max_baseline_(2, sistrip::invalid_),
75  max_smearing_(2, sistrip::invalid_),
76  max_chi2_(2, sistrip::invalid_),
77  spread_amplitude_(2, sistrip::invalid_),
78  spread_tail_(2, sistrip::invalid_),
79  spread_riseTime_(2, sistrip::invalid_),
80  spread_decayTime_(2, sistrip::invalid_),
81  spread_turnOn_(2, sistrip::invalid_),
82  spread_peakTime_(2, sistrip::invalid_),
83  spread_undershoot_(2, sistrip::invalid_),
84  spread_baseline_(2, sistrip::invalid_),
85  spread_smearing_(2, sistrip::invalid_),
86  spread_chi2_(2, sistrip::invalid_),
87  deconv_(deconv),
88  calChan_(0) {
89  ;
90 }
91 
92 // ----------------------------------------------------------------------------
93 //
94 CalibrationAnalysis::CalibrationAnalysis(const bool& deconv)
95  : CommissioningAnalysis("CalibrationAnalysis"),
96  amplitude_(2, VFloat(128, sistrip::invalid_)),
97  tail_(2, VFloat(128, sistrip::invalid_)),
98  riseTime_(2, VFloat(128, sistrip::invalid_)),
99  decayTime_(2, VFloat(128, sistrip::invalid_)),
100  turnOn_(2, VFloat(128, sistrip::invalid_)),
101  peakTime_(2, VFloat(128, sistrip::invalid_)),
102  undershoot_(2, VFloat(128, sistrip::invalid_)),
103  baseline_(2, VFloat(128, sistrip::invalid_)),
104  smearing_(2, VFloat(128, sistrip::invalid_)),
105  chi2_(2, VFloat(128, sistrip::invalid_)),
106  isvalid_(2, VBool(128, sistrip::invalid_)),
107  mean_amplitude_(2, sistrip::invalid_),
108  mean_tail_(2, sistrip::invalid_),
109  mean_riseTime_(2, sistrip::invalid_),
110  mean_decayTime_(2, sistrip::invalid_),
111  mean_turnOn_(2, sistrip::invalid_),
112  mean_peakTime_(2, sistrip::invalid_),
113  mean_undershoot_(2, sistrip::invalid_),
114  mean_baseline_(2, sistrip::invalid_),
115  mean_smearing_(2, sistrip::invalid_),
116  mean_chi2_(2, sistrip::invalid_),
117  min_amplitude_(2, sistrip::invalid_),
118  min_tail_(2, sistrip::invalid_),
119  min_riseTime_(2, sistrip::invalid_),
120  min_decayTime_(2, sistrip::invalid_),
121  min_turnOn_(2, sistrip::invalid_),
122  min_peakTime_(2, sistrip::invalid_),
123  min_undershoot_(2, sistrip::invalid_),
124  min_baseline_(2, sistrip::invalid_),
125  min_smearing_(2, sistrip::invalid_),
126  min_chi2_(2, sistrip::invalid_),
127  max_amplitude_(2, sistrip::invalid_),
128  max_tail_(2, sistrip::invalid_),
129  max_riseTime_(2, sistrip::invalid_),
130  max_decayTime_(2, sistrip::invalid_),
131  max_turnOn_(2, sistrip::invalid_),
132  max_peakTime_(2, sistrip::invalid_),
133  max_undershoot_(2, sistrip::invalid_),
134  max_baseline_(2, sistrip::invalid_),
135  max_smearing_(2, sistrip::invalid_),
136  max_chi2_(2, sistrip::invalid_),
137  spread_amplitude_(2, sistrip::invalid_),
138  spread_tail_(2, sistrip::invalid_),
139  spread_riseTime_(2, sistrip::invalid_),
140  spread_decayTime_(2, sistrip::invalid_),
141  spread_turnOn_(2, sistrip::invalid_),
142  spread_peakTime_(2, sistrip::invalid_),
143  spread_undershoot_(2, sistrip::invalid_),
144  spread_baseline_(2, sistrip::invalid_),
145  spread_smearing_(2, sistrip::invalid_),
146  spread_chi2_(2, sistrip::invalid_),
147  deconv_(deconv),
148  calChan_(0) {
149  ;
150 }
151 
152 // ----------------------------------------------------------------------------
153 //
154 void CalibrationAnalysis::reset() {
155  calChan_ = 0;
156 
157  amplitude_ = VVFloat(2, VFloat(128, sistrip::invalid_));
158  tail_ = VVFloat(2, VFloat(128, sistrip::invalid_));
159  riseTime_ = VVFloat(2, VFloat(128, sistrip::invalid_));
160  decayTime_ = VVFloat(2, VFloat(128, sistrip::invalid_));
161  turnOn_ = VVFloat(2, VFloat(128, sistrip::invalid_));
162  peakTime_ = VVFloat(2, VFloat(128, sistrip::invalid_));
163  undershoot_ = VVFloat(2, VFloat(128, sistrip::invalid_));
164  baseline_ = VVFloat(2, VFloat(128, sistrip::invalid_));
165  smearing_ = VVFloat(2, VFloat(128, sistrip::invalid_));
166  chi2_ = VVFloat(2, VFloat(128, sistrip::invalid_));
167  isvalid_ = VVBool(2, VBool(128, sistrip::invalid_));
168 
169  mean_amplitude_ = VFloat(2, sistrip::invalid_);
170  mean_tail_ = VFloat(2, sistrip::invalid_);
171  mean_riseTime_ = VFloat(2, sistrip::invalid_);
172  mean_decayTime_ = VFloat(2, sistrip::invalid_);
173  mean_turnOn_ = VFloat(2, sistrip::invalid_);
174  mean_peakTime_ = VFloat(2, sistrip::invalid_);
175  mean_undershoot_ = VFloat(2, sistrip::invalid_);
176  mean_baseline_ = VFloat(2, sistrip::invalid_);
177  mean_smearing_ = VFloat(2, sistrip::invalid_);
178  mean_chi2_ = VFloat(2, sistrip::invalid_);
179 
180  min_amplitude_ = VFloat(2, sistrip::invalid_);
181  min_tail_ = VFloat(2, sistrip::invalid_);
182  min_riseTime_ = VFloat(2, sistrip::invalid_);
183  min_decayTime_ = VFloat(2, sistrip::invalid_);
184  min_turnOn_ = VFloat(2, sistrip::invalid_);
185  min_peakTime_ = VFloat(2, sistrip::invalid_);
186  min_undershoot_ = VFloat(2, sistrip::invalid_);
187  min_baseline_ = VFloat(2, sistrip::invalid_);
188  min_smearing_ = VFloat(2, sistrip::invalid_);
189  min_chi2_ = VFloat(2, sistrip::invalid_);
190 
191  max_amplitude_ = VFloat(2, sistrip::invalid_);
192  max_tail_ = VFloat(2, sistrip::invalid_);
193  max_riseTime_ = VFloat(2, sistrip::invalid_);
194  max_decayTime_ = VFloat(2, sistrip::invalid_);
195  max_turnOn_ = VFloat(2, sistrip::invalid_);
196  max_peakTime_ = VFloat(2, sistrip::invalid_);
197  max_undershoot_ = VFloat(2, sistrip::invalid_);
198  max_baseline_ = VFloat(2, sistrip::invalid_);
199  max_smearing_ = VFloat(2, sistrip::invalid_);
200  max_chi2_ = VFloat(2, sistrip::invalid_);
201 
202  spread_amplitude_ = VFloat(2, sistrip::invalid_);
203  spread_tail_ = VFloat(2, sistrip::invalid_);
204  spread_riseTime_ = VFloat(2, sistrip::invalid_);
205  spread_decayTime_ = VFloat(2, sistrip::invalid_);
206  spread_turnOn_ = VFloat(2, sistrip::invalid_);
207  spread_peakTime_ = VFloat(2, sistrip::invalid_);
208  spread_undershoot_ = VFloat(2, sistrip::invalid_);
209  spread_baseline_ = VFloat(2, sistrip::invalid_);
210  spread_smearing_ = VFloat(2, sistrip::invalid_);
211  spread_chi2_ = VFloat(2, sistrip::invalid_);
212 }
213 
214 // ----------------------------------------------------------------------------
215 //
216 void CalibrationAnalysis::print(std::stringstream& ss, uint32_t iapv) {
217  header(ss);
218  ss << " Monitorables for APV number : " << iapv;
219  if (iapv == 0) {
220  ss << " (first of pair)";
221  } else if (iapv == 1) {
222  ss << " (second of pair)";
223  }
224  ss << std::endl;
225  ss << " Mean Amplitude of the pulse : " << mean_amplitude_[iapv] << std::endl
226  << " Mean Tail amplitude after 150ns : " << mean_tail_[iapv] << std::endl
227  << " Mean Rise time : " << mean_riseTime_[iapv] << std::endl
228  << " Mean Time constant : " << mean_decayTime_[iapv] << std::endl
229  << " Mean Turn on time : " << mean_turnOn_[iapv] << std::endl
230  << " Mean peak time : " << mean_peakTime_[iapv] << std::endl
231  << " Mean undershoot amplitude : " << mean_undershoot_[iapv] << std::endl
232  << " Mean baseline amplitude : " << mean_baseline_[iapv] << std::endl
233  << " Mean Smearing parameter : " << mean_smearing_[iapv] << std::endl
234  << " Mean Chi2 of the fit : " << mean_chi2_[iapv] << std::endl;
235  if (deconvMode()) {
236  ss << "Data obtained in deconvolution mode." << std::endl;
237  } else {
238  ss << "Data obtained in peak mode." << std::endl;
239  }
240 }
241 
242 bool CalibrationAnalysis::isValid() const { return true; }
CalibrationAnalysis::minRiseTimeThresholdDeco_
static const float minRiseTimeThresholdDeco_
Definition: CalibrationAnalysis.h:119
CalibrationAnalysis::maxTurnOnThreshold_
static const float maxTurnOnThreshold_
Definition: CalibrationAnalysis.h:113
CalibrationAnalysis::maxRiseTimeThresholdDeco_
static const float maxRiseTimeThresholdDeco_
Definition: CalibrationAnalysis.h:120
CalibrationAnalysis::minPeakTimeThreshold_
static const float minPeakTimeThreshold_
Definition: CalibrationAnalysis.h:108
CalibrationAnalysis::CalibrationAnalysis
CalibrationAnalysis(const uint32_t &key, const bool &deconv)
Definition: CalibrationAnalysis.cc:34
CalibrationAnalysis::maxDecayTimeThresholdDeco_
static const float maxDecayTimeThresholdDeco_
Definition: CalibrationAnalysis.h:116
CalibrationAnalysis::maxRiseTimeThreshold_
static const float maxRiseTimeThreshold_
Definition: CalibrationAnalysis.h:111
CalibrationAnalysis::maxBaselineThreshold_
static const float maxBaselineThreshold_
Definition: CalibrationAnalysis.h:103
CalibrationAnalysis::print
void print(std::stringstream &, uint32_t not_used=0) override
Definition: CalibrationAnalysis.cc:216
sistrip
sistrip classes
Definition: SiStripQualityHelpers.h:14
CalibrationAnalysis::minDecayTimeThresholdDeco_
static const float minDecayTimeThresholdDeco_
Definition: CalibrationAnalysis.h:115
CalibrationAnalysis::minAmplitudeThreshold_
static const float minAmplitudeThreshold_
Definition: CalibrationAnalysis.h:101
CalibrationAnalysis::maxDecayTimeThreshold_
static const float maxDecayTimeThreshold_
Definition: CalibrationAnalysis.h:107
CalibrationAnalysis::minTurnOnThresholdDeco_
static const float minTurnOnThresholdDeco_
Definition: CalibrationAnalysis.h:121
CalibrationAnalysis::minDecayTimeThreshold_
static const float minDecayTimeThreshold_
Definition: CalibrationAnalysis.h:106
CommissioningAnalysis::VFloat
std::vector< float > VFloat
Definition: CommissioningAnalysis.h:30
CommissioningAnalysis::VBool
std::vector< bool > VBool
Definition: CommissioningAnalysis.h:38
CalibrationAnalysis::minPeakTimeThresholdDeco_
static const float minPeakTimeThresholdDeco_
Definition: CalibrationAnalysis.h:117
CalibrationAnalysis::minBaselineThreshold_
static const float minBaselineThreshold_
Definition: CalibrationAnalysis.h:102
CalibrationAnalysis::maxTurnOnThresholdDeco_
static const float maxTurnOnThresholdDeco_
Definition: CalibrationAnalysis.h:122
CalibrationAnalysis::minRiseTimeThreshold_
static const float minRiseTimeThreshold_
Definition: CalibrationAnalysis.h:110
CommissioningAnalysis::header
virtual void header(std::stringstream &) const
Definition: CommissioningAnalysis.cc:33
CalibrationAnalysis::isValid
bool isValid() const override
Definition: CalibrationAnalysis.cc:242
sistrip::invalid_
static const uint16_t invalid_
Definition: Constants.h:16
CommissioningAnalysis::VVBool
std::vector< VBool > VVBool
Definition: CommissioningAnalysis.h:40
CalibrationAnalysis::minTurnOnThreshold_
static const float minTurnOnThreshold_
Definition: CalibrationAnalysis.h:112
CalibrationAnalysis::maxChi2Threshold_
static const float maxChi2Threshold_
Definition: CalibrationAnalysis.h:104
CalibrationAnalysis::maxPeakTimeThresholdDeco_
static const float maxPeakTimeThresholdDeco_
Definition: CalibrationAnalysis.h:118
CalibrationAnalysis::maxPeakTimeThreshold_
static const float maxPeakTimeThreshold_
Definition: CalibrationAnalysis.h:109
CommissioningAnalysis::VVFloat
std::vector< VFloat > VVFloat
Definition: CommissioningAnalysis.h:32
CommissioningAnalysis
Abstract base for derived classes that provide analysis of commissioning histograms.
Definition: CommissioningAnalysis.h:18
Generated for CMSSW Reference Manual by doxygen 1.8.11