TofEPID Class Reference

#include <TofEPID.h>

Inheritance diagram for TofEPID:

Public Member Functions
	~TofEPID ()
void	init ()
void	calculate ()
bool	IsPidInfoValid () const
double	chi (int n) const
double	prob (int n) const
double	sigma (int n) const
double	offset (int n) const
int	ndof () const
double	mass2 () const
int	part () const
double	rhit () const
int	neuronPID () const
	~TofEPID ()
void	init ()
void	calculate ()
bool	IsPidInfoValid () const
double	chi (int n) const
double	prob (int n) const
double	sigma (int n) const
double	offset (int n) const
int	ndof () const
double	mass2 () const
int	part () const
double	rhit () const
int	neuronPID () const
	~TofEPID ()
void	init ()
void	calculate ()
bool	IsPidInfoValid () const
double	chi (int n) const
double	prob (int n) const
double	sigma (int n) const
double	offset (int n) const
int	ndof () const
double	mass2 () const
int	part () const
double	rhit () const
int	neuronPID () const
Public Member Functions inherited from ParticleIDBase
	ParticleIDBase ()
virtual	~ParticleIDBase ()
EvtRecTrack *	PidTrk () const
void	setRecTrack (EvtRecTrack *trk)
double	chiMinCut () const
void	setChiMinCut (const double chi=4)
double	chiMaxCut () const
void	setChiMaxCut (const double chi=6)
double	pdfMinSigmaCut () const
void	setPdfMinSigmaCut (const double pdf=4)
double	getRunNo () const
void	setRunNo (const double runh=8093)
double	p ()
double	pt ()
double	charge ()
double	xmass (int n)
double	velc ()
double	probCalculate (double chi2, int n)
double	pdfCalculate (double offset, double sigma)
double	interpolation (double *x, double *y, double x1)
double	pol2 (double x, double *par)
double	pol3 (double x, double *par)
double	pol4 (double x, double *par)
int	useDedx () const
int	useTof () const
int	useTof1 () const
int	useTof2 () const
int	useTofE () const
int	useTofQ () const
int	useTofC () const
int	useTofCorr () const
int	useEmc () const
int	useMuc () const
int	onlyPionKaon () const
int	onlyPionKaonProton () const
int	onlyPionKaonElectron () const
int	all () const
int	onlyElectron () const
int	onlyMuon () const
int	onlyPion () const
int	onlyKaon () const
int	onlyProton () const
int	methodLikelihood () const
int	methodProbability () const
int	methodNeuronNetwork () const
int	dedxValid () const
int	tofValid () const
int	tofeValid () const
int	tofqValid () const
int	tofcValid () const
int	tofcorrValid () const
int	emcValid () const
int	mucValid () const
void	set_path (const char *s_path=0)
void	set_path (std::string s_path)
	ParticleIDBase ()
virtual	~ParticleIDBase ()
EvtRecTrack *	PidTrk () const
void	setRecTrack (EvtRecTrack *trk)
double	chiMinCut () const
void	setChiMinCut (const double chi=4)
double	chiMaxCut () const
void	setChiMaxCut (const double chi=6)
double	pdfMinSigmaCut () const
void	setPdfMinSigmaCut (const double pdf=4)
double	getRunNo () const
void	setRunNo (const double runh=8093)
double	p ()
double	pt ()
double	charge ()
double	xmass (int n)
double	velc ()
double	probCalculate (double chi2, int n)
double	pdfCalculate (double offset, double sigma)
double	interpolation (double *x, double *y, double x1)
double	pol2 (double x, double *par)
double	pol3 (double x, double *par)
double	pol4 (double x, double *par)
int	useDedx () const
int	useTof () const
int	useTof1 () const
int	useTof2 () const
int	useTofE () const
int	useTofQ () const
int	useTofC () const
int	useTofCorr () const
int	useEmc () const
int	useMuc () const
int	onlyPionKaon () const
int	onlyPionKaonProton () const
int	onlyPionKaonElectron () const
int	all () const
int	onlyElectron () const
int	onlyMuon () const
int	onlyPion () const
int	onlyKaon () const
int	onlyProton () const
int	methodLikelihood () const
int	methodProbability () const
int	methodNeuronNetwork () const
int	dedxValid () const
int	tofValid () const
int	tofeValid () const
int	tofqValid () const
int	tofcValid () const
int	tofcorrValid () const
int	emcValid () const
int	mucValid () const
void	set_path (const char *s_path=0)
void	set_path (std::string s_path)
	ParticleIDBase ()
