Message: Re: Wrong optical photon parameters interpolation Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Warning Re: Wrong optical photon parameters interpolation 

Forum: Processes Involving Optical Photons
Re: None Wrong optical photon parameters interpolation (Porosev Viacheslav)
Re: Feedback Re: Wrong optical photon parameters interpolation (Gumplinger Peter)
Date: 04 Oct, 2016
From: Porosev Viacheslav <Porosev Viacheslav>

I am sorry, but it looks like BUG...

in "G4PhysicsOrderedFreeVector.cc" there is the function to get an energy value:

G4double G4PhysicsOrderedFreeVector::GetEnergy(G4double aValue)
{
        G4double e;
        if (aValue <= GetMinValue()) {
          e = edgeMin;
        } else if (aValue >= GetMaxValue()) {
          e = edgeMax;
        } else { 
          size_t closestBin = FindValueBinLocation(aValue);
          e = LinearInterpolationOfEnergy(aValue, closestBin);
	}
        return e;
}

But instead of:

inline
G4double G4PhysicsOrderedFreeVector::GetMaxLowEdgeEnergy()
{
	return binVector.back();
}

inline
G4double G4PhysicsOrderedFreeVector::GetMinLowEdgeEnergy()
{
	return binVector.front();
}

"GetMinValue()" and "GetMaxValue()" are used !!!! that return function values at given energy point!!!

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Sad: Re: Wrong optical photon parameters interpolation   (Porosev Viacheslav - 04 Oct, 2016)
 Add Message Add Message
to: "Re: Wrong optical photon parameters interpolation"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews