L1TkMuCorrDynamicWindows Class Reference

#include <L1TkMuCorrDynamicWindows.h>

Public Types
typedef std::vector< L1TTTrackType >	L1TTTrackCollectionType
typedef TTTrack< Ref_Phase2TrackerDigi_ >	L1TTTrackType

Public Member Functions
std::vector< int >	find_match (const EMTFTrackCollection &l1mus, const L1TTTrackCollectionType &l1trks)
std::vector< int >	find_match_stub (const EMTFHitCollection &l1mus, const L1TTTrackCollectionType &l1trks, const int &station, bool requireBX0=true)
	L1TkMuCorrDynamicWindows (const std::vector< double > &bounds, TFile *fIn_theta, TFile *fIn_phi)
	L1TkMuCorrDynamicWindows (const std::vector< double > &bounds, TFile *fIn_theta, TFile *fIn_phi, TFile *fIn_theta_S1, TFile *fIn_phi_S1)
const float	max_trk_aeta ()
const float	max_trk_chi2 ()
const int	min_trk_nstubs ()
const float	min_trk_p ()
const int	n_trk_par ()
void	set_do_relax_factor (bool val)
void	set_do_trk_qual_presel (bool val)
void	set_max_trk_aeta (float val)
void	set_max_trk_chi2 (float val)
void	set_min_trk_nstubs (int val)
void	set_min_trk_p (float val)
void	set_n_trk_par (int val)
void	set_relaxation_pattern (float pt_start, float pt_end)
void	set_safety_factor (float sf_l, float sf_h)
void	set_safety_factor (float sf)
void	set_sf_initialrelax (float sf_l, float sf_h)
void	set_sf_initialrelax (float sf)
	~L1TkMuCorrDynamicWindows ()

Static Public Member Functions
static std::vector< double >	prepare_corr_bounds (const string &fname, const string &hname)

Private Member Functions
double	eta_to_theta (double x)
int	findBin (double val)
std::vector< int >	make_unique_coll (const unsigned int &l1musSize, const L1TTTrackCollectionType &l1trks, const std::vector< int > &matches)
double	sf_progressive (double x, double xstart, double xstop, double ystart, double ystop)
double	to_mpio2_pio2 (double x)

Private Attributes
std::vector< double >	bounds_
bool	do_relax_factor_
float	initial_sf_h_
float	initial_sf_l_
float	max_trk_aeta_
float	max_trk_chi2_
int	min_trk_nstubs_
float	min_trk_p_
int	nbins_
int	nTrkPars_
float	pt_end_
float	pt_start_
float	safety_factor_h_
float	safety_factor_l_
bool	track_qual_presel_
std::vector< MuMatchWindow >	wdws_phi_
std::vector< MuMatchWindow >	wdws_phi_S1_
std::vector< MuMatchWindow >	wdws_theta_
std::vector< MuMatchWindow >	wdws_theta_S1_

Detailed Description

Definition at line 24 of file L1TkMuCorrDynamicWindows.h.

Member Typedef Documentation

L1TTTrackCollectionType

typedef std::vector<L1TTTrackType> L1TkMuCorrDynamicWindows::L1TTTrackCollectionType

Definition at line 27 of file L1TkMuCorrDynamicWindows.h.

L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TkMuCorrDynamicWindows::L1TTTrackType

Definition at line 26 of file L1TkMuCorrDynamicWindows.h.

Constructor & Destructor Documentation

L1TkMuCorrDynamicWindows() [1/2]

L1TkMuCorrDynamicWindows::L1TkMuCorrDynamicWindows	(	const std::vector< double > &	bounds,
		TFile *	fIn_theta,
		TFile *	fIn_phi
	)

Definition at line 9 of file L1TkMuCorrDynamicWindows.cc.

