/******************************************************************************* * McStas instrument definition URL=http://mcstas.risoe.dk * * Instrument: ILL_H5 * * %Identification * Written by: FARHI Emmanuel (farhi@ill.fr) * Date: May, 2011 * Origin:ILL * Release: McStas 2.0 * Version: $Revision: 1.0 $ * %INSTRUMENT_SITE: ILL * * The full H5 new cold guide at the ILL, with ThALES, D16, Super-Adam, IN15 * (sample only), D22WASP and CryoEDM * This is a proposed geometry for major H5 guide hall upgrade (since 2013). * The exact adopted H5 geometry may be actually different. * * %Description * * This model decribes the full new H5 cold guide at the ILL, with ThALES, WASP, D16, * Super-Adam, IN15, D22. * * H511:The IN15 Spin-echo spectrometer * is simulated with an incoming polarized beam, but not with a full spin-echo * description. Sample is Vanadium * H512:The D22 Large dynamic range small-angle diffractometer * is fully simulated * H521:The D16 Small momentum transfer diffractometer D16 is realistic * H521:The SuperADAM reflectometer is used in low angle diffraction mode. * H522: The WASP Spin-echo spectrometer * is simulated with an incoming polarized beam, but not with a full spin-echo * description. Sample is Vanadium * H523:The CryoEDM with its polarized beam * H53: The ThALES Cold neutron three-axis spectrometer IN14 * * For each instrument, a sample can be specified (liquid/powder/amorphous), * with monitoring of the scattering in angular (diffraction) and energy modes * (for spectroscopy). * * %Example: lambda=5 Detector: H5_I=1.03234e+14 * * %Parameters * lambda: central wavelength band for guide illumination [Angs] * dlambda: half width of guide wavelength band [Angs] * ThALES_lambda:ThALES monochromator setting wavelength. Usual 2.4 and 4.2 [Angs] * WASP_lambda: IN16 monochromator setting wavelength. Usual 3.3 and 6.3 [Angs] * D16_lambda: D16 monochromator setting wavelength. Usual 4.7 and 5.6 [Angs] * SADAM_lambda: SuperADAM monochromator setting wavelength. Usual 4.4 [Angs] * IN15_lambda: IN15 velocity selector setting wavelength [Angs] * D22_lambda: D22 velocity selector setting wavelength [Angs] * ThALES_sample:ThALES liquid/powder/amorphous sample [string] * D16_sample: D16 liquid/powder/amorphous sample [string] * SADAM_sample: SuperADAM liquid/powder/amorphous sample [string] * D22_sample: D22 liquid/powder/amorphous sample [string] * D22_collimation: D22 collimation length and sample-detector distance [m] * * %Link * The NoteDPT11 at the ILL * %L * The DPT/SMAE 11/070 WASP design report * %L * The DPT/SMAE 10/271 H5 design report * %L * Daily notes from K. Andersen about the H5 project * %L * Mirotron drawing MR-0656-000 for the IN15 V-mirror geometry * %End *******************************************************************************/ DEFINE INSTRUMENT ILL_H5_new(lambda=5, dlambda=4.5, ThALES_lambda=4.2, WASP_lambda=6.3, D16_lambda=5.6, SADAM_lambda=4.4, IN15_lambda=6.5, D22_lambda=4.5, D22_collimation=2, string ThALES_sample="Rb_liq_coh.sqw", string WASP_sample="Rb_liq_coh.sqw", string D16_sample="H2O_liq.qSq", string SADAM_sample="SiO2_quartza.laz", string D22_sample="H2O_liq.qSq", ThALES_RMV=-1,D16_RMV=-1,SADAM_RMV=-1,ThALES_RMH=-1) DECLARE %{ /* HCS (H5) source 'Measured' parameters */ double sT3=413.5,sI3=10.22e12; double sT2=145.8,sI2=3.44e13; double sT1=40.1 ,sI1=2.78e13; /* guide coating parameters */ double gR0 = 1; double gQc = 0.0216; double gAlpha = 4.07; double gW = 1.0/300.0; double Al_Thickness = 0.001; double gGap = 0.001; double ThALES_DM = 3.355; /* PG002 */ double ThALES_A1 = 0; double ThALES_L = 2.12; double WASP_L = 2.75; double D16_DM = 3.355; /* PG002 */ double D16_A1 = 0; double D16_L = 2.8; double SADAM_DM = 3.355; /* PG002 */ double SADAM_A1 = 0; double SADAM_L = 3.36; double CryoEDM_L = 0.3; double IN15_L = 3.7; double IN15_nu = 0; double D22_nu = 0; double WASP_nu = 0; %} USERVARS %{ double flag; %} INITIALIZE %{ double L; /* the effective guide divergence is 1.2*0.1*lambda*pi/180 for m=1.2 */ /* the effective inguide length for focusing is thus (0.12/2)/(1.2*0.1*lambda*PI/180) /* compute H53 parameters: ThALES */ ThALES_A1 = asin(ThALES_lambda/2/ThALES_DM)*RAD2DEG; L = 1/(1/ThALES_L+1/( (0.12/2)/(3*0.1*ThALES_lambda*PI/180) )); if (ThALES_RMV<0) ThALES_RMV = 2*L*sin(DEG2RAD*fabs(ThALES_A1)); if (ThALES_RMH<0) ThALES_RMH = 2*L/sin(DEG2RAD*fabs(ThALES_A1)); printf("%s: ThALES: A1=%g [deg] RMV=%g [m] RMH=%g [m] lambda=%g [Angs] sample=%s\n", NAME_INSTRUMENT, ThALES_A1, ThALES_RMV, ThALES_RMH, ThALES_lambda, ThALES_sample); /* compute H521 parameters: D16 */ D16_A1 = asin(D16_lambda/2/D16_DM)*RAD2DEG; L = 1/(1/D16_L+1/( (0.12/2)/(2*0.1*D16_lambda*PI/180) )); if (D16_RMV<0) D16_RMV = 2*L*sin(DEG2RAD*fabs(D16_A1)); printf("%s: D16: A1=%g [deg] RMV=%g [m] lambda=%g [Angs] sample=%s\n", NAME_INSTRUMENT, D16_A1, D16_RMV, D16_lambda, D16_sample); /* compute H521 parameters: SuperADAM */ SADAM_A1 = asin(SADAM_lambda/2/SADAM_DM)*RAD2DEG; L = 1/(1/SADAM_L+1/( (0.12/2)/(2*0.1*SADAM_lambda*PI/180) )); if (SADAM_RMV<0) SADAM_RMV= 2*L*sin(DEG2RAD*fabs(SADAM_A1)); printf("%s: SuperADAM: A1=%g [deg] RMV=%g [m] lambda=%g [Angs] sample=%s\n", NAME_INSTRUMENT, SADAM_A1, SADAM_RMV, SADAM_lambda, SADAM_sample); /* compute H522 parameters: WASP Velocity Selector */ WASP_nu = 3956*14.2*DEG2RAD/2/PI/WASP_lambda/0.5; printf("%s: WASP: nu=%g [rpm] lambda=%g [Angs] sample=%s\n", NAME_INSTRUMENT, WASP_nu*60, WASP_lambda, WASP_sample); /* compute H511 parameters: IN15 Velocity Selector */ IN15_nu = 3956*64.4*DEG2RAD/2/PI/IN15_lambda/0.3; printf("%s: IN15: nu=%g [rmp] lambda=%g [Angs] sample=%s\n", NAME_INSTRUMENT, IN15_nu*60, IN15_lambda, "Vanadium (polarized)"); /* compute H512 parameters: D22 Velocity Selector */ D22_nu = 3956*48.3*DEG2RAD/2/PI/D22_lambda/0.25; printf("%s: D22: nu=%g [rpm] lambda=%g [Angs] sample=%s\n", NAME_INSTRUMENT, D22_nu*60, D22_lambda, D22_sample); %} TRACE /* ================================ H5 common =============================== */ COMPONENT Origin = Progress_bar() AT(0,0,0) ABSOLUTE COMPONENT HCS = Source_gen( radius = 0.21/2, zdepth = 0.15, dist = 2.155, focus_xw = 0.170, focus_yh = 0.120, lambda0 = lambda, dlambda = dlambda, T1 = sT1, I1 = sI1, T2 = sT2, I2 = sI2, T3 = sT3, I3 = sI3, verbose = 1) AT (0, 0, 0) RELATIVE Origin COMPONENT HCS_Al = PowderN( reflections="Al.laz", xwidth=0.17, yheight=0.12, zdepth=Al_Thickness, p_interact=0.05, p_inc=0.005, p_transmit=0.9450, d_phi=10) AT (0,0,0.21) RELATIVE HCS COMPONENT COPY(HCS_Al) = COPY(HCS_Al) AT (0,0,0.61) RELATIVE HCS COMPONENT COPY(HCS_Al) = COPY(HCS_Al) AT (0,0,0.78) RELATIVE HCS COMPONENT COPY(HCS_Al) = COPY(HCS_Al) AT (0,0,0.92) RELATIVE HCS COMPONENT COPY(HCS_Al) = COPY(HCS_Al) AT (0,0,2.15) RELATIVE HCS /* m=2 common rectangular large guide. MAY USE left/right as m=0.65 (glass, cheaper) */ COMPONENT H5 = Monitor_nD(xwidth=0.17, yheight=0.12, options="x y, slit", bins=50) AT(0,0,2.155) RELATIVE HCS COMPONENT H5_rect = Guide_gravity(w1=0.170, h1=0.12, l=1.0, R0=gR0, Qc=gQc, alpha=gAlpha, m = 2, W=gW) AT(0,0,2.155) RELATIVE HCS /* distribute the beam along channels */ COMPONENT H53_origin = Monitor_nD(xwidth=0.06, yheight=0.120, options="dx limits=[-2 2] dy limits=[-2 2]", bins=25, restore_neutron=1) AT (0.17/2-0.06/2,0,1+gGap) RELATIVE H5 ROTATED (0,1.5,0) RELATIVE H5 JUMP H53_start WHEN SCATTERED COMPONENT H52_origin = COPY(H53_origin)(xwidth=0.062) AT (0.17/2-0.06/2-0.004-0.062,0,1+gGap) RELATIVE H5 ROTATED (0,0,0) RELATIVE H5 JUMP H52_start WHEN SCATTERED COMPONENT H51_origin = COPY(H53_origin)(xwidth=0.04) AT (-0.17/2+0.04/2,0,1+gGap) RELATIVE H5 ROTATED (0,-1.3269,0) RELATIVE H5 JUMP H51_start WHEN SCATTERED /* start of guide sections */ /* ================================ H53 ===================================== */ /* feeds: ThALES */ COMPONENT H53_start = Arm() AT (0,0,0) RELATIVE H53_origin COMPONENT H53_inpile = COPY(H5_rect)(w1=0.06, m=3, l=4.930-3.155) AT (0,0,0) RELATIVE PREVIOUS COMPONENT H53_Al = COPY(HCS_Al)(xwidth=0.06) AT (0,0,COMP_GETPAR(PREVIOUS, l)+0.01) RELATIVE PREVIOUS /* Straight guide H53 m=3 length=14968-4930=10 m */ COMPONENT H53_Obt = COPY(H5_rect)(w1=0.06, m=3, l=3) AT (0,0,Al_Thickness+0.015) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=COMP_GETPAR(PREVIOUS, w1)) AT (0,0,COMP_GETPAR(PREVIOUS, l)+0.01) RELATIVE PREVIOUS /* Gap 75 mm after BOG (barillet obturateur general): VS */ COMPONENT H53_Obt_Out = COPY(H5)(xwidth=0.06, restore_neutron=1) AT (0,0,0.04) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.06) AT (0,0,0.01) RELATIVE PREVIOUS COMPONENT H53_VSComC1 = COPY(H53_inpile)(l=7,nelements=7) AT (0,0,3+0.075) RELATIVE H53_Obt COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.06) AT (0,0,7+0.01) RELATIVE PREVIOUS /* Parabolic taper 60x120 -> 30x120 m=3 */ COMPONENT H53_Nose = Guide_tapering( w1 = 0.06, h1=0.12, linw = 0.0, loutw = 0.7, l=2.0, linh=0.0, louth = 0.0, option="parabolical", R0 = gR0, Qcx = gQc, Qcy = gQc,alphax = gAlpha, alphay = gAlpha, W = gW, mx = 3, my = 3, segno=20) AT (0,0,0.01) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.04) AT (0,0,2+0.01) RELATIVE PREVIOUS /* Monochromator at 18998-16698=2.3 m from nose exit window */ COMPONENT H53_ThALES_Monochromator_Cradle = COPY(H5)(options="x y", restore_neutron=1) AT (0,0,2+2.3) RELATIVE H53_Nose EXTEND %{ flag=0; %} /* ThALES monochromator PG002 w15*h12 7 blades vertically focusing. Sample at 2.12 m */ COMPONENT H53_ThALES_Monochromator = Monochromator_curved( width = 0.15, height = 0.12, NH=7, NV=7, RV=ThALES_RMV, RH=ThALES_RMH, DM=ThALES_DM, mosaich = 30, mosaicv = 30, r0 = 1, reflect="HOPG.rfl", transmit="HOPG.trm") AT (0,0,0) RELATIVE H53_ThALES_Monochromator_Cradle ROTATED (0,ThALES_A1,0) RELATIVE H53_ThALES_Monochromator_Cradle EXTEND %{ flag=SCATTERED; %} COMPONENT H53_ThALES_Monochromator_Jumper = Arm() AT (0,0,0) RELATIVE H53_ThALES_Monochromator_Cradle JUMP H53_ThALES_Transmit WHEN !flag COMPONENT H53_ThALES_Monochromator_Out = Arm() AT (0,0,0) RELATIVE H53_ThALES_Monochromator_Cradle ROTATED (0,2*ThALES_A1,0) RELATIVE H53_ThALES_Monochromator_Cradle COMPONENT H53_ThALES_Sample_Div = Monitor_nD(xwidth=0.05, yheight=0.05, options="dx limits=[-2 2], dy limits=[-2 2]", bins=50, restore_neutron=1) AT (0,0,ThALES_L) RELATIVE H53_ThALES_Monochromator_Out COMPONENT H53_ThALES_Sample_XY = Monitor_nD(xwidth=0.10, yheight=0.10, options="x y", bins=50, restore_neutron=1) AT (0,0,ThALES_L) RELATIVE H53_ThALES_Monochromator_Out COMPONENT H53_ThALES_Sample_L = Monitor_nD(xwidth=0.10, yheight=0.10, options="lambda limits=[1 10]", bins=50, restore_neutron=1) AT (0,0,ThALES_L) RELATIVE H53_ThALES_Monochromator_Out COMPONENT H53_ThALES_Sample = Isotropic_Sqw(radius=0.005,yheight=0.05, Sqw_coh=ThALES_sample, Sqw_inc=NULL, d_phi=30) AT (0,0,ThALES_L) RELATIVE H53_ThALES_Monochromator_Out EXTEND %{ if (!SCATTERED) ABSORB; %} COMPONENT H53_ThALES_Spectrometer = Monitor_nD(radius=0.5, yheight=0.3, options="angle limits=[-150 150] bins=50, energy limits=[0 30], banana", restore_neutron=1) AT (0,0,0) RELATIVE H53_ThALES_Sample COMPONENT H53_ThALES_Diffractometer = Monitor_nD(radius=1, yheight=0.3,options="theta bins=100 limits=[-150 150], y, banana") AT (0,0,0) RELATIVE H53_ThALES_Sample EXTEND %{ ABSORB; /* neutron ends on the ThALES sample */ %} /* For fun we add the ThALES 2nd spectrometer */ COMPONENT Sample_Out = Arm() /* this is the sample-ana axis */ AT (0,0,0) RELATIVE H53_ThALES_Sample ROTATED (0, -42, 0) RELATIVE H53_ThALES_Monochromator_Out COMPONENT SC3 =Collimator_linear( xmin =-0.06/2, ymin =-0.12/2, xmax = 0.06/2, ymax = 0.12/2, length = 0.40) AT (0, 0, 0.5) RELATIVE Sample_Out COMPONENT Ana_Cradle = Arm() AT (0, 0, 1.37) RELATIVE Sample_Out COMPONENT PG2Xtal = Monochromator_curved( width = .2, height = .1, NH=11, NV=1, RV=0, RH=0, DM=3.355, mosaich = 30, mosaicv = 30, r0 = 0.7) AT (0, 0, 0) RELATIVE Ana_Cradle ROTATED (0, 21, 0) RELATIVE Ana_Cradle COMPONENT Ana_Out = Arm() /* this is the sample-ana axis */ AT (0,0,0) RELATIVE Ana_Cradle ROTATED (0, 42, 0) RELATIVE Ana_Cradle COMPONENT SC4 =Collimator_linear( xmin =-0.06/2, ymin =-0.12/2, xmax = 0.06/2, ymax = 0.12/2, length = 0.24) AT (0, 0, 0.25) RELATIVE Ana_Out /* vertical 3He Detector */ COMPONENT He3H = Monitor( xmin = -0.025400, xmax = 0.025400, ymin = -0.042850, ymax = 0.042850) AT (0, 0, 0.7) RELATIVE Ana_Out /* end of ThALES */ COMPONENT H53_ThALES_Transmit = Arm() AT (0,0,0) RELATIVE H53_ThALES_Monochromator_Cradle COMPONENT H53_ThALES_TransmitXY = Monitor_nD( xwidth=0.06, yheight=0.12, options="x lambda limits=[1 10]", bins=50) AT (0,0,0.2) RELATIVE H53_ThALES_Monochromator_Cradle EXTEND %{ ABSORB; %} /* ================================ H52 ===================================== */ /* feeds: H523:CryoEDM, H521: D16, Super-Adam, H522: Wasp */ COMPONENT H52_start = Arm() AT (0,0,0) RELATIVE H52_origin COMPONENT H52_inpile = COPY(H5_rect)(w1=0.062, m=3, l=4.930-3.155, w2=0.068, h2=0.132) AT (0,0,0) RELATIVE PREVIOUS /* divergent section 60x120 -> w120 x h240 length = 20998-4930=16.068 m R=4000 m */ COMPONENT H52_Common = COPY(H52_inpile)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16, l=16.068/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS /* in the copy, the parameters are passed as is, and the PREVIOUS in the cexp is not updated */ COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT H52_Common_3 = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=COMP_GETPAR(H52_Common_3, w2)) AT (0,0,COMP_GETPAR(H52_Common_3, l)+0.01) RELATIVE PREVIOUS /* Gap 75 mm after BOG (barillet obturateur general): VS */ COMPONENT H52_Obt_Out = COPY(H5)(xwidth=COMP_GETPAR(H52_Common_3, w2), restore_neutron=1) AT (0,0,0.04) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=COMP_GETPAR(H52_Common_3, w2)) AT (0,0,0.01) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(H52_Common_3)( w1=COMP_GETPAR(H52_Common_3, w2), h1=COMP_GETPAR(H52_Common_3, h2), w2=COMP_GETPAR(H52_Common_3, w2)+0.052/16, h2=COMP_GETPAR(H52_Common_3, h2)+0.108/16, l=16.068/16) AT (0,0,COMP_GETPAR(H52_Common_3, l)+gGap+0.075) RELATIVE H52_Common_3 ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H52_Common) = COPY(PREVIOUS)( w1=COMP_GETPAR(PREVIOUS, w2), h1=COMP_GETPAR(PREVIOUS, h2), w2=COMP_GETPAR(PREVIOUS, w2)+0.052/16, h2=COMP_GETPAR(PREVIOUS, h2)+0.108/16) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,16.068/16/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.12,yheight=0.24) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS COMPONENT H52_Common_Out = Monitor_nD(xwidth=0.120, yheight=0.240) AT (0,0,0.01+gGap) RELATIVE PREVIOUS /* distribute the H52 beam along channels 1 2 3 */ COMPONENT H523_origin = Monitor_nD(xwidth=0.06, yheight=0.120, options="dx limits=[-2 2] dy limits=[-2 2] ", bins=25, restore_neutron=1) AT (0.03,-.06,0) RELATIVE H52_Common_Out JUMP H523_start WHEN SCATTERED COMPONENT H521_origin = COPY(H523_origin) AT (-0.03,-.06,0) RELATIVE H52_Common_Out JUMP H521_start WHEN SCATTERED COMPONENT H522_origin = COPY(H53_origin)(xwidth=0.12,yheight=0.12) AT (0,0.06,0) RELATIVE H52_Common_Out JUMP H522_start WHEN SCATTERED /* ================================ H521 ==================================== */ /* bottom right (closer to H5 axis) */ /* feeds: D16 and Super-ADAM */ /* rotated downwards by 0.05 deg for 17 m to get 15 mm separation and avoid */ /* collision with H522 */ COMPONENT H521_start = Arm() AT (0,0,0) RELATIVE H521_origin /* curved m=2/2.5 R=1500 m l=56000-20998=35 m 60x120->60x150 */ COMPONENT H521_Curved = COPY(H5_rect)(w1=0.06, h1=0.12, w2=0.06, h2=0.12+0.03/35, l=35.0/35.0, m=2, mright=2.5) AT (0,0,gGap) RELATIVE PREVIOUS ROTATED (0.05,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS /* downwards */ COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35, l=35.0/35.0) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT H521_Curved_preVTE = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.06,yheight=0.15) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap+0.001) RELATIVE PREVIOUS /* VTE space+OT H53: 108 mm + 2 mm Al window */ COMPONENT H521_VTE = COPY(H5)(xwidth=0.06, yheight=0.15, restore_neutron=1) AT (0,0,0.108/2-0.001) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(zdepth=0.002,xwidth=0.06,yheight=0.15) AT (0,0,0.108/2-0.05) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(H521_Curved_preVTE)( h1=COMP_GETPAR(H521_Curved_preVTE, h2), h2=COMP_GETPAR(H521_Curved_preVTE, h2)+0.03/35, l=35.0/35.0) AT (0,0,0.003+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS) /* back to horizontal */ AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (-0.05,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0.05,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H521_Curved) = COPY(PREVIOUS)( h1=COMP_GETPAR(PREVIOUS, h2), h2=COMP_GETPAR(PREVIOUS, h2)+0.03/35) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35.0/1500*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.06,yheight=0.15) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS COMPONENT H521_Curved_Out = Monitor_nD(xwidth=0.06,yheight=0.15,options="x y") AT (0,0,0.01+gGap) RELATIVE PREVIOUS /* straight m=2 l=4 m 60x150 */ COMPONENT H521_Straight = Guide_gravity(w1=0.06, h1=0.15,m=2, l=4,nelements=4) AT (0,0,gGap) RELATIVE PREVIOUS COMPONENT H521_Straight_Out = Monitor_nD(xwidth=0.06,yheight=0.15,options="x y") AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.06,yheight=0.15) AT (0,0,gGap+0.001) RELATIVE PREVIOUS /* D16 */ COMPONENT H521_D16_Monochromator_Cradle = COPY(H5)(xwidth=0.06, restore_neutron=1) AT (0,0,0.200/2-0.001) RELATIVE PREVIOUS EXTEND %{ flag=0; %} /* D16 Monochromator PG002 w6*h12 7 blades vertically focusing. 'Sample' at 2.8 m */ COMPONENT H521_D16_Monochromator = Monochromator_curved( width = 0.06, height = 0.12, NH=1, NV=7, RV=D16_RMV, RH=0, DM=D16_DM, mosaich = 40, mosaicv = 40, r0 = 1, reflect="HOPG.rfl") AT (0,0,0) RELATIVE H521_D16_Monochromator_Cradle ROTATED (0,D16_A1,0) RELATIVE H521_D16_Monochromator_Cradle EXTEND %{ flag=SCATTERED; %} COMPONENT H521_D16_Monochromator_Jumper = Arm() AT (0,0,0) RELATIVE H521_D16_Monochromator_Cradle JUMP H521_D16_Transmit WHEN !flag COMPONENT H521_D16_Monochromator_Out = Arm() AT (0,0,0) RELATIVE H521_D16_Monochromator_Cradle ROTATED (0,2*D16_A1,0) RELATIVE H521_D16_Monochromator_Cradle COMPONENT H521_D16_Sample_Div = Monitor_nD(xwidth=0.05, yheight=0.05, options="dx limits=[-2 2], dy limits=[-2 2]", bins=50, restore_neutron=1) AT (0,0,D16_L) RELATIVE H521_D16_Monochromator_Out COMPONENT H521_D16_Sample_XY = Monitor_nD(xwidth=0.1, yheight=0.1, options="x y", bins=50, restore_neutron=1) AT (0,0,D16_L) RELATIVE H521_D16_Monochromator_Out COMPONENT H521_D16_Sample_L = Monitor_nD(xwidth=0.1, yheight=0.1, options="lambda limits=[1 10]", bins=50, restore_neutron=1) AT (0,0,D16_L) RELATIVE H521_D16_Monochromator_Out COMPONENT H521_D16_Sample = Isotropic_Sqw(radius=0.005,yheight=0.05, Sqw_coh=D16_sample, Sqw_inc=NULL, d_phi=30) AT (0,0,D16_L) RELATIVE H521_D16_Monochromator_Out EXTEND %{ if (!SCATTERED) ABSORB; %} COMPONENT H521_D16_Detector = Monitor_nD( xwidth=0.32, yheight=0.32, options="x y", bins=50, restore_neutron=1) AT (0,0,1.1) RELATIVE H521_D16_Sample COMPONENT H521_D16_Spectrometer = Monitor_nD(radius=0.5, yheight=0.3, options="angle limits=[-150 150] bins=50, energy limits=[0 30], banana", restore_neutron=1) AT (0,0,0) RELATIVE H521_D16_Sample COMPONENT H521_D16_Diffractometer = Monitor_nD(radius=1, yheight=0.3,options="theta bins=100 limits=[-150 150], y, banana") AT (0,0,0) RELATIVE H521_D16_Sample EXTEND %{ ABSORB; /* neutron ends on the D16 sample position */ %} COMPONENT H521_D16_Transmit = Arm() AT (0,0,0) RELATIVE H521_D16_Monochromator_Cradle COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.06,yheight=0.15) AT (0,0,0.200/2-0.01) RELATIVE H521_D16_Transmit COMPONENT H521_P6_Start = Arm() AT (0,0,0.200/2) RELATIVE H521_D16_Transmit /* Carter Poutre 6 */ COMPONENT H521_P6 = COPY(H53_inpile)(l=19.321, m=3,w1=0.06,h1=0.15,h2=0.12) AT (0,0,0) RELATIVE H521_P6_Start COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.06,yheight=0.12) AT (0,0,19.321+0.001) RELATIVE PREVIOUS /* Gap 524 mm: SuperADAM Monochromator */ /* SuperADAM */ COMPONENT H521_SADAM_Monochromator_Cradle = COPY(H5)(xwidth=0.06, restore_neutron=1) AT (0,0,0.524/2-0.001) RELATIVE PREVIOUS EXTEND %{ flag=0; %} /* D16 Monochromator PG002 w6*h12 7 blades vertically focusing. 'Sample' at 2.8 m */ COMPONENT H521_SADAM_Monochromator = Monochromator_curved( width = 0.06, height = 0.12, NH=1, NV=7, RV=SADAM_RMV, RH=0, DM=SADAM_DM, mosaich = 40, mosaicv = 40, r0 = 1, reflect="HOPG.rfl") AT (0,0,0) RELATIVE H521_SADAM_Monochromator_Cradle ROTATED (0,SADAM_A1,0) RELATIVE H521_SADAM_Monochromator_Cradle EXTEND %{ if (!SCATTERED) ABSORB; %} COMPONENT H521_SADAM_Monochromator_Out = Arm() AT (0,0,0) RELATIVE H521_SADAM_Monochromator_Cradle ROTATED (0,2*SADAM_A1,0) RELATIVE H521_SADAM_Monochromator_Cradle COMPONENT H521_SADAM_Sample_Div = Monitor_nD(xwidth=0.05, yheight=0.05, options="dx limits=[-2 2], dy limits=[-2 2]", bins=50, restore_neutron=1) AT (0,0,SADAM_L) RELATIVE H521_SADAM_Monochromator_Out COMPONENT H521_SADAM_Sample_XY = Monitor_nD(xwidth=0.1, yheight=0.1, options="x y", bins=50, restore_neutron=1) AT (0,0,SADAM_L) RELATIVE H521_SADAM_Monochromator_Out COMPONENT H521_SADAM_Sample_L = Monitor_nD(xwidth=0.1, yheight=0.1, options="lambda limits=[1 10]", bins=50, restore_neutron=1) AT (0,0,SADAM_L) RELATIVE H521_SADAM_Monochromator_Out COMPONENT H521_SADAM_Sample = Isotropic_Sqw(radius=0.005,yheight=0.05, Sqw_coh=SADAM_sample, Sqw_inc=NULL, d_phi=30) AT (0,0,SADAM_L) RELATIVE H521_SADAM_Monochromator_Out EXTEND %{ if (!SCATTERED) ABSORB; %} COMPONENT H521_SADAM_Detector = Monitor_nD( xwidth=0.32, yheight=0.32, options="x y", bins=50, restore_neutron=1) AT (0,0,1.1) RELATIVE H521_SADAM_Sample COMPONENT H521_SADAM_Spectrometer = Monitor_nD(radius=0.5, yheight=0.3, options="angle limits=[-150 150] bins=50, energy limits=[0 30], banana", restore_neutron=1) AT (0,0,0) RELATIVE H521_SADAM_Sample COMPONENT H521_SADAM_Diffractometer = Monitor_nD(radius=1, yheight=0.3,options="theta bins=100 limits=[-150 150], y, banana") AT (0,0,0) RELATIVE H521_SADAM_Sample EXTEND %{ ABSORB; /* neutron ends on the SuperADAM sample position */ %} /* ================================ H522 ==================================== */ /* feeds: WASP */ COMPONENT H522_start = Arm() AT (0,0,0) RELATIVE H522_origin /* curved m=2 R=400 m l=35 m 120x120 */ COMPONENT H522_Curved = COPY(H5_rect)(w1=0.12, h1=0.12,l=35.0/35, m=2) AT (0,0,gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.12,yheight=0.12) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap+0.001) RELATIVE PREVIOUS /* VTE space+OT H53: 108 mm + 2 mm Al window */ COMPONENT H522_VTE = COPY(H5)(xwidth=0.12, restore_neutron=1) AT (0,0,0.108/2-0.001) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(zdepth=0.002,xwidth=0.12,yheight=0.12) AT (0,0,0.108/2-0.05) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(H5_rect)(w1=0.12, h1=0.12,l=35.0/35, m=2) AT (0,0,0.003+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H522_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,35.0/35/4000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.12,yheight=0.12) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap+0.001) RELATIVE PREVIOUS COMPONENT H522_Curved_Out_XY = Monitor_nD(xwidth=0.12,yheight=0.12,options="x y", restore_neutron=1) AT (0,0,0.01+gGap) RELATIVE PREVIOUS COMPONENT H522_Curved_Out_XL = COPY(PREVIOUS)(options="x, lambda limits=[1 10]", bins=50) AT (0,0,0) RELATIVE PREVIOUS /* WASP Velocity Selector: NVS 030 in 60 cm gap */ COMPONENT WASP_Vselector = V_selector(xwidth=0.12, yheight=0.085, zdepth=0.60, radius=0.43/2, alpha=14.2, length=0.5, d=0.0004, nu=WASP_nu, nslit=144) AT (0,0,0.31) RELATIVE H522_Curved_Out_XL ROTATED (0,0,-90) RELATIVE H522_Curved_Out_XL COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.12,yheight=0.12) AT (0,0,0.62) RELATIVE H522_Curved_Out_XL /* straight section l=25.77 m=2, including polariser+Vpol */ COMPONENT H522_Straight = Guide_gravity(w1=0.075,h1=0.12,l=25.77) AT (0,0,0.01) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.075,yheight=0.12) AT (0,0,25.77) RELATIVE PREVIOUS COMPONENT H522_Vsel_L = Monitor_nD(xwidth=0.075, yheight=0.12, options="slit lambda limits=[1 10]", bins=50) AT (0,0,0.01) RELATIVE PREVIOUS /* WASP polariser l=4.8 m */ COMPONENT H522_mirror_centre = Arm() AT (0,0,4.8/2) RELATIVE PREVIOUS EXTEND %{ flag=0; %} /* polarising splitter: an iterating assembly of polarising */ /* The single polariser is enough and simpler to handle. We ignore the surrouding m=2 guide */ COMPONENT WASP_polariser = Pol_mirror(zwidth = 4.8, yheight = 0.12, rUpPar ={1.0, 0.0219, 4.07, 3.2, 0.003}, rDownPar={0.1, 0.0219, 4.07, 3.2, 0.003}, p_reflect=0.95) AT (0, 0, 0) RELATIVE H522_mirror_centre ROTATED (0, -0.89, 0) RELATIVE H522_mirror_centre EXTEND %{ flag=SCATTERED; %} COMPONENT H522_mirror_outdir = Arm() AT (0,0,0) RELATIVE H522_mirror_centre ROTATED (0, -0.89*2, 0) RELATIVE H522_mirror_centre COMPONENT H522_mirror_ReflectXY = Monitor_nD(xwidth=0.075, restore_neutron=1, yheight=0.12, bins=50) WHEN flag==1 AT (0,0,4.8/2) RELATIVE H522_mirror_outdir COMPONENT H522_mirror_TransmitXY = COPY(H522_mirror_ReflectXY) WHEN flag==2 AT (0,0,4.8/2) RELATIVE H522_mirror_centre EXTEND %{ if (flag == 2) ABSORB; %} COMPONENT H522_mirror_outdXY = COPY(H522_mirror_ReflectXY)(options="dx limits=[-2 2], dy limits=[-2 2]") AT (0,0,0) RELATIVE H522_mirror_ReflectXY /* WASP polarising V-mirror */ COMPONENT WASP_Vpolariser = Pol_guide_vmirror(xwidth = 0.075, yheight = 0.12, length = 1.210, rPar ={1, 2.0*0.0219, 4.07, 1, 0.003}, rUpPar ={1, 3.2*0.0219, 4.07, 1, 0.003}, rDownPar={1, 0.4*0.0219, 4.07, 1, 0.003}) AT (0,0,0.01) RELATIVE PREVIOUS /* last straight section 1.76 m */ COMPONENT H522_Straight_Last = Guide_gravity(w1=0.075,h1=0.12,l=1.76) AT (0,0,1.21+0.01) RELATIVE PREVIOUS /* sample at 2.75 m */ COMPONENT H522_WASP_Sample_Div = Monitor_nD(xwidth=0.02, yheight=0.05, options="dx limits=[-2 2], dy limits=[-2 2]", bins=100, restore_neutron=1) AT (0,0,WASP_L+1.76) RELATIVE H522_Straight_Last COMPONENT H522_WASP_Sample_XY = Monitor_nD(xwidth=0.02, yheight=0.05, options="x y", bins=50, restore_neutron=1) AT (0,0,WASP_L+1.76) RELATIVE H522_Straight_Last COMPONENT H522_WASP_Sample_L = Monitor_nD(xwidth=0.02, yheight=0.05, options="wavelength limits=[1 10]", bins=50, restore_neutron=1) AT (0,0,WASP_L+1.76) RELATIVE H522_Straight_Last COMPONENT H522_WASP_Sample = Incoherent(radius=0.02,yheight=0.04, focus_aw=180, focus_ah=30, target_index=+1) AT (0,0,WASP_L+1.76) RELATIVE H522_Straight_Last EXTEND %{ if (!SCATTERED) ABSORB; %} COMPONENT H522_WASP_Detector = Monitor_nD(xwidth=0.32, yheight=0.32, options="x y", restore_neutron=1, bins=32) AT (0,0,4.6) RELATIVE H522_WASP_Sample COMPONENT H522_WASP_Spectrometer = Monitor_nD(radius=3, yheight=0.3, options="angle limits=[-170 170] bins=50, energy limits=[0 5], banana", restore_neutron=1) AT (0,0,0) RELATIVE H522_WASP_Sample COMPONENT H522_WASP_Diffractometer = Monitor_nD(radius=3.5,yheight=0.3, options="theta bins=100 limits=[-170 170], y, banana") AT (0,0,0) RELATIVE H522_WASP_Sample EXTEND %{ ABSORB; /* neutron ends on the SuperADAM sample position */ %} /* ================================ H523 ==================================== */ /* feeds: Cryo-EDM */ COMPONENT H523_start = Arm() AT (0,0,0) RELATIVE H523_origin /* curved m=2/2.5 R=800 m l=50 m 60x120 */ COMPONENT H523_Curved = COPY(H5_rect)(w1=0.06, h1=0.12,l=50.0/25, m=1.2, mright=1.5) AT (0,0,gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.06,yheight=0.12) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap+0.001) RELATIVE PREVIOUS /* VTE space+OT H53: 108 mm + 2 mm Al window */ COMPONENT H523_VTE = COPY(H5)(xwidth=0.06, restore_neutron=1) AT (0,0,0.108/2-0.001) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(zdepth=0.002,xwidth=0.06,yheight=0.12) AT (0,0,0.108/2-0.05) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(H5_rect)(w1=0.06, h1=0.12,l=50.0/25, m=1.2, mright=1.5) AT (0,0,0.003+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H523_Curved) = COPY(PREVIOUS) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS ROTATED (0,50.0/25.0/800*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.06,yheight=0.12) AT (0,0,COMP_GETPAR(PREVIOUS, l)+gGap) RELATIVE PREVIOUS /* Cryo-EDM */ COMPONENT H523_CryoEDM_In = Arm() AT (0,0,CryoEDM_L) RELATIVE PREVIOUS COMPONENT H523_CryoEDM_In_L = Monitor_nD(xwidth=0.06, yheight=0.12, options="lambda limits=[1 10]", bins=100, restore_neutron=1) AT (0,0,0) RELATIVE H523_CryoEDM_In COMPONENT H523_CryoEDM_mirror_center = Arm() AT (0,0,2.45/2) RELATIVE H523_CryoEDM_In COMPONENT H523_CryoEDM_polariser = Pol_mirror(zwidth = 2.45, yheight = 0.055, rUpPar ={1.0, 0.0219, 4.07, 3.2, 0.003}, rDownPar={0.1, 0.0219, 4.07, 3.2, 0.003}, p_reflect=0.95) AT (0, 0, 0) RELATIVE H523_CryoEDM_mirror_center ROTATED (0, 2, 0) RELATIVE H523_CryoEDM_mirror_center EXTEND %{ flag=SCATTERED; %} COMPONENT H523_CryoEDM_mirror_outdir = Arm() AT (0,0,0) RELATIVE H523_CryoEDM_mirror_center ROTATED (0, 4, 0) RELATIVE H523_CryoEDM_mirror_center COMPONENT H523_CryoEDM_mirror_ReflectXY = Monitor_nD(xwidth=0.04, restore_neutron=1, yheight=0.055, bins=50) WHEN flag==1 AT (0,0,2.5/2) RELATIVE H523_CryoEDM_mirror_outdir COMPONENT H523_CryoEDM_mirror_TransmitXY = COPY(H523_CryoEDM_mirror_ReflectXY) WHEN flag==2 AT (0,0,2.5/2) RELATIVE H523_CryoEDM_mirror_center EXTEND %{ if (flag == 2) ABSORB; %} COMPONENT H523_CryoEDM_Sample_L = Monitor_nD(xwidth=0.06, yheight=0.12, options="lambda limits=[1 10]", bins=100, restore_neutron=1) AT (0,0,2.45/2+0.2) RELATIVE H523_CryoEDM_mirror_outdir COMPONENT H523_CryoEDM_Sample_Div = COPY(H523_CryoEDM_Sample_L)(options="dx limits=[-2 2], dy limits=[-2 2]") AT (0,0,0) RELATIVE H523_CryoEDM_Sample_L COMPONENT H523_CryoEDM_Sample_XY = COPY(H523_CryoEDM_Sample_L)(options="x y") AT (0,0,0) RELATIVE H523_CryoEDM_Sample_L EXTEND %{ ABSORB; /* end of H53 */ %} /* ================================ H51 ===================================== */ /* feeds: IN15, D22 */ COMPONENT H51_start = Arm() AT (0,0,0) RELATIVE H51_origin COMPONENT H51_inpile = COPY(H5_rect)(w1=0.04, m=1.2, l=4.930-3.155) AT (0,0,0) RELATIVE PREVIOUS COMPONENT H51_Obt_1 = COPY(PREVIOUS)(l=0.863) AT (0,0,4.930-3.155+0.025) RELATIVE PREVIOUS COMPONENT H51_Obt_2 = COPY(PREVIOUS)(l=1.433) AT (0,0,0.863+0.012) RELATIVE PREVIOUS COMPONENT H51_Obt_3 = COPY(PREVIOUS)(l=0.745) AT (0,0,1.433+0.012) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.04) AT (0,0,0.745+0.01) RELATIVE PREVIOUS /* Gap 75 mm after BOG (barillet obturateur general): VS */ COMPONENT H51_Obt_Out = COPY(H5)(xwidth=0.04, restore_neutron=1) AT (0,0,0.06) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.04) AT (0,0,0) RELATIVE PREVIOUS /* Sections 2,3 in NoteDPT11/H51 - curved R=3000 m. Length=2.5 m down to splitter */ /* Carter VS+Movable+Common */ COMPONENT H51_S2 = COPY(H53_inpile)(w1=0.04, l=2.5/3, m=1.2) AT (0,0,0.05) RELATIVE PREVIOUS COMPONENT COPY(H51_S2) = COPY(H51_S2) AT (0,0,2.5/3) RELATIVE PREVIOUS ROTATED (0,-2.5/3/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H51_S2) = COPY(H51_S2) AT (0,0,2.5/3) RELATIVE PREVIOUS ROTATED (0,-2.5/3/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT HCS_Al_H51 = COPY(HCS_Al)(xwidth=0.04) AT (0,0,2.5/3+0.01) RELATIVE PREVIOUS /* Splitter (polarising "Deviateur") H511/H512 L=1.5 */ /* D22 in transmission, IN15 in reflection: 79.70-76.45 = 3.25 deg angle */ /* bottom part direct to d22 height=55mm */ /* top part part in reflection to in15 height=55mm: transmitted part lost */ /* split the beam to either D22 or IN15 (monitors) */ COMPONENT H51_S2_Out = COPY(H5)(xwidth=0.04, yheight=0.12) AT (0,0,0.05) RELATIVE PREVIOUS COMPONENT H51_split_2 = COPY(H5)(xwidth=0.04, restore_neutron=1, yheight=0.055) AT (0,0.12/2-0.055/2,0) RELATIVE H51_S2_Out JUMP H512_Start WHEN SCATTERED COMPONENT H51_split_1 = COPY(H5)(xwidth=0.04, restore_neutron=1, yheight=0.055) AT (0,-0.12/2+0.055/2,0) RELATIVE H51_S2_Out JUMP H511_Start WHEN SCATTERED /* H511 bottom to IN15 ====================================================== */ COMPONENT H511_Start = Arm() AT (0,0,0) RELATIVE H51_split_1 COMPONENT H511_mirror_centre = Arm() AT (0,0,0.75) RELATIVE H511_Start EXTEND %{ flag=0; %} /* polarising splitter: an iterating assembly of polarising */ /* The single polariser is enough and simpler to handle. We ignore the surrouding m=1.2 guide */ COMPONENT IN15_polariser = Pol_mirror(zwidth = 1.5, yheight = 0.055, rUpPar ={1.0, 0.0219, 4.07, 3.2, 0.003}, rDownPar={0.1, 0.0219, 4.07, 3.2, 0.003}, p_reflect=0.95) AT (0, 0, 0) RELATIVE H511_mirror_centre ROTATED (0, -1.5, 0) RELATIVE H511_mirror_centre EXTEND %{ flag=SCATTERED; %} COMPONENT H511_mirror_outdir = Arm() AT (0,0,0) RELATIVE H511_mirror_centre ROTATED (0, -3, 0) RELATIVE H511_mirror_centre COMPONENT H511_mirror_ReflectXY = Monitor_nD(xwidth=0.04, restore_neutron=1, yheight=0.055, bins=50) WHEN flag==1 AT (0,0,0.755) RELATIVE H511_mirror_outdir COMPONENT H511_mirror_TransmitXY = COPY(H511_mirror_ReflectXY) WHEN flag==2 AT (0,0,0.755) RELATIVE H511_mirror_centre EXTEND %{ if (flag == 2) ABSORB; %} COMPONENT H511_mirror_outdXY = COPY(H511_mirror_ReflectXY)(options="dx limits=[-2 2], dy limits=[-2 2]") AT (0,0,0) RELATIVE H511_mirror_ReflectXY COMPONENT H511_S3 = COPY(H51_S2)(l=3, nelements=2, h1=0.055) AT (0,0,0.755) RELATIVE H511_mirror_outdir COMPONENT H511_S9 = COPY(H51_S2)(l=16.483-0.05, nelements=13, m=1, h1=0.055) AT (0,0,3.005) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.04) AT (0,0,16.483+0.001-0.05) RELATIVE PREVIOUS /* VTE space+OT H511: 108 mm + 2 mm Al window */ COMPONENT H511_VTE = COPY(H5)(xwidth=0.06, restore_neutron=1) AT (0,0,0.108/2-0.001) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(thickness=0.002,xwidth=0.04) AT (0,0,0.108/2-0.05) RELATIVE PREVIOUS COMPONENT H511_S10 = COPY(H51_S2)(l=1, m=1, h1=0.055) AT (0,0,0.05) RELATIVE PREVIOUS /* IN15 Velocity Selector EADS/Astrium NVS 024 in 60 cm gap */ COMPONENT IN15_Vselector = V_selector(xwidth=0.04, yheight=0.055, zdepth=0.30, radius=0.12, alpha=64.4, length=0.25, d=0.0004, nu=IN15_nu, nslit=36) AT (0,0,1+0.3) RELATIVE H511_S10 /* H511/IN15 guide L=9.5 m, m=1 */ COMPONENT H511_BeforeV = COPY(H511_S10)(l=9.5,nelements=6) AT (0,0,0.3) RELATIVE PREVIOUS /* IN15 polarising V-mirror */ COMPONENT IN15_Vpolariser = Pol_guide_vmirror(xwidth = 0.04, yheight = 0.055, length = 0.43, rPar ={1, 1.0*0.0219, 4.07, 1, 0.003}, rUpPar ={1, 3.2*0.0219, 4.07, 1, 0.003}, rDownPar={1, 0.4*0.0219, 4.07, 1, 0.003}) AT (0,0,9.5+0.01) RELATIVE PREVIOUS COMPONENT H511_AfterV = COPY(H511_S10)(l=3, nelements=2) AT (0,0,0.44) RELATIVE PREVIOUS /* sample at 3.7 m */ COMPONENT H511_IN15_Sample_Div = Monitor_nD(xwidth=0.02, yheight=0.05, options="dx limits=[-2 2], dy limits=[-2 2]", bins=100, restore_neutron=1) AT (0,0,IN15_L) RELATIVE H511_AfterV COMPONENT H511_IN15_Sample_XY = Monitor_nD(xwidth=0.02, yheight=0.05, options="x y", bins=50, restore_neutron=1) AT (0,0,IN15_L) RELATIVE H511_AfterV COMPONENT H511_IN15_Sample_L = Monitor_nD(xwidth=0.02, yheight=0.05, options="lambda limits=[1 10]", bins=50, restore_neutron=1) AT (0,0,IN15_L) RELATIVE H511_AfterV COMPONENT H511_IN15_Sample = Incoherent(radius=0.005,yheight=0.05, focus_aw=180, focus_ah=30, target_index=+1) AT (0,0,IN15_L) RELATIVE H511_AfterV EXTEND %{ if (!SCATTERED) ABSORB; %} COMPONENT H511_IN15_Detector = Monitor_nD(xwidth=0.32, yheight=0.32, options="x y", restore_neutron=1, bins=32) AT (0,0,4.6) RELATIVE H511_IN15_Sample COMPONENT H511_IN15_Spectrometer = Monitor_nD(radius=0.5, yheight=0.3, options="angle limits=[-150 150] bins=50, energy limits=[0 5], banana", restore_neutron=1) AT (0,0,0) RELATIVE H511_IN15_Sample COMPONENT H511_IN15_Diffractometer = Monitor_nD(radius=1,yheight=0.3, options="theta bins=100 limits=[-150 150], y, banana") AT (0,0,0) RELATIVE H511_IN15_Sample EXTEND %{ ABSORB; /* neutron ends on the SuperADAM sample position */ %} /* H512 top to D22 ========================================================== */ COMPONENT H512_Start = Arm() AT (0,0,0) RELATIVE H51_split_2 /* Sections 3,4,5,6 in NoteDPT11/H51 - curved R=3000 m. Length=2.5 m down to VTE */ /* Carter C1-C4 */ COMPONENT H512_S36 = COPY(H53_inpile)(w1=0.04, l=20.85/16, h1=0.055,m=1.2) AT (0,0,0.0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H512_S36) = COPY(H512_S36) AT (0,0,20.85/16) RELATIVE PREVIOUS ROTATED (0,-20.85/16/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT HCS_Al_H512 = COPY(HCS_Al)(xwidth=0.04) AT (0,0,20.85/16+0.01) RELATIVE PREVIOUS /* VTE space+OT H512: 108 mm + 2 mm Al window */ COMPONENT H512_VTE = COPY(H5)(xwidth=0.06, restore_neutron=1) AT (0,0,0.108/2-0.001) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(thickness=0.002,xwidth=0.04) AT (0,0,0.108/2-0.05) RELATIVE PREVIOUS COMPONENT H51_S78 = COPY(H53_inpile)(w1=0.04, l=6.0/4,h1=0.055,m=1.2) AT (0,0,0.05) RELATIVE PREVIOUS COMPONENT COPY(H51_S78) = COPY(H51_S78) AT (0,0,6.0/4) RELATIVE PREVIOUS ROTATED (0,-6.0/4/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H51_S78) = COPY(H51_S78) AT (0,0,6.0/4) RELATIVE PREVIOUS ROTATED (0,-6.0/4/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(H51_S78) = COPY(H51_S78) AT (0,0,6.0/4) RELATIVE PREVIOUS ROTATED (0,-6.0/4/3000*RAD2DEG,0) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(thickness=0.002, xwidth=0.04) AT (0,0,6.0/4+0.01) RELATIVE PREVIOUS /* 3 m to V selector, straight */ COMPONENT H512_S11 = COPY(H53_inpile)(w1=0.04, l=3.0, nelements=2, h1=0.055, m=1.2) AT (0,0,0.01) RELATIVE PREVIOUS COMPONENT COPY(HCS_Al) = COPY(HCS_Al)(xwidth=0.04) AT (0,0,3+0.01) RELATIVE PREVIOUS COMPONENT D22_PreV = Monitor_nD(xwidth=0.04, yheight=0.05, options="x y", bins=50, restore_neutron=1) AT (0,0,0.01) RELATIVE PREVIOUS /* D22 Velocity Selector: just like D11 EADS/Astrium NVS 023 in 60 cm gap */ COMPONENT D22_Vselector = V_selector(xwidth=0.04, yheight=0.05, zdepth=0.30, radius=0.12, alpha=48.3, length=0.25, d=0.0004, nu=D22_nu, nslit=72) AT (0,0,3+0.01+0.3+0.01) RELATIVE H512_S11 COMPONENT D22_Collimation = COPY(H53_inpile)(w1=0.04, h1=0.055, l=20-D22_collimation, nelements=ceil(D22_collimation/1.5),m=1.2) AT (0,0,0.3) RELATIVE D22_Vselector COMPONENT D22_Sample_Pos = Arm() AT (0,0,20+0.3) RELATIVE D22_Vselector COMPONENT H51_D22_Sample_Div = Monitor_nD(xwidth=0.02, yheight=0.05, options="dx limits=[-2 2], dy limits=[-2 2]", bins=100, restore_neutron=1) AT (0,0,0) RELATIVE D22_Sample_Pos COMPONENT H51_D22_Sample_XY = Monitor_nD(xwidth=0.02, yheight=0.05, options="x y", bins=50, restore_neutron=1) AT (0,0,0) RELATIVE D22_Sample_Pos COMPONENT H51_D22_Sample_L = Monitor_nD(xwidth=0.02, yheight=0.05, options="lambda limits=[1 10]", bins=50, restore_neutron=1) AT (0,0,0) RELATIVE D22_Sample_Pos COMPONENT H51_D22_Sample = Isotropic_Sqw(radius=0.005,yheight=0.05, Sqw_coh=D22_sample, Sqw_inc=NULL, d_phi=20) AT (0,0,0) RELATIVE D22_Sample_Pos EXTEND %{ if (!SCATTERED) ABSORB; %} COMPONENT D22_Detector = Monitor_nD(xwidth=1, yheight=1, options="x y", bins=128) AT (0,0,D22_collimation) RELATIVE D22_Sample_Pos EXTEND %{ ABSORB; /* H512 end here */ %} COMPONENT D22_Detector_vessel = Shape(radius=1.5, yheight=20) AT (0,0,11) RELATIVE D22_Sample_Pos ROTATED (90,0,0) RELATIVE D22_Sample_Pos /* NOTE: wasp 75*120 V selector /* ========================= end of simulation ============================== */ /* inactive components for rendering purposes */ COMPONENT BOG = Shape(radius=0.735, yheight=0.3) AT (0,0,6.55) RELATIVE HCS COMPONENT H53_ThALES_Drum = Shape(radius=0.45, yheight=0.3) AT (0,0,0) RELATIVE H53_ThALES_Monochromator_Cradle COMPONENT H521_D16_Drum = Shape(radius=0.2, yheight=0.3) AT (0,0,0) RELATIVE H521_D16_Monochromator_Cradle COMPONENT H521_SADAM_Drum = Shape(radius=0.2, yheight=0.3) AT (0,0,0) RELATIVE H521_SADAM_Monochromator_Cradle COMPONENT Traversee = Shape(xwidth=2, yheight=0.5) AT (0,0,31.6) RELATIVE HCS /* core at 0.7 m, 44.8 deg */ COMPONENT Core_Centre = Arm() AT (0,0,0) RELATIVE HCS ROTATED (0,44.8+180,0) RELATIVE HCS /* core radius=0.3 m */ COMPONENT Core = Shape (radius=0.3, yheight=0.7) AT (0,0,0.7) RELATIVE Core_Centre ROTATED (0,0,0) RELATIVE Origin /* D2O vessel radius=1.25 m */ COMPONENT D2O_Vessel = Shape (radius=1.25, yheight=1) AT (0,0,0) RELATIVE Core /* H2O vessel radius=1.25 m */ COMPONENT H2O_Vessel = Shape (radius=4.9, yheight=1.2) AT (0,0,0) RELATIVE Core END