PlotCalibrationGains Namespace Reference

Classes
class	EmPartitionPlots
class	GainReader
class	HadPartitionPlots
class	L1CaloGeometryConvertor
class	L1CaloMap

Functions
def	PlotCalibrationGains (input_file_name="", reference_file_name="", isInputXml=False, isInputSqlite=False, isRefXml=False, isRefSqlite=False, isRefOracle=False)

Variables
	parser
	action
	type
	dest
	help
	options
	args

Function Documentation

PlotCalibrationGains()

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

Definition at line 610 of file PlotCalibrationGains.py.

def PlotCalibrationGains(input_file_name="",reference_file_name="",isInputXml=False,isInputSqlite=False,isRefXml=False,isRefSqlite=False,isRefOracle=False):
 
  ROOT.gStyle.SetPalette(1)
  ROOT.gStyle.SetOptStat(111111)
  ROOT.gStyle.SetCanvasColor(10)
  
  c1 = ROOT.TCanvas('c1','Example',200,10,700,500)
  c2 = ROOT.TCanvas('c2','Example Partitions',200,10,700,500)
  c2.Divide(3,2)
 
  h_gains_em  = L1CaloMap("Eta-phi map of EM gains","#eta bin","#phi bin")
  h_gains_had = L1CaloMap("Eta-phi map of HAD gains","#eta bin","#phi bin")
 
  h_gains_em_fselect  = L1CaloMap("Eta-phi map of EM gains that failed selection","#eta bin","#phi bin")
  h_gains_had_fselect = L1CaloMap("Eta-phi map of HAD gains that failed selection","#eta bin","#phi bin")
 
  h_chi2_em  = L1CaloMap("Eta-phi map of EM Chi2","#eta bin","#phi bin")
  h_chi2_had = L1CaloMap("Eta-phi map of HAD Chi2","#eta bin","#phi bin")
 
  h_offset_em  = L1CaloMap("Eta-phi map of EM offsets","#eta bin","#phi bin")
  h_offset_had = L1CaloMap("Eta-phi map of HAD offsets","#eta bin","#phi bin")
 
 
  h_unfitted_em  = L1CaloMap("Eta-phi map of EM failed fits","#eta bin","#phi bin")
  h_unfitted_had = L1CaloMap("Eta-phi map of HAD failed fits","#eta bin","#phi bin")
 
  h_drifted_em  = L1CaloMap("EM TTs that drifted more then 10 %","#eta bin","#phi bin")
  h_drifted_had = L1CaloMap("HAD TTs that drifted more then 10 %","#eta bin","#phi bin")
 
 
  h_gains_em_reference  = L1CaloMap("Eta-phi map of EM gains (gain-reference)","#eta bin","#phi bin")
  h_gains_had_reference = L1CaloMap("Eta-phi map of HAD gains (gain-reference)","#eta bin","#phi bin")
 
  h_gains_em_reference_rel  = L1CaloMap("Eta-phi map of EM gains: (gain-reference)/reference","#eta bin","#phi bin")
  h_gains_had_reference_rel = L1CaloMap("Eta-phi map of HAD gains: (gain-reference)/reference","#eta bin","#phi bin")
 
 
  em_partition_gains = EmPartitionPlots(" EM gains",40,0.6,1.4,"gain","N")
  had_partition_gains = HadPartitionPlots(" HAD gains",40,0.5,2.5,"gain","N")
 
  em_partition_gains_ref = EmPartitionPlots(" EM gains - reference",40,-0.2,0.2,"gain-reference","N")
  had_partition_gains_ref = HadPartitionPlots(" HAD gains - reference",40,-1.,1.,"gain-reference","N")
 
  em_partition_gains_ref_rel = EmPartitionPlots(" EM gains-reference / reference",40,-0.2,0.2,"(gain-reference)/reference","N")
  had_partition_gains_ref_rel = HadPartitionPlots(" HAD gains-reference / reference",40,-1.,1.,"(gain-reference)/reference","N")
 
  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))
 
 
 
#print (measured_gains   )
  c1.cd()
  ROOT.gPad.SetLogy(0)
 
  h_gains_em.SetMinimum(0.6)
  h_gains_em.SetMaximum(1.4)
#  h_gains_em.SetMaximum(2.1)
  h_gains_em.Draw()
  c1.Print("Gains.ps(")
 
  h_gains_had.SetMinimum(0.6)
  h_gains_had.SetMaximum(1.4)
  h_gains_had.Draw()
  c1.Print("Gains.ps")
 
  h_drifted_em.Draw()
  c1.Print("Gains.ps")
  
  h_drifted_had.Draw()
  c1.Print("Gains.ps")
 
 
  c1.cd()
  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()
  c1.Print("Gains.ps")
 
  h_gains_had_reference_rel.SetMinimum(-0.5)
  h_gains_had_reference_rel.SetMaximum(0.5)
  h_gains_had_reference_rel.Draw()
  c1.Print("Gains.ps")
 
#  ROOT.gPad.SetRightMargin(0.01)
  h_gains_em_reference_rel.SetMinimum(-0.1)
  h_gains_em_reference_rel.SetMaximum(0.1)
  h_gains_em_reference_rel.Draw()
  c1.Update()
  palette = h_gains_em_reference_rel.h_1.GetListOfFunctions().FindObject("palette")
  if palette:
    palette.SetLabelSize(0.025)
 
  h_gains_em_reference_rel.Draw()
  c1.Print("Gains.ps")
 
  h_gains_had_reference_rel.SetMinimum(-0.1)
  h_gains_had_reference_rel.SetMaximum(0.1)
  h_gains_had_reference_rel.Draw()
  c1.Update()
 
  palette = h_gains_had_reference_rel.h_1.GetListOfFunctions().FindObject("palette")
  if palette:
    palette.SetLabelSize(0.025)
    
  h_gains_had_reference_rel.Draw()
  c1.Print("Gains.ps")
 
 
  c1.cd()
  ROOT.gPad.SetLogy(0)
 
  h_gains_em_reference.SetMinimum(-0.5)
  h_gains_em_reference.SetMaximum(0.5)
 
  h_gains_em_reference.Draw()
  c1.Print("Gains.ps")
 
  h_gains_had_reference.SetMinimum(-0.5)
  h_gains_had_reference.SetMaximum(0.5)
 
