#include <Isolation.h>

Inheritance diagram for Phase2L1GMT::Isolation:

Public Member Functions
unsigned	compute_trk_iso (l1t::TrackerMuon &in_mu, ConvertedTTTrack &in_trk)

	Isolation (const edm::ParameterSet &iConfig)

	Isolation (const Isolation &cpy)

void	isolation_allmu_alltrk (std::vector< l1t::TrackerMuon > &trkMus, std::vector< ConvertedTTTrack > &convertedTracks)

	~Isolation ()

Public Member Functions inherited from Phase2L1GMT::TopoAlgo
int	deltaEta (const int eta1, const int eta2)

int	deltaPhi (int phi1, int phi2)

int	deltaZ0 (const int Z01, const int Z02)

void	DumpInputs ()

void	load (std::vector< l1t::TrackerMuon > &trkMus, std::vector< ConvertedTTTrack > &convertedTracks)

	TopoAlgo ()

	TopoAlgo (const TopoAlgo &cpy)

	~TopoAlgo ()

Private Types
typedef ap_ufixed< 9, 9, AP_TRN, AP_SAT >	iso_accum_t

typedef ap_ufixed< 9, 0 >	reliso_thresh_t

Private Member Functions
void	DumpOutputs (std::vector< l1t::TrackerMuon > &trkMus)

int	OverlapRemoval (unsigned &ovrl, std::vector< unsigned > &overlaps)

int	SetAbsIsolationBits (int accum)

int	SetRelIsolationBits (int accum, int mupt)

Private Attributes
int	absiso_thrL

int	absiso_thrM

int	absiso_thrT

bool	dumpForHLS_

std::ofstream	dumpOutput

double	reliso_thrL

double	reliso_thrM

double	reliso_thrT

bool	verbose_

Static Private Attributes
static const int	c_iso_dangle_max = 260

static const int	c_iso_dz_max = 17

static const int	c_iso_pt_min = 120

Additional Inherited Members
Protected Attributes inherited from Phase2L1GMT::TopoAlgo
std::vector< ConvertedTTTrack > *	convertedTracks

std::ofstream	dumpInput

std::vector< l1t::TrackerMuon > *	trkMus

Detailed Description

Definition at line 25 of file Isolation.h.

Member Typedef Documentation

◆ iso_accum_t

typedef ap_ufixed<9, 9, AP_TRN, AP_SAT> Phase2L1GMT::Isolation::iso_accum_t

private

Definition at line 56 of file Isolation.h.

◆ reliso_thresh_t

typedef ap_ufixed<9, 0> Phase2L1GMT::Isolation::reliso_thresh_t

private

Definition at line 57 of file Isolation.h.

Constructor & Destructor Documentation

◆ Isolation() [1/2]

Phase2L1GMT::Isolation::Isolation ( const edm::ParameterSet & iConfig )

inline

Definition at line 60 of file Isolation.h.

References dumpForHLS_, Phase2L1GMT::TopoAlgo::dumpInput, dumpOutput, and MillePedeFileConverter_cfg::out.

       : absiso_thrL(iConfig.getParameter<int>("AbsIsoThresholdL")),
         absiso_thrM(iConfig.getParameter<int>("AbsIsoThresholdM")),
         absiso_thrT(iConfig.getParameter<int>("AbsIsoThresholdT")),
         reliso_thrL(iConfig.getParameter<double>("RelIsoThresholdL")),
         reliso_thrM(iConfig.getParameter<double>("RelIsoThresholdM")),
         reliso_thrT(iConfig.getParameter<double>("RelIsoThresholdT")),
         verbose_(iConfig.getParameter<int>("verbose")),
         dumpForHLS_(iConfig.getParameter<int>("IsodumpForHLS")) {
     if (dumpForHLS_) {
       dumpInput.open("Isolation_Mu_Track_infolist.txt", std::ofstream::out);
       dumpOutput.open("Isolation_Mu_Isolation.txt", std::ofstream::out);
     }
   }

◆ ~Isolation()

Phase2L1GMT::Isolation::~Isolation ( )

inline

Definition at line 75 of file Isolation.h.

References dumpForHLS_, Phase2L1GMT::TopoAlgo::dumpInput, and dumpOutput.

                                {
     if (dumpForHLS_) {
       dumpInput.close();
       dumpOutput.close();
     }
   }

◆ Isolation() [2/2]

Phase2L1GMT::Isolation::Isolation ( const Isolation & cpy )

inline

Definition at line 82 of file Isolation.h.

82 : TopoAlgo(cpy) {}

Phase2L1GMT::TopoAlgo::TopoAlgo

TopoAlgo()

