CMS 3D CMS Logo

List of all members | Public Member Functions | Private Member Functions | Private Attributes
SiStripTkMaps Class Reference

#include <SiStripTkMaps.h>

Public Member Functions

void bookMap (const std::string mapTitle, const std::string zAxisTitle)
 
void drawMap (TCanvas &canvas, std::string option="")
 
void fill (long rawid, double val)
 
const std::vector< unsigned int > & getTheFilledIds ()
 
const std::vector< double > & getTheFilledValues ()
 
const TH2Poly * getTheMap ()
 
const std::string & getTheMapTitle ()
 
const std::string & getTheZAxisTitle ()
 
void setZAxisRange (double xmin, double xmax)
 
 SiStripTkMaps (const char *option)
 
 ~SiStripTkMaps ()=default
 

Private Member Functions

void adjustCanvasMargins (TVirtualPad *pad, const float top, const float bottom, const float left, const float right)
 
void drawArrows (const float x_X1, const float x_X2, const float x_Y1, const float y_Y2, const char *x_label, const char *y_label)
 
void dressMap (TCanvas &canv)
 
void readVertices (double &minx, double &maxx, double &miny, double &maxy)
 

Private Attributes

double m_axmax
 
double m_axmin
 
std::map< long, std::shared_ptr< TGraph > > m_bins
 
std::vector< unsigned int > m_detIdVector
 
std::string m_mapTitle = ""
 
Option_t * m_option
 
TH2Poly * m_trackerMap {nullptr}
 
TrackerTopology m_trackerTopo
 
std::vector< double > m_values
 
std::string m_zAxisTitle = ""
 

Detailed Description

Definition at line 40 of file SiStripTkMaps.h.

Constructor & Destructor Documentation

◆ SiStripTkMaps()

SiStripTkMaps::SiStripTkMaps ( const char *  option)
inline

Definition at line 42 of file SiStripTkMaps.h.

References fileinputsource_cfi::option.

43  : m_option{option},
45  edm::FileInPath("Geometry/TrackerCommonData/data/PhaseI/trackerParameters.xml").fullPath())} {}
Option_t * m_option
Definition: SiStripTkMaps.h:73
TrackerTopology fromTrackerParametersXMLFile(const std::string &xmlFileName)
TrackerTopology m_trackerTopo
Definition: SiStripTkMaps.h:80

◆ ~SiStripTkMaps()

SiStripTkMaps::~SiStripTkMaps ( )
default

Member Function Documentation

◆ adjustCanvasMargins()

void SiStripTkMaps::adjustCanvasMargins ( TVirtualPad *  pad,
const float  top,
const float  bottom,
const float  left,
const float  right 
)
private

Definition at line 207 of file SiStripTkMaps.cc.

Referenced by drawMap().

208  {
209  if (top > 0) {
210  pad->SetTopMargin(top);
211  }
212  if (bottom > 0) {
213  pad->SetBottomMargin(bottom);
214  }
215  if (left > 0) {
216  pad->SetLeftMargin(left);
217  }
218  if (right > 0) {
219  pad->SetRightMargin(right);
220  }
221 }

◆ bookMap()

void SiStripTkMaps::bookMap ( const std::string  mapTitle,
const std::string  zAxisTitle 
)

Definition at line 36 of file SiStripTkMaps.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), m_bins, m_mapTitle, m_option, m_trackerMap, m_zAxisTitle, margin, and readVertices().

36  {
37  double minx = 0xFFFFFF, maxx = -0xFFFFFF, miny = 0xFFFFFF, maxy = -0xFFFFFFF;
38  readVertices(minx, maxx, miny, maxy);
39 
40  // set the titles
41  m_zAxisTitle = zAxisTitle;
42  m_mapTitle = mapTitle;
43 
44  TGaxis::SetMaxDigits(2);
45 
46  // margin of the box
47  static constexpr int margin = 5;
48  m_trackerMap =
49  new TH2Poly("Tracker Map", m_mapTitle.c_str(), minx - margin, maxx + margin, miny - margin, maxy + margin);
50  m_trackerMap->SetFloat();
51  m_trackerMap->SetOption(m_option);
52  m_trackerMap->SetStats(false);
53  m_trackerMap->GetZaxis()->SetLabelSize(0.03);
54  m_trackerMap->GetZaxis()->SetTitleOffset(0.5);
55  m_trackerMap->GetZaxis()->SetTitleSize(0.05);
56  m_trackerMap->GetZaxis()->SetTitle(m_zAxisTitle.c_str());
57  m_trackerMap->GetZaxis()->CenterTitle();
58 
59  for (const auto& pair : m_bins) {
60  m_trackerMap->AddBin(pair.second->Clone());
61  }
62 }
void readVertices(double &minx, double &maxx, double &miny, double &maxy)
Option_t * m_option
Definition: SiStripTkMaps.h:73
std::string m_zAxisTitle
Definition: SiStripTkMaps.h:75
Double_t margin
std::string m_mapTitle
Definition: SiStripTkMaps.h:74
TH2Poly * m_trackerMap
Definition: SiStripTkMaps.h:81
std::map< long, std::shared_ptr< TGraph > > m_bins
Definition: SiStripTkMaps.h:77

◆ drawArrows()

void SiStripTkMaps::drawArrows ( const float  x_X1,
const float  x_X2,
const float  x_Y1,
const float  y_Y2,
const char *  x_label,
const char *  y_label 
)
private

Definition at line 177 of file SiStripTkMaps.cc.

Referenced by dressMap().

178  {
179  auto arrow_X = TArrow();
180  arrow_X.SetLineColor(kBlue);
181  arrow_X.SetLineWidth(2);
182  arrow_X.SetOption("|>");
183  arrow_X.SetArrowSize(10);
184  arrow_X.DrawLineNDC(x_X1, x_Y1, x_X2, x_Y1);
185 
186  auto arrow_Y = TArrow();
187  arrow_Y.SetLineColor(kBlue);
188  arrow_Y.SetLineWidth(2);
189  arrow_Y.SetOption("|>");
190  arrow_Y.SetArrowSize(10);
191  arrow_Y.DrawLineNDC(x_X2, x_Y1, x_X2, y_Y2);
192 
193  auto text_X = TLatex();
194  text_X.SetTextSize(0.04);
195  text_X.SetTextAlign(11);
196  text_X.SetTextColor(kBlue);
197  text_X.DrawLatexNDC(x_X1, x_Y1 - 0.03, x_label);
198 
199  auto text_Y = TLatex();
200  text_Y.SetTextSize(0.04);
201  text_Y.SetTextAlign(11);
202  text_Y.SetTextColor(kBlue);
203  text_Y.DrawLatexNDC(x_X2 + 0.005, y_Y2 - 0.01, y_label);
204 }

◆ drawMap()

void SiStripTkMaps::drawMap ( TCanvas &  canvas,
std::string  option = "" 
)

Definition at line 71 of file SiStripTkMaps.cc.

References adjustCanvasMargins(), svgfig::canvas(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), dressMap(), m_trackerMap, and fileinputsource_cfi::option.

71  {
72  // margins
73  static constexpr float tmargin_ = 0.08;
74  static constexpr float bmargin_ = 0.02;
75  static constexpr float lmargin_ = 0.02;
76  static constexpr float rmargin_ = 0.08;
77 
78  // window size
79  static constexpr int wH_ = 3000;
80  static constexpr int hH_ = 700;
81 
82  canvas.cd();
83  adjustCanvasMargins(canvas.cd(), tmargin_, bmargin_, lmargin_, rmargin_);
84  canvas.Update();
85 
86  m_trackerMap->SetTitle("");
87  if (!option.empty()) {
88  m_trackerMap->Draw(option.c_str());
89  } else {
90  m_trackerMap->Draw();
91  }
92 
93  canvas.SetFrameLineColor(0);
94  gPad->Update();
95  TPaletteAxis* palette = (TPaletteAxis*)m_trackerMap->GetListOfFunctions()->FindObject("palette");
96  if (palette != nullptr) {
97  palette->SetLabelSize(0.02);
98  palette->SetX1NDC(1 - rmargin_);
99  palette->SetX2NDC(1 - rmargin_ + lmargin_);
100  }
101 
102  // if not right size, and not drawn in same mode
103  if (canvas.GetWindowHeight() != hH_ && canvas.GetWindowWidth() != wH_ && option.find("same") == std::string::npos) {
104  canvas.SetWindowSize(wH_, hH_);
105  }
106 
107  // call the map dressing
108  dressMap(canvas);
109 }
void dressMap(TCanvas &canv)
void adjustCanvasMargins(TVirtualPad *pad, const float top, const float bottom, const float left, const float right)
def canvas(sub, attr)
Definition: svgfig.py:482
TH2Poly * m_trackerMap
Definition: SiStripTkMaps.h:81

