Fix for focusing problem when processes use focus to calculate cross section.

Processes can now report the need to do so, and in that case a point will be chosen.
If that process is selected, that point is used, and a weight transformation corrects for the probability of that position.

Also introduced control parameters to disable / enable refraction and reflection just like the Refractor component.

Author: mads-bertelsen

mcstas-comps/share/union-lib.c

@@ -33,6 +33,7 @@ enum process {
   Texture,
   IncoherentPhonon,
   NCrystal,
+  Non,
   Template
 };
 
@@ -465,6 +466,7 @@ struct physics_struct
 char name[256]; // User defined material name
 int interact_control;
 int is_vacuum;
+int any_process_needs_cross_section_focus;
 double my_a;
 int number_of_processes;
 // pointer to array of pointers to physics_sub structures that each describe a scattering process
@@ -502,6 +504,7 @@ char name[256];                // User defined process name
 enum process eProcess;         // enum value corresponding to this process GPU
 double process_p_interact;     // double between 0 and 1 that describes the fraction of events forced to undergo this process. -1 for disable
 int non_isotropic_rot_index;   // -1 if process is isotrpic, otherwise is the index of the process rotation matrix in the volume
+int needs_cross_section_focus; // 1 if physics_my needs to call focus functions, otherwise -1
 Rotation rotation_matrix;      // rotation matrix of process, reported by component in local frame, transformed and moved to volume struct in main
 
 union data_transfer_union data_transfer; // The way to reach the storage space allocated for this process (see examples in process.comp files)

mcstas-comps/share/union-suffix.c

@@ -50,6 +50,11 @@ int physics_my(enum process choice, double *my,double *k_initial, union data_tra
         case NCrystal:
             output = NCrystal_physics_my(my, k_initial, data_transfer, focus_data, _particle);
             break;
+        #endif	
+        #ifdef PROCESS_NON_DETECTOR
+        case Non:
+            output = Non_physics_my(my, k_initial, data_transfer, focus_data, _particle);
+            break;
         #endif
         #ifdef PROCESS_TEMPLATE_DETECTOR
         case Template:
@@ -110,9 +115,14 @@ int physics_scattering(enum process choice, double *k_final, double *k_initial,
             output = NCrystal_physics_scattering(k_final, k_initial, weight, data_transfer, focus_data, _particle);
             break;
         #endif
+        #ifdef PROCESS_NON_DETECTOR
+        case Non:
+            output = Non_physics_scattering(k_final, k_initial, weight, data_transfer, focus_data, _particle);
+            break;
+        #endif			
         #ifdef PROCESS_TEMPLATE_DETECTOR
         case Template:
-            output = Incoherent_physics_scattering(k_final, k_initial, weight, data_transfer, focus_data, _particle);
+            output = Template_physics_scattering(k_final, k_initial, weight, data_transfer, focus_data, _particle);
             break;
         #endif
         default: printf("physics_scattering: No scattering process matches input!\n");

mcstas-comps/union/AF_HB_1D_process.comp

@@ -316,6 +316,10 @@ INITIALIZE
   // Need to specify if this process is isotropic
   //This_process.non_isotropic_rot_index = -1; // Yes (powder)
   This_process.non_isotropic_rot_index =  1;  // No (single crystal)
+  
+  // Need to specify if this process need to use focusing in calculation of inverse penetration depth (physics_my)
+  //This_process.needs_cross_section_focus = 1; // Yes
+  This_process.needs_cross_section_focus =  -1;  // No
 
   // The type of the process must be saved in the global enum process
   This_process.eProcess = AF_HB_1D;

mcstas-comps/union/IncoherentPhonon_process.comp

@@ -4649,6 +4649,10 @@ INITIALIZE
   // Need to specify if this process is isotropic
   This_process.non_isotropic_rot_index = -1; // Yes (powder)
   //This_process.non_isotropic_rot_index =  1;  // No (single crystal)
+  
+  // Need to specify if this process need to use focusing in calculation of inverse penetration depth (physics_my)
+  //This_process.needs_cross_section_focus = 1; // Yes
+  This_process.needs_cross_section_focus =  -1;  // No
 
   // The type of the process must be saved in the global enum process
   This_process.eProcess = IncoherentPhonon;

mcstas-comps/union/Incoherent_process.comp

@@ -146,6 +146,10 @@ INITIALIZE
   // Need to specify if this process is isotropic
   This_process.non_isotropic_rot_index = -1; // Yes (powder)
   //This_process.non_isotropic_rot_index =  1;  // No (single crystal)
+  
+  // Need to specify if this process need to use focusing in calculation of inverse penetration depth (physics_my)
+  //This_process.needs_cross_section_focus = 1; // Yes
+  This_process.needs_cross_section_focus =  -1;  // No
 
   // The type of the process must be saved in the global enum process
   This_process.eProcess = Incoherent;

mcstas-comps/union/NCrystal_process.comp

@@ -372,6 +372,10 @@ INITIALIZE
   // Need to specify if this process is isotropic
   //This_process.non_isotropic_rot_index = -1; // Yes (powder)
   This_process.non_isotropic_rot_index =  1;  // No (single crystal)
+  
+  // Need to specify if this process need to use focusing in calculation of inverse penetration depth (physics_my)
+  //This_process.needs_cross_section_focus = 1; // Yes
+  This_process.needs_cross_section_focus =  -1;  // No
 
   // The type of the process must be saved in the global enum process
   This_process.eProcess = NCrystal;

mcstas-comps/union/Non_process.comp

@@ -101,11 +101,13 @@ int Non_physics_scattering(double *k_final, double *k_initial, double *weight, u
 // These lines help with future error correction, and tell other Union components
 //  that at least one process have been defined.
 #ifndef PROCESS_DETECTOR
-    // Obsolete
-    //struct pointer_to_global_process_list global_process_list = {0,NULL};
     #define PROCESS_DETECTOR dummy
 #endif
 
+#ifndef PROCESS_NON_DETECTOR
+    #define PROCESS_NON_DETECTOR dummy
+#endif
+
 %}
 
 DECLARE
@@ -130,6 +132,13 @@ INITIALIZE
   // Need to specify if this process is isotropic
   This_process.non_isotropic_rot_index = -1; // Yes (powder)
   //This_process.non_isotropic_rot_index =  1;  // No (single crystal)
+  
+  // Need to specify if this process need to use focusing in calculation of inverse penetration depth (physics_my)
+  //This_process.needs_cross_section_focus = 1; // Yes
+  This_process.needs_cross_section_focus =  -1;  // No
+  
+  // The type of the process must be saved in the global enum process
+  This_process.eProcess = Non;
 
   // Packing the data into a structure that is transported to the main component
   This_process.data_transfer.pointer_to_a_Non_physics_storage_struct = &Non_storage;

mcstas-comps/union/PhononSimple_process.comp

@@ -515,6 +515,10 @@ INITIALIZE
   // Need to specify if this process is isotropic
   //This_process.non_isotropic_rot_index = -1; // Yes (powder)
   This_process.non_isotropic_rot_index =  1;  // No (single crystal)
+  
+  // Need to specify if this process need to use focusing in calculation of inverse penetration depth (physics_my)
+  This_process.needs_cross_section_focus = 1; // Yes
+  //This_process.needs_cross_section_focus =  -1;  // No
 
   // The type of the process must be saved in the global enum process
   This_process.eProcess = PhononSimple;

mcstas-comps/union/Powder_process.comp

@@ -749,6 +749,10 @@ INITIALIZE
   // Need to specify if this process is isotropic
   This_process.non_isotropic_rot_index = -1; // Yes (powder)
   //This_process.non_isotropic_rot_index =  1;  // No (single crystal)
+  
+  // Need to specify if this process need to use focusing in calculation of inverse penetration depth (physics_my)
+  //This_process.needs_cross_section_focus = 1; // Yes
+  This_process.needs_cross_section_focus =  -1;  // No
 
   // The type of the process must be saved in the global enum process
   This_process.eProcess = Powder;

mcstas-comps/union/Single_crystal_process.comp

@@ -1077,6 +1077,10 @@ INITIALIZE
   // Need to specify if this process is isotropic
   //This_process.non_isotropic_rot_index = -1; // Yes (powder)
   This_process.non_isotropic_rot_index =  1;  // No (single crystal)
+  
+  // Need to specify if this process need to use focusing in calculation of inverse penetration depth (physics_my)
+  //This_process.needs_cross_section_focus = 1; // Yes
+  This_process.needs_cross_section_focus =  -1;  // No
 
   // The type of the process must be saved in the global enum process
   This_process.eProcess = Single_crystal;