12 if input_file_name
is None:
13 print (
"No input file name, assuming output.root")
14 input_file_name =
"output.root"
19 gStyle.SetCanvasColor(10)
20 gStyle.SetFrameFillColor(10)
25 c1 = TCanvas(
'c1',
'Example',200,10,700,500)
28 graphs = TFile(input_file_name)
29 graphs.cd(
'ADC/FineTime')
30 key_list = gDirectory.GetListOfKeys()
32 run_directory_name =
None
37 line = iii_str.split(
"\"")
40 run_directory_name = name
41 if run_directory_name
is None:
44 print (
"ERROR, the file doesn't contain run# directory!")
48 print (
"Entering directory: ", run_directory_name)
49 gDirectory.cd(run_directory_name)
52 c1.Print(
"CalibrationTimingPlots.ps[")
54 ppm_em_2d_profile_etaPhi_adc_fineTime.SetMinimum(-10.)
55 ppm_em_2d_profile_etaPhi_adc_fineTime.SetMaximum(10.)
56 ppm_em_2d_profile_etaPhi_adc_fineTime.Draw(
"colz")
57 c1.Print(
"CalibrationTimingPlots.ps")
60 ppm_had_2d_profile_etaPhi_adc_fineTime.SetMinimum(-10.)
61 ppm_had_2d_profile_etaPhi_adc_fineTime.SetMaximum(10.)
62 ppm_had_2d_profile_etaPhi_adc_fineTime.Draw(
"colz")
63 c1.Print(
"CalibrationTimingPlots.ps")
66 directory_content = gDirectory.GetListOfKeys()
69 gStyle.SetOptStat(111111)
71 for key
in directory_content:
72 print (key,
" Class name= ", key.GetClassName(),
" Name=", key.GetName())
73 if key.GetClassName() ==
'TDirectoryFile':
74 print (
"OK, I will create new histogram for this partition!")
75 print (
"Name of partition is: ", key.GetName())
76 partition_histos.append(
TH1F(key.GetName()+
'_his',
'Timing distribution for partition: '+key.GetName(),100,-10.,10.))
78 gDirectory.cd(
"/ADC/FineTime/"+run_directory_name+
"/"+key.GetName())
80 list_of_histos = gDirectory.GetListOfKeys()
81 for his_key
in list_of_histos:
83 timing_value=his_key.ReadObj().GetBinContent(2)
84 partition_histos[-1].Fill(timing_value)
87 for iHisto
in partition_histos:
89 c1.Print(
"CalibrationTimingPlots.ps")
91 c1.Print(
"CalibrationTimingPlots.ps]")
93 os.system(
"ps2pdf CalibrationTimingPlots.ps")