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 DisplacedVertexTag_(iConfig.getParameter<
edm::InputTag>(
"DisplacedVertexTag")),
46 beamSpotTag_ (iConfig.getParameter<
edm::InputTag> (
"BeamSpotTag")),
48 MuonTag_ (iConfig.getParameter<
edm::InputTag>(
"MuonTag")),
59 desc.
add<
bool>(
"FastAccept",
false);
60 desc.
add<
double>(
"MinLxySignificance",0.0);
61 desc.
add<
double>(
"MaxLxySignificance",0.0);
62 desc.
add<
double>(
"MaxNormalisedChi2",10.0);
63 desc.
add<
double>(
"MinVtxProbability",0.0);
64 desc.
add<
double>(
"MinCosinePointingAngle",-2.0);
68 descriptions.
add(
"hltDisplacedmumuFilter", desc);
92 vertexBeamSpot = *recoBeamSpotHandle;
99 if(foundVertexColl) displacedVertexColl = *displacedVertexCollHandle;
117 bool triggered =
false;
120 for(
auto displacedVertex : displacedVertexColl){
122 if(displacedVertex.tracksSize() != 2)
throw cms::Exception(
"BadLogic") <<
"HLTDisplacedmumuFilter: ERROR: the Jpsi vertex must have exactly two muons by definition. It now has n muons = " 123 << displacedVertex.tracksSize() << std::endl;
125 float normChi2 = displacedVertex.normalizedChi2();
128 double vtxProb = 0.0;
129 if( (displacedVertex.chi2()>=0.0) && (displacedVertex.ndof()>0) ) vtxProb = TMath::Prob(displacedVertex.chi2(), displacedVertex.ndof() );
133 auto trackIt = displacedVertex.tracks_begin();
140 reco::RecoChargedCandidateCollection::const_iterator cand1;
141 reco::RecoChargedCandidateCollection::const_iterator cand2;
144 for (
auto cand=mucands->begin();
cand!=mucands->end();
cand++) {
146 if(tkRef == vertextkRef1) {cand1 =
cand; iFoundRefs++;}
147 if(tkRef == vertextkRef2) {cand2 =
cand; iFoundRefs++;}
149 if(iFoundRefs < 2)
throw cms::Exception(
"BadLogic") <<
"HLTDisplacedmumuFilter: ERROR: the muons matched with the Jpsi vertex tracks should be at least two by definition." << std::endl;
153 cand1->py() + cand2->py(),
159 GlobalPoint secondaryVertex (vpoint.x(), vpoint.y(), vpoint.z());
163 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) );
165 GlobalPoint displacementFromBeamspot( -1*((vertexBeamSpot.
x0() - secondaryVertex.x()) + (secondaryVertex.z() - vertexBeamSpot.
z0()) * vertexBeamSpot.
dxdz()),
166 -1*((vertexBeamSpot.
y0() - secondaryVertex.y())+ (secondaryVertex.z() - vertexBeamSpot.
z0()) * vertexBeamSpot.
dydz()), 0);
168 float lxy = displacementFromBeamspot.perp();
169 float lxyerr =
sqrt(err.rerr(displacementFromBeamspot));
175 float cosAlpha = vperp.Dot(pperp)/(vperp.R()*pperp.R());
192 LogDebug(
"HLTDisplacedMumuFilter") <<
" >>>>> Result of HLTDisplacedMumuFilter is "<< triggered <<std::endl;
HLTDisplacedmumuFilter(const edm::ParameterSet &)
double z0() const
z coordinate
double minVtxProbability_
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)
std::vector< Vertex > VertexCollection
collection of Vertex objects
double maxLxySignificance_
void addObject(int id, const reco::RecoEcalCandidateRef &ref)
setters for L3 collections: (id=physics type, and Ref<C>)
double minCosinePointingAngle_
double dydz() const
dydz slope
~HLTDisplacedmumuFilter() override
math::XYZPoint Point
point in the space
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
double maxNormalisedChi2_
ParameterDescriptionBase * add(U const &iLabel, T const &value)
double dxdz() const
dxdz slope
double minLxySignificance_
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
std::vector< RecoChargedCandidate > RecoChargedCandidateCollection
collectin of RecoChargedCandidate objects
static void makeHLTFilterDescription(edm::ParameterSetDescription &desc)
XYZVectorD XYZVector
spatial vector with cartesian internal representation
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > MuonToken_
void addCollectionTag(const edm::InputTag &collectionTag)
collectionTags
void add(std::string const &label, ParameterSetDescription const &psetDescription)
double y0() const
y coordinate
bool hltFilter(edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct) const override
edm::EDGetTokenT< reco::VertexCollection > DisplacedVertexToken_
double x0() const
x coordinate