diff options
| -rw-r--r-- | psr/psr_vs_detuning_combo.m | 50 |
1 files changed, 46 insertions, 4 deletions
diff --git a/psr/psr_vs_detuning_combo.m b/psr/psr_vs_detuning_combo.m index 1016a82..18ad5cc 100644 --- a/psr/psr_vs_detuning_combo.m +++ b/psr/psr_vs_detuning_combo.m @@ -1,5 +1,6 @@ 1; +data_dir='results/'; N_detun_steps=1000; detuning_p_min=-200.0; %detuning_p_max=-detuning_p_min; @@ -27,7 +28,9 @@ zeeman_splitting=+0.000; 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); +%[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) +load(fname) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % 0.1 G magnetic field,, 30 degree ellipticity @@ -36,7 +39,9 @@ zeeman_splitting=+0.070; 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); +%[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) +load(fname) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % 0.0001 G magnetic field,, 30 degree ellipticity @@ -45,7 +50,9 @@ zeeman_splitting=+0.000070; 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); +%[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) +load(fname) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 % 1.0 G magnetic field,, 30 degree ellipticity @@ -54,5 +61,40 @@ zeeman_splitting=+0.70; 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); +%[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) +load(fname) + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 +% -0.1 G magnetic field,, 30 degree ellipticity + +zeeman_splitting=-0.070; +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) +load(fname) + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 +% -0.0001 G magnetic field,, 30 degree ellipticity + +zeeman_splitting=-0.000070; +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) +load(fname) + +%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5 +% -1.0 G magnetic field,, 30 degree ellipticity + +zeeman_splitting=-0.70; +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) +load(fname) |