AliGenPythia::AliGenPythia():
AliGenMC(),
fProcess(kPyCharm),
+ fItune(-1),
fStrucFunc(kCTEQ5L),
fKineBias(0.),
fTrials(0),
fQuench(0),
fQhat(0.),
fLength(0.),
+ fImpact(0.),
fPtKick(1.),
fFullEvent(kTRUE),
fDecayer(new AliDecayerPythia()),
fTriggerEta(0.9),
fTriggerMultiplicity(0),
fTriggerMultiplicityEta(0),
+ fTriggerMultiplicityPtMin(0),
fCountMode(kCountAll),
fHeader(0),
fRL(0),
fFragPhotonInCalo(kFALSE),
fPi0InCalo(kFALSE) ,
fPhotonInCalo(kFALSE),
+ fEleInEMCAL(kFALSE),
fCheckEMCAL(kFALSE),
fCheckPHOS(kFALSE),
fCheckPHOSeta(kFALSE),
fFragPhotonOrPi0MinPt(0),
fPhotonMinPt(0),
+ fElectronMinPt(0),
fPHOSMinPhi(219.),
fPHOSMaxPhi(321.),
fPHOSEta(0.13),
{
// Default Constructor
fEnergyCMS = 5500.;
- SetNuclei(0,0);
if (!AliPythiaRndm::GetPythiaRandom())
AliPythiaRndm::SetPythiaRandom(GetRandom());
}
AliGenPythia::AliGenPythia(Int_t npart)
:AliGenMC(npart),
fProcess(kPyCharm),
+ fItune(-1),
fStrucFunc(kCTEQ5L),
fKineBias(0.),
fTrials(0),
fQuench(kFALSE),
fQhat(0.),
fLength(0.),
+ fImpact(0.),
fPtKick(1.),
fFullEvent(kTRUE),
fDecayer(new AliDecayerPythia()),
fTriggerEta(0.9),
fTriggerMultiplicity(0),
fTriggerMultiplicityEta(0),
+ fTriggerMultiplicityPtMin(0),
fCountMode(kCountAll),
fHeader(0),
fRL(0),
fFragPhotonInCalo(kFALSE),
fPi0InCalo(kFALSE) ,
fPhotonInCalo(kFALSE),
+ fEleInEMCAL(kFALSE),
fCheckEMCAL(kFALSE),
fCheckPHOS(kFALSE),
fCheckPHOSeta(kFALSE),
fFragPhotonOrPi0MinPt(0),
fPhotonMinPt(0),
+ fElectronMinPt(0),
fPHOSMinPhi(219.),
fPHOSMaxPhi(321.),
fPHOSEta(0.13),
// Set random number generator
if (!AliPythiaRndm::GetPythiaRandom())
AliPythiaRndm::SetPythiaRandom(GetRandom());
- SetNuclei(0,0);
}
AliGenPythia::~AliGenPythia()
fRL = 0x0;
}
//
- fPythia->ProcInit(fProcess,fEnergyCMS,fStrucFunc);
+ fPythia->ProcInit(fProcess,fEnergyCMS,fStrucFunc, fItune);
// Forward Paramters to the AliPythia object
fDecayer->SetForceDecay(fForceDecay);
// Switch off Heavy Flavors on request
fFlavorSelect = 4;
break;
case kPyBeauty:
+ case kPyBeautyJets:
case kPyBeautyPbPbMNR:
case kPyBeautypPbMNR:
case kPyBeautyppMNR:
fParentSelect[0] = 443;
break;
case kPyMbDefault:
+ case kPyMbAtlasTuneMC09:
case kPyMb:
case kPyMbWithDirectPhoton:
case kPyMbNonDiffr:
glauber->Init(2);
//read and store transverse almonds corresponding to differnt
//impact parameters.
- glauber->StoreAlmonds();
glauber->SetCentralityClass(0.,0.1);
fPythia->SetPARJ(200, 1.);
fPythia->SetPARJ(198, fQhat);
fPythia->SetPARJ(199, fLength);
-
- fPythia->SetMSTJ(41, 1); // QCD radiation only
fPythia->SetMSTJ(42, 2); // angular ordering
fPythia->SetMSTJ(44, 2); // option to run alpha_s
- //fPythia->SetMSTJ(47, 0); // No correction back to hard scattering element
- //fPythia->SetMSTJ(50, 0); // No coherence in first branching
fPythia->SetPARJ(82, 1.); // Cut off for parton showers
- // MSTJ(41) must NOT be 11 or 12, as then FSR may go through PYPTFS
- // (kt-ordered cascade) in which medium effects have not been introduced.
}
}
Double_t bimp;
// Quenching comes through medium-modified splitting functions.
AliFastGlauber::Instance()->GetRandomBHard(bimp);
- fPythia->SetPARJ(197,bimp);
+ fPythia->SetPARJ(197, bimp);
+ fImpact = bimp;
}
//
// Either produce new event or read partons from file
}
fTrials++;
fPythia->ImportParticles(&fParticles,"All");
- Boost();
+ if (TMath::Abs(fDyBoost) > 1.e-4) Boost();
//
//
//
fProcess != kPyW &&
fProcess != kPyZ &&
fProcess != kPyCharmppMNRwmi &&
- fProcess != kPyBeautyppMNRwmi) {
+ fProcess != kPyBeautyppMNRwmi &&
+ fProcess != kPyBeautyJets) {
for (i = 0; i < np; i++) {
TParticle* iparticle = (TParticle *) fParticles.At(i);
Int_t* pParent = new Int_t[np];
for (i=0; i< np; i++) pParent[i] = -1;
- if (fProcess == kPyJets || fProcess == kPyDirectGamma) {
+ if (fProcess == kPyJets || fProcess == kPyDirectGamma || fProcess == kPyBeautyJets || fProcess == kPyBeautyppMNRwmi) {
TParticle* jet1 = (TParticle *) fParticles.At(6);
TParticle* jet2 = (TParticle *) fParticles.At(7);
if (!CheckTrigger(jet1, jet2)) {
if(!ok)
return 0;
}
+
+ // Select beauty jets with electron in EMCAL
+ if (fProcess == kPyBeautyJets && fEleInEMCAL) {
+
+ Bool_t ok = kFALSE;
+
+ Int_t pdg = 11; //electron
+
+ Float_t pt = 0.;
+ Float_t eta = 0.;
+ Float_t phi = 0.;
+ for (i=0; i< np; i++) {
+ TParticle* iparticle = (TParticle *) fParticles.At(i);
+ if(iparticle->GetStatusCode()==1 && TMath::Abs(iparticle->GetPdgCode())==pdg &&
+ iparticle->Pt() > fElectronMinPt){
+ pt = iparticle->Pt();
+ phi = iparticle->Phi()*180./TMath::Pi(); //Convert to degrees
+ eta =TMath::Abs(iparticle->Eta()); //in calos etamin=-etamax
+ if(IsInEMCAL(phi,eta))
+ ok =kTRUE;
+ }
+ }
+ if(!ok)
+ 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 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;
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
// Check if there is a ccbar or bbbar pair with at least one of the two
// in fYMin < y < fYMax
- if (fProcess == kPyCharmppMNRwmi || fProcess == kPyBeautyppMNRwmi) {
+
+ if (fProcess == kPyCharmppMNRwmi || fProcess == kPyBeautyppMNRwmi || fProcess == kPyBeautyJets) {
TParticle *partCheck;
TParticle *mother;
Bool_t theQ=kFALSE,theQbar=kFALSE,inYcut=kFALSE;
Int_t km = iparticle->GetFirstMother();
if ((ks == 1 && kf!=0 && KinematicSelection(iparticle, 0)) ||
(ks != 1) ||
- (fProcess == kPyJets && ks == 21 && km == 0 && i>1)) {
+ ((fProcess == kPyJets || fProcess == kPyBeautyJets) && ks == 21 && km == 0 && i>1)) {
nc++;
if (ks == 1) trackIt = 1;
Int_t ipa = iparticle->GetFirstMother()-1;
//
// Jets that have triggered
- if (fProcess == kPyJets || fProcess == kPyDirectGamma)
+ //Need to store jets for b-jet studies too!
+ if (fProcess == kPyJets || fProcess == kPyDirectGamma || fProcess == kPyBeautyJets || fProcess == kPyBeautyppMNRwmi)
{
Int_t ntrig, njet;
Float_t jets[4][10];
} else if (fQuench == 4) {
// QPythia
Double_t xy[2];
- AliFastGlauber::Instance()->GetLastXY(xy);
+ Double_t i0i1[2];
+ AliFastGlauber::Instance()->GetSavedXY(xy);
+ AliFastGlauber::Instance()->GetSavedI0I1(i0i1);
xp = xy[0];
yp = xy[1];
- printf("%13.3e %13.3e \n", xp, yp);
+ ((AliGenPythiaEventHeader*) fHeader)->SetImpactParameter(fImpact);
}
((AliGenPythiaEventHeader*) fHeader)->SetXYJet(xp, yp);
pdg[1] = jet2->GetPdgCode();
Bool_t triggered = kFALSE;
- if (fProcess == kPyJets) {
+ if (fProcess == kPyJets || fProcess == kPyBeautyJets || fProcess == kPyBeautyppMNRwmi) {
Int_t njets = 0;
Int_t ntrig = 0;
Float_t jets[4][10];