|
|
Here is the script.
%%%%%%%%%%%%%
obj = chi3D('tag','(chi2D->chi3D (C)T.Lang)');
obj('Nx',256,'xWindow',0.007,...
'Ny',1, 'yWindow',0.007,... %use Ny=1 for 2+1D, Ny=2^N for 3+1d
'Nt',1024, 'tWindow',1.2e-12,'shiftFreqWindow',8.58e+13,...
'Nz',500, 'propagationLength', 0.003,...
'GPUOn',1, 'singlePrecision',0, 'abortIntensity', inf,...
'showResultEach',10, 'showProgress',1);
obj( 'material','BBO','theta',0.407185,'phi',0,'plane','XZ','temperature',22);
obj('deff1',1.7e-12,'eoo',1,'oeo',0,'ooo',0,...
'deff2',1.4e-12,'oee',0,'eeo',0,'eee',0,...
'n2',4e-20,'SPMOn',0,'selfFocusingOn',0);
obj('walkoffON',1,'noiseOn',1,'noiseFactor',1);
obj({'pulse2','EorI',1.61e-05,'polarization','o','centerWavelength',5.15e-07,...
'pulseDuration',2.2e-13,'GD',0,'GDD',0,'TOD',0,'Phase',0,'beamShape_t','gauss',...
'beamRadius_x',0.0014,'beamRadius_y',0.0014,... %please choose a reasonable values in y
'shift_x',0,'shift_y',0,...
'alpha_x',0,'alpha_y',0,...
'slant_x',0,'slant_y',0,...
'radiusOfCurvature_x',0,'radiusOfCurvature_y',0,...
'beamShape_x',{'gauss'},'beamShape_y',{'gauss'}});
obj.removeDetectors('all'); %delete default detectors
obj.detectors.obs1.plotIntegratedProfiles=1;
obj.detectors.obs1.logScale=0;
obj.detectors.obs1.polarization='o';
obj.detectors.obs1.ftlim=[5.78e+14,5.86e+14];
obj.detectors.obs1.fxlim=[-1000,1000];
obj.detectors.obs1.fylim=[-inf,inf];
obj.detectors.obs1.spectralData=1;
obj.detectors.obs1.spacetimeData=1;
obj.detectors.obs1.peakpower=1;
obj.detectors.obs2.plotIntegratedProfiles=1;
obj.detectors.obs2.logScale=0;
obj.detectors.obs2.polarization='e';
obj.detectors.obs2.ftlim=[1.15e+15,1.20e+15];
obj.detectors.obs2.fxlim=[-1000,1000];
obj.detectors.obs2.fylim=[-inf,inf];
obj.detectors.obs2.spectralData=1;
obj.detectors.obs2.spacetimeData=1;
obj.detectors.obs2.peakpower=1;
[E_ftfxfy_o,E_ftfxfy_e]=run(obj,{'pulse2'});
Results=obj.simResults
|
Locked
