Example #1
0
#!/usr/bin/env python

import sys, os
from CompareTreeVars import getHisto, CompareTreeVars
import ROOT as r
from ROOT import TTree, TH2D, TH1D, TString, gStyle, TMath, gROOT

gROOT.SetMacroPath("~/phi2KsKs/phi2KsKs/preliminaryStudies/");


from Style import rootStyle, printLHCb
rootStyle(False)

for track in ["LL","DD","DD-DTF", "LL-DTF","DeltaTauLL","DeltaTauLD","DeltaTauLL-DTF"]:
	f = r.TFile.Open("/afs/cern.ch/work/s/sbartkow/files/phi2KsKs_incl.root", "read")
	t = f.Get("TuplePhi2KsKs/DecayTree")
	gStyle.SetOptFit(111)
	#gStyle.SetOptTitle(0)
	gStyle.SetOptTitle(1)
	gStyle.SetOptStat("e")

	canvas = r.TCanvas("canvas","Canvas Title")
	canvas.cd()
	if track == "LL-DTF":
		t.Draw("(phi_DTF_CTAU_Ks1/299.792458-Ks1_TRUETAU)*1000. >>histo(50,0,.5)","phi_DTF_CTAU_Ks1 >= 0 && phi_BKGCAT==0 && pi1_TRACK_Type ==3 && pi2_TRACK_Type ==3")
		histo=r.gROOT.FindObject("histo")
		histo.SetDirectory(0)
		t.Draw("(phi_DTF_CTAU_Ks2/299.792458-Ks2_TRUETAU)*1000. >>histo2(50,0,.5)","phi_DTF_CTAU_Ks2 >= 0 && phi_BKGCAT==0 && pi3_TRACK_Type ==3 && pi3_TRACK_Type ==3")
		histo2=r.gROOT.FindObject("histo2")
		histo2.SetDirectory(0)
		histo.Add(histo2)
Example #2
0
#!/usr/bin/env python

import sys, os
from CompareTreeVars import getHisto, CompareTreeVars
import ROOT as r
from ROOT import TTree, TH2D, TH1D, TString, gStyle, TMath, TLatex, gROOT
import LHCbStyle
from Style import rootStyle

gROOT.SetMacroPath("~/phi2KsKs/phi2KsKs/preliminaryStudies/");
# gROOT.ProcessLine(".x LHCbStyle.C")
# # LHCbStyle()

from Style import rootStyle, printLHCb
lhcbStyle = rootStyle(False)


f = r.TFile.Open("/afs/cern.ch/work/s/sbartkow/files/data2012.root", "read")
t = f.Get("TuplePhi2KsKs/DecayTree")
# gStyle.SetOptFit(111)
# gStyle.SetOptTitle(0)

lhcbStyle.SetOptTitle(0)


t.Draw("Ks1_PT/1000. >>histo(100,0,10)","Ks1_TAU>0.005 && Ks1_TAU <0.1")#,"pi1_TRACK_Type==3 && pi2_TRACK_Type==3")
histo=r.gROOT.FindObject("histo")
histo.SetDirectory(0)
t.Draw("Ks2_PT/1000.>>histo2(100,0,10)","Ks2_TAU>0.005 && Ks2_TAU <0.1")#,"pi3_TRACK_Type==3 && pi4_TRACK_Type==3")
histo2=r.gROOT.FindObject("histo2")
histo2.SetDirectory(0)
Example #3
0
#!/usr/bin/env python

import sys, os
from CompareTreeVars import getHisto, CompareTreeVars
import ROOT as r
from ROOT import TTree, TH2D, TH1D, TString, gStyle, TMath, TLatex, gROOT
# from LHCbStyle import printLHCb

gROOT.SetMacroPath("~/phi2KsKs/phi2KsKs/preliminaryStudies/");


from Style import rootStyle, printLHCb
rootStyle(True)

for track in ["DD"]:#["DD","LL","DD-DTF", "LL-DTF","DeltaTauLL","DeltaTauLD", "DeltaTauLL-DTF"]:
	f = r.TFile.Open("/afs/cern.ch/work/s/sbartkow/files/Ds_Phi2KsKs_Ds.root", "read")
	t = f.Get("TuplePhi2KsKs/DecayTree")
	gStyle.SetOptFit(111)
	#gStyle.SetOptTitle(0)
	gStyle.SetOptTitle(1)
	gStyle.SetOptStat("e")

	canvas = r.TCanvas("canvas","Canvas Title")
	canvas.cd()
	if track == "LL-DTF":
		t.Draw("(phi_DTF_CTAU_Ks1_PV/299.792458-Ks1_TRUETAU)*1000. >>histo(50,0,5)","phi_DTF_CTAU_Ks1_PV >= 0 && Ds_DTF_CHI2_PV >0 && phi_BKGCAT==0 && pi1_TRACK_Type ==3 && pi2_TRACK_Type ==3")
		histo=r.gROOT.FindObject("histo")
		histo.SetDirectory(0)
		t.Draw("(phi_DTF_CTAU_Ks2_PV/299.792458-Ks2_TRUETAU)*1000. >>histo2(50,0,5)","phi_DTF_CTAU_Ks2_PV >= 0 && Ds_DTF_CHI2_PV >0 && phi_BKGCAT==0 && pi3_TRACK_Type ==3 && pi3_TRACK_Type ==3")
		histo2=r.gROOT.FindObject("histo2")
		histo2.SetDirectory(0)