◆ dressMap()

void SiStripTkMaps::dressMap ( TCanvas &  canv)
private

Definition at line 112 of file SiStripTkMaps.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), drawArrows(), m_mapTitle, and Skims_PA_cff::name.

Referenced by drawMap().

112  {
113  std::array<std::string, 12> barrelNames = {
114  {"TIB L2", "TIB L1", "TIB L4", "TIB L3", "TOB L2", "TOB L1", "TOB L4", " TOB L3", "TOB L6", "TOB L5"}};
115  std::array<std::string, 4> endcapNames = {{"TID", "TEC", "TID", "TEC"}};
116  std::array<std::string, 24> disknumbering = {{"+1", "+2", "+3", "+1", "+2", "+3", "+4", "+5",
117  "+6", "+7", "+8", "+9", "-1", "-2", "-3", "-1",
118  "-2", "-3", "-4", "-5", "-6", "-7", "-8", "-9"}};
119 
120  static constexpr std::array<float, 12> b_coordx = {
121  {0.1, 0.1, 0.26, 0.26, 0.41, 0.41, 0.56, 0.56, 0.725, 0.725, 0.05, 0.17}};
122  static constexpr std::array<float, 12> b_coordy = {
123  {0.70, 0.45, 0.70, 0.45, 0.70, 0.46, 0.70, 0.46, 0.70, 0.46, 0.85, 0.85}};
124 
125  static constexpr std::array<float, 4> e_coordx = {{0.01, 0.21, 0.01, 0.21}};
126  static constexpr std::array<float, 4> e_coordy = {{0.89, 0.89, 0.17, 0.17}};
127 
128  static constexpr std::array<float, 24> n_coordx = {{0.01, 0.087, 0.165, 0.227, 0.305, 0.383, 0.461, 0.539,
129  0.616, 0.694, 0.772, 0.850, 0.01, 0.087, 0.165, 0.227,
130  0.305, 0.383, 0.461, 0.539, 0.617, 0.695, 0.773, 0.851}};
131 
132  static constexpr std::array<float, 24> n_coordy = {{0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85, 0.85,
133  0.85, 0.85, 0.85, 0.85, 0.13, 0.13, 0.13, 0.13,
134  0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13}};
135 
136  canv.cd();
137  for (const auto& name : barrelNames | boost::adaptors::indexed(0)) {
138  auto ltx = TLatex();
139  ltx.SetTextFont(62);
140  ltx.SetTextSize(0.035);
141  ltx.SetTextAlign(11);
142  ltx.DrawLatexNDC(b_coordx[name.index()], b_coordy[name.index()], name.value().c_str());
143  }
144 
145  for (const auto& name : endcapNames | boost::adaptors::indexed(0)) {
146  auto ltx = TLatex();
147  ltx.SetTextFont(62);
148  ltx.SetTextSize(0.05);
149  ltx.SetTextAlign(11);
150  ltx.DrawLatexNDC(e_coordx[name.index()], e_coordy[name.index()], name.value().c_str());
151  }
152 
153  for (const auto& name : disknumbering | boost::adaptors::indexed(0)) {
154  auto ltx = TLatex();
155  ltx.SetTextFont(62);
156  ltx.SetTextSize(0.035);
157  ltx.SetTextAlign(11);
158  ltx.DrawLatexNDC(n_coordx[name.index()], n_coordy[name.index()], name.value().c_str());
159  }
160 
161  auto ltx = TLatex();
162  ltx.SetTextFont(62);
163  ltx.SetTextSize(0.045);
164  ltx.SetTextAlign(11);
165  ltx.DrawLatexNDC(gPad->GetLeftMargin(), 1 - gPad->GetTopMargin() + 0.03, m_mapTitle.c_str());
166 
167  // barrel axes
168  drawArrows(0.09, 0.23, 0.24, 0.45, "#phi", "z");
169  // endcap axes
170  drawArrows(0.85, 0.89, 0.83, 0.95, "x", "y");
171 
172  canv.Modified();
173  canv.Update(); // make sure it's really (re)drawn
174 }
void drawArrows(const float x_X1, const float x_X2, const float x_Y1, const float y_Y2, const char *x_label, const char *y_label)
std::string m_mapTitle
Definition: SiStripTkMaps.h:74

