Include Lambda(1520) decay
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 23 Aug 2013 11:07:44 +0000 (11:07 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 23 Aug 2013 11:07:44 +0000 (11:07 +0000)
HIJING/hipyset1_35/lucomp_hijing.F
THijing/CMakelibTHijing.pkg
THijing/THijing.cxx
THijing/THijing.h
THijing/openDecayTableHijing.F [new file with mode: 0644]

index 4dae8f8..c1e89f5 100644 (file)
@@ -95,6 +95,9 @@ C...Spin 3/2 baryons.
       ELSEIF(KFLR.EQ.0.AND.KFLS.EQ.4) THEN  
         IF(KFLA.EQ.9.OR.KFLB.EQ.0.OR.KFLB.EQ.9.OR.KFLC.EQ.9) THEN   
         ELSEIF(KFLA.LT.KFLB.OR.KFLB.LT.KFLC) THEN   
+           IF (KFA .EQ. 3124) THEN
+              lucomp_hijing = 387 
+           endif
         ELSEIF(KFLA.GE.6.OR.KFLB.GE.4) THEN 
           LUCOMP_HIJING=80+KFLA    
         ELSE    
index ddb9b63..ecb1aaa 100644 (file)
@@ -28,6 +28,8 @@
 
 set ( SRCS  THijing.cxx AliGenHijing.cxx AliGenBeamGasNew.cxx)
 
+set ( FSRCS  openDecayTableHijing.F)
+
 string(REPLACE ".cxx" ".h" HDRS  "${SRCS}")
 
 set ( DHDR THijingLinkDef.h)
index 257649e..eabf0cf 100644 (file)
@@ -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);
+  
 }
 
 //______________________________________________________________________________
@@ -1315,3 +1333,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);
+}
index 361c606..b9c9b3f 100644 (file)
@@ -196,6 +196,7 @@ public:
    virtual void         Rluget(Int_t lfn, Int_t move=0);
    virtual void         Rluset(Int_t lfn, Int_t move=0);   
    virtual void         Pylist(Int_t flag);
+   virtual void         Luupda(Int_t mupda, Int_t lfn);
    protected:
 
     Float_t      fEfrm;     // Energy in the centre of mass (CMS) or lab-frame (LAB)
diff --git a/THijing/openDecayTableHijing.F b/THijing/openDecayTableHijing.F
new file mode 100644 (file)
index 0000000..751eeb2
--- /dev/null
@@ -0,0 +1,43 @@
+#ifdef __APPLE__
+#ifndef __INTEL_COMPILER
+#define stop CALL EXIT !
+#define STOP CALL EXIT !
+#endif
+#endif
+      SUBROUTINE OPENDECAYTABLE_HIJING(LUN)
+      LOGICAL OPENED,EXISTS
+      LOGICAL LRMSCH, LRD1O2, LTRASP
+      INTEGER LUN
+      CHARACTER*1000 FILNAM
+      CHARACTER*1000 CHROOT
+      CHROOT=' '
+      CALL GETENVF('ALICE_ROOT',CHROOT)
+      LNROOT = LNBLNK(CHROOT)
+      IF(LNROOT.LE.0) THEN
+         FILNAM='decaytable.dat'
+      ELSE
+         FILNAM=CHROOT(1:LNROOT)//'/data/hijingdecaytable.dat'
+      ENDIF
+      INQUIRE(FILE=FILNAM,EXIST=EXISTS)
+      IF(.NOT.EXISTS) THEN
+         PRINT*,'***********************************'
+         PRINT*,'*        OpenDecayTable           *'
+         PRINT*,'*        ---------------          *'
+         PRINT*,'*   File DECAYTABLE.DAT not found *'
+         PRINT*,'*         Program STOP            *'
+         PRINT*,'*   Check ALICE_ROOT environment  *'
+         PRINT*,'*           variable              *'
+         PRINT*,'***********************************'
+         STOP
+      ENDIF
+#ifdef __GFORTRAN__
+      OPEN(LUN,FILE=FILNAM,STATUS='OLD',ACTION='READ')
+#else
+      OPEN(LUN,FILE=FILNAM,STATUS='OLD')
+#endif
+      RETURN
+      END
+      SUBROUTINE CLOSEDECAYTABLE_HIJING(LUN)
+      CLOSE(LUN)
+      RETURN
+      END