TrigTrackSelector Class Reference

L2 tracks. More...

#include <TrigTrackSelector.h>

Inheritance diagram for TrigTrackSelector:

Collaboration diagram for TrigTrackSelector:

Public Types
typedef std::vector< ElementLink< xAOD::TrackParticleContainer > >	TrackParticleLinks_t

Public Member Functions
	TrigTrackSelector (TrackFilter *selector)
	use a radius of 47 mm corresponding to the Run 1 pixel inner radius For the IBL it should be 32 mm, but this was kept as 47 mm for consistency of the definition.
	TrigTrackSelector (TrackFilter *selector, double radius)
	TrigTrackSelector (TrackFilter *selector, int selectPdgId, int selectParentPdgId)
	TrigTrackSelector (TrackFilter *selector, double radius, int selectPdgId, int selectParentPdgId)
template<typename T>
	TrigTrackSelector (T &tm, TrackFilter *selector, double radius=default_radius(), int selectPdgId=0, int selectParentPdgId=0)
template<typename T>
	TrigTrackSelector (T tbegin, T tend, TrackFilter *selector, double radius=default_radius(), int selectPdgId=0, int selectParentPdgId=0)
	~TrigTrackSelector ()
virtual TrackSelector *	clone () override
void	setBeamline (double x, double y, double z=0)
double	getBeamX () const
double	getBeamY () const
double	getBeamZ () const
void	correctTracks (bool b=true)
virtual void	clear () override
bool	selectTrack (const TrigInDetTrack *track, const TrigInDetTrackTruthMap *truthMap=0)
	neater code to make use of vector function also for a single ancestor pdgid, instead of the full code duplication, but less efficienct as it then needs to create a single element vector for each particle, to avoid the code duplication.
void	selectTracks (const TrigInDetTrackCollection *trigtracks, const TrigInDetTrackTruthMap *truthMap=0)
bool	selectTrack (const Rec::TrackParticle *track)
void	selectTracks (const Rec::TrackParticleContainer *trigtracks)
void	selectTracks (const TruthParticleContainer *truthtracks)
void	selectTracks (const xAOD::TruthParticleContainer *truthtracks)
void	truthBeamline (const xAOD::TruthParticleContainer *truthtracks, double &x0, double &y0)
	extract all the tracks from a xAOD::TruthParticle collection and histogram the x and y production coordinates to determine the event "beamline" position
bool	selectTrack (HepMC::ConstGenParticlePtr track)
bool	selectTrack (const TruthParticle &track)
bool	selectTrack (const TruthParticle *track)
bool	selectTrack (const xAOD::TruthParticle *track, double x0=0, double y0=0)
TIDA::Track *	makeTrack (HepMC::ConstGenParticlePtr track)
TIDA::Track *	makeTrack (const TruthParticle *track, unsigned long tid=0)
bool	selectTrack (const Trk::Track *track)
void	selectTracks (const TrackCollection *trigtracks)
bool	selectTrack (const xAOD::TrackParticle *track, void *=0)
	legacy run 2 selector access
void	selectTracks (const xAOD::TrackParticleContainer *tracks, void *=0)
void	selectTracks (xAOD::TrackParticleContainer::const_iterator trackitr, xAOD::TrackParticleContainer::const_iterator trackend, void *=0)
	MT selector access.
void	selectTracks (const TrackParticleLinks_t &tracks)
void	correctToBeamline (double &z0, double &dz0, double &d0, double &dd0, double theta, double phi)
virtual bool	addTrack (TIDA::Track *t, bool(*f)(const TIDA::Track *)=0)
virtual void	addTracks (std::vector< TIDA::Track * > &t, bool(*f)(const TIDA::Track *)=0)
const std::vector< TIDA::Track * > &	tracks () const
std::vector< TIDA::Track * >	tracks (TrackFilter *selector) const
unsigned	size () const
void	delete_track (TIDA::Track *t)

Static Public Member Functions
static double	default_radius ()

Protected Member Functions
const xAOD::TruthParticle *	fromAncestor (const int pdg_id, const xAOD::TruthParticle *p) const
	recursive functions to identify whether a particle comes from some particle of a specific PDG ID, or some numb er of specified PDG IDs.
const xAOD::TruthParticle *	fromAncestor (const std::vector< int > &ids, const xAOD::TruthParticle *p) const
	recursive function to identify whether a particle comes from some other number of ancestor particles, with the pdgids passed in as a vector.
virtual void	cleanup (TIDA::Track *)

Protected Attributes
std::vector< TIDA::Track * >	m_tracks
TrackFilter *	m_selector

Private Attributes
unsigned long	m_id
double	m_xBeam
double	m_yBeam
double	m_zBeam
bool	m_correctTrkTracks
double	m_radius
int	m_selectPdgId
int	m_selectParentPdgId

Static Private Attributes
static const double	s_default_radius = 47
	NB: This was 47 for Run 2, but with the addition of the IBL it should be 32 It was kept at 47 for all Run 2 and migration to MT, but for Run 3 we really want it changed to be 32.

Detailed Description

L2 tracks.

TruthMap offline and EF MC truth FrameWork includes #include "GaudiKernel/ServiceHandle.h" #include "GaudiKernel/IPartPropSvc.h" absolutely pathetic!! write my own class TrigInDetAnalysis classes

Definition at line 58 of file TrigTrackSelector.h.

Member Typedef Documentation

TrackParticleLinks_t

typedef std::vector< ElementLink< xAOD::TrackParticleContainer > > TrigTrackSelector::TrackParticleLinks_t

Definition at line 62 of file TrigTrackSelector.h.

Constructor & Destructor Documentation

TrigTrackSelector() [1/6]

TrigTrackSelector::TrigTrackSelector ( TrackFilter * selector )

inline

use a radius of 47 mm corresponding to the Run 1 pixel inner radius For the IBL it should be 32 mm, but this was kept as 47 mm for consistency of the definition.

this should be changed to 32 mm for Run 3

Definition at line 70 of file TrigTrackSelector.h.

: TrigTrackSelector( selector, default_radius(), 0, 0 ) { } 

TrigTrackSelector() [2/6]

TrigTrackSelector::TrigTrackSelector ( TrackFilter * selector, double radius )

inline

Definition at line 72 of file TrigTrackSelector.h.

: TrigTrackSelector( selector, radius, 0, 0 ) { } 

TrigTrackSelector() [3/6]

TrigTrackSelector::TrigTrackSelector ( TrackFilter * selector, int selectPdgId, int selectParentPdgId )

inline

Definition at line 74 of file TrigTrackSelector.h.

: TrigTrackSelector( selector, default_radius(), selectPdgId, selectParentPdgId ) { } 

TrigTrackSelector() [4/6]

TrigTrackSelector::TrigTrackSelector	(	TrackFilter *	selector,
		double	radius,
		int	selectPdgId,
		int	selectParentPdgId )

Definition at line 24 of file TrigTrackSelector.cxx.

                                                                                                                   : 
    TrackSelector(selector), m_id(0), m_xBeam(0), m_yBeam(0), m_zBeam(0),
    m_correctTrkTracks(false), 
    m_radius(radius), m_selectPdgId(selectPdgId), m_selectParentPdgId(selectParentPdgId) { 
} 

TrigTrackSelector() [5/6]

template<typename T>

TrigTrackSelector::TrigTrackSelector ( T & tm, TrackFilter * selector, double radius = default_radius(), int selectPdgId = 0, int selectParentPdgId = 0 )

inline

Definition at line 79 of file TrigTrackSelector.h.

                                                                                                                                 : 
    TrigTrackSelector( tm.begin(), tm.end(), selector, radius, selectPdgId, selectParentPdgId ) { 
  }

TrigTrackSelector() [6/6]

template<typename T>

TrigTrackSelector::TrigTrackSelector ( T tbegin, T tend, TrackFilter * selector, double radius = default_radius(), int selectPdgId = 0, int selectParentPdgId = 0 )

inline

Definition at line 84 of file TrigTrackSelector.h.

                                                                                                                                            : 
    TrigTrackSelector( selector, radius, selectPdgId, selectParentPdgId ) { 
    selectTracks( tbegin, tend );
  }

