| 1 |
bendavid |
1.1 |
//--------------------------------------------------------------------------------------------------
|
| 2 |
loizides |
1.6 |
// $Id: HitDropper.h,v 1.5 2009/07/09 17:15:57 bendavid Exp $
|
| 3 |
bendavid |
1.1 |
//
|
| 4 |
|
|
// HitDropper
|
| 5 |
|
|
//
|
| 6 |
loizides |
1.4 |
// Utility to remove (or rather flag) hits from a track which are inconsistant with a vertex fit.
|
| 7 |
bendavid |
1.1 |
//
|
| 8 |
|
|
// Authors: J.Bendavid
|
| 9 |
|
|
//--------------------------------------------------------------------------------------------------
|
| 10 |
|
|
|
| 11 |
|
|
#ifndef MITEDM_PRODUCERS_HITDROPPER_H
|
| 12 |
|
|
#define MITEDM_PRODUCERS_HITDROPPER_H
|
| 13 |
|
|
|
| 14 |
|
|
#include "TrackingTools/TransientTrack/interface/TransientTrackBuilder.h"
|
| 15 |
|
|
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
|
| 16 |
|
|
#include "DataFormats/TrackReco/interface/Track.h"
|
| 17 |
|
|
#include "TrackingTools/TransientTrack/interface/TransientTrack.h"
|
| 18 |
bendavid |
1.2 |
#include "RecoTracker/TkDetLayers/interface/GeometricSearchTracker.h"
|
| 19 |
|
|
#include "Geometry/CommonDetUnit/interface/GeomDet.h"
|
| 20 |
bendavid |
1.1 |
#include "MitEdm/DataFormats/interface/Types.h"
|
| 21 |
|
|
|
| 22 |
|
|
namespace mitedm
|
| 23 |
|
|
{
|
| 24 |
|
|
class HitDropper
|
| 25 |
|
|
{
|
| 26 |
|
|
public:
|
| 27 |
loizides |
1.6 |
HitDropper(const TrackerGeometry *geo, const TransientTrackBuilder *builder,
|
| 28 |
|
|
const GeometricSearchTracker *search) :
|
| 29 |
|
|
trackerGeo_(geo),
|
| 30 |
|
|
builder_(builder),
|
| 31 |
|
|
trackerGeoSearch_(search) {}
|
| 32 |
bendavid |
1.1 |
~HitDropper() {}
|
| 33 |
|
|
|
| 34 |
loizides |
1.4 |
reco::HitPattern CorrectedHits(const reco::TransientTrack *tTrack,
|
| 35 |
|
|
const ThreeVector &vtxPos) const;
|
| 36 |
|
|
reco::HitPattern CorrectedHits(const reco::Track *track,
|
| 37 |
|
|
const ThreeVector &vtxPos) const;
|
| 38 |
|
|
reco::HitPattern CorrectedHits(const reco::Track *track,
|
| 39 |
|
|
const ThreeVector &vtxPos,
|
| 40 |
|
|
const ThreeVector &trkMom,
|
| 41 |
|
|
Double_t lxyError,
|
| 42 |
|
|
Double_t lzError,
|
| 43 |
|
|
Double_t sigmaTolerance=3.0) const;
|
| 44 |
bendavid |
1.5 |
reco::HitPattern CorrectedHitsAOD(const reco::Track *track,
|
| 45 |
|
|
const ThreeVector &vtxPos,
|
| 46 |
|
|
const ThreeVector &trkMom,
|
| 47 |
|
|
Double_t lxyError,
|
| 48 |
|
|
Double_t lzError,
|
| 49 |
|
|
Double_t sigmaTolerance=3.0) const;
|
| 50 |
loizides |
1.6 |
bool IsBarrel(const GeomDet *det) const;
|
| 51 |
|
|
bool IsDisk(const GeomDet *det) const;
|
| 52 |
bendavid |
1.1 |
|
| 53 |
|
|
protected:
|
| 54 |
loizides |
1.6 |
const TrackerGeometry *trackerGeo_; //tracker geometry
|
| 55 |
|
|
const TransientTrackBuilder *builder_; //track builder
|
| 56 |
|
|
const GeometricSearchTracker *trackerGeoSearch_; //geometrical track searcher
|
| 57 |
bendavid |
1.1 |
};
|
| 58 |
|
|
}
|
| 59 |
bendavid |
1.2 |
|
| 60 |
|
|
//-------------------------------------------------------------------------------------------------
|
| 61 |
loizides |
1.6 |
inline bool mitedm::HitDropper::IsBarrel(const GeomDet *det) const
|
| 62 |
bendavid |
1.2 |
{
|
| 63 |
loizides |
1.4 |
// Check if detector element is in a barrel layer.
|
| 64 |
|
|
|
| 65 |
bendavid |
1.2 |
GeomDetEnumerators::SubDetector subDet = det->subDetector();
|
| 66 |
|
|
|
| 67 |
|
|
if (subDet==GeomDetEnumerators::PixelBarrel || subDet==GeomDetEnumerators::TIB
|
| 68 |
|
|
|| subDet==GeomDetEnumerators::TOB || subDet==GeomDetEnumerators::DT
|
| 69 |
|
|
|| subDet==GeomDetEnumerators::RPCBarrel)
|
| 70 |
|
|
return true;
|
| 71 |
|
|
else
|
| 72 |
|
|
return false;
|
| 73 |
|
|
}
|
| 74 |
|
|
|
| 75 |
|
|
//-------------------------------------------------------------------------------------------------
|
| 76 |
loizides |
1.6 |
inline bool mitedm::HitDropper::IsDisk(const GeomDet *det) const
|
| 77 |
bendavid |
1.2 |
{
|
| 78 |
loizides |
1.4 |
// Check if detector element is in a barrel layer.
|
| 79 |
|
|
|
| 80 |
bendavid |
1.2 |
GeomDetEnumerators::SubDetector subDet = det->subDetector();
|
| 81 |
|
|
|
| 82 |
|
|
if (subDet==GeomDetEnumerators::PixelEndcap || subDet==GeomDetEnumerators::TID
|
| 83 |
|
|
|| subDet==GeomDetEnumerators::TEC || subDet==GeomDetEnumerators::CSC
|
| 84 |
|
|
|| subDet==GeomDetEnumerators::RPCEndcap)
|
| 85 |
|
|
return true;
|
| 86 |
|
|
else
|
| 87 |
|
|
return false;
|
| 88 |
|
|
}
|
| 89 |
bendavid |
1.1 |
#endif
|
