#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::NAStrangePotential ( )

inlineprivate

Definition at line 34 of file StrangePotential.h.

34 {};

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 ( )

inline

Definition at line 47 of file StrangePotential.h.

47 {};

Member Function Documentation

double NAStrangePotential::CalculateStrangeDensity ( const double strangePotential )

private

Definition at line 46 of file StrangePotential.cc.

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

Referenced by operator()().

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

double NAStrangePotential::CalculateStrangePotential ( )

Definition at line 11 of file StrangePotential.cc.

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

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

                                                      {
   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;
 }

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

inline

Definition at line 49 of file StrangePotential.h.

References CalculateStrangeDensity(), and fStrangeDensity.

Referenced by CalculateStrangePotential().

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