From 9ceb95d86582c18c25b5b99e12f87de501d8df72 Mon Sep 17 00:00:00 2001
From: Andrew Strong <aws@olga4.opt.rzg.mpg.de>
Date: Tue, 14 Jul 2020 11:06:04 +0200
Subject: [PATCH] SourcePopulation.cc: plot model above and below sensivitity
limit, plot_control=3, gen_source_population.cc use this
---
source/SourcePopulation.cc | 29 +++++++++++++++++++++++++++++
source/gen_source_population.cc | 4 ++--
2 files changed, 31 insertions(+), 2 deletions(-)
diff --git a/source/SourcePopulation.cc b/source/SourcePopulation.cc
index e3eba4a..650342f 100644
--- a/source/SourcePopulation.cc
+++ b/source/SourcePopulation.cc
@@ -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;
diff --git a/source/gen_source_population.cc b/source/gen_source_population.cc
index 35f6b15..b14bf22 100644
--- a/source/gen_source_population.cc
+++ b/source/gen_source_population.cc
@@ -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();
--
GitLab