1) decouple the rapidity cut used in AliGenPythia to select the events with c or...

[u/mrichter/AliRoot.git] / PYTHIA6 / AliGenPythia.cxx

diff --git a/PYTHIA6/AliGenPythia.cxx b/PYTHIA6/AliGenPythia.cxx
index f2a4023d2c33c55de511ebc6ffd65d7550268307..5ad46d20b2de6d422cfe75427f25e3d0c286d389 100644 (file)
--- a/PYTHIA6/AliGenPythia.cxx
+++ b/PYTHIA6/AliGenPythia.cxx
@@ -109,6 +109,9 @@ AliGenPythia::AliGenPythia():
     fEtaMaxGamma(20.),      
     fPhiMinGamma(0.),      
     fPhiMaxGamma(2. * TMath::Pi()),      
+    fUseYCutHQ(kFALSE),
+    fYMinHQ(-20.),
+    fYMaxHQ(20.),
     fPycellEtaMax(2.),     
     fPycellNEta(274),       
     fPycellNPhi(432),       
@@ -232,6 +235,9 @@ AliGenPythia::AliGenPythia(Int_t npart)
      fEtaMaxGamma(20.),      
      fPhiMinGamma(0.),      
      fPhiMaxGamma(2. * TMath::Pi()),      
+     fUseYCutHQ(kFALSE),
+     fYMinHQ(-20.),
+     fYMaxHQ(20.),
      fPycellEtaMax(2.),     
      fPycellNEta(274),       
      fPycellNPhi(432),       
@@ -1128,7 +1134,8 @@ Int_t  AliGenPythia::GenerateMB()
          if(pdg>0) { theQ=kTRUE; } else { theQbar=kTRUE; }
          y = 0.5*TMath::Log((partCheck->Energy()+partCheck->Pz()+1.e-13)/
                             (partCheck->Energy()-partCheck->Pz()+1.e-13));
-         if(y>fYMin && y<fYMax) inYcut=kTRUE;
+         if(fUseYCutHQ && y>fYMinHQ && y<fYMaxHQ) inYcut=kTRUE;
+         if(!fUseYCutHQ && y>fYMin && y<fYMax) inYcut=kTRUE;
        }
        if(fTriggerParticle) {
          if(TMath::Abs(pdg)==fTriggerParticle) triggered=kTRUE;
@@ -1285,6 +1292,7 @@ void AliGenPythia::MakeHeader()
 // Builds the event header, to be called after each event
     if (fHeader) delete fHeader;
     fHeader = new AliGenPythiaEventHeader("Pythia");
+    fHeader->SetTitle(GetTitle());
 //
 // Event type  
     if(fProcDiff>0){