Add Hidden cylinder test instrument

Author: Lomholy

mcstas-comps/examples/Tests_union/Hidden_Cylinder/Hidden_Cylinder.instr

@@ -0,0 +1,164 @@
+/********************************************************************************
+* 
+* McStas, neutron ray-tracing package
+*         Copyright (C) 1997-2008, All rights reserved
+*         Risoe National Laboratory, Roskilde, Denmark
+*         Institut Laue Langevin, Grenoble, France
+* 
+* This file was written by McStasScript, which is a 
+* python based McStas instrument generator written by 
+* Mads Bertelsen in 2019 while employed at the 
+* European Spallation Source Data Management and 
+* Software Centre
+* 
+* Instrument: Hidden_Cylinder
+* 
+* %Identification
+* Written by: Daniel Lomholt Christensen
+* Date: 15:54:48 on January 20, 2026
+* Origin: UCPH@NBI, Funded by ACTNXT
+* %INSTRUMENT_SITE: Tests_union
+* 
+* 
+* %Description
+* A small test instrument testing the placement of one cylinder with a low priority
+* fully inside another cylinder with a higher priority. This was added to test 
+* for an mcdisplay error that used to occur, when a cylinder is completely 
+* enclosed in another cylinder.
+* To test this simply do
+* a "mcdisplay Hidden_Cylinder.instr -c -y"
+* 
+* %Example: pin_rad=0.0025 Detector: det_I=3.42884e-06
+* 
+* %Parameters
+* pin_rad: [m] Radius of source
+* d: [m] Distance between source and sample center
+* detector_x: [m] Detector Width
+* detector_y: [m] Detector Height
+* det_dist: [m] Distance between detector and sample center
+* 
+* %Link 
+* 
+* %End 
+********************************************************************************/
+
+DEFINE INSTRUMENT imaging (
+double pin_rad = 0.0025, // Radius of source
+double d = 10, // Distance between source and sample center
+double detector_x = 0.01, // Detector Width
+double detector_y = 0.04, // Detector Height
+double det_dist = 0.005 // Distance between detector and sample center
+)
+
+DECLARE 
+%{
+%}
+
+INITIALIZE 
+%{
+%}
+
+TRACE 
+COMPONENT Origin = Progress_bar()
+AT (0, 0, 0) ABSOLUTE
+
+COMPONENT source = Source_simple(
+ radius = pin_rad, yheight = 0,
+ xwidth = 0, dist = d,
+ focus_xw = detector_x, focus_yh = detector_y,
+ lambda0 = 4, dlambda = 3)
+AT (0, 0, 0) RELATIVE Origin
+ROTATED (0.0, 0.0, 0.0) RELATIVE Origin
+
+COMPONENT entry_mon = PSD_monitor(
+ nx = 200, ny = 200,
+ xwidth = detector_x, yheight = detector_x,
+ restore_neutron = 1)
+AT (0, 0, 0.0001) RELATIVE source
+ROTATED (0.0, 0.0, 0.0) RELATIVE source
+
+COMPONENT arm_sample_position = Arm()
+AT (0, 0, d) RELATIVE source
+ROTATED (0, 0, 0) RELATIVE Origin
+
+COMPONENT init = Union_init()
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT Sample_inc = Incoherent_process(
+ sigma = 0.4 * 2, unit_cell_volume = 24.04)
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT Sample_pow = Powder_process(
+ reflections = "Fe.laz")
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT Sample = Union_make_material(
+ process_string = "Sample_inc,Sample_pow", my_absorption = 10*2.56 * 2 / 24.04 * 100)
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT gas_process = Incoherent_process(
+ sigma = 0.4 * 2, unit_cell_volume = 240000.04)
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT sample_gas = Union_make_material(
+ process_string = "gas_process", my_absorption = 21.29783693843594)
+AT (0, 0, 0) RELATIVE Sample
+ROTATED (0.0, 0.0, 0.0) RELATIVE Sample
+
+COMPONENT sample_cyl_metal = Union_cylinder(
+ material_string = "Sample", priority = 50,
+ radius = 0.005 / 2, yheight = 0.034)
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT sample_cyl_gas = Union_cylinder(
+ material_string = "sample_gas", priority = 51,
+ radius = 0.009 / 2, yheight = 0.053,
+ visualize = 1)
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+  COMPONENT logger_space_zx = Union_logger_2D_space(
+ D_direction_1 = "z", D1_min = -0.05,
+ D1_max = 0.05, n1 = 300,
+ D_direction_2 = "x", D2_min = -0.05,
+ D2_max = 0.05, n2 = 300,
+ filename = "logger_zx.dat")
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT logger_space_zy = Union_logger_2D_space(
+ D_direction_1 = "z", D1_min = -0.05,
+ D1_max = 0.05, n1 = 300,
+ D_direction_2 = "y", D2_min = -0.05,
+ D2_max = 0.05, n2 = 300,
+ filename = "logger_zy.dat")
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT Sample_environment = Union_master(
+ verbal = 1)
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT stop = Union_stop()
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+COMPONENT det = PSD_monitor(
+ nx = 1000, ny = 1000,
+ filename = "det", xwidth = detector_x,
+ yheight = detector_y)
+AT (0, 0, 0.034 / 2 + det_dist) RELATIVE arm_sample_position
+ROTATED (0.0, 0.0, 0.0) RELATIVE arm_sample_position
+
+FINALLY 
+%{
+%}
+
+END