emtf Namespace Reference

Namespaces
	phase2

Classes
class	AbsoluteDeviation
struct	CPPFTag
struct	CSCTag
struct	DTTag
struct	Event
class	Forest
struct	GEMTag
class	Huber
struct	IRPCTag
class	LeastSquares
class	LossFunction
struct	ME0Tag
class	Node
class	PercentErrorSquared
struct	RPCTag
class	Tree

Typedefs
using	Feature = std::array< float, NUM_FEATURES >
using	Prediction = std::array< float, NUM_PREDICTIONS >
template<typename T >
using	sector_array = std::array< T, NUM_SECTORS >
template<typename T >
using	zone_array = std::array< T, NUM_ZONES >

Functions
int	calc_chamber (int station, int sector, int subsector, int ring, int csc_ID)
double	calc_eta (int bits)
double	calc_eta_from_theta_deg (double theta_deg, int endcap)
double	calc_eta_from_theta_rad (double theta_rad)
int	calc_eta_GMT (double val)
double	calc_phi_glob_deg (double loc, int sector)
double	calc_phi_glob_rad (double loc, int sector)
double	calc_phi_GMT_deg (int bits)
int	calc_phi_GMT_int (double val)
double	calc_phi_GMT_rad (int bits)
double	calc_phi_loc_deg (int bits)
double	calc_phi_loc_deg_from_glob (double glob, int sector)
int	calc_phi_loc_int (double glob, int sector)
int	calc_phi_loc_int_rpc (double glob, int sector)
double	calc_phi_loc_rad (int bits)
double	calc_pt (int bits)
int	calc_pt_GMT (double val)
int	calc_ring (int station, int csc_ID, int strip)
double	calc_theta_deg (double eta)
double	calc_theta_deg_from_eta (double eta)
double	calc_theta_deg_from_int (int theta_int)
int	calc_theta_int (double theta, int endcap)
int	calc_theta_int_rpc (double theta, int endcap)
double	calc_theta_rad (double eta)
double	calc_theta_rad_from_eta (double eta)
double	calc_theta_rad_from_int (int theta_int)
int	calc_uGMT_chamber (int csc_ID, int subsector, int neighbor, int station)
double	deg_to_rad (double deg)
void	dump_fw_raw_input (const l1t::EMTFHitCollection &out_hits, const l1t::EMTFTrackCollection &out_tracks)
std::pair< int, int >	get_csc_max_pattern_and_quality (int station, int ring)
int	get_csc_max_slope (int station, int ring, bool useRun3CCLUT_OTMB, bool useRun3CCLUT_TMB)
std::pair< int, int >	get_csc_max_strip_and_wire (int station, int ring)
int	get_trigger_csc_ID (int ring, int station, int chamber)
int	get_trigger_sector (int ring, int station, int chamber)
void	mergeNtuples (const char *ntuplename, const char *filestomerge, const char *outputfile)
template<typename T >
std::string	numToStr (T num)
float	processPrediction (float BDTPt, int Quality, float PrelimFit)
double	rad_to_deg (double rad)
template<class bidiiter >
bidiiter	shuffle (bidiiter begin, bidiiter end, size_t num_random)
void	sort_uGMT_muons (l1t::RegionalMuonCandBxCollection &cands)
void	sortNtupleByEvent (const char *ntuplename, const char *filenametosort, const char *outputfile)
double	wrap_phi_deg (double deg)
double	wrap_phi_rad (double rad)
double	wrap_theta_deg (double deg)
double	wrap_theta_rad (double rad)

Variables
constexpr int	MAX_ENDCAP = 2
constexpr int	MAX_TRIGSECTOR = 6
constexpr int	MIN_ENDCAP = 1
constexpr int	MIN_TRIGSECTOR = 1
constexpr int	NUM_FEATURES = 29
constexpr int	NUM_PREDICTIONS = 2
constexpr int	NUM_SECTORS = 12
constexpr int	NUM_STATION_PAIRS = 6
constexpr int	NUM_STATIONS = 4
constexpr int	NUM_ZONE_HITS = 160
constexpr int	NUM_ZONES = 4
const double	ptscale [31]
const std::vector< double >	ptScale
const double	twoJets_scale [16]
const std::vector< double >	twoJetsScale

Typedef Documentation

Feature

using emtf::Feature = typedef std::array<float, NUM_FEATURES>

Definition at line 71 of file Common.h.

Prediction

using emtf::Prediction = typedef std::array<float, NUM_PREDICTIONS>

Definition at line 72 of file Common.h.

sector_array

template<typename T >

using emtf::sector_array = typedef std::array<T, NUM_SECTORS>

Definition at line 63 of file Common.h.

zone_array

template<typename T >

using emtf::zone_array = typedef std::array<T, NUM_ZONES>

Definition at line 65 of file Common.h.

Function Documentation

calc_chamber()

int emtf::calc_chamber	(	int	station,
		int	sector,
		int	subsector,
		int	ring,
		int	csc_ID
	)

Definition at line 28 of file TrackTools.cc.

References relativeConstraints::chamber, relativeConstraints::ring, nano_mu_digi_cff::sector, relativeConstraints::station, and nano_mu_digi_cff::subsector.

