Commit 36d8503b authored by Reza  ANSARI's avatar Reza ANSARI
Browse files

Suite Debug impression summary des fits , Reza 04/04/2019

parent e04fa54e
......@@ -147,14 +147,6 @@ int main (int narg, char* arg[])
v_trk.push_back(tks);
}
if (do_cxfit) { // Resume resultat fit cross-cor
cout<<" ---------- Summary phase fitting for cross-correlations ---------- " << endl;
for(size_t i=0; i<v_acxd.size(); i++) {
cout << "================================= TrackSet["<<i+1<<"] ================================="<<endl;
v_acxd[i].PrintCxPhaseFitSummary(cout);
}
cout<<"=================================================================================="<<endl;
}
if (do_baselinefit || do_baselinesimplex) { // fit simultane des 6 cross-cor
string cxbofile = "cxb6_"+outfilename;
string cxbckfile = "cxb6_"+checkfilename;
......
......@@ -205,6 +205,7 @@ AcxDataSet & AcxDataSet::operator = (AcxDataSet const & a)
v_acbeams=a.v_acbeams; v_cxbeams=a.v_cxbeams;
v_RcFit_ac=a.v_RcFit_ac; v_xi2red_ac=a.v_xi2red_ac;
v_Ddish=a.v_Ddish; v_thetaant=a.v_thetaant; v_phiant=a.v_phiant;
v_err_Ddish=a.v_err_Ddish; v_err_thetaant=a.v_err_thetaant; v_err_phiant=a.v_err_phiant;
v_RcFit_cx=a.v_RcFit_cx; v_xi2red_cx=a.v_xi2red_cx;
v_phase=a.v_phase; v_phi_0=a.v_phi_0; v_a_phi=a.v_a_phi;
v_err_phi_0=a.v_err_phi_0; v_err_a_phi=a.v_err_a_phi;
......@@ -553,8 +554,9 @@ int ACxSetFitter::doACfit(string outfilename)
//DBG mFit.PrintFit();
// cout << "do_p4_trkfit 2."<<ii+1<<" Performing the fit for AutoCor Antenna= " << ii+1 << endl;
int rcfit = mFit.doFit();
double xi2red=mFit.GetChi2Red();
if (_prtlevel_>1) mFit.PrintFit();
v_RcFit_ac[ii]=rcfit; v_xi2red_ac[ii]=-9999.;
v_RcFit_ac[ii]=rcfit; v_xi2red_ac[ii]=xi2red;
if(rcfit>0) {
cout<< "------- Fit result for Antenna No="<<ii+1<<" Reduce_Chisquare = " << mFit.GetChi2Red()
<< " nstep="<<mFit.GetNStep() << " rc="<<rcfit<<endl;
......@@ -566,7 +568,6 @@ int ACxSetFitter::doACfit(string outfilename)
}
ofr <<setw(4)<<ii+1<<" "<<setw(8)<<mFit.GetChi2Red()<<" ";
v_xi2red_ac[ii]=mFit.GetChi2Red();
double Dfit=mFit.GetParm(0); double err_Dfit=mFit.GetParmErr(0);
cout <<setw(16)<<"DishDiameter= "<<setw(10)<<Dfit<<" +/- "<<setw(10)<<err_Dfit<<" m."<<endl;
ofr <<setw(5)<<rcfit<<" "<<setw(8)<<Dfit<<" "<<setw(8)<<err_Dfit<<" ";
......@@ -820,10 +821,9 @@ int ACxSetFitter::doCxfit(string outfilenamecx, bool useAac, bool fg_B0, bool fg
if (_prtlevel_>1)
cout << " 3."<<ii+1<<" Performing the fit for CrossCor " << ii << " FxF= " << Anum1[ii]+1<<"x"<<Anum2[ii]+1<<endl;
int rcfit = mFit.doFit();
v_RcFit_cx[ii]=rcfit; v_xi2red_cx[ii]=-99999.;
v_RcFit_cx[ii]=rcfit; v_xi2red_cx[ii]=mFit.GetChi2Red();
if (_prtlevel_>1) mFit.PrintFit();
if(rcfit>0) {
v_xi2red_cx[ii]=mFit.GetChi2Red();
// cout<< "-------------------------- Result for Cross No " << ii << endl;
cout<< "------ Fit result for Cross No "<<ii+1<<" Reduce_Chisquare = " << mFit.GetChi2Red()
<< " nstep="<<mFit.GetNStep() << " rc="<<rcfit<<endl;
......
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