Inheritance diagram for L1TkMuonProducer:

Classes
struct	PropState

Public Types
enum	AlgoType { kTP = 1, kDynamicWindows = 2, kMantra = 3 }

typedef std::vector< L1TTTrackType >	L1TTTrackCollectionType

typedef TTTrack< Ref_Phase2TrackerDigi_ >	L1TTTrackType

Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Member Functions
	L1TkMuonProducer (const edm::ParameterSet &)

	~L1TkMuonProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	build_tkMuons_from_idxs (TkMuonCollection &tkMuons, const std::vector< int > &matches, const edm::Handle< L1TTTrackCollectionType > &l1tksH, const edm::Handle< RegionalMuonCandBxCollection > &muonH, int detector) const

void	produce (edm::Event &, const edm::EventSetup &) override

std::vector< L1TkMuMantraDF::muon_df >	product_to_muvec (const EMTFTrackCollection &l1mus) const

std::vector< L1TkMuMantraDF::muon_df >	product_to_muvec (const RegionalMuonCandBxCollection &l1mtfs) const

std::vector< L1TkMuMantraDF::track_df >	product_to_trkvec (const L1TTTrackCollectionType &l1tks) const

PropState	propagateToGMT (const L1TTTrackType &l1tk) const

void	runOnMTFCollection_v1 (const edm::Handle< RegionalMuonCandBxCollection > &, const edm::Handle< L1TTTrackCollectionType > &, TkMuonCollection &tkMuons, const int detector) const

void	runOnMTFCollection_v2 (const edm::Handle< EMTFTrackCollection > &, const edm::Handle< L1TTTrackCollectionType > &, TkMuonCollection &tkMuons) const

double	sigmaEtaTP (const RegionalMuonCand &mu) const

double	sigmaPhiTP (const RegionalMuonCand &mu) const

Private Attributes
AlgoType	bmtfMatchAlgoVersion_

const edm::EDGetTokenT< RegionalMuonCandBxCollection >	bmtfToken_

float	chi2Max_

bool	correctGMTPropForTkZ_

float	dRMax_

std::unique_ptr< L1TkMuCorrDynamicWindows >	dwcorr_

AlgoType	emtfMatchAlgoVersion_

const edm::EDGetTokenT< EMTFTrackCollection >	emtfTCToken_

const edm::EDGetTokenT< RegionalMuonCandBxCollection >	emtfToken_

float	etaMax_

float	etaMin_

int	mantra_n_trk_par_

std::unique_ptr< L1TkMuMantra >	mantracorr_barr_

std::unique_ptr< L1TkMuMantra >	mantracorr_endc_

std::unique_ptr< L1TkMuMantra >	mantracorr_ovrl_

int	nStubsmin_

AlgoType	omtfMatchAlgoVersion_

const edm::EDGetTokenT< RegionalMuonCandBxCollection >	omtfToken_

float	pTMinTra_

const edm::EDGetTokenT< std::vector< TTTrack< Ref_Phase2TrackerDigi_ > > >	trackToken_

bool	use5ParameterFit_

bool	useTPMatchWindows_

float	zMax_

Detailed Description

Definition at line 45 of file L1TkMuonProducer.cc.

Member Typedef Documentation

◆ L1TTTrackCollectionType

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

Definition at line 48 of file L1TkMuonProducer.cc.

◆ L1TTTrackType

typedef TTTrack<Ref_Phase2TrackerDigi_> L1TkMuonProducer::L1TTTrackType

Definition at line 47 of file L1TkMuonProducer.cc.

Member Enumeration Documentation

◆ AlgoType

enum L1TkMuonProducer::AlgoType

Enumerator
kTP
kDynamicWindows
kMantra

Definition at line 61 of file L1TkMuonProducer.cc.

61 { kTP = 1, kDynamicWindows = 2, kMantra = 3 };

Constructor & Destructor Documentation

◆ L1TkMuonProducer()

L1TkMuonProducer::L1TkMuonProducer ( const edm::ParameterSet & iConfig )

explicit

