Hi Sabine,
Look at the two lines modified below
Rene Brun
{
gROOT->Reset();
const Int_t np = 49;
Float_t x[np] = {1.913521, 1.953769, 2.347435, 2.883654,
3.493567,4.047560,4.337210, 4.364347, 4.563004, 5.054247,5.194183,
5.380521, 5.303213, 5.384578, 5.563983,5.728500, 5.685752, 5.080029,
4.251809, 3.372246,2.207432, 1.227541, 0.8597788,0.8220503,
0.8046592,0.7684097,0.7469761,0.8019787,0.8362375,0.8744895,0.9143721,0.9462768,0.9285364,0.8954604,0.8410891,0.7853871,0.7100883,
0.6938808,0.7363682,0.7032954,0.6029015,0.5600163,0.7477068,1.188785,
1.938228,2.602717, 3.472962, 4.465014, 5.177035};
h = new TH1F("g1","Example of several fits in subranges",np,85,134);
h->SetMaximum(7);
for (int i=0;i<np;i++) {
h->SetBinContent(i+1,x[i]);
}
Double_t par[9];
g1 = new TF1("g1","gaus",85,95);
g2 = new TF1("g2","gaus",98,108);
TCanvas *myc = new TCanvas("myc","myc_titre",0,0,360,520);
myc->Clear();
myc->Divide(1,2);
myc->cd(1);
h->Fit("g1","","0",85,95); //<=========
h->DrawCopy(); //<=========
myc->cd(2);
h->Fit("g2","","",98,108);
myc->Update();
myc->Print("fit.eps");
}
On Tue, 6 Aug 2002, Sabine Crepe-Renaudin wrote:
>
>
> Hello,
>
>
> I'm using something equivalent to the following routine
> which simply fit an histo in two ranges and display each fit in 2 pads.
> Everything is OK till I try to print the result.
>
> You can try the routine below to see the effect. If the print
> command is commented, you see the 2 fits, if you uncomment the last line,
> then the first fit is overwritten with the second and this is what is
> saved in the postscript file ...
>
> Have I done something wrong ? is it possible to print the results
> correctly without copying the histo in an other one and fit 2 different
> histos ???
>
> Thanks in advance for any suggestion,
>
> Sabine.
>
>
> {
> gROOT->Reset();
> const Int_t np = 49;
> Float_t x[np] = {1.913521, 1.953769, 2.347435, 2.883654, 3.493567,4.047560,4.337210, 4.364347, 4.563004, 5.054247,5.194183, 5.380521, 5.303213, 5.384578, 5.563983,5.728500, 5.685752, 5.080029, 4.251809, 3.372246,2.207432, 1.227541, 0.8597788,0.8220503,
> 0.8046592,0.7684097,0.7469761,0.8019787,0.8362375,0.8744895,0.9143721,0.9462768,0.9285364,0.8954604,0.8410891,0.7853871,0.7100883,0.6938808,0.7363682,0.7032954,0.6029015,0.5600163,0.7477068,1.188785, 1.938228,2.602717, 3.472962, 4.465014, 5.177035};
>
> h = new TH1F("g1","Example of several fits in subranges",np,85,134);
> h->SetMaximum(7);
>
> for (int i=0;i<np;i++) {
> h->SetBinContent(i+1,x[i]);
> }
>
> Double_t par[9];
> g1 = new TF1("g1","gaus",85,95);
> g2 = new TF1("g2","gaus",98,108);
>
>
> TCanvas *myc = new TCanvas("myc","myc_titre",0,0,360,520);
> myc->Clear();
> myc->Divide(1,2);
> myc->cd(1);
> h->Fit("g1","","",85,95);
> myc->Update();
> myc->cd(2);
> h->Fit("g2","","",98,108);
> myc->Update();
>
> // myc->Print("fit.eps");
> }
>
>
>
> --------------------------------------------------------------------
> Sabine Crépé-Renaudin
> Expérience D0 (Fermilab)
> ISN
> --------------------------------------------------------------------
>
>
This archive was generated by hypermail 2b29 : Sat Jan 04 2003 - 23:51:03 MET