Introduce and adapt parameters ala arXiv:1101.2231v1
authorloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Jan 2011 13:16:33 +0000 (13:16 +0000)
committerloizides <loizides@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 17 Jan 2011 13:16:33 +0000 (13:16 +0000)
TAmpt/AMPT/amptsetdef.f
TAmpt/AliGenAmpt.cxx
TAmpt/AliGenAmpt.h
TAmpt/TAmpt.cxx
TAmpt/TAmpt.h
TAmpt/macros/createAmptMC.C
TAmpt/macros/rootlogon.C

index a31a5ce..578e1bb 100644 (file)
@@ -45,8 +45,8 @@ c     read initialization value for hadron cascade:
       NTMAX    = 150           ! NTMAX: number of timesteps (D=150), see below
       DT       = 0.2           ! DT: timestep in fm (hadron cascade time= DT*NTMAX) (D=0.2)
 c     parj(41) and (42) are a and b parameters in Lund string fragmentation:
-      PARJ(41) = 2.2           ! PARJ(41): parameter a in Lund symmetric splitting function
-      PARJ(42) = 0.5           ! PARJ(42): parameter b in Lund symmetric splitting function
+      PARJ(41) = 0.5           ! PARJ(41): parameter a in Lund symmetric splitting function
+      PARJ(42) = 0.9           ! PARJ(42): parameter b in Lund symmetric splitting function
 c     IHPR2(11)=3 (or 2) allows the popcorn mechanism in PYTHIA and 
 c     increase the net-baryon stopping in rapidity (value HIJING is 1):
       ipop      = 1            ! (D=1,yes;0,no) flag for popcorn mechanism(netbaryon stopping)
@@ -63,7 +63,7 @@ c     Minimum pt of hard or semihard scatterings in HIJING: D=2.0 GeV.
 c     read initialization value for parton cascade:
       xmu      = 3.2264d0      ! parton screening mass in fm^(-1) (D=3.2264d0), see below
       izpc     = 0             ! IZPC: (D=0 forward-angle parton scatterings; 100,isotropic)
-      alpha    = 0.47140452d0  ! alpha in parton cascade
+      alpha    = 0.3333d0      ! alpha in parton cascade
 c     quark coalescence radii in momentum and space for string melting:
       dpcoal   = 1d6           ! dpcoal in GeV
       drcoal   = 1d6           ! drcoal in fm
index 68d59d5..20b5083 100644 (file)
@@ -76,6 +76,9 @@ AliGenAmpt::AliGenAmpt()
     fNtMax(150),
     fIpop(1),
     fXmu(3.2264),