Definition at line 129 of file L1TkMuonProducer.cc.

     : etaMin_((float)iConfig.getParameter<double>("ETAMIN")),
       etaMax_((float)iConfig.getParameter<double>("ETAMAX")),
       zMax_((float)iConfig.getParameter<double>("ZMAX")),
       chi2Max_((float)iConfig.getParameter<double>("CHI2MAX")),
       pTMinTra_((float)iConfig.getParameter<double>("PTMINTRA")),
       dRMax_((float)iConfig.getParameter<double>("DRmax")),
       nStubsmin_(iConfig.getParameter<int>("nStubsmin")),
       // --- mantra corr params
       mantra_n_trk_par_(iConfig.getParameter<int>("mantra_n_trk_par")),
       bmtfToken_(consumes<RegionalMuonCandBxCollection>(iConfig.getParameter<edm::InputTag>("L1BMTFInputTag"))),
       omtfToken_(consumes<RegionalMuonCandBxCollection>(iConfig.getParameter<edm::InputTag>("L1OMTFInputTag"))),
       emtfToken_(consumes<RegionalMuonCandBxCollection>(iConfig.getParameter<edm::InputTag>("L1EMTFInputTag"))),
       emtfTCToken_(consumes<EMTFTrackCollection>(iConfig.getParameter<edm::InputTag>("L1EMTFTrackCollectionInputTag"))),
       trackToken_(consumes<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > >(
           iConfig.getParameter<edm::InputTag>("L1TrackInputTag"))) {
   // --------------------- configuration of the muon algorithm type
  
   std::string bmtfMatchAlgoVersionString = iConfig.getParameter<std::string>("bmtfMatchAlgoVersion");
   std::transform(bmtfMatchAlgoVersionString.begin(),
                  bmtfMatchAlgoVersionString.end(),
                  bmtfMatchAlgoVersionString.begin(),
                  ::tolower);  // make lowercase
  
   std::string omtfMatchAlgoVersionString = iConfig.getParameter<std::string>("omtfMatchAlgoVersion");
   std::transform(omtfMatchAlgoVersionString.begin(),
                  omtfMatchAlgoVersionString.end(),
                  omtfMatchAlgoVersionString.begin(),
                  ::tolower);  // make lowercase
  
   std::string emtfMatchAlgoVersionString = iConfig.getParameter<std::string>("emtfMatchAlgoVersion");
   std::transform(emtfMatchAlgoVersionString.begin(),
                  emtfMatchAlgoVersionString.end(),
                  emtfMatchAlgoVersionString.begin(),
                  ::tolower);  // make lowercase
  
   if (bmtfMatchAlgoVersionString == "tp")
     bmtfMatchAlgoVersion_ = kTP;
   else if (bmtfMatchAlgoVersionString == "mantra")
     bmtfMatchAlgoVersion_ = kMantra;
   else
     throw cms::Exception("TkMuAlgoConfig")
         << "the ID " << bmtfMatchAlgoVersionString << " of the BMTF algo matcher passed is invalid\n";
  
   if (omtfMatchAlgoVersionString == "tp")
     omtfMatchAlgoVersion_ = kTP;
   else if (omtfMatchAlgoVersionString == "mantra")
     omtfMatchAlgoVersion_ = kMantra;
   else
     throw cms::Exception("TkMuAlgoConfig")
         << "the ID " << omtfMatchAlgoVersionString << " of the OMTF algo matcher passed is invalid\n";
  
   if (emtfMatchAlgoVersionString == "tp")
     emtfMatchAlgoVersion_ = kTP;
   else if (emtfMatchAlgoVersionString == "dynamicwindows")
     emtfMatchAlgoVersion_ = kDynamicWindows;
   else if (emtfMatchAlgoVersionString == "mantra")
     emtfMatchAlgoVersion_ = kMantra;
   else
     throw cms::Exception("TkMuAlgoConfig")
         << "the ID " << emtfMatchAlgoVersionString << " of the EMTF algo matcher passed is invalid\n";
  
   correctGMTPropForTkZ_ = iConfig.getParameter<bool>("correctGMTPropForTkZ");
  
   use5ParameterFit_ = iConfig.getParameter<bool>("use5ParameterFit");
   useTPMatchWindows_ = iConfig.getParameter<bool>("useTPMatchWindows");
   produces<TkMuonCollection>();
  
   // initializations
   if (emtfMatchAlgoVersion_ == kDynamicWindows) {
     // FIXME: to merge eventually into an unique file with bith phi and theta boundaries
     std::string fIn_bounds_name = iConfig.getParameter<edm::FileInPath>("emtfcorr_boundaries").fullPath();
     std::string fIn_theta_name = iConfig.getParameter<edm::FileInPath>("emtfcorr_theta_windows").fullPath();
     std::string fIn_phi_name = iConfig.getParameter<edm::FileInPath>("emtfcorr_phi_windows").fullPath();
     const auto& bounds = L1TkMuCorrDynamicWindows::prepare_corr_bounds(fIn_bounds_name, "h_dphi_l");
     TFile* fIn_theta = TFile::Open(fIn_theta_name.c_str());
     TFile* fIn_phi = TFile::Open(fIn_phi_name.c_str());
     dwcorr_ = std::make_unique<L1TkMuCorrDynamicWindows>(bounds, fIn_theta, fIn_phi);
  
     // files can be closed since the correlator code clones the TF1s
     fIn_theta->Close();
     fIn_phi->Close();
  
     // FIXME: more initialisation using the parameters passed from the cfg
     dwcorr_->set_safety_factor(iConfig.getParameter<double>("final_window_factor"));
     dwcorr_->set_sf_initialrelax(iConfig.getParameter<double>("initial_window_factor"));
  
     dwcorr_->set_relaxation_pattern(iConfig.getParameter<double>("pt_start_relax"),
                                     iConfig.getParameter<double>("pt_end_relax"));
  
     dwcorr_->set_do_relax_factor(iConfig.getParameter<bool>("do_relax_factors"));
  
     dwcorr_->set_n_trk_par(iConfig.getParameter<int>("n_trk_par"));
     dwcorr_->set_min_trk_p(iConfig.getParameter<double>("min_trk_p"));
     dwcorr_->set_max_trk_aeta(iConfig.getParameter<double>("max_trk_aeta"));
     dwcorr_->set_max_trk_chi2(iConfig.getParameter<double>("max_trk_chi2"));
     dwcorr_->set_min_trk_nstubs(iConfig.getParameter<int>("min_trk_nstubs"));
     dwcorr_->set_do_trk_qual_presel(true);
   }
  
   if (bmtfMatchAlgoVersion_ == kMantra) {
     std::string fIn_bounds_name = iConfig.getParameter<edm::FileInPath>("mantra_bmtfcorr_boundaries").fullPath();
     std::string fIn_theta_name = iConfig.getParameter<edm::FileInPath>("mantra_bmtfcorr_theta_windows").fullPath();
     std::string fIn_phi_name = iConfig.getParameter<edm::FileInPath>("mantra_bmtfcorr_phi_windows").fullPath();
  
     const auto& bounds = L1TkMuMantra::prepare_corr_bounds(fIn_bounds_name, "h_dphi_l");
  
     TFile* fIn_theta = TFile::Open(fIn_theta_name.c_str());
     TFile* fIn_phi = TFile::Open(fIn_phi_name.c_str());
  
     mantracorr_barr_ = std::make_unique<L1TkMuMantra>(bounds, fIn_theta, fIn_phi, "mantra_barrel");
  
     fIn_theta->Close();
     fIn_phi->Close();
  
     // settings : NB : now hardcoded, could be read from cfg
     mantracorr_barr_->set_safety_factor(0.5, 0.5);
     mantracorr_barr_->setArbitrationType("MaxPt");
   }
  
   if (omtfMatchAlgoVersion_ == kMantra) {
     std::string fIn_bounds_name = iConfig.getParameter<edm::FileInPath>("mantra_omtfcorr_boundaries").fullPath();
     std::string fIn_theta_name = iConfig.getParameter<edm::FileInPath>("mantra_omtfcorr_theta_windows").fullPath();
     std::string fIn_phi_name = iConfig.getParameter<edm::FileInPath>("mantra_omtfcorr_phi_windows").fullPath();
  
     const auto& bounds = L1TkMuMantra::prepare_corr_bounds(fIn_bounds_name, "h_dphi_l");
  
     TFile* fIn_theta = TFile::Open(fIn_theta_name.c_str());
     TFile* fIn_phi = TFile::Open(fIn_phi_name.c_str());
  
     mantracorr_ovrl_ = std::make_unique<L1TkMuMantra>(bounds, fIn_theta, fIn_phi, "mantra_overlap");
  
     fIn_theta->Close();
     fIn_phi->Close();
  
     // settings : NB : now hardcoded, could be read from cfg
     mantracorr_ovrl_->set_safety_factor(0.5, 0.5);
     mantracorr_ovrl_->setArbitrationType("MaxPt");
   }
  
   if (emtfMatchAlgoVersion_ == kMantra) {
     std::string fIn_bounds_name = iConfig.getParameter<edm::FileInPath>("mantra_emtfcorr_boundaries").fullPath();
     std::string fIn_theta_name = iConfig.getParameter<edm::FileInPath>("mantra_emtfcorr_theta_windows").fullPath();
     std::string fIn_phi_name = iConfig.getParameter<edm::FileInPath>("mantra_emtfcorr_phi_windows").fullPath();
  
     const auto& bounds = L1TkMuMantra::prepare_corr_bounds(fIn_bounds_name, "h_dphi_l");
  
     TFile* fIn_theta = TFile::Open(fIn_theta_name.c_str());
     TFile* fIn_phi = TFile::Open(fIn_phi_name.c_str());
  
     mantracorr_endc_ = std::make_unique<L1TkMuMantra>(bounds, fIn_theta, fIn_phi, "mantra_endcap");
  
     fIn_theta->Close();
     fIn_phi->Close();
  
     // settings : NB : now hardcoded, could be read from cfg
     mantracorr_endc_->set_safety_factor(0.5, 0.5);
     mantracorr_endc_->setArbitrationType("MaxPt");
   }
 }

