CMS 3D CMS Logo
Functions
plotT0FromHistos Namespace Reference

Functions

def binNumber (station, sl)
 
def compareDiff (fileNames, sl, ymin=-15., ymax=15.)
 
def plot (fileName, sl, run, ymin=-5, ymax=5, option="HISTOP", draw=True)
 
def SLcompare (fileName, sls, run, ymin=-5, ymax=5, labels=[])
 

Function Documentation

def plotT0FromHistos.binNumber (   station,
  sl 
)

Definition at line 4 of file plotT0FromHistos.py.

Referenced by DTnoiseDBValidation.beginRun(), DTnoiseDBValidation.bookHisto(), plot(), and npstat::HistoND< Numeric, Axis >.setAxisLabel().

4 def binNumber(station,sl):
5  if sl is 3 and station is 4: return 11
6  if sl is 2 and station is 4: return None
7  start = (station - 1)*3
8  return start + sl
9 
plotT0FromHistos.binNumber
def binNumber(station, sl)
Definition: plotT0FromHistos.py:4
def plotT0FromHistos.compareDiff (   fileNames,
  sl,
  ymin = -15.,
  ymax = 15. 
)

Definition at line 110 of file plotT0FromHistos.py.

References drawHistoAllChambers.drawHisto(), HcalObjRepresent.Fill(), plot(), and HcalObjRepresent.Reset().

