/******************************************************************************* * McStas instrument definition URL=http://mcstas.risoe.dk * * Instrument: ILL_H15_D11 * * %Identification * Written by: FARHI Emmanuel (farhi@ill.fr) * Date: 2012 * Origin:ILL * Release: McStas 2.5 * Version: $Revision: 1.0 $ * %INSTRUMENT_SITE: ILL * * D11 at the ILL: Lowest momentum transfer, lowest background small-angle neutron scattering instrument * * %Description * * D11 - LOWEST MOMENTUM TRANSFER and LOWEST BACKGROUND SANS * * D11 is the archetype of a long, pinhole geometry instrument for small angle * neutron scattering (SANS), designed for the study of large scale structures in * soft matter systems, chemistry, biology, solid state physics and materials science. * * This model allows to specify the collimation length, i.e. the free path between * the end of the guide and the sample location. When Lc=1.5 m, the intensity is * maximal, but Q-resolution (divergence) is degraded. The Lc=40.5 m configuration * is the low-Q one. * The collimation can also be given as predefined configuration 'iLc' starting from 0 * * iLc ={ 1.5, 2.5, 4, 5.5, 8, 10.5, 13.5, 16.5, 20.5, 28, 34, 40.5} * * so that e.g. iLc=5 corresponds with a 10.5 m collimation * * The quality of the collimation section can be specified among a set of glass * pre-set specifications. The available configurations are: * * | Borkron_1972 Borkron_2003 Borofloat_2001 Borofloat_2003 * Chamfers [mm] | 0.2 0.2 0.8 0.2 * Wavyness [rad]| 2.5e-5 1e-4 8e-4 2e-4 * * %Example: Lambda=4.51 Detector: SampleC_I=2.1e7 * * %Parameters * Lambda: [Angs] neutron wavelength reaching D11 sample * Config: [string] guide configuration. Must be one of Borkron_1972 Borkron_2003 Borofloat_2001 Borofloat_2003 * Lc: [m] guide collimation, i.e. free path between end of guide and sample location. Must be one of 40.5, 34, 28, 20.5, 16.5 13.5 10.5, 8, 5.5, 4, 2.5, 1.5 * iLc: [index] guide collimation predefined length * Chamfers: [m] chamfers width (input/output movable guide section+longitudinal top/bottom sides) * Waviness: [rad] movable guide Waviness * * %End *******************************************************************************/ DEFINE INSTRUMENT ILL_H15_D11(Lambda=4.51, string Config="Borkron_1972", Lc=0, iLc=5, Chamfers=0.0002, Waviness=2.5e-5) /* H15_D11@ILL simulation for McStas. E. Farhi August 2003 Lambda : wavelength [Angs] Lc : collimation length [m] among 40.5, 34, 28, 20.5, 16.5 13.5 10.5, 8, 5.5, 4, 2.5 1.5 Config : Borkron_1972 Borkron_2003 Borofloat_2001 Borofloat_2003 0.2 0.2 0.8 0.2 Chamfers [mm] 2.5e-5 1e-4 8e-4 2e-4 Waviness [rad] Chamfers: chamfers width (input/output movable guide section+longitudinal top/bottom sides) [m] Waviness: movable guide Waviness [rad] setenv MCSTAS_FORMAT "Matlab binary" setenv MCSTAS_CFLAGS "" mcrun h15_d11.instr --ncount=1e6 --dir=t1_5 Lambda=4.51 Config=Borkron_1972 iLc=5 Lc=0 Chamfers=0.0002 Waviness=2.5e-5 */ DECLARE %{ #define VERSION "H15_D11@ILL Farhi 10/08/2003" /* all variables defined here can be used, and modified, in the whole instrument */ /* variables used in instrument definition and for user computations */ /* Source specifications */ double sLambda = 4; double sDeltaLambdaPercent= 0.08; /* Lambda spread in % (full width) */ double sDeltaLambda; double sT1=216.8,sI1=1.24e+13; /* VCS */ double sT2=33.9, sI2=1.02e+13; double sT3=16.7, sI3=3.0423e+12; double sH = 0.22; double sW = 0.10; /* Guide geometry specifications, default values */ double gInputWidth = 0.038; double gObtWidth = 0.032; double gGuideHeight1 = 0.2; double gGuideHeight2 = 0.05; double gSourceInput = 2.33; double gGuideWidth = 0.03; double gLength1 = 25.5; /* 25 pieces */ double gCurvatureRadius = 2700.0; double gLength2 = 22.284; /* 22 pieces */ double gLength3 = 5.4; /* 5 pieces */ double gToSampleLength = 1.5; double gElementGap = 0.004; double gElmtLength1; double gElmtLength2; double gElmtLength3; double gElmtRot1; /* Guide coatings specifications */ double mInPile = 1; double mGuide1 = 1; double mGuide2 = 0.65; /* 0.953075 */ double mAlpha = 4.07; /* 0 */ double mQc = 0.0218; /* 0.0175455 */ double mR0 = 0.995; /* 0.953075 */ double mW = 1.0/300.0; /* 0.00299565 */ double Vrpm; double Config_l[20]; double Config_n[20]; double Config_g[20]; /* geometry configurations */ double Borkron_1972_l[] = {20, 0,0,0, 0,4,3,3,2.5,0, 2.5,0, 1.5,0,1.5,0,0}; /* length 17 */ double Borkron_2003_l[] = {6.5,0,6,7.5,0,4,3,3,2.5,0, 2.5,0, 1.5,0,1.5,0,1}; double Borofloat_2001_l[]={0.5,6,6,0.5,7,4,3,3,2, 0.5,2, 0.5,0.5,1,0.5,1,1}; double Borkron_1972_n[] = {40,0,0, 0, 0,8,6,6,5,0,5,0,3,0,3,0,0}; double Borkron_2003_n[] = {13,0,12,15,0,8,6,6,5,0,5,0,3,0,3,0,2}; double Borofloat_2001_n[]={1, 6,6, 1, 7,4,3,3,2,1,2,1,1,1,1,1,1}; double Borkron_1972_g[] = {0 ,0 }; double Borkron_2003_g[] = {0.002,0.002}; double Borofloat_2001_g[]={0.002,0.002}; double Length_Lc[] ={1.5,2.5,4,5.5,8,10.5,13.5,16.5,20.5,28,34,40.5}; double microGap = 0.0001; char xlmonopts[128]; %} /* end of DECLARE */ INITIALIZE %{ double PhiC, GammaC, ThetaC, LambdaC, LStar, Zc; double Ki; int i=0; double sum_length=0; double n_elements=0; /* curved section total 25.5 m in 25 elements */ gElmtLength1 = fabs(gLength1)/25; gElmtRot1 = gElmtLength1/gCurvatureRadius*180/PI; gElmtLength2 = fabs(gLength2)/22; gElmtLength3 = fabs(gLength3)/5; /* INFO ******************************************************************* */ printf("Simulation : %s on %s [pid %i]\n",VERSION,getenv("HOSTNAME"), getpid()); /* Source and guide parameters */ sLambda = Lambda; sDeltaLambda = sDeltaLambdaPercent * sLambda / 2; /* First guide Phase-Space parameters */ PhiC = sqrt(2*fabs(gGuideWidth/gCurvatureRadius)); GammaC = mGuide1*mQc/4/PI; ThetaC = sLambda*GammaC; LambdaC = PhiC/GammaC; LStar = sqrt(8*fabs(gGuideWidth*gCurvatureRadius)); Zc = gGuideWidth*(0.5-ThetaC*ThetaC/PhiC/PhiC); printf("Guide: L =%.4f m, rho=%.4f m, m=%.2f\n", gLength1, gCurvatureRadius, mGuide1); printf("Guide: L*=%.2f m, Phi_c = %.2f deg, Gamma_c = %.2f deg/Angs\n", LStar, PhiC*RAD2DEG, GammaC*RAD2DEG); printf("Guide: Zc = %.4f m, Lambda_c = %.2f Angs\n", Zc, LambdaC); if (sLambda < LambdaC) printf("Guide: Garland reflections only\n"); else printf("Guide: Garland and Zig-Zag reflections\n"); if (gLength1 < LStar) printf("Guide: Guide is too short\n"); if (Lambda < 0.12) exit(printf("** D11 Error: Unreachable wavelength for D11: %.4g\n", Lambda)); Ki = 2*PI/Lambda; Vrpm = 124270/(Lambda-0.11832); if (Vrpm < 3100 || Vrpm > 28300) exit(printf("D11 Error: Unreachable V_selector speed for D11: %.4g\n", Vrpm)); if (Vrpm > 9000 && Vrpm < 11900) printf("D11 Warning: Resonance V_selector speed for D11: %.4g\n", Vrpm); printf("D11: Ki = %.4g [Angs-1], Lambda = %.4g [Angs], V_Dolores=%.4g [rpm]\n", Ki, Lambda, Vrpm); if (!Config) exit(printf("** D11 Error: No movable guide configuration selected for D11\n")); if (strstr(Config, "Borkron_1972")) { for (i=0; i<17; i++) { Config_l[i] = Borkron_1972_l[i]; Config_n[i] = Borkron_1972_n[i]; if (i<2) Config_g[i] = Borkron_1972_g[i]; } } if (strstr(Config, "Borkron_2003")) { for (i=0; i<17; i++) { Config_l[i] = Borkron_2003_l[i]; Config_n[i] = Borkron_2003_n[i]; if (i<2) Config_g[i] = Borkron_2003_g[i]; } } if (strstr(Config, "Borofloat_2001")) { for (i=0; i<17; i++) { Config_l[i] = Borofloat_2001_l[i]; Config_n[i] = Borofloat_2001_n[i]; if (i<2) Config_g[i] = Borofloat_2001_g[i]; } } if (strstr(Config, "Borofloat_2003")) { for (i=0; i<17; i++) { Config_l[i] = Borofloat_2001_l[i]; Config_n[i] = Borofloat_2001_n[i]; if (i<2) Config_g[i] = Borofloat_2001_g[i]; } } if (!i) exit(printf("** D11 Error: No movable guide configuration selected for D11\n")); if (strstr(Config,"nogap")) for (i=0; i<2; Config_g[i++]=0); if (Lc == 0) if (iLc >= 0 && iLc < 12) Lc = Length_Lc[(int)iLc]; else exit(printf("** D11 Error: invalid collimator length specification for D11\n")); for (i=0; i< 17; i++) { if (40.5 - Lc <= sum_length) Config_n[i]=0; else { sum_length += Config_l[i]; n_elements += Config_n[i]; } } printf("D11: Coll. L_c=%8.3g [m], Configuration='%s', \n", Lc, Config); printf("D11: Length L_g=%8.3g [m], n_elements=%8.2g\n", sum_length, n_elements); printf("D11: Chamfers =%8.3g [m], Waviness =%8.3g [rad] (%.3g [deg])\n", Chamfers, Waviness, Waviness*RAD2DEG); Waviness *= RAD2DEG; sprintf(xlmonopts,"x bins=100 lambda limits=[%g %g] bins=100",sLambda-sDeltaLambda,sLambda+sDeltaLambda); %} /* end of INITIALIZE */ TRACE /* Source description ******************************************************* */ /* a flat rectangular source with Maxwellian distribution */ /* Reactor source 22x16 cm */ COMPONENT PG = Progress_bar() AT (0,0,0) ABSOLUTE COMPONENT Source = Source_gen( yheight = sH, xwidth = sW, zdepth = 0.1, dist = gSourceInput, focus_xw = gInputWidth, focus_yh = gGuideHeight1, Lmin = sLambda-sDeltaLambda, Lmax = sLambda+sDeltaLambda, T1 = sT1, I1 = sI1, T2 = sT2, I2 = sI2, T3 = sT3, I3 = sI3, verbose = 1 ) AT (0,0,0) ABSOLUTE COMPONENT AlWindow1 = Al_window(thickness=0.002) AT (0,0,gSourceInput-0.01) RELATIVE Source COMPONENT PinkCarter = Guide_gravity( w1=gInputWidth, h1=gGuideHeight1, w2=gObtWidth, h2=gGuideHeight1, l=3.170, R0=mR0, Qc=mQc, alpha=mAlpha, m=mInPile, W=mW) AT (0,0,gSourceInput) RELATIVE Source /* 2.33 from Source */ COMPONENT AlWindow2 = Al_window(thickness=0.002) AT (0,0,gSourceInput+3.522+0.001) RELATIVE Source /* Followed by a 228 mm guide (for H1) in Obturator, after at 10 mm gap, ends at 5.952 m from source after 80 mm gap */ COMPONENT LeadShutter = Guide_gravity( w1=gObtWidth, h1=gGuideHeight1, w2=gObtWidth, h2=gGuideHeight1, l=0.228, R0=mR0, Qc=mQc, alpha=mAlpha, m=mInPile, W=mW) AT (0,0,gSourceInput+3.522+0.01) RELATIVE Source COMPONENT AlWindow3 = Al_window(thickness=0.002) AT (0,0,gSourceInput+3.522+0.01+0.228+0.005) RELATIVE Source /*-------------------------*/ /* Curved Guide ("MAN") */ /*-------------------------*/ /* curvature radius is 2700 m , start at XR=5.913 m from source */ COMPONENT CurvedGuideStart = Arm() AT (0,0,gSourceInput+3.522+0.01+0.228+0.08) RELATIVE Source /* curved part 1, 25 elements : should be 25.5 m long, rho=2700 */ COMPONENT cg1 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0,0) RELATIVE CurvedGuideStart ROTATED (0, gElmtRot1,0) RELATIVE CurvedGuideStart COMPONENT cg2 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg1 ROTATED (0, gElmtRot1,0) RELATIVE cg1 COMPONENT cg3 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg2 ROTATED (0, gElmtRot1,0) RELATIVE cg2 COMPONENT cg4 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg3 ROTATED (0, gElmtRot1,0) RELATIVE cg3 COMPONENT cg5 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg4 ROTATED (0, gElmtRot1,0) RELATIVE cg4 COMPONENT cg6 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg5 ROTATED (0, gElmtRot1,0) RELATIVE cg5 COMPONENT cg7 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg6 ROTATED (0, gElmtRot1,0) RELATIVE cg6 COMPONENT cg8 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg7 ROTATED (0, gElmtRot1,0) RELATIVE cg7 COMPONENT cg9 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg8 ROTATED (0, gElmtRot1,0) RELATIVE cg8 COMPONENT cg10 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg9 ROTATED (0, gElmtRot1,0) RELATIVE cg9 COMPONENT cg11 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg10 ROTATED (0, gElmtRot1,0) RELATIVE cg10 COMPONENT cg12 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg11 ROTATED (0, gElmtRot1,0) RELATIVE cg11 COMPONENT cg13 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg12 ROTATED (0, gElmtRot1,0) RELATIVE cg12 COMPONENT cg14 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg13 ROTATED (0, gElmtRot1,0) RELATIVE cg13 COMPONENT cg15 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg14 ROTATED (0, gElmtRot1,0) RELATIVE cg14 COMPONENT cg16 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg15 ROTATED (0, gElmtRot1,0) RELATIVE cg15 COMPONENT cg17 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg16 ROTATED (0, gElmtRot1,0) RELATIVE cg16 COMPONENT cg18 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg17 ROTATED (0, gElmtRot1,0) RELATIVE cg17 COMPONENT cg19 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg18 ROTATED (0, gElmtRot1,0) RELATIVE cg18 COMPONENT cg20 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg19 ROTATED (0, gElmtRot1,0) RELATIVE cg19 COMPONENT cg21 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg20 ROTATED (0, gElmtRot1,0) RELATIVE cg20 COMPONENT cg22 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg21 ROTATED (0, gElmtRot1,0) RELATIVE cg21 COMPONENT cg23 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg22 ROTATED (0, gElmtRot1,0) RELATIVE cg22 COMPONENT cg24 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg23 ROTATED (0, gElmtRot1,0) RELATIVE cg23 COMPONENT cg25 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength1-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength1) RELATIVE cg24 ROTATED (0, gElmtRot1,0) RELATIVE cg24 /* gap V.T.E. 0.260 at 27408 */ COMPONENT AlWindow4 = Al_window(thickness=0.002) AT (0,0,gElmtLength1+0.001) RELATIVE cg25 COMPONENT PSD_VTE = Monitor_nD(xwidth=gGuideWidth, yheight=gGuideHeight1, options="x y") AT (0,0,gElmtLength1+0.13) RELATIVE cg25 COMPONENT AlWindow5 = Al_window(thickness=0.002) AT (0,0,gElmtLength1+0.26-0.003) RELATIVE cg25 COMPONENT VTEtoIN6GuideStart = Arm() AT (0,0,gElmtLength1+0.26) RELATIVE cg25 /* Ni guide -> IN6: 22.284 m long at 27.588, 22 elements */ COMPONENT sg1 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0,0) RELATIVE VTEtoIN6GuideStart COMPONENT sg2 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg1 COMPONENT sg3 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg2 COMPONENT sg4 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg3 COMPONENT sg5 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg4 COMPONENT sg6 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg5 COMPONENT sg7 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg6 COMPONENT sg8 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg7 COMPONENT sg9 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg8 COMPONENT sg10 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg9 COMPONENT sg11 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg10 COMPONENT sg12 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg11 COMPONENT sg13 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg12 COMPONENT sg14 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg13 COMPONENT sg15 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg14 COMPONENT sg16 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg15 COMPONENT sg17 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg16 COMPONENT sg18 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg17 COMPONENT sg19 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg18 COMPONENT sg20 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg19 COMPONENT sg21 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg20 COMPONENT sg22 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength2-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength2) RELATIVE sg21 /* gap 0.3 m OT, OS IN6 H15 */ COMPONENT AlWindow6 = Al_window(thickness=0.002) AT (0,0,gElmtLength2+0.001) RELATIVE sg22 COMPONENT PSD_IN6 = Monitor_nD(xwidth=gGuideWidth, yheight=gGuideHeight1, options="x y") AT (0,0,gElmtLength2+0.15) RELATIVE sg22 COMPONENT AlWindow7 = Al_window(thickness=0.002) AT (0,0,gElmtLength2+0.290) RELATIVE sg22 COMPONENT IN6toD7GuideStart = Arm() AT (0,0,gElmtLength2+0.30) RELATIVE sg22 /* Ni guide -> D7 Carter Man 2 L=5.4 */ COMPONENT sg23 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength3-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, 0) RELATIVE IN6toD7GuideStart COMPONENT sg24 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength3-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength3) RELATIVE sg23 COMPONENT sg25 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength3-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength3) RELATIVE sg24 COMPONENT sg26 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength3-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength3) RELATIVE sg25 COMPONENT sg27 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight1,w2=gGuideWidth,h2=gGuideHeight1, l=(gElmtLength3-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide1,W=mW) AT (0,0, gElmtLength3) RELATIVE sg26 /* gap 0.3 m OS D7/D11 at 55572 */ COMPONENT AlWindow8 = Al_window(thickness=0.002) AT (0,0,gElmtLength3+0.001) RELATIVE sg27 COMPONENT PSD_D7 = Monitor_nD(xwidth=gGuideWidth, yheight=gGuideHeight1, options="x y") AT (0,0,gElmtLength3+0.15) RELATIVE sg27 COMPONENT AlWindow9 = Al_window(thickness=0.002) AT (0,0,gElmtLength3+0.290) RELATIVE sg27 COMPONENT D7toD11GuideStart = Arm() AT (0,-.065,gElmtLength3+0.30) RELATIVE sg27 /* Glass guide SPRI 1.25 m h=0.05 AT (0,-0.065,0.3) */ COMPONENT sg28 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2,w2=gGuideWidth,h2=gGuideHeight2, l=(1.25-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide2,W=mW) AT (0,0, 0) RELATIVE D7toD11GuideStart /* Glass guide 0.68 h=0.05 */ COMPONENT sg29 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2,w2=gGuideWidth,h2=gGuideHeight2, l=(0.5-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide2,W=mW) AT (0,0, 1.25) RELATIVE sg28 /* Velocity selector 0.3. Path in atm is 7 cm */ COMPONENT AlWindow10 = Al_window(thickness=0.004) AT (0,0,0.5+0.001) RELATIVE sg29 COMPONENT Mon_D11_In = Monitor_nD(xwidth=gGuideWidth, yheight=gGuideHeight2, options=xlmonopts) AT (0,0,0.5+0.01) RELATIVE sg29 COMPONENT Dolores = V_selector(xwidth=0.03, yheight=0.05, zdepth=0.30, radius=0.12, alpha=48.298, length=0.250, d=0.0004, nu=Vrpm/60, nslit=72) AT (0,0, 0.5+0.025+0.30/2) RELATIVE sg29 COMPONENT AlWindow11 = Al_window(thickness=0.004) AT (0,0,0.15+0.01) RELATIVE Dolores /* Glass guide 0.50 h=0.05 */ COMPONENT sg30 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2,w2=gGuideWidth,h2=gGuideHeight2, l=(0.5-gElementGap), R0=mR0,Qc=mQc,alpha=mAlpha,m=mGuide2,W=mW) AT (0,0, 0.15+0.02) RELATIVE Dolores /* Gap 16 cm, start of movable guide */ COMPONENT AlWindow12 = Al_window(thickness=0.002) AT (0,0,0.5+0.001) RELATIVE sg30 COMPONENT Mon_D11_Out = Monitor_nD(xwidth=gGuideWidth, yheight=gGuideHeight2, options=xlmonopts) AT (0,0,0.5+0.01) RELATIVE sg30 COMPONENT AlWindow13 = Al_window(thickness=0.002) AT (0,0,0.5+0.03-0.003) RELATIVE sg30 COMPONENT MovableGuideStart = Arm() AT (0,0, 0.66) RELATIVE sg30 /* D11 Movable guide start */ /* Configurations [c] Borkron 1972 Borkron 2003 Borofloat 2001 Borofloat 2003 0.2 0.2 0.8 0.2 Chamfers [mm] mg 2.5e-5 1e-4 8e-4 2e-4 Waviness [rad] 0 40*500 13*500 1*500 idem Borofloat 2001 1 6*1000 Gap 2 Gap 2 (34 m) 2 12*500 6*1000 Gap 2 Gap 2 (28 m) 3 15*500 1*500 4 7*1000 Gap 17 Gap 17 Gap 17 (20.5 m *) 5 8*500 8*500 4*1000 Gap 2 Gap 2 Gap 2 (16.5 m) 6 6*500 6*500 3*1000 Gap 2 Gap 2 Gap 2 (13.5 m) 7 6*500 6*500 3*1000 Gap 17 Gap 17 Gap 17 (10.5 m *) 8 5*500 5*500 2*1000 9 1*500 Gap 2 Gap 2 Gap 2 (8 m) 10 5*500 5*500 2*1000 11 1*500 Gap 17 Gap 17 Gap 17 (5.5 m *) 12 3*500 3*500 1*500 13 1*1000 Gap 2 Gap 2 Gap 2 (4 m) 14 3*500 3*500 1*500 15 1*1000 Gap 2 Gap 2 (2.5 m *) 16 void 1 m 2*500 1*1000 */ COMPONENT mg0 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[0]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=Config_n[0], wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,0) RELATIVE MovableGuideStart COMPONENT mg1 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[1]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[1]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[0])) RELATIVE mg0 COMPONENT mg2 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[2]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[2]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[1]+Config_g[0])) RELATIVE mg1 COMPONENT mg3 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[3]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[3]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[2]+Config_g[1])) RELATIVE mg2 COMPONENT mg4 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[4]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[4]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[3])) RELATIVE mg3 /* Gap 17 mm at 20.5 m collimation */ COMPONENT mg5 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[5]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[5]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[4]+0.017)) RELATIVE mg4 COMPONENT mg6 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[6]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[6]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[5]+0.002)) RELATIVE mg5 COMPONENT mg7 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[7]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[7]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[6]+0.002)) RELATIVE mg6 COMPONENT mg8 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[8]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[8]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[7]+0.017)) RELATIVE mg7 COMPONENT mg9 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[9]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[9]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[8])) RELATIVE mg8 COMPONENT mg10 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[10]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[10]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[9]+0.002)) RELATIVE mg9 COMPONENT mg11 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[11]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[11]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[10])) RELATIVE mg10 COMPONENT mg12 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[12]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[12]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[11]+0.017)) RELATIVE mg11 COMPONENT mg13 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[13]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[13]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[12])) RELATIVE mg12 COMPONENT mg14 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[14]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[14]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[13]+0.002)) RELATIVE mg13 COMPONENT mg15 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[15]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[15]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[14])) RELATIVE mg14 COMPONENT mg16 = Guide_gravity(w1=gGuideWidth,h1=gGuideHeight2, w2=gGuideWidth,h2=gGuideHeight2, l=(Config_l[16]), chamfers_tb=Chamfers, chamfers_z=Chamfers, nelements=(Config_n[16]), wavy=Waviness, m=mGuide2, W=mW, alpha=mAlpha) AT (0,0,(microGap+Config_l[15]+0.002)) RELATIVE mg15 /* Sample is 1.5 m ahead */ COMPONENT SampleF = Monitor_nD(xwidth=gGuideWidth, yheight=gGuideHeight2, options="x y, per cm2") AT (0,0,(Config_l[15]+2.5)) RELATIVE mg15 COMPONENT SampleC = Monitor_nD(xwidth=gGuideWidth, yheight=gGuideHeight2, options="x y, capture, per cm2") AT (0,0,0.001) RELATIVE SampleF END