AtlasTuning();
break;
+
+ case kPyMbWithDirectPhoton:
+// Minimum Bias pp-Collisions with direct photon processes added
+//
+//
+// select Pythia min. bias model
+ SetMSEL(0);
+ 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
+
+ SetMSUB(14,1); //
+ SetMSUB(18,1); //
+ SetMSUB(29,1); //
+ SetMSUB(114,1); //
+ SetMSUB(115,1); //
+
+
+ AtlasTuning();
+ break;
+
case kPyMbDefault:
// Minimum Bias pp-Collisions
//
// (2) The nuclear geometry using the Glauber Model
//
- fGlauber = new AliFastGlauber();
+ fGlauber = AliFastGlauber::Instance();
fGlauber->Init(2);
fGlauber->SetCentralityClass(cMin, cMax);
//
// Isotropic decay ????
Double_t cost = 2. * gRandom->Rndm() - 1.;
- Double_t sint = TMath::Sqrt(1. - cost * cost);
+ Double_t sint = TMath::Sqrt((1.-cost)*(1.+cost));
Double_t phis = 2. * TMath::Pi() * gRandom->Rndm();
Double_t pz1 = pst * cost;