CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
PedsFullNoiseAnalysis.cc
Go to the documentation of this file.
1 #include "CondFormats/SiStripObjects/interface/PedsFullNoiseAnalysis.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 // ----------------------------------------------------------------------------
12 //
13 PedsFullNoiseAnalysis::PedsFullNoiseAnalysis( const uint32_t& key )
14  : CommissioningAnalysis(key,"PedsFullNoiseAnalysis"),
15  peds_(2,VFloat(128,sistrip::invalid_)),
16  noise_(2,VFloat(128,sistrip::invalid_)),
17  raw_(2,VFloat(128,sistrip::invalid_)),
18  ksProb_(2,VFloat(0,sistrip::invalid_)),
19  chi2Prob_(2,VFloat(0,sistrip::invalid_)),
20  noiseGaus_(2,VFloat(0,sistrip::invalid_)),
21  noiseBin84_(2,VFloat(0,sistrip::invalid_)),
22  noiseRMS_(2,VFloat(0,sistrip::invalid_)),
23  noiseSignif_(2,VFloat(0,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  legacy_(false)
39 {
40  dead_[0].reserve(128); dead_[1].reserve(128);
41  noisy_[0].reserve(128); noisy_[1].reserve(128);
42  ksProb_[0].reserve(128); ksProb_[1].reserve(128);
43  chi2Prob_[0].reserve(128); chi2Prob_[1].reserve(128);
44  noiseGaus_[0].reserve(128); noiseGaus_[1].reserve(128);
45  noiseBin84_[0].reserve(128); noiseBin84_[1].reserve(128);
46  noiseRMS_[0].reserve(128); noiseRMS_[1].reserve(128);
47  noiseSignif_[0].reserve(128); noiseSignif_[1].reserve(128);
48 }
49 
50 // ----------------------------------------------------------------------------
51 //
52 PedsFullNoiseAnalysis::PedsFullNoiseAnalysis()
53  : CommissioningAnalysis("PedsFullNoiseAnalysis"),
54  peds_(2,VFloat(128,sistrip::invalid_)),
55  noise_(2,VFloat(128,sistrip::invalid_)),
56  raw_(2,VFloat(128,sistrip::invalid_)),
57  ksProb_(2,VFloat(0,sistrip::invalid_)),
58  chi2Prob_(2,VFloat(0,sistrip::invalid_)),
59  noiseGaus_(2,VFloat(0,sistrip::invalid_)),
60  noiseBin84_(2,VFloat(0,sistrip::invalid_)),
61  noiseRMS_(2,VFloat(0,sistrip::invalid_)),
62  noiseSignif_(2,VFloat(0,sistrip::invalid_)),
63  dead_(2,VInt(0,sistrip::invalid_)),
64  noisy_(2,VInt(0,sistrip::invalid_)),
65  pedsMean_(2,sistrip::invalid_),
66  pedsSpread_(2,sistrip::invalid_),
67  noiseMean_(2,sistrip::invalid_),
68  noiseSpread_(2,sistrip::invalid_),
69  rawMean_(2,sistrip::invalid_),
70  rawSpread_(2,sistrip::invalid_),
71  pedsMax_(2,sistrip::invalid_),
72  pedsMin_(2,sistrip::invalid_),
73  noiseMax_(2,sistrip::invalid_),
74  noiseMin_(2,sistrip::invalid_),
75  rawMax_(2,sistrip::invalid_),
76  rawMin_(2,sistrip::invalid_),
77  legacy_(false)
78 {
79  dead_[0].reserve(128); dead_[1].reserve(128);
80  noisy_[0].reserve(128); noisy_[1].reserve(128);
81  ksProb_[0].reserve(128); ksProb_[1].reserve(128);
82  chi2Prob_[0].reserve(128); chi2Prob_[1].reserve(128);
83  noiseGaus_[0].reserve(128); noiseGaus_[1].reserve(128);
84  noiseBin84_[0].reserve(128); noiseBin84_[1].reserve(128);
85  noiseRMS_[0].reserve(128); noiseRMS_[1].reserve(128);
86  noiseSignif_[0].reserve(128); noiseSignif_[1].reserve(128);
87 }
88 
89 // ----------------------------------------------------------------------------
90 //
91 void PedsFullNoiseAnalysis::reset() {
92  peds_ = VVFloat(2,VFloat(128,sistrip::invalid_));
93  noise_ = VVFloat(2,VFloat(128,sistrip::invalid_));
94  raw_ = VVFloat(2,VFloat(128,sistrip::invalid_));
95  ksProb_ = VVFloat(2,VFloat(0,sistrip::invalid_));
96  chi2Prob_ = VVFloat(2,VFloat(0,sistrip::invalid_));
97  noiseGaus_ = VVFloat(2,VFloat(0,sistrip::invalid_));
98  noiseBin84_ = VVFloat(2,VFloat(0,sistrip::invalid_));
99  noiseRMS_ = VVFloat(2,VFloat(0,sistrip::invalid_));
100  noiseSignif_ = VVFloat(2,VFloat(0,sistrip::invalid_));
101  dead_ = VVInt(2,VInt(0,sistrip::invalid_));
102  noisy_ = VVInt(2,VInt(0,sistrip::invalid_));
103  pedsMean_ = VFloat(2,sistrip::invalid_);
104  pedsSpread_ = VFloat(2,sistrip::invalid_);
105  noiseMean_ = VFloat(2,sistrip::invalid_);
106  noiseSpread_ = VFloat(2,sistrip::invalid_);
107  rawMean_ = VFloat(2,sistrip::invalid_);
108  rawSpread_ = VFloat(2,sistrip::invalid_);
109  pedsMax_ = VFloat(2,sistrip::invalid_);
110  pedsMin_ = VFloat(2,sistrip::invalid_);
111  noiseMax_ = VFloat(2,sistrip::invalid_);
112  noiseMin_ = VFloat(2,sistrip::invalid_);
113  rawMax_ = VFloat(2,sistrip::invalid_);
114  rawMin_ = VFloat(2,sistrip::invalid_);
115  dead_[0].reserve(128); dead_[1].reserve(128);
116  noisy_[0].reserve(128); noisy_[1].reserve(128);
117  ksProb_[0].reserve(128); ksProb_[1].reserve(128);
118  chi2Prob_[0].reserve(128); chi2Prob_[1].reserve(128);
119  noiseGaus_[0].reserve(128); noiseGaus_[1].reserve(128);
120  noiseBin84_[0].reserve(128); noiseBin84_[1].reserve(128);
121  noiseRMS_[0].reserve(128); noiseRMS_[1].reserve(128);
122  noiseSignif_[0].reserve(128); noiseSignif_[1].reserve(128);
123  legacy_ = false;
124 }
125 
126 // ----------------------------------------------------------------------------
127 //
128 bool PedsFullNoiseAnalysis::isValid() const {
129  return ( pedsMean_[0] < sistrip::maximum_ &&
130  pedsMean_[1] < sistrip::maximum_ &&
131  pedsSpread_[0] < sistrip::maximum_ &&
132  pedsSpread_[1] < sistrip::maximum_ &&
133  noiseMean_[0] < sistrip::maximum_ &&
134  noiseMean_[1] < sistrip::maximum_ &&
135  noiseSpread_[0] < sistrip::maximum_ &&
136  noiseSpread_[1] < sistrip::maximum_ &&
137  rawMean_[0] < sistrip::maximum_ &&
138  rawMean_[1] < sistrip::maximum_ &&
139  rawSpread_[0] < sistrip::maximum_ &&
140  rawSpread_[1] < sistrip::maximum_ &&
141  pedsMax_[0] < sistrip::maximum_ &&
142  pedsMax_[1] < sistrip::maximum_ &&
143  pedsMin_[0] < sistrip::maximum_ &&
144  pedsMin_[1] < sistrip::maximum_ &&
145  noiseMax_[0] < sistrip::maximum_ &&
146  noiseMax_[1] < sistrip::maximum_ &&
147  noiseMin_[0] < sistrip::maximum_ &&
148  noiseMin_[1] < sistrip::maximum_ &&
149  rawMax_[0] < sistrip::maximum_ &&
150  rawMax_[1] < sistrip::maximum_ &&
151  rawMin_[0] < sistrip::maximum_ &&
152  rawMin_[1] < sistrip::maximum_ &&
153  getErrorCodes().empty() );
154  //noiseMean_[0] <= rawMean_[0] && //@@ temp
155  //noiseMean_[1] <= rawMean_[1] ); //@@ temp
156 }
157 
158 // ----------------------------------------------------------------------------
159 //
160 void PedsFullNoiseAnalysis::summary( std::stringstream& ss ) const {
161 
162  SiStripFecKey fec_key( fecKey() );
163  SiStripFedKey fed_key( fedKey() );
164 
165  sistrip::RunType type = SiStripEnumsAndStrings::runType( myName() );
166 
167  std::stringstream extra1,extra2,extra3;
168  if ( legacy_ ) {
169  extra1 << sistrip::extrainfo::pedsAndRawNoise_;
170  extra2 << sistrip::extrainfo::pedsAndCmSubNoise_;
171  extra3 << sistrip::extrainfo::commonMode_;
172  } else {
173  extra1 << sistrip::extrainfo::pedestals_;
174  extra2 << sistrip::extrainfo::rawNoise_;
175  extra3 << sistrip::extrainfo::commonMode_;
176  }
177 
178  std::string title1 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
179  type,
180  sistrip::FED_KEY,
181  fed_key.key(),
182  sistrip::LLD_CHAN,
183  fec_key.lldChan(),
184  extra1.str() ).title();
185  std::string title2 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
186  type,
187  sistrip::FED_KEY,
188  fed_key.key(),
189  sistrip::LLD_CHAN,
190  fec_key.lldChan(),
191  extra2.str() ).title();
192  std::string title3 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
193  type,
194  sistrip::FED_KEY,
195  fed_key.key(),
196  sistrip::APV,
197  SiStripFecKey::i2cAddr( fec_key.lldChan(), true ),
198  extra3.str() ).title();
199  std::string title4 = SiStripHistoTitle( sistrip::EXPERT_HISTO,
200  type,
201  sistrip::FED_KEY,
202  fed_key.key(),
203  sistrip::APV,
204  SiStripFecKey::i2cAddr( fec_key.lldChan(), false ),
205  extra3.str() ).title();
206 
207 
208  ss << " Summary"
209  << ":"
210  << ( isValid() ? "Valid" : "Invalid" )
211  << ":"
212  << sistrip::controlView_ << ":"
213  << fec_key.fecCrate() << "/"
214  << fec_key.fecSlot() << "/"
215  << fec_key.fecRing() << "/"
216  << fec_key.ccuAddr() << "/"
217  << fec_key.ccuChan()
218  << ":"
219  << sistrip::dqmRoot_ << sistrip::dir_
220  << "Collate" << sistrip::dir_
221  << SiStripFecKey( fec_key.fecCrate(),
222  fec_key.fecSlot(),
223  fec_key.fecRing(),
224  fec_key.ccuAddr(),
225  fec_key.ccuChan() ).path()
226  << ":"
227  << title1 << ";" << title2 << ";" << title3 << ";" << title4
228  << std::endl;
229 
230 }
231 
232 // ----------------------------------------------------------------------------
233 //
234 void PedsFullNoiseAnalysis::print( std::stringstream& ss, uint32_t iapv ) {
235 
236  if ( iapv == 1 || iapv == 2 ) { iapv--; }
237  else { iapv = 0; }
238 
239  if ( peds_[iapv].size() < 128 ||
240  noise_[iapv].size() < 128 ||
241  raw_[iapv].size() < 128 ) {
242  edm::LogWarning(mlCommissioning_)
243  << "[" << myName() << "::" << __func__ << "]"
244  << " Unexpected number of pedestal/noise values: "
245  << peds_[iapv].size() << ", "
246  << noise_[iapv].size() << ", "
247  << raw_[iapv].size();
248  return;
249  }
250 
251  header( ss );
252  ss << " Monitorables for APV number : " << iapv;
253  if ( iapv == 0 ) { ss << " (first of pair)"; }
254  else if ( iapv == 1 ) { ss << " (second of pair)"; }
255  ss << std::endl;
256  ss << std::fixed << std::setprecision(2);
257  ss << " Example peds/noise for strips : "
258  << " 0, 31, 63, 127" << std::endl
259  << " Peds [ADC] : "
260  << std::setw(6) << peds_[iapv][0] << ", "
261  << std::setw(6) << peds_[iapv][31] << ", "
262  << std::setw(6) << peds_[iapv][63] << ", "
263  << std::setw(6) << peds_[iapv][127] << std::endl
264  << " Noise [ADC] : "
265  << std::setw(6) << noise_[iapv][0] << ", "
266  << std::setw(6) << noise_[iapv][31] << ", "
267  << std::setw(6) << noise_[iapv][63] << ", "
268  << std::setw(6) << noise_[iapv][127] << std::endl
269  << " Raw noise [ADC] : "
270  << std::setw(6) << raw_[iapv][0] << ", "
271  << std::setw(6) << raw_[iapv][31] << ", "
272  << std::setw(6) << raw_[iapv][63] << ", "
273  << std::setw(6) << raw_[iapv][127] << std::endl
274  << " Dead strips (<5s) [strip] : (" << dead_[iapv].size() << " in total) ";
275  for ( uint16_t ii = 0; ii < dead_[iapv].size(); ii++ ) {
276  ss << dead_[iapv][ii] << " "; }
277 
278  ss << std::endl;
279  ss << " Noisy strips (>5s) [strip] : (" << noisy_[iapv].size() << " in total) ";
280  for ( uint16_t ii = 0; ii < noisy_[iapv].size(); ii++ ) {
281  ss << noisy_[iapv][ii] << " ";
282  }
283  ss << std::endl;
284  ss << " Mean peds +/- spread [ADC] : " << pedsMean_[iapv] << " +/- " << pedsSpread_[iapv] << std::endl
285  << " Min/Max pedestal [ADC] : " << pedsMin_[iapv] << " <-> " << pedsMax_[iapv] << std::endl
286  << " Mean noise +/- spread [ADC] : " << noiseMean_[iapv] << " +/- " << noiseSpread_[iapv] << std::endl
287  << " Min/Max noise [ADC] : " << noiseMin_[iapv] << " <-> " << noiseMax_[iapv] << std::endl
288  << " Mean raw noise +/- spread [ADC] : " << rawMean_[iapv] << " +/- " << rawSpread_[iapv] << std::endl
289  << " Min/Max raw noise [ADC] : " << rawMin_[iapv] << " <-> " << rawMax_[iapv] << std::endl
290  << " Normalised noise : " << "(yet to be implemented...)" << std::endl
291  << std::boolalpha
292  << " isValid : " << isValid() << std::endl
293  << std::noboolalpha
294  << " Error codes (found "
295  << std::setw(2) << std::setfill(' ') << getErrorCodes().size()
296  << ") : ";
297  if ( getErrorCodes().empty() ) { ss << "(none)"; }
298  else {
299  VString::const_iterator istr = getErrorCodes().begin();
300  VString::const_iterator jstr = getErrorCodes().end();
301  for ( ; istr != jstr; ++istr ) { ss << *istr << " "; }
302  }
303  ss << std::endl;
304 }
305 
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:136
SiStripFecKey::fecRing
const uint16_t & fecRing() const
Definition: SiStripFecKey.h:206
sistrip::extrainfo::pedsAndRawNoise_
static const char pedsAndRawNoise_[]
Definition: ConstantsForDqm.h:53
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:84
sistrip::noiseMax_
static const char noiseMax_[]
Definition: ConstantsForMonitorable.h:86
SiStripFecKey::lldChan
const uint16_t & lldChan() const
Definition: SiStripFecKey.h:209
sistrip::noiseMin_
static const char noiseMin_[]
Definition: ConstantsForMonitorable.h:87
SiStripFecKey::fecSlot
const uint16_t & fecSlot() const
Definition: SiStripFecKey.h:205
cuy.ii
int ii
Definition: cuy.py:588
SiStripFecKey::i2cAddr
const uint16_t & i2cAddr() const
Definition: SiStripFecKey.h:210
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
SiStripKey::key
const uint32_t & key() const
Definition: SiStripKey.h:125
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:27
CommissioningAnalysis::VFloat
std::vector< float > VFloat
Definition: CommissioningAnalysis.h:33
cmsHarvester.path
tuple path
else: Piece not in the list, fine.
Definition: cmsHarvester.py:1591
sistrip::mlCommissioning_
static const char mlCommissioning_[]
Definition: ConstantsForLogger.h:15
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:135
sistrip::noiseSpread_
static const char noiseSpread_[]
Definition: ConstantsForMonitorable.h:85
sistrip::extrainfo::commonMode_
static const char commonMode_[]
Definition: ConstantsForDqm.h:48
relval_steps.key
string key
FastSim: produces sample of signal events, overlayed with premixed minbias events.
Definition: relval_steps.py:700
CommissioningAnalysis::VInt
std::vector< uint16_t > VInt
Definition: CommissioningAnalysis.h:37
SiStripFecKey::fecCrate
const uint16_t & fecCrate() const
Definition: SiStripFecKey.h:204
sistrip::extrainfo::pedestals_
static const char pedestals_[]
Definition: ConstantsForDqm.h:45
PedsFullNoiseAnalysis::summary
void summary(std::stringstream &) const
Definition: PedsFullNoiseAnalysis.cc:160
sistrip::extrainfo::pedsAndCmSubNoise_
static const char pedsAndCmSubNoise_[]
Definition: ConstantsForDqm.h:55
CommissioningAnalysis::VVInt
std::vector< VInt > VVInt
Definition: CommissioningAnalysis.h:39
PedsFullNoiseAnalysis::print
void print(std::stringstream &, uint32_t apv_number=0)
Definition: PedsFullNoiseAnalysis.cc:234
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
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
edm::false
volatile std::atomic< bool > shutdown_flag false
Definition: UnixSignalHandlers.cc:22
CommissioningAnalysis::getErrorCodes
const VString & getErrorCodes() const
Definition: CommissioningAnalysis.h:133
CommissioningAnalysis::myName
const std::string & myName() const
Definition: CommissioningAnalysis.h:141
findQualityFiles.size
tuple size
Write out results.
Definition: findQualityFiles.py:442
sistrip::noiseSignif_
static const char noiseSignif_[]
Definition: ConstantsForMonitorable.h:95
Generated for CMSSW Reference Manual by doxygen 1.8.5