PlotCalibrationGains Namespace Reference

Classes
class	EmPartitionPlots
class	GainReader
class	HadPartitionPlots
class	L1CaloGeometryConvertor
class	L1CaloMap

Functions
	setAxisSizes (histos, tsize=0.05, lsize=0.05, zsize=0.05)
	setAxisTitles (histos, xtit, ytit, ztit="foo", xoff=0, yoff=0, zoff=0)
	PlotCalibrationGains (input_file_name="", reference_file_name="", isInputXml=False, isInputSqlite=False, isRefXml=False, isRefSqlite=False, isRefOracle=False, isTileOnly=False)

Variables
	parser = OptionParser()
	action
	type
	dest
	help
	default
	options
	args

Function Documentation

PlotCalibrationGains()

PlotCalibrationGains.PlotCalibrationGains	(		input_file_name = "",
			reference_file_name = "",
			isInputXml = False,
			isInputSqlite = False,
			isRefXml = False,
			isRefSqlite = False,
			isRefOracle = False,
			isTileOnly = False )

Definition at line 608 of file PlotCalibrationGains.py.

                         isTileOnly=False):
 
    ROOT.gROOT.SetBatch( True )    
    ROOT.gStyle.SetPalette(1)
    ROOT.gStyle.SetOptStat(111111)
    ROOT.gStyle.SetCanvasColor(10)
 
    ROOT.gStyle.SetPadTopMargin(0.12)
    ROOT.gStyle.SetPadBottomMargin(0.12)
    ROOT.gStyle.SetPadRightMargin(0.12)
    ROOT.gStyle.SetPadLeftMargin(0.12)
    
    canvas = ROOT.TCanvas('canvas','Gains',0,0,600,800)
    canvas.SetBatch(True)
  
    h_gains_em  = L1CaloMap("EM receiver gains","#eta bin","#phi bin")
    h_gains_had = L1CaloMap("HAD receiver gains","#eta bin","#phi bin")
 
    h_gains_em_fselect  = L1CaloMap("EM gains that failed selection","#eta bin","#phi bin")
    h_gains_had_fselect = L1CaloMap("HAD gains that failed selection","#eta bin","#phi bin")
 
    h_chi2_em  = L1CaloMap("EM fit Chi2","#eta bin","#phi bin")
    h_chi2_had = L1CaloMap("HAD fit Chi2","#eta bin","#phi bin")
 
    h_offset_em  = L1CaloMap("EM fit offsets","#eta bin","#phi bin")
    h_offset_had = L1CaloMap("HAD fit offsets","#eta bin","#phi bin")
 
    h_unfitted_em  = L1CaloMap("EM failed fits","#eta bin","#phi bin")
    h_unfitted_had = L1CaloMap("HAD failed fits","#eta bin","#phi bin")
 
    h_drifted_em  = L1CaloMap("EM TTs that drifted by more than 10%","#eta bin","#phi bin")
    h_drifted_had = L1CaloMap("HAD TTs that drifted by more than 10%","#eta bin","#phi bin")
 
    h_gains_em_reference  = L1CaloMap("EM (gain-reference)","#eta bin","#phi bin")
    h_gains_had_reference = L1CaloMap("HAD (gain-reference)","#eta bin","#phi bin")
 
    h_gains_em_reference_rel  = L1CaloMap("EM (gain-reference)/reference","#eta bin","#phi bin")
    h_gains_had_reference_rel = L1CaloMap("HAD (gain-reference)/reference","#eta bin","#phi bin")
    em_partition_gains = EmPartitionPlots(" EM gains",40,0.6,1.4,"gain","Entries")
    had_partition_gains = HadPartitionPlots(" HAD gains",40,0.5,2.5,"gain","Entries")
 
    em_partition_gains_ref = EmPartitionPlots(" EM (gains-reference)",40,-0.2,0.2,"gain-reference","Entries")
    had_partition_gains_ref = HadPartitionPlots(" HAD (gains-reference)",40,-1.,1.,"gain-reference","Entries")
 
    em_partition_gains_ref_rel = EmPartitionPlots(" EM (gains-reference)/reference",40,-0.2,0.2,"(gain-reference)/reference","Entries")
    had_partition_gains_ref_rel = HadPartitionPlots(" HAD (gains-reference)/reference",40,-1.,1.,"(gain-reference)/reference","Entries")
 
    threshold_change = 0.1
 
    geometry_convertor = L1CaloGeometryConvertor()
    receiver_gains    = GainReader()
 
    bad_gain_file = open('bad_gains.txt','w')
    drifted_towers_file = open('drifted_towers.txt','w')
  
    if isInputXml is True:
        print ("Taking input from xml file: ", input_file_name)
        receiver_gains.LoadGainsXml(input_file_name)
    elif isInputSqlite is True:
        print ("Taking input from Sqlite file: ", input_file_name)
        receiver_gains.LoadGainsSqlite(input_file_name)
    else:
        print ("No option for input file selected, assuming sqlite file energyscanresults.sqlite")
        receiver_gains.LoadGainsSqlite("energyscanresults.sqlite")
        
    if isRefXml is True:
        print ("Taking reference from Xml file: ",reference_file_name)
        receiver_gains.LoadReferenceXml(reference_file_name)
    elif isRefSqlite is True:
        print ("Taking reference from Sqlite file: ",reference_file_name)
        receiver_gains.LoadReferenceSqlite(reference_file_name)
    elif isRefOracle is True:
        print ("Taking reference from Oracle")
        geometry_convertor.LoadReceiverPPMMap()
        receiver_gains.LoadReferenceOracle(geometry_convertor)
    else:
        print (" No option for reference file, assuming Oracle")
        geometry_convertor.LoadReceiverPPMMap()
        receiver_gains.LoadReferenceOracle(geometry_convertor)
                
    for i_eta in range(-49,45):
        for i_phi in range(0,64):
            coolEm  = geometry_convertor.getCoolEm(i_eta,i_phi)
            coolHad = geometry_convertor.getCoolHad(i_eta,i_phi)
      
            if not coolEm == '':                      # there is a channel for this eta-phi      
                gain   = receiver_gains.getGain(coolEm)
                chi2   = receiver_gains.getChi2(coolEm)
                offset = receiver_gains.getOffset(coolEm)
                reference_gain = receiver_gains.getReferenceGain(coolEm)
                passes_selection = receiver_gains.passesSelection(coolEm)
 
                if (not gain == '') and (not reference_gain == ''):  # both  gains should be available
 
                    if gain == -1. :
                        h_unfitted_em.Fill(i_eta,i_phi)  
                        bad_gain_file.write('%i %i %s EM gain= %.3f \n' % (i_eta,i_phi,coolEm,float(gain))) 
                        h_gains_em_reference.Fill(i_eta,i_phi,-100.)
                        h_gains_em_reference_rel.Fill(i_eta,i_phi,-100.)
                    else:
                        h_gains_em.Fill(i_eta,i_phi,gain)
                        h_chi2_em.Fill(i_eta,i_phi,chi2)
                        h_offset_em.Fill(i_eta,i_phi,offset)   
                        em_partition_gains.Fill(i_eta,gain)         
                        em_partition_gains_ref.Fill(i_eta,gain-reference_gain)  
                        h_gains_em_reference.Fill(i_eta,i_phi,gain-reference_gain)
 
                        if passes_selection is False:
                            h_gains_em_fselect.Fill(i_eta,i_phi)
                            bad_gain_file.write('%i %i %s  EM gain= %.3f  chi2= %.3f offset= %.3f \n' %
                                                (i_eta,i_phi,coolEm,float(gain),float(chi2),float(offset))) 
                            #h_gains_em_reference.Fill(i_eta,i_phi,-100.)
                            #h_gains_em_reference_rel.Fill(i_eta,i_phi,-100.)
 
                        if reference_gain > 0:
                            em_partition_gains_ref_rel.Fill(i_eta,(gain-reference_gain)/reference_gain)  
                            h_gains_em_reference_rel.Fill(i_eta,i_phi,(gain-reference_gain)/reference_gain)
                            if fabs((gain-reference_gain)/reference_gain) > threshold_change:
                                h_drifted_em.Fill(i_eta,i_phi)
                                drifted_towers_file.write('%i %i %s  EM gain= %.3f   refGain= %.3f   (%.3f %%) \n' %
                                                          (i_eta,i_phi,coolEm,float(gain),float(reference_gain),
                                                           (gain-reference_gain)*100/reference_gain)) 
 
 
            if not coolHad == '':                    # there is a channel for this eta-phi
                gain = receiver_gains.getGain(coolHad)
                chi2   = receiver_gains.getChi2(coolHad)
                offset = receiver_gains.getOffset(coolHad)
                reference_gain = receiver_gains.getReferenceGain(coolHad)
                passes_selection = receiver_gains.passesSelection(coolHad)
 
                if (not gain == '') and (not reference_gain == ''):      # both gains should be available
 
                    if gain == -1. :
                        h_unfitted_had.Fill(i_eta,i_phi)  
                        bad_gain_file.write('%i %i %s HAD gain= %.3f \n' % (i_eta,i_phi,coolHad,float(gain)))
                        h_gains_had_reference.Fill(i_eta,i_phi,-100.)
                        h_gains_had_reference_rel.Fill(i_eta,i_phi,-100.)
                    else:
                        h_gains_had.Fill(i_eta,i_phi,gain)
                        h_chi2_had.Fill(i_eta,i_phi,chi2)
                        h_offset_had.Fill(i_eta,i_phi,offset)
                        had_partition_gains.Fill(i_eta,gain)
                        had_partition_gains_ref.Fill(i_eta,gain-reference_gain)
                        h_gains_had_reference.Fill(i_eta,i_phi,gain-reference_gain)
 
                        if passes_selection is False:
                            h_gains_had_fselect.Fill(i_eta,i_phi)
                            bad_gain_file.write( '%i %i %s  HAD gain= %.3f  chi2= %.3f offset= %.3f \n' %
                                                 (i_eta,i_phi,coolHad,float(gain),float(chi2),float(offset))) 
                            #h_gains_had_reference.Fill(i_eta,i_phi,-100.)
                            #h_gains_had_reference_rel.Fill(i_eta,i_phi,-100.)
 
                        if reference_gain > 0:
                            had_partition_gains_ref_rel.Fill(i_eta,(gain-reference_gain)/reference_gain)
                            h_gains_had_reference_rel.Fill(i_eta,i_phi,(gain-reference_gain)/reference_gain)
                            if fabs((gain-reference_gain)/reference_gain) > threshold_change:
                                h_drifted_had.Fill(i_eta,i_phi)
                                drifted_towers_file.write('%i %i %s  HAD gain= %.3f   refGain= %.3f   (%.3f %%) \n' %
                                                          (i_eta,i_phi,coolHad,float(gain),float(reference_gain),
                                                           (gain-reference_gain)*100/reference_gain))
 
    # start the plotting
    pdfFileName = "Gains.pdf"
    canvas.Print( pdfFileName + "[" )
 
    # new page: EM gains and drifts
    if not isTileOnly:
        canvas.Clear()
        canvas.Divide(1,2)
        canvas.cd(1)
        ROOT.gPad.SetLogy(0)
        h_gains_em.SetMinimum(0.6)
        h_gains_em.SetMaximum(1.4)
        h_gains_em.Draw()
        canvas.cd(2)
        h_drifted_em.Draw()
        canvas.Print( pdfFileName )
 
    # new page: EM gains that failed selection and failed fits
    if not isTileOnly:
        canvas.Clear()
        canvas.Divide(1,2)
        canvas.cd(1)
        ROOT.gPad.SetLogy(0)
        h_gains_em_fselect.Draw()
        canvas.cd(2)
        h_unfitted_em.Draw()
        canvas.Print( pdfFileName )
 
    # new page: HAD gains and drifts
    canvas.Clear()
    canvas.Divide(1,2)
    canvas.cd(1)
    ROOT.gPad.SetLogy(0)
    h_gains_had.SetMinimum(0.6)
    h_gains_had.SetMaximum(1.4)
    h_gains_had.Draw()
    canvas.cd(2)
    h_drifted_had.Draw()
    canvas.Print( pdfFileName )
 
    # new page: HAD gains that failed selection and failed fits
    canvas.Clear()
    canvas.Divide(1,2)
    canvas.cd(1)
    h_gains_had_fselect.Draw()
    canvas.cd(2)
    h_unfitted_had.Draw()
    canvas.Print( pdfFileName )
    
    # new page: EM gains relative references
    if not isTileOnly:
        canvas.Clear()
        canvas.Divide(1,2)
        canvas.cd(1)
        ROOT.gPad.SetLogy(0)
        h_gains_em_reference_rel.SetMinimum(-0.5)
        h_gains_em_reference_rel.SetMaximum(0.5)
        h_gains_em_reference_rel.Draw()
        canvas.cd(2)
        h_gains_em_reference_rel.SetMinimum(-0.1)
        h_gains_em_reference_rel.SetMaximum(0.1)
        h_gains_em_reference_rel.Draw()
        canvas.Print( pdfFileName )
 
    # new page: HAD gains relative references
    canvas.Clear()
    canvas.Divide(1,2)
    canvas.cd(1)
    ROOT.gPad.SetLogy(0)
    h_gains_had_reference_rel.SetMinimum(-0.5)
    h_gains_had_reference_rel.SetMaximum(0.5)
    h_gains_had_reference_rel.Draw()
    canvas.cd(2)    
    h_gains_had_reference_rel.SetMinimum(-0.1)
    h_gains_had_reference_rel.SetMaximum(0.1)
    h_gains_had_reference_rel.Draw()
    canvas.Print( pdfFileName )
 
    # new page: EM gains absolute references
    if not isTileOnly:
        canvas.Clear()
        canvas.Divide(1,2)
        canvas.cd(1)
        ROOT.gPad.SetLogy(0)
        h_gains_em_reference.SetMinimum(-0.5)
        h_gains_em_reference.SetMaximum(0.5)
        h_gains_em_reference.Draw()
        canvas.cd(2)
        ROOT.gPad.SetLogy(0)
        h_gains_em_reference.SetMinimum(-0.1)
        h_gains_em_reference.SetMaximum(0.1)
        h_gains_em_reference.Draw()
        canvas.Print( pdfFileName )
 
    # new page: HAD gains absolute references
    canvas.Clear()
    canvas.Divide(1,2)
    canvas.cd(1)
    ROOT.gPad.SetLogy(0)
    h_gains_had_reference.SetMinimum(-0.5)
    h_gains_had_reference.SetMaximum(0.5)
    h_gains_had_reference.Draw()
    canvas.cd(2)
    ROOT.gPad.SetLogy(0)
    h_gains_had_reference.SetMinimum(-0.1)
    h_gains_had_reference.SetMaximum(0.1)
    h_gains_had_reference.Draw()
    canvas.Print( pdfFileName )
 
    # new page: EM chi2 and offsets
    if not isTileOnly:
        canvas.Clear()
        canvas.Divide(1,2)
        canvas.cd(1)
        ROOT.gPad.SetLogy(0)
        h_chi2_em.SetMinimum(0.1)
        h_chi2_em.SetMaximum(100)
        h_chi2_em.Draw()
        canvas.cd(2)
        ROOT.gPad.SetLogy(0)
        h_offset_em.SetMinimum(-1.)
        h_offset_em.SetMaximum(1.)
        h_offset_em.Draw()
        canvas.Print( pdfFileName )
 
    # new page: HAD chi2 and offsets
    canvas.Clear()
    canvas.Divide(1,2)
    canvas.cd(1)
    ROOT.gPad.SetLogy(0)
    h_chi2_had.SetMinimum(0.1)
    h_chi2_had.SetMaximum(100)
    h_chi2_had.Draw()
    canvas.cd(2)
    ROOT.gPad.SetLogy(0)
    h_offset_had.SetMinimum(-1.)
    h_offset_had.SetMaximum(1.)
    h_offset_had.Draw()    
    canvas.Print( pdfFileName )
 
    # new page: EM 1D gains per layer
    if not isTileOnly:
        canvas.Clear()
        canvas.Divide(2,3)
        for i_p in range(0,em_partition_gains.nPartitions):
            canvas.cd(i_p+1)
            ROOT.gPad.SetLogy() if em_partition_gains.his_partitions[i_p].GetEntries()>0 else ROOT.gPad.SetLogy(0)
            em_partition_gains.his_partitions[i_p].Draw()
        canvas.Print( pdfFileName )
 
    # new page: HAD 1D gains per layer
    canvas.Clear()
    canvas.Divide(2,3)
    for i_p in range(0,had_partition_gains.nPartitions):
        if (isTileOnly and i_p>2): break
        canvas.cd(i_p+1)
        ROOT.gPad.SetLogy() if  had_partition_gains.his_partitions[i_p].GetEntries()>0 else ROOT.gPad.SetLogy(0)
        had_partition_gains.his_partitions[i_p].Draw()
    canvas.Print( pdfFileName )
 
    # new page: EM 1D gains absolute references per layer
    if not isTileOnly:
        canvas.Clear()
        canvas.Divide(2,3)
        for i_p in range(0,em_partition_gains_ref.nPartitions):
            canvas.cd(i_p+1)
            ROOT.gPad.SetLogy() if em_partition_gains_ref.his_partitions[i_p].GetEntries()>0 else ROOT.gPad.SetLogy(0)
            em_partition_gains_ref.his_partitions[i_p].Draw()  
        canvas.Print( pdfFileName )
    
    # new page: HAD 1D gains absolute references per layer
    canvas.Clear()
    canvas.Divide(2,3)
    for i_p in range(0,had_partition_gains_ref.nPartitions):
        if (isTileOnly and i_p>2): break
        canvas.cd(i_p+1)
        ROOT.gPad.SetLogy() if had_partition_gains_ref.his_partitions[i_p].GetEntries()>0 else ROOT.gPad.SetLogy(0)
        had_partition_gains_ref.his_partitions[i_p].Draw()  
    canvas.Print( pdfFileName )
 
    # new page: EM 1D gains relative references per layer
    if not isTileOnly:
        canvas.Clear()
        canvas.Divide(2,3)
        for i_p in range(0,em_partition_gains_ref_rel.nPartitions):
            canvas.cd(i_p+1)
            ROOT.gPad.SetLogy() if em_partition_gains_ref_rel.his_partitions[i_p].GetEntries()>0 else ROOT.gPad.SetLogy(0)
            em_partition_gains_ref_rel.his_partitions[i_p].Draw()  
        canvas.Print( pdfFileName )
 
    # new page: HAD 1D gains relative references per layer
    canvas.Clear()
    canvas.Divide(2,3)
    for i_p in range(0,had_partition_gains_ref_rel.nPartitions):
        if (isTileOnly and i_p>2): break
        canvas.cd(i_p+1)
        ROOT.gPad.SetLogy() if had_partition_gains_ref_rel.his_partitions[i_p].GetEntries()>0 else ROOT.gPad.SetLogy(0)
        had_partition_gains_ref_rel.his_partitions[i_p].Draw()  
    canvas.Print( pdfFileName )
 
    # Closing files
    canvas.Print( pdfFileName + "]" )
    bad_gain_file.close()
    drifted_towers_file.close()  
    print ("finished!")
 
    

