BmnGlobalTrack.h

Go to the documentation of this file.

/**  BmnGlobalTrack.h
 *@author S.Merts
 **
 ** Data class for Global BMN track. Data level RECO.
 **
 **/
 
#ifndef BmnGlobalTrack_H_
#define BmnGlobalTrack_H_ 1
 
#include "BmnTrack.h"
 
class BmnGlobalTrack : public BmnTrack
{
  public:
    /** Default constructor **/
    BmnGlobalTrack();
 
    /** Destructor **/
    virtual ~BmnGlobalTrack();
 
    enum ParticleType
    {
        kCommon = 0,
        kPion = 1,
        kProton = 2,
        kDeutron = 3,
        kTriton = 4,
        kNumParticleTypes = 5
    };
    TString GetParticleName(ParticleType type);
 
    /** Accessors **/
    Int_t GetGemTrackIndex() const { return fInnerTrack; }
    Int_t GetInnerTrackIndex() const { return fInnerTrack; }
    Int_t GetTof1HitIndex() const { return fTof400Hit[0]; }
    Int_t GetTof2HitIndex() const { return fTof700Hit[0]; }
    Int_t GetTof1HitIndex(ParticleType type) const { return fTof400Hit[type]; }
    Int_t GetTof2HitIndex(ParticleType type) const { return fTof700Hit[type]; }
    Int_t GetNdetHitIndex() const { return fNDetHitIndex; }
    Int_t GetScWallCellId() const { return fScWallCellId; }
 
    Float_t GetMass2(Int_t tofID) const
    {
        return (tofID == 1) ? fMass2Tof400[0] : (tofID == 2) ? fMass2Tof700[0] : -1000.0;
    };
    Float_t GetMass2(Int_t tofID, ParticleType type)
    {
        return (tofID == 1) ? fMass2Tof400[type] : (tofID == 2) ? fMass2Tof700[type] : -1000.0;
    };
    Float_t GetResX(Int_t tofID) const
    {
        return (tofID == 1) ? fTof400ResX[0] : (tofID == 2) ? fTof700ResX[0] : -1000.0;
    };
    Float_t GetResX(Int_t tofID, ParticleType type)
    {
        return (tofID == 1) ? fTof400ResX[type] : (tofID == 2) ? fTof700ResX[type] : -1000.0;
    };
    Float_t GetResY(Int_t tofID) const
    {
        return (tofID == 1) ? fTof400ResY[0] : (tofID == 2) ? fTof700ResY[0] : -1000.0;
    };
    Float_t GetResY(Int_t tofID, ParticleType type)
    {
        return (tofID == 1) ? fTof400ResY[type] : (tofID == 2) ? fTof700ResY[type] : -1000.0;
    };
    Float_t GetScWallSignal() const { return fScWallSignal; }
    Float_t GetDedx() const { return fDedx; }
    Int_t GetQ() { return fQ; }
    Float_t GetPx();
    Float_t GetPy();
    Float_t GetPz();
    Float_t GetPt();
    Float_t GetEta();
    Float_t GetPhi();
    Float_t GetTheta();
 
    Int_t GetRefIndex() { return -1; }   // save  for compatibility
    Float_t GetChi2InVertex() { return fChi2InVertex; }
    Float_t GetDCAInVertex() { return fDCAInVertex; }
 
    /** Modifiers **/
    void SetGemTrackIndex(Int_t iGem) { fInnerTrack = iGem; }
    void SetTof1HitIndex(Int_t iTof1Hit) { fTof400Hit[0] = iTof1Hit; }
    void SetTof2HitIndex(Int_t iTof2Hit) { fTof700Hit[0] = iTof2Hit; }
    void SetTof1HitIndex(Int_t iTof1Hit, ParticleType type) { fTof400Hit[type] = iTof1Hit; }
    void SetTof2HitIndex(Int_t iTof2Hit, ParticleType type) { fTof700Hit[type] = iTof2Hit; }
    void SetNdetHitIndex(Int_t iNdet) { fNDetHitIndex = iNdet; }
    void SetScWallCellId(Int_t iScWall) { fScWallCellId = iScWall; }
    void SetScWallSignal(Float_t sig) { fScWallSignal = sig; }
 
    void SetQ(Int_t q) { fQ = q; }
    void SetDedx(Float_t dedx) { fDedx = dedx; }
    void SetMass2(Int_t tofID, Float_t m2)
    {
        if (tofID == 1)   // tof400
            fMass2Tof400[0] = m2;
        else if (tofID == 2)   // tof700
            fMass2Tof700[0] = m2;
    }
    void SetMass2(Int_t tofID, ParticleType type, Float_t m2)
    {
        if (tofID == 1)   // tof400
            fMass2Tof400[type] = m2;
        else if (tofID == 2)   // tof700
            fMass2Tof700[type] = m2;
    }
    void SetResX(Int_t tofID, Float_t res)
    {
        if (tofID == 1)   // tof400
            fTof400ResX[0] = res;
        else if (tofID == 2)   // tof700
            fTof700ResX[0] = res;
    }
    void SetResX(Int_t tofID, ParticleType type, Float_t res)
    {
        if (tofID == 1)   // tof400
            fTof400ResX[type] = res;
        else if (tofID == 2)   // tof700
            fTof700ResX[type] = res;
    }
    void SetResY(Int_t tofID, Float_t res)
    {
        if (tofID == 1)   // tof400
            fTof400ResY[0] = res;
        else if (tofID == 2)   // tof700
            fTof700ResY[0] = res;
    }
    void SetResY(Int_t tofID, ParticleType type, Float_t res)
    {
        if (tofID == 1)   // tof400
            fTof400ResY[type] = res;
        else if (tofID == 2)   // tof700
            fTof700ResY[type] = res;
    }
 
    void SetChi2InVertex(Float_t chi) { fChi2InVertex = chi; }
    void SetDCAInVertex(Float_t chi) { fDCAInVertex = chi; }
 
    void Print();
 
  private:
    /** Indices of local parts of global track **/
    Int_t fInnerTrack;
    Int_t fScWallCellId;
    Int_t fNDetHitIndex;
 
    // vectors of indices of tof hits for different hypothesis:
    vector<Int_t> fTof400Hit;
    vector<Int_t> fTof700Hit;
 
    // vectors of residuals of tof hits for different hypothesis:
    vector<Float_t> fTof400ResX;
    vector<Float_t> fTof700ResX;
    vector<Float_t> fTof400ResY;
    vector<Float_t> fTof700ResY;
 
    // vectors of mass for different hypothesis:
    vector<Float_t> fMass2Tof400;
    vector<Float_t> fMass2Tof700;
 
    Float_t fScWallSignal;
    Float_t fDedx;
    Short_t fQ;
    Float_t fChi2InVertex;
    Float_t fDCAInVertex;
 
    ClassDef(BmnGlobalTrack, 7);
};
 
#endif