References bounds_, Exception, cuy::ib, nbins_, set_do_relax_factor(), set_relaxation_pattern(), set_safety_factor(), set_sf_initialrelax(), AlCaHLTBitMon_QueryRunRegistry::string, cond::impl::to_string(), track_qual_presel_, wdws_phi_, and wdws_theta_.

    : wdws_theta_(bounds.size() - 1, MuMatchWindow()),
      wdws_phi_(bounds.size() - 1, MuMatchWindow()),
      wdws_theta_S1_(bounds.size() - 1, MuMatchWindow()),
      wdws_phi_S1_(bounds.size() - 1, MuMatchWindow()) {
  set_safety_factor(0.5);
  set_sf_initialrelax(0.0);
  set_relaxation_pattern(2.0, 6.0);
  set_do_relax_factor(true);

  track_qual_presel_ = true;

  nbins_ = bounds.size() - 1;
  bounds_ = bounds;

  // now load in memory the TF1 fits

  for (int ib = 0; ib < nbins_; ++ib) {
    std::string wdn;
    std::string nml;
    std::string nmh;
    TF1* fl;
    TF1* fh;

    // Station 2
    wdn = std::string("wdw_theta_") + std::to_string(ib + 1);
    nml = std::string("fit_low_") + std::to_string(ib + 1);
    nmh = std::string("fit_high_") + std::to_string(ib + 1);
    fl = (TF1*)fIn_theta->Get(nml.c_str());
    fh = (TF1*)fIn_theta->Get(nmh.c_str());
    if (fl == nullptr || fh == nullptr)
      throw cms::Exception("L1TkMuCorrDynamicWindows")
          << "TF1 named " << nml << " or " << nmh << " not found in file " << fIn_theta->GetName() << ".\n";
    wdws_theta_.at(ib).SetName(wdn);
    wdws_theta_.at(ib).SetLower(fl);
    wdws_theta_.at(ib).SetUpper(fh);

    wdn = std::string("wdw_phi_") + std::to_string(ib + 1);
    nml = std::string("fit_low_") + std::to_string(ib + 1);
    nmh = std::string("fit_high_") + std::to_string(ib + 1);
    fl = (TF1*)fIn_phi->Get(nml.c_str());
    fh = (TF1*)fIn_phi->Get(nmh.c_str());
    if (fl == nullptr || fh == nullptr)
      throw cms::Exception("L1TkMuCorrDynamicWindows")
          << "TF1 named " << nml << " or " << nmh << " not found in file " << fIn_theta->GetName() << ".\n";
    wdws_phi_.at(ib).SetName(wdn);
    wdws_phi_.at(ib).SetLower(fl);
    wdws_phi_.at(ib).SetUpper(fh);
  }
}

L1TkMuCorrDynamicWindows() [2/2]

L1TkMuCorrDynamicWindows::L1TkMuCorrDynamicWindows	(	const std::vector< double > &	bounds,
		TFile *	fIn_theta,
		TFile *	fIn_phi,
		TFile *	fIn_theta_S1,
		TFile *	fIn_phi_S1
	)

Definition at line 60 of file L1TkMuCorrDynamicWindows.cc.

References bounds_, Exception, cuy::ib, nbins_, set_do_relax_factor(), set_relaxation_pattern(), set_safety_factor(), set_sf_initialrelax(), AlCaHLTBitMon_QueryRunRegistry::string, cond::impl::to_string(), track_qual_presel_, wdws_phi_, wdws_phi_S1_, wdws_theta_, and wdws_theta_S1_.

    : wdws_theta_(bounds.size() - 1, MuMatchWindow()),
      wdws_phi_(bounds.size() - 1, MuMatchWindow()),
      wdws_theta_S1_(bounds.size() - 1, MuMatchWindow()),
      wdws_phi_S1_(bounds.size() - 1, MuMatchWindow()) {
  set_safety_factor(0.0);
  set_sf_initialrelax(0.0);
  set_relaxation_pattern(2.0, 6.0);
  set_do_relax_factor(true);

  track_qual_presel_ = true;

  nbins_ = bounds.size() - 1;
  bounds_ = bounds;

  // now load in memory the TF1 fits

  for (int ib = 0; ib < nbins_; ++ib) {
    std::string wdn;
    std::string nml;
    std::string nmh;
    TF1* fl;
    TF1* fh;

    // Station 2
    wdn = std::string("wdw_theta_") + std::to_string(ib + 1);
    nml = std::string("fit_low_") + std::to_string(ib + 1);
    nmh = std::string("fit_high_") + std::to_string(ib + 1);
    fl = (TF1*)fIn_theta->Get(nml.c_str());
    fh = (TF1*)fIn_theta->Get(nmh.c_str());
    if (fl == nullptr || fh == nullptr)
      throw cms::Exception("L1TkMuCorrDynamicWindows")
          << "TF1 named " << nml << " or " << nmh << " not found in file " << fIn_theta->GetName() << ".\n";
    wdws_theta_.at(ib).SetName(wdn);
    wdws_theta_.at(ib).SetLower(fl);
    wdws_theta_.at(ib).SetUpper(fh);

    wdn = std::string("wdw_phi_") + std::to_string(ib + 1);
    nml = std::string("fit_low_") + std::to_string(ib + 1);
    nmh = std::string("fit_high_") + std::to_string(ib + 1);
    fl = (TF1*)fIn_phi->Get(nml.c_str());
    fh = (TF1*)fIn_phi->Get(nmh.c_str());
    if (fl == nullptr || fh == nullptr)
      throw cms::Exception("L1TkMuCorrDynamicWindows")
          << "TF1 named " << nml << " or " << nmh << " not found in file " << fIn_theta->GetName() << ".\n";
    wdws_phi_.at(ib).SetName(wdn);
    wdws_phi_.at(ib).SetLower(fl);
    wdws_phi_.at(ib).SetUpper(fh);

    // Station 1 - MW's don't have to exist for TkMuon correlator
    // It is only needed for TkMuStub
    wdn = std::string("wdw_theta_") + std::to_string(ib + 1);
    nml = std::string("fit_low_") + std::to_string(ib + 1);
    nmh = std::string("fit_high_") + std::to_string(ib + 1);
    fl = (TF1*)fIn_theta_S1->Get(nml.c_str());
    fh = (TF1*)fIn_theta_S1->Get(nmh.c_str());
    if (fl == nullptr || fh == nullptr)
      throw cms::Exception("L1TkMuCorrDynamicWindows")
          << "TF1 named " << nml << " or " << nmh << " not found in file " << fIn_theta->GetName() << ".\n";
    wdws_theta_S1_.at(ib).SetName(wdn);
    wdws_theta_S1_.at(ib).SetLower(fl);
    wdws_theta_S1_.at(ib).SetUpper(fh);

    wdn = std::string("wdw_phi_") + std::to_string(ib + 1);
    nml = std::string("fit_low_") + std::to_string(ib + 1);
    nmh = std::string("fit_high_") + std::to_string(ib + 1);
    fl = (TF1*)fIn_phi_S1->Get(nml.c_str());
    fh = (TF1*)fIn_phi_S1->Get(nmh.c_str());
    if (fl == nullptr || fh == nullptr)
      throw cms::Exception("L1TkMuCorrDynamicWindows")
          << "TF1 named " << nml << " or " << nmh << " not found in file " << fIn_theta->GetName() << ".\n";
    wdws_phi_S1_.at(ib).SetName(wdn);
    wdws_phi_S1_.at(ib).SetLower(fl);
    wdws_phi_S1_.at(ib).SetUpper(fh);
  }
}