~TrigTrackSelector()

TrigTrackSelector::~TrigTrackSelector ( )

inline

Definition at line 90 of file TrigTrackSelector.h.

{ clear(); }

Member Function Documentation

addTrack()

virtual bool TrackSelector::addTrack ( TIDA::Track * t, bool(* f )(const TIDA::Track *) = 0 )

inlinevirtualinherited

Definition at line 35 of file Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h.

                                                                      {
    //std::cout << "addtrack()  before f: t  " << *t << " " << size() << "\t f  " << f << std::endl;
    if ( f==0 ) { 
      if ( m_selector && m_selector->select(t) )  { m_tracks.push_back(t); return true; }
      else                                      { cleanup(t); }
    }
    else { 
      if ( f(t) )  { m_tracks.push_back(t); return true; } 
      else         { cleanup(t); } 
    }
    return false;
  }

addTracks()

virtual void TrackSelector::addTracks ( std::vector< TIDA::Track * > & t, bool(* f )(const TIDA::Track *) = 0 )

inlinevirtualinherited

Definition at line 48 of file Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h.

                                                                                   {
    for ( size_t i=0 ; i<t.size() ; i++ ) addTrack( t[i], f ); 
  }  

cleanup()

virtual void TrackSelector::cleanup ( TIDA::Track * )

inlineprotectedvirtualinherited

Definition at line 77 of file Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h.

{ } 

clear()

virtual void TrigTrackSelector::clear ( )

inlineoverridevirtual

Reimplemented from TrackSelector.

Definition at line 102 of file TrigTrackSelector.h.

{ for ( size_t i=m_tracks.size() ; i-- ; ) delete m_tracks[i]; m_tracks.clear(); }   

clone()

virtual TrackSelector * TrigTrackSelector::clone ( )

inlineoverridevirtual

Implements TrackSelector.

Definition at line 92 of file TrigTrackSelector.h.

{ return new TrigTrackSelector(*this); }

correctToBeamline()

void TrigTrackSelector::correctToBeamline	(	double &	z0,
		double &	dz0,
		double &	d0,
		double &	dd0,
		double	theta,
		double	phi )

make sure that users have set the beamline parameters

error estimates

Definition at line 1276 of file TrigTrackSelector.cxx.

                                                    {
    
    
    //   if ( m_first ) { 
    //     if ( m_xBeam==0 && m_yBeam==0 ) { 
    //      std::cerr << "TrigTrackSelector::correctToBeamline() WARNING -- Beamline set to (0,0) -- WARNING" << std::endl;  
    //     }
    //     else { 
    //     std::cout << "TrigTrackSelector::correctToBeamline() Beamline set to " << m_xBeam << " " << m_yBeam << std::endl;  
    //     }
    //     m_first = false;
    //   }
    
    
    //      double theta = 2*std::atan( exp( (-1)*eta ) );
    double z0t = z0 + ((std::cos(phi)*m_xBeam + std::sin(phi)*m_yBeam)/std::tan(theta));    
    double a0t = d0 +   std::sin(phi)*m_xBeam - std::cos(phi)*m_yBeam; 
    
    double dz0t = dz0 + ((std::cos(phi)*m_xBeam + std::sin(phi)*m_yBeam)/std::tan(theta));
    double da0t = dd0 +   std::sin(phi)*m_xBeam - std::cos(phi)*m_yBeam;
    
    z0  = z0t;
    d0  = a0t;
    
    dz0 = dz0t;
    dd0 = da0t;
}

correctTracks()

void TrigTrackSelector::correctTracks ( bool b = true )

inline

Definition at line 100 of file TrigTrackSelector.h.

{ m_correctTrkTracks = b; } 

default_radius()

double TrigTrackSelector::default_radius ( )

inlinestatic

Definition at line 188 of file TrigTrackSelector.h.

{ return s_default_radius; }

delete_track()

void TrackSelector::delete_track ( TIDA::Track * t )

inlineinherited

Definition at line 69 of file Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h.

                                    { 
    m_tracks.erase( std::remove( m_tracks.begin(), m_tracks.end(), t ), m_tracks.end() );
  }

fromAncestor() [1/2]

const xAOD::TruthParticle * TrigTrackSelector::fromAncestor ( const int pdg_id, const xAOD::TruthParticle * p ) const

protected

recursive functions to identify whether a particle comes from some particle of a specific PDG ID, or some numb er of specified PDG IDs.

recursive function to identify whether a particle comes from some other ancestor particle at any point in its history - checks its parents and then function calls itself to check their parents etc., until it either finds the pdgid it is looking for, or it has no more ancestors if it finds an appropriate ancestor it returns the pointer to it, otherwise it returns a nullptr

finds an appropriate ancestor and returns the pointer to it, otherwise it just returns a nullptr

Definition at line 39 of file TrigTrackSelector.cxx.

                                                                                                            { 
  if ( p==nullptr ) return nullptr;
  if ( MC::isElectron(p) || MC::isMuon(p) ) return nullptr; //don't want light leptons from eg tau decays - they are found directly
  if ( p->absPdgId()==pdg_id ) {
    return p;   // recursive stopping conditions
  }
  auto vertex = p->prodVtx();
  if ( vertex == nullptr ) {
    return nullptr; // has no production vertex !!!
  }
  if ( vertex->nIncomingParticles() < 1 ) {
    return nullptr;  // recursive stopping conditions
  }
  for( unsigned ip = 0; ip < vertex->nIncomingParticles(); ip++ ) {
    auto* in = vertex->incomingParticle(ip);
    auto parent = fromAncestor( pdg_id, in);
    if ( parent!=nullptr ) { 
      if (parent->absPdgId()==pdg_id) return parent;
    }
  }
  
  return nullptr;
}

fromAncestor() [2/2]

const xAOD::TruthParticle * TrigTrackSelector::fromAncestor ( const std::vector< int > & ids, const xAOD::TruthParticle * p ) const

protected

recursive function to identify whether a particle comes from some other number of ancestor particles, with the pdgids passed in as a vector.

Any such ancestors at point in it's history will do - checks it's parents and then calls the function calls itself to check the parent's parents etc, until it either finds one of the pdgids it is looking for, or it has no more ancestors if it finds an appropriate ancestor it returns the pointer to it, otherwise it returns a nullptr

Definition at line 73 of file TrigTrackSelector.cxx.

                                                                                                                       { 
  if ( p==nullptr ) return nullptr;
  if (MC::isElectron(p) || MC::isMuon(p)) return nullptr; //don't want light leptons from eg tau decays - they are found directly
  for ( size_t i=ids.size() ; i-- ; ) { 
    if ( p->absPdgId()==ids[i] ) return p; // recursive stopping conditions
  }
 
  auto vertex = p->prodVtx();
  if ( vertex == nullptr ) return nullptr; // has no production vertex !!!
 
  if ( vertex->nIncomingParticles()<1 ) return nullptr;  // recursive stopping conditions
 
  for( unsigned ip = 0; ip < vertex->nIncomingParticles(); ip++ ) {
    auto* in = vertex->incomingParticle(ip);
    auto parent = fromAncestor( ids, in);
    if ( parent!=nullptr ) { 
      for ( size_t i=ids.size() ; i-- ; ) { 
    if ( parent->absPdgId()==ids[i] ) return parent; 
      }
    }
  }  
 
  return nullptr;
}

getBeamX()

double TrigTrackSelector::getBeamX ( ) const

inline

Definition at line 96 of file TrigTrackSelector.h.

{ return m_xBeam; } 

getBeamY()

double TrigTrackSelector::getBeamY ( ) const

inline

Definition at line 97 of file TrigTrackSelector.h.

{ return m_yBeam; } 

getBeamZ()

double TrigTrackSelector::getBeamZ ( ) const

inline

Definition at line 98 of file TrigTrackSelector.h.

{ return m_zBeam; } 

makeTrack() [1/2]

TIDA::Track * TrigTrackSelector::makeTrack	(	const TruthParticle *	track,
		unsigned long	tid = 0 )

check for final state

