corrected TRD/TOF MV position

[u/mrichter/AliRoot.git] / THijing / THijing.cxx

diff --git a/THijing/THijing.cxx b/THijing/THijing.cxx
index a0f2fe13187669cb0095a776c5ac5b91eaefc1da..0ba5b14b2499bdcf21e603176043199b6e961527 100644 (file)
--- a/THijing/THijing.cxx
+++ b/THijing/THijing.cxx
@@ -13,7 +13,7 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
-/* $Id$ */
+/* $Id: THijing.cxx 64005 2013-08-28 12:23:17Z hristov $ */
 
 ////////////////////////////////////////////////////////////////////////////////
 //                                                                            //
@@ -56,6 +56,9 @@
 # define rluget_hijing rluget_hijing_
 # define rluset_hijing rluset_hijing_
 # define lulist_hijing lulist_hijing_
+# define luupda_hijing luupda_hijing_
+# define opendecaytable_hijing opendecaytable_hijing_
+# define closedecaytable_hijing closedecaytable_hijing_
 # define type_of_call
 #else
 # define hijset HIJSET
@@ -64,6 +67,9 @@
 # define rluget_hijing RLUGET_HIJING
 # define rluset_hijing RLUSET_HIJING
 # define lulist_hijing LULIST_HIJING
+# define luupda_hijing LUUPDA_HIJING
+# define opendecaytable_hijing OPENDECAYTABLE_HIJING
+# define closedecaytable_hijing CLOSEDECAYTABLE_HIJING
 # define type_of_call _stdcall
 #endif
 
@@ -80,6 +86,11 @@ extern "C" void type_of_call rluget_hijing(Int_t & lfn, Int_t & move);
 
 extern "C" void type_of_call rluset_hijing(Int_t & lfn, Int_t & move);
 extern "C" void type_of_call lulist_hijing(Int_t &);
+extern "C" void type_of_call luupda_hijing(Int_t &, Int_t &);
+extern "C" void type_of_call 
+          opendecaytable_hijing(Int_t&);
+extern "C" void type_of_call 
+          closedecaytable_hijing(Int_t&);
 
 #else
 #endif
@@ -125,6 +136,13 @@ THijing::THijing(Float_t efrm, const char *frame="CMS",
 // THijing constructor: 
 // Note that there may be only one functional THijing object
 // at a time, so it's not use to create more than one instance of it.
+//  
+// Read decaytable
+  Int_t lun = 15;
+  opendecaytable_hijing(lun);
+  Luupda(2,lun);
+  closedecaytable_hijing(lun);
+  
 }
 
 //______________________________________________________________________________
@@ -556,6 +574,8 @@ Int_t  THijing::GetN0() const
    return HIMAIN1.n0;
 
 }
+
+
 //______________________________________________________________________________
 Int_t  THijing::GetN01() const
 {
@@ -592,6 +612,7 @@ Float_t  THijing::GetBB() const
 
 }
 
+
 //====================== access to common HIMAIN2 ===============================
 
 //______________________________________________________________________________
@@ -1161,19 +1182,47 @@ void  THijing::SetMDCY(Int_t key1, Int_t key2, Int_t   parm)
   // Set value of array MDCY
   if ( key1 < 1 || key1 > 500) {
     printf("ERROR in THijing::SetMDCY(key1, key2, parm):\n");
-    printf("      key1=%i is out of range [1..200]\n", key1);
+    printf("      key1=%i is out of range [1..500]\n", key1);
   } else if ( key2 < 1 || key2 > 3) {
     printf("ERROR in THijing::SetMDCY(key1, key2, parm):\n");
-    printf("      key2=%i is out of range [1..200]\n", key2);
+    printf("      key2=%i is out of range [1..3]\n", key2);
   } else {
     LUDAT3_HIJING.mdcy[key2-1][key1-1] = parm;
   }
   
 }
 
+void  THijing::SetMDME(Int_t key1, Int_t key2, Int_t   parm)
+{
+  // Set value of array MDME
+  if ( key1 < 1 || key1 > 2000) {
+    printf("ERROR in THijing::SetMDME(key1, key2, parm):\n");
+    printf("      key1=%i is out of range [1..2000]\n", key1);
+  } else if ( key2 < 1 || key2 > 2) {
+    printf("ERROR in THijing::SetMDME(key1, key2, parm):\n");
+    printf("      key2=%i is out of range [1..3]\n", key2);
+  } else {
+    LUDAT3_HIJING.mdme[key2-1][key1-1] = parm;
+  }
+  
+}
+
+void  THijing::SetBRAT(Int_t key, Float_t   parm)
+{
+  // Set value of array MDME
+  if ( key < 1 || key > 2000) {
+    printf("ERROR in THijing::SetBRAT(key, parm):\n");
+    printf("      key=%i is out of range [1..2000]\n", key);
+  } else {
+    LUDAT3_HIJING.brat[key-1] = parm;
+  }
+  
+}
+
 Int_t THijing::GetMDCY(Int_t key1, Int_t key2) 
 {
   // Get value of array MDCY
+
   if ( key1 < 1 || key1 > 500) {
     printf("ERROR in THijing::GetMDCY(key1, key2, parm):\n");
     printf("      key1=%i is out of range [1..200]\n", key1);
@@ -1186,6 +1235,7 @@ Int_t THijing::GetMDCY(Int_t key1, Int_t key2)
     return (LUDAT3_HIJING.mdcy[key2-1][key1-1]);
   }
 
+
 }
 //====================== access to Hijing subroutines =========================
 
@@ -1298,3 +1348,9 @@ void  THijing::Pylist(Int_t flag)
 // call lulist
   lulist_hijing(flag);
 }
+
+void  THijing::Luupda(Int_t mupda, Int_t lfn)
+{
+// call luupda
+  luupda_hijing(mupda, lfn);
+}