PFBlockElementGsfTrack.cc

Go to the documentation of this file.

#include "DataFormats/ParticleFlowReco/interface/PFBlockElementGsfTrack.h"
#include "DataFormats/ParticleFlowReco/interface/PFRecTrack.h"
#include "DataFormats/Common/interface/Ref.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/ParticleFlowReco/interface/PFTrajectoryPoint.h"

#include <iomanip>

using namespace reco;
using namespace std;

PFBlockElementGsfTrack::PFBlockElementGsfTrack(const GsfPFRecTrackRef& gsfref,
                                               const math::XYZTLorentzVector& Pin,
                                               const math::XYZTLorentzVector& Pout)
    : PFBlockElement(GSF),
      GsftrackRefPF_(gsfref),
      GsftrackRef_(gsfref->gsfTrackRef()),
      Pin_(Pin),
      Pout_(Pout),
      trackType_(0) {
  if (gsfref.isNull())
    throw cms::Exception("NullRef") << " PFBlockElementGsfTrack constructed from a null reference to PFGsfRecTrack.";
  const reco::PFTrajectoryPoint& atECAL = gsfref->extrapolatedPoint(reco::PFTrajectoryPoint::ECALEntrance);
  if (atECAL.isValid())
    positionAtECALEntrance_.SetCoordinates(atECAL.position().x(), atECAL.position().y(), atECAL.position().z());

  setTrackType(DEFAULT, true);
}

void PFBlockElementGsfTrack::Dump(ostream& out, const char* tab) const {
  if (!out)
    return;

  if (!GsftrackRefPF_.isNull()) {
    //    double charge = trackPF().charge;
    double charge = GsftrackRefPF_->charge();
    math::XYZTLorentzVector pin = Pin_;
    math::XYZTLorentzVector pout = Pout_;
    double ptin = pin.pt();
    double etain = pin.eta();
    double phiin = pin.phi();
    double ptout = pout.pt();
    double etaout = pout.eta();
    double phiout = pout.phi();
    out << setprecision(0);
    out << tab << setw(7) << "charge=" << setw(3) << charge;
    out << setprecision(3);
    out << setiosflags(ios::right);
    out << setiosflags(ios::fixed);
    out << ", Inner pT  =" << setw(7) << ptin;
    out << " Inner (eta,phi)= (";
    out << etain << ",";
    out << phiin << ")";
    out << ", Outer pT  =" << setw(7) << ptout;
    out << " Outer (eta,phi)= (";
    out << etaout << ",";
    out << phiout << ")";
    out << resetiosflags(ios::right | ios::fixed);
  }
}