the is our new "final state" requirement the inner and outer radii are in some sense arbitrary - these correspond to an envelope around the pixel detector, so the track must pass through the entire pixel detector NB: In order to ensure you don't miss any tracks they really need to be the same ie if a track brems in your "volume" then you will miss that track, and also the resulting track, even if it is a high et track

keep anything over 10 GeV with the old requirement

avoid default (unset) TruthParticle charge

only use charged tracks

what a faff - refuse to mess about with the classes to swim tracks etc - why can't they just encode this sort of information in the class!! It's not as if it doesn't actually have members for anything else useless!! Classes should be designed for ease of use !!!!

is there a sign issue here ?

correct back to the beamline

static HepPDT::ParticleDataTable* m_pdt = new ParticleDataTable();

// Get the Particle Properties Service if ( m_pdt==0 ) {
ServiceHandle<IPartPropSvc> partPropSvc("PartPropSvc", "TrigTestMonToolAC"); // , name()); if ( !partPropSvc.retrieve().isSuccess() ) { m_pdt = partPropSvc->PDT(); } else { std::cerr << " Could not initialize Particle Properties Service" << std::endl; return; // StatusCode::FAILURE; } }

how about storing uniqueID/status/pidg info?

this isn't good!! but it will do for testing

probably won't work either

get the production vertex for the z0 and d0 but should we store the z0 of the production vertex? or swim the track to the perigee point with respect to 0,0 and use the parameters there?

what to do with these???

Definition at line 700 of file TrigTrackSelector.cxx.

                                                                                       { 
  
    if ( track==0 ) return 0; 
    if ( !MC::isStable(track)) return 0;   
 
 
    double phi = track->phi();
    double eta = track->eta();

    double xp[3] = { 0, 0, 0 };
 
    if ( track->genParticle()->production_vertex() ) { 
      xp[0] = track->genParticle()->production_vertex()->position().x();
      xp[1] = track->genParticle()->production_vertex()->position().y();
      xp[2] = track->genParticle()->production_vertex()->position().z();
    }
 
    // CHANGED BY JK - z0 with respect to (0,0)
    //    double z0 = xp[2];
    double theta = 2*std::atan( exp( (-1)*eta ) );
    double z0 = xp[2] - (xp[0]*std::cos(phi) + xp[1]*std::sin(phi))/std::tan(theta);
    
    double xd[3] = { 0, 0, 0 };
 
    if ( track->genParticle()->end_vertex() ) { 
      xd[0] = track->genParticle()->end_vertex()->position().x();
      xd[1] = track->genParticle()->end_vertex()->position().y();
      xd[2] = track->genParticle()->end_vertex()->position().z();
    }
 
    double rp = std::sqrt( xp[0]*xp[0] + xp[1]*xp[1] ); 
    double rd = std::sqrt( xd[0]*xd[0] + xd[1]*xd[1] ); 
 
    
    bool final_state = false; 

    const double inner_radius = m_radius;
    const double outer_radius = m_radius;
    if ( ( track->genParticle()->production_vertex() && rp<=inner_radius ) && 
     ( track->genParticle()->end_vertex()==0 || rd>outer_radius ) ) final_state = true; 
     
    
    if ( !final_state ) return 0; 
    
 
    double q = track->charge();
    
    static const particleType ptype;

    if ( q==-999 ) q = ptype.charge( track->pdgId() );

    if ( q==0 ) return 0;
    
    double pT  = q*track->pt();
 
 
    double d0  = 0;


 
 
    // CHANGED BY JK - d0 with respect to (0,0)
    // d0 = q*rp*std::sin(phi);
    d0  =  xp[1]*std::cos(phi) -  xp[0]*std::sin(phi);
 

 
    double dz0 = 0;
    double dd0 = 0;
   
    correctToBeamline( z0, dz0, d0, dd0, theta, phi );
 
 

 

    int author  = track->pdgId();   
    int uniqueID = HepMC::uniqueID(track); 
 
 
    unsigned long id = (unsigned long)track;
    if ( tid!=0 ) id = tid;

 

        
    //    std::cout << "\t\t\tSUTT Truth track" 
    //              << "\teta=" << eta  // << " +- " << (*trackitr)->params()->deta()
    //              << "\tphi=" << phi  // << " +- " << (*trackitr)->params()->dphi()
    //              << "\tz0="  << z0 
    //              << "\tpT="  << pT // << "\t( " << 1/qoverp << ")"
    //              << "\td0="  << d0
    //              << "\tauthor=" << author
    //              << std::endl;
    
 
 
    TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, 0, 0,
                                     0, 0, 0, 0, 0,
                                     0, 0, 0, 0,
                                     0, 0, 0, 0,
                                     author, false, uniqueID, -1,
                                     false, 
                                     id ) ;  
 
    return t;
 
}

makeTrack() [2/2]

TIDA::Track * TrigTrackSelector::makeTrack

(

HepMC::ConstGenParticlePtr

track

)

Definition at line 688 of file TrigTrackSelector.cxx.

                                                                       { 
//AV Using memory to get some value is not a good idea. This is not a repruducible/portable way, but I leave it as is.
#ifdef HEPMC3
    unsigned long id = (unsigned long)(track.get());
#else
    unsigned long id = (unsigned long)track;
#endif
    TruthParticle t = TruthParticle(track); 
    return  makeTrack( &t, id );
}

selectTrack() [1/8]

bool TrigTrackSelector::selectTrack

(

const Rec::TrackParticle *

track

)

temporarily remove interaction mean z position

<< check this <<–

check this !!!

now some ridiculous code to get the author of the TrackParticle (offline) tracks

