Skip to content
GitLab
Commit da6c69ac authored by Andrew Strong's avatar Andrew Strong
Browse files

gen_source_population.cc: output summary table for sources above and below... 

gen_source_population.cc: output summary table for sources above and below threshold from threshold file
parent 8e61542c
Hide whitespace changes
Inline Side-by-side
source/gen_source_population.cc
View file @ da6c69ac
......@@ -445,7 +445,7 @@ int Galplot::gen_source_population()
//---------------------------------------------------------
txt_stream<<endl<<"Comparison source synthesis : Fermi catalogue 1"<<endl; //AWS20081216
//sourcepop1 replaced by sourcepop5 AWS20111011
txt_stream<<"source counts for region "<<sourcepop5.long_min1 <<" < l < "<<sourcepop5.long_max1 <<", "<<sourcepop5.long_min2 <<" < l < "<<sourcepop5.long_max2<<" / "
<<sourcepop5. lat_min1 <<" < b < "<<sourcepop5.lat_max1 <<", "<< sourcepop5. lat_min2 <<" < b < "<<sourcepop5. lat_max2<<endl; // AWS20111011
......@@ -462,6 +462,8 @@ int Galplot::gen_source_population()
txt_stream << "sourcepop5.n_dlnN_dlnS="<< sourcepop5.n_dlnN_dlnS<<endl;//AWS20111011
txt_stream << "sourcepop4.n_dlnN_dlnS="<< sourcepop4.n_dlnN_dlnS<<endl;//AWS20111011
txt_stream << "sourcepop5 model all sources, Fermi catalogue 1 "<<endl;
for (ii=0;ii< sourcepop5.n_dlnN_dlnS;ii++)
{
if( sourcepop5.dlnN_dlnS[ii]>0 || sourcepop4.dlnN_dlnS[ii]>0 )
......@@ -476,6 +478,44 @@ int Galplot::gen_source_population()
<<endl;
}
txt_stream <<endl;
txt_stream << "sourcepop5: model sources above threshold, sourcepop4: Fermi catalogue 1"<<endl;
for (ii=0;ii< sourcepop5.n_dlnN_dlnS;ii++)
{
if( sourcepop5.dlnN_dlnS[ii]>0 || sourcepop4.dlnN_dlnS[ii]>0 )
txt_stream <<"ii="<<ii
<<" lnS="<< sourcepop5.lnS_min+ii*sourcepop5.dlnS <<"/"<< sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS
<<" S=" <<pow(10, sourcepop5.lnS_min+ii*sourcepop5.dlnS) <<"/"<<pow(10, sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS)
<<" N(S)= "<< sourcepop5.dlnN_dlnS_soplimit [ii]<<":" << sourcepop4.dlnN_dlnS [ii]
<<" N(>S)= "<< sourcepop5.dlnN_dlnS_int_soplimit[ii]<<":" << sourcepop4.dlnN_dlnS_int[ii]
<<" F(S)= "<< sourcepop5.FS_soplimit [ii]<<":" << sourcepop4.FS [ii]
<<" F(>S)= "<< sourcepop5.FS_int_soplimit [ii]<<":" << sourcepop4.FS_int [ii]
<<endl;
}
txt_stream <<endl;
txt_stream << "sourcepop5: model sources below threshold, sourcepop4: Fermi catalogue 1"<<endl;
for (ii=0;ii< sourcepop5.n_dlnN_dlnS;ii++)
{
if( sourcepop5.dlnN_dlnS[ii]>0 || sourcepop4.dlnN_dlnS[ii]>0 )
txt_stream <<"ii="<<ii
<<" lnS="<< sourcepop5.lnS_min+ii*sourcepop5.dlnS <<"/"<< sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS
<<" S=" <<pow(10, sourcepop5.lnS_min+ii*sourcepop5.dlnS) <<"/"<<pow(10, sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS)
<<" N(S)= "<< sourcepop5.dlnN_dlnS_sublimit [ii]<<":" << sourcepop4.dlnN_dlnS [ii]
<<" N(>S)= "<< sourcepop5.dlnN_dlnS_int_sublimit[ii]<<":" << sourcepop4.dlnN_dlnS_int[ii]
<<" F(S)= "<< sourcepop5.FS_sublimit [ii]<<":" << sourcepop4.FS [ii]
<<" F(>S)= "<< sourcepop5.FS_int_sublimit [ii]<<":" << sourcepop4.FS_int [ii]
<<endl;
}
// one-line summary statistics
txt_stream << " galdef_ID = " << galplotdef.galdef_ID;
txt_stream << " psfile_tag = " << galplotdef.psfile_tag;
......@@ -554,6 +594,45 @@ int Galplot::gen_source_population()
<<endl;
}
txt_stream <<endl;
txt_stream << "sourcepop5: model sources above threshold, sourcepop6: Fermi catalogue 2"<<endl;
for (ii=0;ii< sourcepop5.n_dlnN_dlnS;ii++)
{
if( sourcepop5.dlnN_dlnS[ii]>0 || sourcepop6.dlnN_dlnS[ii]>0 )
txt_stream <<"ii="<<ii
<<" lnS="<< sourcepop5.lnS_min+ii*sourcepop5.dlnS <<"/"<< sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS
<<" S=" <<pow(10, sourcepop5.lnS_min+ii*sourcepop5.dlnS) <<"/"<<pow(10, sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS)
<<" N(S)= "<< sourcepop5.dlnN_dlnS_soplimit [ii]<<":" << sourcepop6.dlnN_dlnS [ii]
<<" N(>S)= "<< sourcepop5.dlnN_dlnS_int_soplimit[ii]<<":" << sourcepop6.dlnN_dlnS_int[ii]
<<" F(S)= "<< sourcepop5.FS_soplimit [ii]<<":" << sourcepop6.FS [ii]
<<" F(>S)= "<< sourcepop5.FS_int_soplimit [ii]<<":" << sourcepop6.FS_int [ii]
<<endl;
}
txt_stream <<endl;
txt_stream << "sourcepop5: model sources below threshold, sourcepop6: Fermi catalogue 2"<<endl;
for (ii=0;ii< sourcepop5.n_dlnN_dlnS;ii++)
{
if( sourcepop5.dlnN_dlnS[ii]>0 || sourcepop6.dlnN_dlnS[ii]>0 )
txt_stream <<"ii="<<ii
<<" lnS="<< sourcepop5.lnS_min+ii*sourcepop5.dlnS <<"/"<< sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS
<<" S=" <<pow(10, sourcepop5.lnS_min+ii*sourcepop5.dlnS) <<"/"<<pow(10, sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS)
<<" N(S)= "<< sourcepop5.dlnN_dlnS_sublimit [ii]<<":" << sourcepop6.dlnN_dlnS [ii]
<<" N(>S)= "<< sourcepop5.dlnN_dlnS_int_sublimit[ii]<<":" << sourcepop6.dlnN_dlnS_int[ii]
<<" F(S)= "<< sourcepop5.FS_sublimit [ii]<<":" << sourcepop6.FS [ii]
<<" F(>S)= "<< sourcepop5.FS_int_sublimit [ii]<<":" << sourcepop6.FS_int [ii]
<<endl;
}
// one-line summary statistics
txt_stream << " galdef_ID = " << galplotdef.galdef_ID;
......@@ -633,6 +712,45 @@ int Galplot::gen_source_population()
}
txt_stream <<endl;
txt_stream << "sourcepop5: model sources above threshold, sourcepop7: Fermi catalogue 3"<<endl;
for (ii=0;ii< sourcepop5.n_dlnN_dlnS;ii++)
{
if( sourcepop5.dlnN_dlnS[ii]>0 || sourcepop7.dlnN_dlnS[ii]>0 )
txt_stream <<"ii="<<ii
<<" lnS="<< sourcepop5.lnS_min+ii*sourcepop5.dlnS <<"/"<< sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS
<<" S=" <<pow(10, sourcepop5.lnS_min+ii*sourcepop5.dlnS) <<"/"<<pow(10, sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS)
<<" N(S)= "<< sourcepop5.dlnN_dlnS_soplimit [ii]<<":" << sourcepop7.dlnN_dlnS [ii]
<<" N(>S)= "<< sourcepop5.dlnN_dlnS_int_soplimit[ii]<<":" << sourcepop7.dlnN_dlnS_int[ii]
<<" F(S)= "<< sourcepop5.FS_soplimit [ii]<<":" << sourcepop7.FS [ii]
<<" F(>S)= "<< sourcepop5.FS_int_soplimit [ii]<<":" << sourcepop7.FS_int [ii]
<<endl;
}
txt_stream <<endl;
txt_stream << "sourcepop5: model sources below threshold, sourcepop7: Fermi catalogue 3"<<endl;
for (ii=0;ii< sourcepop5.n_dlnN_dlnS;ii++)
{
if( sourcepop5.dlnN_dlnS[ii]>0 || sourcepop7.dlnN_dlnS[ii]>0 )
txt_stream <<"ii="<<ii
<<" lnS="<< sourcepop5.lnS_min+ii*sourcepop5.dlnS <<"/"<< sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS
<<" S=" <<pow(10, sourcepop5.lnS_min+ii*sourcepop5.dlnS) <<"/"<<pow(10, sourcepop5.lnS_min+(ii+1)*sourcepop5.dlnS)
<<" N(S)= "<< sourcepop5.dlnN_dlnS_sublimit [ii]<<":" << sourcepop7.dlnN_dlnS [ii]
<<" N(>S)= "<< sourcepop5.dlnN_dlnS_int_sublimit[ii]<<":" << sourcepop7.dlnN_dlnS_int[ii]
<<" F(S)= "<< sourcepop5.FS_sublimit [ii]<<":" << sourcepop7.FS [ii]
<<" F(>S)= "<< sourcepop5.FS_int_sublimit [ii]<<":" << sourcepop7.FS_int [ii]
<<endl;
}
// one-line summary statistics
txt_stream << " galdef_ID = " << galplotdef.galdef_ID;
......
Supports Markdown
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