336 {
338 assert(cascadeinfoContainer!=nullptr && cascadeinfoContainer_noConstr!=nullptr);
339
340
343
344 std::vector<const xAOD::TrackParticle*> tracksJpsi;
345 std::vector<const xAOD::TrackParticle*> tracksDiTrk;
346 std::vector<const xAOD::TrackParticle*> tracksPsi;
347 std::vector<const xAOD::TrackParticle*> tracksJpsi2;
348 std::vector<double> massesPsi;
354
355
358
359
362
363
364 std::vector<const xAOD::Vertex*> selectedJpsiCandidates;
365 for(auto vxcItr=jpsiContainer.cptr()->cbegin(); vxcItr!=jpsiContainer.cptr()->cend(); ++vxcItr) {
366
371 if(flagAcc.isAvailable(*vtx) && flagAcc(*vtx)) {
373 }
374 }
376
377
378 double mass_jpsi2 =
m_V0Tools->invariantMass(*vxcItr, massesJpsi2);
379 if (mass_jpsi2 < m_jpsi2MassLower || mass_jpsi2 >
m_jpsi2MassUpper)
continue;
380
381 double chi2DOF = (*vxcItr)->chiSquared()/(*vxcItr)->numberDoF();
383
384 selectedJpsiCandidates.push_back(*vxcItr);
385 }
386 if(selectedJpsiCandidates.size()==0) return StatusCode::SUCCESS;
387
388
389 std::vector<const xAOD::Vertex*> selectedPsiCandidates;
390 for(auto vxcItr=psiContainer.cptr()->cbegin(); vxcItr!=psiContainer.cptr()->cend(); ++vxcItr) {
391
396 if(flagAcc.isAvailable(*vtx) && flagAcc(*vtx)) {
398 }
399 }
401
402
403 double mass_psi =
m_V0Tools->invariantMass(*vxcItr,massesPsi);
404 if(mass_psi < m_psiMassLower || mass_psi >
m_psiMassUpper)
continue;
405
406
407 TLorentzVector p4_mu1, p4_mu2;
408 p4_mu1.SetPtEtaPhiM( (*vxcItr)->trackParticle(0)->pt(),
409 (*vxcItr)->trackParticle(0)->eta(),
411 p4_mu2.SetPtEtaPhiM( (*vxcItr)->trackParticle(1)->pt(),
412 (*vxcItr)->trackParticle(1)->eta(),
414 double mass_jpsi = (p4_mu1 + p4_mu2).M();
415 if (mass_jpsi < m_jpsiMassLower || mass_jpsi >
m_jpsiMassUpper)
continue;
416
418 TLorentzVector p4_trk1, p4_trk2;
419 p4_trk1.SetPtEtaPhiM( (*vxcItr)->trackParticle(2)->pt(),
420 (*vxcItr)->trackParticle(2)->eta(),
422 p4_trk2.SetPtEtaPhiM( (*vxcItr)->trackParticle(3)->pt(),
423 (*vxcItr)->trackParticle(3)->eta(),
425 double mass_diTrk = (p4_trk1 + p4_trk2).M();
427 }
428
429 double chi2DOF = (*vxcItr)->chiSquared()/(*vxcItr)->numberDoF();
431
432 selectedPsiCandidates.push_back(*vxcItr);
433 }
434 if(selectedPsiCandidates.size()==0) return StatusCode::SUCCESS;
435
436 std::vector<std::pair<const xAOD::Vertex*, const xAOD::Vertex*> > candidatePairs;
437 for(auto jpsiItr=selectedJpsiCandidates.cbegin(); jpsiItr!=selectedJpsiCandidates.cend(); ++jpsiItr) {
438 tracksJpsi2.clear();
439 for(
size_t i=0;
i<(*jpsiItr)->nTrackParticles();
i++) tracksJpsi2.push_back((*jpsiItr)->trackParticle(i));
440 for(auto psiItr=selectedPsiCandidates.cbegin(); psiItr!=selectedPsiCandidates.cend(); ++psiItr) {
442 for(
size_t j=0;
j<(*psiItr)->nTrackParticles();
j++) {
443 if(std::find(tracksJpsi2.cbegin(), tracksJpsi2.cend(), (*psiItr)->trackParticle(j)) != tracksJpsi2.cend()) {
skip =
true;
break; }
444 }
446 candidatePairs.push_back(std::pair<const xAOD::Vertex*, const xAOD::Vertex*>(*jpsiItr,*psiItr));
447 }
448 }
449
450 std::sort( candidatePairs.begin(), candidatePairs.end(), [](std::pair<const xAOD::Vertex*, const xAOD::Vertex*>
a, std::pair<const xAOD::Vertex*, const xAOD::Vertex*> b) { return a.first->chiSquared()/a.first->numberDoF()+a.second->chiSquared()/a.second->numberDoF() < b.first->chiSquared()/b.first->numberDoF()+b.second->chiSquared()/b.second->numberDoF(); } );
452 candidatePairs.erase(candidatePairs.begin()+
m_maxCandidates, candidatePairs.end());
453 }
454
455 for(
size_t ic=0;
ic<candidatePairs.size();
ic++) {
458
459 tracksJpsi2.clear();
461 if (tracksJpsi2.size() != 2 || massesJpsi2.size() != 2) {
462 ATH_MSG_ERROR(
"Problems with Jpsi input: number of tracks or track mass inputs is not 2!");
463 }
464 tracksPsi.clear();
466 if (tracksPsi.size() != massesPsi.size()) {
467 ATH_MSG_ERROR(
"Problems with Psi input: number of tracks or track mass inputs is not correct!");
468 }
469
470 tracksJpsi.clear();
473 tracksDiTrk.clear();
477 }
478
479 TLorentzVector p4_moth;
484 }
488 }
490
491
492 std::unique_ptr<Trk::IVKalState> state =
m_iVertexFitter->makeState(ctx);
493
494 int robustness = 0;
496
497
498 std::vector<Trk::VertexID> vrtList;
499
501
504 } else {
506 }
507 vrtList.push_back(vID1);
508
512 } else {
514 }
515 vrtList.push_back(vID2);
516
517 std::vector<const xAOD::TrackParticle*> tp; tp.clear();
518 std::vector<double> tp_masses; tp_masses.clear();
521 std::vector<Trk::VertexID> cnstV; cnstV.clear();
524 }
525 }
527 std::vector<Trk::VertexID> cnstV; cnstV.clear();
530 }
531 }
532
534
535 bool pass = false;
536 if (result != nullptr) {
537 for(
auto v :
result->vertices()) {
538 if(
v->nTrackParticles()==0) {
539 std::vector<ElementLink<xAOD::TrackParticleContainer> > nullLinkVector;
540 v->setTrackParticleLinks(nullLinkVector);
541 }
542 }
543
545
546
547 result->setSVOwnership(
true);
548
549
552
553 if(chi2CutPassed) {
554 cascadeinfoContainer->push_back(
result.release());
555 pass = true;
556 }
557 }
558
559
560 if(pass) {
562 std::unique_ptr<Trk::IVKalState> state (
m_iVertexFitter->makeState(ctx));
564 std::vector<Trk::VertexID> vrtList_nc;
565
567 vrtList_nc.push_back(vID1_nc);
569 vrtList_nc.push_back(vID2_nc);
570
571 std::vector<const xAOD::TrackParticle*> tp; tp.clear();
572 std::vector<double> tp_masses; tp_masses.clear();
574
575 std::unique_ptr<Trk::VxCascadeInfo> result_nc(
m_iVertexFitter->fitCascade(*state));
576
577 if (result_nc != nullptr) {
578 for(auto v : result_nc->vertices()) {
579 if(
v->nTrackParticles()==0) {
580 std::vector<ElementLink<xAOD::TrackParticleContainer> > nullLinkVector;
581 v->setTrackParticleLinks(nullLinkVector);
582 }
583 }
584
586
587
588 result_nc->setSVOwnership(true);
589 cascadeinfoContainer_noConstr->push_back(result_nc.release());
590 }
591 else cascadeinfoContainer_noConstr->push_back(0);
592 }
593 else cascadeinfoContainer_noConstr->push_back(0);
594 }
595 }
596
597 return StatusCode::SUCCESS;
598 }
virtual double phi() const override final
The azimuthal angle ( ) of the particle (has range to .).
virtual double pt() const override final
The transverse momentum ( ) of the particle.
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
size_t nTrackParticles() const
Get the number of tracks associated with this vertex.
const TrackParticle * trackParticle(size_t i) const
Get the pointer to a given track that was used in vertex reco.
float j(const xAOD::IParticle &, const xAOD::TrackMeasurementValidation &hit, const Eigen::Matrix3d &jab_inv)
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.
1.16.1