diff --git a/macro/mcbm/pl_lambda.C b/macro/mcbm/pl_lambda.C index 50797a37e64bc31d80f2ec3bce4d4eef0f117a83..207f856d8d2f5b5e507a5f4d4d9470e272e9eaf0 100755 --- a/macro/mcbm/pl_lambda.C +++ b/macro/mcbm/pl_lambda.C @@ -199,10 +199,10 @@ can->cd(10); if(NULL != fSignal) fSignal->Delete(); fSignal=new TF1("fSignal","gaus",1.090,1.140); //h1->Fit("gaus","","H",1.1,1.2); - Double_t MMin=1.09; - Double_t MMax=1.142; + Double_t MMin=1.1; + Double_t MMax=1.14; Double_t param[3]; - param[0]=h1->GetMaximum(); + param[0]=h1->GetMaximum()*0.5; param[1]=1.1156; param[2]=0.01; fSignal->SetParameters(param); @@ -214,7 +214,10 @@ can->cd(10); Double_t dFRange=2.5; MMin=MMEAN-dFRange*MSIG; MMax=MMEAN+dFRange*MSIG; - cout<<" Fit results: "<<MMEAN<<", "<<MSIG<<", new range: "<<MMin<<" - "<<MMax<<endl; + cout<<" 1. fit results: "<<MMEAN<<", "<<MSIG<<", new range: "<<MMin<<" - "<<MMax<<endl; + param[0]=h1->GetMaximum(); + param[1]=1.1156; + fSignal->SetParameters(param); h1->Fit("fSignal","","",MMin,MMax); fSignal->GetParameters(param); MMEAN=param[1];