ERSupport.h

  

#ifndef ERSupport_H
#define ERSupport_H

#include <map>
#include <limits>

#include "TMatrixD.h"
#include "TString.h"

#include "FairLogger.h"

using ERChannel = ushort;
namespace consts {
  const unsigned short undefined_channel = std::numeric_limits<ERChannel>::max();
  // Averaged numbers of objects in events. Used as TClonesArray capacitance.
  // Сan affect the efficiency of calculations in case of an obvious 
  // mismatch with the real numbers of objects.
  const int approx_beamdet_point_number = 100;
  const int approx_beamdet_tof_digi_number = 10;
  const int approx_beamdet_mwpc_digi_number = 100;
  const int approx_beamdet_track_number = 10;
  const int approx_beamdet_particle_number = 10;
  const int approx_telescope_point_number = 100;
  const int approx_telescope_digi_number = 10;
  const int approx_telescope_track_number = 10;
  const int approx_telescope_particle_number = 10;
}
enum ERDataObjectType {Point, Digi, Track, Particle};
inline TString ERDataObjectTypeStr(ERDataObjectType object) {
    switch(object) {
        case Point : return "Point";
        case Digi  : return "Digi";
        case Track : return "Track";
        case Particle : return "Particle";
        default : LOG(FATAL) << "Unknown ERDataObject type\n";
    }
    return "";
}
enum ChannelSide {First, Second, None};
enum SensetiveType {Si, CsI, LaBr};
enum OrientationAroundZ {X, Y, Default};
inline TString SensetiveTypeStr(SensetiveType type) {
  switch(type) {
    case Si : return "Si";
    case CsI : return "CsI";
    case LaBr : return "LaBr";
    default : LOG(FATAL) << "Unknown telescope sensetive type\n"; 
  }
  return "";
}
inline TString ChannelSideStr(ChannelSide side) {
  switch(side) {
    case First : return "1";
    case Second : return "2";
    case None : LOG(FATAL) 
        << "String representation for None channel side is not available\n"; 
    default : LOG(FATAL) << "Unknown telescope channel side\n"; 
  }
  return "";
}
inline TString OrientationAroundZStr(OrientationAroundZ type) {
  switch(type) {
    case X : return "X";
    case Y : return "Y";
    case Default : LOG(FATAL) 
        << "String representation for Default orientation around Z is not available\n"; 
    default : LOG(FATAL) << "Unknown telescope orientation around z\n"; 
  }
  return "";
}

class G4ParticleDefinition;
class G4Material;

#define TempNorm 293.16 /* Normal temperature (K) */
#define AMU 931.494028 /* atomic mass unit (MeV) */

TMatrixD* ReadCalFile(TString fileName);
double EiEo(double tableER[][105],double Tp,double Rp);
int ReadRint(char* Fname,double Ranges[][105]);
int intrp4(double* x,double* y, double* c);

Double_t CalcElossIntegralVolStep(Double_t T, const G4ParticleDefinition& ion, 
                                   const G4Material& mat, const Double_t range, 
                                   const Double_t intStep=1e-4);

Double_t ComputeElossIntegralVolStep(Double_t kineticEnergy, const G4ParticleDefinition& particle, 
                                     const G4Material& material, const TString& process_name,
                                     const Double_t range);
class ElossCalculator {
    public:
        static double CalcDeDx_long(const TString& material, int pdg, double kineticEnergy,const Double_t range);
        static double CalcDeDx(const TString& material, int pdg, double kineticEnergy);
    ClassDef(ElossCalculator,1)
};

#endif