Util.h

Go to the documentation of this file.

#ifndef L1Trigger_TrackFindingTracklet_interface_Util_h
#define L1Trigger_TrackFindingTracklet_interface_Util_h
 
#include <sstream>
#include <cassert>
#include <cmath>
#include <algorithm>
 
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"
 
namespace trklet {
 
  //Converts string in binary to hex (used in writing out memory content)
  inline std::string hexFormat(const std::string& binary) {
    std::stringstream ss;
 
    unsigned int radix = 1, value = 0;
    for (int i = binary.length() - 1; i >= 0; i--) {
      if (binary.at(i) != '0' && binary.at(i) != '1')
        continue;
      value += (binary.at(i) - '0') * radix;
      if (radix == 8) {
        ss << std::hex << value;
        radix = 1;
        value = 0;
      } else
        radix <<= 1;
    }
    if (radix != 1)
      ss << std::hex << value;
 
    std::string str = ss.str() + "x0";
    std::reverse(str.begin(), str.end());
    return str;
  }
 
  //Should be optimized by layer - now first implementation to make sure it works OK
  inline int bendencode(double bend, bool isPS) {
    int ibend = 2.0 * bend;
 
    assert(std::abs(ibend - 2.0 * bend) < 0.1);
 
    if (isPS) {
      if (ibend == 0 || ibend == 1)
        return 0;
      if (ibend == 2 || ibend == 3)
        return 1;
      if (ibend == 4 || ibend == 5)
        return 2;
      if (ibend >= 6)
        return 3;
      if (ibend == -1 || ibend == -2)
        return 4;
      if (ibend == -3 || ibend == -4)
        return 5;
      if (ibend == -5 || ibend == -6)
        return 6;
      if (ibend <= -7)
        return 7;
 
      throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__
                                        << " Unknown bendencode for PS module for bend = " << bend
                                        << " ibend = " << ibend;
    }
 
    if (ibend == 0 || ibend == 1)
      return 0;
    if (ibend == 2 || ibend == 3)
      return 1;
    if (ibend == 4 || ibend == 5)
      return 2;
    if (ibend == 6 || ibend == 7)
      return 3;
    if (ibend == 8 || ibend == 9)
      return 4;
    if (ibend == 10 || ibend == 11)
      return 5;
    if (ibend == 12 || ibend == 13)
      return 6;
    if (ibend >= 14)
      return 7;
    if (ibend == -1 || ibend == -2)
      return 8;
    if (ibend == -3 || ibend == -4)
      return 9;
    if (ibend == -5 || ibend == -6)
      return 10;
    if (ibend == -7 || ibend == -8)
      return 11;
    if (ibend == -9 || ibend == -10)
      return 12;
    if (ibend == -11 || ibend == -12)
      return 13;
    if (ibend == -13 || ibend == -14)
      return 14;
    if (ibend <= -15)
      return 15;
 
    throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__
                                      << " Unknown bendencode for 2S module for bend = " << bend
                                      << " ibend = " << ibend;
  }
 
  //Should be optimized by layer - now first implementation to make sure it works OK
  inline double benddecode(int ibend, bool isPS) {
    if (isPS) {
      if (ibend == 0)
        return 0.25;
      if (ibend == 1)
        return 1.25;
      if (ibend == 2)
        return 2.25;
      if (ibend == 3)
        return 3.25;
      if (ibend == 4)
        return -0.75;
      if (ibend == 5)
        return -1.75;
      if (ibend == 6)
        return -2.75;
      if (ibend == 7)
        return -3.75;
 
      throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__
                                        << " Unknown benddecode for PS module for ibend = " << ibend;
    }
 
    if (ibend == 0)
      return 0.25;
    if (ibend == 1)
      return 1.25;
    if (ibend == 2)
      return 2.25;
    if (ibend == 3)
      return 3.25;
    if (ibend == 4)
      return 4.25;
    if (ibend == 5)
      return 5.25;
    if (ibend == 6)
      return 6.25;
    if (ibend == 7)
      return 7.25;
    if (ibend == 8)
      return -0.75;
    if (ibend == 9)
      return -1.75;
    if (ibend == 10)
      return -2.75;
    if (ibend == 11)
      return -3.75;
    if (ibend == 12)
      return -4.75;
    if (ibend == 13)
      return -5.75;
    if (ibend == 14)
      return -6.75;
    if (ibend == 15)
      return -7.75;
 
    throw cms::Exception("BadConfig") << __FILE__ << " " << __LINE__
                                      << " Unknown benddecode for 2S module for ibend = " << ibend;
  }
 
  inline double bend(double r, double rinv, double stripPitch) {
    constexpr double dr = 0.18;
    double delta = r * dr * 0.5 * rinv;
    double bend = -delta / stripPitch;
    return bend;
  }
 
  inline double rinv(double phi1, double phi2, double r1, double r2) {
    if (r2 <= r1) {  //can not form tracklet
      return 20.0;
    }
 
    double dphi = phi2 - phi1;
    double dr = r2 - r1;
 
    return 2.0 * sin(dphi) / dr / sqrt(1.0 + 2 * r1 * r2 * (1.0 - cos(dphi)) / (dr * dr));
  }
 
};  // namespace trklet
#endif