setAxisSizes()

PlotCalibrationGains.setAxisSizes	(		histos,
			tsize = 0.05,
			lsize = 0.05,
			zsize = 0.05 )

Definition at line 14 of file PlotCalibrationGains.py.

def setAxisSizes(histos, tsize=0.05, lsize=0.05, zsize=0.05):
    for h in histos:
        h.GetXaxis().SetTitleSize(tsize)
        h.GetXaxis().SetLabelSize(lsize)
        h.GetYaxis().SetTitleSize(tsize)
        h.GetYaxis().SetLabelSize(lsize)
 
        if h.GetZaxis():
            h.GetZaxis().SetTitleSize(zsize)
            h.GetZaxis().SetLabelSize(zsize)
 

setAxisTitles()

PlotCalibrationGains.setAxisTitles	(		histos,
			xtit,
			ytit,
			ztit = "foo",
			xoff = 0,
			yoff = 0,
			zoff = 0 )

Definition at line 25 of file PlotCalibrationGains.py.

def setAxisTitles(histos, xtit, ytit, ztit="foo", xoff=0, yoff=0, zoff=0):
    for h in histos:
        h.GetXaxis().SetTitle(xtit)
        h.GetYaxis().SetTitle(ytit)
        if ztit != "foo":
            h.GetZaxis().SetTitle(ztit)
            
        if xoff!=0: h.GetXaxis().SetTitleOffset(xoff)
        if yoff!=0: h.GetYaxis().SetTitleOffset(yoff)
        if zoff!=0: h.GetZaxis().SetTitleOffset(zoff)
 
 

Variable Documentation

action

PlotCalibrationGains.action

Definition at line 985 of file PlotCalibrationGains.py.

args

PlotCalibrationGains.args

Definition at line 994 of file PlotCalibrationGains.py.

default

PlotCalibrationGains.default

Definition at line 992 of file PlotCalibrationGains.py.

dest

PlotCalibrationGains.dest

Definition at line 985 of file PlotCalibrationGains.py.

help

PlotCalibrationGains.help

Definition at line 985 of file PlotCalibrationGains.py.

options

PlotCalibrationGains.options

Definition at line 994 of file PlotCalibrationGains.py.

parser

PlotCalibrationGains.parser = OptionParser()

Definition at line 984 of file PlotCalibrationGains.py.

type

PlotCalibrationGains.type

Definition at line 985 of file PlotCalibrationGains.py.