virtual	~ParticleIDBase ()
EvtRecTrack *	PidTrk () const
void	setRecTrack (EvtRecTrack *trk)
double	chiMinCut () const
void	setChiMinCut (const double chi=4)
double	chiMaxCut () const
void	setChiMaxCut (const double chi=6)
double	pdfMinSigmaCut () const
void	setPdfMinSigmaCut (const double pdf=4)
double	getRunNo () const
void	setRunNo (const double runh=8093)
double	p ()
double	pt ()
double	charge ()
double	xmass (int n)
double	velc ()
double	probCalculate (double chi2, int n)
double	pdfCalculate (double offset, double sigma)
double	interpolation (double *x, double *y, double x1)
double	pol2 (double x, double *par)
double	pol3 (double x, double *par)
double	pol4 (double x, double *par)
int	useDedx () const
int	useTof () const
int	useTof1 () const
int	useTof2 () const
int	useTofE () const
int	useTofQ () const
int	useTofC () const
int	useTofCorr () const
int	useEmc () const
int	useMuc () const
int	onlyPionKaon () const
int	onlyPionKaonProton () const
int	onlyPionKaonElectron () const
int	all () const
int	onlyElectron () const
int	onlyMuon () const
int	onlyPion () const
int	onlyKaon () const
int	onlyProton () const
int	methodLikelihood () const
int	methodProbability () const
int	methodNeuronNetwork () const
int	dedxValid () const
int	tofValid () const
int	tofeValid () const
int	tofqValid () const
int	tofcValid () const
int	tofcorrValid () const
int	emcValid () const
int	mucValid () const
void	set_path (const char *s_path=0)
void	set_path (std::string s_path)

Static Public Member Functions
static TofEPID *	instance ()
static TofEPID *	instance ()
static TofEPID *	instance ()

Protected Member Functions
int	neuronPIDCalculation ()
int	particleIDCalculation ()
int	LikelihoodCalculation ()
double	offsetTofE (int n, int cntr, double ptrk, double rtof, double ph, double charge)
double	sigmaTofE (int n, int cntr, double ptrk, double rtof, double ph, double charge)
int	neuronPIDCalculation ()
int	particleIDCalculation ()
int	LikelihoodCalculation ()
double	offsetTofE (int n, int cntr, double ptrk, double rtof, double ph, double charge)
double	sigmaTofE (int n, int cntr, double ptrk, double rtof, double ph, double charge)
int	neuronPIDCalculation ()
int	particleIDCalculation ()
int	LikelihoodCalculation ()
double	offsetTofE (int n, int cntr, double ptrk, double rtof, double ph, double charge)
double	sigmaTofE (int n, int cntr, double ptrk, double rtof, double ph, double charge)

Additional Inherited Members
Static Protected Attributes inherited from ParticleIDBase
static std::string	path = ""

Detailed Description

Definition at line 13 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

Constructor & Destructor Documentation

~TofEPID() [1/3]

TofEPID::~TofEPID ( )

inline

Definition at line 17 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ ; }

~TofEPID() [2/3]

TofEPID::~TofEPID ( )

inline

Definition at line 17 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ ; }

~TofEPID() [3/3]

TofEPID::~TofEPID ( )

inline

Definition at line 17 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ ; }

Member Function Documentation

calculate() [1/3]

void TofEPID::calculate ( )

virtual

Implements ParticleIDBase.

Definition at line 40 of file TofEPID.cxx.

                        {
  if ( particleIDCalculation() == 0 ) m_ndof = 1;
}

calculate() [2/3]

void TofEPID::calculate ( )

virtual

Implements ParticleIDBase.

calculate() [3/3]

