82 const double MuMass(0.106);
83 const double MuMass2(MuMass * MuMass);
117 double e1, e2, e3_m3, e3_m4, e4_m3, e4_m4;
120 if (mucands->size() < 2)
122 if (trkcands->size() < 2)
125 RecoChargedCandidateCollection::const_iterator mucand1;
126 RecoChargedCandidateCollection::const_iterator mucand2;
127 RecoChargedCandidateCollection::const_iterator trkcand1;
128 RecoChargedCandidateCollection::const_iterator trkcand2;
134 vector<RecoChargedCandidateRef> vPrevCands;
137 for (mucand1 = mucands->begin(); mucand1 != mucands->end(); ++mucand1) {
139 LogDebug(
"HLTmumutktkVtxProducer") <<
" 1st muon: q*pt= " << trk1->charge() * trk1->pt() <<
", eta= " << trk1->eta()
140 <<
", hits= " << trk1->numberOfValidHits();
146 if (fabs(trk1->eta()) >
maxEta_)
153 for (; mucand2 != mucands->end(); mucand2++) {
158 LogDebug(
"HLTDisplacedMumukFilter") <<
" 2nd muon: q*pt= " << trk2->charge() * trk2->pt()
159 <<
", eta= " << trk2->eta() <<
", hits= " << trk2->numberOfValidHits();
165 if (fabs(trk2->eta()) >
maxEta_)
171 for (trkcand1 = trkcands->begin(); trkcand1 != trkcands->end(); ++trkcand1) {
179 LogDebug(
"HLTDisplacedMumukFilter") <<
" 3rd track: q*pt= " << trk3->charge() * trk3->pt()
180 <<
", eta= " << trk3->eta() <<
", hits= " << trk3->numberOfValidHits();
183 if (fabs(trk3->eta()) >
maxEta_)
190 double d0sigTrk3 = tscb_Trk3.transverseImpactParameter().significance();
195 for (trkcand2 = trkcands->begin(); trkcand2 != trkcands->end(); ++trkcand2) {
199 if (trk3->charge() * trk4->charge() != -1)
209 LogDebug(
"HLTDisplacedMumukFilter") <<
" 4th track: q*pt= " << trk4->charge() * trk4->pt()
210 <<
", eta= " << trk4->eta() <<
", hits= " << trk4->numberOfValidHits();
213 if (fabs(trk4->eta()) >
maxEta_)
220 double d0sigTrk4 = tscb_Trk4.transverseImpactParameter().significance();
225 e1 =
sqrt(trk1->momentum().Mag2() + MuMass2);
226 e2 =
sqrt(trk2->momentum().Mag2() + MuMass2);
227 e3_m3 =
sqrt(trk3->momentum().Mag2() + thirdTrackMass2);
228 e3_m4 =
sqrt(trk3->momentum().Mag2() + fourthTrackMass2);
229 e4_m3 =
sqrt(trk4->momentum().Mag2() + thirdTrackMass2);
230 e4_m4 =
sqrt(trk4->momentum().Mag2() + fourthTrackMass2);
239 p = p1 + p2 + p3_m3 + p4_m4;
240 pBar = p1 + p2 + p3_m4 + p4_m3;
241 p_m3m4 = p3_m3 + p4_m4;
242 p_m4m3 = p3_m4 + p4_m3;
253 vector<TransientTrack> t_tks;
254 t_tks.push_back((*theB).build(&trk1));
255 t_tks.push_back((*theB).build(&trk2));
256 t_tks.push_back((*theB).build(&trk3));
257 t_tks.push_back((*theB).build(&trk4));
258 if (t_tks.size() != 4)
bool checkPreviousCand(const reco::TrackRef &trackref, const std::vector< reco::RecoChargedCandidateRef > &ref2) const
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
bool overlap(const reco::TrackRef &trackref1, const reco::TrackRef &trackref2)
std::vector< Vertex > VertexCollection
collection of Vertex objects
static FreeTrajectoryState initialFreeState(const reco::Track &, const MagneticField *)
const edm::EDGetTokenT< reco::RecoChargedCandidateCollection > trkCandToken_
const edm::EDGetTokenT< reco::RecoChargedCandidateCollection > muCandToken_
const edm::EDGetTokenT< reco::BeamSpot > beamSpotToken_
CachingVertex< 5 > vertex(const std::vector< reco::TransientTrack > &tracks) const override
T const * get() const
Returns C++ pointer to the item.
const std::string mfName_
const double minTrkTrkMass_
const double fourthTrackMass_
const double thirdTrackMass_
const edm::EDGetTokenT< trigger::TriggerFilterObjectWithRefs > previousCandToken_
const double minD0Significance_
T const * product() const
const double maxTrkTrkMass_
math::PtEtaPhiELorentzVectorF LorentzVector