L1Triplet.h

Go to the documentation of this file.

#ifndef L1Triplet_H
#define L1Triplet_H
 
class L1Triplet
{
 
 private:
 
  THitI w0; // left hit (16b)  index in vStsHits array
  THitI w1; // middle hit(16b)
  THitI w2; // right hit(16b)
  unsigned char      b0; // chi2(5b) + level(3b)
  unsigned char      bl; // level
  unsigned char      b1; // qp (8b)
  unsigned char      st; // staL (4b), staM-1-staL (2b), staR-2-staL (2b)
  float chi2;
  
 public:
 
  unsigned char Cqp;
 
  std::vector<unsigned int> neighbours;
 
 L1Triplet(): w0(0), w1(0), w2(0), b0(0), bl(0), b1(0), st(0), chi2(0.),
    Cqp(0), neighbours() {
    neighbours.clear();
  };
  
  void Set( unsigned int iHitL, unsigned int iHitM, unsigned int iHitR,
            unsigned int iStaL, unsigned int iStaM, unsigned int iStaR,
      unsigned char Level, unsigned char  Qp, 
      float Chi2  
      ){ 
    chi2 = Chi2;
    
    w0 = iHitL;
    w1 = iHitM;
    w2 = iHitR;
 
    Chi2 = sqrt(fabs(Chi2))*31./3.5 ; // 3.5 == 31
    if( Chi2>31 || !finite(Chi2) ) Chi2 = 31;
    b0 = ( (static_cast<unsigned char>( Chi2 ))<<3 ) + (Level%8);
    b1 = Qp;
    bl = Level;
 
    st = (iStaL<<4) + ((iStaM-iStaL-1)<<2) + (iStaR-iStaL-2);
  }
 
  THitI GetLHit() const { 
    return w0; 
  }
  THitI GetMHit() const { 
    return w1; 
  }
  THitI GetRHit() const { 
    return w2; 
  }
 
  // unsigned int GetFirstNeighbour() const { 
  //   return w1; 
  // }
  // unsigned int GetNNeighbours() const { 
  //   return w2; 
  // }
  
  unsigned char GetLevel() const {
//    return b0%8;
      return bl;
  }
 
  unsigned char GetQp() const {
    return b1;
  }
 
  float GetChi2() const {
    // float x = (b0>>3)*3.5/31.;
    // return x*x;
    return chi2;
  }
 
  float GetQpOrig( float MaxInvMom ){
    return(GetQp()*(MaxInvMom*2)/255. - MaxInvMom);
  }
 
  int GetLSta() const { 
    return st>>4; 
  }
  
  int GetMSta() const { 
    return ((st%16)>>2) + GetLSta() + 1; 
  }
  
  int GetRSta() const {
    return (st%4) + GetLSta() + 2; 
  }
 
};
 
#endif