void TofEPID::calculate ( )

virtual

Implements ParticleIDBase.

chi() [1/3]

double TofEPID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 22 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return m_chi[n]; }

chi() [2/3]

double TofEPID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 22 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return m_chi[n]; }

chi() [3/3]

double TofEPID::chi ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 22 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return m_chi[n]; }

init() [1/3]

void TofEPID::init ( )

virtual

Implements ParticleIDBase.

Definition at line 25 of file TofEPID.cxx.

                   {
  for ( int i = 0; i < 5; i++ )
  {
    m_chi[i]    = 99.0;
    m_prob[i]   = -1.0;
    m_offset[i] = 99.0;
    m_sigma[i]  = 1.0;
  }
  m_chimin = 99.;
  m_pdfmin = 99.;
  m_ndof   = 0;
  m_mass2  = -999;
  m_rhit   = -99;
}

init() [2/3]

void TofEPID::init ( )

virtual

Implements ParticleIDBase.

init() [3/3]

void TofEPID::init ( )

virtual

Implements ParticleIDBase.

instance() [1/3]

TofEPID * TofEPID::instance ( )

static

Definition at line 16 of file TofEPID.cxx.

                           {
  if ( !m_pointer ) m_pointer = new TofEPID();
  return m_pointer;
}

instance() [2/3]

TofEPID * TofEPID::instance ( )

static

instance() [3/3]

TofEPID * TofEPID::instance ( )

static

IsPidInfoValid() [1/3]

bool TofEPID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 21 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return ( m_ndof > 0 ); }

IsPidInfoValid() [2/3]

bool TofEPID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 21 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return ( m_ndof > 0 ); }

IsPidInfoValid() [3/3]

bool TofEPID::IsPidInfoValid ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 21 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return ( m_ndof > 0 ); }

LikelihoodCalculation() [1/3]

int TofEPID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 35 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return -1; }

LikelihoodCalculation() [2/3]

int TofEPID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 35 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return -1; }

LikelihoodCalculation() [3/3]

int TofEPID::LikelihoodCalculation ( )

inlineprotectedvirtual

Implements ParticleIDBase.

Definition at line 35 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return -1; }

mass2() [1/3]

double TofEPID::mass2 ( ) const

inline

Definition at line 27 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return m_mass2; }

mass2() [2/3]

double TofEPID::mass2 ( ) const

inline

Definition at line 27 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return m_mass2; }

mass2() [3/3]

double TofEPID::mass2 ( ) const

inline

Definition at line 27 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return m_mass2; }

ndof() [1/3]

int TofEPID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 26 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return m_ndof; }

ndof() [2/3]

int TofEPID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 26 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return m_ndof; }

ndof() [3/3]

int TofEPID::ndof ( ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 26 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return m_ndof; }

neuronPID() [1/3]

int TofEPID::neuronPID ( ) const

inline

Definition at line 30 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return -1; }

neuronPID() [2/3]

int TofEPID::neuronPID ( ) const

inline

Definition at line 30 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return -1; }

neuronPID() [3/3]

int TofEPID::neuronPID ( ) const

inline

Definition at line 30 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return -1; }

neuronPIDCalculation() [1/3]

int TofEPID::neuronPIDCalculation ( )

inlineprotected

Definition at line 33 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return -1; }

neuronPIDCalculation() [2/3]

int TofEPID::neuronPIDCalculation ( )

inlineprotected

Definition at line 33 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return -1; }

neuronPIDCalculation() [3/3]

int TofEPID::neuronPIDCalculation ( )

inlineprotected

Definition at line 33 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return -1; }

offset() [1/3]

double TofEPID::offset ( int n ) const

inline

Definition at line 25 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return m_offset[n]; }

Referenced by offsetTofE().

offset() [2/3]

double TofEPID::offset ( int n ) const

inline

Definition at line 25 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return m_offset[n]; }

offset() [3/3]

double TofEPID::offset ( int n ) const

inline

Definition at line 25 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return m_offset[n]; }

offsetTofE() [1/3]

double TofEPID::offsetTofE ( int n, int cntr, double ptrk, double rtof, double ph, double charge )

