From 91c052beddf58aabcce3a30f2860b7c893e9f8d2 Mon Sep 17 00:00:00 2001
From: shiroka <shiroka>
Date: Wed, 29 Oct 2008 16:37:52 +0000
Subject: [PATCH] Added two missing header files.

---
 geant4/LEMuSR/include/meyer.h  | 65 ++++++++++++++++++++++++++++++++++
 geant4/LEMuSR/include/yields.h | 39 ++++++++++++++++++++
 2 files changed, 104 insertions(+)
 create mode 100644 geant4/LEMuSR/include/meyer.h
 create mode 100644 geant4/LEMuSR/include/yields.h

diff --git a/geant4/LEMuSR/include/meyer.h b/geant4/LEMuSR/include/meyer.h
new file mode 100644
index 0000000..d9e52ce
--- /dev/null
+++ b/geant4/LEMuSR/include/meyer.h
@@ -0,0 +1,65 @@
+//$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$//*
+//             LOW ENERGY MUON SPIN RELAXATION, ROTATION, RADIATION   Geant4 SIMULATION
+//  ID    : MEYER.hh , v 1.0
+//  AUTHOR: Taofiq PARAISO
+//  DATE  : 2005-04
+//                   
+//$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$//
+//
+//  &           &&&&&&&&&&                            &&&&&&&         &&&&&&&&
+//  &           &                                   &&      &&       &       &&
+//  &           &                   &      &       &                &       &&
+//  &           &&&&&&&            &      &         &&&&&&         &&&&&&&&  
+//  &           &                 &      &&                &      &      &&
+//  &           &                &&     &  &     &&      &&      &        &
+//  &&&&&&&&&&  &&&&&&&&&&      &  &&&&&    &&    &&&&&&&       &        &&   
+//                             &
+//                            &
+//                           &
+//                         &  
+//                             MEYER
+/*
+  fIRST IMPLEMENTATION BY ANLSEM,H. IN FORTRAN
+  C++ CONVERSION T.K.PARAISO 04-2005
+
+  !!! IMPORTANT  !!!
+ 
+  Notice: 
+  Tables definition changes between FORTRAN and C++:
+  1/ Fortran indices start at 1 and C++ indices start at 0
+  2/ Tables are defined as table[column][row] in Fortran
+  table[row][column] in c++
+
+  usefull reference 
+  http://gershwin.ens.fr/vdaniel/Doc-Locale/Langages-Program-Scientific/Fortran/Tutorial/arrays.htm
+
+*/
+//$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$// 
+
+
+#ifndef meyer_h
+#define meyer_h 1
+
+#include <iomanip>
+#include <stdlib.h>
+#include <iostream>
+#include <sstream>
+#include <string>
+#include <fstream>
+#include <ios>
+#include "globals.hh"
+
+class meyer
+{
+ public:
+  meyer();
+ ~meyer();
+
+ void GFunctions(double*, double*, double);
+ void Get_F_Function_Meyer(double tau, double Ekin, double Z1, double Z2, double m1, double m2);
+ void F_Functions_Meyer( double tau,double thetaSchlange,double *f1,double *f2);
+ void Get_F_Function(double tau,double theta, double Ekin, double Z1, double Z2, double m1, double m2, double* F);
+
+};
+
+#endif
diff --git a/geant4/LEMuSR/include/yields.h b/geant4/LEMuSR/include/yields.h
new file mode 100644
index 0000000..bbbc9b1
--- /dev/null
+++ b/geant4/LEMuSR/include/yields.h
@@ -0,0 +1,39 @@
+//$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
+//  Muonium yields as a function of initial mu+ energies.
+//          The method GetYields is used by MuFormation.
+//  Id    : yields.cc, v 1.1
+//  Author: Taofiq PARAISO, T. Shiroka
+//  Date  : 2007-12
+//  Notes : First implemented in Fortran by A. Hofer
+//          C++ conversion by T.K. Paraiso 04-2005
+//          Slight modifications by T. Shiroka
+//$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
+
+#ifndef Yields_h
+#define Yield_h 1
+
+#include "globals.hh"
+
+/*  - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+ The Muonium Yield function as well as the parameters are taken from:
+ M. Gonin, R. Kallenbach, P. Bochsler: "Charge exchange of hydrogen atoms
+ in carbon foils at 0.4 - 120 keV", Rev.Sci.Instrum. 65 (3), March 1994
+- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - */
+
+class Yields
+{
+
+public:
+  Yields(); // Class constructor
+ ~Yields(); // Class destructor
+
+  void GetYields(double E, double mass, double yvector[]);
+
+private:   // Some internal variables
+  double Q_zero, Q_minus, D;
+  double Yield_minus, Yield_zero, Yield_plus;
+
+  double aux1, aux2, aux3; // Auxiliary variables
+};
+
+#endif