Fix the generation status code for d

diff --git a/EVGEN/AliGenDeuteron.cxx b/EVGEN/AliGenDeuteron.cxx
index 03b9794daa3c148144274e1c30343216e3cc347e..124ee13ab002daf710e0a9fd5cc10874b38bce48 100644 (file)
--- a/EVGEN/AliGenDeuteron.cxx
+++ b/EVGEN/AliGenDeuteron.cxx
@@ -67,7 +67,6 @@ AliGenDeuteron::AliGenDeuteron(Int_t sign, Double_t pmax, Double_t rmax, Int_t c
  ,fR(0)
  ,fPsiR(0)
  ,fCurStack(0)
- ,fNtrk(0)
 {
 //
 // constructor
@@ -349,11 +348,19 @@ void AliGenDeuteron::PushDeuteron(TParticle* parent1, TParticle* parent2)
        // E^2 = p^2 + m^2
        Double_t energy = TMath::Sqrt(pN.Mag2() + kDeuteronMass*kDeuteronMass);
        
+       Int_t ntrk = 0;
+       Double_t weight = 1;
+       Int_t is = 1; // final state particle
+       
        // Add a new (anti)deuteron to current event stack
        fCurStack->PushTrack(1, -1, fSign*kDeuteronPdg,
                         pN.X(), pN.Y(), pN.Z(), energy,
                         vN.X(), vN.Y(), vN.Z(), parent1->T(),
-                        0., 0., 0., kPNCapture, fNtrk, 1., 0);
+                        0., 0., 0., kPNCapture, ntrk, weight, is);
+       
+       // change the status code of the parents
+       parent1->SetStatusCode(kCluster);
+       parent2->SetStatusCode(kCluster);
        
        // Set kDoneBit for the parents
        parent1->SetBit(kDoneBit);

diff --git a/EVGEN/AliGenDeuteron.h b/EVGEN/AliGenDeuteron.h
index ff2f21a426a1a7d666ee204d3f4f9b038947d00f..5333aee1667d04ac09ff0db289efbb1bdbe6e9f3 100644 (file)
--- a/EVGEN/AliGenDeuteron.h
+++ b/EVGEN/AliGenDeuteron.h
@@ -43,6 +43,7 @@ class AliGenDeuteron: public AliGenerator
 
        enum { kFirstPartner=0, kLowestMomentum, kLowestDistance, kBoth };
        enum { kNone=0, kThermal, kExpansion };
+       enum { kCluster=77 };
        
  private:
  
@@ -71,7 +72,6 @@ class AliGenDeuteron: public AliGenerator
        Double_t fR;          // Projectile/Target nuclear radius (fm)
        Double_t fPsiR;       // Reaction plane angle
        AliStack* fCurStack;  //! current event stack
-       Int_t fNtrk;          //! number of the stored track
        
        ClassDef(AliGenDeuteron,2)
 };