Fix the generation status code for d
authoreserradi <eulogio.serradilla@cern.ch>
Fri, 7 Nov 2014 14:53:49 +0000 (15:53 +0100)
committermorsch <andreas.morsch@cern.ch>
Thu, 13 Nov 2014 12:38:05 +0000 (13:38 +0100)
EVGEN/AliGenDeuteron.cxx
EVGEN/AliGenDeuteron.h

index 03b9794..124ee13 100644 (file)
@@ -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);
index ff2f21a..5333aee 100644 (file)
@@ -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)
 };