L1Trigger/TrackFindingTracklet/src/Track.cc - master - githubmirror/cmssw

L1Trigger/TrackFindingTracklet/src/Track.cc

54 строки1 KB

Manos Vourliotis

Rename reduceRange to reducePhiRange

21 янв 2025, 01:04

21 янв 2025, 01:043139210

#include "L1Trigger/TrackFindingTracklet/interface/Track.h"

#include "DataFormats/Math/interface/deltaPhi.h"

#include <algorithm>

using namespace std;
using namespace trklet;

Track::Track(TrackPars<int> ipars,
             int ichisqrphi,
             int ichisqrz,
             double chisqrphi,
             double chisqrz,
             int hitpattern,
             std::map<int, int> stubID,
             const std::vector<L1TStub>& l1stub,
             int seed) {
  ipars_ = ipars;
  ichisqrphi_ = ichisqrphi;
  ichisqrz_ = ichisqrz;

  chisqrphi_ = chisqrphi;
  chisqrz_ = chisqrz;

  hitpattern_ = hitpattern;

  nstubs_ = std::max((int)l1stub.size(), (int)N_FITSTUB);

  stubID_ = stubID;
  l1stub_ = l1stub;

  seed_ = seed;
  duplicate_ = false;
  sector_ = -1;
}

double Track::phi0(Settings const& settings) const {
  double dphi = 2 * M_PI / N_SECTOR;
  double dphiHG = 0.5 * settings.dphisectorHG() - M_PI / N_SECTOR;
  double phimin = sector_ * dphi - dphiHG;
  double phimax = phimin + dphi + 2 * dphiHG;
  phimin -= M_PI / N_SECTOR;
  phimax -= M_PI / N_SECTOR;
  phimin = reco::reducePhiRange(phimin);
  phimax = reco::reducePhiRange(phimax);
  if (phimin > phimax)
    phimin -= 2 * M_PI;
  double phioffset = phimin;

  double phi0 = ipars_.phi0() * settings.kphi0pars() + phioffset;
  phi0 = reco::reducePhiRange(phi0);
  return phi0;
}