histo.GetXaxis().SetTitle("#Delta t(reconstructed)- #Delta t(true) [ps]");
elif track == "DeltaTauLD" :
histo.GetYaxis().SetTitle("MC events/ps")
histo.Fit(f,"I","",-200,200);
histo.GetXaxis().SetTitle("#Delta t(reconstructed)- #Delta t(true) [ps]");
histo.Draw()
#f.SetParameters(1200,0.001,1,0.000)
f.SetLineColor(2)
f.Draw("same")
canvas.Update()
printLHCb(x=0.35)
canvas.Print("time_res_incl/timeResolution-"+track+".pdf")
del t
#del f
del canvas
del histo
if not ("DeltaTauLL" in track or track == "DeltaTauLD"):
del histo2
dictionary = {
"LL1" : "&& pi1_TRACK_Type ==3 && pi2_TRACK_Type == 3",
"LL2" : "&& pi3_TRACK_Type ==3 && pi4_TRACK_Type ==3",
"DD1" : "&& pi1_TRACK_Type ==5 && pi2_TRACK_Type == 5",
"DD2" : "&& pi3_TRACK_Type ==5 && pi4_TRACK_Type ==5",
histo2=r.gROOT.FindObject("histo2")
histo2.SetDirectory(0)
histo.Add(histo2)
canvas = r.TCanvas("canvas","Canvas Title")
canvas.cd()
histo.Draw()
histo.GetXaxis().SetTitle("p_{T}(K_{S}) [GeV]");
histo.GetYaxis().SetTitle("[a.u.]");
# funct1 = r.TF1("f1","[0]*(pow(x*1000,[1])*exp(-1000*[2]*x))",0,10)
funct1 = r.TF1("f1","30*(pow(x*1000,1.52)*exp(-1000*2.19e-3*x))",0,10)
funct1.SetLineColor(2)
funct1.SetParLimits(0,0,50)
funct1.SetParLimits(1,0,10)
funct1.SetParLimits(2,0,10)
funct1.SetParameters(1,1.52)
funct1.SetParameters(2,2.19e-3)
funct1.SetParameters(0,.3)
# funct1.Draw("same")
# funct2 = r.TF1("f2",".2*(pow(x*1000,2.19)*exp(-1.52e-3*1000*x)) ",0,10);
# funct2.SetLineColor(2)
# funct2.Draw("same")
# histo.Fit(funct1,"I0","",0,10);
funct1.Draw("same")
printLHCb()
canvas.Print("mommodel.pdf")
c_m_phi_Ds = TCanvas("c_m_phi_Ds","m_phi",2);
l_m_phi_Ds = TLegend(0.6, 0.70, .85, .85);
c_m_Ds = TCanvas("c_m_Ds","m_Ds",2);
l_m_Ds_Ds = TLegend(0.6, 0.70, .85, .85);
c_m_phi_incl.cd();
h_m_phi_mc_incl.GetXaxis().SetTitle("m(#phi)[MeV]");
h_m_phi_mc_incl.SetLineColor(kRed);
h_m_phi_mc_incl.GetYaxis().SetTitle("Entries normalized to unit area");
h_m_phi_mc_incl.GetYaxis().SetTitleOffset(1.3);
h_m_phi_mc_incl.DrawNormalized("");
l_m_phi_incl.AddEntry(h_m_phi_mc_incl, "MC prompt #phi", "l");
h_m_phi_data_incl.SetLineColor(kBlue);
h_m_phi_data_incl.DrawNormalized("same");
l_m_phi_incl.AddEntry(h_m_phi_data_incl, "2012 data", "l");
l_m_phi_incl.Draw();
printLHCb(x=0.35)
c_m_phi_incl.Print("m_phi_incl.pdf");
c_m_phi_Ds.cd();
h_m_phi_mc_Ds.GetXaxis().SetTitle("m(#phi)[MeV]");
h_m_phi_mc_Ds.SetLineColor(kRed);
h_m_phi_mc_Ds.GetYaxis().SetTitle("Entries normalized to unit area");
h_m_phi_mc_Ds.GetYaxis().SetTitleOffset(1.5);
h_m_phi_mc_Ds.DrawNormalized("");
l_m_phi_Ds.AddEntry(h_m_phi_mc_Ds, "MC D_{S} #rightarrow #phi #pi", "l");
h_m_phi_data_Ds.SetLineColor(kBlue);
h_m_phi_data_Ds.DrawNormalized("same");
l_m_phi_Ds.AddEntry(h_m_phi_data_Ds, "2012 data", "l");
l_m_phi_Ds.Draw();
printLHCb(x=0.35)
c_m_phi_Ds.Print("m_phi_Ds.pdf");
c_m_Ds.cd();
