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
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)
void	get_csc_max_strip_and_wire (int station, int ring, int &max_strip, int &max_wire)
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_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

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 91 of file TrackTools.h.

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

References bits, and 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 109 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 104 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 116 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 162 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 168 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 213 of file TrackTools.h.

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

References bits.

Referenced by calc_phi_GMT_rad().

calc_phi_GMT_int()

int emtf::calc_phi_GMT_int ( double  val )

inline

Definition at line 225 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 223 of file TrackTools.h.

{ return deg_to_rad(calc_phi_GMT_deg(bits)); }

References bits, calc_phi_GMT_deg(), and deg_to_rad().

calc_phi_loc_deg()

double emtf::calc_phi_loc_deg ( int  bits )

inline

Definition at line 172 of file TrackTools.h.

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

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

calc_phi_loc_deg_from_glob()

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

inline

Definition at line 191 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 197 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 205 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 178 of file TrackTools.h.

{ return deg_to_rad(calc_phi_loc_deg(bits)); }

References bits, calc_phi_loc_deg(), and deg_to_rad().

calc_pt()

double emtf::calc_pt ( int  bits )

inline

Definition at line 76 of file TrackTools.h.

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

References bits, and DiDispStaMuonMonitor_cfi::pt.

calc_pt_GMT()

int emtf::calc_pt_GMT ( double  val )

inline

Definition at line 82 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 137 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 158 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 124 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 139 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 146 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 132 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 153 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 130 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 64 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 5 of file DebugTools.cc.

                                                                                                         {
    // from interface/Common.h
    typedef L1TMuon::TriggerPrimitive TriggerPrimitive;
    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() == TriggerPrimitive::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() == TriggerPrimitive::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::TriggerPrimitive::kCSC, L1TMuon::TriggerPrimitive::kRPC, MAX_ENDCAP, MAX_TRIGSECTOR, MIN_ENDCAP, MIN_TRIGSECTOR, relativeConstraints::station, digitizers_cfi::strip, OrderedSet::t, and validateGeometry_cfg::valid.

Referenced by TrackFinder::process().

get_csc_max_strip_and_wire()

void emtf::get_csc_max_strip_and_wire	(	int	station,
		int	ring,
		int &	max_strip,
		int &	max_wire
	)

Definition at line 141 of file TrackTools.cc.

                                                                                        {
    max_strip = 0;                    // halfstrip
    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;
  }

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 69 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 begin, HLT_2018_cff::distance, end, alignCSCRings::r, rand(), and std::swap().

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

sort_uGMT_muons()

void emtf::sort_uGMT_muons

(

l1t::RegionalMuonCandBxCollection &

cands

)

Definition at line 109 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_2018_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 28 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 36 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 44 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 53 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_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 CSCTFPtLUT::CSCTFPtLUT(), 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.