110 def compareDiff(fileNames,sl,ymin=-15.,ymax=15.):
111  option = "HISTOP"
112  colors = (2,4,9,12,38,44,46,55)
113  markers = (24,25,26,27,28,30,32,5)
114 
115  idx = 0
116  canvases = [None,None]
117  objects = None
118  histoRef = None
119  histos = []
120  histosDist = []
121  for fileName in fileNames:
122  objs = plot(fileName,sl,300,360,'',False)
123  histos.append( objs[1].Clone(objs[1].GetName() + "_diff") )
124  histos[-1].SetName( "%s_%d" % (histos[-1].GetName(),idx) )
125  if not idx:
126  histoRef = objs[1]
127  histos[-1].Reset()
128  else:
129  histos[-1].Add(histoRef,-1.)
130 
131  draw = False
132  if not idx: draw = True
133 
134  objs = drawHisto(histos[-1],
135  title="t_{Trig} difference (ns)",
136  ymin=ymin,ymax=ymax,option=option,draw=draw)
137 
138  if not idx:
139  canvases[0] = objs[0]
140  objects = objs[2]
141 
142  if idx:
143  canvases[0].cd()
144  histos[-1].SetLineColor(colors[ (idx - 1) % len(colors) ])
145  histos[-1].SetMarkerColor(colors[ (idx - 1) % len(colors) ])
146  histos[-1].SetMarkerStyle(markers[ (idx - 1) % len(markers) ])
147 
148  histos[-1].Draw(option + "SAME")
149 
150  histosDist.append( ROOT.TH1F(histos[-1].GetName() + "_dist","tTrig distribution",200,ymin,ymax) )
151  for ibin in range(1,histos[-1].GetNbinsX()+1):
152  histosDist[-1].Fill( histos[-1].GetBinContent(ibin) )
153 
154  histosDist[-1].SetLineColor(colors[ (idx - 1) % len(colors) ])
155  histosDist[-1].SetMarkerColor(colors[ (idx - 1) % len(colors) ])
156  histosDist[-1].SetMarkerStyle(markers[ (idx - 1) % len(markers) ])
157 
158  idx += 1
159 
160 
161  canvases[1] = ROOT.TCanvas("c_tTrigDist")
162  canvases[1].SetGridy()
163  canvases[1].SetFillColor(0)
164  canvases[1].cd()
165  option = "HISTO"
166  idx = 0
167  for histo in histosDist:
168  if not idx:
169  histo.GetXaxis().SetTitle("t_{Trig} difference (ns)")
170  histo.GetYaxis().SetTitle("Number of chambers")
171  histo.Draw(option)
172  else:
173  histo.Draw(option + "SAME")
174  idx += 1
175 
176  return (canvases,histos,histosDist,objects)
177 
plotT0FromHistos.plot
def plot(fileName, sl, run, ymin=-5, ymax=5, option="HISTOP", draw=True)
Definition: plotT0FromHistos.py:10
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
plotT0FromHistos.compareDiff
def compareDiff(fileNames, sl, ymin=-15., ymax=15.)
Definition: plotT0FromHistos.py:110
HcalObjRepresent::Reset
void Reset(std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:620
drawHistoAllChambers.drawHisto
def drawHisto(histo, title, ymin, ymax, option="HISTOP", draw=True)
Definition: drawHistoAllChambers.py:3
def plotT0FromHistos.plot (   fileName,
  sl,
  run,
  ymin = -5,
  ymax = 5,
  option = "HISTOP",
  draw = True 
)

Definition at line 10 of file plotT0FromHistos.py.

References binNumber(), drawHistoAllChambers.drawHisto(), and str.

Referenced by compareDiff(), and SLcompare().

10 def plot(fileName,sl,run,ymin=-5,ymax=5,option="HISTOP",draw=True):
11  path = "DQMData/Run "+str(run)+"/DT/Run summary/DtCalib/TTrigDBValidation/"
12  slType = sl
13  slStr = "SL%d" % slType
14  verbose = True #False
15 
16  ROOT.TH1.AddDirectory(False)
17 
18  file = ROOT.TFile(fileName,'read')
19 
20  wheels = (-2,-1,0,1,2)
21  stations = (1,2,3,4)
22 
23  #histosWheel = {}
24  #for wh in wheels:
25  # histoName = path+'Wheel%d_%s_TTrig' % (wh,slStr)
26  # print "Accessing",histoName
27  # histosWheel[wh] = file.Get(histoName)
28 
29  # (Wh-2 MB1 Sec1 ... Wh-2 MB1 Sec12 ... Wh-1 MB1 Sec1 ... Wh-1 MB1 Sec12 ...)
30  # (Wh-2 MB2 Sec1 ... Wh-2 MB2 Sec12 ... Wh-1 MB2 Sec1 ... Wh-1 MB1 Sec12 ...) ...
31  nBins = 250
32  #if slType == 2: nBins = 180
33  histo = ROOT.TH1F("h_TTrigAll","TTrig Run "+str(run),nBins,0,nBins)
34  for st in stations:
35  nSectors = 12
36  if st == 4: nSectors = 14
37  if st == 4 and sl == 2: continue
38  if verbose: print "Station",st
39  for wh in wheels:
40  if verbose: print "Wheel",wh
41  for sec in range(1,nSectors+1):
42  if verbose: print "Sector",sec
43  binHisto = binNumber(st,sl)
44  if verbose: print "Bin from histos:",binHisto
45  histoName = path+'Wheel%d/TTrigDifference_W%d_Sec%d' % (wh,wh,sec)
46  value = file.Get(histoName).GetBinContent(binHisto)
47 
48  binHistoNew = (st - 1)*60 + (wh + 2)*nSectors + sec
49  if verbose: print "Bin final",binHistoNew,value
50  histo.SetBinContent(binHistoNew,value)
51 
52  if sec == 1:
53  label = "Wheel %d" % wh
54  if wh == -2: label += " MB%d" % st
55  histo.GetXaxis().SetBinLabel(binHistoNew,label)
56 
57  objects = drawHisto(histo,
58  title="#Deltat_{Trig} (ns)",
59  ymin=ymin,ymax=ymax,option=option,draw=draw)
60 
61  return objects
62 
plotT0FromHistos.plot
def plot(fileName, sl, run, ymin=-5, ymax=5, option="HISTOP", draw=True)
Definition: plotT0FromHistos.py:10
plotT0FromHistos.binNumber
def binNumber(station, sl)
Definition: plotT0FromHistos.py:4
str
#define str(s)
Definition: TestProcessor.cc:47
drawHistoAllChambers.drawHisto
def drawHisto(histo, title, ymin, ymax, option="HISTOP", draw=True)
Definition: drawHistoAllChambers.py:3
def plotT0FromHistos.SLcompare (   fileName,
  sls,
  run,
  ymin = -5,
  ymax = 5,
  labels = [] 
)

Definition at line 63 of file plotT0FromHistos.py.

References plot().

63 def SLcompare(fileName,sls,run,ymin=-5,ymax=5,labels=[]):
64  option = "HISTOP"
65  colors = (2,4,12,44,55,38,27,46)
66  markers = (24,25,26,27,28,30,32,5)
67 
68  idx = 0
69  canvas = None
70  objects = None
71  histos = []
72  for sl in sls:
73  draw = False
74  if not idx: draw = True
75 
76  objs = plot(fileName,sl,run,ymin,ymax,option,draw)
77  histos.append(objs[1])
78  histos[-1].SetName( "SL %d" % sl)
79  if not idx:
80  canvas = objs[0]
81  objects = objs[2]
82 
83  canvas.cd()
84  if idx:
85  histos[-1].SetLineColor(colors[ (idx - 1) % len(colors) ])
86  histos[-1].SetMarkerColor(colors[ (idx - 1) % len(colors) ])
87  histos[-1].SetMarkerStyle(markers[ (idx - 1) % len(markers) ])
88 
89  histos[-1].Draw(option + "SAME")
90 
91  idx += 1
92 
93  legend = ROOT.TLegend(0.4,0.7,0.95,0.8)
94  for idx in range( len(histos) ):
95  histo = histos[idx]
96  label = histo.GetName()
97  if len(labels): label = labels[idx]
98  legend.AddEntry(histo,label,"LP")
99 
100  idx += 1
101 
102  canvas.cd()
103  legend.SetFillColor( canvas.GetFillColor() )
104  legend.Draw("SAME")
105 
106  objects.append(legend)
107 
108  return (canvas,histos,objects)
109 
plotT0FromHistos.plot
def plot(fileName, sl, run, ymin=-5, ymax=5, option="HISTOP", draw=True)
Definition: plotT0FromHistos.py:10
plotT0FromHistos.SLcompare
def SLcompare(fileName, sls, run, ymin=-5, ymax=5, labels=[])
Definition: plotT0FromHistos.py:63
Generated for CMSSW Reference Manual by doxygen 1.8.11