IsolatedPixelTrackCandidate.cc

Go to the documentation of this file.

#include "DataFormats/HcalIsolatedTrack/interface/IsolatedPixelTrackCandidate.h"

using namespace reco;

IsolatedPixelTrackCandidate::IsolatedPixelTrackCandidate(const IsolatedPixelTrackCandidate& right) : RecoCandidate(right), track_(right.track_), l1tauJet_(right.l1tauJet_) {
  maxPtPxl_  = right.maxPtPxl_;
  sumPtPxl_  = right.sumPtPxl_;
  enIn_      = right.enIn_;
  enOut_     = right.enOut_;
  nhitIn_    = right.nhitIn_;
  nhitOut_   = right.nhitOut_;
  etaPhiEcal_= right.etaPhiEcal_;
  etaEcal_   = right.etaEcal_; 
  phiEcal_   = right.phiEcal_;  
}

IsolatedPixelTrackCandidate::~IsolatedPixelTrackCandidate() { }


IsolatedPixelTrackCandidate * IsolatedPixelTrackCandidate::clone() const { 
  return new IsolatedPixelTrackCandidate( * this ); 
}

TrackRef IsolatedPixelTrackCandidate::track() const {
  return track_;
}

l1extra::L1JetParticleRef IsolatedPixelTrackCandidate::l1tau() const {
  return l1tauJet_;
}

bool IsolatedPixelTrackCandidate::overlap( const Candidate & c ) const {
  const RecoCandidate * o = dynamic_cast<const RecoCandidate *>( & c );
  return ( o != 0 &&  checkOverlap( track(), o->track() ) );
}

std::pair<int,int> IsolatedPixelTrackCandidate::towerIndex() const
{
  int ieta=0, iphi=0;
  for (int i=1; i<21; i++)
    {
      if (fabs(track_->eta())<(i*0.087)&&fabs(track_->eta())>(i-1)*0.087) ieta=int(fabs(track_->eta())/track_->eta())*i;
    }
  if (fabs(track_->eta())>1.740&&fabs(track_->eta())<1.830) ieta=int(fabs(track_->eta())/track_->eta())*21;
  if (fabs(track_->eta())>1.830&&fabs(track_->eta())<1.930) ieta=int(fabs(track_->eta())/track_->eta())*22;
  if (fabs(track_->eta())>1.930&&fabs(track_->eta())<2.043) ieta=int(fabs(track_->eta())/track_->eta())*23;

  double delta=track_->phi()+0.174532925;
  if (delta<0) delta=delta+2*acos(-1);
  if (fabs(track_->eta())<1.740) 
    {
      for (int i=0; i<72; i++)
    {
      if (delta<(i+1)*0.087266462&&delta>i*0.087266462) iphi=i;
    }
    }
  else 
    {
      for (int i=0; i<36; i++)
    {
      if (delta<2*(i+1)*0.087266462&&delta>2*i*0.087266462) iphi=2*i;
    }
    }

  return std::pair<int,int>(ieta,iphi);

}