Definition: TopologicalAlgorithm.h:46

Member Function Documentation

◆ compute_trk_iso()

unsigned Phase2L1GMT::Isolation::compute_trk_iso	(	l1t::TrackerMuon &	in_mu,
		ConvertedTTTrack &	in_trk
	)

inline

Definition at line 198 of file Isolation.h.

References c_iso_dangle_max, c_iso_dz_max, c_iso_pt_min, Phase2L1GMT::TopoAlgo::deltaEta(), Phase2L1GMT::TopoAlgo::deltaPhi(), Phase2L1GMT::TopoAlgo::deltaZ0(), Phase2L1GMT::ConvertedTTTrack::eta(), l1t::L1Candidate::hwEta(), l1t::L1Candidate::hwPhi(), l1t::L1Candidate::hwPt(), l1t::TrackerMuon::hwZ0(), Phase2L1GMT::ConvertedTTTrack::phi(), Phase2L1GMT::ConvertedTTTrack::pt(), funct::true, verbose_, and Phase2L1GMT::ConvertedTTTrack::z0().

Referenced by isolation_allmu_alltrk().

                                                                                             {
     int dphi = deltaPhi(in_mu.hwPhi(), in_trk.phi());
     int deta = deltaEta(in_mu.hwEta(), in_trk.eta());
     int dz0 = deltaZ0(in_mu.hwZ0(), in_trk.z0());
 
     bool pass_deta = (deta < c_iso_dangle_max ? true : false);
     bool pass_dphi = (dphi < c_iso_dangle_max ? true : false);
     bool pass_dz0 = (dz0 < c_iso_dz_max ? true : false);
     bool pass_trkpt = (in_trk.pt() >= c_iso_pt_min ? true : false);
     bool pass_ovrl = (deta > 0 || dphi > 0 ? true : false);
 
     if (verbose_) {
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : Start of debug msg for compute_trk_iso";
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : incoming muon (pt / eta / phi / z0 / isvalid)";
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : MU  =  " << in_mu.hwPt() << " / " << in_mu.hwEta()
                                 << " / " << in_mu.hwPhi() << " / " << in_mu.hwZ0() << " / " << 1;
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : incoming track (pt / eta / phi / z0 / isvalid)";
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : TRK =  " << in_trk.pt() << " / " << in_trk.eta() << " / "
                                 << in_trk.phi() << " / " << in_trk.z0() << " / " << 1;
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : Delta phi : " << dphi;
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : Delta eta : " << deta;
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : Delta z0  : " << dz0;
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : pass_deta      : " << pass_deta;
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : pass_dphi      : " << pass_dphi;
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : pass_dz0       : " << pass_dz0;
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : pass_trkpt     : " << pass_trkpt;
       edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : pass_ovrl      : " << pass_ovrl;
     }
     // match conditions
     if (pass_deta && pass_dphi && pass_dz0 && pass_trkpt && pass_ovrl) {
       if (verbose_) {
         edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : THE TRACK WAS MATCHED";
         edm::LogInfo("Isolation") << " [DEBUG compute_trk_iso] : RETURN         : " << in_trk.pt();
       }
 
       //return in_trk.pt();
       // Return fixed bit output for duplication removal.
       // dZ0(8bis) + deta(10bits)+dphi(10bits)
       unsigned int retbits = 0;
       retbits |= (dz0 & ((1 << 9) - 1)) << 20;
       retbits |= (deta & ((1 << 11) - 1)) << 10;
       retbits |= (dphi & ((1 << 11) - 1));
       return retbits;
     } else {
       return 0;
     }
   }

◆ DumpOutputs()

void Phase2L1GMT::Isolation::DumpOutputs ( std::vector< l1t::TrackerMuon > & trkMus )

inlineprivate

Definition at line 84 of file Isolation.h.

References dumpOutput, mps_fire::i, Phase2L1GMT::LSBeta, Phase2L1GMT::LSBGTz0, Phase2L1GMT::LSBphi, Phase2L1GMT::LSBpt, M_PI, amptDefaultParameters_cff::mu, and Phase2L1GMT::TopoAlgo::trkMus.

Referenced by isolation_allmu_alltrk().

                                                                       {
     static std::atomic<int> nevto = 0;
     auto evto = nevto++;
     for (unsigned int i = 0; i < trkMus.size(); ++i) {
       auto mu = trkMus.at(i);
       if (mu.hwPt() != 0) {
         double convertphi = mu.hwPhi() * LSBphi;
         if (convertphi > M_PI) {
           convertphi -= 2 * M_PI;
         }
         dumpOutput << evto << " " << i << " " << mu.hwPt() * LSBpt << " " << mu.hwEta() * LSBeta << " " << convertphi
                    << " " << mu.hwZ0() * LSBGTz0 << " " << mu.hwIso() << endl;
       }
     }
   }

