redefine p0 as 1/2*DeltaP_max
authoreserradi <eulogio.serradilla@cern.ch>
Tue, 16 Dec 2014 09:38:49 +0000 (10:38 +0100)
committermorsch <andreas.morsch@cern.ch>
Fri, 19 Dec 2014 15:27:53 +0000 (16:27 +0100)
EVGEN/AliGenDeuteron.cxx
EVGEN/AliGenDeuteron.h
EVGEN/AliGenLightNuclei.cxx

index 124ee13..327f5cf 100644 (file)
@@ -195,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)
index 5333aee..646f0a0 100644 (file)
@@ -17,7 +17,7 @@ class AliGenDeuteron: public AliGenerator
 
  public:
 
-       AliGenDeuteron(Int_t sign=1, Double_t pmax=0.2, Double_t rmax=2.1, Int_t cluster=0 );
+       AliGenDeuteron(Int_t sign=1, Double_t pmax=0.1, Double_t rmax=2.1, Int_t cluster=0 );
        virtual ~AliGenDeuteron();
 
        virtual void Init();
index 19ffaf1..04d8d62 100644 (file)
@@ -1,5 +1,5 @@
 /**************************************************************************
- * Copyright(c) 2009-2010, ALICE Experiment at CERN, All rights reserved. *
+ * Copyright(c) 2009-2014, ALICE Experiment at CERN, All rights reserved. *
  *                                                                        *
  * Author: The ALICE Off-line Project.                                    *
  * Contributors are mentioned in the code where appropriate.              *
@@ -23,7 +23,7 @@
 // same momentum.
 //
 // By default it starts with He4 nuclei which are the most stable,
-// then He3 nuclei and tritons and finally deuterons. It can also generate
+// then He3 nuclei, tritons and finally deuterons. It can also generate
 // a single nucleus species by disabling the others.
 //
 // Sample code for PYTHIA:
@@ -36,7 +36,7 @@
 //
 //    gener->UsePerEventRates();
 //    gener->AddGenerator(pythia, "PYTHIA", 1);
-//    gener->SetCoalescenceMomentum(0.200); // default (GeV/c)
+//    gener->SetCoalescenceMomentum(0.100); // default (GeV/c)
 //
 //////////////////////////////////////////////////////////////////////
 
@@ -60,7 +60,7 @@ ClassImp(AliGenLightNuclei)
 
 AliGenLightNuclei::AliGenLightNuclei()
 :AliGenCocktail()
-,fP0(0.200)
+,fP0(0.100)
 ,fGenDeuterons(kTRUE)
 ,fGenTritons(kTRUE)
 ,fGenHe3Nuclei(kTRUE)
@@ -173,8 +173,8 @@ Bool_t AliGenLightNuclei::Coalescence(const TParticle* n1, const TParticle* n2)
 // returns true if the nucleons are inside of an sphere of radius p0
 // (assume the nucleons are in the same place e.g. PYTHIA, PHOJET,...)
 //
-       Double_t deltaP = 2.*this->GetPcm(n1->Px(), n1->Py(), n1->Pz(), n1->GetMass(), 
-                                         n2->Px(), n2->Py(), n2->Pz(), n2->GetMass());
+       Double_t deltaP = this->GetPcm(n1->Px(), n1->Py(), n1->Pz(), n1->GetMass(), 
+                                      n2->Px(), n2->Py(), n2->Pz(), n2->GetMass());
        
        return ( deltaP < fP0);
 }