3 Nikolai Amelin, Ludmila Malinina, Timur Pocheptsov (C) JINR/Dubna
4 amelin@sunhe.jinr.ru, malinina@sunhe.jinr.ru, pocheptsov@sunhe.jinr.ru
9 //This class is used to calculate strange potential from
10 //the known initial strange density = 0 at given temperature and baryon potential.
12 #ifndef NAStrangePotential_h
13 #define NAStrangePotential_h 1
15 #ifndef NAStrangeDensity_h
16 #include "StrangeDensity.h"
18 #ifndef NAEquationSolver_h
19 #include "EquationSolver.h"
22 #include "DatabasePDG.h"
25 class NAStrangePotential {
27 NAStrangePotential(const Double_t initialStrangeDensity=0, DatabasePDG* database=0x0) :
30 fStrangeDensity(initialStrangeDensity),
31 fMinStrangePotential(0.0001*GeV),
32 fMaxStrangePotential(0.9*GeV),
34 fNSolverIteration(100),
40 ~NAStrangePotential() {};
42 Double_t operator()(const Double_t strangePotential) {
43 return (fStrangeDensity - this->CalculateStrangeDensity(strangePotential))/fStrangeDensity;
46 void SetTemperature(Double_t value) {fTemperature = value;}
47 void SetBaryonPotential(Double_t value) {fBaryonPotential = value;}
48 void SetMinStrangePotential(Double_t value) {fMinStrangePotential = value;}
49 void SetMaxStrangePotential(Double_t value) {fMaxStrangePotential = value;}
50 Double_t CalculateStrangePotential();
53 NAStrangePotential(const NAStrangePotential&);
54 NAStrangePotential& operator=(const NAStrangePotential&);
56 Double_t fTemperature;
57 Double_t fBaryonPotential;
58 Double_t fStrangeDensity;
59 Double_t fMinStrangePotential;//initial min value of strange potential
60 Double_t fMaxStrangePotential;//initial max value of strange potential
61 Int_t fNIteration; //to find proper [minStrangePotential, maxStrangePotential] interval
62 Int_t fNSolverIteration; //to find root in [minStrangePotential,maxStrangePotential] interval
63 Double_t fTolerance;//to find root
64 DatabasePDG* fDatabase;
66 //compute hadron system strange density through strange potential
67 Double_t CalculateStrangeDensity(const Double_t strangePotential);
68 //default constructor is not accesible