Event/DstEvent/include/DstEvent/DstMdcDedx.h

Go to the documentation of this file.

#ifndef DST_DSTMDCDEDX_H
#define DST_DSTMDCDEDX_H
#include "EventModel/EventModel.h"
#include "GaudiKernel/ContainedObject.h"
#include "GaudiKernel/ObjectVector.h"
 
using namespace EventModel;
extern const CLID& CLID_DstMdcDedx;
enum pid_dedx { electron, muon, kaon, pion, proton };
 
class DstMdcDedx : virtual public ContainedObject {
public:
  // ctor and dector
  DstMdcDedx() {}
  virtual ~DstMdcDedx() {}
 
  virtual const CLID& clID() const { return DstMdcDedx::classID(); }
  static const CLID&  classID() { return CLID_DstMdcDedx; }
 
  // extractors
  // int     getTrackId()    const {return m_trackId ;}
  // int    getParticleId()      const { return m_particleId; }
  // int    getStatus()     const { return m_status; }
  // int    geTruncAlg() const { return m_trunc_alg;  }
  // double getNumChidedx(int i)        const { return m_numChidedx[i]; }
  // int    getNumGoodHits()      const { return m_numGoodHits; }
  // int    getNumTotalHits()     const { return m_numTotalHits; }
  // double getProbPH()           const { return m_probPH; }
  // double getNormPH()        const { return m_normPH; }
  // extractors
  int      trackId() const { return m_trackId; }
  int      particleId() const { return m_particleId; } // will converge to particleType() later
  pid_dedx particleType() const {
    switch ( particleId() )
    {
    case 0: return electron; break;
    case 1: return muon; break;
    case 2: return kaon; break;
    case 3: return pion; break;
    case 4: return proton; break;
    default:
      std::cerr << "the particleType receive an incorrect input!" << std::endl;
      exit( 1 );
    }
  }
  int    status() const { return m_status; }
  int    truncAlg() const { return m_trunc_alg; }
  double chi( int i ) const { return m_chi[i]; }
  double chiE() const { return m_chi[0]; }
  double chiMu() const { return m_chi[1]; }
  double chiPi() const { return m_chi[2]; }
  double chiK() const { return m_chi[3]; }
  double chiP() const { return m_chi[4]; }
  int    numGoodHits() const { return m_numGoodHits; }
  int    numTotalHits() const { return m_numTotalHits; }
  double probPH() const { return m_probPH; }
  double normPH() const { return m_normPH; }
  double errorPH() const { return m_errorPH; }
  double twentyPH() const { return m_twentyPH; }
 
  // modifiers
  void setTrackId( int trackId ) { m_trackId = trackId; }
  void setParticleId( int particleId ) { m_particleId = particleId; }
  void setStatus( int status ) { m_status = status; }
  void setTruncAlg( int trunc_alg ) { m_trunc_alg = trunc_alg; }
 
  void setChi( double* chi ) {
    for ( int i = 0; i < 5; i++ ) m_chi[i] = chi[i];
  }
 
  void setNumGoodHits( int numGoodHits ) { m_numGoodHits = numGoodHits; }
  void setNumTotalHits( int numTotalHits ) { m_numTotalHits = numTotalHits; }
  void setProbPH( double probPH ) { m_probPH = probPH; }
  void setNormPH( double normPH ) { m_normPH = normPH; }
  void setErrorPH( double errorPH ) { m_errorPH = errorPH; }
  void setTwentyPH( double twentyPH ) { m_twentyPH = twentyPH; }
 
protected:
  int    m_trackId{ 0 };      // Track ID Add 2005-10-18
  int    m_particleId{ 0 };   // Particle ID from De/Dx
  int    m_status{ 0 };       // Status
  int    m_trunc_alg{ 0 };    // truncate method
  double m_chi[5]{};          // Number of Chi_dEdx for different particles
                              // 0: e  1: mu   2: pi  3: K  4:p
  int    m_numGoodHits{ 0 };  // No. of good de/dx hits(exclude overflow)
  int    m_numTotalHits{ 0 }; // No. of good de/dx hits(include overflow)
  double m_probPH{ 0.0 };     // Most probable pulse height from truncated mean
  double m_normPH{ 0.0 };     // normalized pulse height
  double m_errorPH{ 0.0 };    // resolution of truncated mean
  double m_twentyPH{ 0.0 };   // de/dx value of bitrunction rate equal to 20%
};
 
typedef ObjectVector<DstMdcDedx> DstMdcDedxCol;
 
#endif // DST_DEDX_H