~L1TkMuCorrDynamicWindows()

L1TkMuCorrDynamicWindows::~L1TkMuCorrDynamicWindows ( )

inline

Definition at line 32 of file L1TkMuCorrDynamicWindows.h.

{}

Member Function Documentation

eta_to_theta()

double L1TkMuCorrDynamicWindows::eta_to_theta ( double  x )

inlineprivate

Definition at line 87 of file L1TkMuCorrDynamicWindows.h.

References JetChargeProducer_cfi::exp, and x.

Referenced by find_match(), and find_match_stub().

                                {
    //  give theta in rad
    return (2. * atan(exp(-1. * x)));
  }

find_match()

std::vector< int > L1TkMuCorrDynamicWindows::find_match	(	const EMTFTrackCollection &	l1mus,
		const L1TTTrackCollectionType &	l1trks
	)

Definition at line 155 of file L1TkMuCorrDynamicWindows.cc.

References funct::abs(), angle_units::operators::convertDegToRad(), reco::deltaPhi(), HLT_2022v12_cff::distance, do_relax_factor_, eta_to_theta(), findBin(), initial_sf_h_, initial_sf_l_, make_unique_coll(), muonTagProbeFilters_cff::matched, max_trk_aeta_, max_trk_chi2_, min_trk_nstubs_, min_trk_p_, MillePedeFileConverter_cfg::out, pt_end_, pt_start_, safety_factor_h_, safety_factor_l_, sf_progressive(), jetUpdater_cfi::sort, to_mpio2_pio2(), track_qual_presel_, wdws_phi_, and wdws_theta_.

                                                                                             {
  std::vector<int> out(l1trks.size());
  for (auto l1trkit = l1trks.begin(); l1trkit != l1trks.end(); ++l1trkit) {
    float trk_pt = l1trkit->momentum().perp();
    float trk_p = l1trkit->momentum().mag();
    float trk_aeta = std::abs(l1trkit->momentum().eta());
    float trk_theta = to_mpio2_pio2(eta_to_theta(l1trkit->momentum().eta()));
    float trk_phi = l1trkit->momentum().phi();
    int trk_charge = (l1trkit->rInv() > 0 ? 1 : -1);

    // porting some selections from the MuonTrackCorr finder
    // https://github.com/cms-l1t-offline/cmssw/blob/l1t-phase2-932-v1.6/L1Trigger/L1TTrackMatch/plugins/L1TkMuonProducer.cc#L264
    // in future, make preselections confiuguable
    bool reject_trk = false;
    if (trk_p < min_trk_p_)
      reject_trk = true;
    if (trk_aeta > max_trk_aeta_)
      reject_trk = true;
    if (track_qual_presel_) {
      float l1tk_chi2 = l1trkit->chi2();
      int l1tk_nstubs = l1trkit->getStubRefs().size();
      if (l1tk_chi2 >= max_trk_chi2_)
        reject_trk = true;
      if (l1tk_nstubs < min_trk_nstubs_)
        reject_trk = true;
    }

    int ibin = findBin(trk_aeta);

    std::vector<std::tuple<float, float, int>> matched;  // dtheta, dphi, idx
    // loop on muons to see which match
    for (auto l1muit = l1mus.begin(); l1muit != l1mus.end(); ++l1muit) {
      // match only muons in the central bx - as the track collection refers anyway to bx 0 only
      if (l1muit->BX() != 0)
        continue;

      // putting everything in rad
      float emtf_theta = to_mpio2_pio2(eta_to_theta(l1muit->Eta()));
      float emtf_phi = angle_units::operators::convertDegToRad(l1muit->Phi_glob());

      float dtheta = std::abs(emtf_theta - trk_theta);
      float dphi = reco::deltaPhi(emtf_phi, trk_phi);
      float adphi = std::abs(dphi);

      double sf_l;
      double sf_h;
      if (do_relax_factor_) {
        sf_l = sf_progressive(trk_pt, pt_start_, pt_end_, initial_sf_l_, safety_factor_l_);
        sf_h = sf_progressive(trk_pt, pt_start_, pt_end_, initial_sf_h_, safety_factor_h_);
      } else {
        sf_l = safety_factor_l_;
        sf_h = safety_factor_h_;
      }

      if (
          // emtf_theta * trk_theta > 0 &&
          dtheta > (1 - sf_l) * wdws_theta_.at(ibin).bound_low(trk_pt) &&
          dtheta <= (1 + sf_h) * wdws_theta_.at(ibin).bound_high(trk_pt) &&
          adphi > (1 - sf_l) * wdws_phi_.at(ibin).bound_low(trk_pt) &&
          adphi <= (1 + sf_h) * wdws_phi_.at(ibin).bound_high(trk_pt) && dphi * trk_charge < 0 &&  // sign requirement
          // rndm > 0.5
          true)
        matched.emplace_back(dtheta, adphi, std::distance(l1mus.begin(), l1muit));
    }

    if (reject_trk)
      matched.clear();  // quick fix - to be optimised to avoid the operations above

    if (matched.empty())
      out.at(std::distance(l1trks.begin(), l1trkit)) = -1;
    else {
      std::sort(matched.begin(), matched.end());  // closest in theta, then in phi
      out.at(std::distance(l1trks.begin(), l1trkit)) = std::get<2>(matched.at(0));
    }
  }

  // now convert out to a unique set
  return make_unique_coll(l1mus.size(), l1trks, out);
}

