emtf Namespace Reference

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)
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 = 23
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 std::vector< double >	ptScale
const double	ptscale [31]
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.

                                                                                 {
    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;
  }

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

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

calc_eta()

double emtf::calc_eta ( int  bits )

inline

Definition at line 95 of file TrackTools.h.

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

References PVValHelper::eta.

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

calc_eta_from_theta_deg()

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

inline

Definition at line 113 of file TrackTools.h.

                                                                      {  // 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;
  }

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().

calc_eta_from_theta_rad()

double emtf::calc_eta_from_theta_rad ( double  theta_rad )

inline

Definition at line 108 of file TrackTools.h.

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

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

Referenced by calc_eta_from_theta_deg().

calc_eta_GMT()

int emtf::calc_eta_GMT ( double  val )

inline

Definition at line 120 of file TrackTools.h.

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

References heppy_batch::val.

calc_phi_glob_deg()

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

inline

Definition at line 166 of file TrackTools.h.

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

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().

calc_phi_glob_rad()

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

inline

Definition at line 172 of file TrackTools.h.

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

References calc_phi_glob_deg(), deg_to_rad(), and rad_to_deg().

calc_phi_GMT_deg()

double emtf::calc_phi_GMT_deg ( int  bits )

inline

Definition at line 217 of file TrackTools.h.

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

Referenced by calc_phi_GMT_rad().

calc_phi_GMT_int()

int emtf::calc_phi_GMT_int ( double  val )

inline

Definition at line 229 of file TrackTools.h.

                                          {  // 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;
  }

References heppy_batch::val, and wrap_phi_deg().

calc_phi_GMT_rad()

double emtf::calc_phi_GMT_rad ( int  bits )

inline

Definition at line 227 of file TrackTools.h.

{ return deg_to_rad(calc_phi_GMT_deg(bits)); }

References calc_phi_GMT_deg(), and deg_to_rad().

calc_phi_loc_deg()

double emtf::calc_phi_loc_deg ( int  bits )

inline

Definition at line 176 of file TrackTools.h.

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

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().

calc_phi_loc_deg_from_glob()

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

inline

Definition at line 195 of file TrackTools.h.

                                                                    {  // 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;
  }

References wrap_phi_deg().

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

calc_phi_loc_int()

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

inline

Definition at line 201 of file TrackTools.h.

                                                       {  // 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;
  }

References calc_phi_loc_deg_from_glob().

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

calc_phi_loc_int_rpc()

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

inline

Definition at line 209 of file TrackTools.h.

                                                           {  // 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;
  }

References calc_phi_loc_deg_from_glob().

Referenced by PrimitiveConversion::convert_rpc().

calc_phi_loc_rad()

double emtf::calc_phi_loc_rad ( int  bits )

inline

Definition at line 182 of file TrackTools.h.

{ return deg_to_rad(calc_phi_loc_deg(bits)); }

References calc_phi_loc_deg(), and deg_to_rad().

calc_pt()

double emtf::calc_pt ( int  bits )

inline

Definition at line 80 of file TrackTools.h.

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

References DiDispStaMuonMonitor_cfi::pt.

calc_pt_GMT()

int emtf::calc_pt_GMT ( double  val )

inline

Definition at line 86 of file TrackTools.h.

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

References heppy_batch::val.

calc_ring()

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

Definition at line 5 of file TrackTools.cc.

                                                    {
    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;
  }

References relativeConstraints::station.

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

calc_theta_deg()

double emtf::calc_theta_deg ( double  eta )

inline

Definition at line 141 of file TrackTools.h.

{ return rad_to_deg(calc_theta_rad(eta)); }

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

calc_theta_deg_from_eta()

double emtf::calc_theta_deg_from_eta ( double  eta )

inline

Definition at line 162 of file TrackTools.h.

{ return rad_to_deg(calc_theta_rad_from_eta(eta)); }

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

calc_theta_deg_from_int()

double emtf::calc_theta_deg_from_int ( int  theta_int )

inline

Definition at line 128 of file TrackTools.h.

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

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().

calc_theta_int()

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

inline

Definition at line 143 of file TrackTools.h.

                                                      {  // 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;
  }

References makeMuonMisalignmentScenario::endcap, and theta().

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

calc_theta_int_rpc()

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

inline

Definition at line 150 of file TrackTools.h.

                                                          {  // 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;
  }

References makeMuonMisalignmentScenario::endcap, and theta().

Referenced by PrimitiveConversion::convert_rpc().

calc_theta_rad()

double emtf::calc_theta_rad ( double  eta )

inline

Definition at line 136 of file TrackTools.h.

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

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

Referenced by calc_theta_deg().

calc_theta_rad_from_eta()

double emtf::calc_theta_rad_from_eta ( double  eta )

inline

Definition at line 157 of file TrackTools.h.

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

References PVValHelper::eta, and theta().

Referenced by calc_theta_deg_from_eta().

calc_theta_rad_from_int()

double emtf::calc_theta_rad_from_int ( int  theta_int )

inline

Definition at line 134 of file TrackTools.h.

