Updated version of deuteron emission in nuclear fragment production algorithm

diff --git a/ZDC/AliZDCDigitizer.cxx b/ZDC/AliZDCDigitizer.cxx
index 18147db6d1dbfd5f7f25437d25a8ef1c82a297c3..8fa35e6524c3b6622504797a08ab29716eabf0ca 100644 (file)
--- a/ZDC/AliZDCDigitizer.cxx
+++ b/ZDC/AliZDCDigitizer.cxx
@@ -330,6 +330,10 @@ void AliZDCDigitizer::Fragmentation(Float_t impPar, Int_t specN, Int_t specP,
   frag.AttachNeutrons(zz, nn, ztot, ntot);
   freeSpecN = specN-ntot-2*nAlpha;
   freeSpecP = specP-ztot-2*nAlpha;
+  // Removing deuterons
+  Int_t ndeu = (Int_t) (frag.DeuteronNumber());
+  freeSpecN -= ndeu;
+  //
   if(freeSpecN<0) freeSpecN=0;
   if(freeSpecP<0) freeSpecP=0;
   AliDebug(2, Form("FreeSpn = %d, FreeSpp = %d", freeSpecN, freeSpecP));

diff --git a/ZDC/AliZDCFragment.cxx b/ZDC/AliZDCFragment.cxx
index 19f128e7bbd8b4c292ccbf62b2fbdc253dced928..174f6ae705df6683d632b6a803806d94a3b17696 100644 (file)
--- a/ZDC/AliZDCFragment.cxx
+++ b/ZDC/AliZDCFragment.cxx
@@ -356,12 +356,12 @@ void AliZDCFragment::AttachNeutrons(Int_t *fZZ, Int_t *fNN, Int_t &fZtot,Int_t &
 }
 
 //_____________________________________________________________________________
-Float_t AliZDCFragment::DeuteronFraction()
+Float_t AliZDCFragment::DeuteronNumber()
 {
     // Calculates the fraction of deuterum nucleus produced
     //
-    Float_t deuteronProdPar[2] = {1.068,0.0385};
-    Float_t deutFrac = deuteronProdPar[0]-deuteronProdPar[1]*fB;
-    if(deutFrac>1.) deutFrac=1.;
-    return deutFrac;
+    Float_t deuteronProdPar[2] = {-0.068,0.0385};
+    Float_t deutNum = deuteronProdPar[0] + deuteronProdPar[1]*fB;
+    if(deutNum<0.) deutNum = 0.;
+    return deutNum;
 }

diff --git a/ZDC/AliZDCFragment.h b/ZDC/AliZDCFragment.h
index 79eafd056faa63640b18ab8becaf380420fea04d..425379d7f5e80754615e3e900ef399dc5a3e7b94 100644 (file)
--- a/ZDC/AliZDCFragment.h
+++ b/ZDC/AliZDCFragment.h
@@ -24,7 +24,7 @@ public:
   virtual      ~AliZDCFragment() {}
   void GenerateIMF(Int_t* fZZ, Int_t &fNalpha);
   void AttachNeutrons(Int_t* fZZ, Int_t* fNN, Int_t &Ztot, Int_t &Ntot);
-  Float_t DeuteronFraction();
+  Float_t DeuteronNumber();
   
   // Setting parameters
   virtual void SetImpactParameter(Float_t b) {fB=b;};