find_match_stub()

std::vector< int > L1TkMuCorrDynamicWindows::find_match_stub	(	const EMTFHitCollection &	l1mus,
		const L1TTTrackCollectionType &	l1trks,
		const int &	station,
		bool	requireBX0 = true
	)

Definition at line 236 of file L1TkMuCorrDynamicWindows.cc.

References funct::abs(), angle_units::operators::convertDegToRad(), reco::deltaPhi(), HLT_2022v12_cff::distance, do_relax_factor_, eta_to_theta(), findBin(), initial_sf_h_, initial_sf_l_, make_unique_coll(), muonTagProbeFilters_cff::matched, max_trk_aeta_, max_trk_chi2_, min_trk_nstubs_, min_trk_p_, MillePedeFileConverter_cfg::out, pt_end_, pt_start_, safety_factor_h_, safety_factor_l_, sf_progressive(), jetUpdater_cfi::sort, relativeConstraints::station, to_mpio2_pio2(), track_qual_presel_, wdws_phi_, wdws_phi_S1_, wdws_theta_, and wdws_theta_S1_.

                                                                            {
  std::vector<int> out(l1trks.size());
  for (auto l1trkit = l1trks.begin(); l1trkit != l1trks.end(); ++l1trkit) {
    float trk_pt = l1trkit->momentum().perp();
    float trk_p = l1trkit->momentum().mag();
    float trk_aeta = std::abs(l1trkit->momentum().eta());
    float trk_theta = to_mpio2_pio2(eta_to_theta(l1trkit->momentum().eta()));
    float trk_phi = l1trkit->momentum().phi();
    int trk_charge = (l1trkit->rInv() > 0 ? 1 : -1);

    // porting some selections from the MuonTrackCorr finder
    // https://github.com/cms-l1t-offline/cmssw/blob/l1t-phase2-932-v1.6/L1Trigger/L1TTrackMatch/plugins/L1TkMuonProducer.cc#L264
    // for future: make preselections confiuguable
    bool reject_trk = false;
    if (trk_p < min_trk_p_)
      reject_trk = true;
    if (trk_aeta > max_trk_aeta_)
      reject_trk = true;
    if (track_qual_presel_) {
      float l1tk_chi2 = l1trkit->chi2();
      int l1tk_nstubs = l1trkit->getStubRefs().size();
      if (l1tk_chi2 >= max_trk_chi2_)
        reject_trk = true;
      if (l1tk_nstubs < min_trk_nstubs_)
        reject_trk = true;
    }

    int ibin = findBin(trk_aeta);

    std::vector<std::tuple<float, float, int>> matched;  // dtheta, dphi, idx
    // loop on stubs to see which match
    for (auto l1muit = l1mus.begin(); l1muit != l1mus.end(); ++l1muit) {
      if (!(l1muit->Is_CSC() || l1muit->Is_RPC()))
        continue;

      int hit_station = l1muit->Station();
      // match only stubs in the central bx - as the track collection refers anyway to bx 0 only
      if (requireBX0 && l1muit->BX() != 0)
        continue;

      // allow only track matching to stubs from the given station, station= 1,2,3,4
      if (station < 5 && hit_station != station)
        continue;
      // in case of station=12 allow track matching to stubs from either station 1 or 2.
      else if (station == 12 && hit_station > 2)  // good for tkMuStub12
        continue;
      // in case of station=123 allow track matching to stubs from either station 1, 2, or 3.
      else if (station == 123 && hit_station > 3)  // good for tkMuStub123
        continue;
      // in case of station=1234 allow track matching to stubs from either station 1, 2, 3, or 4.
      else if (station == 1234 && hit_station > 4)  // good for tkMuStub1234
        continue;

      float emtf_theta = to_mpio2_pio2(eta_to_theta(l1muit->Eta_sim()));
      float emtf_phi = angle_units::operators::convertDegToRad(l1muit->Phi_sim());

      float dtheta = std::abs(emtf_theta - trk_theta);
      float dphi = reco::deltaPhi(emtf_phi, trk_phi);
      float adphi = std::abs(dphi);

      double sf_l;
      double sf_h;
      if (do_relax_factor_) {
        sf_l = sf_progressive(trk_pt, pt_start_, pt_end_, initial_sf_l_, safety_factor_l_);
        sf_h = sf_progressive(trk_pt, pt_start_, pt_end_, initial_sf_h_, safety_factor_h_);
      } else {
        sf_l = safety_factor_l_;
        sf_h = safety_factor_h_;
      }

      if (hit_station == 1 &&
          //if hit in station 1 use these  matching windows for checking

          dtheta > (1 - sf_l) * wdws_theta_S1_.at(ibin).bound_low(trk_pt) &&
          dtheta <= (1 + sf_h) * wdws_theta_S1_.at(ibin).bound_high(trk_pt) &&
          adphi > (1 - sf_l) * wdws_phi_S1_.at(ibin).bound_low(trk_pt) &&
          adphi <= (1 + sf_h) * wdws_phi_S1_.at(ibin).bound_high(trk_pt) &&
          dphi * trk_charge < 0 &&  // sign requirement
          true)
        matched.emplace_back(dtheta, adphi, std::distance(l1mus.begin(), l1muit));

      if (hit_station == 2 && dtheta > (1 - sf_l) * wdws_theta_.at(ibin).bound_low(trk_pt) &&
          dtheta <= (1 + sf_h) * wdws_theta_.at(ibin).bound_high(trk_pt) &&
          adphi > (1 - sf_l) * wdws_phi_.at(ibin).bound_low(trk_pt) &&
          adphi <= (1 + sf_h) * wdws_phi_.at(ibin).bound_high(trk_pt) && dphi * trk_charge < 0 &&  // sign requirement
          // rndm > 0.5
          true)
        matched.emplace_back(dtheta, adphi, std::distance(l1mus.begin(), l1muit));
    }

    if (reject_trk)
      matched.clear();  // quick fix - to be optimised to avoid the operations above

    if (matched.empty())
      out.at(std::distance(l1trks.begin(), l1trkit)) = -1;
    else {
      std::sort(matched.begin(), matched.end());  // closest in theta, then in phi
      out.at(std::distance(l1trks.begin(), l1trkit)) = std::get<2>(matched.at(0));
    }
  }

  // return out;

  // now convert out to a unique set
  auto unique_out = make_unique_coll(l1mus.size(), l1trks, out);

  return unique_out;
}