protected

Definition at line 215 of file TofEPID.cxx.

                                            {
  double offset=0.0;
  //   double gb = ptrk/xmass(n);
  switch ( n )
  {
  case 0: { // Electron
    double ptemp = ptrk;
    if ( ptrk < 0.2 ) ptemp = 0.2;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    offset      = 0.001 * ( -28.8481 + 138.159 * plog - 249.334 * plog * plog );
    break;
  }
 
  case 1: { // Muon
    double ptemp = ptrk;
    if ( ptrk < 0.2 ) ptemp = 0.2;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    offset      = 0.001 * ( -33.6966 + 1.91915 * plog - 0.592320 * plog * plog );
    break;
  }
  case 2: { // Pion
    double ptemp = ptrk;
    if ( ptrk < 0.2 ) ptemp = 0.2;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    offset      = 0.001 * ( -27.9965 + 1.213 * plog - 2.02740 * plog * plog );
    break;
  }
  case 3: { // Kaon
    double ptemp = ptrk;
    if ( ptrk < 0.3 ) ptemp = 0.3;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    offset      = 0.001 * ( -23.4842 - 28.7569 * plog + 78.21 * plog * plog );
    break;
  }
 
  case 4: { // Proton
    double ptemp = ptrk;
    if ( ptrk < 0.4 ) ptemp = 0.4;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    if ( charge > 0 ) offset = 0.001 * ( -4.854 - 110.540 * plog + 99.8732 * plog * plog );
    if ( charge < 0 ) offset = 0.001 * ( 27.047 - 145.120 * plog + 167.014 * plog * plog );
    break;
  }
 
  default: offset = 0.0; break;
  }
  //  offset = 0.0;
  return offset;
}

offsetTofE() [2/3]

double TofEPID::offsetTofE ( int n, int cntr, double ptrk, double rtof, double ph, double charge )

protected

offsetTofE() [3/3]

double TofEPID::offsetTofE ( int n, int cntr, double ptrk, double rtof, double ph, double charge )

protected

part() [1/3]

int TofEPID::part ( ) const

inline

Definition at line 28 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return m_part; }

part() [2/3]

int TofEPID::part ( ) const

inline

Definition at line 28 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return m_part; }

part() [3/3]

int TofEPID::part ( ) const

inline

Definition at line 28 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return m_part; }

particleIDCalculation() [1/3]

int TofEPID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

