Commit fc449d32 authored by Andrew Strong's avatar Andrew Strong

SourcePopulation.h add print_ctl, SourcePopulation.cc use print_ctl for cout...

SourcePopulation.h add print_ctl, SourcePopulation.cc use print_ctl for cout of sources, init print_crl=0, not yet in Galplotdef
parent 7f2f793c
...@@ -45,6 +45,8 @@ int SourcePopulation::init() ...@@ -45,6 +45,8 @@ int SourcePopulation::init()
IDL_plot_control = 1; //AWS20200623 IDL_plot_control = 1; //AWS20200623
IDL_psym = -6; //AWS20200623 IDL_psym = -6; //AWS20200623
print_ctl = 0; //AWS20200720
strcpy(galdef_ID, "sourcepop_needs_galdef_ID" ); //AWS20200701 strcpy(galdef_ID, "sourcepop_needs_galdef_ID" ); //AWS20200701
strcpy(psfile_tag,"sourcepop_needs_psfile_tag"); //AWS20200701 strcpy(psfile_tag,"sourcepop_needs_psfile_tag"); //AWS20200701
...@@ -790,12 +792,12 @@ int SourcePopulation::analyse_sample() ...@@ -790,12 +792,12 @@ int SourcePopulation::analyse_sample()
pointing pointing_=pointing(theta,phi); pointing pointing_=pointing(theta,phi);
if(verbose!=4)cout<<"source statistics sample i="<<i<<" Fermi sensitivity l="<<sample_l[i] <<" b="<<sample_b[i] if(print_ctl==1)cout<<"source statistics sample i="<<i<<" Fermi sensitivity l="<<sample_l[i] <<" b="<<sample_b[i]
<<" theta="<<theta<<" phi="<<phi; //AWS20200716 <<" theta="<<theta<<" phi="<<phi; //AWS20200720
// cout<<" interpolated value="<<healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_)<<endl; // cout<<" interpolated value="<<healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_)<<endl;
flux_detection_limit_sensitivity_file = healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_); //AWS20200713 was flux_detection_limit flux_detection_limit_sensitivity_file = healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_); //AWS20200713 was flux_detection_limit
if(verbose!=4)cout<<" flux detection limit_sensitivity_file="<<flux_detection_limit_sensitivity_file<<" sample flux="<<sample_flux[i]<<endl;//AWS20200716 if(print_ctl==1)cout<<" flux detection limit_sensitivity_file="<<flux_detection_limit_sensitivity_file<<" sample flux="<<sample_flux[i]<<endl;//AWS20200720
} }
dlnN_dlnS [ii]++; dlnN_dlnS [ii]++;
...@@ -1065,12 +1067,12 @@ int SourcePopulation::analyse_sample() ...@@ -1065,12 +1067,12 @@ int SourcePopulation::analyse_sample()
double phi = sample_l[i] *dtr; double phi = sample_l[i] *dtr;
pointing pointing_=pointing(theta,phi); pointing pointing_=pointing(theta,phi);
if(verbose!=4)cout<<"skymap generation source sample i="<<i<<" Fermi sensitivity l="<<sample_l[i] <<" b="<<sample_b[i] if(print_ctl==1)cout<<"skymap generation source sample i="<<i<<" Fermi sensitivity l="<<sample_l[i] <<" b="<<sample_b[i]
<<" theta="<<theta<<" phi="<<phi; //AWS20200716 <<" theta="<<theta<<" phi="<<phi; //AWS20200720
// cout<<" interpolated value="<<healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_)<<endl; // cout<<" interpolated value="<<healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_)<<endl;
flux_detection_limit_sensitivity_file = healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_); //AWS20200713 flux_detection_limit_sensitivity_file = healpix_skymap_Fermi_sensitivity.interpolated_value(pointing_); //AWS20200713
if(verbose!=4)cout<<" flux detection limit_sensitivity_file="<<flux_detection_limit_sensitivity_file<<" sample flux="<<sample_flux[i]<<endl;//AWS20200716 if(print_ctl==1)cout<<" flux detection limit_sensitivity_file="<<flux_detection_limit_sensitivity_file<<" sample flux="<<sample_flux[i]<<endl;//AWS20200720
} }
ii=int( (log10(sample_flux[i])-lnS_min)/dlnS ); ii=int( (log10(sample_flux[i])-lnS_min)/dlnS );
......
...@@ -162,6 +162,7 @@ class SourcePopulation ...@@ -162,6 +162,7 @@ class SourcePopulation
char* Fermi_sensitivity_file; //AWS20170112 char* Fermi_sensitivity_file; //AWS20170112
int verbose; int verbose;
int print_ctl; //AWS20200720
ofstream *txt_stream; ofstream *txt_stream;
......
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