Baur MC

Combining all Quarks

.L ntread.C+
TNtuple* botnt = ntread("results/bot_1_0.out")
TH1F* bottom = new TH1F("bottom","Pt of first Z; events [2 GeV]",250,0,500)
botnt->Draw("PtZ1>>bottom","W")
TNtuple* chmnt = ntread("results/chm_1_0.out")
TH1F* charm = new TH1F("charm","Pt of first Z; events [2 GeV]",250,0,500)
chmnt->Draw("PtZ1>>charm","W")
TNtuple* strnt = ntread("results/strv_1_0")
TH1F* strange = new TH1F("strange","Pt of first Z; events [2 GeV]",250,0,500)
strnt->Draw("PtZ1>>strange","W")
TNtuple* dwnnt = ntread("results/dsum_1_0.out")
TH1F* down = new TH1F("down","Pt of first Z; events [2 GeV]",250,0,500)    
dwnnt->Draw("PtZ1>>down","W")
TNtuple* dseant = ntread("results/dsea_1_0.out")
TH1F* downsea = new TH1F("downsea","Pt of first Z; events [2 GeV]",250,0,500)
dseant->Draw("PtZ1>>downsea","W")
TNtuple* upnt = ntread("results/usum_1_0.out")
TH1F* up = new TH1F("up","Pt of first Z; events [2 GeV]",250,0,500)    
upnt->Draw("PtZ1>>up","W")
TNtuple* upvnt = ntread("results/upv_1_0.out")
TH1F* upv = new TH1F("upv","Pt of first Z; events [2 GeV]",250,0,500)
upvnt->Draw("PtZ1>>upv","W")
TNtuple* useant = ntread("results/usea_1_0.out")
TH1F* upsea = new TH1F("upsea","Pt of first Z; events [2 GeV]",250,0,500)
useant->Draw("PtZ1>>upsea","W")
up->Scale(0.00540/up->Integral())
down->Scale(0.00128/down->Integral())
strange->Scale(0.000032/strange->Integral())
charm->Scale(0.00000427/charm->Integral())
bottom->Scale(0.00000163/bottom->Integral())
TH1F* Total = new TH1F("Total","Pt of first Z; events [2 GeV]",250,0,500)
Total->Add(up)
Total->Add(down)
Total->Add(charm)
Total->Add(strange)
Total->Add(bottom)
TNtuple* fullnt = ntread("results/full_1_0.out")
TH1F* full = new TH1F("full","Pt of first Z; events [2 GeV]",250,0,500)
fullnt->Draw("PtZ1>>full","W")
full->Scale(0.0067/full->Integral())

-- MatthewNorman - 08 Dec 2006