Add SetERAN function

diff --git a/TGeant3/TGeant3.cxx b/TGeant3/TGeant3.cxx
index 40a99ad9bde67a8bf6bf551d900418df116d66b2..6397d405e07238574607c5aa2816c66e71598ae6 100644 (file)
--- a/TGeant3/TGeant3.cxx
+++ b/TGeant3/TGeant3.cxx
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.20  2000/01/12 11:29:27  fca
+Close material file
+
 Revision 1.19  1999/12/17 09:03:12  fca
 Introduce a names array
 
@@ -550,6 +553,7 @@ void TGeant3::LoadAddress()
   gcomad(PASSCHARD("GCBANK"),(int*&) fGcbank  PASSCHARL("GCBANK"));
   gcomad(PASSCHARD("GCLINK"),(int*&) fGclink  PASSCHARL("GCLINK"));
   gcomad(PASSCHARD("GCCUTS"),(int*&) fGccuts  PASSCHARL("GCCUTS"));
+  gcomad(PASSCHARD("GCMULO"),(int*&) fGcmulo  PASSCHARL("GCMULO"));
   gcomad(PASSCHARD("GCFLAG"),(int*&) fGcflag  PASSCHARL("GCFLAG"));
   gcomad(PASSCHARD("GCKINE"),(int*&) fGckine  PASSCHARL("GCKINE"));
   gcomad(PASSCHARD("GCKING"),(int*&) fGcking  PASSCHARL("GCKING"));
@@ -3179,6 +3183,20 @@ void TGeant3::SetDRAY(Int_t par)
   fGcphys->idray = par;
 }
  
+//_____________________________________________________________________________
+void TGeant3::SetERAN(Float_t ekmin, Float_t ekmax, Int_t nekbin)
+{
+  //
+  //  To control cross section tabulations
+  //   ekmin = minimum kinetic energy in GeV
+  //   ekmax = maximum kinetic energy in GeV
+  //   nekbin = number of logatithmic bins (<200)
+  //  
+  fGcmulo->ekmin = ekmin;
+  fGcmulo->ekmax = ekmax;
+  fGcmulo->nekbin = nekbin;
+}
+ 
 //_____________________________________________________________________________
 void TGeant3::SetHADR(Int_t par)
 {

diff --git a/TGeant3/TGeant3.h b/TGeant3/TGeant3.h
index be6882cb8d49543f72bb90f1399ce1e8e5808be7..77b96e0e400a12a33c8c9c2e63c6ed789c4e7c3f 100644 (file)
--- a/TGeant3/TGeant3.h
+++ b/TGeant3/TGeant3.h
@@ -274,6 +274,26 @@ typedef struct {
   Float_t gcuts[5]; 
 } Gccuts_t; 
 
+//----------GCMULO
+//      COMMON/GCMULO/SINMUL(101),COSMUL(101),SQRMUL(101),OMCMOL,CHCMOL
+//     +  ,EKMIN,EKMAX,NEKBIN,NEK1,EKINV,GEKA,GEKB,EKBIN(200),ELOW(200)
+typedef struct {
+  Float_t sinmul[101];
+  Float_t cosmul[101];
+  Float_t sqrmul[101];
+  Float_t omcmol;
+  Float_t chcmol;
+  Float_t ekmin;
+  Float_t ekmax;
+  Int_t   nekbin;
+  Int_t   nek1;
+  Float_t ekinv;
+  Float_t geka;
+  Float_t gekb;
+  Float_t ekbin[200];
+  Float_t elow[200];
+} Gcmulo_t;
+
 //----------GCPHYS
 //      COMMON/GCPHYS/IPAIR,SPAIR,SLPAIR,ZINTPA,STEPPA
 //     +             ,ICOMP,SCOMP,SLCOMP,ZINTCO,STEPCO
@@ -503,6 +523,7 @@ private:
   Gcbank_t *fGcbank;
   Gclink_t *fGclink; 
   Gccuts_t *fGccuts; 
+  Gcmulo_t *fGcmulo; 
   Gcmate_t *fGcmate; 
   Gctpol_t *fGctpol; 
   Gcnum_t  *fGcnum; 
@@ -621,6 +642,7 @@ public:
   virtual Gcbank_t* Gcbank() const {return fGcbank;}
   virtual Gclink_t* Gclink() const {return fGclink;}
   virtual Gccuts_t* Gccuts() const {return fGccuts;}
+  virtual Gcmulo_t* Gcmulo() const {return fGcmulo;}
   virtual Gcmate_t* Gcmate() const {return fGcmate;}
   virtual Gctpol_t* Gctpol() const {return fGctpol;}
   virtual Gcnum_t* Gcnum() const {return fGcnum;}
@@ -782,6 +804,8 @@ public:
    virtual  void  SetDCAY(Int_t par=1);
    virtual  void  SetDEBU(Int_t emin=1, Int_t emax=999, Int_t emod=1);
    virtual  void  SetDRAY(Int_t par=1);
+   virtual  void  SetERAN(Float_t ekmin=1.e-5, Float_t ekmax=1.e4,
+                         Int_t nekbin=90);
    virtual  void  SetHADR(Int_t par=1);
    virtual  void  SetKINE(Int_t kine, Float_t xk1=0, Float_t xk2=0, Float_t xk3=0, Float_t xk4=0,
                          Float_t xk5=0, Float_t xk6=0, Float_t xk7=0, Float_t xk8=0, Float_t xk9=0,

diff --git a/TGeant3/TGeant3Dummy.cxx b/TGeant3/TGeant3Dummy.cxx
index 2421d1633647788cf8586ddee12f85ba99552c95..fb528c7127cbe3224e13402926b1c727874c52cd 100644 (file)
--- a/TGeant3/TGeant3Dummy.cxx
+++ b/TGeant3/TGeant3Dummy.cxx
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.14  1999/11/03 16:31:36  fca
+Dummy ProdProcess should return a value
+
 Revision 1.13  1999/11/03 13:17:08  fca
 Have ProdProcess return const char*
 
@@ -253,6 +256,7 @@ void TGeant3::SetCUTS(Float_t,Float_t,Float_t,Float_t,
 void TGeant3::SetDCAY(Int_t)  {}
 void TGeant3::SetDEBU(Int_t, Int_t, Int_t)  {}
 void TGeant3::SetDRAY(Int_t)  {}
+void TGeant3::SetERAN(Float_t, Float_t, Int_t) {}
 void TGeant3::SetHADR(Int_t)  {}
 void TGeant3::SetKINE(Int_t, Float_t, Float_t, Float_t, Float_t,
                      Float_t, Float_t, Float_t, Float_t, Float_t,