Definition at line 200 of file TrigTrackSelector.cxx.

                                                                   { 
        
    // do the track extraction stuff here....
 
    static const int hpmap[20] = { 0, 1, 2,  7, 8, 9,  3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
 
    if ( track ) { 
        
#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
      const Trk::MeasuredPerigee* measPer = track->measuredPerigee();
#else
      const Trk::Perigee* measPer = track->measuredPerigee();
#endif
      //     CLHEP::HepVector perigeeParams = measPer->parameters();
      
      double pT    = measPer->pT(); 
      double eta   = measPer->eta();
      double phi   = measPer->parameters()[Trk::phi0];
      double z0    = measPer->parameters()[Trk::z0] + m_zBeam; 
      double d0    = measPer->parameters()[Trk::d0];
 
      double theta = measPer->parameters()[Trk::theta];
      double  p    = 1/measPer->parameters()[Trk::qOverP];
  
      // AAARCH!!!!! the TrackParticle pT is NOT SIGNED!!!! ( I ask you! ) 
      if ( measPer->parameters()[Trk::qOverP]<0 && pT>0 ) pT *= -1;
 
#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
      const Trk::ErrorMatrix err = measPer->localErrorMatrix();
      double dtheta = err.error(Trk::theta);
      double dqovp  = err.error(Trk::qOverP); 
      double covthetaOvP  = err.covValue(Trk::qOverP,Trk::theta);  
#else
      double dtheta      = std::sqrt((*measPer->covariance())(Trk::theta,Trk::theta));
      double dqovp       = std::sqrt((*measPer->covariance())(Trk::qOverP,Trk::qOverP)); 
      double covthetaOvP = (*measPer->covariance())(Trk::qOverP,Trk::theta);  
#endif
 
 
      double deta = 0.5*dtheta/(std::cos(0.5*theta)*std::cos(0.5*theta)*std::tan(0.5*theta));   
 
#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
      double dphi = err.error(Trk::phi0);
      double dz0  = err.error(Trk::z0);
      double dd0  = err.error(Trk::d0);
#else
      double dphi = std::sqrt((*measPer->covariance())(Trk::phi0,Trk::phi0));
      double dz0  = std::sqrt((*measPer->covariance())(Trk::z0,Trk::z0));
      double dd0  = std::sqrt((*measPer->covariance())(Trk::d0,Trk::d0));
#endif
 
      double dpT  = 0; 
 
 
      double sintheta = std::sin(theta);
      double costheta = std::cos(theta);
      double dpt2 = (p*p*sintheta)*(p*p*sintheta)*dqovp*dqovp + (p*costheta)*(p*costheta)*dtheta*dtheta - 2*(p*p*sintheta)*(p*costheta)*covthetaOvP;   
 
      if ( dpt2>0 ) dpT = std::sqrt( dpt2 );
 
      // Check number of hits
      // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really 
      // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersecting
      // stereo clusters making a spacepoint are two "hits"
      const Trk::TrackSummary *summary = track->trackSummary();
      int nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits); 
      int nPixelHits  = 2*summary->get(Trk::numberOfPixelHits);  
      int nSctHits    = summary->get(Trk::numberOfSCTHits); 
      int nStrawHits  = summary->get(Trk::numberOfTRTHits);
      int nTrHits     = summary->get(Trk::numberOfTRTHighThresholdHits);
 
      int nSiHits     = nPixelHits + nSctHits;
      bool expectBL   = false; // Not stored for Rec::TrackParticle
 
      const Trk::FitQuality *quality   = track->fitQuality();
      double chi2 = quality->chiSquared();
      double dof  = quality->numberDoF();
      
      unsigned bitmap = 0;
 
      
      unsigned long id = (unsigned long)track;
 
      for ( int ih=0 ; ih<20 ; ih++ ) {
    if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );    
      }

      
      //      std::cout << "fetching author info :" << track->info().trackFitter() << ":" 
      //        << track->info().dumpInfo() << ": bm 0x" << std::hex << bitmap << std::dec << std::endl;
      
      int fitter = track->info().trackFitter();
      std::string dumpinfo = track->info().dumpInfo();
      
      int trackAuthor = -1; 
      if ( fitter>0 && fitter<Trk::TrackInfo::NumberOfTrackFitters ) {
    if      ( dumpinfo.find("TRTStandalone")!=std::string::npos)        trackAuthor = 2;
    else if ( dumpinfo.find("TRTSeededTrackFinder")!=std::string::npos) trackAuthor = 1;
    else                                                                trackAuthor = 0;
      }
      
#if 0
      std::cout << "\t\t\tSUTT TP track" 
        << "\teta=" << eta  // << " +- " << (*trackitr)->params()->deta()
        << "\tphi=" << phi  // << " +- " << (*trackitr)->params()->dphi()
        << "\tz0="  << z0 
        << "\tpT="  << pT // << "\t( " << 1/qoverp << ")"
        << "\td0="  << d0
        << "\tNsi=" << nSiHits 
        << "\tNtrt=" << nTrHits 
        << "\tNstr=" << nStrawHits
        << "\tauthor=" << trackAuthor
        << std::endl;
#endif  
 
      // Create and save Track
      
      TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
                                       deta, dphi, dz0, dd0, dpT,
                                       nBlayerHits, nPixelHits, nSctHits, nSiHits,
                                       nStrawHits, nTrHits, bitmap, 0,
                                       trackAuthor, false, -1, -1,  
                                       expectBL, id) ;  
 
      //      std::cout << "SUTT TP track " << *t << "\t0x" << std::hex << bitmap << std::dec << std::endl; 
      
      if ( !addTrack( t ) ){
    delete t;
    return false;
      }
      return true;
      
    }
    return false;
}

selectTrack() [2/8]

bool TrigTrackSelector::selectTrack	(	const TrigInDetTrack *	track,
		const TrigInDetTrackTruthMap *	truthMap = 0 )

neater code to make use of vector function also for a single ancestor pdgid, instead of the full code duplication, but less efficienct as it then needs to create a single element vector for each particle, to avoid the code duplication.

Perhaps something can be done for the compiler optimisation to realise this and somehow
reduce the opverhead of creating the single int vector - need to investiogate this more thoroghly before including this, as it would be a much neater solution

not definied ofr TrigInDetTracks

Definition at line 116 of file TrigTrackSelector.cxx.

                                                                                                         {     
    // do the track extraction stuff here....
    if ( track ) { 
 
    double eta    = track->param()->eta();
    double phi    = track->param()->phi0();
    double z0     = track->param()->z0(); 
    double pT     = track->param()->pT(); 
    double d0     = track->param()->a0();
 
    double deta    = track->param()->eeta();
    double dphi    = track->param()->ephi0();
    double dz0     = track->param()->ez0(); 
    double dpT     = track->param()->epT(); 
    double dd0     = track->param()->ea0();
 
    double theta = 2*std::atan( exp( (-1)*eta ) );
    correctToBeamline( z0, dz0, d0, dd0, theta, phi );
    
    int   algoid  = track->algorithmId();        
 
    int nBlayerHits = (track->HitPattern() & 0x1);
    int nPixelHits  = 2 * track->NPixelSpacePoints();  // NB: for comparison with offline 
    int nSctHits    = 2 * track->NSCT_SpacePoints();   //     a spacepoint is 2 "hits"
    int nStrawHits  =  track->NStrawHits();
    int nTrHits     =  track->NTRHits();
 
    int nSiHits     = nPixelHits + nSctHits;
 
    bool expectBL   = false;                           //not filled in 
 
    unsigned long id = (unsigned long)track;
 
    unsigned hitPattern = track->HitPattern();
    unsigned multiPattern = 0;
 
    double chi2    = track->chi2();
    double dof     = 0; 
 
    bool truth = false;
    int match_uniqueID = HepMC::INVALID_PARTICLE_ID;
    
    if ( truthMap ) { 
      const TrigInDetTrackTruth* trackTruth = truthMap->truth(track);
      if (trackTruth!=0 && trackTruth->nrMatches() > 0) {
        match_uniqueID = HepMC::uniqueID(trackTruth->bestSiMatch());
        truth = true;
      }
    }
    
    
    TIDA::Track* t = new TIDA::Track(  eta,  phi,  z0,  d0,  pT, chi2, dof, 
                                     deta, dphi, dz0, dd0, dpT, 
                                     nBlayerHits, nPixelHits, nSctHits, nSiHits, 
                                     nStrawHits, nTrHits, 
                                     hitPattern, multiPattern, 
                                     algoid, truth, -1, match_uniqueID,
                                     expectBL, id) ; 
    
    //  std::cout << "SUTT ID track " << *t << "\t0x" << std::hex << track->HitPattern() << std::dec << std::endl;
    
    if ( !addTrack( t ) ){
      delete t;
      return false;
    }   
    return true;
    }
    return false;
}

selectTrack() [3/8]

bool TrigTrackSelector::selectTrack

(

const Trk::Track *

track

)

a covariance!

check this <<–

now some ridiculous code to get the author of the TrackParticle (offline) tracks

