CMS 3D CMS Logo
RunInfoPayloadInspectoHelper.h
Go to the documentation of this file.
1 #ifndef CONDCORE_RUNINFOPLUGINS_RUNINFOPAYLOADINSPECTORHELPER_H
2 #define CONDCORE_RUNINFOPLUGINS_RUNINFOPAYLOADINSPECTORHELPER_H
3 
4 #include <vector>
5 #include <string>
6 #include <ctime>
7 #include "TH1.h"
8 #include "TH2.h"
9 #include "TStyle.h"
10 #include "TPaveText.h"
11 #include "TStyle.h"
12 #include "FWCore/MessageLogger/interface/MessageLogger.h"
13 
14 // for the FED intervals
15 #include "DataFormats/FEDRawData/interface/FEDNumbering.h"
16 
17 #define MK_P(a, b) std::make_pair(a, b)
18 
19 namespace RunInfoPI {
20 
21  enum state { k_fake = 0, k_valid = 1, k_invalid = 2 };
22 
23  // values are taken from https://github.com/cms-sw/cmssw/blob/master/MagneticField/GeomBuilder/plugins/VolumeBasedMagneticFieldESProducerFromDB.cc#L74-L75
24  constexpr std::array<int, 7> nominalCurrents{{-1, 0, 9558, 14416, 16819, 18268, 19262}};
25  constexpr std::array<float, 7> nominalFields{{3.8, 0., 2., 3., 3.5, 3.8, 4.}};
26 
27  // all parameter than can be displayed
28  enum parameters {
29  m_run, // int
30  m_start_time_ll, // long long;
31  m_stop_time_ll, // long long
32  m_start_current, // float
33  m_stop_current, // float
34  m_avg_current, // float
35  m_max_current, // float
36  m_min_current, // float
37  m_run_interval_seconds, // float
38  m_fedIN, // unsigned int
39  m_BField, // float
40  END_OF_TYPES
41  };
42 
43  /************************************************/
44  inline float theBField(const float current) {
45  // logic is taken from https://github.com/cms-sw/cmssw/blob/master/MagneticField/GeomBuilder/plugins/VolumeBasedMagneticFieldESProducerFromDB.cc#L156
46 
47  int i = 0;
48  for (; i < (int)nominalFields.size() - 1; i++) {
49  if (2 * current < nominalCurrents[i] + nominalCurrents[i + 1]) {
50  return nominalFields[i];
51  }
52  }
53  return nominalFields[i];
54  }
55 
56  /************************************************/
57  inline float runDuration(const std::shared_ptr<RunInfo>& payload) {
58  // calculation of the run duration in seconds
59  time_t start_time = payload->m_start_time_ll;
60  ctime(&start_time);
61  time_t end_time = payload->m_stop_time_ll;
62  ctime(&end_time);
63  return difftime(end_time, start_time) / 1.0e+6;
64  }
65 
66  /************************************************/
67  inline std::string runStartTime(const std::shared_ptr<RunInfo>& payload) {
68  const time_t start_time = payload->m_start_time_ll / 1.0e+6;
69  return std::asctime(std::gmtime(&start_time));
70  }
71 
72  /************************************************/
73  inline std::string runEndTime(const std::shared_ptr<RunInfo>& payload) {
74  const time_t end_time = payload->m_stop_time_ll / 1.0e+6;
75  return std::asctime(std::gmtime(&end_time));
76  }
77 
78  /************************************************/
79  inline std::string getStringFromTypeEnum(const parameters& parameter) {
80  switch (parameter) {
81  case m_run:
82  return "run number";
83  case m_start_time_ll:
84  return "start time";
85  case m_stop_time_ll:
86  return "stop time";
87  case m_start_current:
88  return "start current [A]";
89  case m_stop_current:
90  return "stop current [A]";
91  case m_avg_current:
92  return "average current [A]";
93  case m_max_current:
94  return "max current [A]";
95  case m_min_current:
96  return "min current [A]";
97  case m_run_interval_seconds:
98  return "run duration [s]";
99  case m_fedIN:
100  return "n. FEDs";
101  case m_BField:
102  return "B-field intensity [T]";
103  default:
104  return "should never be here";
105  }
106  }
107 
108  /************************************************/
109  inline void reportSummaryMapPalette(TH2* obj) {
110  static int pcol[20];
111 
112  float rgb[20][3];
113 
114  for (int i = 0; i < 20; i++) {
115  if (i < 17) {
116  rgb[i][0] = 0.80 + 0.01 * i;
117  rgb[i][1] = 0.00 + 0.03 * i;
118  rgb[i][2] = 0.00;
119  } else if (i < 19) {
120  rgb[i][0] = 0.80 + 0.01 * i;
121  rgb[i][1] = 0.00 + 0.03 * i + 0.15 + 0.10 * (i - 17);
122  rgb[i][2] = 0.00;
123  } else if (i == 19) {
124  rgb[i][0] = 0.00;
125  rgb[i][1] = 0.80;
126  rgb[i][2] = 0.00;
127  }
128  pcol[i] = TColor::GetColor(rgb[i][0], rgb[i][1], rgb[i][2]);
129  }
130 
131  gStyle->SetPalette(20, pcol);
132 
133  if (obj) {
134  obj->SetMinimum(-1.e-15);
135  obj->SetMaximum(+1.0);
136  obj->SetOption("colz");
137  }
138  }
139 
140  // FED bounds
141  enum DET {
142  SIPIXEL,
143  SISTRIP,
144  PRESHOWER,
145  TOTEMRP_H,
146  CTPPSDIAMOND,
147  TOTEMRP_V,
148  TOTEMRP_T,
149  ECAL,
150  CASTOR,
151  HCAL,
152  CSC,
153  RPC,
154  HCALPHASE1,
155  SIPIXELPHASE1,
156  GEM
157  };
158 
159  using FEDMAP_T = std::map<RunInfoPI::DET, std::pair<int, int> >;
160 
161  inline FEDMAP_T buildFEDBounds() {
162  static RunInfoPI::FEDMAP_T fb;
163  fb.insert(MK_P(SIPIXEL, MK_P(FEDNumbering::MINSiPixelFEDID, FEDNumbering::MAXSiPixelFEDID)));
164  fb.insert(MK_P(SISTRIP, MK_P(FEDNumbering::MINSiStripFEDID, FEDNumbering::MAXSiStripFEDID)));
165  fb.insert(MK_P(PRESHOWER, MK_P(FEDNumbering::MINPreShowerFEDID, FEDNumbering::MAXPreShowerFEDID)));
166  fb.insert(MK_P(TOTEMRP_H, MK_P(FEDNumbering::MINTotemRPHorizontalFEDID, FEDNumbering::MAXTotemRPHorizontalFEDID)));
167  fb.insert(MK_P(CTPPSDIAMOND, MK_P(FEDNumbering::MINCTPPSDiamondFEDID, FEDNumbering::MAXCTPPSDiamondFEDID)));
168  fb.insert(MK_P(TOTEMRP_V, MK_P(FEDNumbering::MINTotemRPVerticalFEDID, FEDNumbering::MAXTotemRPVerticalFEDID)));
169  fb.insert(MK_P(TOTEMRP_T,
170  MK_P(FEDNumbering::MINTotemRPTimingVerticalFEDID, FEDNumbering::MAXTotemRPTimingVerticalFEDID)));
171  fb.insert(MK_P(ECAL, MK_P(FEDNumbering::MINECALFEDID, FEDNumbering::MAXECALFEDID)));
172  fb.insert(MK_P(CASTOR, MK_P(FEDNumbering::MINCASTORFEDID, FEDNumbering::MAXCASTORFEDID)));
173  fb.insert(MK_P(HCAL, MK_P(FEDNumbering::MINHCALFEDID, FEDNumbering::MAXHCALFEDID)));
174  fb.insert(MK_P(CSC, MK_P(FEDNumbering::MINCSCFEDID, FEDNumbering::MAXCSCFEDID)));
175  fb.insert(MK_P(RPC, MK_P(FEDNumbering::MINRPCFEDID, FEDNumbering::MAXRPCFEDID)));
176  fb.insert(MK_P(HCALPHASE1, MK_P(FEDNumbering::MINHCALuTCAFEDID, FEDNumbering::MAXHCALuTCAFEDID)));
177  fb.insert(MK_P(SIPIXELPHASE1, MK_P(FEDNumbering::MINSiPixeluTCAFEDID, FEDNumbering::MAXSiPixeluTCAFEDID)));
178  fb.insert(MK_P(GEM, MK_P(FEDNumbering::MINGEMFEDID, FEDNumbering::MAXGEMFEDID)));
179  return fb;
180  }
181 }; // namespace RunInfoPI
182 #endif
mps_fire.i
i
Definition: mps_fire.py:429
RunInfoPI::reportSummaryMapPalette
void reportSummaryMapPalette(TH2 *obj)
Definition: RunInfoPayloadInspectoHelper.h:109
RunInfoPI::runDuration
float runDuration(const std::shared_ptr< RunInfo > &payload)
Definition: RunInfoPayloadInspectoHelper.h:57
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
RunInfoPI::runEndTime
std::string runEndTime(const std::shared_ptr< RunInfo > &payload)
Definition: RunInfoPayloadInspectoHelper.h:73
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
RunInfoPI::FEDMAP_T
std::map< RunInfoPI::DET, std::pair< int, int > > FEDMAP_T
Definition: RunInfoPayloadInspectoHelper.h:159
RunInfoPI::getStringFromTypeEnum
std::string getStringFromTypeEnum(const parameters &parameter)
Definition: RunInfoPayloadInspectoHelper.h:79
jetsAK4_Puppi_cff.payload
payload
Definition: jetsAK4_Puppi_cff.py:17
HCAL
Definition: HCAL.py:1
pfMETsysShiftCorrections_cfi.parameter
parameter
Definition: pfMETsysShiftCorrections_cfi.py:118
RunInfoPI::nominalFields
constexpr std::array< float, 7 > nominalFields
Definition: RunInfoPayloadInspectoHelper.h:25
RunInfoPI::theBField
float theBField(const float current)
Definition: RunInfoPayloadInspectoHelper.h:44
RunInfoPI::nominalCurrents
constexpr std::array< int, 7 > nominalCurrents
Definition: RunInfoPayloadInspectoHelper.h:24
RunInfoPI::runStartTime
std::string runStartTime(const std::shared_ptr< RunInfo > &payload)
Definition: RunInfoPayloadInspectoHelper.h:67
svgfig.rgb
def rgb(r, g, b, maximum=1.)
Definition: svgfig.py:40
MK_P
#define MK_P(a, b)
Definition: RunInfoPayloadInspectoHelper.h:17
Generated for CMSSW Reference Manual by doxygen 1.8.14