From 589380c6dbe34deac8523a5fbcb62addf6ebca33 Mon Sep 17 00:00:00 2001
From: morsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Fri, 4 Mar 2005 07:55:51 +0000
Subject: [PATCH] Option kPyW for W-production added. (Zaida Conesa del Valle
 <zconesa@subatech.in2p3.fr>)

---
 PYTHIA6/AliGenPythia.cxx |  6 +++++-
 PYTHIA6/AliPythia.cxx    | 23 +++++++++++++++++++++++
 PYTHIA6/AliPythia.h      |  3 ++-
 3 files changed, 30 insertions(+), 2 deletions(-)

diff --git a/PYTHIA6/AliGenPythia.cxx b/PYTHIA6/AliGenPythia.cxx
index 2683ca5e4e4..afb5b8540ba 100644
--- 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 43ab97e06d8..7dc79c6737d 100644
--- 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 03aae7f95dd..c71ffd6fee7 100644
--- 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
 
 
 
+
-- 
2.43.0