findBin()

int L1TkMuCorrDynamicWindows::findBin ( double  val )

private

Definition at line 138 of file L1TkMuCorrDynamicWindows.cc.

References bounds_, Exception, cuy::ib, nbins_, parallelization::uint, and heppy_batch::val.

Referenced by find_match(), and find_match_stub().

                                                {
  // not the most efficient, nor the most elegant implementation for now
  if (val < bounds_.at(0))
    return 0;
  if (val >= bounds_.back())
    return (nbins_ - 1);  // i.e. bounds_size() -2

  for (uint ib = 0; ib < bounds_.size() - 1; ++ib) {
    if (val >= bounds_.at(ib) && val < bounds_.at(ib + 1))
      return ib;
  }

  //"Something strange happened at val.
  throw cms::Exception("L1TkMuCorrDynamicWindows") << "Can't find bin.\n";
  return 0;
}

make_unique_coll()

std::vector< int > L1TkMuCorrDynamicWindows::make_unique_coll ( const unsigned int &  l1musSize, const L1TTTrackCollectionType &  l1trks, const std::vector< int > &  matches  )

private

Definition at line 348 of file L1TkMuCorrDynamicWindows.cc.

References oniaPATMuonsWithTrigger_cff::matches, nTrkPars_, MillePedeFileConverter_cfg::out, HLT_2022v12_cff::pt1, HLT_2022v12_cff::pt2, and jetUpdater_cfi::sort.