◆ fill()

void SiStripTkMaps::fill ( long  rawid,
double  val 
)

Definition at line 65 of file SiStripTkMaps.cc.

References m_trackerMap, m_values, and heppy_batch::val.

65  {
66  m_trackerMap->Fill(TString::Format("%ld", rawid), val);
67  m_values.push_back(val);
68 }
std::vector< double > m_values
Definition: SiStripTkMaps.h:79
TH2Poly * m_trackerMap
Definition: SiStripTkMaps.h:81

◆ getTheFilledIds()

const std::vector<unsigned int>& SiStripTkMaps::getTheFilledIds ( )
inline

Definition at line 63 of file SiStripTkMaps.h.

References m_detIdVector.

63 { return m_detIdVector; }
std::vector< unsigned int > m_detIdVector
Definition: SiStripTkMaps.h:78

◆ getTheFilledValues()

const std::vector<double>& SiStripTkMaps::getTheFilledValues ( )
inline

Definition at line 66 of file SiStripTkMaps.h.

References m_values.

66 { return m_values; }
std::vector< double > m_values
Definition: SiStripTkMaps.h:79

◆ getTheMap()

const TH2Poly* SiStripTkMaps::getTheMap ( )
inline

Definition at line 54 of file SiStripTkMaps.h.

References m_trackerMap.

54 { return m_trackerMap; }
TH2Poly * m_trackerMap
Definition: SiStripTkMaps.h:81

◆ getTheMapTitle()

const std::string& SiStripTkMaps::getTheMapTitle ( )
inline

Definition at line 57 of file SiStripTkMaps.h.

References m_mapTitle.

57 { return m_mapTitle; }
std::string m_mapTitle
Definition: SiStripTkMaps.h:74

◆ getTheZAxisTitle()

const std::string& SiStripTkMaps::getTheZAxisTitle ( )
inline

Definition at line 60 of file SiStripTkMaps.h.

References m_zAxisTitle.

60 { return m_zAxisTitle; }
std::string m_zAxisTitle
Definition: SiStripTkMaps.h:75

◆ readVertices()

void SiStripTkMaps::readVertices ( double &  minx,
double &  maxx,
double &  miny,
double &  maxy 
)
private

Definition at line 224 of file SiStripTkMaps.cc.

