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 CalibrationAnalysis::CalibrationAnalysis(const bool& deconv)
94  : CommissioningAnalysis("CalibrationAnalysis"),
95  amplitude_(2,VFloat(128,sistrip::invalid_)),
96  tail_(2,VFloat(128,sistrip::invalid_)),
97  riseTime_(2,VFloat(128,sistrip::invalid_)),
98  decayTime_(2,VFloat(128,sistrip::invalid_)),
99  turnOn_(2,VFloat(128,sistrip::invalid_)),
100  peakTime_(2,VFloat(128,sistrip::invalid_)),
101  undershoot_(2,VFloat(128,sistrip::invalid_)),
102  baseline_(2,VFloat(128,sistrip::invalid_)),
103  smearing_(2,VFloat(128,sistrip::invalid_)),
104  chi2_(2,VFloat(128,sistrip::invalid_)),
105  isvalid_(2,VBool(128,sistrip::invalid_)),
106  mean_amplitude_(2,sistrip::invalid_),
107  mean_tail_(2,sistrip::invalid_),
108  mean_riseTime_(2,sistrip::invalid_),
109  mean_decayTime_(2,sistrip::invalid_),
110  mean_turnOn_(2,sistrip::invalid_),
111  mean_peakTime_(2,sistrip::invalid_),
112  mean_undershoot_(2,sistrip::invalid_),
113  mean_baseline_(2,sistrip::invalid_),
114  mean_smearing_(2,sistrip::invalid_),
115  mean_chi2_(2,sistrip::invalid_),
116  min_amplitude_(2,sistrip::invalid_),
117  min_tail_(2,sistrip::invalid_),
118  min_riseTime_(2,sistrip::invalid_),
119  min_decayTime_(2,sistrip::invalid_),
120  min_turnOn_(2,sistrip::invalid_),
121  min_peakTime_(2,sistrip::invalid_),
122  min_undershoot_(2,sistrip::invalid_),
123  min_baseline_(2,sistrip::invalid_),
124  min_smearing_(2,sistrip::invalid_),
125  min_chi2_(2,sistrip::invalid_),
126  max_amplitude_(2,sistrip::invalid_),
127  max_tail_(2,sistrip::invalid_),
128  max_riseTime_(2,sistrip::invalid_),
129  max_decayTime_(2,sistrip::invalid_),
130  max_turnOn_(2,sistrip::invalid_),
131  max_peakTime_(2,sistrip::invalid_),
132  max_undershoot_(2,sistrip::invalid_),
133  max_baseline_(2,sistrip::invalid_),
134  max_smearing_(2,sistrip::invalid_),
135  max_chi2_(2,sistrip::invalid_),
136  spread_amplitude_(2,sistrip::invalid_),
137  spread_tail_(2,sistrip::invalid_),
138  spread_riseTime_(2,sistrip::invalid_),
139  spread_decayTime_(2,sistrip::invalid_),
140  spread_turnOn_(2,sistrip::invalid_),
141  spread_peakTime_(2,sistrip::invalid_),
142  spread_undershoot_(2,sistrip::invalid_),
143  spread_baseline_(2,sistrip::invalid_),
144  spread_smearing_(2,sistrip::invalid_),
145  spread_chi2_(2,sistrip::invalid_),
146  deconv_(deconv),
147  calChan_(0)
148 {;}
149 
150 // ----------------------------------------------------------------------------
151 //
152 void CalibrationAnalysis::reset() {
153 
154  calChan_ = 0;
155 
156  amplitude_ = VVFloat(2,VFloat(128,sistrip::invalid_));
157  tail_ = VVFloat(2,VFloat(128,sistrip::invalid_));
158  riseTime_ = VVFloat(2,VFloat(128,sistrip::invalid_));
159  decayTime_ = VVFloat(2,VFloat(128,sistrip::invalid_));
160  turnOn_ = VVFloat(2,VFloat(128,sistrip::invalid_));
161  peakTime_ = VVFloat(2,VFloat(128,sistrip::invalid_));
162  undershoot_ = VVFloat(2,VFloat(128,sistrip::invalid_));
163  baseline_ = VVFloat(2,VFloat(128,sistrip::invalid_));
164  smearing_ = VVFloat(2,VFloat(128,sistrip::invalid_));
165  chi2_ = VVFloat(2,VFloat(128,sistrip::invalid_));
166  isvalid_ = VVBool(2,VBool(128,sistrip::invalid_));
167 
168  mean_amplitude_ = VFloat(2,sistrip::invalid_);
169  mean_tail_ = VFloat(2,sistrip::invalid_);
170  mean_riseTime_ = VFloat(2,sistrip::invalid_);
171  mean_decayTime_ = VFloat(2,sistrip::invalid_);
172  mean_turnOn_ = VFloat(2,sistrip::invalid_);
173  mean_peakTime_ = VFloat(2,sistrip::invalid_);
174  mean_undershoot_ = VFloat(2,sistrip::invalid_);
175  mean_baseline_ = VFloat(2,sistrip::invalid_);
176  mean_smearing_ = VFloat(2,sistrip::invalid_);
177  mean_chi2_ = VFloat(2,sistrip::invalid_);
178 
179  min_amplitude_ = VFloat(2,sistrip::invalid_);
180  min_tail_ = VFloat(2,sistrip::invalid_);
181  min_riseTime_ = VFloat(2,sistrip::invalid_);
182  min_decayTime_ = VFloat(2,sistrip::invalid_);
183  min_turnOn_ = VFloat(2,sistrip::invalid_);
184  min_peakTime_ = VFloat(2,sistrip::invalid_);
185  min_undershoot_ = VFloat(2,sistrip::invalid_);
186  min_baseline_ = VFloat(2,sistrip::invalid_);
187  min_smearing_ = VFloat(2,sistrip::invalid_);
188  min_chi2_ = VFloat(2,sistrip::invalid_);
189 
190  max_amplitude_ = VFloat(2,sistrip::invalid_);
191  max_tail_ = VFloat(2,sistrip::invalid_);
192  max_riseTime_ = VFloat(2,sistrip::invalid_);
193  max_decayTime_ = VFloat(2,sistrip::invalid_);
194  max_turnOn_ = VFloat(2,sistrip::invalid_);
195  max_peakTime_ = VFloat(2,sistrip::invalid_);
196  max_undershoot_ = VFloat(2,sistrip::invalid_);
197  max_baseline_ = VFloat(2,sistrip::invalid_);
198  max_smearing_ = VFloat(2,sistrip::invalid_);
199  max_chi2_ = VFloat(2,sistrip::invalid_);
200 
201  spread_amplitude_ = VFloat(2,sistrip::invalid_);
202  spread_tail_ = VFloat(2,sistrip::invalid_);
203  spread_riseTime_ = VFloat(2,sistrip::invalid_);
204  spread_decayTime_ = VFloat(2,sistrip::invalid_);
205  spread_turnOn_ = VFloat(2,sistrip::invalid_);
206  spread_peakTime_ = VFloat(2,sistrip::invalid_);
207  spread_undershoot_ = VFloat(2,sistrip::invalid_);
208  spread_baseline_ = VFloat(2,sistrip::invalid_);
209  spread_smearing_ = VFloat(2,sistrip::invalid_);
210  spread_chi2_ = VFloat(2,sistrip::invalid_);
211 }
212 
213 // ----------------------------------------------------------------------------
214 //
215 void CalibrationAnalysis::print( std::stringstream& ss, uint32_t iapv ) {
216  header( ss );
217  ss << " Monitorables for APV number : " << iapv;
218  if ( iapv == 0 ) { ss << " (first of pair)"; }
219  else if ( iapv == 1 ) { ss << " (second of pair)"; }
220  ss << std::endl;
221  ss << " Mean Amplitude of the pulse : " << mean_amplitude_[iapv] << std::endl
222  << " Mean Tail amplitude after 150ns : " << mean_tail_[iapv] << std::endl
223  << " Mean Rise time : " << mean_riseTime_[iapv] << std::endl
224  << " Mean Time constant : " << mean_decayTime_[iapv] << std::endl
225  << " Mean Turn on time : " << mean_turnOn_[iapv] << std::endl
226  << " Mean peak time : " << mean_peakTime_[iapv] << std::endl
227  << " Mean undershoot amplitude : " << mean_undershoot_[iapv] << std::endl
228  << " Mean baseline amplitude : " << mean_baseline_[iapv] << std::endl
229  << " Mean Smearing parameter : " << mean_smearing_[iapv] << std::endl
230  << " Mean Chi2 of the fit : " << mean_chi2_[iapv] << std::endl;
231  if(deconvMode()) {
232  ss << "Data obtained in deconvolution mode." << std::endl;
233  } else {
234  ss << "Data obtained in peak mode." << std::endl;
235  }
236 }
237 
238 bool CalibrationAnalysis::isValid() const { return true;}
CalibrationAnalysis::minRiseTimeThresholdDeco_
static const float minRiseTimeThresholdDeco_
Definition: CalibrationAnalysis.h:122
CalibrationAnalysis::maxTurnOnThreshold_
static const float maxTurnOnThreshold_
Definition: CalibrationAnalysis.h:116
CalibrationAnalysis::maxRiseTimeThresholdDeco_
static const float maxRiseTimeThresholdDeco_
Definition: CalibrationAnalysis.h:123
CalibrationAnalysis::minPeakTimeThreshold_
static const float minPeakTimeThreshold_
Definition: CalibrationAnalysis.h:111
CalibrationAnalysis::CalibrationAnalysis
CalibrationAnalysis(const uint32_t &key, const bool &deconv)
Definition: CalibrationAnalysis.cc:34
CalibrationAnalysis::maxDecayTimeThresholdDeco_
static const float maxDecayTimeThresholdDeco_
Definition: CalibrationAnalysis.h:119
CalibrationAnalysis::maxRiseTimeThreshold_
static const float maxRiseTimeThreshold_
Definition: CalibrationAnalysis.h:114
CalibrationAnalysis::maxBaselineThreshold_
static const float maxBaselineThreshold_
Definition: CalibrationAnalysis.h:106
CalibrationAnalysis::print
void print(std::stringstream &, uint32_t not_used=0) override
Definition: CalibrationAnalysis.cc:215
sistrip
sistrip classes
Definition: EnsembleCalibrationLA.cc:11
CalibrationAnalysis::minDecayTimeThresholdDeco_
static const float minDecayTimeThresholdDeco_
Definition: CalibrationAnalysis.h:118
CalibrationAnalysis::minAmplitudeThreshold_
static const float minAmplitudeThreshold_
Definition: CalibrationAnalysis.h:104
CalibrationAnalysis::maxDecayTimeThreshold_
static const float maxDecayTimeThreshold_
Definition: CalibrationAnalysis.h:110
CalibrationAnalysis::minTurnOnThresholdDeco_
static const float minTurnOnThresholdDeco_
Definition: CalibrationAnalysis.h:124
CalibrationAnalysis::minDecayTimeThreshold_
static const float minDecayTimeThreshold_
Definition: CalibrationAnalysis.h:109
CommissioningAnalysis::VFloat
std::vector< float > VFloat
Definition: CommissioningAnalysis.h:33
CommissioningAnalysis::VBool
std::vector< bool > VBool
Definition: CommissioningAnalysis.h:41
CalibrationAnalysis::minPeakTimeThresholdDeco_
static const float minPeakTimeThresholdDeco_
Definition: CalibrationAnalysis.h:120
CalibrationAnalysis::minBaselineThreshold_
static const float minBaselineThreshold_
Definition: CalibrationAnalysis.h:105
CalibrationAnalysis::maxTurnOnThresholdDeco_
static const float maxTurnOnThresholdDeco_
Definition: CalibrationAnalysis.h:125
CalibrationAnalysis::minRiseTimeThreshold_
static const float minRiseTimeThreshold_
Definition: CalibrationAnalysis.h:113
CommissioningAnalysis::header
virtual void header(std::stringstream &) const
Definition: CommissioningAnalysis.cc:32
CalibrationAnalysis::isValid
bool isValid() const override
Definition: CalibrationAnalysis.cc:238
sistrip::invalid_
static const uint16_t invalid_
Definition: Constants.h:16
CommissioningAnalysis::VVBool
std::vector< VBool > VVBool
Definition: CommissioningAnalysis.h:43
CalibrationAnalysis::minTurnOnThreshold_
static const float minTurnOnThreshold_
Definition: CalibrationAnalysis.h:115
CalibrationAnalysis::maxChi2Threshold_
static const float maxChi2Threshold_
Definition: CalibrationAnalysis.h:107
CalibrationAnalysis::maxPeakTimeThresholdDeco_
static const float maxPeakTimeThresholdDeco_
Definition: CalibrationAnalysis.h:121
CalibrationAnalysis::maxPeakTimeThreshold_
static const float maxPeakTimeThreshold_
Definition: CalibrationAnalysis.h:112
CommissioningAnalysis::VVFloat
std::vector< VFloat > VVFloat
Definition: CommissioningAnalysis.h:35
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