TCgemKalTrack.h

Go to the documentation of this file.

#ifndef RootEventData_TCgemKalTrack_H
#define RootEventData_TCgemKalTrack_H 1
 
#include "TObject.h"
#include "TString.h"
 
class TCgemKalTrack : public TObject {
    
public:
    
    TCgemKalTrack();
    
    ~TCgemKalTrack ();
    
    //extractors
    Int_t    getTrackId() const { return m_trackId;  }
    Int_t    getStat(const Int_t pid) const { return m_stat[pid]; }        
    Int_t    getStat2(const Int_t pid) const { return m_stat2[pid]; }        
    Double_t getChisq(const Int_t pid) const { return m_chisq[pid];}
    Int_t    getNdf(const Int_t pid) const { return m_ndf[pid]; }
    Int_t    getFirstLayer(const Int_t pid) const { return m_firstLayer[pid]; }
    Int_t    getLastLayer(const Int_t pid) const {return m_lastLayer[pid];}
    
    Double_t getZHelix(Int_t i) const {return m_zhelix[i];}
    Double_t getZError(Int_t i, Int_t j) const {return m_zerror[i][j];}
    //Double_t getPoca(Int_t i) const {return m_poca[i];} 
    Double_t getZHelixE(Int_t i) const {return m_zhelix_e[i];}            
    Double_t getZErrorE(Int_t i, Int_t j) const {return m_zerror_e[i][j];}
    //Double_t getPocaE(Int_t i) const {return m_poca_e[i];}
    Double_t getZHelixMu(Int_t i) const {return m_zhelix_mu[i];}            
    Double_t getZErrorMu(Int_t i, Int_t j) const {return m_zerror_mu[i][j];}
    //Double_t getPocaMu(Int_t i) const {return m_poca_mu[i];}
    Double_t getZHelixK(Int_t i) const {return m_zhelix_k[i];}            
    Double_t getZErrorK(Int_t i, Int_t j) const {return m_zerror_k[i][j];}
    //Double_t getPocaK(Int_t i) const {return m_poca_k[i];}
    Double_t getZHelixP(Int_t i) const {return m_zhelix_p[i];}            
    Double_t getZErrorP(Int_t i, Int_t j) const {return m_zerror_p[i][j];}     
    //Double_t getPocaP(Int_t i) const {return m_poca_p[i];}    
         
    Double_t getFHelix(Int_t i) const {return m_fhelix[i];}
    Double_t getFError(Int_t i, Int_t j) const {return m_ferror[i][j];}
    Double_t getFHelixE(Int_t i) const {return m_fhelix_e[i];}            
    Double_t getFErrorE(Int_t i, Int_t j) const {return m_ferror_e[i][j];}
    Double_t getFHelixMu(Int_t i) const {return m_fhelix_mu[i];}            
    Double_t getFErrorMu(Int_t i, Int_t j) const {return m_ferror_mu[i][j];}
    Double_t getFHelixK(Int_t i) const {return m_fhelix_k[i];}            
    Double_t getFErrorK(Int_t i, Int_t j) const {return m_ferror_k[i][j];}
    Double_t getFHelixP(Int_t i) const {return m_fhelix_p[i];}            
    Double_t getFErrorP(Int_t i, Int_t j) const {return m_ferror_p[i][j];}     
    //modifiers
    void setTrackId (const Int_t trackId) { m_trackId = trackId; }
    void setStat(const Int_t stat, const Int_t pid) {m_stat[pid] = stat;}
    void setStat2(const Int_t stat, const Int_t pid) {m_stat2[pid] = stat;}
    void setChisq(const Double_t chisq, const Int_t pid) {m_chisq[pid] = chisq;}
    void setNdf(const Int_t ndf, const Int_t pid) {m_ndf[pid] = ndf;}
    void setFirstLayer(const Int_t fL, const Int_t pid) { m_firstLayer[pid] = fL; }
    void setLastLayer(const Int_t lL, const Int_t pid){ m_lastLayer[pid] = lL;}
    
