G4DNARevertProbability Class Reference

#include <G4DNARevertProbability.hh>

Public Member Functions
	G4DNARevertProbability (double toleranceY)
	G4DNARevertProbability (double xmin, double xmax, double toleranceY)
void	Reset ()
void	SetBoundaries (double xmin, double xmax)
void	SetTolerance (double toleranceY)
bool	Propose (double proposedXValue, double proposedProba, double nextProba, double &returnedValue)
double	Revert (double probaForSearchedTime, std::function< double(double)> &funct)

Public Attributes
double	fXmax
double	fXmin
double	fXmaxDef
double	fXminDef
double	fToleranceProba
double	fIncreasingCumulativeFunction

Detailed Description

Definition at line 43 of file G4DNARevertProbability.hh.

Constructor & Destructor Documentation

G4DNARevertProbability::G4DNARevertProbability ( double  toleranceY )

inline

Definition at line 53 of file G4DNARevertProbability.hh.

                                            :
      fToleranceProba(toleranceY)
  {
    fXmax = 0.;
    fXmin = 0.;
    fXminDef = 0.;
    fXmaxDef = 0.;
    fIncreasingCumulativeFunction = 0.;
  }

G4DNARevertProbability::G4DNARevertProbability ( double  xmin, double  xmax, double  toleranceY  )

inline

Definition at line 63 of file G4DNARevertProbability.hh.

                                                                      :
      fXmax(xmax), fXmin(xmin), fToleranceProba(toleranceY)
  {
    if(fXmax < fXmin)
    {
      double tmp = fXmin;
      fXmin = fXmax;
      fXmax = tmp;
    }

    fXminDef = fXmin;
    fXmaxDef = fXmax;
    fIncreasingCumulativeFunction = 0;
  }

Member Function Documentation

bool G4DNARevertProbability::Propose ( double  proposedXValue, double  proposedProba, double  nextProba, double &  returnedValue  )

inline

Definition at line 104 of file G4DNARevertProbability.hh.

  {
    bool returnFlag = false;

    if(proposedProba < nextProba - fToleranceProba) // proba trop petite ==> augmente
    {
      if(fIncreasingCumulativeFunction > 0) // croissant
      {
        if(proposedXValue > fXmin) fXmin = proposedXValue;
      }
      else if(fIncreasingCumulativeFunction < 0) // decroissant
      {
        if(proposedXValue < fXmax) fXmax = proposedXValue;
      }

      returnedValue = (fXmax + fXmin) / 2.;
      returnFlag = false;
    }
    else if(proposedProba > nextProba + fToleranceProba) // proba trop grande
    {
      if(fIncreasingCumulativeFunction > 0)
      {
        if(proposedXValue < fXmax)
        {
          fXmax = proposedXValue;
        }
      }
      else if(fIncreasingCumulativeFunction < 0)
      {
        if(proposedXValue > fXmin)
        {
          fXmin = proposedXValue;
        }
      }

      returnedValue = (fXmax + fXmin) / 2.;
      returnFlag = false;
    }
    else
    {
      // Assuming search for next proba is increasing
      if(fIncreasingCumulativeFunction < 0)
      {
        fXmin = fXminDef;
        fXmax = proposedXValue;
      }
      else if(fIncreasingCumulativeFunction > 0)
      {
        fXmin = proposedXValue;
        fXmax = fXmaxDef;
      }
      returnFlag = true;
    }

    return returnFlag;
  }

Here is the caller graph for this function:

void G4DNARevertProbability::Reset ( )

inline

Definition at line 78 of file G4DNARevertProbability.hh.

  {
    fXmin = fXminDef;
    fXmax = fXmaxDef;
  }

Here is the caller graph for this function:

double G4DNARevertProbability::Revert ( double  probaForSearchedTime, std::function< double(double)> &  funct  )

inline

Definition at line 168 of file G4DNARevertProbability.hh.

  {
    Reset();
    bool notFound = true;
    double proposedX;
    double x = fXmax;

    fIncreasingCumulativeFunction = (funct(fXmax) - funct(fXmin))
        / (fXmax - fXmin);

    while(notFound)
    {
      double newProba = funct(x);

      if(Propose(x, newProba, probaForSearchedTime, proposedX))
      {
        notFound = false;
      }
      else
      {
        if(x == proposedX)
        {
          return x;
//          G4cout << "BREAK : x= " << x << G4endl;
//          G4cout << "nextProba= " << probaForSearchedTime << G4endl;
//          G4cout << "newProba= " << newProba << G4endl;
//          abort();
        }
        x = proposedX;
      }
    }
    return x;
  }

Here is the call graph for this function:

void G4DNARevertProbability::SetBoundaries ( double  xmin, double  xmax  )

inline

Definition at line 84 of file G4DNARevertProbability.hh.

  {
    fXmax = xmax;
    fXmin = xmin;
    if(fXmax < fXmin)
    {
      double tmp = fXmin;
      fXmin = fXmax;
      fXmax = tmp;
    }

    fXminDef = fXmin;
    fXmaxDef = fXmax;
  }

void G4DNARevertProbability::SetTolerance ( double  toleranceY )

inline

Definition at line 99 of file G4DNARevertProbability.hh.

  {
    fToleranceProba = toleranceY;
  }

Member Data Documentation

double G4DNARevertProbability::fIncreasingCumulativeFunction

Definition at line 51 of file G4DNARevertProbability.hh.

double G4DNARevertProbability::fToleranceProba

Definition at line 50 of file G4DNARevertProbability.hh.

double G4DNARevertProbability::fXmax

Definition at line 46 of file G4DNARevertProbability.hh.

double G4DNARevertProbability::fXmaxDef

Definition at line 48 of file G4DNARevertProbability.hh.

double G4DNARevertProbability::fXmin

Definition at line 47 of file G4DNARevertProbability.hh.

double G4DNARevertProbability::fXminDef

Definition at line 49 of file G4DNARevertProbability.hh.

---

The documentation for this class was generated from the following file:

• source/geant4.10.03.p03/source/processes/electromagnetic/dna/utils/include/G4DNARevertProbability.hh