fGinit(1),
fGfinal(1),
fHadronisation(1),
+ fPatchOmegaDalitz(0),
fNpartons(0),
fReadFromFile(0),
fQuench(0),
fTriggerEta(0.9),
fTriggerMultiplicity(0),
fTriggerMultiplicityEta(0),
+ fTriggerMultiplicityPtMin(0),
fCountMode(kCountAll),
fHeader(0),
fRL(0),
- fFileName(0),
+ fkFileName(0),
fFragPhotonInCalo(kFALSE),
fPi0InCalo(kFALSE) ,
fPhotonInCalo(kFALSE),
{
// Default Constructor
fEnergyCMS = 5500.;
- SetNuclei(0,0);
if (!AliPythiaRndm::GetPythiaRandom())
AliPythiaRndm::SetPythiaRandom(GetRandom());
}
fGinit(kTRUE),
fGfinal(kTRUE),
fHadronisation(kTRUE),
+ fPatchOmegaDalitz(0),
fNpartons(0),
fReadFromFile(kFALSE),
fQuench(kFALSE),
fTriggerEta(0.9),
fTriggerMultiplicity(0),
fTriggerMultiplicityEta(0),
+ fTriggerMultiplicityPtMin(0),
fCountMode(kCountAll),
fHeader(0),
fRL(0),
- fFileName(0),
+ fkFileName(0),
fFragPhotonInCalo(kFALSE),
fPi0InCalo(kFALSE) ,
fPhotonInCalo(kFALSE),
// Set random number generator
if (!AliPythiaRndm::GetPythiaRandom())
AliPythiaRndm::SetPythiaRandom(GetRandom());
- SetNuclei(0,0);
}
AliGenPythia::~AliGenPythia()
if (fReadFromFile) {
- fRL = AliRunLoader::Open(fFileName, "Partons");
+ fRL = AliRunLoader::Open(fkFileName, "Partons");
fRL->LoadKinematics();
fRL->LoadHeader();
} else {
fParentSelect[0] = 431;
fFlavorSelect = 4;
break;
+ case kPyLambdacppMNR:
+ fParentSelect[0] = 4122;
+ fFlavorSelect = 4;
+ break;
case kPyBeauty:
case kPyBeautyJets:
case kPyBeautyPbPbMNR:
fParentSelect[0] = 443;
break;
case kPyMbDefault:
+ case kPyMbAtlasTuneMC09:
case kPyMb:
case kPyMbWithDirectPhoton:
case kPyMbNonDiffr:
AliFastGlauber::Instance()->GetRandomBHard(bimp);
fPythia->SetPARJ(197, bimp);
fImpact = bimp;
+ fPythia->Qpygin0();
}
//
// Either produce new event or read partons from file
//
// .. and perform hadronisation
// printf("Calling hadronisation %d\n", fPythia->GetN());
- fPythia->Pyexec();
+
+ if (fPatchOmegaDalitz) {
+ fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 0);
+ fPythia->Pyexec();
+ fPythia->DalitzDecays();
+ fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1);
+ }
+ fPythia->Pyexec();
}
fTrials++;
fPythia->ImportParticles(&fParticles,"All");
- Boost();
+
+ if (TMath::Abs(fDyBoost) > 1.e-4) Boost();
+ if(TMath::Abs(fXingAngleY) > 1.e-10) BeamCrossAngle();
+
//
//
//
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;
Int_t statusCode = iparticle->GetStatusCode();
// Initial state particle
- if (statusCode > 20)
- continue;
-
- // skip quarks and gluons
- Int_t pdgCode = TMath::Abs(iparticle->GetPdgCode());
- if (pdgCode <= 10 || pdgCode == 21)
+ if (statusCode != 1)
continue;
-
+ // eta cut
if (fTriggerMultiplicityEta > 0 && TMath::Abs(iparticle->Eta()) > fTriggerMultiplicityEta)
continue;
-
+ // pt cut
+ if (iparticle->Pt() < fTriggerMultiplicityPtMin)
+ continue;
+
TParticlePDG* pdgPart = iparticle->GetPDG();
if (pdgPart && pdgPart->Charge() == 0)
continue;
delete [] pParent;
return 0;
}
-
- Printf("Triggered on event with multiplicity of %d > %d", multiplicity, fTriggerMultiplicity);
+ Printf("Triggered on event with multiplicity of %d >= %d", multiplicity, fTriggerMultiplicity);
}
// Select events with a photon pt > min pt going to PHOS eta acceptance or exactly PHOS eta phi
fProcess == kPyZ ||
fProcess == kPyMbDefault ||
fProcess == kPyMb ||
+ fProcess == kPyMbAtlasTuneMC09 ||
fProcess == kPyMbWithDirectPhoton ||
fProcess == kPyMbNonDiffr)
&& (fCutOnChild == 1) ) {
//
// Event Vertex
fHeader->SetPrimaryVertex(fVertex);
-
+ fHeader->SetInteractionTime(fEventTime);
//
// Number of primaries
fHeader->SetNProduced(fNprimaries);
//loop for all particles and produce the phi rotation
Int_t np = (fHadronisation) ? fParticles.GetEntriesFast() : fNpartons;
Double_t oldphi, newphi;
- Double_t newVx, newVy, R, Vz, time;
- Double_t newPx, newPy, pt, Pz, e;
+ Double_t newVx, newVy, r, vZ, time;
+ Double_t newPx, newPy, pt, pz, e;
for(Int_t i=0; i< np; i++) {
TParticle* iparticle = (TParticle *) fParticles.At(i);
oldphi = iparticle->Phi();
if(newphi < 0) newphi = 2*TMath::Pi() + newphi; // correct angle
if(newphi > 2*TMath::Pi()) newphi = newphi - 2*TMath::Pi(); // correct angle
- R = iparticle->R();
- newVx = R*TMath::Cos(newphi);
- newVy = R*TMath::Sin(newphi);
- Vz = iparticle->Vz(); // don't transform
+ r = iparticle->R();
+ newVx = r * TMath::Cos(newphi);
+ newVy = r * TMath::Sin(newphi);
+ vZ = iparticle->Vz(); // don't transform
time = iparticle->T(); // don't transform
pt = iparticle->Pt();
- newPx = pt*TMath::Cos(newphi);
- newPy = pt*TMath::Sin(newphi);
- Pz = iparticle->Pz(); // don't transform
- e = iparticle->Energy(); // don't transform
+ newPx = pt * TMath::Cos(newphi);
+ newPy = pt * TMath::Sin(newphi);
+ pz = iparticle->Pz(); // don't transform
+ e = iparticle->Energy(); // don't transform
// apply rotation
- iparticle->SetProductionVertex(newVx, newVy, Vz, time);
- iparticle->SetMomentum(newPx, newPy, Pz, e);
+ iparticle->SetProductionVertex(newVx, newVy, vZ, time);
+ iparticle->SetMomentum(newPx, newPy, pz, e);
} //end particle loop
}
-#ifdef never
-void AliGenPythia::Streamer(TBuffer &R__b)
-{
- // Stream an object of class AliGenPythia.
-
- if (R__b.IsReading()) {
- Version_t R__v = R__b.ReadVersion(); if (R__v) { }
- AliGenerator::Streamer(R__b);
- R__b >> (Int_t&)fProcess;
- R__b >> (Int_t&)fStrucFunc;
- R__b >> (Int_t&)fForceDecay;
- R__b >> fEnergyCMS;
- R__b >> fKineBias;
- R__b >> fTrials;
- fParentSelect.Streamer(R__b);
- fChildSelect.Streamer(R__b);
- R__b >> fXsection;
-// (AliPythia::Instance())->Streamer(R__b);
- R__b >> fPtHardMin;
- R__b >> fPtHardMax;
-// if (fDecayer) fDecayer->Streamer(R__b);
- } else {
- R__b.WriteVersion(AliGenPythia::IsA());
- AliGenerator::Streamer(R__b);
- R__b << (Int_t)fProcess;
- R__b << (Int_t)fStrucFunc;
- R__b << (Int_t)fForceDecay;
- R__b << fEnergyCMS;
- R__b << fKineBias;
- R__b << fTrials;
- fParentSelect.Streamer(R__b);
- fChildSelect.Streamer(R__b);
- R__b << fXsection;
-// R__b << fPythia;
- R__b << fPtHardMin;
- R__b << fPtHardMax;
- // fDecayer->Streamer(R__b);
- }
-}
-#endif
-
-