    void setZHelix(const Double_t zhelix[5]){
      for (int i=0; i<5; i++)
    m_zhelix[i] = zhelix[i];
    }
    void setZError(const Double_t zerror[5][5]){
      for (int i=0 ; i<5 ; i++)
    for (int j=0; j<=i; j++){
      m_zerror[i][j] = zerror[i][j];
      m_zerror[j][i] = zerror[i][j];
    }
    }
    
    //void setPoca(const Double_t poca[3]){
    //  for(int i=0; i<3; i++) m_poca[i] = poca[i];
    //}
    
    void setZHelixE(const Double_t zhelix_e[5]){   
      for (int i = 0 ; i<5 ; i++)             
    m_zhelix_e[i] = zhelix_e[i];             
    }                                          
    void setZErrorE(const Double_t zerror_e[5][5]){
      for (int i= 0 ; i<5 ; i++)               
    for (int j=0; j<=i; j++){              
      m_zerror_e[i][j] = zerror_e[i][j];       
      m_zerror_e[j][i] = zerror_e[i][j];       
    }
    }   
    //void setPocaE(const Double_t poca_e[3]){
    //  for(int i=0; i<3; i++) m_poca_e[i] = poca_e[i];
    //}
    
    void setZHelixMu(const Double_t zhelix_mu[5]){   
      for (int i = 0 ; i<5 ; i++)             
    m_zhelix_mu[i] = zhelix_mu[i];             
    }                                          
    void setZErrorMu(const Double_t zerror_mu[5][5]){
      for (int i= 0 ; i<5 ; i++)               
    for (int j=0; j<=i; j++){               
      m_zerror_mu[i][j] = zerror_mu[i][j];       
      m_zerror_mu[j][i] = zerror_mu[i][j];                                          }             
    }
    //void setPocaMu(const Double_t poca_mu[3]){
    //  for(int i=0; i<3; i++) m_poca_mu[i] = poca_mu[i];
    //}
    
    void setZHelixK(const Double_t zhelix_k[5]){             
      for (int i = 0 ; i<5 ; i++)                  
    m_zhelix_k[i] = zhelix_k[i];            
    }                                               
    void setZErrorK(const Double_t zerror_k[5][5]){
      for (int i= 0 ; i<5 ; i++)                    
    for (int j=0; j<=i; j++){                    
      m_zerror_k[i][j] = zerror_k[i][j];      
      m_zerror_k[j][i] = zerror_k[i][j];
    }   
    }               
    //void setPocaK(const Double_t poca_k[3]){
    //  for(int i=0; i<3; i++) m_poca_k[i] = poca_k[i];
    //}
    
    void setZHelixP(const Double_t zhelix_p[5]){   
      for (int i = 0 ; i<5 ; i++)                  
    m_zhelix_p[i] = zhelix_p[i];            
    }                                               
    void setZErrorP(const Double_t zerror_p[5][5]){
      for (int i=0; i<5; i++)                     
    for (int j=0; j<=i; j++){                    
      m_zerror_p[i][j] = zerror_p[i][j];      
      m_zerror_p[j][i] = zerror_p[i][j];
    }   
    }                                               
    //void setPocaP(const Double_t poca_p[3]){
    //  for(int i=0; i<3; i++) m_poca_p[i] = poca_p[i];
    //}
 
    void setFHelix(const Double_t fhelix[5]){
      for (int i=0; i<5; i++)
    m_fhelix[i] = fhelix[i];
    }
    void setFError(const Double_t ferror[5][5]){
      for (int i=0 ; i<5 ; i++)
    for (int j=0; j<=i; j++){
      m_ferror[i][j] = ferror[i][j];
      m_ferror[j][i] = ferror[i][j];
    }
    }
    
