FindEtaPhi.cc

Go to the documentation of this file.

#include "Calibration/IsolatedParticles/interface/FindEtaPhi.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

namespace spr {

  spr::EtaPhi getEtaPhi(int ieta, int iphi, bool debug) {
    int ietal = (ieta - 1) / 2;
    int ietar = ieta - ietal - 1;
    int iphil = (iphi - 1) / 2;
    int iphir = iphi - iphil - 1;
    spr::EtaPhi etaphi;
    etaphi.ietaE[0] = ietal;
    etaphi.ietaW[0] = ietar;
    etaphi.iphiN[0] = iphil;
    etaphi.iphiS[0] = iphir;
    if (ietal == ietar && iphil == iphir) {
      etaphi.ntrys = 1;
    } else if (ietal == ietar || iphil == iphir) {
      etaphi.ntrys = 2;
      etaphi.ietaE[1] = ietar;
      etaphi.ietaW[1] = ietal;
      etaphi.iphiN[1] = iphir;
      etaphi.iphiS[1] = iphil;
    } else {
      etaphi.ntrys = 4;
      etaphi.ietaE[1] = ietar;
      etaphi.ietaW[1] = ietal;
      etaphi.iphiN[1] = iphil;
      etaphi.iphiS[1] = iphir;
      etaphi.ietaE[2] = ietal;
      etaphi.ietaW[1] = ietar;
      etaphi.iphiN[2] = iphir;
      etaphi.iphiS[1] = iphil;
      etaphi.ietaE[3] = ietar;
      etaphi.ietaW[1] = ietal;
      etaphi.iphiN[3] = iphir;
      etaphi.iphiS[1] = iphil;
    }

    if (debug) {
      edm::LogVerbatim("IsoTrack") << "getEtaPhi:: Has " << etaphi.ntrys << " possibilites for " << ieta << "X" << iphi
                                   << " matrix";
      for (int itry = 0; itry < etaphi.ntrys; itry++) {
        edm::LogVerbatim("IsoTrack") << "Trial " << itry << " with etaE|etaW " << etaphi.ietaE[itry] << "|"
                                     << etaphi.ietaW[itry] << " and phiN|PhiS " << etaphi.iphiN[itry] << "|"
                                     << etaphi.iphiS[itry];
      }
    }
    return etaphi;
  }
}  // namespace spr