AliGenGSIlib with parametersations for GSI physics simulation added (YF, MI)

diff --git a/EVGEN/AliGenGSIlib.cxx b/EVGEN/AliGenGSIlib.cxx
new file mode 100644 (file)
index 0000000..fb923d3
--- /dev/null
+++ b/EVGEN/AliGenGSIlib.cxx
@@ -0,0 +1,274 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+/*
+$Log$
+Revision 1.7  2000/05/02 08:12:13  morsch
+Coding rule violations corrected.
+
+Revision 1.6  1999/09/29 09:24:14  fca
+Introduction of the Copyright and cvs Log
+
+*/
+
+#include "AliGenGSIlib.h"
+#include "AliRun.h"
+#include "iostream.h"
+
+ClassImp(AliGenGSIlib)
+
+  Bool_t AliGenGSIlib::fgDebug =kFALSE;
+//                      Upsilon
+//
+//
+//                  pt-distribution
+//____________________________________________________________
+Double_t AliGenGSIlib::PtUpsilonRitman( Double_t *px, Double_t *dummy )
+{
+  // Upsilon pT
+  /*   AliGenMUONlib parametrisation
+  const Double_t kpt0 = 5.3;
+  const Double_t kxn  = 2.5;
+  Double_t x=*px;
+  //
+  Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+  return x/TMath::Power(pass1,kxn);
+  */
+  if (fgDebug) cout<<"Ritman Pt paramtrisation\n";
+  const Double_t kpt0 = 4.7;
+  const Double_t kxn  = 3.5;
+  Double_t x=*px;
+  //
+  Double_t pass1 = 1.+((x*x)/(kpt0*kpt0));
+  return x/TMath::Power(pass1,kxn);
+   
+}
+//
+//                    y-distribution
+//
+//____________________________________________________________
+Double_t AliGenGSIlib::YUpsilonRitman(Double_t *py, Double_t *dummy)
+{
+
+    // Upsilon y
+
+/*   original from AliGenMUON
+    const Double_t ky0 = 3.;
+    const Double_t kb=1.;
+    Double_t yu;
+    Double_t y=TMath::Abs(*py);
+    //
+    if (y < ky0)
+    yu=kb;
+    else
+    yu=kb*TMath::Exp(-(y-ky0)*(y-ky0)/2);
+    return yu;
+    */
+  if (fgDebug) cout<<"Ritman Y paramtrisation\n";
+  return 0.003;  //GSI parametrisation 
+}
+//                 particle composition
+//
+Int_t AliGenGSIlib::IpUpsilonRitman()
+{
+// y composition
+  if (fgDebug) cout<<"Ritman Ip paramtrisation\n";
+  return 553;     
+}
+
+
+Double_t AliGenGSIlib::PtUpsilonKarel( Double_t *px, Double_t *dummy )
+{
+  // Upsilon pT
+  //to implement
+  if (fgDebug) cout<<"Karel Pt paramtrisation\n";
+
+  return 0.;   
+}
+//
+//                    y-distribution
+//
+//____________________________________________________________
+Double_t AliGenGSIlib::YUpsilonKarel(Double_t *py, Double_t *dummy)
+{
+  
+  // Upsilon y
+//to implement
+  if (fgDebug) cout<<"Karel Y paramtrisation\n";
+  return 0.003;  //Karel parametrisation 
+}
+
+//                 particle composition
+//
+
+Int_t AliGenGSIlib::IpUpsilonKarel()
+{
+  // y composition//
+  //to implement
+  if (fgDebug) cout<<"Karel Ip paramtrisation\n";
+  return 553;     
+}
+
+
+
+Double_t AliGenGSIlib::PtUpsilonMUON( Double_t *px, Double_t *dummy )
+{
+// Upsilon pT
+  const Double_t kpt0 = 5.3;
+  const Double_t kxn  = 2.5;
+  Double_t x=*px;
+  //
+  Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+  return x/TMath::Power(pass1,kxn);
+}
+//
+//                    y-distribution
+//
+//____________________________________________________________
+Double_t AliGenGSIlib::YUpsilonMUON(Double_t *py, Double_t *dummy)
+{
+// Upsilon y
+  const Double_t ky0 = 3.;
+  const Double_t kb=1.;
+  Double_t yu;
+  Double_t y=TMath::Abs(*py);
+  //
+  if (y < ky0)
+    yu=kb;
+  else
+    yu=kb*TMath::Exp(-(y-ky0)*(y-ky0)/2);
+  return yu;
+}
+//                 particle composition
+//
+Int_t AliGenGSIlib::IpUpsilonMUON()
+{
+// y composition
+    return 553;
+}
+
+
+typedef Double_t (*GenFunc) (Double_t*,  Double_t*);
+
+typedef Int_t (*GenFuncIp) ();
+
+
+
+GenFunc AliGenGSIlib::GetPt(Param_t param, const char * tname)
+{
+// Return pointer to pT parameterisation
+  GenFunc func=0;
+  TString sname(tname);
+  switch (param) 
+    {
+    case upsilon_p:
+      if (sname=="MUON"){
+       func= PtUpsilonMUON;
+       break;
+      }
+      if (sname=="RITMAN"){
+       func=PtUpsilonRitman;
+       break;
+      }
+      if (sname=="KAREL"){
+       func=PtUpsilonKarel;
+       break;
+      }
+      break;
+       func=0;
+        printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
+    default:
+        func=0;
+        printf("<AliGenGSIlib::GetPt> unknown parametrisation\n");
+    }
+    return func;
+}
+
+
+
+GenFunc AliGenGSIlib::GetY(Param_t param, const char * tname)
+{
+  // Return pointer to y- parameterisation
+   GenFunc func=0;
+    TString sname(tname);
+    switch (param) 
+    {
+    case upsilon_p:
+      if (sname=="MUON"){
+       func= YUpsilonMUON;
+       break;
+      }
+      if (sname=="RITMAN"){
+       func=YUpsilonRitman;
+       break;
+      }
+      if (sname=="KAREL"){
+       func=YUpsilonKarel;
+       break;
+      }
+      func=0;
+      printf("<AliGenGSIlib::GetY> unknown parametrisation\n");
+      break;
+    default:
+        func=0;
+        printf("<AliGenGSIlib::GetY> unknown parametrisation\n");
+    }
+    return func;
+}
+
+
+
+GenFuncIp AliGenGSIlib::GetIp(Param_t param, const char * tname)
+{
+// Return pointer to particle type parameterisation
+    GenFuncIp func=0;
+    TString sname(tname);
+    switch (param) 
+    {
+    case upsilon_p:
+      if (sname=="MUON"){
+       func=IpUpsilonMUON;
+       break;
+      }
+      if (sname=="RITMAN"){
+       func=IpUpsilonRitman;
+       break;
+      }
+      if (sname=="KAREL"){
+       func=IpUpsilonKarel;
+       break;
+      }
+      func=0;
+      printf("<AliGenGSIlib::GetIP> unknown parametrisation\n");
+      break;
+    default:
+        func=0;
+        printf("<AliGenGSIlib::GetIp> unknown parametrisation\n");
+    }
+    return func;
+}
+
+
+
+
+
+
+
+
+
+
+
+
+

