Added DptDpt QA Task

[u/mrichter/AliRoot.git] / TEvtGen / AliGenEvtGen.cxx

diff --git a/TEvtGen/AliGenEvtGen.cxx b/TEvtGen/AliGenEvtGen.cxx
index 08d2c58c9633c6f357c4f2b7e768171897033be7..c1d70ad1bb43f86baa93f038260b745ae520358d 100644 (file)
--- a/TEvtGen/AliGenEvtGen.cxx
+++ b/TEvtGen/AliGenEvtGen.cxx
@@ -26,6 +26,7 @@
 #include "AliStack.h"
 #include "AliGenEvtGen.h"
 #include "AliRun.h"
+#include "AliLog.h"
 #include <TParticle.h>
 
 ClassImp(AliGenEvtGen)
@@ -129,12 +130,14 @@ void AliGenEvtGen::Generate()
   //check if particle is already decayed by Pythia  
   if(part->GetStatusCode() != 1 || part->GetNDaughters()>0) 
     {
-    Info("AliGenEvtGen","Attention: particle %d is already decayed by Pythia!",pdg); 
+    AliDebug(1,Form("Attention: particle %d is already decayed by Pythia!",pdg));
     continue;
     }
   
   part->SetStatusCode(11); //Set particle as decayed : change the status code
-  
+  part->SetBit(kDoneBit);
+  part->ResetBit(kTransportBit);
+
   mom->SetPxPyPzE(part->Px(),part->Py(),part->Pz(),part->Energy());
   Int_t np;
 
@@ -153,7 +156,7 @@ void AliGenEvtGen::Generate()
          }
   //select trackable particle  
   if (np >1) {
-  TParticle* iparticle =  (TParticle *) particles->At(0);//parent particle
+    TParticle* iparticle =  0;//(TParticle *) particles->At(0);//parent particle
    for (int i = 1; i<np ; i++) { 
      iparticle = (TParticle*) particles->At(i);
      Int_t ks = iparticle->GetStatusCode();
@@ -178,13 +181,13 @@ void AliGenEvtGen::Generate()
          Int_t jpa  = iparticle->GetFirstMother()-1; //jpa = 0 for daughters of beauty particles
          Int_t iparent = (jpa > 0) ? pParent[jpa] : iTrack;      
 
-         och[0] = origin0[0]+iparticle->Vx()/10; //[cm]
-         och[1] = origin0[1]+iparticle->Vy()/10; //[cm]
-         och[2] = origin0[2]+iparticle->Vz()/10; //[cm]
+         och[0] = origin0[0]+iparticle->Vx(); //[cm]
+         och[1] = origin0[1]+iparticle->Vy(); //[cm]
+         och[2] = origin0[2]+iparticle->Vz(); //[cm]
          pc[0]  = iparticle->Px(); //[GeV/c]
          pc[1]  = iparticle->Py(); //[GeV/c]
          pc[2]  = iparticle->Pz(); //[GeV/c]
-         tof = part->T()+kconv*iparticle->T(); 
+         tof = part->T()+iparticle->T(); 
 
          AliDebug(1,Form("FirstMother = %d e indicePart = %d e pdg = %d \n",jpa,i,kf));