◆ isolation_allmu_alltrk()

void Phase2L1GMT::Isolation::isolation_allmu_alltrk	(	std::vector< l1t::TrackerMuon > &	trkMus,
		std::vector< ConvertedTTTrack > &	convertedTracks
	)

inline

Definition at line 139 of file Isolation.h.

References compute_trk_iso(), Phase2L1GMT::TopoAlgo::convertedTracks, dumpForHLS_, Phase2L1GMT::TopoAlgo::DumpInputs(), DumpOutputs(), Phase2L1GMT::TopoAlgo::load(), amptDefaultParameters_cff::mu, OverlapRemoval(), analyzePatCleaning_cfg::overlaps, SetAbsIsolationBits(), SetRelIsolationBits(), submitPVValidationJobs::t, groupFilesInBlocks::temp, Phase2L1GMT::TopoAlgo::trkMus, and verbose_.

                                                                                               {
     load(trkMus, convertedTracks);
     if (dumpForHLS_) {
       DumpInputs();
     }
 
     static std::atomic<int> itest = 0;
     if (verbose_) {
       edm::LogInfo("Isolation") << "........ RUNNING TEST NUMBER .......... " << itest++;
     }
 
     for (auto &mu : trkMus) {
       int accum = 0;
       int iso_ = 0;
       std::vector<unsigned> overlaps;
       for (auto t : convertedTracks) {
         unsigned ovrl = compute_trk_iso(mu, t);
         if (ovrl != 0) {
           accum += OverlapRemoval(ovrl, overlaps) * t.pt();
         }
       }
 
       // Only 8 bit for accumation?
       mu.setHwIsoSum(accum);
 
       iso_accum_t temp(accum);
       accum = temp.to_int();
 
       mu.setHwIsoSumAp(accum);
 
       iso_ |= SetAbsIsolationBits(accum);
       iso_ |= SetRelIsolationBits(accum, mu.hwPt());
 
       mu.setHwIso(iso_);
     }
 
     if (dumpForHLS_) {
       DumpOutputs(trkMus);
     }
   }

◆ OverlapRemoval()

int Phase2L1GMT::Isolation::OverlapRemoval	(	unsigned &	ovrl,
		std::vector< unsigned > &	overlaps
	)

inlineprivate

Definition at line 185 of file Isolation.h.

References change_name::diff, mps_fire::i, and analyzePatCleaning_cfg::overlaps.

Referenced by isolation_allmu_alltrk().

                                                                                     {
     for (auto i : overlaps) {
       // same tracks with Phi can be off by 1 LSB
       unsigned diff = ovrl - i;
       if (diff <= 1 || diff == unsigned(-1)) {
         // When Overlap, return 0 so that this track won't be consider
         return 0;
       }
     }
     overlaps.push_back(ovrl);
     return 1;
   }  // -----  end of function Isolation::OverlapRemoval  -----

◆ SetAbsIsolationBits()

int Phase2L1GMT::Isolation::SetAbsIsolationBits ( int accum )

inlineprivate

Definition at line 100 of file Isolation.h.

References absiso_thrL, absiso_thrM, absiso_thrT, genparticles_cff::iso, and verbose_.

Referenced by isolation_allmu_alltrk().

                                                      {
     int iso = (accum <= absiso_thrT ? 3 : accum <= absiso_thrM ? 2 : accum <= absiso_thrL ? 1 : 0);
 
     if (verbose_) {
       edm::LogInfo("Isolation") << " [DEBUG Isolation] : absiso_threshold L : " << absiso_thrL << " accum " << accum
                                 << " bit set : " << (accum < absiso_thrL);
       edm::LogInfo("Isolation") << " [DEBUG Isolation] : absiso_threshold M : " << absiso_thrM << " accum " << accum
                                 << " bit set : " << (accum < absiso_thrM);
       edm::LogInfo("Isolation") << " [DEBUG Isolation] : absiso_threshold T : " << absiso_thrT << " accum " << accum
                                 << " bit set : " << (accum < absiso_thrT);
       edm::LogInfo("Isolation") << " [DEBUG Isolation] : absiso : " << (iso);
     }
     return iso;
   }

◆ SetRelIsolationBits()

int Phase2L1GMT::Isolation::SetRelIsolationBits	(	int	accum,
		int	mupt
	)

inlineprivate

Definition at line 115 of file Isolation.h.

References genparticles_cff::iso, reliso_thrL, reliso_thrM, reliso_thrT, and verbose_.

