For quenching option 4 obtain jet production point from AliFastGlauber.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 5 May 2009 12:19:27 +0000 (12:19 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 5 May 2009 12:19:27 +0000 (12:19 +0000)
PYTHIA6/AliGenPythia.cxx

index 4313444..db281c4 100644 (file)
@@ -1236,8 +1236,11 @@ void AliGenPythia::MakeHeader()
            
        } else if (fQuench == 4) {
            // QPythia
-           xp = fPythia->GetPARI(199);
-           yp = fPythia->GetPARI(200);
+           Double_t xy[2];
+           AliFastGlauber::Instance()->GetLastXY(xy);
+           xp = xy[0];
+           yp = xy[1];
+           printf("%13.3e %13.3e \n", xp, yp);
        }
        
            ((AliGenPythiaEventHeader*) fHeader)->SetXYJet(xp, yp);