StrangePotential.h

Go to the documentation of this file.

/*                                                                            
                                                                            
Nikolai Amelin, Ludmila Malinina, Timur Pocheptsov (C) JINR/Dubna
amelin@sunhe.jinr.ru, malinina@sunhe.jinr.ru, pocheptsov@sunhe.jinr.ru 
November. 2, 2005                                

*/

//This class is used to calculate strange potential from 
//the known initial strange density = 0 at given temperature and baryon potential.

#ifndef NAStrangePotential_h
#define NAStrangePotential_h 1

#include "StrangeDensity.h"
#include "EquationSolver.h"
#include "DatabasePDG.h"

class NAStrangePotential {
 private:
  double fTemperature;
  double fBaryonPotential;
  double fStrangeDensity;
  double fMinStrangePotential;//initial min value of strange potential 
  double fMaxStrangePotential;//initial max value of strange potential
  int fNIteration; //to find proper [minStrangePotential, maxStrangePotential] interval
  int fNSolverIteration; //to find root in [minStrangePotential,maxStrangePotential] interval
  double fTolerance;//to find root 
  DatabasePDG* fDatabase;
  NAStrangeDensity fGc;
  //compute hadron  system strange density through strange potential
  double CalculateStrangeDensity(const double strangePotential);
  //default constructor is not accesible
  NAStrangePotential(){};

 public:
  NAStrangePotential(const double initialStrangeDensity, DatabasePDG* database) :
    fStrangeDensity(initialStrangeDensity),
    fMinStrangePotential(0.0001*GeV),
    fMaxStrangePotential(0.9*GeV),
    fNIteration(100),
    fNSolverIteration(100),
    fTolerance(1.e-8),
    fDatabase(database)
      {};

    ~NAStrangePotential() {};
   
    double operator()(const double strangePotential) { 
      return (fStrangeDensity - this->CalculateStrangeDensity(strangePotential))/fStrangeDensity; 
    }   

    void SetTemperature(double value) {fTemperature = value;}
    void SetBaryonPotential(double value) {fBaryonPotential = value;}
    void SetMinStrangePotential(double value) {fMinStrangePotential = value;}
    void SetMaxStrangePotential(double value) {fMaxStrangePotential = value;}
    double CalculateStrangePotential();
};

#endif