Skip to content
Snippets Groups Projects
Commit dc7b32f2 authored by Frost, Matthew's avatar Frost, Matthew
Browse files

Added IRP PostPPU Source, readied for comparison to moderator measurments.

parent b84e021c
No related branches found
No related tags found
No related merge requests found
......@@ -127,14 +127,15 @@ COMPONENT init = Union_init()
E0=5) //Because sometimes you just need a simple source...
AT (0,0,0) RELATIVE Origin */
/*
COMPONENT Moderator = SNS_source_xy_McStas3(S_filename="ETD_BL10.dat", //Added on 2023jul28 after BOAT comparison
cview_x=0.0, cview_y=0.0,
dview_x=0.1, dview_y=0.12,
dist=2.5,
xw=0.06, yh=0.06,
xfoc = -0.58865,
Emin=3.0, Emax=400.0,
Tmin=0.0, Tmax=1000.0,
Emin=0.07, Emax=100000,
//Tmin=0.0, Tmax=10.0,
sample_E=2,
sample_t=0,
sample_xy=1,
......@@ -142,12 +143,17 @@ COMPONENT Moderator = SNS_source_xy_McStas3(S_filename="ETD_BL10.dat", //Ad
debug=0)
AT (0.006,0,0) RELATIVE Origin
ROTATED (0,13.25,0) RELATIVE Origin
*/
COMPONENT Moderator = SNS_source(filename="a1G3AD-11-f5.dat",dist=2.5,focus_xw=0.06, focus_yh=0.06,
Emin=0.07, Emax=100000)
AT (0.006,0,0) RELATIVE Origin
EXTEND
%{
p*=59.9; //Per Pulse scaling to time averaged per second
p*=59.9/2.0; //Per Pulse scaling to time averaged per MW-second
%}
//FIXED APERTURE
//PER DRAWING VENUS1020M8U8710A120
COMPONENT CV_A1 = Slit(radius=0.028067)
......@@ -201,21 +207,21 @@ COMPONENT Coll1US = Slit(xwidth=0.105, yheight=0.105)//radius=1.0)
COMPONENT Coll1US20cmOFF = Guide_anyshape(geometry="collimatorData/venus1020m8u8710a410Surface.off", R0=0.6, m=0.5, center=0, transmit=0)
WHEN (CollSetting==1)
AT ( 0.0, 0.0, 7.430287) RELATIVE Origin
EXTEND
%{
if (SCATTERED>0) BNCollScatter++;
%}
//EXTEND
//%{
// if (SCATTERED>0) BNCollScatter++;
//%}
//Collimator 1 4cm FOV- Upstream Aperture
//PER DRAWING VENUS1020M8U8710A411
COMPONENT Coll1US04cmOFF = Guide_anyshape(geometry="collimatorData/venus1020m8u8710a411Surface.off", R0=0.6, m=0.5, center=0, transmit=0)
COMPONENT Coll1US04cmOFF = Guide_anyshape(geometry="collimatorData/venus1020m8u8710a411Surface.off", R0=0.6, m=0, center=0, transmit=0)
WHEN (CollSetting==2)
AT ( 0.0, 0.0, 7.430287) RELATIVE Origin
EXTEND
%{
if (SCATTERED>0) BNCollScatter++;
%}
//EXTEND
//%{
// if (SCATTERED>0) BNCollScatter++;
//%}
//%include "Coll1AUS.code" //Replaced by COMPONENT Coll1US20cmOFF
......@@ -227,24 +233,24 @@ COMPONENT Coll1DS = Slit(xwidth=0.105, yheight=0.105)//radius=1.0)
//Collimator 1 20cm FOV- Downstream Aperture
//PER DRAWING VENUS1020M8U8710A412
COMPONENT Coll1DS20cmOFF = Guide_anyshape(geometry="collimatorData/venus1020m8u8710a412Surface.off", R0=0.6, m=0.5, center=0, transmit=0)
COMPONENT Coll1DS20cmOFF = Guide_anyshape(geometry="collimatorData/venus1020m8u8710a412Surface.off", R0=0.6, m=0, center=0, transmit=0)
WHEN (CollSetting==1)
AT ( 0.0, 0.0, 8.353679) RELATIVE Origin
EXTEND
%{
if (SCATTERED>0) BNCollScatter++;
%}
//EXTEND
//%{
// if (SCATTERED>0) BNCollScatter++;
//%}
//Collimator 1 4cm FOV- Downstream Aperture
//PER DRAWING VENUS1020M8U8710A413
COMPONENT Coll1DS04cmOFF = Guide_anyshape(geometry="collimatorData/venus1020m8u8710a413Surface.off", R0=0.6, m=0.5, center=0, transmit=0)
COMPONENT Coll1DS04cmOFF = Guide_anyshape(geometry="collimatorData/venus1020m8u8710a413Surface.off", R0=0.6, m=0, center=0, transmit=0)
WHEN (CollSetting==2)
AT ( 0.0, 0.0, 8.353679) RELATIVE Origin
EXTEND
%{
if (SCATTERED>0) BNCollScatter++;
%}
//EXTEND
//%{
// if (SCATTERED>0) BNCollScatter++;
//%}
//%include "Coll1ADS.code" //Replaced by COMPONENT Coll1DS20cmOFF
......@@ -346,7 +352,7 @@ COMPONENT CaveFTScraper04 = Slit(xwidth = 0.193015, yheight = 0.193015)
//VARIABLE Sample Slit Package - PER DRAWING VENUS1020M8U8710A600
COMPONENT VariableSampleSlitCSYS = Arm()
AT (0, 0, 24.25) RELATIVE Origin
AT (0, 0, 23.75) RELATIVE Origin
//EACH SLIT BLADE IS 1" THICK SINTERED B4C NATURAL ENRICHMENT
COMPONENT SlitBladeTOPUS = Beamstop(xwidth=8.25*0.0254, yheight=8.25*0.0254)
......@@ -373,43 +379,43 @@ COMPONENT SlitBladeLEFTUS = Beamstop(xwidth=8.25*0.0254, yheight=8.25*0.0254)
COMPONENT SlitBladeLEFTDS = Beamstop(xwidth=8.25*0.0254, yheight=8.25*0.0254)
AT (0.2, 0, (1.636+1.0)*0.0254) RELATIVE VariableSampleSlitCSYS
//COMPONENT VariableSampleSlit = Slit(xwidth=SampleSlitWidth, yheight=SampleSlitHeight)
// AT (0,0,24.250000) RELATIVE Origin
COMPONENT VariableSampleSlit = Slit(xwidth=SampleSlitWidth, yheight=SampleSlitHeight)
AT (0,0,24.250000) RELATIVE Origin
COMPONENT SamplePosition = Arm()
AT (0,0,24.50000) ABSOLUTE
//%include "VENUS.Sample"
COMPONENT Sample = PowderN(reflections="Fe_Gamma.laz", zdepth=0.005, yheight=0.02, xwidth=0.02, p_transmit=1.0)
WHEN (SAMPLE)
AT (0,0,0) RELATIVE SamplePosition
//COMPONENT Sample = PowderN(reflections="Fe_Gamma.laz", zdepth=0.005, yheight=0.02, xwidth=0.02, p_transmit=1.0)
//WHEN (SAMPLE)
//AT (0,0,0) RELATIVE SamplePosition
COMPONENT PhosScreenWeighting = Arm()
WHEN (ScintScreenWeighting==1)
AT (0,0,24.999) RELATIVE Origin
EXTEND
%{
lambdaAngstroms = 3956.0/sqrt(vx*vx+vy*vy+vz*vz);
p*=1-exp(-0.3*lambdaAngstroms);
// lambdaAngstroms = 3956.0/sqrt(vx*vx+vy*vy+vz*vz);
// p*=1-exp(-0.3*lambdaAngstroms);
%}
// Detector
COMPONENT Image_Detector = PSD_monitor(nx = 2201, ny = 2201, filename = "PSD_output_file.sim",
COMPONENT Image_Detector = PSD_monitor(nx = 221, ny = 221, filename = "PSD_output_file.sim",
xwidth = DetectorWidth, yheight = DetectorWidth, restore_neutron=1)
AT (0, 0, 25.00) RELATIVE Origin
//COMPONENT Image_Detector_BN = PSD_monitor(nx = 2201, ny = 2201, filename = "PSD_output_file_BN.sim",
// xwidth = DetectorWidth, yheight = DetectorWidth, restore_neutron=1)
// WHEN (BNCollScatter>0)
// AT (0, 0, 25.00) RELATIVE Origin
COMPONENT Image_Detector_BN = PSD_monitor(nx = 221, ny = 221, filename = "PSD_output_file_BN.sim",
xwidth = DetectorWidth, yheight = DetectorWidth, restore_neutron=1)
WHEN (BNCollScatter>0)
AT (0, 0, 25.00) RELATIVE Origin
//COMPONENT Penumbra_Detector = PSD_monitor(nx = 1501, ny = 1501, filename = "Penumbra_output_file.sim",
// xwidth = PenumbraDetWidth, yheight = PenumbraDetWidth, restore_neutron=1)
//COMPONENT Penumbra_Detector = Monitor_nD(filename = "Penumbra_output_file", xwidth = PenumbraDetWidth, yheight = PenumbraDetWidth,
COMPONENT Penumbra_Detector = Monitor_nD(filename = "Penumbra_output_file", xwidth = 0.02, yheight = 0.02,
options = "t bins=10000 limits=[0,0.033333]")
COMPONENT Penumbra_Detector = Monitor_nD(filename = "Penumbra_output_file", xwidth = 0.2, yheight = 0.2,
options = "t bins=1000 limits=[0,0.2]")
AT (0, 0, 25.00+1e-9) RELATIVE Origin
//COMPONENT Penumbra_Detector_BN = PSD_monitor(nx = 1501, ny = 1501, filename = "Penumbra_output_file_BN.sim",
......
This diff is collapsed.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment