case kPyJets:
case kPyDirectGamma:
break;
+ case kPyW:
+ break;
}
//
//
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);
#endif
+
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
}
+
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
+