Definition at line 44 of file TofEPID.cxx.

                                   {
  int          irc    = -1;
  EvtRecTrack* recTrk = PidTrk();
  if ( !( recTrk->isMdcTrackValid() ) ) return irc;
  if ( !( recTrk->isTofTrackValid() ) ) return irc;
 
#ifndef BEAN
  SmartRefVector<RecTofTrack>           tofTrackCol = recTrk->tofTrack();
  SmartRefVector<RecTofTrack>::iterator it;
#else
  const std::vector<TTofTrack*>&          tofTrackCol = recTrk->tofTrack();
  std::vector<TTofTrack*>::const_iterator it;
#endif
 
  bool isMRPC = false;
  if ( tofTrackCol.size() != 1 && tofTrackCol.size() != 3 ) { return irc; }
  else
  {
    TofHitStatus* hitst = new TofHitStatus;
    if ( tofTrackCol.size() == 1 )
    {
      it = tofTrackCol.begin();
      hitst->setStatus( ( *it )->status() );
      bool isRaw = hitst->is_raw();
      if ( isRaw ) { return irc; }
      isMRPC         = hitst->is_mrpc();
      bool isReadout = hitst->is_readout();
      bool isCounter = hitst->is_counter();
      bool isCluster = hitst->is_cluster();
      if ( !isReadout || !isCounter || !isCluster ) { return irc; }
    }
    else if ( tofTrackCol.size() == 3 )
    {
      unsigned int icounter = -1;
      unsigned int inumber  = 0;
      for ( it = tofTrackCol.begin(); it != tofTrackCol.end(); it++, inumber++ )
      {
        hitst->setStatus( ( *it )->status() );
        bool isRaw = hitst->is_raw();
        if ( isRaw ) continue;
        isMRPC = hitst->is_mrpc();
        if ( !isMRPC ) { return irc; } // MRPC is possible to have 3 record.
        bool isReadout = hitst->is_readout();
        bool isCounter = hitst->is_counter();
        bool isCluster = hitst->is_cluster();
        if ( !isReadout && isCounter && isCluster ) { icounter = inumber; }
      }
      if ( icounter == -1 ) { return irc; }
      it = tofTrackCol.begin() + icounter;
    }
  }
 
  double tof = ( *it )->tof();
  if ( tof <= 0 ) { return irc; }
  double path         = ( *it )->path();
  m_rhit              = ( *it )->zrhit();
  double       beta   = path / velc() / tof;
  double       beta2  = beta * beta;
  RecMdcTrack* mdcTrk = recTrk->mdcTrack();
  double       ptrk   = mdcTrk->p();
  m_mass2             = ptrk * ptrk * ( 1.0 / beta2 - 1.0 );
 
  double chitemp = 99.;
  double pdftemp = 0;
  for ( unsigned int i = 0; i < 5; i++ )
  {
    double texp      = ( *it )->texp( i );
    m_offset[i]      = tof - texp - ( *it )->toffset( i );
    double sigma_tmp = ( *it )->sigma( 0 );
    if ( !isMRPC )
    {
      if ( sigma_tmp > 0 )
      {
        if ( i < 2 ) { m_sigma[i] = sigma_tmp; }
        else { m_sigma[i] = 1.2 * sigma_tmp; }
      }
      else { m_sigma[i] = 0.12; }
    }
    else
    {
      //       m_sigma[i] = 0.065;
      int strip = ( *it )->strip();
      if ( strip < 0 || strip > 11 ) { m_sigma[i] = 0.065; }
      else
      {
        double p0, p1;
        if ( strip == 0 )
        {
          p0 = 0.16;
          p1 = 0.0;
        }
        else if ( strip == 1 )
        {
          p0 = 0.051094;
          p1 = 0.019467;
        }
        else if ( strip == 2 )
        {
          p0 = 0.056019;
          p1 = 0.012964;
        }
        else if ( strip == 3 )
        {
          p0 = 0.043901;
          p1 = 0.015933;
        }
        else if ( strip == 4 )
        {
          p0 = 0.049750;
          p1 = 0.010473;
        }
        else if ( strip == 5 )
        {
          p0 = 0.048345;
          p1 = 0.008545;
        }
        else if ( strip == 6 )
        {
          p0 = 0.046518;
          p1 = 0.009038;
        }
        else if ( strip == 7 )
        {
          p0 = 0.048803;
          p1 = 0.007251;
        }
        else if ( strip == 8 )
        {
          p0 = 0.047523;
          p1 = 0.008434;
        }
        else if ( strip == 9 )
        {
          p0 = 0.042187;
          p1 = 0.010307;
        }
        else if ( strip == 10 )
        {
          p0 = 0.050337;
          p1 = 0.005951;
        }
        else if ( strip == 11 )
        {
          p0 = 0.054740;
          p1 = 0.005629;
        }
        if ( ptrk < 0.05 ) { m_sigma[i] = p0 + p1 / 0.05; }
        else { m_sigma[i] = p0 + p1 / ptrk; }
      }
    }
    m_chi[i] = m_offset[i] / m_sigma[i];
 
    if ( fabs( m_chi[i] ) < chitemp ) { chitemp = fabs( m_chi[i] ); }
    double ppp = pdfCalculate( m_chi[i], 1 );
    if ( fabs( ppp ) > pdftemp ) { pdftemp = fabs( ppp ); }
  }
  m_chimin = chitemp;
  //   if( pdftemp < pdfCalculate(pdfMinSigmaCut(),1) ) { return irc; }
  //   if( fabs(m_chimin) > chiMinCut() ) { return irc; }
 
  for ( unsigned int i = 0; i < 5; i++ )
  { m_prob[i] = probCalculate( m_chi[i] * m_chi[i], 1 ); }
 
  irc = 0;
  return irc;
}