References bmtfMatchAlgoVersion_, correctGMTPropForTkZ_, dwcorr_, emtfMatchAlgoVersion_, Exception, contentValuesFiles::fullPath, edm::ParameterSet::getParameter(), kDynamicWindows, kMantra, kTP, mantracorr_barr_, mantracorr_endc_, mantracorr_ovrl_, omtfMatchAlgoVersion_, L1TkMuCorrDynamicWindows::prepare_corr_bounds(), L1TkMuMantra::prepare_corr_bounds(), AlCaHLTBitMon_QueryRunRegistry::string, HcalDetIdTransform::transform(), use5ParameterFit_, and useTPMatchWindows_.

◆ ~L1TkMuonProducer()

L1TkMuonProducer::~L1TkMuonProducer ( )

override

Definition at line 290 of file L1TkMuonProducer.cc.

290 {}

Member Function Documentation

◆ build_tkMuons_from_idxs()

void L1TkMuonProducer::build_tkMuons_from_idxs	(	TkMuonCollection &	tkMuons,
		const std::vector< int > &	matches,
		const edm::Handle< L1TTTrackCollectionType > &	l1tksH,
		const edm::Handle< RegionalMuonCandBxCollection > &	muonH,
		int	detector
	)		const

private

Definition at line 647 of file L1TkMuonProducer.cc.

                                                                    {
   for (uint imatch = 0; imatch < matches.size(); ++imatch) {
     int match_trk_idx = matches[imatch];
     if (match_trk_idx < 0)
       continue;  // this muon was not matched to any candidate
  
     // take properties of the track
     const L1TTTrackType& matchTk = (*l1tksH.product())[match_trk_idx];
     const auto& p3 = matchTk.momentum();
     const auto& tkv3 = matchTk.POCA();
     float p4e = sqrt(mu_mass * mu_mass + p3.mag2());
     math::XYZTLorentzVector l1tkp4(p3.x(), p3.y(), p3.z(), p4e);
  
     edm::Ptr<L1TTTrackType> l1tkPtr(l1tksH, match_trk_idx);
     auto l1muRef = muonH.isValid() ? edm::Ref<RegionalMuonCandBxCollection>(muonH, imatch)
                                    : edm::Ref<RegionalMuonCandBxCollection>();
  
     float trkisol = -999;
     TkMuon l1tkmu(l1tkp4, l1muRef, l1tkPtr, trkisol);
     l1tkmu.setTrackCurvature(matchTk.rInv());
     l1tkmu.setTrkzVtx((float)tkv3.z());
     l1tkmu.setMuonDetector(detector);
     tkMuons.push_back(l1tkmu);
   }
   return;
 }

References hgcalTestNeighbor_cfi::detector, edm::HandleBase::isValid(), patCandidatesForDimuonsSequences_cff::matches, TTTrack< T >::momentum(), mu_mass, p3, TTTrack< T >::POCA(), edm::Handle< T >::product(), TTTrack< T >::rInv(), l1t::TkMuon::setMuonDetector(), l1t::TkMuon::setTrackCurvature(), l1t::TkMuon::setTrkzVtx(), mathSSE::sqrt(), l1pfProducer_cfi::tkMuons, and parallelization::uint.

