TSU::Trigo Struct Reference

#include <Trigo.h>

Collaboration diagram for TSU::Trigo:

Static Public Member Functions
static int	atan2leg (TSU::L1TopoDataTypes< 16, 0 > x, TSU::L1TopoDataTypes< 16, 0 > y)
static int	atan2 (TSU::L1TopoDataTypes< 16, 0 > x, TSU::L1TopoDataTypes< 16, 0 > y)

Static Public Attributes
static const std::vector< std::string >	Cosleg
static const std::vector< std::string >	Sinleg
static const std::vector< std::string >	Cos
static const std::vector< std::string >	Sin
static const std::vector< int >	CoslegInt = convertVector(TSU::Trigo::Cosleg)
static const std::vector< int >	SinlegInt = convertVector(TSU::Trigo::Sinleg)
static const std::vector< int >	CosInt = convertVector(TSU::Trigo::Cos)
static const std::vector< int >	SinInt = convertVector(TSU::Trigo::Sin)

Detailed Description

Definition at line 28 of file Trigo.h.

Member Function Documentation

atan2()

int TSU::Trigo::atan2 ( TSU::L1TopoDataTypes< 16, 0 >  x, TSU::L1TopoDataTypes< 16, 0 >  y  )

static

Definition at line 512 of file Trigo.cxx.

                                                                         {
  short int octant=0;
  if((x.value()&(1<<16))&&(y.value()&(1<<16))){ // Ex and Ey negative
    if(std::abs(y.to_float()) < std::abs(x.to_float())) octant = 4; //100
    else octant = 5; // 101
  } else if(x.value()&(1<<16)){
    if(std::abs(y.to_float()) < std::abs(x.to_float())) octant = 3; //011
    else octant = 2; // 010
  } else if(y.value()&(1<<16)){
    if(std::abs(y.to_float()) < std::abs(x.to_float())) octant = 7; //111
    else octant = 6; // 110
  } else { // Ex and Ey both positive
    if(std::abs(y.to_float()) < std::abs(x.to_float())) octant = 0; //000
    else octant = 1; // 001
  }
 
  auto signed_x_f = x;
  auto signed_y_f = y;
  if(octant==0||octant==3||octant==4||octant==7){
    //signed_x_f = (abs(int(x)));
    signed_x_f = x.abs();
    signed_y_f = y.abs();
  } else {
    //signed_x_f = (abs(int(y)));
    signed_x_f = y.abs();
    signed_y_f = x.abs();
  }
  
  // For phi in the first sector (0.8 < phi < 1.6) store the possible values of the tan
  const std::vector<std::string> atan2_thresholds =
    {
     "00000010010001101", // tan(0.85)=1.1383327132284393 rounded to 1.1376953125
     "00000010100001010", // tan(0.9)=1.2601582175503392 rounded to 1.259765625
     "00000010110010111", // tan(0.95)=1.398382589287699 rounded to 1.3974609375
     "00000011000111010", // tan(1.0)=1.5574077246549023 rounded to 1.556640625
     "00000011011111001", // tan(1.05)=1.7433153099831704 rounded to 1.7431640625
     "00000011111011011", // tan(1.1)=1.9647596572486523 rounded to 1.9638671875
     "00000100011110000", // tan(1.15)=2.2344969487553255 rounded to 2.234375
     "00000101001001001", // tan(1.2)=2.5721516221263188 rounded to 2.5712890625
     "00000110000001001", // tan(1.25)=3.0095696738628313 rounded to 3.0087890625
     "00000111001101000", // tan(1.3)=3.6021024479679786 rounded to 3.6015625
     "00001000111010010", // tan(1.35)=4.455221759562705 rounded to 4.455078125
     "00001011100110001", // tan(1.4)=5.797883715482887 rounded to 5.7978515625
     "00010000011110011", // tan(1.4500000000000002)=8.238092752965619 rounded to 8.2373046875
     "00011100001100111", // tan(1.5)=14.101419947171719 rounded to 14.1005859375
     "01100000001010000" // tan(1.55)=48.07848247921907 rounded to 48.078125
    };
  
  unsigned int thresholdsPassed = atan2_thresholds.size();
  for(size_t j=0;j<atan2_thresholds.size(); ++j){
    auto thresh =  TSU::L1TopoDataTypes<17,10>(atan2_thresholds.at(j));
    auto yy = signed_y_f*thresh;
    if(! (signed_x_f.to_float() > yy.to_float()) ) {
      //if the current threshold is not passed, then the previous one must have been the last one passed 
      //(no) off-by-one consideration: j starts at 0, so if the first threshold is already failed, that's just the number we're looking for
      thresholdsPassed = j;
      break;
    }
  }
  
  //for first half of a quadrant fine steps are counted in downward 
  // (range: 15 downto 0)
  unsigned int octant_fine = 16 - (thresholdsPassed + 1);
 
  //for second half of a quadrant count fine steps upwards
  // (range: 0 to 15)
  unsigned int octant_fine2 = thresholdsPassed;
 
  int intphi = octant << 4;
  if(octant==0||octant==2||octant==4||octant==6){
    intphi |= octant_fine;
  } else {
    intphi |= octant_fine2;
  }
 
  return intphi;
}

