CMS 3D CMS Logo
DaqScopeModeAnalysis.cc
Go to the documentation of this file.
1 #include "CondFormats/SiStripObjects/interface/DaqScopeModeAnalysis.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 <iomanip>
7 #include <cmath>
8 
9 using namespace sistrip;
10 
11 const float DaqScopeModeAnalysis::tickMarkHeightThreshold_ = 50.; // [ADC]
12 const float DaqScopeModeAnalysis::frameFindingThreshold_ = (2./3.); // fraction of tick mark height
13 
14 // ----------------------------------------------------------------------------
15 //
16 DaqScopeModeAnalysis::DaqScopeModeAnalysis( const uint32_t& key )
17  : CommissioningAnalysis(key,"DaqScopeModeAnalysis"),
18  height_(1.*sistrip::invalid_),
19  base_(1.*sistrip::invalid_),
20  peak_(1.*sistrip::invalid_),
21  peds_(2,VFloat(128,sistrip::invalid_)),
22  noise_(2,VFloat(128,sistrip::invalid_)),
23  raw_(2,VFloat(128,sistrip::invalid_)),
24  dead_(2,VInt(0,sistrip::invalid_)),
25  noisy_(2,VInt(0,sistrip::invalid_)),
26  pedsMean_(2,sistrip::invalid_),
27  pedsSpread_(2,sistrip::invalid_),
28  noiseMean_(2,sistrip::invalid_),
29  noiseSpread_(2,sistrip::invalid_),
30  rawMean_(2,sistrip::invalid_),
31  rawSpread_(2,sistrip::invalid_),
32  pedsMax_(2,sistrip::invalid_),
33  pedsMin_(2,sistrip::invalid_),
34  noiseMax_(2,sistrip::invalid_),
35  noiseMin_(2,sistrip::invalid_),
36  rawMax_(2,sistrip::invalid_),
37  rawMin_(2,sistrip::invalid_)
38 {
39  dead_[0].reserve(256); dead_[1].reserve(256);
40  noisy_[0].reserve(256); noisy_[1].reserve(256);
41 }
42 
43 // ----------------------------------------------------------------------------
44 //
45 DaqScopeModeAnalysis::DaqScopeModeAnalysis()
46  : CommissioningAnalysis("DaqScopeModeAnalysis"),
47  height_(1.*sistrip::invalid_),
48  base_(1.*sistrip::invalid_),
49  peak_(1.*sistrip::invalid_),
50  peds_(2,VFloat(128,sistrip::invalid_)),
51  noise_(2,VFloat(128,sistrip::invalid_)),
52  raw_(2,VFloat(128,sistrip::invalid_)),
53  dead_(2,VInt(0,sistrip::invalid_)),
54  noisy_(2,VInt(0,sistrip::invalid_)),
55  pedsMean_(2,sistrip::invalid_),
56  pedsSpread_(2,sistrip::invalid_),
57  noiseMean_(2,sistrip::invalid_),
58  noiseSpread_(2,sistrip::invalid_),
59  rawMean_(2,sistrip::invalid_),
60  rawSpread_(2,sistrip::invalid_),
61  pedsMax_(2,sistrip::invalid_),
62  pedsMin_(2,sistrip::invalid_),
63  noiseMax_(2,sistrip::invalid_),
64  noiseMin_(2,sistrip::invalid_),
65  rawMax_(2,sistrip::invalid_),
66  rawMin_(2,sistrip::invalid_)
67 {
68  dead_[0].reserve(256); dead_[1].reserve(256);
69  noisy_[0].reserve(256); noisy_[1].reserve(256);
70 }
71 
72 // ----------------------------------------------------------------------------
73 //
74 void DaqScopeModeAnalysis::reset() {
75  height_ = 1.*sistrip::invalid_;
76  base_ = 1.*sistrip::invalid_;
77  peak_ = 1.*sistrip::invalid_;
78 
79  peds_ = VVFloat(2,VFloat(128,sistrip::invalid_));
80  noise_ = VVFloat(2,VFloat(128,sistrip::invalid_));
81  raw_ = VVFloat(2,VFloat(128,sistrip::invalid_));
82  dead_ = VVInt(2,VInt(0,sistrip::invalid_));
83  noisy_ = VVInt(2,VInt(0,sistrip::invalid_));
84  pedsMean_ = VFloat(2,sistrip::invalid_);
85  pedsSpread_ = VFloat(2,sistrip::invalid_);
86  noiseMean_ = VFloat(2,sistrip::invalid_);
87  noiseSpread_ = VFloat(2,sistrip::invalid_);
88  rawMean_ = VFloat(2,sistrip::invalid_);
89  rawSpread_ = VFloat(2,sistrip::invalid_);
90  pedsMax_ = VFloat(2,sistrip::invalid_);
91  pedsMin_ = VFloat(2,sistrip::invalid_);
92  noiseMax_ = VFloat(2,sistrip::invalid_);
93  noiseMin_ = VFloat(2,sistrip::invalid_);
94  rawMax_ = VFloat(2,sistrip::invalid_);
95  rawMin_ = VFloat(2,sistrip::invalid_);
96  dead_[0].reserve(256);
97  dead_[1].reserve(256);
98  noisy_[0].reserve(256);
99  noisy_[1].reserve(256);
100 }
101 
102 
103 // ----------------------------------------------------------------------------
104 //
105 uint16_t DaqScopeModeAnalysis::frameFindingThreshold() const {
106  if ( ( getErrorCodes().empty() || getErrorCodes()[0] == "TickMarkRecovered" ) &&
107  base_ < sistrip::valid_ &&
108  peak_ < sistrip::valid_ &&
109  height_ < sistrip::valid_ &&
110  height_ > tickMarkHeightThreshold_ ) {
111  return ( ( static_cast<uint16_t>( base_ +
112  height_ *
113  DaqScopeModeAnalysis::frameFindingThreshold_ ) / 32 ) * 32 );
114  } else { return sistrip::invalid_; }
115 }
116 
117 // ----------------------------------------------------------------------------
118 //
119 bool DaqScopeModeAnalysis::foundTickMark() const {
120  return ( ( getErrorCodes().empty() || getErrorCodes()[0] == "TickMarkRecovered") &&
121  base_ < sistrip::valid_ &&
122  peak_ < sistrip::valid_ &&
123  height_ < sistrip::valid_ &&
124  frameFindingThreshold() < sistrip::valid_ );
125 }
126 
127 // ----------------------------------------------------------------------------
128 //
129 bool DaqScopeModeAnalysis::isValid() const {
130 
131  return ( getErrorCodes().empty() &&
132  base_ < sistrip::valid_ &&
133  peak_ < sistrip::valid_ &&
134  height_ < sistrip::valid_ &&
135  frameFindingThreshold() < sistrip::valid_ &&
136  pedsMean_[0] < sistrip::maximum_ &&
137  pedsMean_[1] < sistrip::maximum_ &&
138  pedsSpread_[0] < sistrip::maximum_ &&
139  pedsSpread_[1] < sistrip::maximum_ &&
140  noiseMean_[0] < sistrip::maximum_ &&
141  noiseMean_[1] < sistrip::maximum_ &&
142  noiseSpread_[0] < sistrip::maximum_ &&
143  noiseSpread_[1] < sistrip::maximum_ &&
144  rawMean_[0] < sistrip::maximum_ &&
145  rawMean_[1] < sistrip::maximum_ &&
146  rawSpread_[0] < sistrip::maximum_ &&
147  rawSpread_[1] < sistrip::maximum_ &&
148  pedsMax_[0] < sistrip::maximum_ &&
149  pedsMax_[1] < sistrip::maximum_ &&
150  pedsMin_[0] < sistrip::maximum_ &&
151  pedsMin_[1] < sistrip::maximum_ &&
152  noiseMax_[0] < sistrip::maximum_ &&
153  noiseMax_[1] < sistrip::maximum_ &&
154  noiseMin_[0] < sistrip::maximum_ &&
155  noiseMin_[1] < sistrip::maximum_ &&
156  rawMax_[0] < sistrip::maximum_ &&
157  rawMax_[1] < sistrip::maximum_ &&
158  rawMin_[0] < sistrip::maximum_ &&
159  rawMin_[1] < sistrip::maximum_ );
160 }
161 
162 // ----------------------------------------------------------------------------
163 //
164 void DaqScopeModeAnalysis::print( std::stringstream& ss, uint32_t iapv ) {
165 
166 
167  if ( iapv == 1 || iapv == 2 ) { iapv--; }
168  else { iapv = 0; }
169 
170  header( ss );
171  ss << std::fixed << std::setprecision(2)
172  << " Tick mark bottom (baseline) [ADC] : " << base_ << std::endl
173  << " Tick mark top [ADC] : " << peak_ << std::endl
174  << " Tick mark height [ADC] : " << height_ << std::endl
175  << " Frame finding threshold [ADC] : " << frameFindingThreshold() << std::endl
176  << std::boolalpha
177  << " Tick mark found : " << foundTickMark() << std::endl
178  << " isValid : " << isValid() << std::endl;
179 
180  if ( peds_[iapv].size() < 128 ||
181  noise_[iapv].size() < 128 ||
182  raw_[iapv].size() < 128 ) {
183 
184  edm::LogWarning(mlCommissioning_)
185  << "[" << myName() << "::" << __func__ << "]"
186  << " Unexpected number of pedestal/noise values: "
187  << peds_[iapv].size() << ", "
188  << noise_[iapv].size() << ", "
189  << raw_[iapv].size();
190  return;
191  }
192 
193  ss << " Monitorables for APV number : " << iapv;
194  if ( iapv == 0 ) { ss << " (first of pair)"; }
195  else if ( iapv == 1 ) { ss << " (second of pair)"; }
196  ss << std::endl;
197  ss << std::fixed << std::setprecision(2);
198  ss << " Example peds/noise for strips : "
199  << " 0, 31, 63, 127" << std::endl
200  << " Peds [ADC] : "
201  << std::setw(6) << peds_[iapv][0] << ", "
202  << std::setw(6) << peds_[iapv][31] << ", "
203  << std::setw(6) << peds_[iapv][63] << ", "
204  << std::setw(6) << peds_[iapv][127] << std::endl
205  << " Noise [ADC] : "
206  << std::setw(6) << noise_[iapv][0] << ", "
207  << std::setw(6) << noise_[iapv][31] << ", "
208  << std::setw(6) << noise_[iapv][63] << ", "
209  << std::setw(6) << noise_[iapv][127] << std::endl
210  << " Raw noise [ADC] : "
211  << std::setw(6) << raw_[iapv][0] << ", "
212  << std::setw(6) << raw_[iapv][31] << ", "
213  << std::setw(6) << raw_[iapv][63] << ", "
214  << std::setw(6) << raw_[iapv][127] << std::endl
215  << " Dead strips (<5s) [strip] : (" << dead_[iapv].size() << " in total) ";
216  for ( uint16_t ii = 0; ii < dead_[iapv].size(); ii++ ) {
217  ss << dead_[iapv][ii] << " "; }
218 
219  ss << std::endl;
220  ss << " Noisy strips (>5s) [strip] : (" << noisy_[iapv].size() << " in total) ";
221  for ( uint16_t ii = 0; ii < noisy_[iapv].size(); ii++ ) {
222  ss << noisy_[iapv][ii] << " ";
223  }
224  ss << std::endl;
225  ss << " Mean peds +/- spread [ADC] : " << pedsMean_[iapv] << " +/- " << pedsSpread_[iapv] << std::endl
226  << " Min/Max pedestal [ADC] : " << pedsMin_[iapv] << " <-> " << pedsMax_[iapv] << std::endl
227  << " Mean noise +/- spread [ADC] : " << noiseMean_[iapv] << " +/- " << noiseSpread_[iapv] << std::endl
228  << " Min/Max noise [ADC] : " << noiseMin_[iapv] << " <-> " << noiseMax_[iapv] << std::endl
229  << " Mean raw noise +/- spread [ADC] : " << rawMean_[iapv] << " +/- " << rawSpread_[iapv] << std::endl
230  << " Min/Max raw noise [ADC] : " << rawMin_[iapv] << " <-> " << rawMax_[iapv] << std::endl
231  << " Normalised noise : " << "(yet to be implemented...)" << std::endl
232  << std::boolalpha
233  << " isValid : " << isValid() << std::endl
234  << std::noboolalpha
235  << " Error codes (found "
236  << std::setw(2) << std::setfill(' ') << getErrorCodes().size()
237  << ") : ";
238  if ( getErrorCodes().empty() ) { ss << "(none)"; }
239  else {
240  VString::const_iterator istr = getErrorCodes().begin();
241  VString::const_iterator jstr = getErrorCodes().end();
242  for ( ; istr != jstr; ++istr ) { ss << *istr << " "; }
243  }
244  ss << std::endl;
245 }
246 
247 
248 
249 // ----------------------------------------------------------------------------
250 //
251 void DaqScopeModeAnalysis::summary( std::stringstream& ss ) const {
252 
253  SiStripFecKey fec_key( fecKey() );
254  SiStripFedKey fed_key( fedKey() );
255 
256  sistrip::RunType type = SiStripEnumsAndStrings::runType( myName() );
257 
258  std::stringstream extra1,extra2,extra3,extra4,extra5,extra6;
259  extra1 << sistrip::extrainfo::pedestals_;
260  extra2 << sistrip::extrainfo::rawNoise_;
261  extra3 << sistrip::extrainfo::commonMode_;
262  extra4 << sistrip::extrainfo::scopeModeFrame_;
263  extra5 << sistrip::extrainfo::scopeModeHeaderLow_;
264  extra6 << sistrip::extrainfo::scopeModeHeaderHigh_;
265 
266  std::string title1 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
267  type,
268  sistrip::FED_KEY,
269  fed_key.key(),
270  sistrip::LLD_CHAN,
271  fec_key.lldChan(),
272  extra1.str() ).title();
273  std::string title2 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
274  type,
275  sistrip::FED_KEY,
276  fed_key.key(),
277  sistrip::LLD_CHAN,
278  fec_key.lldChan(),
279  extra2.str() ).title();
280  std::string title3 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
281  type,
282  sistrip::FED_KEY,
283  fed_key.key(),
284  sistrip::APV,
285  SiStripFecKey::i2cAddr( fec_key.lldChan(), true ),
286  extra3.str() ).title();
287  std::string title4 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
288  type,
289  sistrip::FED_KEY,
290  fed_key.key(),
291  sistrip::APV,
292  SiStripFecKey::i2cAddr( fec_key.lldChan(), false ),
293  extra3.str() ).title();
294 
295  std::string title5 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
296  type,
297  sistrip::FED_KEY,
298  fed_key.key(),
299  sistrip::LLD_CHAN,
300  SiStripFecKey::i2cAddr( fec_key.lldChan()),
301  extra4.str() ).title();
302 
303  std::string title6 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
304  type,
305  sistrip::FED_KEY,
306  fed_key.key(),
307  sistrip::LLD_CHAN,
308  SiStripFecKey::i2cAddr( fec_key.lldChan()),
309  extra5.str() ).title();
310 
311  std::string title7 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
312  type,
313  sistrip::FED_KEY,
314  fed_key.key(),
315  sistrip::LLD_CHAN,
316  SiStripFecKey::i2cAddr( fec_key.lldChan()),
317  extra6.str() ).title();
318 
319  ss << " Summary"
320  << ":"
321  << ( isValid() ? "Valid" : "Invalid" )
322  << ":"
323  << sistrip::controlView_ << ":"
324  << fec_key.fecCrate() << "/"
325  << fec_key.fecSlot() << "/"
326  << fec_key.fecRing() << "/"
327  << fec_key.ccuAddr() << "/"
328  << fec_key.ccuChan()
329  << ":"
330  << sistrip::dqmRoot_ << sistrip::dir_
331  << "Collate" << sistrip::dir_
332  << SiStripFecKey( fec_key.fecCrate(),
333  fec_key.fecSlot(),
334  fec_key.fecRing(),
335  fec_key.ccuAddr(),
336  fec_key.ccuChan() ).path()
337  << ":"
338  << title1 << ";" << title2 << ";" << title3 << ";" << title4 << ";" << title5 << ";" << title6 << ";" << title7
339  << std::endl;
340 
341 }
alignBH_cfg.fixed
fixed
Definition: alignBH_cfg.py:54
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
sistrip::extrainfo::noise_
static const char noise_[]
Definition: ConstantsForDqm.h:47
type
type
Definition: HCALResponse.h:21
sistrip::extrainfo::rawNoise_
static const char rawNoise_[]
Definition: ConstantsForDqm.h:46
CommissioningAnalysis::fedKey
const uint32_t & fedKey() const
Definition: CommissioningAnalysis.h:141
SiStripFecKey::fecRing
const uint16_t & fecRing() const
Definition: SiStripFecKey.h:206
sistrip::dir_
static const char dir_[]
Definition: ConstantsForDqm.h:17
SiStripHistoTitle
Utility class that holds histogram title.
Definition: SiStripHistoTitle.h:20
SiStripFedKey
A container class for generic run and event-related info, information required by the commissioning a...
Definition: SiStripFedKey.h:56
sistrip::noiseMean_
static const char noiseMean_[]
Definition: ConstantsForMonitorable.h:83
DaqScopeModeAnalysis::isValid
bool isValid() const override
Definition: DaqScopeModeAnalysis.cc:129
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
DaqScopeModeAnalysis::frameFindingThreshold
uint16_t frameFindingThreshold() const
Definition: DaqScopeModeAnalysis.cc:105
sistrip::noiseMax_
static const char noiseMax_[]
Definition: ConstantsForMonitorable.h:85
SiStripFecKey::lldChan
const uint16_t & lldChan() const
Definition: SiStripFecKey.h:209
sistrip::noiseMin_
static const char noiseMin_[]
Definition: ConstantsForMonitorable.h:86
SiStripFecKey::fecSlot
const uint16_t & fecSlot() const
Definition: SiStripFecKey.h:205
sistrip::valid_
static const uint16_t valid_
Definition: Constants.h:17
SiStripFecKey::i2cAddr
const uint16_t & i2cAddr() const
Definition: SiStripFecKey.h:210
sistrip
sistrip classes
Definition: EnsembleCalibrationLA.cc:11
SiStripKey::key
const uint32_t & key() const
Definition: SiStripKey.h:125
DaqScopeModeAnalysis::tickMarkHeightThreshold_
static const float tickMarkHeightThreshold_
Definition: DaqScopeModeAnalysis.h:80
SiStripFecKey
Utility class that identifies a position within the strip tracker control structure, down to the level of an APV25.
Definition: SiStripFecKey.h:45
SiStripEnumsAndStrings::runType
static std::string runType(const sistrip::RunType &)
Definition: SiStripEnumsAndStrings.cc:26
CommissioningAnalysis::VFloat
std::vector< float > VFloat
Definition: CommissioningAnalysis.h:33
sistrip::mlCommissioning_
static const char mlCommissioning_[]
Definition: ConstantsForLogger.h:15
sistrip::extrainfo::scopeModeHeaderLow_
static const char scopeModeHeaderLow_[]
Definition: ConstantsForDqm.h:64
DaqScopeModeAnalysis::frameFindingThreshold_
static const float frameFindingThreshold_
Definition: DaqScopeModeAnalysis.h:83
sistrip::controlView_
static const char controlView_[]
Definition: ConstantsForView.h:21
sistrip::maximum_
static const uint16_t maximum_
Definition: Constants.h:20
sistrip::dqmRoot_
static const char dqmRoot_[]
Definition: ConstantsForDqm.h:28
CommissioningAnalysis::fecKey
const uint32_t & fecKey() const
Definition: CommissioningAnalysis.h:140
sistrip::noiseSpread_
static const char noiseSpread_[]
Definition: ConstantsForMonitorable.h:84
sistrip::extrainfo::commonMode_
static const char commonMode_[]
Definition: ConstantsForDqm.h:48
DaqScopeModeAnalysis::print
void print(std::stringstream &, uint32_t apv_number=0) override
Definition: DaqScopeModeAnalysis.cc:164
CommissioningAnalysis::VInt
std::vector< uint16_t > VInt
Definition: CommissioningAnalysis.h:37
cuy.ii
ii
Definition: cuy.py:590
SiStripFecKey::fecCrate
const uint16_t & fecCrate() const
Definition: SiStripFecKey.h:204
sistrip::extrainfo::pedestals_
static const char pedestals_[]
Definition: ConstantsForDqm.h:45
sistrip::extrainfo::scopeModeHeaderHigh_
static const char scopeModeHeaderHigh_[]
Definition: ConstantsForDqm.h:65
CommissioningAnalysis::VVInt
std::vector< VInt > VVInt
Definition: CommissioningAnalysis.h:39
CommissioningAnalysis::header
virtual void header(std::stringstream &) const
Definition: CommissioningAnalysis.cc:32
SiStripFecKey::ccuAddr
const uint16_t & ccuAddr() const
Definition: SiStripFecKey.h:207
sistrip::invalid_
static const uint16_t invalid_
Definition: Constants.h:16
sistrip::extrainfo::scopeModeFrame_
static const char scopeModeFrame_[]
Definition: ConstantsForDqm.h:63
SiStripFecKey::ccuChan
const uint16_t & ccuChan() const
Definition: SiStripFecKey.h:208
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
CommissioningAnalysis::getErrorCodes
const VString & getErrorCodes() const
Definition: CommissioningAnalysis.h:138
CommissioningAnalysis::myName
const std::string & myName() const
Definition: CommissioningAnalysis.h:146
DaqScopeModeAnalysis::summary
void summary(std::stringstream &) const override
Definition: DaqScopeModeAnalysis.cc:251
Generated for CMSSW Reference Manual by doxygen 1.8.11