#include <ConfigurableTertiaryTracksVertexFinder.h>

Inheritance diagram for ConfigurableTertiaryTracksVertexFinder:

Public Member Functions
virtual ConfigurableTertiaryTracksVertexFinder *	clone () const
	ConfigurableTertiaryTracksVertexFinder (const VertexFitter< 5 > vf, const VertexUpdator< 5 > vu, const VertexTrackCompatibilityEstimator< 5 > *ve)
AddTvTrack::TrackInfoVector	getTrackInfo ()
virtual std::vector < TransientVertex >	vertices (const std::vector< reco::TransientTrack > &) const
virtual std::vector < TransientVertex >	vertices (const std::vector< reco::TransientTrack > &, const TransientVertex &pv) const
virtual	~ConfigurableTertiaryTracksVertexFinder ()
Private Member Functions
std::vector< TransientVertex >	reconstruct (const std::vector< reco::TransientTrack > &tracks, const TransientVertex &pv) const
Private Attributes
Flight2DSvFilter *	theFlight2DSvFilter
double	theK0sMassWindow
double	theMaxDist2D
double	theMaxInPvFrac
double	theMaxSigOnDistTrackToB
double	theMaxVtxMass
double	theMinDist2D
double	theMinSign2D
double	theMinTrackPt
int	theMinTracks
ConfigurableTrimmedVertexFinder *	theTKVF
AddTvTrack::TrackInfoVector	theTrackInfoVector
V0SvFilter *	theV0SvFilter
Static Private Attributes
static const bool	debug = false

Detailed Description

Definition at line 21 of file ConfigurableTertiaryTracksVertexFinder.h.

Constructor & Destructor Documentation

ConfigurableTertiaryTracksVertexFinder::ConfigurableTertiaryTracksVertexFinder	(	const VertexFitter< 5 > *	vf,
		const VertexUpdator< 5 > *	vu,
		const VertexTrackCompatibilityEstimator< 5 > *	ve
	)

Definition at line 16 of file ConfigurableTertiaryTracksVertexFinder.cc.

Referenced by clone().

{
  theTKVF = new ConfigurableTrimmedVertexFinder(vf,vu,ve);

  theMinTrackPt = 1.0;
  theMaxVtxMass = 6.5;
  theMaxSigOnDistTrackToB = 3.0; // this is being overwritten to 10 in AddTvTtrack (why?)

  theMaxInPvFrac = 0.65;

  // set up V0SvFilter
  theK0sMassWindow = 0.05; // mass window around K0s
  theV0SvFilter = new V0SvFilter(theK0sMassWindow);

  // set up Flight2DSvFilter
  theMaxDist2D = 2.5;  // max transv. dist to beam line
  theMinDist2D = 0.01; // min transv. dist to beam line
  theMinSign2D = 3.0;  // min transverse distance significance
  theMinTracks = 2;    // min number of tracks
  theFlight2DSvFilter= new Flight2DSvFilter(theMaxDist2D,theMinDist2D,
                                            theMinSign2D,theMinTracks);

  //  thePrimaryVertex = new TransientVertex;
  // FIXME this is incomplete!? -> get real primary vertex!

  //theNewTrackInfoVector = new NewTrackInfoVector;

}

ConfigurableTertiaryTracksVertexFinder::~ConfigurableTertiaryTracksVertexFinder ( ) [virtual]

Definition at line 50 of file ConfigurableTertiaryTracksVertexFinder.cc.

{
  delete theTKVF;
  delete theV0SvFilter;
  delete theFlight2DSvFilter;
}

Member Function Documentation

virtual ConfigurableTertiaryTracksVertexFinder* ConfigurableTertiaryTracksVertexFinder::clone ( void ) const [inline, virtual]

Implements VertexReconstructor.

Definition at line 36 of file ConfigurableTertiaryTracksVertexFinder.h.

References ConfigurableTertiaryTracksVertexFinder().

                                                                 {
    return new ConfigurableTertiaryTracksVertexFinder(*this);
  }

AddTvTrack::TrackInfoVector ConfigurableTertiaryTracksVertexFinder::getTrackInfo ( ) [inline]

Definition at line 53 of file ConfigurableTertiaryTracksVertexFinder.h.

References theTrackInfoVector.

Referenced by TertiaryTracksVertexFinder::getTrackInfo().

{ return theTrackInfoVector; }

std::vector< TransientVertex > ConfigurableTertiaryTracksVertexFinder::reconstruct	(	const std::vector< reco::TransientTrack > &	tracks,
		const TransientVertex &	pv
	)		const `[private]`

Definition at line 87 of file ConfigurableTertiaryTracksVertexFinder.cc.