Referenced by find_match(), and find_match_stub().

                                                                                           {
  std::vector<int> out(matches.size(), -1);

  std::vector<std::vector<int>> macthed_to_emtf(l1musSize,
                                                std::vector<int>(0));  // one vector of matched trk idx per EMTF

  for (unsigned int itrack = 0; itrack < matches.size(); ++itrack) {
    int iemtf = matches.at(itrack);
    if (iemtf < 0)
      continue;
    macthed_to_emtf.at(iemtf).push_back(itrack);
  }

  std::function<bool(int, int, const L1TTTrackCollectionType&, int)> track_less_than_proto =
      [](int idx1, int idx2, const L1TTTrackCollectionType& l1trkcoll, int nTrackParams) {
        float pt1 = l1trkcoll.at(idx1).momentum().perp();
        float pt2 = l1trkcoll.at(idx2).momentum().perp();
        return (pt1 < pt2);
      };

  // // and binds to accept only 2 params
  std::function<bool(int, int)> track_less_than =
      std::bind(track_less_than_proto, std::placeholders::_1, std::placeholders::_2, l1trks, nTrkPars_);

  for (unsigned int iemtf = 0; iemtf < macthed_to_emtf.size(); ++iemtf) {
    std::vector<int>& thisv = macthed_to_emtf.at(iemtf);
    if (thisv.empty())
      continue;

    std::sort(thisv.begin(), thisv.end(), track_less_than);

    // copy to the output
    int best_trk = thisv.back();
    out.at(best_trk) = iemtf;
  }

  return out;
}

max_trk_aeta()

const float L1TkMuCorrDynamicWindows::max_trk_aeta ( )

inline

Definition at line 73 of file L1TkMuCorrDynamicWindows.h.

References max_trk_aeta_.

{ return max_trk_aeta_; }

max_trk_chi2()

const float L1TkMuCorrDynamicWindows::max_trk_chi2 ( )

inline

Definition at line 74 of file L1TkMuCorrDynamicWindows.h.

References max_trk_chi2_.

{ return max_trk_chi2_; }

min_trk_nstubs()

const int L1TkMuCorrDynamicWindows::min_trk_nstubs ( )

inline

Definition at line 75 of file L1TkMuCorrDynamicWindows.h.

References min_trk_nstubs_.

{ return min_trk_nstubs_; }

min_trk_p()

const float L1TkMuCorrDynamicWindows::min_trk_p ( )

inline

Definition at line 72 of file L1TkMuCorrDynamicWindows.h.

References min_trk_p_.

{ return min_trk_p_; }

n_trk_par()

const int L1TkMuCorrDynamicWindows::n_trk_par ( )

inline

Definition at line 71 of file L1TkMuCorrDynamicWindows.h.

References nTrkPars_.

{ return nTrkPars_; }

prepare_corr_bounds()

std::vector< double > L1TkMuCorrDynamicWindows::prepare_corr_bounds ( const string &  fname, const string &  hname  )

static

Definition at line 389 of file L1TkMuCorrDynamicWindows.cc.

References Exception, alignmentValidation::fname, and cuy::ib.

