diff options
Diffstat (limited to 'psr')
| -rw-r--r-- | psr/make_representative_psr_vs_detuning_for_given_B_and_psi_el.m | 25 | ||||
| -rw-r--r-- | psr/psr_vs_detuning_combo.m | 65 |
2 files changed, 58 insertions, 32 deletions
diff --git a/psr/make_representative_psr_vs_detuning_for_given_B_and_psi_el.m b/psr/make_representative_psr_vs_detuning_for_given_B_and_psi_el.m index 127eba6..46d8262 100644 --- a/psr/make_representative_psr_vs_detuning_for_given_B_and_psi_el.m +++ b/psr/make_representative_psr_vs_detuning_for_given_B_and_psi_el.m @@ -24,30 +24,7 @@ function ... [psr_rad_grEp_neg_el]=psr_vs_detuning(detuning_freq, Ep,-psi_el, B_field, theta, phi) ; - figure(6); - tn_cl='blue'; - sm_cl='green'; - lg_cl='magenta'; - gt_cl='red'; - - plot ( ... - detuning_freq, psr_rad_tnEp_pos_el, '-;tn,pos el;', "linewidth",3, 'color',tn_cl - , detuning_freq, psr_rad_tnEp_neg_el, '-;tn,neg el;', "linewidth",1, 'color',tn_cl - , detuning_freq, psr_rad_smEp_pos_el, '-;sm,pos el;', "linewidth",3, 'color',sm_cl - , detuning_freq, psr_rad_smEp_neg_el, '-;sm,neg el;', "linewidth",1, 'color',sm_cl - , detuning_freq, psr_rad_lgEp_pos_el, '-;lg,pos el;', "linewidth",3, 'color',lg_cl - , detuning_freq, psr_rad_lgEp_neg_el, '-;lg,neg el;', "linewidth",1, 'color',lg_cl - , detuning_freq, psr_rad_grEp_pos_el, '-;gr,pos el;', "linewidth",3, 'color',gt_cl - , detuning_freq, psr_rad_grEp_neg_el, '-;gr,neg el;', "linewidth",1, 'color',gt_cl - ); - - str_title=sprintf("PSR. B field=%.5f Gauss, ellipticity=%.2f rad, theta=%.2f, phi=%.2f", B_field, psi_el, theta, phi); - title(str_title); - xlabel('detuning, MHz'); - ylabel('PSR, radians'); - fname=strcat('results/',str_title, '.pdf'); - print(fname); - fname=strcat('results/',str_title, '.mat'); + fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) save(fname, ... 'detuning_freq', 'B_field', 'theta', 'phi', 'psi_el', ... 'psr_rad_tnEp_pos_el', 'psr_rad_tnEp_neg_el', ... diff --git a/psr/psr_vs_detuning_combo.m b/psr/psr_vs_detuning_combo.m index 18ad5cc..b2ea9d3 100644 --- a/psr/psr_vs_detuning_combo.m +++ b/psr/psr_vs_detuning_combo.m @@ -1,6 +1,7 @@ 1; data_dir='results/'; +output_dir='results/'; N_detun_steps=1000; detuning_p_min=-200.0; %detuning_p_max=-detuning_p_min; @@ -29,8 +30,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, psi_el, theta, phi); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % 0.1 G magnetic field,, 30 degree ellipticity @@ -40,8 +48,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % 0.0001 G magnetic field,, 30 degree ellipticity @@ -51,8 +66,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % 1.0 G magnetic field,, 30 degree ellipticity @@ -62,8 +84,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % -0.1 G magnetic field,, 30 degree ellipticity @@ -73,8 +102,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % -0.0001 G magnetic field,, 30 degree ellipticity @@ -84,8 +120,15 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % -1.0 G magnetic field,, 30 degree ellipticity @@ -95,6 +138,12 @@ B_field=zeeman_splitting/gmg; psi_el=30/180*pi; %[psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, psr_rad_grEp_pos_el, psr_rad_grEp_neg_el] =make_representative_psr_vs_detuning_for_given_B_and_psi_el(detuning_freq, B_field, theta, phi, psi_el); -fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el) +fname=data_file_name('results/', 'PSR.','mat', B_field, theta,phi,psi_el); load(fname) - +ret=ouput_psr_vs_detuning_combo( ... + psr_rad_tnEp_pos_el, psr_rad_tnEp_neg_el, ... + psr_rad_smEp_pos_el, psr_rad_smEp_neg_el, ... + psr_rad_lgEp_pos_el, psr_rad_lgEp_neg_el, ... + psr_rad_grEp_pos_el, psr_rad_grEp_neg_el ... + , detuning_freq, B_field, theta, phi, psi_el ... + , output_dir ); |