doxy: emit comment when complete

[u/mrichter/AliRoot.git] / EVGEN / AliGenDeuteron.cxx

diff --git a/EVGEN/AliGenDeuteron.cxx b/EVGEN/AliGenDeuteron.cxx
index 03b9794daa3c148144274e1c30343216e3cc347e..327f5cf1f877fe48cb8af3f76a65b1aaf5142a54 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
@@ -196,7 +195,7 @@ Double_t AliGenDeuteron::GetCoalescenceProbability(const TParticle* nucleon1, co
        TVector3 v2(nucleon2->Vx(), nucleon2->Vy(), nucleon2->Vz());
        TVector3 p2(nucleon2->Px(), nucleon2->Py(), nucleon2->Pz());
        
-       Double_t deltaP = 2.*this->GetPcm(p1, kProtonMass, p2, kNeutronMass); // relative momentum in CM frame
+       Double_t deltaP = this->GetPcm(p1, kProtonMass, p2, kNeutronMass); // relative momentum in CM frame
        if( deltaP >= fPmax) return -1.;
        
        Double_t deltaR = (v2-v1).Mag();       // relative distance (cm)
@@ -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);