Reference Manual

Functions
def	acustompalette

def	drawHalfEtaValues

def	TwikiPrintout

Variables
dictionary	_LABELS2COMPS

tuple	COMPOUNDS = OrderedDict()

tuple	dEdx = OrderedDict()

tuple	DETECTORS = OrderedDict()

tuple	hist_label_to_num = OrderedDict()

tuple	MatXo = OrderedDict()

dictionary	plots = {}

tuple	sDETS = OrderedDict()

Function Documentation

def plot_hgcal_utils.acustompalette ( )

Definition at line 192 of file plot_hgcal_utils.py.

References mps_check.array, and sistrip::SpyUtilities.range().

 
 def acustompalette(): 
     NRGBs = 7
     NCont = 100
     ncolors = array('i', [])
     gStyle.SetNumberContours(NCont);
     stops   = [ 0.00, 0.10, 0.25, 0.45, 0.60, 0.75, 1.00 ]
     red     = [ 1.00, 0.00, 0.00, 0.00, 0.97, 0.97, 0.10 ]
     green   = [ 1.00, 0.97, 0.30, 0.40, 0.97, 0.00, 0.00 ]
     blue    = [ 1.00, 0.97, 0.97, 0.00, 0.00, 0.00, 0.00 ]
     stopsArray = array('d', stops)
     redArray = array('d', red)
     greenArray = array('d', green)
     blueArray = array('d', blue)
     first_color_number = TColor.CreateGradientColorTable(NRGBs, stopsArray, redArray, greenArray, blueArray, NCont);
     gStyle.SetNumberContours(NCont)
 
 
     palsize = NCont
     palette = []
     for i in range(palsize):
         palette.append(first_color_number+i)
         palarray = array('i',palette)
 
     gStyle.SetPalette(palsize,palarray)
 
 
#In MeV/mm

mps_check.array

tuple array

Definition: mps_check.py:216

sistrip::SpyUtilities::range

const uint16_t range(const Frame &aFrame)

Definition: SiStripSpyUtilities.cc:76

plot_hgcal_utils.acustompalette

def acustompalette

Definition: plot_hgcal_utils.py:192

def plot_hgcal_utils.drawHalfEtaValues ( )

Function to draw the eta.
Function to draw the eta references on top of an already existing
TCanvas. The lines and labels drawn are collected inside a list and
the list is returned to the user to extend the live of the objects
contained, otherwise no lines and labels will be drawn, since they
will be garbage-collected as soon as this function returns.

Definition at line 270 of file plot_hgcal_utils.py.

References funct.cos(), funct.exp(), sistrip::SpyUtilities.range(), funct.sin(), and funct.tan().

 
 def drawHalfEtaValues():
     """Function to draw the eta.
     Function to draw the eta references on top of an already existing
     TCanvas. The lines and labels drawn are collected inside a list and
     the list is returned to the user to extend the live of the objects
     contained, otherwise no lines and labels will be drawn, since they
     will be garbage-collected as soon as this function returns.
     """
 
     # Add eta labels
     keep_alive = []
     etas = [ 0.2*i for i in range(0,18) ]
 
     etax = 2850.#6850.
     etay = 1240.#5200.
     lineL = 110#8600.
     offT = 10.
 
     for ieta in etas:
         th = 2*atan(exp(-ieta))
         talign = 21
 
         #IP
         lineh = TLine(-20.,0.,20.,0.)
         lineh.Draw()
         linev = TLine(0.,-10.,0.,10.)
         linev.Draw()
         keep_alive.append(lineh)
         keep_alive.append(linev)
 
         x1 = 0
         y1 = 0
         if ieta>-1.6 and ieta<1.6:
             x1 = etay/tan(th)
             y1 = etay
         elif ieta <=-1.6:
             x1 = -etax
             y1 = -etax*tan(th)
             talign = 11
         elif ieta>=1.6:
             x1 = etax
             y1 = etax*tan(th)
             talign = 31
         x2 = x1+lineL*cos(th)
         y2 = y1+lineL*sin(th)
         xt = x2
         yt = y2+offT
 
         line1 = TLine(x1,y1,x2,y2)
         line1.Draw()
         keep_alive.append(line1)
 
         text = "%3.1f" % ieta
         t1 = TLatex(xt, yt, '%s' % ('#eta = 0' if ieta == 0 else text))
         t1.SetTextSize(0.03)
         t1.SetTextAlign(talign)
         t1.Draw()
         keep_alive.append(t1)
     return keep_alive