#  h_gains_had_reference.SetMinimum(-0.2)
#  h_gains_had_reference.SetMaximum(0.2)
 
  h_gains_had_reference.Draw()
  c1.Print("Gains.ps")
 
  c1.cd()
  ROOT.gPad.SetLogy(0)
  h_chi2_em.SetMinimum(0.1)
  h_chi2_em.SetMaximum(100)
  h_chi2_em.Draw()
  c1.Print("Gains.ps")
 
  c1.cd()
  ROOT.gPad.SetLogy(0)
  h_chi2_had.SetMinimum(0.1)
  h_chi2_had.SetMaximum(100)
  h_chi2_had.Draw()
  c1.Print("Gains.ps")
 
  c1.cd()
  ROOT.gPad.SetLogy(0)
  h_offset_em.SetMinimum(-1.)
  h_offset_em.SetMaximum(1.)
  h_offset_em.Draw()
  c1.Print("Gains.ps")
 
  c1.cd()
  ROOT.gPad.SetLogy(0)
  h_offset_had.SetMinimum(-1.)
  h_offset_had.SetMaximum(1.)
  h_offset_had.Draw()
  c1.Print("Gains.ps")
 
 
  #c2.cd()
  c2.Clear()
  c2.Divide(3,2)
  for i_p in range(0,em_partition_gains.nPartitions):
    c2.cd(i_p+1)
    if em_partition_gains.his_partitions[i_p].GetEntries() > 0:
      ROOT.gPad.SetLogy()
    else:
      ROOT.gPad.SetLogy(0)
    em_partition_gains.his_partitions[i_p].Draw()
  
  c2.Print("Gains.ps")
  
  #c2.cd()
  c2.Clear()
  c2.Divide(3,2)
  for i_p in range(0,had_partition_gains.nPartitions):
    c2.cd(i_p+1)
    if had_partition_gains.his_partitions[i_p].GetEntries() > 0:
      ROOT.gPad.SetLogy()
    else:
      ROOT.gPad.SetLogy(0)
    had_partition_gains.his_partitions[i_p].Draw()
  
  c2.Print("Gains.ps")
  
  
  #c2.cd()
  c2.Clear()
  c2.Divide(3,2)
  for i_p in range(0,em_partition_gains_ref.nPartitions):
    c2.cd(i_p+1)
    if em_partition_gains_ref.his_partitions[i_p].GetEntries() > 0:
      ROOT.gPad.SetLogy()
    else:
      ROOT.gPad.SetLogy(0)
    em_partition_gains_ref.his_partitions[i_p].Draw()
  
  c2.Print("Gains.ps")
  
  #c2.cd()
  c2.Clear()
  c2.Divide(3,2)
  for i_p in range(0,had_partition_gains_ref.nPartitions):
    c2.cd(i_p+1)
    if had_partition_gains_ref.his_partitions[i_p].GetEntries() > 0:
      ROOT.gPad.SetLogy()
    else:
      ROOT.gPad.SetLogy(0)
    had_partition_gains_ref.his_partitions[i_p].Draw()
  
  c2.Print("Gains.ps")
 
  #c2.cd()
  c2.Clear()
  c2.Divide(3,2)
  for i_p in range(0,em_partition_gains_ref_rel.nPartitions):
    c2.cd(i_p+1)
    if em_partition_gains_ref_rel.his_partitions[i_p].GetEntries() > 0:
      ROOT.gPad.SetLogy()
    else:
      ROOT.gPad.SetLogy(0)
    em_partition_gains_ref_rel.his_partitions[i_p].Draw()
  
  c2.Print("Gains.ps")
  
  #c2.cd()
  c2.Clear()
  c2.Divide(3,2)
  for i_p in range(0,had_partition_gains_ref_rel.nPartitions):
    c2.cd(i_p+1)
    if had_partition_gains_ref_rel.his_partitions[i_p].GetEntries() > 0:
      ROOT.gPad.SetLogy()
    else:
      ROOT.gPad.SetLogy(0)
    had_partition_gains_ref_rel.his_partitions[i_p].Draw()
  
  c2.Print("Gains.ps")
 
 
  c1.cd()
  ROOT.gPad.SetLogy(0)
 
  h_gains_em_fselect.Draw()
  c1.Print("Gains.ps")
 
  h_gains_had_fselect.Draw()
  c1.Print("Gains.ps")
 
 
  h_unfitted_em.Draw()
  c1.Print("Gains.ps")
 
  h_unfitted_had.Draw()
  c1.Print("Gains.ps)")
 
  os.system("ps2pdf Gains.ps")
  
  bad_gain_file.close()
  drifted_towers_file.close()  
 
  print ("finished!")
 

Variable Documentation

action

PlotCalibrationGains.action

Definition at line 990 of file PlotCalibrationGains.py.

args

PlotCalibrationGains.args

Definition at line 1000 of file PlotCalibrationGains.py.

dest

PlotCalibrationGains.dest

Definition at line 990 of file PlotCalibrationGains.py.

help

PlotCalibrationGains.help

Definition at line 990 of file PlotCalibrationGains.py.

options

PlotCalibrationGains.options

Definition at line 1000 of file PlotCalibrationGains.py.

parser

PlotCalibrationGains.parser

Definition at line 988 of file PlotCalibrationGains.py.

type

PlotCalibrationGains.type

Definition at line 990 of file PlotCalibrationGains.py.