Referenced by l1t::stage2::emtf::ImportME().

                                                                                 {
    int chamber = -999;
    if (station == 1) {
      chamber = ((sector - 1) * 6) + csc_ID + 2;  // Chamber offset of 2: First chamber in sector 1 is chamber 3
      if (ring == 2)
        chamber -= 3;
      if (ring == 3)
        chamber -= 6;
      if (subsector == 2)
        chamber += 3;
      if (chamber > 36)
        chamber -= 36;
    } else if (ring == 1) {
      chamber = ((sector - 1) * 3) + csc_ID + 1;  // Chamber offset of 1: First chamber in sector 1 is chamber 2
      if (chamber > 18)
        chamber -= 18;
    } else if (ring == 2) {
      chamber = ((sector - 1) * 6) + csc_ID - 3 + 2;  // Chamber offset of 2: First chamber in sector 1 is chamber 3
      if (chamber > 36)
        chamber -= 36;
    }
    return chamber;
  }

calc_eta()

double emtf::calc_eta ( int  bits )

inline

Definition at line 99 of file TrackTools.h.

References ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::bits, and PVValHelper::eta.

Referenced by l1t::stage2::emtf::ImportSP().

                                   {
    double eta = static_cast<double>(bits);
    eta *= 0.010875;
    return eta;
  }

calc_eta_from_theta_deg()

double emtf::calc_eta_from_theta_deg ( double  theta_deg, int  endcap  )

inline

Definition at line 117 of file TrackTools.h.

References calc_eta_from_theta_rad(), deg_to_rad(), makeMuonMisalignmentScenario::endcap, PVValHelper::eta, and wrap_theta_deg().

Referenced by AngleCalculation::calculate_angles(), PrimitiveConversion::convert_csc_details(), PrimitiveConversion::convert_other_details(), PrimitiveConversion::convert_rpc_details(), and l1t::stage2::emtf::ImportRPC().

                                                                      {  // endcap [-1,+1]
    double theta_rad = deg_to_rad(wrap_theta_deg(theta_deg));            // put theta in [0, 90] range
    double eta = calc_eta_from_theta_rad(theta_rad);
    eta = (endcap == -1) ? -eta : eta;
    return eta;
  }

calc_eta_from_theta_rad()

double emtf::calc_eta_from_theta_rad ( double  theta_rad )

inline

Definition at line 112 of file TrackTools.h.

References PVValHelper::eta, dqm-mbProfile::log, and funct::tan().

Referenced by calc_eta_from_theta_deg().

                                                          {
    double eta = -1. * std::log(std::tan(theta_rad / 2.));
    return eta;
  }

calc_eta_GMT()

int emtf::calc_eta_GMT ( double  val )

inline

Definition at line 124 of file TrackTools.h.

References heppy_batch::val.

                                      {
    val /= 0.010875;
    int gmt_eta = static_cast<int>(std::round(val));
    return gmt_eta;
  }

calc_phi_glob_deg()

double emtf::calc_phi_glob_deg ( double  loc, int  sector  )

inline

Definition at line 170 of file TrackTools.h.

References nano_mu_digi_cff::sector.

Referenced by calc_phi_glob_rad(), AngleCalculation::calculate_angles(), PrimitiveConversion::convert_csc_details(), PrimitiveConversion::convert_other_details(), PrimitiveConversion::convert_rpc_details(), l1t::stage2::emtf::ImportRPC(), and l1t::stage2::emtf::ImportSP().

                                                          {  // loc in deg, sector [1-6]
    double glob = loc + 15. + (60. * (sector - 1));
    glob = (glob < 180.) ? glob : glob - 360.;
    return glob;
  }

calc_phi_glob_rad()

double emtf::calc_phi_glob_rad ( double  loc, int  sector  )

inline

Definition at line 176 of file TrackTools.h.

References calc_phi_glob_deg(), deg_to_rad(), rad_to_deg(), and nano_mu_digi_cff::sector.

                                                          {  // loc in rad, sector [1-6]
    return deg_to_rad(calc_phi_glob_deg(rad_to_deg(loc), sector));
  }

calc_phi_GMT_deg()

double emtf::calc_phi_GMT_deg ( int  bits )

inline

Definition at line 221 of file TrackTools.h.

References ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::bits.

Referenced by calc_phi_GMT_rad().

                                           {
    double phi = static_cast<double>(bits);
    phi = (phi * 360. / 576.) + (180. / 576.);
    return phi;
  }

calc_phi_GMT_int()

int emtf::calc_phi_GMT_int ( double  val )

inline

Definition at line 233 of file TrackTools.h.

References heppy_batch::val, and wrap_phi_deg().

                                          {  // phi in deg
    val = wrap_phi_deg(val);                 // put phi in [-180,180] range
    val = (val - 180. / 576.) / (360. / 576.);
    int gmt_phi = static_cast<int>(std::round(val));
    return gmt_phi;
  }

calc_phi_GMT_rad()

double emtf::calc_phi_GMT_rad ( int  bits )

inline

Definition at line 231 of file TrackTools.h.

References ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::bits, calc_phi_GMT_deg(), and deg_to_rad().

{ return deg_to_rad(calc_phi_GMT_deg(bits)); }

calc_phi_loc_deg()

double emtf::calc_phi_loc_deg ( int  bits )

inline

Definition at line 180 of file TrackTools.h.

References ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::bits.

Referenced by calc_phi_loc_rad(), AngleCalculation::calculate_angles(), PrimitiveConversion::convert_csc_details(), PrimitiveConversion::convert_other_details(), PrimitiveConversion::convert_rpc_details(), l1t::stage2::emtf::ImportRPC(), and l1t::stage2::emtf::ImportSP().

                                           {
    double loc = static_cast<double>(bits);
    loc = (loc / 60.) - 22.;
    return loc;
  }