Referenced by L1TkMuonProducer::L1TkMuonProducer().

                                                                                                        {
  // find the boundaries of the match windoww
  TFile* fIn = TFile::Open(fname.c_str());
  if (fIn == nullptr) {
    throw cms::Exception("L1TkMuMantra") << "Can't find file " << fname << " to derive bounds.\n";
  }
  TH2* h_test = (TH2*)fIn->Get(hname.c_str());
  if (h_test == nullptr) {
    throw cms::Exception("L1TkMuCorrDynamicWindows")
        << "Can't find histo " << hname << " in file " << fname << " to derive bounds.\n";
  }

  int nbds = h_test->GetNbinsY() + 1;
  vector<double> bounds(nbds);
  for (int ib = 0; ib < nbds; ++ib) {
    bounds.at(ib) = h_test->GetYaxis()->GetBinLowEdge(ib + 1);
  }
  fIn->Close();
  return bounds;
}

set_do_relax_factor()

void L1TkMuCorrDynamicWindows::set_do_relax_factor ( bool  val )

inline

Definition at line 59 of file L1TkMuCorrDynamicWindows.h.

References do_relax_factor_, and heppy_batch::val.

Referenced by L1TkMuCorrDynamicWindows().

{ do_relax_factor_ = val; }

set_do_trk_qual_presel()

void L1TkMuCorrDynamicWindows::set_do_trk_qual_presel ( bool  val )

inline

Definition at line 61 of file L1TkMuCorrDynamicWindows.h.

References track_qual_presel_, and heppy_batch::val.

{ track_qual_presel_ = val; }

set_max_trk_aeta()

void L1TkMuCorrDynamicWindows::set_max_trk_aeta ( float  val )

inline

Definition at line 66 of file L1TkMuCorrDynamicWindows.h.

References max_trk_aeta_, and heppy_batch::val.

{ max_trk_aeta_ = val; }

set_max_trk_chi2()

void L1TkMuCorrDynamicWindows::set_max_trk_chi2 ( float  val )

inline

Definition at line 67 of file L1TkMuCorrDynamicWindows.h.

References max_trk_chi2_, and heppy_batch::val.

{ max_trk_chi2_ = val; }

set_min_trk_nstubs()

void L1TkMuCorrDynamicWindows::set_min_trk_nstubs ( int  val )

inline

Definition at line 68 of file L1TkMuCorrDynamicWindows.h.

References min_trk_nstubs_, and heppy_batch::val.

{ min_trk_nstubs_ = val; }

set_min_trk_p()

void L1TkMuCorrDynamicWindows::set_min_trk_p ( float  val )

inline

Definition at line 65 of file L1TkMuCorrDynamicWindows.h.

References min_trk_p_, and heppy_batch::val.

{ min_trk_p_ = val; }

set_n_trk_par()

void L1TkMuCorrDynamicWindows::set_n_trk_par ( int  val )

inline

Definition at line 64 of file L1TkMuCorrDynamicWindows.h.

References nTrkPars_, and heppy_batch::val.

{ nTrkPars_ = val; }

set_relaxation_pattern()

void L1TkMuCorrDynamicWindows::set_relaxation_pattern ( float  pt_start, float  pt_end  )

inline

Definition at line 53 of file L1TkMuCorrDynamicWindows.h.

References pt_end_, and pt_start_.

Referenced by L1TkMuCorrDynamicWindows().

                                                            {
    pt_start_ = pt_start;
    pt_end_ = pt_end;
  }

set_safety_factor() [1/2]

void L1TkMuCorrDynamicWindows::set_safety_factor ( float  sf_l, float  sf_h  )

inline

Definition at line 45 of file L1TkMuCorrDynamicWindows.h.

References safety_factor_h_, and safety_factor_l_.

Referenced by L1TkMuCorrDynamicWindows().

                                                 {
    safety_factor_l_ = sf_l;
    safety_factor_h_ = sf_h;
  }

set_safety_factor() [2/2]

void L1TkMuCorrDynamicWindows::set_safety_factor ( float  sf )

inline

Definition at line 57 of file L1TkMuCorrDynamicWindows.h.

References set_safety_factor().

Referenced by set_safety_factor().

{ set_safety_factor(sf, sf); }

set_sf_initialrelax() [1/2]

void L1TkMuCorrDynamicWindows::set_sf_initialrelax ( float  sf_l, float  sf_h  )

inline

Definition at line 49 of file L1TkMuCorrDynamicWindows.h.

References initial_sf_h_, and initial_sf_l_.

Referenced by L1TkMuCorrDynamicWindows().

                                                   {
    initial_sf_l_ = sf_l;
    initial_sf_h_ = sf_h;
  }

set_sf_initialrelax() [2/2]

void L1TkMuCorrDynamicWindows::set_sf_initialrelax ( float  sf )

inline

Definition at line 58 of file L1TkMuCorrDynamicWindows.h.

References set_sf_initialrelax().

Referenced by set_sf_initialrelax().

{ set_sf_initialrelax(sf, sf); }

