KVChIo.cpp

#include "KVChIo.h"
 
ClassImp(KVChIo)
 
 
 
Double_t KVChIo::GetELossMylar(UInt_t z, UInt_t a, Double_t egas, Bool_t stopped)
{
   //Based on energy loss in gas, calculates sum of energy losses in mylar windows
   //from energy loss tables.
   //If argument 'egas' not given, KVChIo::GetEnergy() is used.
   //if stopped=kTRUE, we give the correction for a particle which stops in the detector
   //(by default we assume the particle continues after the detector)
   //
   // WARNING: if stopped=kFALSE, and if the residual energy after the detector
   //   is known (i.e. measured in a detector placed after this one), you should
   //   first call
   //       SetEResAfterDetector(Eres);
   //   before calling this method. Otherwise, especially for heavy ions, the
   //   correction may be false for particles which are just above the punch-through energy.
 
   egas = ((egas < 0.) ? GetEnergy() : egas);
   if (egas <= 0.)
      return 0.0;               //check measured (calibrated) energy in gas is reasonable (>0)
 
   KVNucleus tmp(z, a);
   Double_t emylar = GetCorrectedEnergy(&tmp, egas, !stopped) - egas;
   return emylar;
}