calc_phi_loc_deg_from_glob()

double emtf::calc_phi_loc_deg_from_glob ( double  glob, int  sector  )

inline

Definition at line 199 of file TrackTools.h.

References nano_mu_digi_cff::sector, and wrap_phi_deg().

Referenced by calc_phi_loc_int(), calc_phi_loc_int_rpc(), and PrimitiveConversion::convert_dt().

                                                                    {  // glob in deg, sector [1-6]
    glob = wrap_phi_deg(glob);                                         // put phi in [-180,180] range
    double loc = glob - 15. - (60. * (sector - 1));
    return loc;
  }

calc_phi_loc_int()

int emtf::calc_phi_loc_int ( double  glob, int  sector  )

inline

Definition at line 205 of file TrackTools.h.

References calc_phi_loc_deg_from_glob(), and nano_mu_digi_cff::sector.

Referenced by PrimitiveConversion::convert_dt(), PrimitiveConversion::convert_gem(), PrimitiveConversion::convert_me0(), and PrimitiveConversion::convert_rpc().

                                                       {  // glob in deg, sector [1-6]
    double loc = calc_phi_loc_deg_from_glob(glob, sector);
    loc = ((loc + 22.) < 0.) ? loc + 360. : loc;
    loc = (loc + 22.) * 60.;
    int phi_int = static_cast<int>(std::round(loc));
    return phi_int;
  }

calc_phi_loc_int_rpc()

int emtf::calc_phi_loc_int_rpc ( double  glob, int  sector  )

inline

Definition at line 213 of file TrackTools.h.

References calc_phi_loc_deg_from_glob(), and nano_mu_digi_cff::sector.

Referenced by PrimitiveConversion::convert_rpc().

                                                           {  // glob in deg, sector [1-6]
    double loc = calc_phi_loc_deg_from_glob(glob, sector);
    loc = ((loc + 22.) < 0.) ? loc + 360. : loc;
    loc = (loc + 22.) * 60. / 4.;  // 4x coarser resolution
    int phi_int = static_cast<int>(std::round(loc));
    return phi_int;
  }

calc_phi_loc_rad()

double emtf::calc_phi_loc_rad ( int  bits )

inline

Definition at line 186 of file TrackTools.h.

References ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::bits, calc_phi_loc_deg(), and deg_to_rad().

{ return deg_to_rad(calc_phi_loc_deg(bits)); }

calc_pt()

double emtf::calc_pt ( int  bits )

inline

Definition at line 84 of file TrackTools.h.

References ALPAKA_ACCELERATOR_NAMESPACE::pixelClustering::pixelStatus::bits, and DiDispStaMuonMonitor_cfi::pt.

                                  {
    double pt = static_cast<double>(bits);
    pt = 0.5 * (pt - 1);
    return pt;
  }

calc_pt_GMT()

int emtf::calc_pt_GMT ( double  val )

inline

Definition at line 90 of file TrackTools.h.

References heppy_batch::val.

                                     {
    val = (val * 2) + 1;
    int gmt_pt = static_cast<int>(std::round(val));
    gmt_pt = (gmt_pt > 511) ? 511 : gmt_pt;
    return gmt_pt;
  }

calc_ring()

int emtf::calc_ring	(	int	station,
		int	csc_ID,
		int	strip
	)

Definition at line 5 of file TrackTools.cc.

References relativeConstraints::station.

Referenced by l1t::stage2::emtf::ImportME(), and l1t::stage2::emtf::MEBlockUnpacker::unpack().

                                                    {
    if (station > 1) {
      if (csc_ID < 4)
        return 1;
      else if (csc_ID < 10)
        return 2;
      else
        return -999;
    } else if (station == 1) {
      if (csc_ID < 4 && strip > 127)
        return 4;
      else if (csc_ID < 4 && strip >= 0)
        return 1;
      else if (csc_ID > 3 && csc_ID < 7)
        return 2;
      else if (csc_ID > 6 && csc_ID < 10)
        return 3;
      else
        return -999;
    } else
      return -999;
  }

calc_theta_deg()

double emtf::calc_theta_deg ( double  eta )

inline

Definition at line 145 of file TrackTools.h.

References calc_theta_rad(), PVValHelper::eta, and rad_to_deg().

{ return rad_to_deg(calc_theta_rad(eta)); }

calc_theta_deg_from_eta()

double emtf::calc_theta_deg_from_eta ( double  eta )

inline

Definition at line 166 of file TrackTools.h.

References calc_theta_rad_from_eta(), PVValHelper::eta, and rad_to_deg().

{ return rad_to_deg(calc_theta_rad_from_eta(eta)); }

calc_theta_deg_from_int()

double emtf::calc_theta_deg_from_int ( int  theta_int )

inline

Definition at line 132 of file TrackTools.h.

References theta().

Referenced by calc_theta_rad_from_int(), AngleCalculation::calculate_angles(), PrimitiveConversion::convert_csc_details(), PrimitiveConversion::convert_other_details(), PrimitiveConversion::convert_rpc_details(), and l1t::stage2::emtf::ImportRPC().

                                                       {
    double theta = static_cast<double>(theta_int);
    theta = theta * (45.0 - 8.5) / 128. + 8.5;
    return theta;
  }

calc_theta_int()

int emtf::calc_theta_int ( double  theta, int  endcap  )

inline

