Format McXtrace union

Author: willend

mcxtrace-comps/union/Incoherent_process.comp

@@ -57,106 +57,107 @@ SETTING PARAMETERS(sigma=5.08,f_QE=0,gamma=0,packing_factor=1,unit_cell_volume=1
 
 SHARE
 %{
-#ifndef Union
-#define Union $Revision: 0.8 $
+  #ifndef Union
+  #define Union $Revision: 0.8 $
 
-#include "Union_functions.c"
-#include "Union_initialization.c"
+  #include "Union_functions.c"
+  #include "Union_initialization.c"
 
-#endif
+  #endif
 
-
-struct Incoherent_physics_storage_struct{
+  struct Incoherent_physics_storage_struct {
     // Variables that needs to be transfered between any of the following places:
     // The initialize in this component
     // The function for calculating my
     // The function for calculating scattering
-    
+
     double my_scattering;
     double QE_sampling_frequency;
     double lorentzian_width;
-    
-};
+  };
 
-// Function for calculating my in Incoherent case
-int Incoherent_physics_my(double *my,double *k_initial, union data_transfer_union data_transfer, struct focus_data_struct *focus_data) {
+  // Function for calculating my in Incoherent case
+  int
+  Incoherent_physics_my (double* my, double* k_initial, union data_transfer_union data_transfer, struct focus_data_struct* focus_data) {
     *my = data_transfer.pointer_to_a_Incoherent_physics_storage_struct->my_scattering;
     return 1;
-};
+  };
 
-// Function for basic incoherent scattering event
-int Incoherent_physics_scattering(double *k_final, double *k_initial, double *weight, union data_transfer_union data_transfer, struct focus_data_struct *focus_data) {
+  // Function for basic incoherent scattering event
+  int
+  Incoherent_physics_scattering (double* k_final, double* k_initial, double* weight, union data_transfer_union data_transfer,
+                                 struct focus_data_struct* focus_data) {
 
-    //New version of incoherent scattering
-    double k_length = sqrt(k_initial[0]*k_initial[0]+k_initial[1]*k_initial[1]+k_initial[2]*k_initial[2]);
+    // New version of incoherent scattering
+    double k_length = sqrt (k_initial[0] * k_initial[0] + k_initial[1] * k_initial[1] + k_initial[2] * k_initial[2]);
 
     Coords k_out;
     // Here is the focusing system in action, get a vector
     double solid_angle;
-    focus_data->focusing_function(&k_out,&solid_angle,focus_data);
-    NORM(k_out.x,k_out.y,k_out.z);
-    *weight *= solid_angle*0.25/PI;
-    
-    double E_i,dE,E_f;
-    
-    if (rand01() < data_transfer.pointer_to_a_Incoherent_physics_storage_struct->QE_sampling_frequency) {
-      E_i = K2E*k_length;
-      dE = data_transfer.pointer_to_a_Incoherent_physics_storage_struct->lorentzian_width*tan(PI/2*randpm1());
+    focus_data->focusing_function (&k_out, &solid_angle, focus_data);
+    NORM (k_out.x, k_out.y, k_out.z);
+    *weight *= solid_angle * 0.25 / PI;
+
+    double E_i, dE, E_f;
+
+    if (rand01 () < data_transfer.pointer_to_a_Incoherent_physics_storage_struct->QE_sampling_frequency) {
+      E_i = K2E * k_length;
+      dE = data_transfer.pointer_to_a_Incoherent_physics_storage_struct->lorentzian_width * tan (PI / 2 * randpm1 ());
       E_f = E_i + dE;
       if (E_f <= 0)
         return 0;
-      k_length = E_f*E2K;
+      k_length = E_f * E2K;
     }
-    
-    k_final[0] = k_out.x*k_length; k_final[1] = k_out.y*k_length; k_final[2] = k_out.z*k_length;
-    return 1;
-};
 
+    k_final[0] = k_out.x * k_length;
+    k_final[1] = k_out.y * k_length;
+    k_final[2] = k_out.z * k_length;
+    return 1;
+  };
 %}
 
 DECLARE
 %{
-// Needed for transport to the main component
-struct global_process_element_struct global_process_element;
-struct scattering_process_struct This_process;
-
-#ifndef PROCESS_DETECTOR
-	//struct pointer_to_global_process_list global_process_list = {0,NULL};
-	#define PROCESS_DETECTOR dummy
-#endif
-
-// Declare for this component, to do calculations on the input / store in the transported data
-struct Incoherent_physics_storage_struct Incoherent_storage;
-double effective_my_scattering;
-
+  // Needed for transport to the main component
+  struct global_process_element_struct global_process_element;
+  struct scattering_process_struct This_process;
+
+  #ifndef PROCESS_DETECTOR
+  // struct pointer_to_global_process_list global_process_list = {0,NULL};
+  #define PROCESS_DETECTOR dummy
+  #endif
+
+  // Declare for this component, to do calculations on the input / store in the transported data
+  struct Incoherent_physics_storage_struct Incoherent_storage;
+  double effective_my_scattering;
 %}
 
 INITIALIZE
 %{
   // Initialize done in the component
-  effective_my_scattering = ((packing_factor/unit_cell_volume) * 100 * sigma);
+  effective_my_scattering = ((packing_factor / unit_cell_volume) * 100 * sigma);
   Incoherent_storage.my_scattering = effective_my_scattering;
-  
+
   Incoherent_storage.QE_sampling_frequency = f_QE;
   Incoherent_storage.lorentzian_width = gamma;
 
   // 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)
+  // This_process.non_isotropic_rot_index =  1;  // No (single crystal)
 
   // Packing the data into a structure that is transported to the main component
-  sprintf(This_process.name,"%s",NAME_CURRENT_COMP);
+  sprintf (This_process.name, "%s", NAME_CURRENT_COMP);
   This_process.process_p_interact = interact_fraction;
   This_process.data_transfer.pointer_to_a_Incoherent_physics_storage_struct = &Incoherent_storage;
-  //This_process.data_transfer.pointer_to_a_Incoherent_physics_storage_struct->my_scattering = effective_my_scattering;
+  // This_process.data_transfer.pointer_to_a_Incoherent_physics_storage_struct->my_scattering = effective_my_scattering;
   This_process.probability_for_scattering_function = &Incoherent_physics_my;
   This_process.scattering_function = &Incoherent_physics_scattering;
 
   // This will be the same for all process's, and can thus be moved to an include.
-  sprintf(global_process_element.name,"%s",NAME_CURRENT_COMP);
+  sprintf (global_process_element.name, "%s", NAME_CURRENT_COMP);
   global_process_element.component_index = INDEX_CURRENT_COMP;
   global_process_element.p_scattering_process = &This_process;
-  add_element_to_process_list(&global_process_list,global_process_element);
+  add_element_to_process_list (&global_process_list, global_process_element);
  %}
 
 TRACE
@@ -165,8 +166,7 @@ TRACE
 
 FINALLY
 %{
-// Since the process and it's storage is a static allocation, there is nothing to deallocate
-
+  // Since the process and it's storage is a static allocation, there is nothing to deallocate
 %}
 
 END