Store jet production point in header.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 4 May 2009 09:03:49 +0000 (09:03 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 4 May 2009 09:03:49 +0000 (09:03 +0000)
PYTHIA6/AliGenPythia.cxx
PYTHIA6/QPYTHIA/q-pyshow.1.0.F

index 7eaa45c..4313444 100644 (file)
@@ -1224,7 +1224,7 @@ void AliGenPythia::MakeHeader()
        if (fQuench == 1) {
            // Pythia::Quench()
            fPythia->GetQuenchingParameters(xp, yp, z);
-       } else {
+       } else if (fQuench == 2){
            // Pyquen
            Double_t r1 = PARIMP.rb1;
            Double_t r2 = PARIMP.rb2;
@@ -1234,10 +1234,15 @@ void AliGenPythia::MakeHeader()
            xp = r * TMath::Cos(phi);
            yp = r * TMath::Sin(phi);
            
+       } else if (fQuench == 4) {
+           // QPythia
+           xp = fPythia->GetPARI(199);
+           yp = fPythia->GetPARI(200);
        }
+       
            ((AliGenPythiaEventHeader*) fHeader)->SetXYJet(xp, yp);
            ((AliGenPythiaEventHeader*) fHeader)->SetZQuench(z);
-       }
+    }
 //
 // Store pt^hard 
     ((AliGenPythiaEventHeader*) fHeader)->SetPtHard(fPythia->GetVINT(47));
index b4e766f..6ca3249 100644 (file)
@@ -2434,11 +2434,14 @@ C     CENTER-OF-MASS OF THE COLLISION TO THE CENTER-OF-MASS OF THE HARD
 C     SCATTERING. THEY ARE THE ENTRIES THREE TO SEVEN IN ROUTINE PYROBO.
       COMMON/QPLT/AA1,AA2,BBX,BBY,BBZ
 cforalice+
+      COMMON/PYPARS/MSTP(200),PARP(200),MSTI(200),PARI(200)
 c     Here the transverse coordinates of the hard scattering are set by
 c     glauber geometry. 
       call GetRandomXY(xrang,yrang) 
       xin=xrang ! fm
       yin=yrang ! fm
+      PARI(199) = xin
+      PARI(200) = yin
 cforalice-
       zin=0.d0 ! fm
       tin=0.d0 ! fm