Referenced by isolation_allmu_alltrk().

                                                                {
     const static reliso_thresh_t relisoL(reliso_thrL);
     const static reliso_thresh_t relisoM(reliso_thrM);
     const static reliso_thresh_t relisoT(reliso_thrT);
 
     iso_accum_t thrL = relisoL * mupt;
     iso_accum_t thrM = relisoM * mupt;
     iso_accum_t thrT = relisoT * mupt;
 
     int iso = (accum <= thrT.to_int() ? 3 : accum <= thrM.to_int() ? 2 : accum <= thrL.to_int() ? 1 : 0);
 
     if (verbose_) {
       edm::LogInfo("Isolation") << " [DEBUG Isolation] : reliso_threshold L : " << thrL << " accum " << accum
                                 << " bit set : " << (accum < thrL.to_int());
       edm::LogInfo("Isolation") << " [DEBUG Isolation] : reliso_threshold M : " << thrM << " accum " << accum
                                 << " bit set : " << (accum < thrM.to_int());
       edm::LogInfo("Isolation") << " [DEBUG Isolation] : reliso_threshold T : " << thrT << " accum " << accum
                                 << " bit set : " << (accum < thrT.to_int());
       edm::LogInfo("Isolation") << " [DEBUG Isolation] : reliso : " << (iso << 2) << " org " << iso;
     }
 
     return iso << 2;
   }

Member Data Documentation

◆ absiso_thrL

int Phase2L1GMT::Isolation::absiso_thrL

private

Definition at line 46 of file Isolation.h.

Referenced by SetAbsIsolationBits().

◆ absiso_thrM

int Phase2L1GMT::Isolation::absiso_thrM

private

Definition at line 47 of file Isolation.h.

Referenced by SetAbsIsolationBits().

◆ absiso_thrT

int Phase2L1GMT::Isolation::absiso_thrT

private

Definition at line 48 of file Isolation.h.

Referenced by SetAbsIsolationBits().

◆ c_iso_dangle_max

const int Phase2L1GMT::Isolation::c_iso_dangle_max = 260

staticprivate

Definition at line 41 of file Isolation.h.

Referenced by compute_trk_iso().

◆ c_iso_dz_max

const int Phase2L1GMT::Isolation::c_iso_dz_max = 17

staticprivate

Definition at line 42 of file Isolation.h.

Referenced by compute_trk_iso().

◆ c_iso_pt_min

const int Phase2L1GMT::Isolation::c_iso_pt_min = 120

staticprivate

Definition at line 43 of file Isolation.h.

Referenced by compute_trk_iso().

◆ dumpForHLS_

bool Phase2L1GMT::Isolation::dumpForHLS_

private

Definition at line 53 of file Isolation.h.

Referenced by Isolation(), isolation_allmu_alltrk(), and ~Isolation().

◆ dumpOutput

std::ofstream Phase2L1GMT::Isolation::dumpOutput

private

Definition at line 54 of file Isolation.h.

Referenced by DumpOutputs(), Isolation(), and ~Isolation().

◆ reliso_thrL

double Phase2L1GMT::Isolation::reliso_thrL

private

Definition at line 49 of file Isolation.h.

Referenced by SetRelIsolationBits().

◆ reliso_thrM

double Phase2L1GMT::Isolation::reliso_thrM

private

Definition at line 50 of file Isolation.h.

Referenced by SetRelIsolationBits().

◆ reliso_thrT

double Phase2L1GMT::Isolation::reliso_thrT

private

Definition at line 51 of file Isolation.h.

Referenced by SetRelIsolationBits().

◆ verbose_

bool Phase2L1GMT::Isolation::verbose_

private

Definition at line 52 of file Isolation.h.

Referenced by compute_trk_iso(), isolation_allmu_alltrk(), SetAbsIsolationBits(), and SetRelIsolationBits().

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ iso_accum_t

◆ reliso_thresh_t

Constructor & Destructor Documentation

◆ Isolation() [1/2]

◆ ~Isolation()

◆ Isolation() [2/2]

Member Function Documentation

◆ compute_trk_iso()

◆ DumpOutputs()

◆ isolation_allmu_alltrk()

◆ OverlapRemoval()

◆ SetAbsIsolationBits()

◆ SetRelIsolationBits()

Member Data Documentation

◆ absiso_thrL

◆ absiso_thrM

◆ absiso_thrT

◆ c_iso_dangle_max

◆ c_iso_dz_max

◆ c_iso_pt_min

◆ dumpForHLS_

◆ dumpOutput

◆ reliso_thrL

◆ reliso_thrM

◆ reliso_thrT

◆ verbose_