98 const double MuMass(0.106);
99 const double MuMass2(MuMass*MuMass);
144 vector<TrackRef> trkMuCands;
147 vector<bool> isUsedCand(trkcands->size(),
false);
152 for (
auto mucand1=mucands->begin(), endCand1=mucands->end(); mucand1!=endCand1; ++mucand1) {
154 if ( mucands->size()<2)
break;
155 if ( trkcands->empty())
break;
158 LogDebug(
"HLTDisplacedMumukFilter") <<
" 1st muon: q*pt= " << trk1->charge()*trk1->pt() <<
", eta= " << trk1->eta() <<
", hits= " << trk1->numberOfValidHits();
161 if (fabs(trk1->eta()) >
maxEta_)
continue;
164 if (trk1->pt() <
minPt_)
continue;
166 auto mucand2 = mucand1; ++mucand2;
168 for (
auto endCand2=mucands->end(); mucand2!=endCand2; ++mucand2) {
172 LogDebug(
"HLTDisplacedMumukFilter") <<
" 2nd muon: q*pt= " << trk2->charge()*trk2->pt() <<
", eta= " << trk2->eta() <<
", hits= " << trk2->numberOfValidHits();
175 if (fabs(trk2->eta()) >
maxEta_)
continue;
178 if (trk2->pt() <
minPt_)
continue;
180 RecoChargedCandidateCollection::const_iterator trkcand, endCandTrk;
182 std::vector<bool>::iterator isUsedIter, endIsUsedCand;
185 for ( trkcand = trkcands->begin(), endCandTrk=trkcands->end(), isUsedIter = isUsedCand.begin(), endIsUsedCand = isUsedCand.end(); trkcand != endCandTrk && isUsedIter != endIsUsedCand; ++trkcand, ++isUsedIter) {
189 if (
overlap(*mucand1,*trkcand)) {
190 trkMuCands.push_back(trk3);
194 else if (
overlap(*mucand2,*trkcand)){
195 trkMuCands.push_back(trk3);
200 if(trkMuCands.size()==2)
break;
207 for ( trkcand = trkcands->begin(), endCandTrk=trkcands->end(), isUsedIter = isUsedCand.begin(), endIsUsedCand = isUsedCand.end(); trkcand != endCandTrk && isUsedIter != endIsUsedCand; ++trkcand, ++isUsedIter) {
211 LogDebug(
"HLTDisplacedMumukFilter") <<
" 3rd track: q*pt= " << trk3->charge()*trk3->pt() <<
", eta= " << trk3->eta() <<
", hits= " << trk3->numberOfValidHits();
215 for (
auto & trkMuCand : trkMuCands)
if(trk3==trkMuCand) skip=
true;
219 if(*isUsedIter)
continue;
222 if (fabs(trk3->eta()) >
maxEta_)
continue;
225 if (trk3->pt() <
minPt_)
continue;
228 e1 =
sqrt(trk1->momentum().Mag2()+MuMass2);
229 e2 =
sqrt(trk2->momentum().Mag2()+MuMass2);
230 e3 =
sqrt(trk3->momentum().Mag2()+thirdTrackMass2);
239 double invmass =
abs(p.mass());
241 LogDebug(
"HLTDisplacedMumukFilter") <<
" Invmass= " << invmass;
246 vector<TransientTrack> t_tks;
247 t_tks.push_back((*theB).build(&trk1));
248 t_tks.push_back((*theB).build(&trk2));
249 t_tks.push_back((*theB).build(&trk3));
251 if (t_tks.size()!=3)
continue;
255 double d0sig = tscb.transverseImpactParameter().significance();
271 GlobalPoint displacementFromBeamspot( -1*((vertexBeamSpot.
x0() -secondaryVertex.
x()) + (secondaryVertex.
z() - vertexBeamSpot.
z0()) * vertexBeamSpot.
dxdz()), -1*((vertexBeamSpot.
y0() - secondaryVertex.
y())+ (secondaryVertex.
z() -vertexBeamSpot.
z0()) * vertexBeamSpot.
dydz()), 0);
273 float lxy = displacementFromBeamspot.perp();
274 float lxyerr =
sqrt(err.
rerr(displacementFromBeamspot));
280 Vertex::Point vperp(displacementFromBeamspot.x(),displacementFromBeamspot.y(),0.);
283 float cosAlpha = vperp.Dot(pperp)/(vperp.R()*pperp.R());
285 LogDebug(
"HLTDisplacedMumukFilter") <<
" vertex fit normalised chi2: " << normChi2 <<
", Lxy significance: " << lxy/lxyerr <<
", cosine pointing angle: " << cosAlpha;
294 LogDebug(
"HLTDisplacedMumukFilter") <<
" Event passed!";
303 vector<RecoChargedCandidateRef> vref;
305 for (
auto &
i : vref) {
310 }
else if (trktmp==trk2) {
312 }
else if (trktmp==trk3) {
315 if (i1done && i2done && i3done)
break;
339 const bool accept (counter >= 1);
341 LogDebug(
"HLTDisplacedMumukFilter") <<
" >>>>> Result of HLTDisplacedMumukFilter is "<<
accept <<
", number of muon pairs passing thresholds= " <<
counter;
GlobalError positionError() const
static int overlap(const reco::Candidate &, const reco::Candidate &)
double z0() const
z coordinate
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > trkCandToken_
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
void getObjects(Vids &ids, VRphoton &refs) const
various physics-level getters:
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::Ref< RecoChargedCandidateCollection > RecoChargedCandidateRef
reference to an object in a collection of RecoChargedCandidate objects
std::vector< Vertex > VertexCollection
collection of Vertex objects
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
void addObject(int id, const reco::RecoEcalCandidateRef &ref)
setters for L3 collections: (id=physics type, and Ref<C>)
double dydz() const
dydz slope
float normalisedChiSquared() const
GlobalPoint position() const
Abs< T >::type abs(const T &t)
CachingVertex< 5 > vertex(const std::vector< reco::TransientTrack > &tracks) const override
math::XYZPoint Point
point in the space
edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
T const * get() const
Returns C++ pointer to the item.
edm::InputTag trkCandTag_
double dxdz() const
dxdz slope
edm::EDGetTokenT< reco::RecoChargedCandidateCollection > muCandToken_
XYZVectorD XYZVector
spatial vector with cartesian internal representation
T rerr(const GlobalPoint &aPoint) const
void addCollectionTag(const edm::InputTag &collectionTag)
collectionTags
const double minCosinePointingAngle_
const double thirdTrackMass_
static std::atomic< unsigned int > counter
double y0() const
y coordinate
const double maxNormalisedChi2_
static FreeTrajectoryState initialFreeState(const reco::Track &, const MagneticField *)
edm::OwnVector< Candidate > CandidateCollection
collection of Candidate objects
const double minD0Significance_
T const * product() const
math::PtEtaPhiELorentzVectorF LorentzVector
const double minLxySignificance_
double x0() const
x coordinate