References hcalRecHitTable_cff::detId, Exception, contentValuesFiles::fullPath, mps_fire::i, recoMuon::in, fastTrackerRecHitType::isPixel(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::ix(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::endcap::iy(), mps_splice::line, m_bins, m_detIdVector, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, AlCaHLTBitMon_QueryRunRegistry::string, submitPVValidationJobs::t, x, and y.

Referenced by bookMap().

224  {
225  std::ifstream in;
226 
227  // TPolyline vertices stored at https://github.com/cms-data/DQM-SiStripMonitorClient
228  in.open(edm::FileInPath("DQM/SiStripMonitorClient/data/Geometry/tracker_map_bare").fullPath().c_str());
229 
230  if (!in.good()) {
231  throw cms::Exception("FileError") << "SiStripTkMaps: problem opening vertices file!!" << std::endl;
232  return;
233  }
234 
235  while (in.good()) {
236  long detid = 0;
237  double x[5], y[5];
238 
240  std::getline(in, line);
241  typedef boost::tokenizer<boost::char_separator<char>> tokenizer;
242  boost::char_separator<char> sep{" "};
243  tokenizer tok{line, sep};
244 
245  int ix{0}, iy{0};
246  bool isPixel{false};
247  for (const auto& t : tok | boost::adaptors::indexed(0)) {
248  int i = t.index();
249  if (i == 0) {
250  detid = atoll((t.value()).c_str());
251 
252  // Drop Pixel Data
253  DetId detId(detid);
254  if (detId.subdetId() == PixelSubdetector::PixelBarrel || detId.subdetId() == PixelSubdetector::PixelEndcap) {
255  isPixel = true;
256  break;
257  }
258  } else {
259  if (i % 2 == 0) {
260  x[ix] = atof((t.value()).c_str());
261  if (x[ix] < minx) {
262  minx = x[ix];
263  }
264  if (x[ix] > maxx) {
265  maxx = x[ix];
266  }
267  ++ix;
268  } else {
269  y[iy] = atof((t.value()).c_str());
270  if (y[iy] < miny) {
271  miny = y[iy];
272  }
273  if (y[iy] > maxy) {
274  maxy = y[iy];
275  }
276  ++iy;
277  } // else
278  } // else
279  } // loop on entries
280 
281  if (isPixel) {
282  continue;
283  }
284 
285  m_bins[detid] = std::make_shared<TGraph>(ix, x, y);
286  m_bins[detid]->SetName(TString::Format("%ld", detid));
287  m_bins[detid]->SetTitle(TString::Format("Module ID=%ld", detid));
288  m_detIdVector.push_back(detid);
289  }
290 }
ALPAKA_FN_ACC ALPAKA_FN_INLINE uint32_t ix(uint32_t id)
Definition: DetId.h:17
std::vector< unsigned int > m_detIdVector
Definition: SiStripTkMaps.h:78
bool isPixel(HitType hitType)
ALPAKA_FN_ACC ALPAKA_FN_INLINE uint32_t iy(uint32_t id)
std::map< long, std::shared_ptr< TGraph > > m_bins
Definition: SiStripTkMaps.h:77

◆ setZAxisRange()

void SiStripTkMaps::setZAxisRange ( double  xmin,
double  xmax 
)
inline

Member Data Documentation

◆ m_axmax

double SiStripTkMaps::m_axmax
private

Definition at line 76 of file SiStripTkMaps.h.

◆ m_axmin

double SiStripTkMaps::m_axmin
private

Definition at line 76 of file SiStripTkMaps.h.

◆ m_bins

std::map<long, std::shared_ptr<TGraph> > SiStripTkMaps::m_bins
private

Definition at line 77 of file SiStripTkMaps.h.

Referenced by bookMap(), and readVertices().

◆ m_detIdVector

std::vector<unsigned int> SiStripTkMaps::m_detIdVector
private

Definition at line 78 of file SiStripTkMaps.h.

Referenced by getTheFilledIds(), and readVertices().

◆ m_mapTitle

std::string SiStripTkMaps::m_mapTitle = ""
private

Definition at line 74 of file SiStripTkMaps.h.

Referenced by bookMap(), dressMap(), and getTheMapTitle().

◆ m_option

Option_t* SiStripTkMaps::m_option
private

Definition at line 73 of file SiStripTkMaps.h.

Referenced by bookMap().

◆ m_trackerMap

TH2Poly* SiStripTkMaps::m_trackerMap {nullptr}
private

Definition at line 81 of file SiStripTkMaps.h.

Referenced by bookMap(), drawMap(), fill(), getTheMap(), and setZAxisRange().

◆ m_trackerTopo

TrackerTopology SiStripTkMaps::m_trackerTopo
private

Definition at line 80 of file SiStripTkMaps.h.

◆ m_values

std::vector<double> SiStripTkMaps::m_values
private

Definition at line 79 of file SiStripTkMaps.h.

Referenced by fill(), and getTheFilledValues().

◆ m_zAxisTitle

std::string SiStripTkMaps::m_zAxisTitle = ""
private

Definition at line 75 of file SiStripTkMaps.h.

Referenced by bookMap(), and getTheZAxisTitle().