Merging the VirtualMC branch to the main development branch (HEAD)
[u/mrichter/AliRoot.git] / EVGEN / AliGenPythia.cxx
index 25da29c..0820dcb 100644 (file)
 
 /*
 $Log$
+Revision 1.52.4.4  2002/10/10 16:40:08  hristov
+Updating VirtualMC to v3-09-02
+
+Revision 1.62  2002/09/24 10:00:01  morsch
+CheckTrigger() corrected.
+
+Revision 1.61  2002/07/30 09:52:38  morsch
+Call SetGammaPhiRange() and SetGammaEtaRange() in the constructor.
+
+Revision 1.60  2002/07/19 14:49:03  morsch
+Typo corrected.
+
+Revision 1.59  2002/07/19 14:35:36  morsch
+Count total number of trials. Print mean Q, x1, x2.
+
 Revision 1.58  2002/07/17 10:04:09  morsch
 SetYHard method added.
 
@@ -203,6 +218,8 @@ AliGenPythia::AliGenPythia()
   SetEventListRange();
   SetJetPhiRange();
   SetJetEtaRange();
+  SetGammaPhiRange();
+  SetGammaEtaRange();
 }
 
 AliGenPythia::AliGenPythia(Int_t npart)
@@ -233,6 +250,8 @@ AliGenPythia::AliGenPythia(Int_t npart)
     SetEventListRange();
     SetJetPhiRange();
     SetJetEtaRange();
+    SetGammaPhiRange();
+    SetGammaEtaRange();
     // Options determining what to keep in the stack (Heavy flavour generation)
     fStackFillOpt = kFlavorSelection; // Keep particle with selected flavor
     fFeedDownOpt = kTRUE;             // allow feed down from higher family
@@ -345,7 +364,7 @@ void AliGenPythia::Init()
     fQ         = 0.;
     fX1        = 0.;
     fX2        = 0.;    
-    fNev       = 0.;
+    fNev       = 0 ;
 //    
     AliGenMC::Init();
 }
@@ -413,7 +432,8 @@ void AliGenPythia::Generate()
        Int_t nc = 0;        // Total n. of selected particles
        Int_t nParents = 0;  // Selected parents
        Int_t nTkbles = 0;   // Trackable particles
-       if (fProcess != kPyMb && fProcess != kPyJets && fProcess != kPyDirectGamma &&
+       if (fProcess != kPyMb && fProcess != kPyJets && 
+           fProcess != kPyDirectGamma &&
            fProcess != kPyMbNonDiffr) {
            
            for (i = 0; i<np; i++) {
@@ -574,9 +594,6 @@ void AliGenPythia::Generate()
            if (jev >= fNpart || fNpart == -1) {
                fKineBias=Float_t(fNpart)/Float_t(fTrials);
                printf("\n Trials: %i %i %i\n",fTrials, fNpart, jev);
-                      fPythia->GetVINT(41), 
-                      fPythia->GetVINT(42),
-                      fPythia->GetVINT(51));
 
                fQ  += fPythia->GetVINT(51);
                fX1 += fPythia->GetVINT(41);
@@ -716,15 +733,17 @@ Bool_t AliGenPythia::CheckTrigger(TParticle* jet1, TParticle* jet2) const
 
     if (fProcess == kPyJets) {
        //Check eta range first...
-       if ((eta[0] < fEtaMaxJet && eta[0] > fEtaMinJet) ||
-           (eta[1] < fEtaMaxJet && eta[1] > fEtaMinJet))
+       if (
+           ((eta[0] < fEtaMaxJet && eta[0] > fEtaMinJet)   &&
+            (phi[0] < fPhiMaxJet && phi[0] > fPhiMinJet)) 
+           
+           ||
+           
+           ((eta[1] < fEtaMaxJet && eta[1] > fEtaMinJet)   &&
+            (phi[1] < fPhiMaxJet && phi[1] > fPhiMinJet))
+           )
        {
-           //Eta is okay, now check phi range
-           if ((phi[0] < fPhiMaxJet && phi[0] > fPhiMinJet) ||
-               (phi[1] < fPhiMaxJet && phi[1] > fPhiMinJet))
-           {
-               triggered = kTRUE;
-           }
+           triggered = kTRUE;
        }
     } else {
        Int_t ij = 0;
@@ -745,8 +764,6 @@ Bool_t AliGenPythia::CheckTrigger(TParticle* jet1, TParticle* jet2) const
            }
        }
     }
-    
-    
     return triggered;
 }