From 3848c666076072bb8c3c398b471785cd9655bcf5 Mon Sep 17 00:00:00 2001 From: coppedis Date: Thu, 21 Jun 2007 12:25:37 +0000 Subject: [PATCH] Updated version of deuteron emission in nuclear fragment production algorithm --- ZDC/AliZDCDigitizer.cxx | 4 ++++ ZDC/AliZDCFragment.cxx | 10 +++++----- ZDC/AliZDCFragment.h | 2 +- 3 files changed, 10 insertions(+), 6 deletions(-) diff --git a/ZDC/AliZDCDigitizer.cxx b/ZDC/AliZDCDigitizer.cxx index 18147db6d1d..8fa35e6524c 100644 --- 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 19f128e7bbd..174f6ae705d 100644 --- 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 79eafd056fa..425379d7f5e 100644 --- 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;}; -- 2.39.3