Commit dd08bcda authored by Reza  ANSARI's avatar Reza ANSARI
Browse files

Extension intervalle permis pour le fit de phase au dela de [0,2Pi] pour...

Extension intervalle permis pour le fit de phase au dela de [0,2Pi] pour rendre le fit plus robuste, Reza 23/12/2018
parent 9d987114
......@@ -557,7 +557,7 @@ int do_p4_cxtrkfit(vector< vector<double> > & v_time_data, vector< vector< vecto
}
}
}
mFit.SetParam(0,"Phase",bestphase,M_PI/720.,0.,2.2*M_PI);
mFit.SetParam(0,"Phase",bestphase,M_PI/720.,-0.5*M_PI,2.5*M_PI);
cout << "2."<<ii+1<<" Scan param bestxi2_red="<<bestxi2/(double)(tot_npoints-(1+NTRK))<<" bestphase="<<Angle(bestphase).ToDegree()
<<" bestnpts="<<bestnpts<<" bestafact="<<bestafact<< " A= ";
v_phase[ii]=bestphase;
......@@ -592,6 +592,8 @@ int do_p4_cxtrkfit(vector< vector<double> > & v_time_data, vector< vector< vecto
<< " nstep="<<mFit.GetNStep() << " rc="<<rcfit<<endl;
ofr <<setw(4)<<ii+1<<" "<<setw(8)<<mFit.GetChi2Red()<<" ";
double phase=mFit.GetParm(0); double err_phase=mFit.GetParmErr(0);
if (phase<0.) phase += 2.*M_PI;
if (phase<2.*M_PI) phase -= 2.*M_PI;
cout <<"Phase= "<<setw(10)<<Angle(phase).ToDegree()<<" +/- "<<setw(10)<<Angle(err_phase).ToDegree()<<" deg."<<endl;
ofr <<setw(8)<<Angle(phase).ToDegree()<<" "<<setw(8)<<Angle(err_phase).ToDegree()<<" ";
v_phase[ii]=phase;
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment