%%% Exemples pour un calcule de précession
%
% $Revision: 15 $
% $Author: lenertz $
% $Date: 2015-07-30 17:36:07 +0200 (jeu., 30 juil. 2015) $
% contact: lenertz@cea.fr


%% Déclaration des paramètres pour une modification plus simple de l'exemple
% voir doc de computePolarization pour plus d'informations sur les
% arguments

% paramètres optionnelles

H = [0.00,0,0.008]; % vecteur Hxyz en T
M = 0.6; % en T (peut aussi être un vecteur dans une direction différente de H)
MagMod = @wireVortexField_vector_V2; % déclaration de la fonction calculant le champ magnétique

% Paramètres optionnelles
longueur = 200e-6; %en m /!\ attention à choisir en accord avec MagMod
x_step = 0.1e-6; % en m
largeur = 400e-6; % en m
y_step = 10e-6; % en m
r = 100e-6; % en m
wavelength = 4e-10; % en m
wl_sigma = 1.e-10; % en m
wl_nbStep = 3; % en m
wl_mesure = 1; % en unité de wl_sigma
beam_width = 30e-6; % en m 
BO_save = false; % booléen
startPolar = 0.9; % entre 1 et -1
    %angPhase = 0:2*pi/10+0.001:2*pi; % en rad
    angPhase = 0; % en rad


Params = [longueur, x_step, largeur, y_step, r, wavelength, wl_sigma,...
          wl_nbStep, wl_mesure, beam_width, BO_save];


%% Execution
% Calcul
p = computePolarization(H,M,MagMod, Params, startPolar, angPhase,1,2);

% % Affichage du résultat en polarisation sur une même figure pour plusieurs
% % largeurs de faisceau:
% fig_polar = p.plotPolar(10e-6);
% p.plotPolar(20e-6, fig_polar);
% p.plotPolar(40e-6, fig_polar);
% p.plotPolar(100e-6, fig_polar);
