Introducing the interaction time into the aliroot generators. In case of gaussian...

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

diff --git a/PYTHIA6/AliGenPythia.cxx b/PYTHIA6/AliGenPythia.cxx
index e0ab292b16993dc3498b174374ad8b5a93ff4af5..6c46455f9ef92c8f159aba4d072f441ca7b999a3 100644 (file)
--- a/PYTHIA6/AliGenPythia.cxx
+++ b/PYTHIA6/AliGenPythia.cxx
@@ -834,7 +834,7 @@ void AliGenPythia::Generate()
                    origin[1] = fVertex[1]+iparticle->Vy()/10; // [cm]
                    origin[2] = fVertex[2]+iparticle->Vz()/10; // [cm]
                    
-                   Float_t tof   = kconv*iparticle->T();
+                   Float_t tof   = fTime + kconv*iparticle->T();
                    Int_t ipa     = iparticle->GetFirstMother()-1;
                    Int_t iparent = (ipa > -1) ? pParent[ipa] : -1;
  
@@ -1174,7 +1174,7 @@ Int_t  AliGenPythia::GenerateMB()
            origin[1] = fVertex[1]+iparticle->Vy()/10; // [cm]
            origin[2] = fVertex[2]+iparticle->Vz()/10; // [cm]
            
-           Float_t tof = fEventTime + kconv * iparticle->T();
+           Float_t tof = fTime + fEventTime + kconv * iparticle->T();
 
            PushTrack(fTrackIt*trackIt, iparent, kf, 
                      p[0], p[1], p[2], p[3], 
@@ -1263,7 +1263,7 @@ void AliGenPythia::MakeHeader()
 //
 // Event Vertex 
     fHeader->SetPrimaryVertex(fVertex);
-    fHeader->SetInteractionTime(fEventTime);
+    fHeader->SetInteractionTime(fTime+fEventTime);
 //
 // Number of primaries
     fHeader->SetNProduced(fNprimaries);
@@ -1898,8 +1898,8 @@ Double_t w[]={
       return kTRUE;
     }
 
-    Int_t ibin=nbin;
-    for(Int_t i=1; i<nbin; i++) 
+    Int_t ibin=nbin-1;
+    for(Int_t i=1; i<=nbin; i++) 
       if(M<=bin[i]) {
        ibin=i-1;
        //      printf("Mi> %f && Mi< %f\n", bin[i-1], bin[i]);