Phase2L1GMT::Isolation Class Reference

#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.

: TopoAlgo(cpy) {}

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, 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;
      std::vector<unsigned> overlaps;
      for (auto t : convertedTracks) {
        unsigned ovrl = compute_trk_iso(mu, t);
        if (ovrl != 0) {
          accum += OverlapRemoval(ovrl, overlaps) * t.pt();
        }
      }

      // Accumation without fixed bit width
      mu.setHwIsoSum(accum);

      // Accumation with fixed bit width
      // Bit shifts with 3 bits with LSB of 0.25GeV
      iso_accum_t temp(accum);
      mu.setHwIsoSumAp(temp.to_int() >> 3);

      //Disable isolation bit, sending isolation sumPT to GT
      //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_.

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

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