+    fAlpha(1./3),
+    fStringA(0.5),
+    fStringB(0.9),
     fHeader(new AliGenAmptEventHeader("Ampt"))
 {
   // Constructor
@@ -125,6 +128,9 @@ AliGenAmpt::AliGenAmpt(Int_t npart)
     fNtMax(150),
     fIpop(1),
     fXmu(3.2264),
+    fAlpha(1./3),
+    fStringA(0.5),
+    fStringB(0.9),
     fHeader(new AliGenAmptEventHeader("Ampt"))
 {
   // Default PbPb collisions at 5.5 TeV
@@ -207,6 +213,8 @@ void AliGenAmpt::Init()
   fAmpt->SetNtMax(fNtMax);
   fAmpt->SetIpop(fIpop);
   fAmpt->SetXmu(fXmu);
+  fAmpt->SetAlpha(fAlpha);
+  fAmpt->SetStringFrag(fStringA, fStringB);
 
   AliGenMC::Init();
     
index d3cd7f7..d160548 100644 (file)
@@ -53,10 +53,12 @@ class AliGenAmpt : public AliGenMC
     virtual void    SetBoostLHC(Int_t flag = 0)       {fLHC        = flag;}
     virtual void    SetRandomPz(Bool_t flag = 0)      {fRandomPz   = flag;}
     virtual void    SwitchOffHeavyQuarks(Bool_t flag = kTRUE) {fNoHeavyQuarks = flag;}
-    virtual void    SetIsoft(Int_t i)                 {fIsoft = i;  }
-    virtual void    SetNtMax(Int_t max)               {fNtMax = max;}
-    virtual void    SetIpop(Int_t pop)                {fIpop  = pop;}
-    virtual void    SetXmu(Float_t m)                 {fXmu   = m;  }
+    virtual void    SetIsoft(Int_t i)                         {fIsoft = i;  }
+    virtual void    SetNtMax(Int_t max)                       {fNtMax = max;}
+    virtual void    SetIpop(Int_t pop)                        {fIpop  = pop;}
+    virtual void    SetXmu(Float_t m)                         {fXmu   = m;  }
+    virtual void    SetAlpha(Float_t alpha)                   {fAlpha = alpha;            }
+    virtual void    SetStringFrag(Float_t a, Float_t b)       {fStringA = a; fStringB = b;}
            
     // Getters
     virtual TString GetReferenceFrame()  const {return fFrame;}
@@ -128,6 +130,9 @@ class AliGenAmpt : public AliGenMC
     Int_t                    fNtMax;           // NTMAX: number of timesteps (D=150)
     Int_t                    fIpop;            // (D=1,yes;0,no) flag for popcorn mechanism(netbaryon stopping)
     Float_t                  fXmu;             // parton screening mass in fm^(-1) (D=3.2264d0)
+    Float_t                  fAlpha;           // alpha running (fixed) coupling
+    Float_t                  fStringA;         // string frag parameter A
+    Float_t                  fStringB;         // string frag parameter B
     AliGenHijingEventHeader *fHeader;          // header
 
   private:
@@ -141,6 +146,6 @@ class AliGenAmpt : public AliGenMC
     // check if stable
     Bool_t Stable(TParticle*  particle) const;
 
-    ClassDef(AliGenAmpt, 1) // AliGenerator interface to Ampt
+    ClassDef(AliGenAmpt, 2) // AliGenerator interface to Ampt
 };
 #endif
index 96067e1..e45d6df 100644 (file)
@@ -1239,6 +1239,21 @@ void TAmpt::SetIpop(Int_t pop)
 //______________________________________________________________________________
 void TAmpt::SetXmu(Float_t m)     
 {
-  // set parton screening mass in fm^
+  // set parton screening mass in fm^-1.
   PARA2.xmu     = m;  
 }
+
+//______________________________________________________________________________
+void TAmpt::SetAlpha(Float_t alpha)     
+{
+  // set parton screening mass in fm^-1.
+  PARA2.alpha     = alpha;  
+}
+
+//______________________________________________________________________________
+void TAmpt::SetStringFrag(Float_t a, Float_t b)
+{
+  // Set string frag parameters, f(z) = 1/z*(1-z)^a*exp(-b*mz^2/z).
+   SetPARJ(41, a);
+   SetPARJ(42, b);
+}
index 92ce1ad..5281488 100644 (file)
@@ -132,10 +132,12 @@ public:
   virtual void        Rluset(Int_t lfn, Int_t move=0);   
 
   // set AMPT specific
-  virtual void       SetIsoft(Int_t i);
-  virtual void       SetNtMax(Int_t max);
-  virtual void       SetIpop(Int_t pop);
-  virtual void       SetXmu(Float_t m);
+  virtual void        SetIsoft(Int_t i);
+  virtual void        SetNtMax(Int_t max);
+  virtual void        SetIpop(Int_t pop);
+  virtual void        SetXmu(Float_t m);
+  virtual void        SetAlpha(Float_t alpha);
+  virtual void        SetStringFrag(Float_t a, Float_t b);
 
 protected:
   Double_t     fEfrm;     // Energy in the centre of mass (CMS) or lab-frame (LAB)
index 7bb1143..3159630 100644 (file)
@@ -119,8 +119,12 @@ void createAmptMC(Int_t nEvents,
   genHi->SetDecaysOff(1);    // neutral pion and heavy particle decays switched off
   genHi->SetSpectators(1);   // track spectators 
   genHi->SetIsoft(1);        // standard AMPT
-  genHi->SetNtMax(3);        // time bins
   genHi->SetXmu(3.2264);     // parton xsection
+  genHi->SetNtMax(3);        // time bins
+  if (0) { //RHIC settings
+    genHi->SetAlpha(0.47140452);
+    genHi->SetStringFrag(2.2,0.5);
+  }
   genHi->Init();
 
   TClonesArray *inucs = new TClonesArray("TParticle",10000);
index c69371b..d11b236 100644 (file)
@@ -22,8 +22,7 @@
   gSystem->Load("libpythia6.so");
   gSystem->Load("libAliPythia6.so");
   gSystem->AddIncludePath("-I$ALICE_ROOT/include");
-  gSystem->AddIncludePath("-I$ALICE_ROOT_ORIG/include");
-  gSystem->AddIncludePath("-I$ALICE_ROOT_ORIG/PYTHIA6");
+  gSystem->AddIncludePath("-I$ALICE_ROOT/PYTHIA6");
   if (gSystem->Getenv("TMPDIR")) 
     gSystem->SetBuildDir(gSystem->Getenv("TMPDIR"));
   if (!gAlice)