From: eserradi Date: Fri, 7 Nov 2014 14:53:49 +0000 (+0100) Subject: Fix the generation status code for d X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=3d87465c79ff9b1f7148ed5b4c8e39ed7d386f08;p=u%2Fmrichter%2FAliRoot.git Fix the generation status code for d --- diff --git a/EVGEN/AliGenDeuteron.cxx b/EVGEN/AliGenDeuteron.cxx index 03b9794daa3..124ee13ab00 100644 --- 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 ff2f21a426a..5333aee1667 100644 --- 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) };