Commit 9ceb95d8 authored by Andrew Strong's avatar Andrew Strong

SourcePopulation.cc: plot model above and below sensivitity limit,...

SourcePopulation.cc: plot model above and below sensivitity limit, plot_control=3, gen_source_population.cc use this
parent cd129a14
......@@ -1579,6 +1579,25 @@ int SourcePopulation::analyse_sample()
*txt_stream<<endl;
*txt_stream<<"NS_soplimit = [";
for (ii=0;ii< n_dlnN_dlnS;ii++)
{ *txt_stream<<dlnN_dlnS_soplimit[ii];
if(ii<n_dlnN_dlnS-1) *txt_stream<<",";
}
*txt_stream<<"]";
*txt_stream<<"; idl format";
*txt_stream<<endl;
*txt_stream<<"NS_sublimit = [";
for (ii=0;ii< n_dlnN_dlnS;ii++)
{ *txt_stream<<dlnN_dlnS_sublimit[ii];
if(ii<n_dlnN_dlnS-1) *txt_stream<<",";
}
*txt_stream<<"]";
*txt_stream<<"; idl format";
*txt_stream<<endl;
// *txt_stream<<"device,file='"<<title<<".ps' ;idl"<<endl;
......@@ -1601,6 +1620,16 @@ int SourcePopulation::analyse_sample()
*txt_stream<<"oplot,S,NS,psym="<<IDL_psym<<" ; idl"<<endl;
if(IDL_plot_control==3) // overplot on start plot
*txt_stream<<"oplot,S,NS,psym="<<IDL_psym<<" ; idl"<<endl;
if(IDL_plot_control==3) // overplot on start plot
*txt_stream<<"oplot,S,NS_soplimit,psym="<<IDL_psym<<" ; idl"<<endl;
if(IDL_plot_control==3 ) // overplot on start plot
*txt_stream<<"oplot,S,NS_sublimit,psym="<<IDL_psym<<" ; idl"<<endl;
// python format
*txt_stream<<"title='" <<title<< "' # python format"<<endl;
......
......@@ -375,8 +375,8 @@ int Galplot::gen_source_population()
sourcepop5.print();
sourcepop5.IDL_plot_control=2; //AWS20200623 oplot on first plot
sourcepop5.IDL_psym =0; //AWS2020623 just line
sourcepop5.IDL_plot_control=3; //AWS20200714 oplot on first plot model NS, soplimit, sublimit
sourcepop5.IDL_psym =0; //AWS20200714 just line
sourcepop5.analyse_sample();
......
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