Referenced by calculate().

particleIDCalculation() [2/3]

int TofEPID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

particleIDCalculation() [3/3]

int TofEPID::particleIDCalculation ( )

protectedvirtual

Implements ParticleIDBase.

prob() [1/3]

double TofEPID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 23 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return m_prob[n]; }

prob() [2/3]

double TofEPID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 23 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return m_prob[n]; }

prob() [3/3]

double TofEPID::prob ( int n ) const

inlinevirtual

Implements ParticleIDBase.

Definition at line 23 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return m_prob[n]; }

rhit() [1/3]

double TofEPID::rhit ( ) const

inline

Definition at line 29 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return m_rhit; }

rhit() [2/3]

double TofEPID::rhit ( ) const

inline

Definition at line 29 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return m_rhit; }

rhit() [3/3]

double TofEPID::rhit ( ) const

inline

Definition at line 29 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return m_rhit; }

sigma() [1/3]

double TofEPID::sigma ( int n ) const

inline

Definition at line 24 of file Analysis/ParticleID/include/ParticleID/TofEPID.h.

{ return m_sigma[n]; }

Referenced by sigmaTofE().

sigma() [2/3]

double TofEPID::sigma ( int n ) const

inline

Definition at line 24 of file InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h.

{ return m_sigma[n]; }

sigma() [3/3]

double TofEPID::sigma ( int n ) const

inline

Definition at line 24 of file InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h.

{ return m_sigma[n]; }

sigmaTofE() [1/3]

double TofEPID::sigmaTofE ( int n, int cntr, double ptrk, double rtof, double ph, double charge )

protected

Definition at line 271 of file TofEPID.cxx.

                                           {
 
  double sigma=0.100;
  //   double gb = ptrk/xmass(n);
  switch ( n )
  {
 
  case 0: { // Electron
    double ptemp = ptrk;
    if ( ptrk < 0.2 ) ptemp = 0.2;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    sigma       = 0.001 * ( 109.974 + 15.2457 * plog + 36.8139 * plog * plog );
 
    break;
  }
 
  case 1: { // Muon
    double ptemp = ptrk;
    if ( ptrk < 0.2 ) ptemp = 0.2;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    sigma       = 0.001 * ( 96.5077 - 2.96232 * plog + 3.12910 * plog * plog );
    break;
  }
 
  case 2: { // pion
    double ptemp = ptrk;
    if ( ptrk < 0.2 ) ptemp = 0.2;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    sigma       = 0.001 * ( 105.447 - 2.08044 * plog + 3.44846 * plog * plog );
    break;
  }
 
  case 3: { // Kaon
    double ptemp = ptrk;
    if ( ptrk < 0.3 ) ptemp = 0.3;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    sigma       = 0.001 * ( 88.8806 - 26.8464 * plog + 113.672 * plog * plog );
    break;
  }
  case 4: { // Proton
    double ptemp = ptrk;
    if ( ptrk < 0.5 ) ptemp = 0.5;
    if ( ptrk > 2.1 ) ptemp = 2.1;
    double plog = log( ptemp );
    if ( charge > 0 ) sigma = 0.001 * ( 96.3534 - 44.1139 * plog + 53.9805 * plog * plog );
    if ( charge < 0 ) sigma = 0.001 * ( 157.345 - 98.7357 * plog + 55.1145 * plog * plog );
    break;
  }
 
  default: sigma = 0.100; break;
  }
  // sigma =1;
  return sigma;
}

sigmaTofE() [2/3]

double TofEPID::sigmaTofE ( int n, int cntr, double ptrk, double rtof, double ph, double charge )

protected

sigmaTofE() [3/3]

double TofEPID::sigmaTofE ( int n, int cntr, double ptrk, double rtof, double ph, double charge )

protected

---

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

• Analysis/ParticleID/include/ParticleID/TofEPID.h
• InstallArea/x86_64-el9-gcc13-dbg/include/ParticleID/TofEPID.h
• InstallArea/x86_64-el9-gcc13-opt/include/ParticleID/TofEPID.h
• TofEPID.cxx