atan2leg()

int TSU::Trigo::atan2leg ( TSU::L1TopoDataTypes< 16, 0 >  x, TSU::L1TopoDataTypes< 16, 0 >  y  )

static

Definition at line 440 of file Trigo.cxx.

                                                                            {
  short int octant=0;
  if((x.value()&(1<<16))&&(y.value()&(1<<16))){ // Ex and Ey negative
    if(std::abs(y.to_float()) < std::abs(x.to_float())) octant = 4; //100
    else octant = 5; // 101
  } else if(x.value()&(1<<16)){
    if(std::abs(y.to_float()) < std::abs(x.to_float())) octant = 3; //011
    else octant = 2; // 010
  } else if(y.value()&(1<<16)){
    if(std::abs(y.to_float()) < std::abs(x.to_float())) octant = 7; //111
    else octant = 6; // 110
  } else { // Ex and Ey both positive
    if(std::abs(y.to_float()) < std::abs(x.to_float())) octant = 0; //000
    else octant = 1; // 001
  }
 
  auto signed_x_f = x;
  auto signed_y_f = y;
  if(octant==0||octant==3||octant==4||octant==7){
    //signed_x_f = (abs(int(x)));
    signed_x_f = x.abs();
    signed_y_f = y.abs();
  } else {
    //signed_x_f = (abs(int(y)));
    signed_x_f = y.abs();
    signed_y_f = x.abs();
  }
 
  const std::vector<std::string> atan2_thresholds = {
                             "000010011101100",  //  1.260 (0.9)
                             "000011000010101", //  1.557 (1.0)
                             "000011110101100", //  1.964 (1.1)
                             "000101000001100", //  2.572 (1.2)
                             "000111000010010", //  3.602 (1.3)
                             "001011010100101", //  5.797 (1.4)
                             "011011100010101"}; // 14.101 (1.5)
  unsigned int flag = 0;
  for(size_t j=0;j<atan2_thresholds.size(); ++j){
    auto thresh =  TSU::L1TopoDataTypes<15,10>(atan2_thresholds.at(j));
    auto yy = signed_y_f*thresh;
    if(signed_x_f.to_float() > yy.to_float()) flag = (flag | (1<<j));
  }
  unsigned int octant_fine = 0;
  if(flag==0) octant_fine = 7;
  else if(flag==1) octant_fine = 6;
  else if(flag==3) octant_fine = 5;
  else if(flag==7) octant_fine = 4;
  else if(flag==15) octant_fine = 3;
  else if(flag==31) octant_fine = 2;
  else if(flag==63) octant_fine = 1;
  else if(flag==127) octant_fine = 0;
 
  unsigned int octant_fine2 = 0;
  if(flag==0) octant_fine2 = 0;
  else if(flag==1) octant_fine2 = 1;
  else if(flag==3) octant_fine2 = 2;
  else if(flag==7) octant_fine2 = 3;
  else if(flag==15) octant_fine2 = 4;
  else if(flag==31) octant_fine2 = 5;
  else if(flag==63) octant_fine2 = 6;
  else if(flag==127) octant_fine2 = 7;
 
  int intphi = octant << 3;
  if(octant==0||octant==2||octant==4||octant==6){
    intphi |= octant_fine;
  } else {
    intphi |= octant_fine2;
  }
 
  return intphi;
}

Member Data Documentation

Cos

const std::vector< std::string > TSU::Trigo::Cos

static

Definition at line 38 of file Trigo.h.

CosInt

const std::vector< int > TSU::Trigo::CosInt = convertVector(TSU::Trigo::Cos)

static

Definition at line 42 of file Trigo.h.

Cosleg

const std::vector< std::string > TSU::Trigo::Cosleg

static

Definition at line 36 of file Trigo.h.

CoslegInt

const std::vector< int > TSU::Trigo::CoslegInt = convertVector(TSU::Trigo::Cosleg)

static

Definition at line 40 of file Trigo.h.

Sin

const std::vector< std::string > TSU::Trigo::Sin

static

Definition at line 39 of file Trigo.h.

SinInt

const std::vector< int > TSU::Trigo::SinInt = convertVector(TSU::Trigo::Sin)

static

Definition at line 43 of file Trigo.h.

Sinleg

const std::vector< std::string > TSU::Trigo::Sinleg

static

Definition at line 37 of file Trigo.h.

SinlegInt

const std::vector< int > TSU::Trigo::SinlegInt = convertVector(TSU::Trigo::Sinleg)

static

Definition at line 41 of file Trigo.h.

---

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

• Trigo.h
• Trigo.cxx