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