Definition at line 147 of file TrackTools.h.

References makeMuonMisalignmentScenario::endcap, and theta().

Referenced by PrimitiveConversion::convert_dt(), PrimitiveConversion::convert_gem(), PrimitiveConversion::convert_me0(), and PrimitiveConversion::convert_rpc().

                                                      {  // theta in deg, endcap [-1,+1]
    theta = (endcap == -1) ? (180. - theta) : theta;
    theta = (theta - 8.5) * 128. / (45.0 - 8.5);
    int theta_int = static_cast<int>(std::round(theta));
    return theta_int;
  }

calc_theta_int_rpc()

int emtf::calc_theta_int_rpc ( double  theta, int  endcap  )

inline

Definition at line 154 of file TrackTools.h.

References makeMuonMisalignmentScenario::endcap, and theta().

Referenced by PrimitiveConversion::convert_rpc().

                                                          {  // theta in deg, endcap [-1,+1]
    theta = (endcap == -1) ? (180. - theta) : theta;
    theta = (theta - 8.5) * (128. / 4.) / (45.0 - 8.5);  // 4x coarser resolution
    int theta_int = static_cast<int>(std::round(theta));
    return theta_int;
  }

calc_theta_rad()

double emtf::calc_theta_rad ( double  eta )

inline

Definition at line 140 of file TrackTools.h.

References PVValHelper::eta, and JetChargeProducer_cfi::exp.

Referenced by calc_theta_deg().

                                           {
    double theta_rad = 2. * std::atan(std::exp(-1. * eta));
    return theta_rad;
  }

calc_theta_rad_from_eta()

double emtf::calc_theta_rad_from_eta ( double  eta )

inline

Definition at line 161 of file TrackTools.h.

References PVValHelper::eta, and theta().

Referenced by calc_theta_deg_from_eta().

                                                    {
    double theta = std::atan2(1.0, std::sinh(eta));  // cot(theta) = sinh(eta)
    return theta;
  }

calc_theta_rad_from_int()

double emtf::calc_theta_rad_from_int ( int  theta_int )

inline

Definition at line 138 of file TrackTools.h.

References calc_theta_deg_from_int(), and deg_to_rad().

{ return deg_to_rad(calc_theta_deg_from_int(theta_int)); }

calc_uGMT_chamber()

int emtf::calc_uGMT_chamber	(	int	csc_ID,
		int	subsector,
		int	neighbor,
		int	station
	)

Definition at line 53 of file TrackTools.cc.

References relativeConstraints::station, and nano_mu_digi_cff::subsector.

Referenced by l1t::stage2::emtf::SPBlockUnpacker::unpack().

                                                                              {
    if (station == 1) {
      if (csc_ID == 3 && neighbor == 1 && subsector == 2)
        return 1;
      else if (csc_ID == 6 && neighbor == 1 && subsector == 2)
        return 2;
      else if (csc_ID == 9 && neighbor == 1 && subsector == 2)
        return 3;
      else if (csc_ID == 3 && neighbor == 0 && subsector == 2)
        return 4;
      else if (csc_ID == 6 && neighbor == 0 && subsector == 2)
        return 5;
      else if (csc_ID == 9 && neighbor == 0 && subsector == 2)
        return 6;
      else
        return 0;
    } else {
      if (csc_ID == 3 && neighbor == 1)
        return 1;
      else if (csc_ID == 9 && neighbor == 1)
        return 2;
      else if (csc_ID == 3 && neighbor == 0)
        return 3;
      else if (csc_ID == 9 && neighbor == 0)
        return 4;
      else
        return 0;
    }
  }

deg_to_rad()

double emtf::deg_to_rad ( double  deg )

inline

Definition at line 72 of file TrackTools.h.

References ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), and M_PI.

Referenced by calc_eta_from_theta_deg(), calc_phi_glob_rad(), calc_phi_GMT_rad(), calc_phi_loc_rad(), and calc_theta_rad_from_int().

                                       {
    constexpr double factor = M_PI / 180.;
    return deg * factor;
  }

dump_fw_raw_input()

void emtf::dump_fw_raw_input	(	const l1t::EMTFHitCollection &	out_hits,
		const l1t::EMTFTrackCollection &	out_tracks
	)

Definition at line 6 of file DebugTools.cc.

References nano_mu_digi_cff::bx, relativeConstraints::chamber, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), gather_cfg::cout, makeMuonMisalignmentScenario::endcap, h, L1TMuon::kCSC, L1TMuon::kRPC, MAX_ENDCAP, MAX_TRIGSECTOR, MIN_ENDCAP, MIN_TRIGSECTOR, nano_mu_digi_cff::sector, relativeConstraints::station, nano_mu_digi_cff::strip, submitPVValidationJobs::t, validateGeometry_cfg::valid, and nano_mu_digi_cff::wire.

