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

News 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: 06 Oct, 2016
From: Porosev Viacheslav <Porosev Viacheslav>

Dear Peter,

I am sorry once again...

but as I understood in "G$" file in "void G4Scintillation::BuildThePhysicsTable()" method, two internal variables are calculated: "fFastIntegralTable" "fSlowIntegralTable". They are commulative integrals of the scintillator emission spectrum. To calculate an integral a standart trapezoidal method is used.

currentCII = 0.5 * (prevIN + currentIN); currentCII = prevCII + (currentPM - prevPM) * currentCII;

and in "PostStepDoIt" to calculate enegy of the photons the following method is used:

// Max Scintillation Integral G4double CIImax = ScintillationIntegral->GetMaxValue();

for (G4int i = 0; i < Num; i++)
      // Determine photon energy
      G4double CIIvalue = G4UniformRand()*CIImax;
      G4double sampledEnergy =  ScintillationIntegral->GetEnergy(CIIvalue);

but "GetEnergy" makes linear interpolation only!!!!

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;

Please, explain how I can change linear interpolation to SPLINE for scintillator emission spectrum!!!!

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

1 Feedback: Re: Wrong optical photon parameters interpolation   (Gumplinger Peter - 06 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