hi,
igot root to fit a gaussian over a triangular distribution, which
comes naturally out of the (ptl+ptr) coplanarity plots. here is a
preview; hopefully by tomorrow it works inside show_proton.C, and then
we will have a much better measure of real (e,ep) yields.
here's the commands:
//- positive triangle function
TH1F b("b","b", 100,-1000,1000)
eep->Draw("ptl+ptr>>b", "ntl==4&&ntr==12")
TF1 t("tri","[2]*((1-abs(x-[0])/[1])+abs(1-abs(x-[0])/[1]))/2")
tri->SetParameters(-100,300,40)
b->Fit("tri")
//- triangle + gaus
c=new TF1("tot","tri+gaus(3)",-1000,1000)
c->SetParameter(3,100)
c->SetParameter(4,-100)
c->SetParameter(5,50)
pv(c->GetParameters(),6)
c->FixParameter(1,316)
b->Fit("tot")
//- plot individual functions
d=c->GetParameters()
ga=new TF1("gaus")
ga=new TF1("ga","gaus")
ga->SetParameters(d+3)
ga->Print()
ga->SetRange(-1000,1000)
ga->Draw("same")
ta=new TF1("ta", "tri",-1000,1000)
ta->SetParameters(d)
ta->Draw("same")
--chris
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:28 EST
