In DAs:

[u/mrichter/AliRoot.git] / macros / Config_gener.C

diff --git a/macros/Config_gener.C b/macros/Config_gener.C
index 1d556da42aaeaf1c99b2833016c96411a6fbc793..8820a33bbd3ba49d08835badd20d6cc95cdfff0c 100644 (file)
--- a/macros/Config_gener.C
+++ b/macros/Config_gener.C
@@ -28,11 +28,11 @@
 #include "STRUCT/AliMAG.h"
 #endif
 
-enum gentype_t {hijing, hijingParam, gun, box, pythia, herwig, 
+enum gentype_t {hijing, hijingParam, gun, box, pythia, herwig, mcatnlo, 
                param1, param2, param3, param4, 
                cocktail, fluka, halo, ntuple, scan};
 
-gentype_t gentype = herwig;
+gentype_t gentype = mcatnlo;
 
 Int_t ntracks=1;
 
@@ -128,7 +128,7 @@ void Config()
        gener->SetPhi(180.);
        gener->SetTheta(5.);
        gener->SetOrigin(0,0,0);        //vertex position
-       gener->SetPart(13);             //GEANT particle type
+       gener->SetPart(kNeutron);
        gGener = gener;
       }
       break;
@@ -144,7 +144,7 @@ void Config()
        gener->SetOrigin(0,0,0);   
        //vertex position
        gener->SetSigma(0,0,0);         //Sigma in (X,Y,Z) (cm) on IP position
-       gener->SetPart(5);              //GEANT particle type
+       gener->SetPart(kMuonPlus);
        gGener = gener;
       }
       break;
@@ -159,7 +159,7 @@ void Config()
        gener->SetThetaRange(9,9);
        //vertex position
        gener->SetSigma(6,6,0);         //Sigma in (X,Y,Z) (cm) on IP position
-       gener->SetPart(5); 
+       gener->SetPart(kMuonPlus); 
        gener->SetRange(20, -100, 100, 20, -100, 100, 1, 500, 500);
        gGener = gener;
       }
@@ -290,6 +290,29 @@ void Config()
       }
       break;              
 
+    case mcatnlo:
+    {
+      AliGenHerwig *gener = new AliGenHerwig(-1);
+      gener->SetMomentumRange(0,7000);
+      gener->SetPhiRange(0. ,360.);
+      gener->SetThetaRange(0., 180.);
+      gener->SetYRange(-10,10);
+      gener->SetPtRange(0,7000);
+      gener->SetOrigin(0,0,0);       // vertex position
+      gener->SetVertexSmear(kPerEvent);
+      gener->SetSigma(0,0,5.6);      // Sigma in (X,Y,Z) (cm) on IP position
+      gener->SetBeamMomenta(7000,7000);
+      gener->SetProjectile("P");
+      gener->SetTarget("P");
+      gener->SetStrucFunc(kCTEQ5M);
+      gener->SetProcess(-1705);
+      gener->SetHardProcessFile("sb.events");
+      gener->SetEventListRange(0,1);
+      gener->SetTrackingFlag(0);
+      gGener = gener;
+    }
+    break;  
+      
     case param1:
 //*******************************************************
 // Example for J/psi  Production from  Parameterisation