Strange error in propagating from point source to detector

[capitanevs]

This code is not working. I did not have time to understand why.

`from LibWiser.EasyGo import * import LibWiser.FermiSource

%% DEFAULT SETTINGS: You can change parameter here or in a similar dictionary

%=============================================================================

BeamlineName = '' SettingsDefault = {#chosen to mimic a plane wave on slits at 70m 'Lambda' : 633e-9, 'Waist0' : 1e-15, 'DetectorSize' : 50e-3, 'NSamples' : 8000, 'DetectorDefocus' : 0, 'UseCustomSampling' : True, 'SlitSize' : 10e-6, 'SlitDistance' : 50, 'DeltaZ' : 100e-3 }

=============================================================================

S = SettingsDefault DetectorSize = S['DetectorSize']
NSamples = S['NSamples'] FresnelNumber = S['SlitSize']**2 / S['Lambda'] / S['DeltaZ'] print("FresnelNumber = %0.1e" % FresnelNumber)

%% LAYOUT: creating Foundation.Optical element objects

%=============================================================================

---SOURCE (H,V)

------------------------------------------------------------

DeltaSource = 0 s = OpticalElement( Name = 's', IsSource = True, CoreOpticsElement = Optics.SourceGaussian( Lambda = S['Lambda'], Waist0 = S['Waist0'], Orientation = Optics.OPTICS_ORIENTATION.ISOTROPIC), PositioningDirectives = Foundation.PositioningDirectives( ReferTo = 'absolute', XYCentre = [0,0], Angle = 0) )

---slits_v

------------------------------------------------------------

slits_v = OpticalElement( Name = 'slits_v', CoreOpticsElement = Optics.Slits( L = S['SlitSize'], Orientation = Optics.OPTICS_ORIENTATION.VERTICAL), PositioningDirectives = Foundation.PositioningDirectives( PlaceWhat = 'centre', PlaceWhere = 'centre', ReferTo = 'source', Distance = S['SlitDistance'] ) )

---Detector

------------------------------------------------------------

det_v= OpticalElement( Name = 'det_v', CoreOpticsElement = Optics.Detector( L = DetectorSize, AngleGrazing = np.deg2rad(90), Orientation = Optics.OPTICS_ORIENTATION.VERTICAL), PositioningDirectives = Foundation.PositioningDirectives( PlaceWhat = 'centre', PlaceWhere = 'centre', ReferTo = 'source', Distance = S['SlitDistance'] + S['DeltaZ']))

Create Beamlinee beamline starting from Beamlinee Beamline Elements Objects

------------------------------------------------------------------------------------------------

Beamline = None Beamline = Foundation.BeamlineElements() Beamline.Append(s)

Beamline.Append(slits_v)

Beamline.Append(det_v) Beamline.RefreshPositions() print(Beamline) # Beamline.Name = BeamlineName

%% Settings: UseFigureError

==========================================================

Beamline.ComputeFields() Beamline.PlotIntensity()

d = S['SlitSize'] z = S['DeltaZ'] Lambda = S['Lambda'] NF = d*2 /Lambda/z x0 = Lambda z / d print('NF = %0.2f' % NF)

R = s.CoreOptics.RCurvature(S['SlitDistance']) z_eff = (1/z + 1/R)*-1 x0_g = Lambda z_eff / d

%% Add constant field

E = det_v.ComputationData.Field plot(abs(E+10)**2)

`