{ return deg_to_rad(calc_theta_deg_from_int(theta_int)); }

References calc_theta_deg_from_int(), and deg_to_rad().

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.

                                                                              {
    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;
    }
  }

References relativeConstraints::station.

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

deg_to_rad()

double emtf::deg_to_rad ( double  deg )

inline

Definition at line 68 of file TrackTools.h.

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

References DQMScaleToClient_cfi::factor, 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().

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.

                                                                                                         {
    // 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
  }

References l1GtPatternGenerator_cfi::bx, relativeConstraints::chamber, gather_cfg::cout, makeMuonMisalignmentScenario::endcap, L1TMuon::kCSC, L1TMuon::kRPC, MAX_ENDCAP, MAX_TRIGSECTOR, MIN_ENDCAP, MIN_TRIGSECTOR, relativeConstraints::station, digitizers_cfi::strip, submitPVValidationJobs::t, and RunInfoPI::valid.

Referenced by TrackFinder::process().

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.

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

Referenced by PrimitiveSelection::select_csc().

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.

                                                                      {
    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);
  }

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

Referenced by PrimitiveSelection::select_csc().

get_trigger_csc_ID()

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

Definition at line 100 of file TrackTools.cc.

                                                             {
    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;
  }

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().

get_trigger_sector()

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

Definition at line 85 of file TrackTools.cc.

                                                             {
    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.
  }

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().

mergeNtuples()

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

Definition at line 88 of file Utilities.cc.

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

References writedatasetfile::outputfile.

numToStr()

template<typename T >

std::string emtf::numToStr

(

T

num

)

Definition at line 43 of file Utilities.h.

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

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

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

processPrediction()

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

Definition at line 46 of file Utilities.cc.

                                                                   {
  // 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;
}

References RPCpg::pts, and ptscale.

rad_to_deg()

double emtf::rad_to_deg ( double  rad )

inline

Definition at line 73 of file TrackTools.h.

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

References DQMScaleToClient_cfi::factor, 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().

shuffle()

template<class bidiiter >

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

Definition at line 27 of file Utilities.h.

                                                                    {
    // 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;
  }

References HLT_FULL_cff::distance, mps_fire::end, alignCSCRings::r, and std::swap().

Referenced by emtf::Forest::prepareRandomSubsample().

sort_uGMT_muons()

void emtf::sort_uGMT_muons

(

l1t::RegionalMuonCandBxCollection &

cands

)

Definition at line 111 of file MicroGMTConverter.cc.

                                                               {
    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();
  }

References HLT_FULL_cff::cands, l1t::emtf_neg, simKBmtfStubs_cfi::maxBX, simKBmtfStubs_cfi::minBX, ValidateTausOnZEEFastSim_cff::proc, and std::swap().

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

sortNtupleByEvent()

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

Definition at line 98 of file Utilities.cc.

                                                                                                   {
  //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;
}

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

wrap_phi_deg()

double emtf::wrap_phi_deg ( double  deg )

inline

Definition at line 32 of file TrackTools.h.

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

Referenced by calc_phi_GMT_int(), and calc_phi_loc_deg_from_glob().

wrap_phi_rad()

double emtf::wrap_phi_rad ( double  rad )

inline

Definition at line 40 of file TrackTools.h.

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

References M_PI.

wrap_theta_deg()

double emtf::wrap_theta_deg ( double  deg )

inline

Definition at line 48 of file TrackTools.h.

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

References funct::abs().

Referenced by calc_eta_from_theta_deg().

wrap_theta_rad()

double emtf::wrap_theta_rad ( double  rad )

inline

Definition at line 57 of file TrackTools.h.

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

References funct::abs(), and M_PI.

Variable Documentation

MAX_ENDCAP

constexpr int emtf::MAX_ENDCAP = 2

constexpr

Definition at line 46 of file Common.h.

Referenced by 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

constexpr

Definition at line 50 of file Common.h.

Referenced by 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

constexpr

Definition at line 45 of file Common.h.

Referenced by 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

constexpr

Definition at line 49 of file Common.h.

Referenced by 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 = 23

constexpr

Definition at line 68 of file Common.h.

Referenced by PtAssignmentEngineDxy::call_tensorflow_dxy().

NUM_PREDICTIONS

constexpr int emtf::NUM_PREDICTIONS = 2

constexpr

Definition at line 69 of file Common.h.

Referenced by PtAssignmentEngineDxy::call_tensorflow_dxy().

NUM_SECTORS

constexpr int emtf::NUM_SECTORS = 12

constexpr

Definition at line 51 of file Common.h.

NUM_STATION_PAIRS

constexpr int emtf::NUM_STATION_PAIRS = 6

constexpr

Definition at line 59 of file Common.h.

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

NUM_STATIONS

constexpr int emtf::NUM_STATIONS = 4

constexpr

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

constexpr

Definition at line 55 of file Common.h.

Referenced by PatternRecognition::process_single_zone().

NUM_ZONES

constexpr int emtf::NUM_ZONES = 4

constexpr

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 std::vector<double> emtf::ptScale

Definition at line 33 of file Utilities.cc.

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

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().

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.