case kPyDirectGamma:
case kPyLhwgMb:
break;
- case kPyWPWHG: /// !!!!!!!!!! Change done for W prod with POWHEG !!!!!!!!!! :)
+ case kPyWPWHG:
case kPyW:
case kPyZ:
+ case kPyZgamma:
case kPyMBRSingleDiffraction:
case kPyMBRDoubleDiffraction:
case kPyMBRCentralDiffraction:
fPythia->SetMSTJ(44, 2); // option to run alpha_s
fPythia->SetPARJ(82, 1.); // Cut off for parton showers
}
+
+ StdoutToAliDebug(1, fPythia->Pystat(4));
+ StdoutToAliDebug(1, fPythia->Pystat(5));
}
void AliGenPythia::Generate()
fProcess != kPyMbMSEL1 &&
fProcess != kPyW &&
fProcess != kPyZ &&
+ fProcess != kPyZgamma &&
fProcess != kPyCharmppMNRwmi &&
fProcess != kPyBeautyppMNRwmi &&
fProcess != kPyBeautyJets &&
- fProcess != kPyWPWHG && /// !!!!!!!!!!!!!!!!! Change done for W with POHWEG !!!!!!!!!!!!!!!!!!! :)
+ fProcess != kPyWPWHG &&
fProcess != kPyJetsPWHG &&
fProcess != kPyCharmPWHG &&
fProcess != kPyBeautyPWHG) {
//Introducing child cuts in case kPyW, kPyZ, kPyMb, and kPyMbNonDiff
if ( (
- fProcess == kPyWPWHG || /// !!!!!!!!!!!!!!!! Added for W with POWHEG !!!!!!!!!! :)
+ fProcess == kPyWPWHG ||
fProcess == kPyW ||
fProcess == kPyZ ||
+ fProcess == kPyZgamma ||
fProcess == kPyMbDefault ||
fProcess == kPyMb ||
fProcess == kPyMbAtlasTuneMC09 ||
case kPyDirectGamma:
case kPyLhwgMb:
break;
- case kPyWPWHG: ///// !!!!!!!!!!!!!!!!!!!1 Change for W with POWEG !!!!!!!!!!!!!!!!!1 :)
+ case kPyWPWHG:
case kPyW:
case kPyZ:
+ case kPyZgamma:
case kPyMBRSingleDiffraction:
case kPyMBRDoubleDiffraction:
case kPyMBRCentralDiffraction:
fProcess != kPyMbMSEL1 &&
fProcess != kPyW &&
fProcess != kPyZ &&
+ fProcess != kPyZgamma &&
fProcess != kPyCharmppMNRwmi &&
fProcess != kPyBeautyppMNRwmi &&
- fProcess != kPyWPWHG && /// !!!!!!!!!!!!!!!!!!! Change for W with POWHEG !!!!!!!!!!!!!!!!1 :)
+ fProcess != kPyWPWHG &&
fProcess != kPyJetsPWHG &&
fProcess != kPyCharmPWHG &&
fProcess != kPyBeautyPWHG) {
//Introducing child cuts in case kPyW, kPyZ, kPyMb, and kPyMbNonDiff
if ( (
- fProcess == kPyWPWHG || // !!!!!!!!!!!!! Change for W with POWEG !!!!!!!!!!!!!!!1 :)
+ fProcess == kPyWPWHG ||
fProcess == kPyW ||
fProcess == kPyZ ||
+ fProcess == kPyZgamma ||
fProcess == kPyMbDefault ||
fProcess == kPyMb ||
fProcess == kPyMbWithDirectPhoton ||
SetMSTP(71,1); //Final QCD & QED showers on
break;
+ case kPyZgamma:
+ //Inclusive production of Z
+ SetMSEL(0);
+ //f fbar -> Z/gamma
+ SetMSUB(1,1);
+ // Initial/final parton shower on (Pythia default)
+ // With parton showers on we are generating "Z inclusive process"
+ SetMSTP(61,1); //Initial QCD & QED showers on
+ SetMSTP(71,1); //Final QCD & QED showers on
+ break;
case kPyMBRSingleDiffraction:
case kPyMBRDoubleDiffraction:
case kPyMBRCentralDiffraction:
case kPyCharmPWHG:
case kPyBeautyPWHG:
- case kPyWPWHG: //!!!! CHANGE done for W with Powheg !!!!!!! :)
+ case kPyWPWHG:
// number of warnings printed on the shell
SetMSTU(26,20);
SetMSTP(61,1); //Initial QCD & QED showers on
SetMSTP(71,1); //Final QCD & QED showers on
break;
+ case kPyZgamma:
+ //Inclusive production of Z
+ SetMSEL(0);
+ //f fbar -> Z/gamma
+ SetMSUB(1,1);
+ // Initial/final parton shower on (Pythia default)
+ // With parton showers on we are generating "Z inclusive process"
+ SetMSTP(61,1); //Initial QCD & QED showers on
+ SetMSTP(71,1); //Final QCD & QED showers on
+ break;
case kPyMBRSingleDiffraction:
case kPyMBRDoubleDiffraction:
case kPyMBRCentralDiffraction:
break;
case kPyCharmPWHG:
case kPyBeautyPWHG:
+ case kPyWPWHG:
// number of warnings printed on the shell
SetMSTU(26,20);
kPyCharmppMNR, kPyCharmppMNRwmi, kPyD0ppMNR, kPyDPlusppMNR, kPyDPlusStrangeppMNR,
kPyBeautyppMNR, kPyBeautyppMNRwmi, kPyBeautyJets, kPyW, kPyZ, kPyLambdacppMNR, kPyMbMSEL1,
kPyOldUEQ2ordered, kPyOldUEQ2ordered2, kPyOldPopcorn,
- kPyLhwgMb, kPyMbDefault, kPyMbAtlasTuneMC09, kPyMBRSingleDiffraction, kPyMBRDoubleDiffraction, kPyMBRCentralDiffraction, kPyJetsPWHG, kPyCharmPWHG, kPyBeautyPWHG, kPyWPWHG
+ kPyLhwgMb, kPyMbDefault, kPyMbAtlasTuneMC09, kPyMBRSingleDiffraction, kPyMBRDoubleDiffraction, kPyMBRCentralDiffraction, kPyJetsPWHG, kPyCharmPWHG, kPyBeautyPWHG, kPyWPWHG, kPyZgamma
}
Process_t;
// With parton showers on we are generating "Z inclusive process"
ReadString("PartonLevel:ISR = on");
ReadString("PartonLevel:FSR = on");
- break;
+ break;
+ case kPyZgamma:
+ //Inclusive production of Z/gamma*
+ //f fbar -> Z/gamma
+ ReadString("WeakSingleBoson:ffbar2gmZ = on");
+ // Initial/final parton shower on (Pythia default)
+ // With parton showers on we are generating "Z inclusive process"
+ ReadString("PartonLevel:ISR = on");
+ ReadString("PartonLevel:FSR = on");
+ break;
case kPyMBRSingleDiffraction:
ReadString("Diffraction:PomFlux = 5");
ReadString("SoftQCD:singleDiffractive = on");