Add many meshes instrument that fails without extra allocation in dummyint

Author: Lomholy

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

@@ -0,0 +1,171 @@
+/********************************************************************************
+* 
+* 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: Testing placing many meshes in the same instrument file
+* 
+* %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 many overlapping 
+* mesh components.
+* Test with:
+* mcdisplay -y -c
+*
+* %Example: pin_rad=0.0025 Detector:det_I=4.63384e-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
+* crack_width: [m] Height masking crack should be raised from the original crack
+* 
+* %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
+double crack_width = 0.003 // Height masking crack should be raised from the original crack
+)
+DEPENDENCY " @NCRYSTALFLAGS@ -DFUNNEL "
+
+DECLARE 
+%{
+%}
+
+INITIALIZE 
+%{
+// Start of initialize for generated imaging
+%}
+
+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 = NCrystal_process(
+ cfg = "gasmix::H2/10bar")
+AT (0, 0, 0) RELATIVE Sample
+ROTATED (0.0, 0.0, 0.0) RELATIVE Sample
+
+COMPONENT sample_gas = Union_make_material(
+ process_string = "gas_process", my_absorption = 21.29783693843594)
+AT (0, 0, 0) RELATIVE gas_process
+ROTATED (0.0, 0.0, 0.0) RELATIVE gas_process
+
+COMPONENT arm_crack = Arm()
+AT (0, 0, 0) RELATIVE arm_sample_position
+ROTATED (90, 0, 0) RELATIVE arm_sample_position
+
+COMPONENT masking_crack = Union_mesh(
+ filename = "test.stl", material_string = "Sample",
+ priority = 140, skip_convex_check = 1)
+AT (0, 0, 0.0005) RELATIVE arm_crack
+ROTATED (0, 0, 0) RELATIVE arm_crack
+
+COMPONENT crack = Union_mesh(
+ filename = "test.stl", material_string = "sample_gas",
+ priority = 139, skip_convex_check = 1,
+ coordinate_scale = 0.0009811453007353597)
+AT (0, 0, 0) RELATIVE arm_crack
+ROTATED (0, 0, 0) RELATIVE arm_crack
+
+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 
+%{
+// Start of finally for generated imaging
+%}
+
+END