sf_progressive()

double L1TkMuCorrDynamicWindows::sf_progressive ( double  x, double  xstart, double  xstop, double  ystart, double  ystop  )

inlineprivate

Definition at line 101 of file L1TkMuCorrDynamicWindows.h.

References x.

Referenced by find_match(), and find_match_stub().

                                                                                            {
    if (x < xstart)
      return ystart;
    if (x >= xstart && x < xstop)
      return ystart + (x - xstart) * (ystop - ystart) / (xstop - xstart);
    return ystop;
  }

to_mpio2_pio2()

double L1TkMuCorrDynamicWindows::to_mpio2_pio2 ( double  x )

inlineprivate

Definition at line 92 of file L1TkMuCorrDynamicWindows.h.

References M_PI, and x.

Referenced by find_match(), and find_match_stub().

                                 {
    //  put the angle in radians between -pi/2 and pi/2
    while (x >= 0.5 * M_PI)
      x -= M_PI;
    while (x < -0.5 * M_PI)
      x += M_PI;
    return x;
  }

Member Data Documentation

bounds_

std::vector<double> L1TkMuCorrDynamicWindows::bounds_

private

Definition at line 110 of file L1TkMuCorrDynamicWindows.h.

Referenced by findBin(), and L1TkMuCorrDynamicWindows().

do_relax_factor_

bool L1TkMuCorrDynamicWindows::do_relax_factor_

private

Definition at line 121 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), and set_do_relax_factor().

initial_sf_h_

float L1TkMuCorrDynamicWindows::initial_sf_h_

private

Definition at line 118 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), and set_sf_initialrelax().

initial_sf_l_

float L1TkMuCorrDynamicWindows::initial_sf_l_

private

Definition at line 117 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), and set_sf_initialrelax().

max_trk_aeta_

float L1TkMuCorrDynamicWindows::max_trk_aeta_

private

Definition at line 127 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), max_trk_aeta(), and set_max_trk_aeta().

max_trk_chi2_

float L1TkMuCorrDynamicWindows::max_trk_chi2_

private

Definition at line 128 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), max_trk_chi2(), and set_max_trk_chi2().

min_trk_nstubs_

int L1TkMuCorrDynamicWindows::min_trk_nstubs_

private

Definition at line 129 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), min_trk_nstubs(), and set_min_trk_nstubs().

min_trk_p_

float L1TkMuCorrDynamicWindows::min_trk_p_

private

Definition at line 126 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), min_trk_p(), and set_min_trk_p().

nbins_

int L1TkMuCorrDynamicWindows::nbins_

private

Definition at line 109 of file L1TkMuCorrDynamicWindows.h.

Referenced by findBin(), and L1TkMuCorrDynamicWindows().

nTrkPars_

int L1TkMuCorrDynamicWindows::nTrkPars_

private

Definition at line 125 of file L1TkMuCorrDynamicWindows.h.

Referenced by make_unique_coll(), n_trk_par(), and set_n_trk_par().

pt_end_

float L1TkMuCorrDynamicWindows::pt_end_

private

Definition at line 120 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), and set_relaxation_pattern().

pt_start_

float L1TkMuCorrDynamicWindows::pt_start_

private

Definition at line 119 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), and set_relaxation_pattern().

safety_factor_h_

float L1TkMuCorrDynamicWindows::safety_factor_h_

private

Definition at line 116 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), and set_safety_factor().

safety_factor_l_

float L1TkMuCorrDynamicWindows::safety_factor_l_

private

Definition at line 115 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), and set_safety_factor().

track_qual_presel_

bool L1TkMuCorrDynamicWindows::track_qual_presel_

private

Definition at line 122 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), L1TkMuCorrDynamicWindows(), and set_do_trk_qual_presel().

wdws_phi_

std::vector<MuMatchWindow> L1TkMuCorrDynamicWindows::wdws_phi_

private

Definition at line 112 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), and L1TkMuCorrDynamicWindows().

wdws_phi_S1_

std::vector<MuMatchWindow> L1TkMuCorrDynamicWindows::wdws_phi_S1_

private

Definition at line 114 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match_stub(), and L1TkMuCorrDynamicWindows().

wdws_theta_

std::vector<MuMatchWindow> L1TkMuCorrDynamicWindows::wdws_theta_

private

Definition at line 111 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match(), find_match_stub(), and L1TkMuCorrDynamicWindows().

wdws_theta_S1_

std::vector<MuMatchWindow> L1TkMuCorrDynamicWindows::wdws_theta_S1_

private

Definition at line 113 of file L1TkMuCorrDynamicWindows.h.

Referenced by find_match_stub(), and L1TkMuCorrDynamicWindows().