Referenced by TrackFinder::process().

                                                                                                         {
    // from interface/Common.h
    constexpr int MIN_ENDCAP = 1;
    constexpr int MAX_ENDCAP = 2;
    constexpr int MIN_TRIGSECTOR = 1;
    constexpr int MAX_TRIGSECTOR = 6;

    for (int endcap = MIN_ENDCAP; endcap <= MAX_ENDCAP; ++endcap) {
      for (int sector = MIN_TRIGSECTOR; sector <= MAX_TRIGSECTOR; ++sector) {
        const int es = (endcap - MIN_ENDCAP) * (MAX_TRIGSECTOR - MIN_TRIGSECTOR + 1) + (sector - MIN_TRIGSECTOR);

        // _____________________________________________________________________
        // This prints the hits as raw text input to the firmware simulator
        // "12345" is the BX separator

        std::cout << "==== Endcap " << endcap << " Sector " << sector << " Hits ====" << std::endl;
        std::cout << "bx e s ss st vf ql cp wg id bd hs" << std::endl;

        bool empty_sector = true;
        for (const auto& h : out_hits) {
          if (h.Sector_idx() != es)
            continue;
          empty_sector = false;
        }

        for (int ibx = -3 - 5; (ibx < +3 + 5 + 5) && !empty_sector; ++ibx) {
          for (const auto& h : out_hits) {
            if (h.Subsystem() == L1TMuon::kCSC) {
              if (h.Sector_idx() != es)
                continue;
              if (h.BX() != ibx)
                continue;

              int bx = 1;
              int endcap = (h.Endcap() == 1) ? 1 : 2;
              int sector = h.PC_sector();
              int station = (h.PC_station() == 0 && h.Subsector() == 1) ? 1 : h.PC_station();
              int chamber = h.PC_chamber() + 1;
              int strip = (h.Station() == 1 && h.Ring() == 4) ? h.Strip() + 128 : h.Strip();  // ME1/1a
              int wire = h.Wire();
              int valid = 1;
              std::cout << bx << " " << endcap << " " << sector << " " << h.Subsector() << " " << station << " "
                        << valid << " " << h.Quality() << " " << h.Pattern() << " " << wire << " " << chamber << " "
                        << h.Bend() << " " << strip << std::endl;

            } else if (h.Subsystem() == L1TMuon::kRPC) {
              if (h.Sector_idx() != es)
                continue;
              if (h.BX() + 6 != ibx)
                continue;  // RPC hits should be supplied 6 BX later relative to CSC hits

              // Assign RPC link index. Code taken from src/PrimitiveSelection.cc
              int rpc_sub = -1;
              int rpc_chm = -1;
              if (!h.Neighbor()) {
                rpc_sub = ((h.Subsector_RPC() + 3) % 6);
              } else {
                rpc_sub = 6;
              }
              if (h.Station() <= 2) {
                rpc_chm = (h.Station() - 1);
              } else {
                rpc_chm = 2 + (h.Station() - 3) * 2 + (h.Ring() - 2);
              }

              int bx = 1;
              int endcap = (h.Endcap() == 1) ? 1 : 2;
              int sector = h.PC_sector();
              int station = rpc_sub;
              int chamber = rpc_chm + 1;
              int strip = (h.Phi_fp() >> 2);
              int wire = (h.Theta_fp() >> 2);
              int valid = 2;  // this marks RPC stub
              std::cout << bx << " " << endcap << " " << sector << " " << 0 << " " << station << " " << valid << " "
                        << 0 << " " << 0 << " " << wire << " " << chamber << " " << 0 << " " << strip << std::endl;
            }
          }  // end loop over hits

          std::cout << "12345" << std::endl;
        }  // end loop over bx

        // _____________________________________________________________________
        // This prints the tracks as raw text output from the firmware simulator

        std::cout << "==== Endcap " << endcap << " Sector " << sector << " Tracks ====" << std::endl;
        std::cout << "bx e s a mo et ph cr q pt" << std::endl;

        for (const auto& t : out_tracks) {
          if (t.Sector_idx() != es)
            continue;

          std::cout << t.BX() << " " << (t.Endcap() == 1 ? 1 : 2) << " " << t.Sector() << " " << t.PtLUT().address
                    << " " << t.Mode() << " " << (t.GMT_eta() >= 0 ? t.GMT_eta() : t.GMT_eta() + 512) << " "
                    << t.GMT_phi() << " " << t.GMT_charge() << " " << t.GMT_quality() << " " << t.Pt() << std::endl;
        }  // end loop over tracks

      }  // end loop over sector
    }    // end loop over endcap
  }

get_csc_max_pattern_and_quality()

std::pair< int, int > emtf::get_csc_max_pattern_and_quality	(	int	station,
		int	ring
	)

Definition at line 169 of file TrackTools.cc.

Referenced by PrimitiveSelection::select_csc().

                                                                           {
    int max_pattern = 11;
    int max_quality = 16;
    return std::make_pair(max_pattern, max_quality);
  }

get_csc_max_slope()

int emtf::get_csc_max_slope	(	int	station,
		int	ring,
		bool	useRun3CCLUT_OTMB,
		bool	useRun3CCLUT_TMB
	)

Definition at line 175 of file TrackTools.cc.

References or, and relativeConstraints::ring.

Referenced by PrimitiveSelection::select_csc().

                                                                                              {
    int max_slope = 65536;  // Uninitialized slope can be 65536. This is expected when CCLUT is not running
    if (useRun3CCLUT_OTMB and (ring == 1 or ring == 4))
      max_slope = 16;
    if (useRun3CCLUT_TMB and (ring == 2 or ring == 3))
      max_slope = 16;
    return max_slope;
  }

get_csc_max_strip_and_wire()

std::pair< int, int > emtf::get_csc_max_strip_and_wire	(	int	station,
		int	ring
	)

Definition at line 141 of file TrackTools.cc.

References relativeConstraints::ring, and relativeConstraints::station.

Referenced by PrimitiveSelection::select_csc().

                                                                      {
    int max_strip = 0;                // halfstrip
    int max_wire = 0;                 // wiregroup
    if (station == 1 && ring == 4) {  // ME1/1a
      max_strip = 96;
      max_wire = 48;
    } else if (station == 1 && ring == 1) {  // ME1/1b
      max_strip = 128;
      max_wire = 48;
    } else if (station == 1 && ring == 2) {  // ME1/2
      max_strip = 160;
      max_wire = 64;
    } else if (station == 1 && ring == 3) {  // ME1/3
      max_strip = 128;
      max_wire = 32;
    } else if (station == 2 && ring == 1) {  // ME2/1
      max_strip = 160;
      max_wire = 112;
    } else if (station >= 3 && ring == 1) {  // ME3/1, ME4/1
      max_strip = 160;
      max_wire = 96;
    } else if (station >= 2 && ring == 2) {  // ME2/2, ME3/2, ME4/2
      max_strip = 160;
      max_wire = 64;
    }
    return std::make_pair(max_strip, max_wire);
  }

get_trigger_csc_ID()

int emtf::get_trigger_csc_ID	(	int	ring,
		int	station,
		int	chamber
	)

Definition at line 100 of file TrackTools.cc.

References relativeConstraints::chamber, mps_fire::result, relativeConstraints::ring, and relativeConstraints::station.

Referenced by PrimitiveConversion::convert_dt(), PrimitiveConversion::convert_gem(), PrimitiveConversion::convert_me0(), PrimitiveConversion::convert_rpc(), PrimitiveSelection::select_dt(), PrimitiveSelection::select_gem(), and PrimitiveSelection::select_me0().

                                                             {
    int result = 0;
    if (station == 1) {
      result = (chamber) % 3 + 1;  // 1,2,3
      switch (ring) {
        case 1:
          break;
        case 2:
          result += 3;  // 4,5,6
          break;
        case 3:
          result += 6;  // 7,8,9
          break;
        case 4:                            // ME0
          result = (chamber + 1) % 3 + 1;  // 1,2,3
          break;
      }
    } else {
      if (ring == 1) {
        result = (chamber + 1) % 3 + 1;  // 1,2,3
      } else {
        result = (chamber + 3) % 6 + 4;  // 4,5,6,7,8,9
      }
    }
    return result;
  }

get_trigger_sector()

int emtf::get_trigger_sector	(	int	ring,
		int	station,
		int	chamber
	)

Definition at line 85 of file TrackTools.cc.

References relativeConstraints::chamber, mps_fire::result, relativeConstraints::ring, and relativeConstraints::station.

Referenced by PrimitiveConversion::convert_dt(), PrimitiveConversion::convert_gem(), PrimitiveConversion::convert_me0(), PrimitiveSelection::select_dt(), PrimitiveSelection::select_gem(), and PrimitiveSelection::select_me0().

                                                             {
    int result = 0;
    if (station > 1 && ring > 1) {
      result = ((static_cast<unsigned>(chamber - 3) & 0x7f) / 6) + 1;  // ch 3-8->1, 9-14->2, ... 1,2 -> 6
    } else if (station == 1) {
      result = ((static_cast<unsigned>(chamber - 3) & 0x7f) / 6) + 1;  // ch 3-8->1, 9-14->2, ... 1,2 -> 6
    } else {
      result = ((static_cast<unsigned>(chamber - 2) & 0x1f) / 3) + 1;  // ch 2-4-> 1, 5-7->2, ...
    }
    return (result <= 6) ? result
                         : 6;  // max sector is 6, some calculations give a value greater than six but this is expected.
  }

mergeNtuples()

void emtf::mergeNtuples	(	const char *	ntuplename,
		const char *	filestomerge,
		const char *	outputfile
	)

Definition at line 88 of file Utilities.cc.

References writedatasetfile::outputfile.

                                                                                            {
  TChain chain(ntuplename);
  chain.Add(filestomerge);
  chain.Merge(outputfile);
}

numToStr()

template<typename T >

std::string emtf::numToStr

(

T

num

)

Definition at line 43 of file Utilities.h.

References EgammaValidation_cff::num, alignCSCRings::s, contentValuesCheck::ss, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by emtf::Tree::addXMLAttributes().

                            {
    // Convert a number to a string.
    std::stringstream ss;
    ss << num;
    std::string s = ss.str();
    return s;
  };

processPrediction()

float emtf::processPrediction	(	float	BDTPt,
		int	Quality,
		float	PrelimFit
	)

Definition at line 46 of file Utilities.cc.

References RPCpg::pts, and ptscale.

                                                                   {
  // Discretize and scale the BDTPt prediction

  // Fix terrible predictions
  if (BDTPt < 0)
    BDTPt = PrelimFit;
  if (BDTPt > 250)
    BDTPt = PrelimFit;

  float BDTPt1 = BDTPt;
  float scaleF = 1.0;

  // Scale based upon quality
  if (Quality == 3)
    scaleF = 1.15;
  if (Quality == 2)
    scaleF = 1.3;
  if (Quality == 1)
    scaleF = 1.7;

  BDTPt1 = scaleF * BDTPt1;

  // Discretize based upon ptscale
  for (int pts = 0; pts < 31; pts++) {
    if (ptscale[pts] <= BDTPt1 && ptscale[pts + 1] > BDTPt1) {
      BDTPt1 = ptscale[pts];
      break;
    }
  }

  if (BDTPt1 > 140)
    BDTPt1 = 140;
  if (BDTPt1 < 0)
    BDTPt1 = 0;

  return BDTPt1;
}