def plot_hgcal_utils.TwikiPrintout	(	plotname,
		label,
		zoom
	)

The plots in the twiki are already too much and to avoid mistakes 
   we will try to automatize the procedure

Definition at line 131 of file plot_hgcal_utils.py.

References print().

 
 def TwikiPrintout(plotname, label, zoom): 
     """The plots in the twiki are already too much and to avoid mistakes 
        we will try to automatize the procedure
     """
 
     #Twiki will strip out spaces
     label = label.replace(" ", "_")
 
     zoomstring = ""
 
     if zoom == "all":
         zoomstring = ""
         zoomtitle = "in all HGCal"
         zoomdir = "%s/" % label
     elif zoom == "zplus":
         zoomstring = "_ZplusZoom"
         zoomtitle = "in Z+ endcap of HGCal"
         zoomdir = "%s/ZPlusZoom/" % label
     elif zoom == "zminus":
         zoomstring = "_ZminusZoom"
         zoomtitle = "in Z- endcap of HGCal"
         zoomdir = "%s/ZMinusZoom/" % label
     else :
         print("WRONG OPTION")
 
 
     #Here for the hide button
     if plotname == "x_vs_z_vs_Rsum":
         print("%%TWISTY{ mode=\"div\" showlink=\"Click to see the %s plots %s \" hidelink=\"Hide %s %s\" showimgright=\"%%ICONURLPATH{toggleopen-small}%%\" hideimgright=\"%%ICONURLPATH{toggleclose-small}%%\"}%%" % (label,zoomtitle, label, zoomtitle))
 
     if "Rsum" in plotname and "x_vs" in plotname and not "cos" in plotname: 
         print("| <img alt=\"HGCal_%s%s%s.png\" height=\"300\" width=\"550\" src=\"http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.png\" /> | The plot on the left shows the 2D profile histogram for *%s* %s that displays the mean value of the material budget in units of radiation length in each R-z cell. R-z cell is 1 cm x 1 mm. The plot depicts the accumulated material budget as seen by the track, as the track travels throughout the detector.[[http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.pdf][Click to enlarge plot]] |" % (plotname,label,zoomstring,zoomdir,plotname,label,zoomstring, label, zoomtitle,zoomdir,plotname,label,zoomstring))
 
     if "Rsum" in plotname and "l_vs" in plotname and not "cos" in plotname: 
         print("| <img alt=\"HGCal_%s%s%s.png\" height=\"300\" width=\"550\" src=\"http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.png\" /> | The plot on the left shows the 2D profile histogram for *%s* %s that displays the mean value of the material budget in units of interaction length in each R-z cell. R-z cell is 1 cm x 1 mm. The plot depicts the accumulated material budget as seen by the track, as the track travels throughout the detector.[[http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.pdf][Click to enlarge plot]] |" % (plotname,label,zoomstring,zoomdir,plotname,label,zoomstring, label, zoomtitle,zoomdir,plotname,label,zoomstring))
 
     if "Rsumcos" in plotname and "x_vs" in plotname: 
         print("| <img alt=\"HGCal_%s%s%s.png\" height=\"300\" width=\"550\" src=\"http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.png\" /> | The plot on the left shows the 2D profile histogram for *%s* %s that displays the mean value of the material budget in units of radiation length in each R-z cell. R-z cell is 1 cm x 1 mm. The plot depicts the orthogonal accumulated material budget, that is cos(theta) what the track sees.[[http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.pdf][Click to enlarge plot]] |" % (plotname,label,zoomstring,zoomdir,plotname,label,zoomstring, label, zoomtitle,zoomdir,plotname,label,zoomstring))
 
     if "Rsumcos" in plotname and "l_vs" in plotname: 
         print("| <img alt=\"HGCal_%s%s%s.png\" height=\"300\" width=\"550\" src=\"http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.png\" /> | The plot on the left shows the 2D profile histogram for *%s* %s that displays the mean value of the material budget in units of interaction length in each R-z cell. R-z cell is 1 cm x 1 mm. The plot depicts the orthogonal accumulated material budget, that is cos(theta) what the track sees.[[http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.pdf][Click to enlarge plot]] |" % (plotname,label,zoomstring,zoomdir,plotname,label,zoomstring, label, zoomtitle,zoomdir,plotname,label,zoomstring))
 
     if "Rloc" in plotname and "x_vs" in plotname and not "cos" in plotname: 
         print("| <img alt=\"HGCal_%s%s%s.png\" height=\"300\" width=\"550\" src=\"http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.png\" /> | The plot on the left shows the 2D profile histogram for *%s* %s that displays the local mean value of the material budget in units of radiation length in each R-z cell. R-z cell is 1 cm x 1 mm. The plot depicts the local material budget as seen by the track, as the track travels throughout the detector.[[http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.pdf][Click to enlarge plot]] |" % (plotname,label,zoomstring,zoomdir,plotname,label,zoomstring, label, zoomtitle,zoomdir,plotname,label,zoomstring))
 
     if "Rloc" in plotname and "l_vs" in plotname and not "cos" in plotname: 
         print("| <img alt=\"HGCal_%s%s%s.png\" height=\"300\" width=\"550\" src=\"http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.png\" /> | The plot on the left shows the 2D profile histogram for *%s* %s that displays the local mean value of the material budget in units of interaction length in each R-z cell. R-z cell is 1 cm x 1 mm. The plot depicts the local material budget as seen by the track, as the track travels throughout the detector.[[http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.pdf][Click to enlarge plot]] |" % (plotname,label,zoomstring,zoomdir,plotname,label,zoomstring, label, zoomtitle,zoomdir,plotname,label,zoomstring))
 
     #Here again for the closing of the hide button
     if plotname == "l_vs_z_vs_Rloc":
         print("%ENDTWISTY%")
 
     """ 
     I won't put the local cos plots for now, only the sum cos above
     if "Rloccos" in plotname and "x_vs" in plotname: 
          print "| <img alt=\"HGCal_%s%s%s.png\" height=\"300\" width=\"550\" src=\"http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.png\" /> | The plot on the left shows the 2D profile histogram for *%s* %s that displays the local mean value of the material budget in units of radiation length in each R-z cell. R-z cell is 1 cm x 1 mm. The plot depicts the orthogonal accumulated material budget, that is cos(theta) what the track sees.[[http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.pdf][Click to enlarge plot]] |" % (plotname,label,zoomstring,zoomdir,plotname,label,zoomstring, label, zoomtitle,zoomdir,plotname,label,zoomstring)
 
     if "Rloccos" in plotname and "l_vs" in plotname: 
          print "| <img alt=\"HGCal_%s%s%s.png\" height=\"300\" width=\"550\" src=\"http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.png\" /> | The plot on the left shows the 2D profile histogram for *%s* %s that displays the local mean value of the material budget in units of interaction length in each R-z cell. R-z cell is 1 cm x 1 mm. The plot depicts the orthogonal accumulated material budget, that is cos(theta) what the track sees.[[http://apsallid.web.cern.ch/apsallid/HGCalMaterial/%sHGCal_%s%s%s.pdf][Click to enlarge plot]] |" % (plotname,label,zoomstring,zoomdir,plotname,label,zoomstring, label, zoomtitle,zoomdir,plotname,label,zoomstring)
     """

void print(TMatrixD &m, const char *label=nullptr, bool mathematicaFormat=false)

Definition: Utilities.cc:47

plot_hgcal_utils.TwikiPrintout

def TwikiPrintout

Definition: plot_hgcal_utils.py:131

Functions

Variables

Function Documentation

Variable Documentation