From 2591bd0ee7f4606df93225e6e7086769dcd74513 Mon Sep 17 00:00:00 2001 From: hristov Date: Tue, 20 Mar 2012 15:00:54 +0000 Subject: [PATCH] Changes for #92223: Modify AliGenPythia.cxx, .h to enable triggering with specific species within a given pt window --- PYTHIA6/AliGenPythia.cxx | 7 ++++++- PYTHIA6/AliGenPythia.h | 7 +++++-- 2 files changed, 11 insertions(+), 3 deletions(-) diff --git a/PYTHIA6/AliGenPythia.cxx b/PYTHIA6/AliGenPythia.cxx index 10843cef809..a6283061154 100644 --- a/PYTHIA6/AliGenPythia.cxx +++ b/PYTHIA6/AliGenPythia.cxx @@ -122,6 +122,8 @@ AliGenPythia::AliGenPythia(): fNucPdf(0), fTriggerParticle(0), fTriggerEta(0.9), + fTriggerMinPt(-1), + fTriggerMaxPt(1000), fTriggerMultiplicity(0), fTriggerMultiplicityEta(0), fTriggerMultiplicityPtMin(0), @@ -224,7 +226,9 @@ AliGenPythia::AliGenPythia(Int_t npart) fHFoff(kFALSE), fNucPdf(0), fTriggerParticle(0), - fTriggerEta(0.9), + fTriggerEta(0.9), + fTriggerMinPt(-1), + fTriggerMaxPt(1000), fTriggerMultiplicity(0), fTriggerMultiplicityEta(0), fTriggerMultiplicityPtMin(0), @@ -1082,6 +1086,7 @@ Int_t AliGenPythia::GenerateMB() if (kf != fTriggerParticle) continue; if (iparticle->Pt() == 0.) continue; if (TMath::Abs(iparticle->Eta()) > fTriggerEta) continue; + if ( iparticle->Pt() > fTriggerMaxPt || iparticle->Pt() < fTriggerMinPt ) continue; triggered = kTRUE; break; } diff --git a/PYTHIA6/AliGenPythia.h b/PYTHIA6/AliGenPythia.h index 4a13eed822f..8ccb53f0abd 100644 --- a/PYTHIA6/AliGenPythia.h +++ b/PYTHIA6/AliGenPythia.h @@ -114,8 +114,9 @@ class AliGenPythia : public AliGenMC // Trigger and rotate event void RotatePhi(Int_t iphcand, Bool_t& okdd); // Trigger on a single particle - virtual void SetTriggerParticle(Int_t particle = 0, Float_t etamax = 0.9) - {fTriggerParticle = particle; fTriggerEta = etamax;} + virtual void SetTriggerParticle(Int_t particle = 0, Float_t etamax = 0.9, Float_t ptmin = -1, Float_t ptmax = 1000) + {fTriggerParticle = particle; fTriggerEta = etamax; fTriggerMinPt = ptmin; fTriggerMaxPt = ptmax;} + // // Heavy flavor options // @@ -276,6 +277,8 @@ class AliGenPythia : public AliGenMC Int_t fNucPdf; // Nuclear pdf 0: EKS98 1: EPS08 Int_t fTriggerParticle; // Trigger on this particle ... Float_t fTriggerEta; // .. within |eta| < fTriggerEta + Float_t fTriggerMinPt; // .. within pt > fTriggerMinPt + Float_t fTriggerMaxPt; // .. within pt < fTriggerMaxPt Int_t fTriggerMultiplicity; // Trigger on events with a minimum charged multiplicity Float_t fTriggerMultiplicityEta; // in a given eta range Float_t fTriggerMultiplicityPtMin; // above this pT -- 2.43.0