PFBlockElementTrack.cc

Go to the documentation of this file.

#include "DataFormats/ParticleFlowReco/interface/PFBlockElementTrack.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;


PFBlockElementTrack::PFBlockElementTrack(const PFRecTrackRef& ref) : 
  PFBlockElement( TRACK ),
  trackRefPF_( ref ), 
  trackRef_( ref->trackRef() ),
  trackType_(0) {
  

  if( ref.isNull() ) 
    throw cms::Exception("NullRef")<<" PFBlockElementTrack constructed from a null reference to PFRecTrack.";
  
  const reco::PFTrajectoryPoint& atECAL 
    = ref->extrapolatedPoint( reco::PFTrajectoryPoint::ECALEntrance );

  if( atECAL.isValid() ) 
    positionAtECALEntrance_.SetCoordinates( atECAL.position().x(),
                        atECAL.position().y(),
                        atECAL.position().z() );
  // if the position at ecal entrance is invalid, 
  // positionAtECALEntrance_ is initialized by default to 0,0,0

  setTrackType( DEFAULT, true );      
} 


void PFBlockElementTrack::Dump(ostream& out, 
                               const char* tab ) const {
  
  if(! out ) return;
  
  if( !trackRef_.isNull() ) {
    
    double charge = trackRef_->charge();
    double pt = trackRef_->pt();
    double p = trackRef_->p();
    string s = "  at vertex";
    double tracketa = trackRef_->eta();
    double trackphi = trackRef_->phi();

    // COLIN
    // the following lines rely on the presence of the PFRecTrack, 
    // which for most people is not there (PFRecTracks are transient) 
    // commented these lines out to remove the spurious error message
    // for the missing PFRecTrack product
    //     const reco::PFTrajectoryPoint& atECAL 
    //       = trackRefPF_->extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax );
    //     // check if  reach ecal Shower max 
    //     if( atECAL.isValid() ) { 
    //       s = "  at ECAL shower max";  
    //       tracketa = atECAL.position().Eta();
    //       trackphi = atECAL.position().Phi();
    //     }
    
    out<<setprecision(0);
    out<<tab<<setw(7)<<"charge="<<setw(3)<<charge;
    out<<setprecision(3);
    out<<setiosflags(ios::right);
    out<<setiosflags(ios::fixed);
    out<<", pT ="<<setw(7)<<pt;
    out<<", p ="<<setw(7)<<p;
    out<<" (eta,phi)= (";
    out<<tracketa<<",";
    out<<trackphi<<")" << s;
    
    out<<resetiosflags(ios::right|ios::fixed);  
  }
}