References gather_cfg::cout, debug, TauTagTools::filteredTracks(), AddTvTrack::getSecondaryVertices(), AddTvTrack::getTrackInfo(), TransientTrackInVertices::isInVertex(), TransientVertex::isValid(), and edm::RefVector< C, T, F >::push_back().

{
  // get  primary vertices;
  std::vector<TransientVertex> primaryVertices;
  if(pv.isValid()) { primaryVertices.push_back(pv); 
    if (debug) cout <<"[TTVF] add PV ...\n";
  }

  VertexMass theVertexMass;

  //filter tracks in pt
  std::vector<reco::TransientTrack> filteredTracks;
  for (std::vector<reco::TransientTrack>::const_iterator it=tracks.begin();
       it!=tracks.end();it++)
    if ((*it).impactPointState().globalMomentum().perp() > theMinTrackPt) 
      filteredTracks.push_back(*it);

  if (debug) cout <<"[TTVF] tracks: " << filteredTracks.size() <<endl;

  // get vertices
  std::vector<TransientVertex> vertices;
  if (filteredTracks.size()>1) vertices = theTKVF->vertices(filteredTracks);

  if (debug) cout <<"[TTVF] found secondary vertices with TKVF: "<<vertices.size()<<endl;

  std::vector<TransientVertex> secondaryVertices;

  for(std::vector<TransientVertex>::const_iterator ivx=vertices.begin();
      ivx!=vertices.end(); ivx++) {
    TransientVertex vtx=*ivx;

    double mass=theVertexMass(vtx);
    if (debug) cout <<"[TTVF] new svx: mass: "<<mass<<endl;

    if ((*theV0SvFilter)(vtx)) {
      if (debug) cout <<"[TTVF] survived V0SvFilter\n";
      if((*theFlight2DSvFilter)(vtx)) {
        if (debug) cout <<"[TTVF] survived 2DSvFilter\n";
        if (mass<theMaxVtxMass) {
          if (!primaryVertices.empty()) {
            PvSvFilter thePvSvFilter(theMaxInPvFrac,primaryVertices[0]);
            if (thePvSvFilter(vtx)) secondaryVertices.push_back(vtx);
            else { if (debug) cout <<"[TTVF] failed PvSvFilter\n";}
          }
          else secondaryVertices.push_back(vtx);
        }
        else {if(debug)cout<<"[TTVF] failed mass cut\n";}
      }
    }

  }

  if (debug) cout<<"[TTVF] remaining svx: "<<secondaryVertices.size()<<endl;

  if (primaryVertices.empty() || secondaryVertices.empty()) 
    return secondaryVertices; // unable to reconstruct b-flight-trajectory

  if (debug) cout<<"[TTVF] still here ...\n";

  // find tracks not used in primaryVertex or in vertices
  vector<TransientTrack> unusedTracks;
  for( vector<TransientTrack>::const_iterator itT = filteredTracks.begin(); 
    itT != filteredTracks.end(); itT++ ) 
    if( (!TransientTrackInVertices::isInVertex((*itT),primaryVertices)) 
     && (!TransientTrackInVertices::isInVertex((*itT),vertices)) )
       unusedTracks.push_back( *itT ); 
  if (debug) cout <<"[TTVF] remaining tracks: "<<unusedTracks.size()<<endl;

  // now add tracks to the SV candidate
  AddTvTrack MyAddTVTrack( &primaryVertices, &secondaryVertices, 
                           theMaxSigOnDistTrackToB);  
  vector<TransientVertex> newVertices =
    MyAddTVTrack.getSecondaryVertices(unusedTracks); 

  // for tdr studies
  theTrackInfoVector = MyAddTVTrack.getTrackInfo();

  //std::vector<pair<reco::TransientTrack,double> > theTrackInfo;
  //std::vector<pair<reco::TransientTrack,double* > > theTrackInfo2;
  //theTrackInfo = MyAddTVTrack.getTrackInfo();
  //theTrackInfo2= MyAddTVTrack.getTrackInfo2();

  //TrackInfo = theTrackInfo;
  //TrackInfo2= theTrackInfo2;

  if (debug) cout <<"[TTVF] vertices found: "<<newVertices.size()<<endl;
  return newVertices;

}

virtual std::vector<TransientVertex> ConfigurableTertiaryTracksVertexFinder::vertices ( const std::vector< reco::TransientTrack > & ) const [virtual]

Reconstruct vertices

Implements VertexReconstructor.

Referenced by TertiaryTracksVertexFinder::vertices().

virtual std::vector<TransientVertex> ConfigurableTertiaryTracksVertexFinder::vertices	(	const std::vector< reco::TransientTrack > &	,
		const TransientVertex &	pv
	)		const `[virtual]`