rad_to_deg()

double emtf::rad_to_deg ( double  rad )

inline

Definition at line 77 of file TrackTools.h.

References ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), and M_PI.

Referenced by calc_phi_glob_rad(), calc_theta_deg(), calc_theta_deg_from_eta(), PrimitiveConversion::convert_csc(), PrimitiveConversion::convert_dt(), PrimitiveConversion::convert_gem(), PrimitiveConversion::convert_me0(), PrimitiveConversion::convert_rpc(), RecHitProcessor::process(), and RecHitProcessor::processLook().

                                       {
    constexpr double factor = 180. / M_PI;
    return rad * factor;
  }

shuffle()

template<class bidiiter >

bidiiter emtf::shuffle	(	bidiiter	begin,
		bidiiter	end,
		size_t	num_random
	)

Definition at line 27 of file Utilities.h.

References HLT_2024v10_cff::distance, and edm::swap().

Referenced by createIOVlist::fillJson(), emtf::Forest::prepareRandomSubsample(), and InputFile::randomizeOrder().

                                                                    {
    // We will end up with the same elements in the collection except that
    // the first num_random elements will be randomized.

    size_t left = std::distance(begin, end);
    while (num_random--) {
      bidiiter r = begin;
      std::advance(r, rand() % left);
      std::swap(*begin, *r);
      ++begin;
      --left;
    }
    return begin;
  }

sort_uGMT_muons()

void emtf::sort_uGMT_muons

(

l1t::RegionalMuonCandBxCollection &

cands

)

Definition at line 111 of file MicroGMTConverter.cc.

References HLT_2024v10_cff::cands, l1t::emtf_neg, hltL1TEGammaFilteredCollectionProducer_cfi::maxBX, hltL1TEGammaFilteredCollectionProducer_cfi::minBX, ValidateTausOnZEEFastSim_cff::proc, and edm::swap().

Referenced by MicroGMTConverter::convert_all(), and l1t::stage2::EMTFCollections::~EMTFCollections().

                                                               {
    int minBX = cands.getFirstBX();
    int maxBX = cands.getLastBX();
    int emtfMinProc = 0;   // ME+ sector 1
    int emtfMaxProc = 11;  // ME- sector 6

    // New collection, sorted by processor to match uGMT unpacked order
    auto sortedCands = std::make_unique<l1t::RegionalMuonCandBxCollection>();
    sortedCands->clear();
    sortedCands->setBXRange(minBX, maxBX);
    for (int iBX = minBX; iBX <= maxBX; ++iBX) {
      for (int proc = emtfMinProc; proc <= emtfMaxProc; proc++) {
        for (l1t::RegionalMuonCandBxCollection::const_iterator cand = cands.begin(iBX); cand != cands.end(iBX);
             ++cand) {
          int cand_proc = cand->processor();
          if (cand->trackFinderType() == l1t::tftype::emtf_neg)
            cand_proc += 6;
          if (cand_proc != proc)
            continue;
          sortedCands->push_back(iBX, *cand);
        }
      }
    }

    // Return sorted collection
    std::swap(cands, *sortedCands);
    sortedCands.reset();
  }

sortNtupleByEvent()

void emtf::sortNtupleByEvent	(	const char *	ntuplename,
		const char *	filenametosort,
		const char *	outputfile
	)

Definition at line 98 of file Utilities.cc.

References f, DeadROC_duringRun::f2, mps_fire::i, createfilelist::int, and writedatasetfile::outputfile.

                                                                                                   {
  //TFile f("../../all_test_redux_post.root");
  TFile f(filenametosort);
  TNtuple* tree = (TNtuple*)f.Get(ntuplename);
  int nentries = (int)tree->GetEntries();
  //Drawing variable pz with no graphics option.
  //variable pz stored in array fV1 (see TTree::Draw)
  tree->Draw("Event", "", "goff");
  int* index = new int[nentries];
  //sort array containing pz in decreasing order
  //The array index contains the entry numbers in decreasing order
  TMath::Sort(nentries, tree->GetV1(), index);

  //open new file to store the sorted Tree
  //TFile f2("../../test_events_sorted.root","recreate");
  TFile f2(outputfile, "recreate");

  //Create an empty clone of the original tree
  TTree* tsorted = (TTree*)tree->CloneTree(0);
  for (int i = 0; i < nentries; i++) {
    tree->GetEntry(index[i]);
    tsorted->Fill();
  }
  tsorted->Write();
  delete[] index;
}

wrap_phi_deg()

double emtf::wrap_phi_deg ( double  deg )

inline

Definition at line 36 of file TrackTools.h.

Referenced by calc_phi_GMT_int(), and calc_phi_loc_deg_from_glob().

                                         {
    while (deg < -180.)
      deg += 360.;
    while (deg >= +180.)
      deg -= 360.;
    return deg;
  }

wrap_phi_rad()

double emtf::wrap_phi_rad ( double  rad )

inline

Definition at line 44 of file TrackTools.h.

References M_PI.

                                         {
    while (rad < -M_PI)
      rad += 2. * M_PI;
    while (rad >= +M_PI)
      rad -= 2. * M_PI;
    return rad;
  }

wrap_theta_deg()

double emtf::wrap_theta_deg ( double  deg )

inline

Definition at line 52 of file TrackTools.h.

References funct::abs().

