Int_t nTkbles = 0; // Trackable particles
if (fProcess != kPyMbDefault &&
fProcess != kPyMb &&
+ fProcess != kPyMbAtlasTuneMC09 &&
fProcess != kPyMbWithDirectPhoton &&
fProcess != kPyJets &&
fProcess != kPyDirectGamma &&
Int_t pdg = 0;
if (fFragPhotonInCalo) pdg = 22 ; // Photon
- else if (fPi0InCalo) pdg = 111 ; // Pi0
+ else if (fPi0InCalo) pdg = 111 ; // Pi0
for (i=0; i< np; i++) {
TParticle* iparticle = (TParticle *) fParticles.At(i);
return 0;
AliDebug(5,Form("Found an electron jet (pt,eta,phi) = (%f,%f,%f)",pt,eta,phi));
}
-
// Check for minimum multiplicity
if (fTriggerMultiplicity > 0) {
Int_t multiplicity = 0;
delete [] pParent;
return 0;
}
-
Printf("Triggered on event with multiplicity of %d >= %d", multiplicity, fTriggerMultiplicity);
}
fProcess == kPyZ ||
fProcess == kPyMbDefault ||
fProcess == kPyMb ||
+ fProcess == kPyMbAtlasTuneMC09 ||
fProcess == kPyMbWithDirectPhoton ||
fProcess == kPyMbNonDiffr)
&& (fCutOnChild == 1) ) {
//
// select Pythia min. bias model
SetMSEL(0);
- SetMSUB(92,0); // single diffraction AB-->XB
- SetMSUB(93,0); // single diffraction AB-->AX
+ SetMSUB(92,1); // single diffraction AB-->XB
+ SetMSUB(93,1); // single diffraction AB-->AX
SetMSUB(94,1); // double diffraction
SetMSUB(95,1); // low pt production
-
- SetPARP(67,2.5); // Regulates Initial State Radiation (value from best fit to D0 dijet analysis)
- SetMSTP(82,4); // Double Gaussian Model
- SetPARP(82,2.0); // [GeV] PT_min at Ref. energy
- SetPARP(84,0.4); // Core radius
- SetPARP(85,0.90) ; // Regulates gluon prod. mechanism
- SetPARP(86,0.95); // Regulates gluon prod. mechanism
- SetPARP(89,1800.); // [GeV] Ref. energy
- SetPARP(90,0.25); // 2*epsilon (exponent in power law)
-
break;
case kPyLhwgMb:
// Les Houches Working Group 05 Minimum Bias pp-Collisions: hep-ph/0604120