Definition at line 863 of file TrigTrackSelector.cxx.

                                                           { 
        
    // do the track extraction stuff here....
 
    static const int hpmap[20] = { 0, 1, 2,  7, 8, 9,  3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
    //    std::cout << __FILE__<< " " <<__LINE__<<std::endl;
 
    if ( track ) { 
        
      // const Trk::Perigee* startPerigee = track->perigeeParameters();
 
#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
      const Trk::MeasuredPerigee* startPerigee = dynamic_cast<const Trk::MeasuredPerigee*>(track->perigeeParameters());
      //      const Trk::MeasuredPerigee* measPer = startPerigee; // just out of laziness
#else
      const Trk::Perigee* startPerigee = track->perigeeParameters();
      const Trk::Perigee* measPer = startPerigee; // just out of laziness
#endif
 
 
      //      CLHEP::HepVector perigeeParams = measPer->parameters();      
      //      double pT    = measPer->pT(); 
      //      double eta   = measPer->eta();
      //      double phi   = perigeeParams[Trk::phi0];
      //      double z0    = perigeeParams[Trk::z0];
      //      double d0    = perigeeParams[Trk::d0];
      //      // AAARCH!!!!! the TrackParticle pT is NOT SIGNED!!!! ( I ask you! ) 
      //      if ( perigeeParams[Trk::qOverP]<0 ) pT *= -1;
      //      std::cout <<pT1<<" pt1vspT "<<pT<<std::endl;
      
      if (startPerigee){
     
    double theta   = startPerigee->parameters()[Trk::theta];  
    double p       = 1/startPerigee->parameters()[Trk::qOverP];
    double qOverPt = startPerigee->parameters()[Trk::qOverP]/std::sin(theta);
    double charge  = startPerigee->charge();
    double eta     = startPerigee->eta();
    double phi     = startPerigee->parameters()[Trk::phi0];
    double z0      = startPerigee->parameters()[Trk::z0];
    double d0      = startPerigee->parameters()[Trk::d0];
    //  double pT = (1./qOverPt)*(charge);   
    double pT      = (1./qOverPt); // always use signed PT
    
    if ( charge<0 && pT>0 ) pT *= -1;
    if ( charge<0 && p>0 )  p  *= -1;
    
 
 
#ifdef TRKPARAMETERS_MEASUREDPERIGEE_H
    const Trk::ErrorMatrix err = startPerigee->localErrorMatrix();
    double dtheta = err.error(Trk::theta);
    double dqovp  = err.error(Trk::qOverP); 
    double covthetaOvP  = err.covValue(Trk::qOverP,Trk::theta);  
    
    double dphi = err.error(Trk::phi0);
    double dz0  = err.error(Trk::z0);
    double dd0  = err.error(Trk::d0);
#else
    double dtheta      = std::sqrt((*measPer->covariance())(Trk::theta,Trk::theta));
    double dqovp       = std::sqrt((*measPer->covariance())(Trk::qOverP,Trk::qOverP)); 
    double covthetaOvP = (*measPer->covariance())(Trk::qOverP,Trk::theta);    
 
    double dphi = std::sqrt((*measPer->covariance())(Trk::phi0,Trk::phi0));
    double dz0  = std::sqrt((*measPer->covariance())(Trk::z0,Trk::z0));
    double dd0  = std::sqrt((*measPer->covariance())(Trk::d0,Trk::d0));
#endif  
    
    double deta = 0.5*dtheta/(std::cos(0.5*theta)*std::cos(0.5*theta)*std::tan(0.5*theta));  
 
 
    if ( m_correctTrkTracks ) correctToBeamline( z0, dz0, d0, dd0, theta, phi );
 
    double dpT  = 0; 
    
    
    double sintheta = std::sin(theta);
    double costheta = std::cos(theta);
    double dpT2 = (p*p*sintheta)*(p*p*sintheta)*dqovp*dqovp + (p*costheta)*(p*costheta)*dtheta*dtheta - 2*(p*p*sintheta)*(p*costheta)*covthetaOvP;
    
    if ( dpT2>0 ) dpT = std::sqrt( dpT2 );
 
    // Check number of hits
    // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really 
    // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
    // stereo clusters making a spacepoint are two "hits"
    // const Trk::TrackSummary *summary = dynamic_cast<const Trk::TrackSummary*>(track->trackSummary());
    //ToolHandle< Trk::ITrackSummaryTool > m_trackSumTool;
    //m_trackSumTool  = ToolHandle<Trk::ITrackSummaryTool>("Trk::TrackSummaryTool/InDetTrackSummaryTool");
    //const Trk::TrackSummary* summary = NULL;       
    //summary = m_trackSumTool->createSummary(*track); 
 
    const Trk::TrackSummary * summary = track->trackSummary();
    int nBlayerHits = 0;
    int nPixelHits  = 0;
    int nSctHits    = 0;
    int nStrawHits  = 0;
    int nTrHits     = 0;
    int nSiHits     = 0;
    bool expectBL = false; // Not stored for Trk::Track
    unsigned bitmap = 0;
 
    if(summary==0){
            std::cerr << "Could not create TrackSummary  - Track will likely fail hits requirements" << std::endl;
    }    
    else{      
            nBlayerHits = 2*summary->get(Trk::numberOfBLayerHits); 
            nPixelHits  = 2*summary->get(Trk::numberOfPixelHits);  
        nSctHits    = summary->get(Trk::numberOfSCTHits); 
            nStrawHits  = summary->get(Trk::numberOfTRTHits);
            nTrHits     = summary->get(Trk::numberOfTRTHighThresholdHits);
        nSiHits     = nPixelHits + nSctHits;
 
        for ( int ih=0 ; ih<20 ; ih++ ) {
          if ( summary->isHit(Trk::DetectorType(ih)) ) bitmap |= ( 1<<hpmap[ih] );  
        }
    }
 
    unsigned long id = (unsigned long)track;
    double chi2 = 0;
    double dof  = 0; 
    //const Trk::FitQuality *quality   = dynamic_cast<const Trk::FitQuality*>(track->fitQuality());
    const Trk::FitQuality *quality   = (track->fitQuality());
    if(quality==0) std::cerr << "Could not create FitQuality  - Track will likely fail hits requirements" << std::endl;
    else{
            chi2 = quality->chiSquared();          
        dof  = quality->numberDoF();
    }
 
    int trackAuthor = -1;

 
    //      std::cout << "fetching author info :" << track->info().trackFitter() << ":" 
    //      << track->info().dumpInfo() << ": bm 0x" << std::hex << bitmap << std::dec << std::endl;
 
    int fitter = track->info().trackFitter();
    //      std::string dumpinfo = track->info().dumpInfo();
 
    if ( fitter>0 && fitter<Trk::TrackInfo::NumberOfTrackFitters ) {
      if     ((track->info().dumpInfo()).find("TRTStandalone")        != std::string::npos) trackAuthor = 2;
      else if   ((track->info().dumpInfo()).find("TRTSeededTrackFinder") != std::string::npos) trackAuthor = 1;
      else                                                                                    trackAuthor = 0;
    }
 
  #if 0
    std::cout << "\t\t\tSUTT TP track" 
          << "\teta=" << eta  // << " +- " << (*trackitr)->params()->deta()
          << "\tphi=" << phi  // << " +- " << (*trackitr)->params()->dphi()
          << "\tz0="  << z0 
          << "\tpT="  << pT // << "\t( " << 1/qoverp << ")"
          << "\td0="  << d0
          << "\tNsi=" << nSiHits 
          << "\tNtrt=" << nTrHits 
          << "\tNstr=" << nStrawHits
          << "\tauthor=" << trackAuthor
          << std::endl;
  #endif    
    // Create and save Track      
    TIDA::Track* t = new TIDA::Track(eta, phi, z0, d0, pT, chi2, dof,
                                         deta, dphi, dz0, dd0, dpT,
                                         nBlayerHits, nPixelHits, nSctHits, nSiHits,
                                         nStrawHits, nTrHits, bitmap, 0,
                                         trackAuthor,  false, -1, -1,  
                                         expectBL, id) ;  
 
    if ( !addTrack( t ) ){
      delete t;
      return false;
    }
    return true;
    
    //std::cout << "SUTT TP track " << *t << "\t0x" << std::hex << bitmap << std::dec << std::endl; 
      }
    }
    
    return false;
}

selectTrack() [4/8]

bool TrigTrackSelector::selectTrack

(

const TruthParticle &

track

)

Definition at line 522 of file TrigTrackSelector.cxx.

                                                                {
 
  return selectTrack( &track );
  
}

selectTrack() [5/8]

bool TrigTrackSelector::selectTrack

(

const TruthParticle *

track

)

Definition at line 530 of file TrigTrackSelector.cxx.

                                                                { 
    TIDA::Track* t = makeTrack( track, m_id );
    if ( t == 0 ) return false;
    if ( !addTrack(t) ) {
      delete t;
      return false;
    }
    return true;
}

selectTrack() [6/8]

bool TrigTrackSelector::selectTrack	(	const xAOD::TrackParticle *	track,
		void *	= 0 )

legacy run 2 selector access

lazy just to avoid a find-replace of measPer to track

Grrrrr remove interaction mean z position!!!

protect against spurious qOverP = 0 BUT !!! this should NEVER happen !!! so shouldn't really have this failing silently, but not allowed to put in std::cerr output leave the original code in as this is mostly for a test

a covariance not an error

don't correct xaod tracks to the beamline

??? <<– check this

holes

cheat !! pack the holes into the hits so that eg int pixelholes() { return npix/1000; } int pixelhits() { return npix%1000; }

get the total number of holes as well

fit quality

Definition at line 1059 of file TrigTrackSelector.cxx.

                                                                           {
    // do the track extraction stuff here....
 
    //    static int hpmap[20] = { 0, 1, 2,  7, 8, 9,  3, 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19 };
 
    if ( track ) { 
        
      const xAOD::TrackParticle* measPer = track;
 
      //     CLHEP::HepVector perigeeParams = measPer->parameters();
      
      double pT    = measPer->pt(); 
      double eta   = measPer->eta();
      double phi   = measPer->phi0();
      double z0    = measPer->z0() + measPer->vz();  
      double d0    = measPer->d0();
 
      double theta = measPer->theta();

      // if ( measPer->qOverP()==0 ) return false;
      if ( measPer->qOverP()==0 ) throw std::runtime_error( "probable corrupted track - this should never happen" );
      double  p    = 1/measPer->qOverP();
  
      // AAARCH!!!!! the TrackParticle pT is NOT SIGNED!!!! ( I ask you! ) 
      if ( measPer->qOverP()<0 && pT>0 ) pT *= -1;
 
      double dtheta      = std::sqrt(measPer->definingParametersCovMatrix()(Trk::theta,Trk::theta));
      double dqovp       = std::sqrt(measPer->definingParametersCovMatrix()(Trk::qOverP,Trk::qOverP)); 
      double covthetaOvP = measPer->definingParametersCovMatrix()(Trk::qOverP,Trk::theta);  
 
      double deta = 0.5*dtheta/(std::cos(0.5*theta)*std::cos(0.5*theta)*std::tan(0.5*theta));  // ???? CHECK THIS <<--
 
      double dphi = std::sqrt(measPer->definingParametersCovMatrix()(Trk::phi0,Trk::phi0));
      double dz0  = std::sqrt(measPer->definingParametersCovMatrix()(Trk::z0,Trk::z0));
      double dd0  = std::sqrt(measPer->definingParametersCovMatrix()(Trk::d0,Trk::d0));
 
      double dpT  = 0; 
 

      //      if ( m_xBeam!=0 || m_yBeam!=0 ) correctToBeamline( z0, dz0, d0, dd0, theta, phi );
      
 
      double sintheta = std::sin(theta);
      double costheta = std::cos(theta);
      double dpt2 = (p*p*sintheta)*(p*p*sintheta)*dqovp*dqovp + (p*costheta)*(p*costheta)*dtheta*dtheta - 2*(p*p*sintheta)*(p*costheta)*covthetaOvP;   
 
      if ( dpt2>0 ) dpT = std::sqrt( dpt2 );
 
      // Check number of hits
      // NB: a spacepoint is two offline "hits", so a pixel spacepoint is really 
      // 2 "hits" and an offline SCT "hit" is really a 1D cluster, so two intersetcting
      // stereo clusters making a spacepoint are two "hits"
      
      uint8_t sum_nBlayerHits = 0;
      track->summaryValue( sum_nBlayerHits, xAOD::numberOfInnermostPixelLayerHits);
      int nBlayerHits = 2*sum_nBlayerHits;
 
      uint8_t  sum_nPixelHits = 0;
      track->summaryValue( sum_nPixelHits, xAOD::numberOfPixelHits);
      int nPixelHits = 2*sum_nPixelHits;
 
      uint8_t  sum_nSctHits = 0;
      track->summaryValue( sum_nSctHits, xAOD::numberOfSCTHits);
      int nSctHits = sum_nSctHits;
 
      uint8_t  sum_nStrawHits = 0;
      track->summaryValue( sum_nStrawHits, xAOD::numberOfTRTHits);
      int nStrawHits  = sum_nStrawHits;
 
      uint8_t  sum_nTrtHits = 0;
      track->summaryValue( sum_nTrtHits, xAOD::numberOfTRTHighThresholdHits);
      int nTrtHits  = sum_nTrtHits;
      
 
      uint8_t sum_expectBL  = 0;
      track->summaryValue( sum_expectBL, xAOD::expectInnermostPixelLayerHit);
      bool expectBL = ( sum_expectBL ? true : false );

 
      uint8_t sum_sctholes  = 0;
      track->summaryValue( sum_sctholes, xAOD::numberOfSCTHoles);
 
      uint8_t sum_pixholes  = 0;
      track->summaryValue( sum_pixholes, xAOD::numberOfPixelHoles);

 
      nSctHits   += 1000*sum_sctholes;
      nPixelHits += 1000*sum_pixholes;

      int nSiHits     = nPixelHits + nSctHits;

 
      double chi2 = track->chiSquared();
      double dof  = track->numberDoF();
      
      unsigned long id = (unsigned long)track;
 
      unsigned bitmap = track->hitPattern();
 
 
 
      double xbeam = track->vx(); 
      double ybeam = track->vy(); 
      double zbeam = track->vz(); 
      
      if ( xbeam!=getBeamX() || ybeam!=getBeamY() || zbeam!=getBeamZ() ) setBeamline( xbeam, ybeam, zbeam );
    
      int trackAuthor = 0;
 
      int fitter = track->trackFitter();
      std::bitset<xAOD::NumberOfTrackRecoInfo>  patternrec = track->patternRecoInfo();
 
      //int icount = 0;<- never used if section below is commented
      for ( unsigned ipr=patternrec.size() ; ipr-- ; ) { 
    if ( patternrec[ipr] ) {
      //icount++; <- never used if section below is commented
      trackAuthor |= (ipr >> 16);
      // static bool first = true;
      // if ( first && icount>1 ) { 
      //   std::cerr << "more than one pattern rec strategy " << ipr << "\t(suppressing further output)" << std::endl;
      //   first = false;
      // }
    }
      }
      
      trackAuthor |= fitter;
 
      //   if ( fitter>0 && fitter<Trk::TrackInfo::NumberOfTrackFitters ) {
      //    if      ( dumpinfo.find("TRTStandalone")!=std::string::npos)        trackAuthor = 2;
      //    else if ( dumpinfo.find("TRTSeededTrackFinder")!=std::string::npos) trackAuthor = 1;
      //    else                                                                trackAuthor = 0;
      //   }
      
#if 0
      std::cout << "\t\t\tSUTT TP track" 
        << "\teta=" << eta  // << " +- " << (*trackitr)->params()->deta()
        << "\tphi=" << phi  // << " +- " << (*trackitr)->params()->dphi()
        << "\tz0="  << z0 
        << "\tpT="  << pT // << "\t( " << 1/qoverp << ")"
        << "\td0="  << d0
        << "\tNsi=" << nSiHits 
    //      << "\tNtrt=" << nTrtHits 
    //      << "\tNstr=" << nStrawHits
        << "\tfitter=" << fitter
        << "\tauthor=" << trackAuthor
        << "\tVTX x " << track->vx() << "\ty " << track->vy() << "\tz " << track->vz() 
        << std::endl;
#endif  
 
      // Create and save Track
      
      TIDA::Track* t = new TIDA::Track( eta,   phi,  z0,  d0,  pT, chi2, dof,
                        deta,  dphi, dz0, dd0, dpT,
                        nBlayerHits, nPixelHits, nSctHits, nSiHits,
                        nStrawHits,  nTrtHits,   bitmap, 0,
                        trackAuthor,  false, -1, -1,  
                        expectBL, id) ;  
 
      //      std::cout << "SUTT TP track " << *t << "\t0x" << std::hex << bitmap << std::dec << std::endl; 
      
      if ( !addTrack( t ) ){
    delete t;
    return false;
      }
      return true;
      
    }
    return false;
     
 
}

selectTrack() [7/8]

bool TrigTrackSelector::selectTrack	(	const xAOD::TruthParticle *	track,
		double	x0 = 0,
		double	y0 = 0 )

lazy just to avoid a find-replace of measPer to track

avoid default (unset) TruthParticle charge

only use charged tracks

need to calculate the origin, and the beamline, and the d0 and z0 with respect to the beamline leave this in until we have checked whether everything is implemented correctly

need to calculate d0 and z0 correctly.

these are the d0 and z at the point of closest approach to x0, y0, the event "beamline", which is the best that we can do for the moment

the is our new "final state" requirement the inner and outer radii are in some sense arbitrary - these correspond to an envelope around the pixel detector, so the track must pass through the entire pixel detector NB: In order to ensure you don't miss any tracks they really need to be the same ie if a track brems in your "volume" then you will miss that track, and also the resulting track, even if it is a high et track

was hardcoded as 47 - now is set from the constructor

get the total number of holes as well

useful debug info - leave in

Definition at line 544 of file TrigTrackSelector.cxx.

                                                                                         { 
  if ( track ) { 
        
    if (!MC::isStable(track) ) return false;

    const xAOD::TruthParticle* measPer = track;
 
    double pT    = measPer->pt(); 
    double eta   = measPer->eta();
    double phi   = measPer->phi();
  
    // AAARCH!!!!! the TrackParticle pT is NOT SIGNED!!!! ( I ask you! ) 
    if ( measPer->charge()<0 && pT>0 ) pT *= -1;
    double q = track->charge();
 
    static const particleType ptype;

    if ( q==-999 ) q = ptype.charge( track->pdgId() );
    
    if ( q==0 ) return 0;

    //      double xbeam = getBeamX(); // track->vx(); 
    //      double ybeam = getBeamY(); // track->vy(); 
    //      double zbeam = getBeamZ(); // track->vz(); 
      
    if ( !track->hasProdVtx() ) return false; 

 
    double xp[3] = { measPer->prodVtx()->x(), measPer->prodVtx()->y(), measPer->prodVtx()->z() };
    double xb[3] = { xp[0]-x0, xp[1]-y0, measPer->prodVtx()->z() };
    double xd[3] = { 0, 0, 0 };
      
    if ( track->hasDecayVtx() ) { 
      xd[0] = track->decayVtx()->x();
      xd[1] = track->decayVtx()->y();
      xd[2] = track->decayVtx()->z();
    }
      
    double rp = std::sqrt( xp[0]*xp[0] + xp[1]*xp[1] ); 
    double rd = std::sqrt( xd[0]*xd[0] + xd[1]*xd[1] ); 
      
    double theta = 2*std::atan( std::exp( -eta ) );
    double z0 = xb[2] - (xb[0]*std::cos(phi) + xb[1]*std::sin(phi))/std::tan(theta);
    double d0 = xb[1]*std::cos(phi) -  xb[0]*std::sin(phi);
 
    bool final_state = false; 
      
    const double inner_radius = m_radius; 
    const double outer_radius = m_radius;
 
    if ( (  track->hasProdVtx() && rp<=inner_radius ) && 
     ( !track->hasDecayVtx() || rd>outer_radius ) ) final_state = true; 
      
    if ( !final_state ) return false; 
    
    double deta = 0;
    double dphi = 0;
    double dz0  = 0;
    double dd0  = 0;
 
    double dpT  = 0; 
 
    int nBlayerHits = 0;
    int nPixelHits  = 0;
    int nSctHits    = 0;
    int nStrawHits  = 0;
    int nTrtHits    = 0;
      
    double chi2 = 0;
    double dof  = 0;
 
    bool expectBL = false;
 
    nSctHits   += 0;
    nPixelHits += 0;

    int nSiHits     = 0;
 
    unsigned long id = (unsigned long)track;
 
    unsigned bitmap = 0;
 
    int trackAuthor = track->pdgId();
    int uniqueID     = HepMC::uniqueID(track);
 
#if 0
    std::cout << "\t\t\tSUTT TP track" 
          << "\teta=" << eta  
          << "\tphi=" << phi  
          << "\tz0="  << z0 
          << "\tpT="  << pT 
          << "\td0="  << d0
          << "\tauthor=" << trackAuthor
          << "\tVTX x " << xp[0]<< "\ty " << xp[1] << "\tz " << xp[2] 
          << std::endl;
#endif
 
    // Create and save Track
      
    TIDA::Track* t = new TIDA::Track( eta,   phi,  z0,  d0,  pT, chi2, dof,
                      deta,  dphi, dz0, dd0, dpT,
                      nBlayerHits, nPixelHits, nSctHits, nSiHits,
                      nStrawHits,  nTrtHits,   bitmap, 0,
                      trackAuthor,  false, uniqueID, -1,  
                      expectBL, id) ;  

    //      std::cout << "SUTT TP track " << *t << "\t0x" << std::hex << bitmap << std::dec << std::endl; 
      
    // addTrack applies additional cuts using the Filter
    if ( !addTrack( t ) ){
      delete t;
      return false;
    }
  }
  return false;
     
}

selectTrack() [8/8]

bool TrigTrackSelector::selectTrack

(

HepMC::ConstGenParticlePtr

track

)

not a "final state" particle

set this so can use it as the identifier - don't forget to reset!!

Definition at line 500 of file TrigTrackSelector.cxx.

                                                                    {
  
    if ( !MC::isStable(track) ) return false;

//AV Using memory to get some value is not a good idea. This is not a repruducible/portable way, but I leave it as is.
#ifdef HEPMC3
    m_id = (unsigned long)(track.get());
#else
    m_id = (unsigned long)track;
#endif
    bool sel;
    sel = selectTrack( TruthParticle(track) );
    m_id = 0;
    
    return sel; 
    
}

selectTracks() [1/8]

void TrigTrackSelector::selectTracks

(

const Rec::TrackParticleContainer *

trigtracks

)

Definition at line 340 of file TrigTrackSelector.cxx.

                                                                                  { 
    
    //    std::cout << "\t\t\tSUTT \tTrackParticleContainer->size() = " << trigtracks->size() << std::endl;
    
    Rec::TrackParticleContainer::const_iterator  trackitr = trigtracks->begin();
    Rec::TrackParticleContainer::const_iterator  trackend = trigtracks->end();
   
    while ( trackitr!=trackend ) { 
    
      selectTrack( *trackitr );
      
      ++trackitr;
 
    } // loop over tracks
    
}

selectTracks() [2/8]

void TrigTrackSelector::selectTracks

(

const TrackCollection *

trigtracks

)

Definition at line 1043 of file TrigTrackSelector.cxx.

                                                                         { 
    
    //    std::cout << "\t\t\tSUTT \tTrackContainer->size() = " << trigtracks->size() << std::endl;
    
    TrackCollection::const_iterator  trackitr = trigtracks->begin();
    TrackCollection::const_iterator  trackend = trigtracks->end();
   
    while ( trackitr!=trackend ) { 
      selectTrack( *trackitr );
      ++trackitr;
    } // loop over tracks
    
}

selectTracks() [3/8]

void TrigTrackSelector::selectTracks

(

const TrackParticleLinks_t &

tracks

)

Definition at line 1268 of file TrigTrackSelector.cxx.

                                                                         {
  for( const auto& track : tracks ) selectTrack( *track );
}

selectTracks() [4/8]

void TrigTrackSelector::selectTracks	(	const TrigInDetTrackCollection *	trigtracks,
		const TrigInDetTrackTruthMap *	truthMap = 0 )

Definition at line 188 of file TrigTrackSelector.cxx.

                                                                                                                         {     
    // do the track extraction stuff here....
    TrigInDetTrackCollection::const_iterator  trackitr = trigtracks->begin();
    TrigInDetTrackCollection::const_iterator  trackend = trigtracks->end();
    while ( trackitr!=trackend ) { 
      selectTrack( *trackitr, truthMap );
      ++trackitr;
    }
}

selectTracks() [5/8]

void TrigTrackSelector::selectTracks

(

const TruthParticleContainer *

truthtracks

)

Definition at line 360 of file TrigTrackSelector.cxx.

                                                                                { 
    //    std::cout << "\t\t\tSUTT \tTrackParticleContainer->size() = " << trigtracks->size() << std::endl;
    
    TruthParticleContainer::const_iterator  trackitr = truthtracks->begin();
    TruthParticleContainer::const_iterator  trackend = truthtracks->end();
   
    while ( trackitr!=trackend ) { 
    
      selectTrack( *trackitr );
      
      ++trackitr;
 
    } // loop over tracks
    
}

selectTracks() [6/8]

void TrigTrackSelector::selectTracks	(	const xAOD::TrackParticleContainer *	tracks,
		void *	= 0 )

Definition at line 1244 of file TrigTrackSelector.cxx.

                                                                                      {
    //    std::cout << "\t\t\tSUTT \tTrackContainer->size() = " << trigtracks->size() << std::endl;
    xAOD::TrackParticleContainer::const_iterator  trackitr = tracks->begin();
    xAOD::TrackParticleContainer::const_iterator  trackend = tracks->end();
    while ( trackitr!=trackend ) { 
      selectTrack( *trackitr );
      ++trackitr;
    } // loop over tracks     
}

selectTracks() [7/8]

void TrigTrackSelector::selectTracks

(

const xAOD::TruthParticleContainer *

truthtracks

)

get truth beamline

fix default (unset) TruthParticle charge

not just yet - save for later ... else ( gotPdgId && m_selectParentPdgIds.size()!=0 ) gotParentPdgId = fromAncestor(m_selectParentPdgIds, (*trackitr))!=nullptr;

Definition at line 455 of file TrigTrackSelector.cxx.

                                                                                   { 
 
  xAOD::TruthParticleContainer::const_iterator  trackitr = truthtracks->begin();
  xAOD::TruthParticleContainer::const_iterator  trackend = truthtracks->end();
  
  double x0 = 0;
  double y0 = 0;
 
  truthBeamline( truthtracks, x0, y0 );
 
 
  for ( ; trackitr!=trackend; ++trackitr) {
 
 
    // Only select charged final state particles
    double q = (*trackitr)->charge();

    static const particleType ptype;
    if ( q==-999 ) q = ptype.charge( (*trackitr)->pdgId() );
 
    if (q == 0 || !MC::isStable(*trackitr) ) continue;
 
    // If looking for tau parents, don't select mu or e children
 
    // select based on the pdg of final state particle 
    bool gotPdgId = true;
    if (m_selectPdgId!=0) gotPdgId = (*trackitr)->absPdgId()==m_selectPdgId;
 
    // select based on the pdg of the parent or ancestor
    bool gotParentPdgId = true;
    if   ( gotPdgId && m_selectParentPdgId!=0 )         gotParentPdgId = fromAncestor(m_selectParentPdgId, (*trackitr))!=nullptr;
 
    if ( gotParentPdgId && gotPdgId ) selectTrack( *trackitr, x0, y0);
 
  } // loop over tracks
    
}

selectTracks() [8/8]

void TrigTrackSelector::selectTracks	(	xAOD::TrackParticleContainer::const_iterator	trackitr,
		xAOD::TrackParticleContainer::const_iterator	trackend,
		void *	= 0 )

MT selector access.

will need this printout during debugging, so leave commented until all the feature access has been properly debugged

Definition at line 1256 of file TrigTrackSelector.cxx.

                                                                                 {
    //    std::cout << "\t\t\tSUTT \tTrackContainer->size() = " << trigtracks->size() << std::endl;
    while ( trackitr!=trackend ) { 
      selectTrack( *trackitr );
      ++trackitr;
    } // loop over tracks     
}

setBeamline()

void TrigTrackSelector::setBeamline ( double x, double y, double z = 0 )

inline

Definition at line 94 of file TrigTrackSelector.h.

{ m_xBeam = x; m_yBeam = y; m_zBeam=z; }

size()

unsigned TrackSelector::size ( ) const

inlineinherited

Definition at line 63 of file Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h.

{ return m_tracks.size(); } 

tracks() [1/2]

const std::vector< TIDA::Track * > & TrackSelector::tracks ( ) const

inlineinherited

Definition at line 53 of file Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h.

{ return m_tracks; }

tracks() [2/2]

std::vector< TIDA::Track * > TrackSelector::tracks ( TrackFilter * selector ) const

inlineinherited

Definition at line 55 of file Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h.

                                                              {
    if ( selector==0 ) return m_tracks; 
    std::vector<TIDA::Track*> t;
    for ( int i=m_tracks.size() ; i-- ; ) if ( selector->select(m_tracks[i]) ) t.push_back(m_tracks[i]);
    return t;
  }

truthBeamline()

void TrigTrackSelector::truthBeamline	(	const xAOD::TruthParticleContainer *	truthtracks,
		double &	x0,
		double &	y0 )

extract all the tracks from a xAOD::TruthParticle collection and histogram the x and y production coordinates to determine the event "beamline" position

histograms and book keeoing

positions

numbers of tracks

fill histograms ...

get track production vertex

add to the correct bins

calculate the most populous bin in x and y

require more than 1 particle for it to be classed a "vertex" therefore, for single particle Monte Carlo, this will not be updated and (0,0) will still be correctly used

Definition at line 381 of file TrigTrackSelector.cxx.

                                                                                                             { 
 
  x0 = 0;
  y0 = 0;
    
 
  int Nx = 300;
  int Ny = 300;

  std::vector<double> xpos(Nx,0);
  std::vector<double> ypos(Ny,0);

  std::vector<int> xn(Nx,0);
  std::vector<int> yn(Ny,0);
 
  int xoffset = Nx/2;; 
  int yoffset = Ny/2; 
 
  double deltax = 3.0/Nx;
  double deltay = 3.0/Ny;

 
  xAOD::TruthParticleContainer::const_iterator  trackitr = truthtracks->begin();
  xAOD::TruthParticleContainer::const_iterator  trackend = truthtracks->end();
    
  for ( ; trackitr!=trackend ; ++trackitr ) { 
    
    const xAOD::TruthParticle* track = (*trackitr);
    
    if ( !MC::isStable(track) || !track->hasProdVtx() ) continue; 

 
    double xp[3] = { track->prodVtx()->x(), track->prodVtx()->y(), track->prodVtx()->z() };

 
    int ix = xp[0]/deltax + xoffset;
    int iy = xp[1]/deltay + yoffset;
 
    if ( ix<0 || ix>=Nx || iy<0 || iy>=Nx ) continue;
 
    xpos[ix] += xp[0];
    ypos[iy] += xp[1];
 
    xn[ix]++;
    yn[iy]++;
 
  } // loop over tracks

 
  int imx = 0; 
  int imy = 0; 
 
  for ( size_t i=0 ; i<xpos.size() ; i++ ) {
    if ( xn[i]>xn[imx] ) imx = i;
    if ( yn[i]>yn[imy] ) imy = i;
  }

  if ( xn[imx]>1 ) x0 = xpos[imx]/xn[imx];
  if ( yn[imy]>1 ) y0 = ypos[imy]/yn[imy];
 
}

Member Data Documentation

m_correctTrkTracks

bool TrigTrackSelector::m_correctTrkTracks

private

Definition at line 198 of file TrigTrackSelector.h.

m_id

unsigned long TrigTrackSelector::m_id

private

Definition at line 192 of file TrigTrackSelector.h.

m_radius

double TrigTrackSelector::m_radius

private

Definition at line 200 of file TrigTrackSelector.h.

m_selector

TrackFilter* TrackSelector::m_selector

protectedinherited

Definition at line 85 of file Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h.

m_selectParentPdgId

int TrigTrackSelector::m_selectParentPdgId

private

Definition at line 203 of file TrigTrackSelector.h.

m_selectPdgId

int TrigTrackSelector::m_selectPdgId

private

Definition at line 202 of file TrigTrackSelector.h.

m_tracks

std::vector<TIDA::Track*> TrackSelector::m_tracks

protectedinherited

Definition at line 82 of file Trigger/TrigAnalysis/TrigInDetAnalysis/TrigInDetAnalysis/TrackSelector.h.

m_xBeam

double TrigTrackSelector::m_xBeam

private

Definition at line 194 of file TrigTrackSelector.h.

m_yBeam

double TrigTrackSelector::m_yBeam

private

Definition at line 195 of file TrigTrackSelector.h.

m_zBeam

double TrigTrackSelector::m_zBeam

private

Definition at line 196 of file TrigTrackSelector.h.

s_default_radius

const double TrigTrackSelector::s_default_radius = 47

staticprivate

NB: This was 47 for Run 2, but with the addition of the IBL it should be 32 It was kept at 47 for all Run 2 and migration to MT, but for Run 3 we really want it changed to be 32.

Definition at line 205 of file TrigTrackSelector.h.

---

The documentation for this class was generated from the following files:

• TrigTrackSelector.h
• TrigTrackSelector.cxx