66 double const MuMass = 0.106;
67 double const MuMass2 = MuMass * MuMass;
83 RecoChargedCandidateCollection::const_iterator cand1;
84 RecoChargedCandidateCollection::const_iterator cand2;
85 RecoChargedCandidateCollection::const_iterator cand3;
91 vector<RecoChargedCandidateRef> vPrevCands;
94 for (cand1 = mucands->begin(); cand1 != mucands->end(); cand1++) {
96 LogDebug(
"HLTDisplacedMumumuFilter") <<
" 1st muon in loop: q*pt= " << cand1->charge() * cand1->pt()
97 <<
", eta= " << cand1->eta() <<
", hits= " << tk1->numberOfValidHits();
104 if (fabs(cand1->eta()) >
maxEta_)
111 for (; cand2 != mucands->end(); cand2++) {
115 LogDebug(
"HLTMuonDimuonFilter") <<
" 2nd muon in loop: q*pt= " << cand2->charge() * cand2->pt()
116 <<
", eta= " << cand2->eta() <<
", hits= " << tk2->numberOfValidHits()
117 <<
", d0= " << tk2->d0();
123 if (fabs(cand2->eta()) >
maxEta_)
130 for (; cand3 != mucands->end(); cand3++) {
134 LogDebug(
"HLTMuonDimuonFilter") <<
" 3rd muon in loop: q*pt= " << cand3->charge() * cand3->pt()
135 <<
", eta= " << cand3->eta() <<
", hits= " << tk3->numberOfValidHits()
136 <<
", d0= " << tk3->d0();
142 if (fabs(cand3->eta()) >
maxEta_)
149 if (fabs(cand1->charge() + cand2->charge() + cand3->charge()) !=
chargeOpt_)
154 e1 =
sqrt(cand1->momentum().Mag2() + MuMass2);
155 e2 =
sqrt(cand2->momentum().Mag2() + MuMass2);
156 e3 =
sqrt(cand3->momentum().Mag2() + MuMass2);
165 double invmass =
abs(p.mass());
166 LogDebug(
"HLTDisplacedMumumuFilter") <<
" ... 1-2 invmass= " << invmass;
174 vector<TransientTrack> t_tks;
178 t_tks.push_back(ttkp1);
179 t_tks.push_back(ttkp2);
180 t_tks.push_back(ttkp3);
182 if (t_tks.size() != 3)
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
const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > previousCandToken_
std::vector< Vertex > VertexCollection
collection of Vertex objects
const edm::EDGetTokenT< reco::RecoChargedCandidateCollection > srcToken_
Abs< T >::type abs(const T &t)
CachingVertex< 5 > vertex(const std::vector< reco::TransientTrack > &tracks) const override
const double minPtTriplet_
T const * get() const
Returns C++ pointer to the item.
bool checkPreviousCand(const reco::TrackRef &trackref, const std::vector< reco::RecoChargedCandidateRef > &ref2) const
math::PtEtaPhiELorentzVectorF LorentzVector