BPHPlusMinusVertex.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  \author Paolo Ronchese INFN Padova
 *
 */

//-----------------------
// This Class' Header --
//-----------------------
#include "HeavyFlavorAnalysis/RecoDecay/interface/BPHPlusMinusVertex.h"

//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "TrackingTools/TransientTrack/interface/TransientTrack.h"

//---------------
// C++ Headers --
//---------------
#include <iostream>

using namespace std;

//-------------------
// Initializations --
//-------------------


//----------------
// Constructors --
//----------------
BPHPlusMinusVertex::BPHPlusMinusVertex( const edm::EventSetup* es ):
 BPHDecayVertex( es ),
 oldA( true ),
 inRPhi( nullptr ) {
}

//--------------
// Destructor --
//--------------
BPHPlusMinusVertex::~BPHPlusMinusVertex() {
  delete inRPhi;
}

//--------------
// Operations --
//--------------
const ClosestApproachInRPhi& BPHPlusMinusVertex::cAppInRPhi() const {
  if ( oldA ) computeApp();
  if ( inRPhi == nullptr ) {
    static const ClosestApproachInRPhi ca;
    return ca;
  }
  return *inRPhi;
}


bool BPHPlusMinusVertex::chkSize( const string& msg ) const {
  return chkSize( daughters(), msg );
}


void BPHPlusMinusVertex::setNotUpdated() const {
  BPHDecayVertex::setNotUpdated();
  oldA = true;
  return;
}


void BPHPlusMinusVertex::computeApp() const {
  static const string msg =
  "BPHPlusMinusVertex::computeApp: incomplete, no closest approach available";
  delete inRPhi;
  if ( !chkSize( msg ) ) {
    inRPhi = nullptr;
    return;
  }
  inRPhi = new ClosestApproachInRPhi;
  const vector<reco::TransientTrack>& ttk = transientTracks();
  const reco::TransientTrack& ttp = ttk[0];
  const reco::TransientTrack& ttn = ttk[1];
  inRPhi->calculate( ttp.impactPointTSCP().theState(),
                     ttn.impactPointTSCP().theState() );
  oldA = false;
  return;
}