    void setFHelixE(const Double_t fhelix_e[5]){   
      for (int i = 0 ; i<5 ; i++)             
    m_fhelix_e[i] = fhelix_e[i];             
    }                                          
    void setFErrorE(const Double_t ferror_e[5][5]){
      for (int i= 0 ; i<5 ; i++)               
    for (int j=0; j<=i; j++){              
      m_ferror_e[i][j] = ferror_e[i][j];       
      m_ferror_e[j][i] = ferror_e[i][j];       
    }
    }   
 
  void setFHelixMu(const Double_t fhelix_mu[5]){
      for (int i = 0 ; i<5 ; i++)
        m_fhelix_mu[i] = fhelix_mu[i];
    }                  
    void setFErrorMu(const Double_t ferror_mu[5][5]){
      for (int i= 0 ; i<5 ; i++)
        for (int j=0; j<=i; j++){
          m_ferror_mu[i][j] = ferror_mu[i][j];
          m_ferror_mu[j][i] = ferror_mu[i][j];
        }
    }   
 
  void setFHelixK(const Double_t fhelix_k[5]){
      for (int i = 0 ; i<5 ; i++)
        m_fhelix_k[i] = fhelix_k[i];
    }                  
    void setFErrorK(const Double_t ferror_k[5][5]){
      for (int i= 0 ; i<5 ; i++)
        for (int j=0; j<=i; j++){
          m_ferror_k[i][j] = ferror_k[i][j];
          m_ferror_k[j][i] = ferror_k[i][j];
        }
    }   
 
  void setFHelixP(const Double_t fhelix_p[5]){
      for (int i = 0 ; i<5 ; i++)
        m_fhelix_p[i] = fhelix_p[i];
    }                  
    void setFErrorP(const Double_t ferror_p[5][5]){
      for (int i= 0 ; i<5 ; i++)
        for (int j=0; j<=i; j++){
          m_ferror_p[i][j] = ferror_p[i][j];
          m_ferror_p[j][i] = ferror_p[i][j];
        }
    }   
 
 
private:
   Int_t    m_trackId;       
   Int_t    m_stat[5];  
   Int_t    m_stat2[5];  
   Double_t m_chisq[5]; 
   Int_t    m_ndf[5]; 
   Int_t    m_firstLayer[5];
   Int_t    m_lastLayer[5];    
 
   //Double_t m_poca[3];
   Double_t m_zhelix[5];      // 5 track parameters at zero point for pi   
   Double_t m_zerror[5][5];   // error matrix at zero point for pion       
 
   //Double_t m_poca_e[3];
   Double_t m_zhelix_e[5];    // 5 track parameters at zero point for el   
   Double_t m_zerror_e[5][5]; // error matrix at zero point for electron   
 
   //Double_t m_poca_mu[3];
   Double_t m_zhelix_mu[5];   // 5 track parameters at zero point for mu   
   Double_t m_zerror_mu[5][5];// error matrix at zero point  for muon   ;  
 
   //Double_t m_poca_k[3];   
   Double_t m_zhelix_k[5];    // 5 track parameters at zero point for ka   
   Double_t m_zerror_k[5][5]; // error matrix at zero point for kaon       
  
   //Double_t m_poca_p[3]; 
   Double_t m_zhelix_p[5];    // 5 track parameters at zero point for pr   
   Double_t m_zerror_p[5][5]; // error matrix at zero point for proton     
 
   Double_t m_fhelix[5];      // 5 track parameters at first Cgemhit  for pi   
   Double_t m_ferror[5][5];   // error matrix at first Cgem hit for pion    
   Double_t m_fhelix_e[5];       
   Double_t m_ferror_e[5][5];    
   Double_t m_fhelix_mu[5];      
   Double_t m_ferror_mu[5][5];
   Double_t m_fhelix_k[5];       
   Double_t m_ferror_k[5][5];        
   Double_t m_fhelix_p[5];      
   Double_t m_ferror_p[5][5];      
   
   ClassDef(TCgemKalTrack,2)
};
 
#endif