Referenced by produce().

◆ fillDescriptions()

void L1TkMuonProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 522 of file L1TkMuonProducer.cc.

                                                                                   {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

References edm::ConfigurationDescriptions::addDefault(), and submitPVResolutionJobs::desc.

◆ produce()

void L1TkMuonProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 293 of file L1TkMuonProducer.cc.

                                                                             {
   // the L1Mu objects
   edm::Handle<RegionalMuonCandBxCollection> l1bmtfH;
   edm::Handle<RegionalMuonCandBxCollection> l1omtfH;
   edm::Handle<RegionalMuonCandBxCollection> l1emtfH;
   edm::Handle<EMTFTrackCollection> l1emtfTCH;
  
   iEvent.getByToken(bmtfToken_, l1bmtfH);
   iEvent.getByToken(omtfToken_, l1omtfH);
   iEvent.getByToken(emtfToken_, l1emtfH);
   iEvent.getByToken(emtfTCToken_, l1emtfTCH);
  
   // the L1Tracks
   edm::Handle<L1TTTrackCollectionType> l1tksH;
   iEvent.getByToken(trackToken_, l1tksH);
  
   TkMuonCollection oc_bmtf_tkmuon;
   TkMuonCollection oc_omtf_tkmuon;
   TkMuonCollection oc_emtf_tkmuon;
  
   std::vector<L1TkMuMantraDF::track_df> mantradf_tracks;  // if needed, just encode once for all trk finders
   if (bmtfMatchAlgoVersion_ == kMantra || omtfMatchAlgoVersion_ == kMantra || emtfMatchAlgoVersion_ == kMantra) {
     mantradf_tracks = product_to_trkvec(*l1tksH.product());
   }
  
   // process each of the MTF collections separately! -- we don't want to filter the muons
  
   // ----------------------------------------------------- barrel
   if (bmtfMatchAlgoVersion_ == kTP)
     runOnMTFCollection_v1(l1bmtfH, l1tksH, oc_bmtf_tkmuon, barrel_MTF_region);
   else if (bmtfMatchAlgoVersion_ == kMantra) {
     const auto& muons = product_to_muvec(*l1bmtfH.product());
     const auto& match_idx = mantracorr_barr_->find_match(mantradf_tracks, muons);
     build_tkMuons_from_idxs(oc_bmtf_tkmuon, match_idx, l1tksH, l1bmtfH, barrel_MTF_region);
   } else
     throw cms::Exception("TkMuAlgoConfig") << " barrel : trying to run an invalid algorithm version "
                                            << bmtfMatchAlgoVersion_ << " (this should never happen)\n";
  
   // ----------------------------------------------------- overlap
   if (omtfMatchAlgoVersion_ == kTP)
     runOnMTFCollection_v1(l1omtfH, l1tksH, oc_omtf_tkmuon, overlap_MTF_region);
   else if (omtfMatchAlgoVersion_ == kMantra) {
     const auto& muons = product_to_muvec(*l1omtfH.product());
     const auto& match_idx = mantracorr_ovrl_->find_match(mantradf_tracks, muons);
     build_tkMuons_from_idxs(oc_omtf_tkmuon, match_idx, l1tksH, l1omtfH, overlap_MTF_region);
   } else
     throw cms::Exception("TkMuAlgoConfig") << " overlap : trying to run an invalid algorithm version "
                                            << omtfMatchAlgoVersion_ << " (this should never happen)\n";
  
   // ----------------------------------------------------- endcap
   if (emtfMatchAlgoVersion_ == kTP)
     runOnMTFCollection_v1(l1emtfH, l1tksH, oc_emtf_tkmuon, endcap_MTF_region);
   else if (emtfMatchAlgoVersion_ == kDynamicWindows)
     runOnMTFCollection_v2(l1emtfTCH, l1tksH, oc_emtf_tkmuon);
   else if (emtfMatchAlgoVersion_ == kMantra) {
     const auto& muons = product_to_muvec(*l1emtfTCH.product());
     const auto& match_idx = mantracorr_endc_->find_match(mantradf_tracks, muons);
     //for the TkMu that were built from a EMTFCollection - do not produce a valid muon ref
     edm::Handle<RegionalMuonCandBxCollection> invalidMuonH;
     build_tkMuons_from_idxs(oc_emtf_tkmuon, match_idx, l1tksH, invalidMuonH, endcap_MTF_region);
   } else
     throw cms::Exception("TkMuAlgoConfig") << "endcap : trying to run an invalid algorithm version "
                                            << emtfMatchAlgoVersion_ << " (this should never happen)\n";
  
   // now combine all trk muons into a single output collection!
   auto oc_tkmuon = std::make_unique<TkMuonCollection>();
   for (const auto& p : {oc_bmtf_tkmuon, oc_omtf_tkmuon, oc_emtf_tkmuon}) {
     oc_tkmuon->insert(oc_tkmuon->end(), p.begin(), p.end());
   }
  
   // put the new track+muon objects in the event!
   iEvent.put(std::move(oc_tkmuon));
 };

References barrel_MTF_region, bmtfMatchAlgoVersion_, bmtfToken_, build_tkMuons_from_idxs(), emtfMatchAlgoVersion_, emtfTCToken_, emtfToken_, endcap_MTF_region, Exception, iEvent, kDynamicWindows, kMantra, kTP, mantracorr_barr_, mantracorr_endc_, mantracorr_ovrl_, eostools::move(), PDWG_BPHSkim_cff::muons, omtfMatchAlgoVersion_, omtfToken_, overlap_MTF_region, AlCaHLTBitMon_ParallelJobs::p, edm::Handle< T >::product(), product_to_muvec(), product_to_trkvec(), runOnMTFCollection_v1(), runOnMTFCollection_v2(), and trackToken_.

◆ product_to_muvec() [1/2]

std::vector< L1TkMuMantraDF::muon_df > L1TkMuonProducer::product_to_muvec ( const EMTFTrackCollection & l1mus ) const

private

Definition at line 634 of file L1TkMuonProducer.cc.

                                                                                                           {
   std::vector<L1TkMuMantraDF::muon_df> result(l1mus.size());
   for (uint imu = 0; imu < l1mus.size(); ++imu) {
     auto& mu = l1mus[imu];
     result[imu].pt = mu.Pt();
     result[imu].eta = mu.Eta();
     result[imu].theta = L1TkMuMantra::to_mpio2_pio2(L1TkMuMantra::eta_to_theta(mu.Eta()));
     result[imu].phi = angle_units::operators::convertDegToRad(mu.Phi_glob());
     result[imu].charge = mu.Charge();
   }
   return result;
 }

References angle_units::operators::convertDegToRad(), L1TkMuMantra::eta_to_theta(), amptDefaultParameters_cff::mu, mps_fire::result, L1TkMuMantra::to_mpio2_pio2(), and parallelization::uint.

◆ product_to_muvec() [2/2]

std::vector< L1TkMuMantraDF::muon_df > L1TkMuonProducer::product_to_muvec ( const RegionalMuonCandBxCollection & l1mtfs ) const

private

Definition at line 617 of file L1TkMuonProducer.cc.

                                                       {
   std::vector<L1TkMuMantraDF::muon_df> result;
   for (auto l1mu = l1mtfs.begin(0); l1mu != l1mtfs.end(0); ++l1mu)  // considering BX = 0 only
   {
     L1TkMuMantraDF::muon_df this_mu;
     this_mu.pt = l1mu->hwPt() * 0.5;
     this_mu.eta = l1mu->hwEta() * eta_scale;
     this_mu.theta = L1TkMuMantra::to_mpio2_pio2(L1TkMuMantra::eta_to_theta(l1mu->hwEta() * eta_scale));
     this_mu.phi =
         MicroGMTConfiguration::calcGlobalPhi(l1mu->hwPhi(), l1mu->trackFinderType(), l1mu->processor()) * phi_scale;
     this_mu.charge = (l1mu->hwSign() == 0 ? 1 : -1);  // charge sign bit (charge = (-1)^(sign))
     result.push_back(this_mu);
   }
   return result;
 }

References BXVector< T >::begin(), L1TkMuMantraDF::muon_df::charge, BXVector< T >::end(), L1TkMuMantraDF::muon_df::eta, eta_scale, L1TkMuMantra::eta_to_theta(), L1TkMuMantraDF::muon_df::phi, phi_scale, L1TkMuMantraDF::muon_df::pt, mps_fire::result, L1TkMuMantraDF::muon_df::theta, and L1TkMuMantra::to_mpio2_pio2().

Referenced by produce().

◆ product_to_trkvec()

std::vector< L1TkMuMantraDF::track_df > L1TkMuonProducer::product_to_trkvec ( const L1TTTrackCollectionType & l1tks ) const

private

Definition at line 600 of file L1TkMuonProducer.cc.

                                                                                                                 {
   std::vector<L1TkMuMantraDF::track_df> result(l1tks.size());
   for (uint itrk = 0; itrk < l1tks.size(); ++itrk) {
     auto& trk = l1tks[itrk];
  
     result[itrk].pt = trk.momentum().perp();
     result[itrk].eta = trk.momentum().eta();
     result[itrk].theta = L1TkMuMantra::to_mpio2_pio2(L1TkMuMantra::eta_to_theta(trk.momentum().eta()));
     result[itrk].phi = trk.momentum().phi();
     result[itrk].nstubs = trk.getStubRefs().size();
     result[itrk].chi2 = trk.chi2();
     result[itrk].charge = (trk.rInv() > 0 ? 1 : -1);
   }
  
   return result;
 }

References L1TkMuMantra::eta_to_theta(), mps_fire::result, L1TkMuMantra::to_mpio2_pio2(), and parallelization::uint.

Referenced by produce().

◆ propagateToGMT()

L1TkMuonProducer::PropState L1TkMuonProducer::propagateToGMT ( const L1TTTrackType & l1tk ) const

private

Definition at line 530 of file L1TkMuonProducer.cc.

                                                                                                         {
   auto p3 = tk.momentum();
   float tk_pt = p3.perp();
   float tk_p = p3.mag();
   float tk_eta = p3.eta();
   float tk_aeta = std::abs(tk_eta);
   float tk_phi = p3.phi();
   float tk_q = tk.rInv() > 0 ? 1. : -1.;
   float tk_z = tk.POCA().z();
   if (!correctGMTPropForTkZ_)
     tk_z = 0;
  
   L1TkMuonProducer::PropState dest;
   if (tk_p < min_mu_propagator_p)
     return dest;
   if (tk_aeta < 1.1 && tk_pt < min_mu_propagator_barrel_pT)
     return dest;
   if (tk_aeta > max_mu_propagator_eta)
     return dest;
  
   //0th order:
   dest.valid = true;
  
   float dzCorrPhi = 1.;
   float deta = 0;
   float etaProp = tk_aeta;
  
   if (tk_aeta < 1.1) {
     etaProp = 1.1;
     deta = tk_z / 550. / cosh(tk_aeta);
   } else {
     float delta = tk_z / 850.;  //roughly scales as distance to 2nd station
     if (tk_eta > 0)
       delta *= -1;
     dzCorrPhi = 1. + delta;
  
     float zOzs = tk_z / 850.;
     if (tk_eta > 0)
       deta = zOzs / (1. - zOzs);
     else
       deta = zOzs / (1. + zOzs);
     deta = deta * tanh(tk_eta);
   }
   float resPhi = tk_phi - 1.464 * tk_q * cosh(1.7) / cosh(etaProp) / tk_pt * dzCorrPhi - M_PI / 144.;
   resPhi = reco::reduceRange(resPhi);
  
   dest.eta = tk_eta + deta;
   dest.phi = resPhi;
   dest.pt = tk_pt;  //not corrected for eloss
  
   dest.sigmaEta = 0.100 / tk_pt;  //multiple scattering term
   dest.sigmaPhi = 0.106 / tk_pt;  //need a better estimate for these
   return dest;
 }

References funct::abs(), correctGMTPropForTkZ_, dumpMFGeometry_cfg::delta, mps_fire::dest, M_PI, max_mu_propagator_eta, min_mu_propagator_barrel_pT, min_mu_propagator_p, TTTrack< T >::momentum(), p3, TTTrack< T >::POCA(), reco::reduceRange(), TTTrack< T >::rInv(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by runOnMTFCollection_v1().

◆ runOnMTFCollection_v1()

void L1TkMuonProducer::runOnMTFCollection_v1	(	const edm::Handle< RegionalMuonCandBxCollection > &	muonH,
		const edm::Handle< L1TTTrackCollectionType > &	l1tksH,
		TkMuonCollection &	tkMuons,
		const int	detector
	)		const

private

Definition at line 367 of file L1TkMuonProducer.cc.

                                                                        {
   const L1TTTrackCollectionType& l1tks = (*l1tksH.product());
   const RegionalMuonCandBxCollection& l1mtfs = (*muonH.product());
  
   int imu = 0;
   for (auto l1mu = l1mtfs.begin(0); l1mu != l1mtfs.end(0); ++l1mu) {  // considering BX = only
  
     edm::Ref<RegionalMuonCandBxCollection> l1muRef(muonH, imu);
     imu++;
  
     float l1mu_eta = l1mu->hwEta() * eta_scale;
     // get the global phi
     float l1mu_phi =
         MicroGMTConfiguration::calcGlobalPhi(l1mu->hwPhi(), l1mu->trackFinderType(), l1mu->processor()) * phi_scale;
  
     float l1mu_feta = std::abs(l1mu_eta);
     if (l1mu_feta < etaMin_)
       continue;
     if (l1mu_feta > etaMax_)
       continue;
  
     float drmin = 999;
  
     PropState matchProp;
     int match_idx = -1;
     int il1tk = -1;
  
     int nTracksMatch = 0;
  
     for (const auto& l1tk : l1tks) {
       il1tk++;
  
       float l1tk_pt = l1tk.momentum().perp();
       if (l1tk_pt < pTMinTra_)
         continue;
  
       float l1tk_z = l1tk.POCA().z();
       if (std::abs(l1tk_z) > zMax_)
         continue;
  
       float l1tk_chi2 = l1tk.chi2();
       if (l1tk_chi2 > chi2Max_)
         continue;
  
       int l1tk_nstubs = l1tk.getStubRefs().size();
       if (l1tk_nstubs < nStubsmin_)
         continue;
  
       float l1tk_eta = l1tk.momentum().eta();
       float l1tk_phi = l1tk.momentum().phi();
  
       float dr2 = reco::deltaR2(l1mu_eta, l1mu_phi, l1tk_eta, l1tk_phi);
       if (dr2 > dr2_cutoff)
         continue;
  
       nTracksMatch++;
  
       const PropState& pstate = propagateToGMT(l1tk);
       if (!pstate.valid)
         continue;
  
       float dr2prop = reco::deltaR2(l1mu_eta, l1mu_phi, pstate.eta, pstate.phi);
       // FIXME: check if this matching procedure can be improved with
       // a pT dependent dR window
       if (dr2prop < drmin) {
         drmin = dr2prop;
         match_idx = il1tk;
         matchProp = pstate;
       }
     }  // over l1tks
  
     LogDebug("L1TkMuonProducer") << "matching index is " << match_idx;
     if (match_idx >= 0) {
       const L1TTTrackType& matchTk = l1tks[match_idx];
  
       float sigmaEta = sigmaEtaTP(*l1mu);
       float sigmaPhi = sigmaPhiTP(*l1mu);
  
       float etaCut = matching_factor_eta * sqrt(sigmaEta * sigmaEta + matchProp.sigmaEta * matchProp.sigmaEta);
       float phiCut = matching_factor_phi * sqrt(sigmaPhi * sigmaPhi + matchProp.sigmaPhi * matchProp.sigmaPhi);
  
       float dEta = std::abs(matchProp.eta - l1mu_eta);
       float dPhi = std::abs(deltaPhi(matchProp.phi, l1mu_phi));
  
       bool matchCondition = useTPMatchWindows_ ? dEta < etaCut && dPhi < phiCut : drmin < dRMax_;
  
       if (matchCondition) {
         edm::Ptr<L1TTTrackType> l1tkPtr(l1tksH, match_idx);
  
         const auto& p3 = matchTk.momentum();
         float p4e = sqrt(mu_mass * mu_mass + p3.mag2());
  
         math::XYZTLorentzVector l1tkp4(p3.x(), p3.y(), p3.z(), p4e);
  
         const auto& tkv3 = matchTk.POCA();
         math::XYZPoint v3(tkv3.x(), tkv3.y(), tkv3.z());  // why is this defined?
  
         float trkisol = -999;
  
         TkMuon l1tkmu(l1tkp4, l1muRef, l1tkPtr, trkisol);
  
         l1tkmu.setTrackCurvature(matchTk.rInv());
         l1tkmu.setTrkzVtx((float)tkv3.z());
         l1tkmu.setdR(drmin);
         l1tkmu.setNTracksMatched(nTracksMatch);
         l1tkmu.setMuonDetector(detector);
         tkMuons.push_back(l1tkmu);
       }
     }
   }  //over l1mus
 }

References funct::abs(), BXVector< T >::begin(), chi2Max_, SiPixelRawToDigiRegional_cfi::deltaPhi, reco::deltaR2(), HLT_FULL_cff::dEta, hgcalTestNeighbor_cfi::detector, HLT_FULL_cff::dPhi, dr2_cutoff, dRMax_, BXVector< T >::end(), L1TkMuonProducer::PropState::eta, eta_scale, fftjetproducer_cfi::etaCut, etaMax_, etaMin_, LogDebug, matching_factor_eta, matching_factor_phi, TTTrack< T >::momentum(), mu_mass, nStubsmin_, p3, L1TkMuonProducer::PropState::phi, phi_scale, TTTrack< T >::POCA(), edm::Handle< T >::product(), propagateToGMT(), pTMinTra_, TTTrack< T >::rInv(), l1t::TkMuon::setdR(), l1t::TkMuon::setMuonDetector(), l1t::TkMuon::setNTracksMatched(), l1t::TkMuon::setTrackCurvature(), l1t::TkMuon::setTrkzVtx(), L1TkMuonProducer::PropState::sigmaEta, sigmaEtaTP(), L1TkMuonProducer::PropState::sigmaPhi, sigmaPhiTP(), mathSSE::sqrt(), l1pfProducer_cfi::tkMuons, useTPMatchWindows_, L1TkMuonProducer::PropState::valid, and zMax_.

Referenced by produce().

◆ runOnMTFCollection_v2()

void L1TkMuonProducer::runOnMTFCollection_v2	(	const edm::Handle< EMTFTrackCollection > &	muonH,
		const edm::Handle< L1TTTrackCollectionType > &	l1tksH,
		TkMuonCollection &	tkMuons
	)		const

private

Definition at line 482 of file L1TkMuonProducer.cc.

                                                                               {
   const EMTFTrackCollection& l1mus = (*muonH.product());
   const L1TTTrackCollectionType& l1trks = (*l1tksH.product());
   const auto& corr_mu_idxs = dwcorr_->find_match(l1mus, l1trks);
   // it's a vector with as many entries as the L1TT vector.
   // >= 0 : the idx in the EMTF vector of matched mu
   // < 0: no match
  
   // sanity check
   if (corr_mu_idxs.size() != l1trks.size())
     throw cms::Exception("TkMuAlgoOutput")
         << "the size of tkmu indices does not match the size of input trk collection\n";
  
   for (uint il1ttrack = 0; il1ttrack < corr_mu_idxs.size(); ++il1ttrack) {
     int emtf_idx = corr_mu_idxs[il1ttrack];
     if (emtf_idx < 0)
       continue;
  
     const L1TTTrackType& matchTk = l1trks[il1ttrack];
     const auto& p3 = matchTk.momentum();
     const auto& tkv3 = matchTk.POCA();
     float p4e = sqrt(mu_mass * mu_mass + p3.mag2());
     math::XYZTLorentzVector l1tkp4(p3.x(), p3.y(), p3.z(), p4e);
  
     edm::Ref<RegionalMuonCandBxCollection> l1muRef;
     edm::Ptr<L1TTTrackType> l1tkPtr(l1tksH, il1ttrack);
     float trkisol = -999;  // now doing as in the TP algo
     TkMuon l1tkmu(l1tkp4, l1muRef, l1tkPtr, trkisol);
     l1tkmu.setTrackCurvature(matchTk.rInv());
     l1tkmu.setTrkzVtx((float)tkv3.z());
     l1tkmu.setMuonDetector(endcap_MTF_region);
     tkMuons.push_back(l1tkmu);
   }
  
   return;
 }

References dwcorr_, endcap_MTF_region, TTTrack< T >::momentum(), mu_mass, p3, TTTrack< T >::POCA(), edm::Handle< T >::product(), TTTrack< T >::rInv(), l1t::TkMuon::setMuonDetector(), l1t::TkMuon::setTrackCurvature(), l1t::TkMuon::setTrkzVtx(), mathSSE::sqrt(), l1pfProducer_cfi::tkMuons, and parallelization::uint.

Referenced by produce().

◆ sigmaEtaTP()

double L1TkMuonProducer::sigmaEtaTP ( const RegionalMuonCand & mu ) const

private

Definition at line 585 of file L1TkMuonProducer.cc.

                                                                       {
   float l1mu_eta = l1mu.hwEta() * eta_scale;
   if (std::abs(l1mu_eta) <= 1.55)
     return 0.0288;
   else if (std::abs(l1mu_eta) > 1.55 && std::abs(l1mu_eta) <= 1.65)
     return 0.025;
   else if (std::abs(l1mu_eta) > 1.65 && std::abs(l1mu_eta) <= 2.4)
     return 0.0144;
   return 0.0288;
 }

References funct::abs(), eta_scale, and l1t::RegionalMuonCand::hwEta().

Referenced by runOnMTFCollection_v1().

◆ sigmaPhiTP()

double L1TkMuonProducer::sigmaPhiTP ( const RegionalMuonCand & mu ) const

private

Definition at line 596 of file L1TkMuonProducer.cc.

596 { return 0.0126; }

Referenced by runOnMTFCollection_v1().

Member Data Documentation

◆ bmtfMatchAlgoVersion_

AlgoType L1TkMuonProducer::bmtfMatchAlgoVersion_

private

Definition at line 110 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer(), and produce().

◆ bmtfToken_

const edm::EDGetTokenT<RegionalMuonCandBxCollection> L1TkMuonProducer::bmtfToken_

private

Definition at line 121 of file L1TkMuonProducer.cc.

Referenced by produce().

◆ chi2Max_

float L1TkMuonProducer::chi2Max_

private

Definition at line 102 of file L1TkMuonProducer.cc.

Referenced by runOnMTFCollection_v1().

◆ correctGMTPropForTkZ_

bool L1TkMuonProducer::correctGMTPropForTkZ_

private

Definition at line 106 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer(), and propagateToGMT().

◆ dRMax_

float L1TkMuonProducer::dRMax_

private

Definition at line 104 of file L1TkMuonProducer.cc.

Referenced by runOnMTFCollection_v1().

◆ dwcorr_

std::unique_ptr<L1TkMuCorrDynamicWindows> L1TkMuonProducer::dwcorr_

private

Definition at line 114 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer(), and runOnMTFCollection_v2().

◆ emtfMatchAlgoVersion_

AlgoType L1TkMuonProducer::emtfMatchAlgoVersion_

private

Definition at line 112 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer(), and produce().

◆ emtfTCToken_

const edm::EDGetTokenT<EMTFTrackCollection> L1TkMuonProducer::emtfTCToken_

private

Definition at line 125 of file L1TkMuonProducer.cc.

Referenced by produce().

◆ emtfToken_

const edm::EDGetTokenT<RegionalMuonCandBxCollection> L1TkMuonProducer::emtfToken_

private

Definition at line 123 of file L1TkMuonProducer.cc.

Referenced by produce().

◆ etaMax_

float L1TkMuonProducer::etaMax_

private

Definition at line 100 of file L1TkMuonProducer.cc.

Referenced by runOnMTFCollection_v1().

◆ etaMin_

float L1TkMuonProducer::etaMin_

private

Definition at line 99 of file L1TkMuonProducer.cc.

Referenced by runOnMTFCollection_v1().

◆ mantra_n_trk_par_

int L1TkMuonProducer::mantra_n_trk_par_

private

Definition at line 119 of file L1TkMuonProducer.cc.

◆ mantracorr_barr_

std::unique_ptr<L1TkMuMantra> L1TkMuonProducer::mantracorr_barr_

private

Definition at line 116 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer(), and produce().

◆ mantracorr_endc_

std::unique_ptr<L1TkMuMantra> L1TkMuonProducer::mantracorr_endc_

private

Definition at line 118 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer(), and produce().

◆ mantracorr_ovrl_

std::unique_ptr<L1TkMuMantra> L1TkMuonProducer::mantracorr_ovrl_

private

Definition at line 117 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer(), and produce().

◆ nStubsmin_

int L1TkMuonProducer::nStubsmin_

private

Definition at line 105 of file L1TkMuonProducer.cc.

Referenced by runOnMTFCollection_v1().

◆ omtfMatchAlgoVersion_

AlgoType L1TkMuonProducer::omtfMatchAlgoVersion_

private

Definition at line 111 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer(), and produce().

◆ omtfToken_

const edm::EDGetTokenT<RegionalMuonCandBxCollection> L1TkMuonProducer::omtfToken_

private

Definition at line 122 of file L1TkMuonProducer.cc.

Referenced by produce().

◆ pTMinTra_

float L1TkMuonProducer::pTMinTra_

private

Definition at line 103 of file L1TkMuonProducer.cc.

Referenced by runOnMTFCollection_v1().

◆ trackToken_

const edm::EDGetTokenT<std::vector<TTTrack<Ref_Phase2TrackerDigi_> > > L1TkMuonProducer::trackToken_

private

Definition at line 126 of file L1TkMuonProducer.cc.

Referenced by produce().

◆ use5ParameterFit_

bool L1TkMuonProducer::use5ParameterFit_

private

Definition at line 107 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer().

◆ useTPMatchWindows_

bool L1TkMuonProducer::useTPMatchWindows_

private

Definition at line 108 of file L1TkMuonProducer.cc.

Referenced by L1TkMuonProducer(), and runOnMTFCollection_v1().

◆ zMax_

float L1TkMuonProducer::zMax_

private

Definition at line 101 of file L1TkMuonProducer.cc.

Referenced by runOnMTFCollection_v1().

Classes

Public Types

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ L1TTTrackCollectionType

◆ L1TTTrackType

Member Enumeration Documentation

◆ AlgoType

Constructor & Destructor Documentation

◆ L1TkMuonProducer()

◆ ~L1TkMuonProducer()

Member Function Documentation

◆ build_tkMuons_from_idxs()

◆ fillDescriptions()

◆ produce()

◆ product_to_muvec() [1/2]

◆ product_to_muvec() [2/2]

◆ product_to_trkvec()

◆ propagateToGMT()

◆ runOnMTFCollection_v1()

◆ runOnMTFCollection_v2()

◆ sigmaEtaTP()

◆ sigmaPhiTP()

Member Data Documentation

◆ bmtfMatchAlgoVersion_

◆ bmtfToken_

◆ chi2Max_

◆ correctGMTPropForTkZ_

◆ dRMax_

◆ dwcorr_

◆ emtfMatchAlgoVersion_

◆ emtfTCToken_

◆ emtfToken_

◆ etaMax_

◆ etaMin_

◆ mantra_n_trk_par_

◆ mantracorr_barr_

◆ mantracorr_endc_

◆ mantracorr_ovrl_

◆ nStubsmin_

◆ omtfMatchAlgoVersion_

◆ omtfToken_

◆ pTMinTra_

◆ trackToken_

◆ use5ParameterFit_

◆ useTPMatchWindows_

◆ zMax_