Referenced by calc_eta_from_theta_deg().

                                           {
    deg = std::abs(deg);
    while (deg >= 180.)
      deg -= 180.;
    if (deg >= 180. / 2.)
      deg = 180. - deg;
    return deg;
  }

wrap_theta_rad()

double emtf::wrap_theta_rad ( double  rad )

inline

Definition at line 61 of file TrackTools.h.

References funct::abs(), and M_PI.

                                           {
    rad = std::abs(rad);
    while (rad >= M_PI)
      rad -= M_PI;
    if (rad >= M_PI / 2.)
      rad = M_PI - rad;
    return rad;
  }

Variable Documentation

MAX_ENDCAP

constexpr int emtf::MAX_ENDCAP = 2

Definition at line 46 of file Common.h.

Referenced by SectorProcessorShower::configure(), SectorProcessor::configure(), dump_fw_raw_input(), TrackFinder::process(), PrimitiveSelection::select_csc(), PrimitiveSelection::select_dt(), PrimitiveSelection::select_gem(), PrimitiveSelection::select_me0(), PrimitiveSelection::select_rpc(), and L1DummyProducer::SimpleDigi().

MAX_TRIGSECTOR

constexpr int emtf::MAX_TRIGSECTOR = 6

Definition at line 50 of file Common.h.

Referenced by SectorProcessorShower::configure(), SectorProcessor::configure(), dump_fw_raw_input(), TrackFinder::process(), PrimitiveSelection::select_csc(), PrimitiveSelection::select_dt(), PrimitiveSelection::select_gem(), PrimitiveSelection::select_me0(), and PrimitiveSelection::select_rpc().

MIN_ENDCAP

constexpr int emtf::MIN_ENDCAP = 1

Definition at line 45 of file Common.h.

Referenced by SectorProcessorShower::configure(), SectorProcessor::configure(), dump_fw_raw_input(), TrackFinder::process(), PrimitiveSelection::select_csc(), PrimitiveSelection::select_dt(), PrimitiveSelection::select_gem(), PrimitiveSelection::select_me0(), PrimitiveSelection::select_rpc(), and L1DummyProducer::SimpleDigi().

MIN_TRIGSECTOR

constexpr int emtf::MIN_TRIGSECTOR = 1

Definition at line 49 of file Common.h.

Referenced by SectorProcessorShower::configure(), SectorProcessor::configure(), dump_fw_raw_input(), TrackFinder::process(), PrimitiveSelection::select_csc(), PrimitiveSelection::select_dt(), PrimitiveSelection::select_gem(), PrimitiveSelection::select_me0(), and PrimitiveSelection::select_rpc().

NUM_FEATURES

constexpr int emtf::NUM_FEATURES = 29

Definition at line 68 of file Common.h.

Referenced by PtAssignmentEngineDxy::call_tensorflow_dxy().

NUM_PREDICTIONS

constexpr int emtf::NUM_PREDICTIONS = 2

Definition at line 69 of file Common.h.

Referenced by PtAssignmentEngineDxy::call_tensorflow_dxy().

NUM_SECTORS

constexpr int emtf::NUM_SECTORS = 12

Definition at line 51 of file Common.h.

NUM_STATION_PAIRS

constexpr int emtf::NUM_STATION_PAIRS = 6

Definition at line 59 of file Common.h.

Referenced by AngleCalculation::calculate_angles(), and PtAssignment::process().

NUM_STATIONS

constexpr int emtf::NUM_STATIONS = 4

Definition at line 58 of file Common.h.

Referenced by AngleCalculation::calculate_angles(), AngleCalculation::erase_tracks(), BestTrackSelection::process(), and PrimitiveMatching::process().

NUM_ZONE_HITS

constexpr int emtf::NUM_ZONE_HITS = 160

Definition at line 55 of file Common.h.

Referenced by PatternRecognition::process_single_zone().

NUM_ZONES

constexpr int emtf::NUM_ZONES = 4

Definition at line 54 of file Common.h.

Referenced by BestTrackSelection::cancel_multi_bx(), BestTrackSelection::cancel_one_bx(), PrimitiveConversion::get_zone_code(), AngleCalculation::process(), PrimitiveMatching::process(), and PatternRecognition::process().

ptscale

const double emtf::ptscale

Initial value:

= {0,    1.5,  2.0,  2.5,  3.0,  3.5,  4.0,   4.5,   5.0,  6.0,  7.0,
                                  8.0,  10.0, 12.0, 14.0, 16.0, 18.0, 20.0,  25.0,  30.0, 35.0, 40.0,
                                  45.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0, 120.0, 140.0}

Definition at line 29 of file Utilities.cc.

Referenced by L1TGMT::bookHistograms(), and processPrediction().

ptScale

const std::vector<double> emtf::ptScale

Definition at line 33 of file Utilities.cc.

Referenced by CSCTFSectorProcessor::CSCTFSectorProcessor(), CSCTFTrackBuilder::CSCTFTrackBuilder(), and CSCTFTrackProducer::produce().

twoJets_scale

const double emtf::twoJets_scale

Initial value:

= {
    0, 0.5, 1.0, 2.0, 3.0, 4.0, 5.0, 10.0, 20.0, 50.0, 100, 500, 1000, 5000, 7500, 50000}

Definition at line 36 of file Utilities.cc.

twoJetsScale

const std::vector< double > emtf::twoJetsScale

Initial value:

=
    std::vector<double>(twoJets_scale, twoJets_scale + sizeof twoJets_scale / sizeof twoJets_scale[0])

Definition at line 39 of file Utilities.cc.