diff --git a/EVGEN/AliGenGSIlib.h b/EVGEN/AliGenGSIlib.h
new file mode 100644 (file)
index 0000000..437f078
--- /dev/null
+++ b/EVGEN/AliGenGSIlib.h
@@ -0,0 +1,46 @@
+#ifndef ALIGENGSILIB_H
+#define ALIGENGSILIB_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+
+/* $Id$ */
+
+#include <TROOT.h>
+#include "GenTypeDefs.h"
+
+class AliGenGSIlib :public TObject{
+ public:
+
+// Upsilon RITMAN   
+    static Double_t PtUpsilonRitman( Double_t *px, Double_t *dummy );
+    static Double_t YUpsilonRitman(Double_t *py, Double_t *dummy);
+    static Int_t    IpUpsilonRitman();
+// Upsilon Karel
+    static Double_t PtUpsilonKarel( Double_t *px, Double_t *dummy );
+    static Double_t YUpsilonKarel(Double_t *py, Double_t *dummy);
+    static Int_t    IpUpsilonKarel();
+// YpsMUON
+    static Double_t PtUpsilonMUON( Double_t *px, Double_t *dummy );
+    static Double_t YUpsilonMUON(Double_t *py, Double_t *dummy);
+    static Int_t    IpUpsilonMUON();
+
+//
+    typedef Double_t (*GenFunc)  (Double_t *, Double_t *);
+    typedef Int_t    (*GenFuncIp)();    
+    static GenFunc   GetPt(Param_t param,const char * tname=0);
+    static GenFunc   GetY(Param_t param,const char * tname=0);
+    static GenFuncIp GetIp(Param_t param,const char *tname=0);    
+    static void SetDebug(Bool_t debug){fDebug=debug;}
+private:
+    static Bool_t fgDebug;  // Debug flag 
+  ClassDef(AliGenGSIlib,1)
+};
+
+#endif
+
+
+
+
+
+
+

diff --git a/EVGEN/EVGENLinkDef.h b/EVGEN/EVGENLinkDef.h
index 49ef80da85074d9a235e87fd560e8a060197fa36..9be3856d6d2339909ba5ee4f039b5ccba7e9ae3d 100644 (file)
--- a/EVGEN/EVGENLinkDef.h
+++ b/EVGEN/EVGENLinkDef.h
@@ -30,6 +30,7 @@
 #pragma link C++ class  AliGenHalo;
 #pragma link C++ class  AliDimuCombinator;
 #pragma link C++ class  AliGenPHOSlib;
+#pragma link C++ class  AliGenGSIlib;
 #endif
 
 

diff --git a/EVGEN/Makefile b/EVGEN/Makefile
index 5939f65c4088184a1068e54622ff0ad5e93b60ff..1cf214f9ca623fe6374be0bb260dc6d5721f8cdf 100644 (file)
--- a/EVGEN/Makefile
+++ b/EVGEN/Makefile
@@ -14,7 +14,8 @@ SRCS          = AliGenHIJINGpara.cxx AliGenBox.cxx AliGenFixed.cxx \
                 AliGenMUONlib.cxx AliGenParam.cxx AliGenPythia.cxx \
                 AliPythia.cxx AliDimuCombinator.cxx AliGenHalo.cxx \
                 AliGenExtFile.cxx AliGenScan.cxx AliGenPHOSlib.cxx \
-               AliGenDoubleScan.cxx AliGenCocktailEntry.cxx
+               AliGenDoubleScan.cxx AliGenCocktailEntry.cxx \
+               AliGenGSIlib.cxx
 # C++ Headers
 
 HDRS          = $(SRCS:.cxx=.h) EVGENLinkDef.h