Option kPyW for W-production added.

(Zaida Conesa del Valle <zconesa@subatech.in2p3.fr>)

diff --git a/PYTHIA6/AliGenPythia.cxx b/PYTHIA6/AliGenPythia.cxx
index 2683ca5e4e4909898350cd7d5029c9f13743e30d..afb5b8540ba178b394a9a48edcc49e7296fe0ee6 100644 (file)
--- a/PYTHIA6/AliGenPythia.cxx
+++ b/PYTHIA6/AliGenPythia.cxx
@@ -263,6 +263,8 @@ void AliGenPythia::Init()
     case kPyJets:
     case kPyDirectGamma:
        break;
+    case kPyW:
+        break;
     }
 //
 //
@@ -396,7 +398,8 @@ void AliGenPythia::Generate()
        Int_t nTkbles = 0;   // Trackable particles
        if (fProcess != kPyMb && fProcess != kPyJets && 
            fProcess != kPyDirectGamma &&
-           fProcess != kPyMbNonDiffr) {
+           fProcess != kPyMbNonDiffr  &&
+           fProcess != kPyW) {
            
            for (i = 0; i < np; i++) {
                TParticle* iparticle = (TParticle *) fParticles->At(i);
@@ -1014,3 +1017,4 @@ void AliGenPythia::Streamer(TBuffer &R__b)
 #endif
 
 
+

diff --git a/PYTHIA6/AliPythia.cxx b/PYTHIA6/AliPythia.cxx
index 43ab97e06d8af9e7a3ec3a34751cfb2baf379881..7dc79c6737de62465141c9c19e7f80a7428779ff 100644 (file)
--- a/PYTHIA6/AliPythia.cxx
+++ b/PYTHIA6/AliPythia.cxx
@@ -549,6 +549,28 @@ void AliPythia::ProcInit(Process_t process, Float_t energy, StrucFunc_t strucfun
       SetPMAS(5,1,4.75);
 
       break;
+ 
+    case kPyW:
+
+      //Inclusive production of W+/-
+      SetMSEL(0);
+      //f fbar -> W+ 
+      SetMSUB(2,1);
+      //       //f fbar -> g W+
+      //       SetMSUB(16,1);
+      //       //f fbar -> gamma W+
+      //       SetMSUB(20,1);
+      //       //f g -> f W+  
+      //       SetMSUB(31,1);
+      //       //f gamma -> f W+
+      //       SetMSUB(36,1);
+      
+      // Initial/final parton shower on (Pythia default)
+      // With parton showers on we are generating "W inclusive process"
+      SetMSTP(61,1); //Initial QCD & QED showers on
+      SetMSTP(71,1); //Final QCD & QED showers on
+      
+      break;  
     }
 //
 //  Initialize PYTHIA
@@ -1223,3 +1245,4 @@ void AliPythia::GetQuenchingParameters(Double_t& xp, Double_t& yp, Double_t z[4]
 
 }
 
+

diff --git a/PYTHIA6/AliPythia.h b/PYTHIA6/AliPythia.h
index 03aae7f95ddf04f84a99578f9b1f89b3539bc04c..c71ffd6fee76271b62984e674b4321be3a98cc2d 100644 (file)
--- a/PYTHIA6/AliPythia.h
+++ b/PYTHIA6/AliPythia.h
@@ -13,7 +13,7 @@ typedef enum
  kPyJpsi, kPyJpsiChi, kPyMb, kPyMbNonDiffr, kPyJets, kPyDirectGamma,
  kPyCharmPbPbMNR, kPyD0PbPbMNR, kPyDPlusPbPbMNR, kPyBeautyPbPbMNR,
  kPyCharmpPbMNR, kPyD0pPbMNR, kPyDPluspPbMNR, kPyBeautypPbMNR,
- kPyCharmppMNR, kPyD0ppMNR, kPyDPlusppMNR, kPyBeautyppMNR}
+ kPyCharmppMNR, kPyD0ppMNR, kPyDPlusppMNR, kPyBeautyppMNR, kPyW}
 Process_t;
 /*
 typedef enum
@@ -87,3 +87,4 @@ class AliPythia : public TPythia6, public AliRndm
 
 
 
+