38 fastAccept_ (iConfig.getParameter<
bool>(
"FastAccept")),
39 minLxySignificance_ (iConfig.getParameter<double>(
"MinLxySignificance")),
40 maxLxySignificance_ (iConfig.getParameter<double>(
"MaxLxySignificance")),
41 maxNormalisedChi2_ (iConfig.getParameter<double>(
"MaxNormalisedChi2")),
42 minVtxProbability_ (iConfig.getParameter<double>(
"MinVtxProbability")),
43 minCosinePointingAngle_ (iConfig.getParameter<double>(
"MinCosinePointingAngle")),
44 triggerTypeDaughters_(iConfig.getParameter<
int>(
"triggerTypeDaughters")),
45 DisplacedVertexTag_(iConfig.getParameter<
edm::InputTag>(
"DisplacedVertexTag")),
47 beamSpotTag_ (iConfig.getParameter<
edm::InputTag> (
"BeamSpotTag")),
49 TrackTag_ (iConfig.getParameter<
edm::InputTag>(
"TrackTag")),
60 desc.
add<
bool>(
"FastAccept",
false);
61 desc.
add<
double>(
"MinLxySignificance",0.0);
62 desc.
add<
double>(
"MaxLxySignificance",0.0);
63 desc.
add<
double>(
"MaxNormalisedChi2",10.0);
64 desc.
add<
double>(
"MinVtxProbability",0.0);
65 desc.
add<
double>(
"MinCosinePointingAngle",-2.0);
66 desc.
add<
int>(
"triggerTypeDaughters",0);
70 descriptions.
add(
"hltDisplacedtktkFilter", desc);
94 vertexBeamSpot = *recoBeamSpotHandle;
101 if(foundVertexColl) displacedVertexColl = *displacedVertexCollHandle;
119 bool triggered =
false;
122 for(
auto displacedVertex : displacedVertexColl){
124 if(displacedVertex.tracksSize() != 2){
125 edm::LogError(
"HLTDisplacedtktkFilter") <<
"HLTDisplacedtktkFilter: ERROR: the Jpsi vertex must have exactly two tracks by definition. It now has n tracks = "<< displacedVertex.tracksSize() << std::endl;
129 float normChi2 = displacedVertex.normalizedChi2();
132 double vtxProb = 0.0;
133 if( (displacedVertex.chi2()>=0.0) && (displacedVertex.ndof()>0) ) vtxProb = TMath::Prob(displacedVertex.chi2(), displacedVertex.ndof() );
137 auto trackIt = displacedVertex.tracks_begin();
144 reco::RecoChargedCandidateCollection::const_iterator cand1;
145 reco::RecoChargedCandidateCollection::const_iterator cand2;
148 for (
auto cand=trackcands->begin();
cand!=trackcands->end();
cand++) {
150 if(tkRef == vertextkRef1) {cand1 =
cand; iFoundRefs++;}
151 if(tkRef == vertextkRef2) {cand2 =
cand; iFoundRefs++;}
155 edm::LogError(
"HLTDisplacedtktkFilter") <<
"HLTDisplacedtktkFilter: ERROR: the tracks matched with the Jpsi vertex tracks should be at least two by definition." << std::endl;
160 cand1->py() + cand2->py(),
166 GlobalPoint secondaryVertex (vpoint.x(), vpoint.y(), vpoint.z());
170 GlobalError err(verr.At(0,0), verr.At(1,0), verr.At(1,1), verr.At(2,0), verr.At(2,1), verr.At(2,2) );
172 GlobalPoint displacementFromBeamspot( -1*((vertexBeamSpot.
x0() - secondaryVertex.x()) + (secondaryVertex.z() - vertexBeamSpot.
z0()) * vertexBeamSpot.
dxdz()),
173 -1*((vertexBeamSpot.
y0() - secondaryVertex.y())+ (secondaryVertex.z() - vertexBeamSpot.
z0()) * vertexBeamSpot.
dydz()), 0);
175 float lxy = displacementFromBeamspot.perp();
176 float lxyerr =
sqrt(err.rerr(displacementFromBeamspot));
182 float cosAlpha = vperp.Dot(pperp)/(vperp.R()*pperp.R());
199 LogDebug(
"HLTDisplacedtktkFilter") <<
" >>>>> Result of HLTDisplacedtktkFilter is "<< triggered <<std::endl;
double z0() const
z coordinate
double minCosinePointingAngle_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::Ref< RecoChargedCandidateCollection > RecoChargedCandidateRef
reference to an object in a collection of RecoChargedCandidate objects
math::Error< dimension >::type Error
covariance error matrix (3x3)
bool hltFilter(edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override
std::vector< Vertex > VertexCollection
collection of Vertex objects
edm::EDGetTokenT< reco::VertexCollection > DisplacedVertexToken_
~HLTDisplacedtktkFilter() override
void addObject(int id, const reco::RecoEcalCandidateRef &ref)
setters for L3 collections: (id=physics type, and Ref<C>)
const int triggerTypeDaughters_
double dydz() const
dydz slope
HLTDisplacedtktkFilter(const edm::ParameterSet &)
double minLxySignificance_
math::XYZPoint Point
point in the space
ParameterDescriptionBase * add(U const &iLabel, T const &value)
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
double dxdz() const
dxdz slope
std::vector< RecoChargedCandidate > RecoChargedCandidateCollection
collectin of RecoChargedCandidate objects
static void makeHLTFilterDescription(edm::ParameterSetDescription &desc)
XYZVectorD XYZVector
spatial vector with cartesian internal representation
void addCollectionTag(const edm::InputTag &collectionTag)
collectionTags
void add(std::string const &label, ParameterSetDescription const &psetDescription)
double y0() const
y coordinate
double minVtxProbability_
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > TrackToken_
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
double maxNormalisedChi2_
double maxLxySignificance_
double x0() const
x coordinate