#include <StrangePotential.h>

Public Member Functions
double	CalculateStrangePotential ()

	NAStrangePotential (const double initialStrangeDensity, DatabasePDG *database)

double	operator() (const double strangePotential)

void	SetBaryonPotential (double value)

void	SetMaxStrangePotential (double value)

void	SetMinStrangePotential (double value)

void	SetTemperature (double value)

	~NAStrangePotential ()

Private Member Functions
double	CalculateStrangeDensity (const double strangePotential)

	NAStrangePotential ()

Private Attributes
double	fBaryonPotential

DatabasePDG *	fDatabase

NAStrangeDensity	fGc

double	fMaxStrangePotential

double	fMinStrangePotential

int	fNIteration

int	fNSolverIteration

double	fStrangeDensity

double	fTemperature

double	fTolerance

Detailed Description

Definition at line 19 of file StrangePotential.h.

Constructor & Destructor Documentation

◆ NAStrangePotential() [1/2]

NAStrangePotential::NAStrangePotential ( )

inlineprivate

Definition at line 34 of file StrangePotential.h.

34 {};

◆ NAStrangePotential() [2/2]

NAStrangePotential::NAStrangePotential	(	const double	initialStrangeDensity,
		DatabasePDG *	database
	)

inline

Definition at line 37 of file StrangePotential.h.

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

◆ ~NAStrangePotential()

NAStrangePotential::~NAStrangePotential ( )

inline

Definition at line 46 of file StrangePotential.h.

46 {};

Member Function Documentation

◆ CalculateStrangeDensity()

double NAStrangePotential::CalculateStrangeDensity ( const double strangePotential )

private

Definition at line 46 of file StrangePotential.cc.

                                                                                 {
   fGc.SetStrangePotential(strangePotential);
   fGc.SetTemperature(fTemperature);
   fGc.SetBaryonPotential(fBaryonPotential);
   return fGc.StrangenessDensity(fDatabase);
 }

References fBaryonPotential, fDatabase, fGc, fTemperature, NAStrangeDensity::SetBaryonPotential(), NAStrangeDensity::SetStrangePotential(), NAStrangeDensity::SetTemperature(), and NAStrangeDensity::StrangenessDensity().

Referenced by operator()().

◆ CalculateStrangePotential()

double NAStrangePotential::CalculateStrangePotential ( )

Definition at line 11 of file StrangePotential.cc.

                                                      {
   double minFunction = this->operator()(fMinStrangePotential);
   double maxFunction = this->operator()(fMaxStrangePotential);
  
   int iter = 0;
   while (minFunction < 0.0 && iter++ < fNIteration) {
     fMinStrangePotential -= 0.5 * fMinStrangePotential;
     minFunction = this->operator()(fMinStrangePotential);
   }
  
   iter = 0;
   while (minFunction * maxFunction > 0.0 && iter++ < fNIteration) {
     fMaxStrangePotential += 1.5 * Abs(fMaxStrangePotential - fMinStrangePotential);
     maxFunction = this->operator()(fMaxStrangePotential);
   }
  
   if (minFunction * maxFunction > 0.0) {
     edm::LogError("StrangePotential") << "CalculateStrangePotential: minFunction*maxFunction is positive!";
     return 0.;
   }
  
   NAEquationSolver<NAStrangePotential>* theSolver =
       new NAEquationSolver<NAStrangePotential>(fNSolverIteration, fTolerance);
  
   theSolver->SetIntervalLimits(fMinStrangePotential, fMaxStrangePotential);
  
   if (!theSolver->Brent(*this))
     edm::LogError("StrangePotential") << "CalculateStrangePotential: the root is not found!";
  
   double strangePotential = theSolver->GetRoot();
   delete theSolver;
   return strangePotential;
 }

References Abs(), NAEquationSolver< Function >::Brent(), fMaxStrangePotential, fMinStrangePotential, fNIteration, fNSolverIteration, fTolerance, NAEquationSolver< Function >::GetRoot(), operator()(), and NAEquationSolver< Function >::SetIntervalLimits().

Referenced by gen::Hydjet2Hadronizer::initializeForInternalPartons().

◆ operator()()

double NAStrangePotential::operator() ( const double strangePotential )

inline

Definition at line 48 of file StrangePotential.h.

                                                    {
     return (fStrangeDensity - this->CalculateStrangeDensity(strangePotential)) / fStrangeDensity;
   }

References CalculateStrangeDensity(), and fStrangeDensity.

Referenced by CalculateStrangePotential().