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 // ----------------------------------------------------------------------------
13 //
14 CalibrationAnalysis::CalibrationAnalysis( const uint32_t& key, const bool& deconv, int calchan )
15  : CommissioningAnalysis(key,"CalibrationAnalysis"),
16  amplitude_(2,VFloat(128,sistrip::invalid_)),
17  tail_(2,VFloat(128,sistrip::invalid_)),
18  riseTime_(2,VFloat(128,sistrip::invalid_)),
19  timeConstant_(2,VFloat(128,sistrip::invalid_)),
20  turnOn_(2,VFloat(128,sistrip::invalid_)),
21  maximum_(2,VFloat(128,sistrip::invalid_)),
22  undershoot_(2,VFloat(128,sistrip::invalid_)),
23  baseline_(2,VFloat(128,sistrip::invalid_)),
24  smearing_(2,VFloat(128,sistrip::invalid_)),
25  chi2_(2,VFloat(128,sistrip::invalid_)),
26  mean_amplitude_(2,sistrip::invalid_),
27  mean_tail_(2,sistrip::invalid_),
28  mean_riseTime_(2,sistrip::invalid_),
29  mean_timeConstant_(2,sistrip::invalid_),
30  mean_turnOn_(2,sistrip::invalid_),
31  mean_maximum_(2,sistrip::invalid_),
32  mean_undershoot_(2,sistrip::invalid_),
33  mean_baseline_(2,sistrip::invalid_),
34  mean_smearing_(2,sistrip::invalid_),
35  mean_chi2_(2,sistrip::invalid_),
36  min_amplitude_(2,sistrip::invalid_),
37  min_tail_(2,sistrip::invalid_),
38  min_riseTime_(2,sistrip::invalid_),
39  min_timeConstant_(2,sistrip::invalid_),
40  min_turnOn_(2,sistrip::invalid_),
41  min_maximum_(2,sistrip::invalid_),
42  min_undershoot_(2,sistrip::invalid_),
43  min_baseline_(2,sistrip::invalid_),
44  min_smearing_(2,sistrip::invalid_),
45  min_chi2_(2,sistrip::invalid_),
46  max_amplitude_(2,sistrip::invalid_),
47  max_tail_(2,sistrip::invalid_),
48  max_riseTime_(2,sistrip::invalid_),
49  max_timeConstant_(2,sistrip::invalid_),
50  max_turnOn_(2,sistrip::invalid_),
51  max_maximum_(2,sistrip::invalid_),
52  max_undershoot_(2,sistrip::invalid_),
53  max_baseline_(2,sistrip::invalid_),
54  max_smearing_(2,sistrip::invalid_),
55  max_chi2_(2,sistrip::invalid_),
56  spread_amplitude_(2,sistrip::invalid_),
57  spread_tail_(2,sistrip::invalid_),
58  spread_riseTime_(2,sistrip::invalid_),
59  spread_timeConstant_(2,sistrip::invalid_),
60  spread_turnOn_(2,sistrip::invalid_),
61  spread_maximum_(2,sistrip::invalid_),
62  spread_undershoot_(2,sistrip::invalid_),
63  spread_baseline_(2,sistrip::invalid_),
64  spread_smearing_(2,sistrip::invalid_),
65  spread_chi2_(2,sistrip::invalid_),
66  deconv_(deconv),
67  calchan_(calchan),
68  isScan_(false)
69 {;}
70 
71 // ----------------------------------------------------------------------------
72 //
73 CalibrationAnalysis::CalibrationAnalysis(const bool& deconv, int calchan)
74  : CommissioningAnalysis("CalibrationAnalysis"),
75  amplitude_(2,VFloat(128,sistrip::invalid_)),
76  tail_(2,VFloat(128,sistrip::invalid_)),
77  riseTime_(2,VFloat(128,sistrip::invalid_)),
78  timeConstant_(2,VFloat(128,sistrip::invalid_)),
79  turnOn_(2,VFloat(128,sistrip::invalid_)),
80  maximum_(2,VFloat(128,sistrip::invalid_)),
81  undershoot_(2,VFloat(128,sistrip::invalid_)),
82  baseline_(2,VFloat(128,sistrip::invalid_)),
83  smearing_(2,VFloat(128,sistrip::invalid_)),
84  chi2_(2,VFloat(128,sistrip::invalid_)),
85  mean_amplitude_(2,sistrip::invalid_),
86  mean_tail_(2,sistrip::invalid_),
87  mean_riseTime_(2,sistrip::invalid_),
88  mean_timeConstant_(2,sistrip::invalid_),
89  mean_turnOn_(2,sistrip::invalid_),
90  mean_maximum_(2,sistrip::invalid_),
91  mean_undershoot_(2,sistrip::invalid_),
92  mean_baseline_(2,sistrip::invalid_),
93  mean_smearing_(2,sistrip::invalid_),
94  mean_chi2_(2,sistrip::invalid_),
95  min_amplitude_(2,sistrip::invalid_),
96  min_tail_(2,sistrip::invalid_),
97  min_riseTime_(2,sistrip::invalid_),
98  min_timeConstant_(2,sistrip::invalid_),
99  min_turnOn_(2,sistrip::invalid_),
100  min_maximum_(2,sistrip::invalid_),
101  min_undershoot_(2,sistrip::invalid_),
102  min_baseline_(2,sistrip::invalid_),
103  min_smearing_(2,sistrip::invalid_),
104  min_chi2_(2,sistrip::invalid_),
105  max_amplitude_(2,sistrip::invalid_),
106  max_tail_(2,sistrip::invalid_),
107  max_riseTime_(2,sistrip::invalid_),
108  max_timeConstant_(2,sistrip::invalid_),
109  max_turnOn_(2,sistrip::invalid_),
110  max_maximum_(2,sistrip::invalid_),
111  max_undershoot_(2,sistrip::invalid_),
112  max_baseline_(2,sistrip::invalid_),
113  max_smearing_(2,sistrip::invalid_),
114  max_chi2_(2,sistrip::invalid_),
115  spread_amplitude_(2,sistrip::invalid_),
116  spread_tail_(2,sistrip::invalid_),
117  spread_riseTime_(2,sistrip::invalid_),
118  spread_timeConstant_(2,sistrip::invalid_),
119  spread_turnOn_(2,sistrip::invalid_),
120  spread_maximum_(2,sistrip::invalid_),
121  spread_undershoot_(2,sistrip::invalid_),
122  spread_baseline_(2,sistrip::invalid_),
123  spread_smearing_(2,sistrip::invalid_),
124  spread_chi2_(2,sistrip::invalid_),
125  deconv_(deconv),
126  calchan_(calchan),
127  isScan_(false)
128 {;}
129 
130 // ----------------------------------------------------------------------------
131 //
132 void CalibrationAnalysis::reset() {
133  amplitude_ = VVFloat(2,VFloat(128,sistrip::invalid_));
134  tail_ = VVFloat(2,VFloat(128,sistrip::invalid_));
135  riseTime_ = VVFloat(2,VFloat(128,sistrip::invalid_));
136  timeConstant_ = VVFloat(2,VFloat(128,sistrip::invalid_));
137  turnOn_ = VVFloat(2,VFloat(128,sistrip::invalid_));
138  maximum_ = VVFloat(2,VFloat(128,sistrip::invalid_));
139  undershoot_ = VVFloat(2,VFloat(128,sistrip::invalid_));
140  baseline_ = VVFloat(2,VFloat(128,sistrip::invalid_));
141  smearing_ = VVFloat(2,VFloat(128,sistrip::invalid_));
142  chi2_ = VVFloat(2,VFloat(128,sistrip::invalid_));
143  mean_amplitude_ = VFloat(2,sistrip::invalid_);
144  mean_tail_ = VFloat(2,sistrip::invalid_);
145  mean_riseTime_ = VFloat(2,sistrip::invalid_);
146  mean_timeConstant_ = VFloat(2,sistrip::invalid_);
147  mean_turnOn_ = VFloat(2,sistrip::invalid_);
148  mean_maximum_ = VFloat(2,sistrip::invalid_);
149  mean_undershoot_ = VFloat(2,sistrip::invalid_);
150  mean_baseline_ = VFloat(2,sistrip::invalid_);
151  mean_smearing_ = VFloat(2,sistrip::invalid_);
152  mean_chi2_ = VFloat(2,sistrip::invalid_);
153  min_amplitude_ = VFloat(2,sistrip::invalid_);
154  min_tail_ = VFloat(2,sistrip::invalid_);
155  min_riseTime_ = VFloat(2,sistrip::invalid_);
156  min_timeConstant_ = VFloat(2,sistrip::invalid_);
157  min_turnOn_ = VFloat(2,sistrip::invalid_);
158  min_maximum_ = VFloat(2,sistrip::invalid_);
159  min_undershoot_ = VFloat(2,sistrip::invalid_);
160  min_baseline_ = VFloat(2,sistrip::invalid_);
161  min_smearing_ = VFloat(2,sistrip::invalid_);
162  min_chi2_ = VFloat(2,sistrip::invalid_);
163  max_amplitude_ = VFloat(2,sistrip::invalid_);
164  max_tail_ = VFloat(2,sistrip::invalid_);
165  max_riseTime_ = VFloat(2,sistrip::invalid_);
166  max_timeConstant_ = VFloat(2,sistrip::invalid_);
167  max_turnOn_ = VFloat(2,sistrip::invalid_);
168  max_maximum_ = VFloat(2,sistrip::invalid_);
169  max_undershoot_ = VFloat(2,sistrip::invalid_);
170  max_baseline_ = VFloat(2,sistrip::invalid_);
171  max_smearing_ = VFloat(2,sistrip::invalid_);
172  max_chi2_ = VFloat(2,sistrip::invalid_);
173  spread_amplitude_ = VFloat(2,sistrip::invalid_);
174  spread_tail_ = VFloat(2,sistrip::invalid_);
175  spread_riseTime_ = VFloat(2,sistrip::invalid_);
176  spread_timeConstant_ = VFloat(2,sistrip::invalid_);
177  spread_turnOn_ = VFloat(2,sistrip::invalid_);
178  spread_maximum_ = VFloat(2,sistrip::invalid_);
179  spread_undershoot_ = VFloat(2,sistrip::invalid_);
180  spread_baseline_ = VFloat(2,sistrip::invalid_);
181  spread_smearing_ = VFloat(2,sistrip::invalid_);
182  spread_chi2_ = VFloat(2,sistrip::invalid_);
183 }
184 
185 // ----------------------------------------------------------------------------
186 //
187 void CalibrationAnalysis::print( std::stringstream& ss, uint32_t iapv ) {
188  header( ss );
189  ss << " Monitorables for APV number : " << iapv;
190  if ( iapv == 0 ) { ss << " (first of pair)"; }
191  else if ( iapv == 1 ) { ss << " (second of pair)"; }
192  ss << std::endl;
193  ss << " Mean Amplitude of the pulse : " << mean_amplitude_[iapv] << std::endl
194  << " Mean Tail amplitude after 150ns : " << mean_tail_[iapv] << std::endl
195  << " Mean Rise time : " << mean_riseTime_[iapv] << std::endl
196  << " Mean Time constant : " << mean_timeConstant_[iapv] << std::endl
197  << " Mean Turn on time : " << mean_turnOn_[iapv] << std::endl
198  << " Mean peak time : " << mean_maximum_[iapv] << std::endl
199  << " Mean undershoot amplitude : " << mean_undershoot_[iapv] << std::endl
200  << " Mean baseline amplitude : " << mean_baseline_[iapv] << std::endl
201  << " Mean Smearing parameter : " << mean_smearing_[iapv] << std::endl
202  << " Mean Chi2 of the fit : " << mean_chi2_[iapv] << std::endl;
203  if(deconvMode()) {
204  ss << "Data obtained in deconvolution mode." << std::endl;
205  } else {
206  ss << "Data obtained in peak mode." << std::endl;
207  }
208 }
209 
CalibrationAnalysis::print
void print(std::stringstream &, uint32_t not_used=0) override
Definition: CalibrationAnalysis.cc:187
sistrip
sistrip classes
Definition: EnsembleCalibrationLA.cc:12
CommissioningAnalysis::VFloat
std::vector< float > VFloat
Definition: CommissioningAnalysis.h:33
sistrip::maximum_
static const uint16_t maximum_
Definition: Constants.h:20
CommissioningAnalysis::header
virtual void header(std::stringstream &) const
Definition: CommissioningAnalysis.cc:32
CalibrationAnalysis::CalibrationAnalysis
CalibrationAnalysis(const uint32_t &key, const bool &deconv, int calchan)
Definition: CalibrationAnalysis.cc:14
sistrip::invalid_
static const uint16_t invalid_
Definition: Constants.h:16
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