366 {
368 assert(cascadeinfoContainer!=nullptr && cascadeinfoContainer_noConstr!=nullptr);
369
370
373
374 std::vector<const xAOD::TrackParticle*> tracksJpsi1;
375 std::vector<const xAOD::TrackParticle*> tracksJpsi2;
376 std::vector<const xAOD::TrackParticle*> tracksDiTrk1;
377 std::vector<const xAOD::TrackParticle*> tracksDiTrk2;
378 std::vector<const xAOD::TrackParticle*> tracksPsi1;
379 std::vector<const xAOD::TrackParticle*> tracksPsi2;
380 std::vector<double> massesPsi1;
385 std::vector<double> massesPsi2;
390
391
394
395
398
399
400 std::vector<const xAOD::Vertex*> selectedPsi2Candidates;
401 for(auto vxcItr=psi2Container.cptr()->cbegin(); vxcItr!=psi2Container.cptr()->cend(); ++vxcItr) {
402
407 if(flagAcc.isAvailable(*vtx) && flagAcc(*vtx)) {
409 }
410 }
412
413
414 double mass_psi2 =
m_V0Tools->invariantMass(*vxcItr, massesPsi2);
415 if (mass_psi2 < m_psi2MassLower || mass_psi2 >
m_psi2MassUpper)
continue;
416
417
418 TLorentzVector p4_mu1, p4_mu2;
419 p4_mu1.SetPtEtaPhiM( (*vxcItr)->trackParticle(0)->pt(),
420 (*vxcItr)->trackParticle(0)->eta(),
422 p4_mu2.SetPtEtaPhiM( (*vxcItr)->trackParticle(1)->pt(),
423 (*vxcItr)->trackParticle(1)->eta(),
425 double mass_jpsi2 = (p4_mu1 + p4_mu2).M();
426 if (mass_jpsi2 < m_jpsi2MassLower || mass_jpsi2 >
m_jpsi2MassUpper)
continue;
427
429 TLorentzVector p4_trk1, p4_trk2;
430 p4_trk1.SetPtEtaPhiM( (*vxcItr)->trackParticle(2)->pt(),
431 (*vxcItr)->trackParticle(2)->eta(),
433 p4_trk2.SetPtEtaPhiM( (*vxcItr)->trackParticle(3)->pt(),
434 (*vxcItr)->trackParticle(3)->eta(),
436 double mass_diTrk2 = (p4_trk1 + p4_trk2).M();
438 }
439
440 double chi2DOF = (*vxcItr)->chiSquared()/(*vxcItr)->numberDoF();
442
443 selectedPsi2Candidates.push_back(*vxcItr);
444 }
445 if(selectedPsi2Candidates.size()==0) return StatusCode::SUCCESS;
446
447
448 std::vector<const xAOD::Vertex*> selectedPsi1Candidates;
449 for(auto vxcItr=psi1Container.cptr()->cbegin(); vxcItr!=psi1Container.cptr()->cend(); ++vxcItr) {
450
455 if(flagAcc.isAvailable(*vtx) && flagAcc(*vtx)) {
457 }
458 }
460
461
462 double mass_psi1 =
m_V0Tools->invariantMass(*vxcItr,massesPsi1);
463 if(mass_psi1 < m_psi1MassLower || mass_psi1 >
m_psi1MassUpper)
continue;
464
465
466 TLorentzVector p4_mu1, p4_mu2;
467 p4_mu1.SetPtEtaPhiM( (*vxcItr)->trackParticle(0)->pt(),
468 (*vxcItr)->trackParticle(0)->eta(),
470 p4_mu2.SetPtEtaPhiM( (*vxcItr)->trackParticle(1)->pt(),
471 (*vxcItr)->trackParticle(1)->eta(),
473 double mass_jpsi1 = (p4_mu1 + p4_mu2).M();
474 if (mass_jpsi1 < m_jpsi1MassLower || mass_jpsi1 >
m_jpsi1MassUpper)
continue;
475
477 TLorentzVector p4_trk1, p4_trk2;
478 p4_trk1.SetPtEtaPhiM( (*vxcItr)->trackParticle(2)->pt(),
479 (*vxcItr)->trackParticle(2)->eta(),
481 p4_trk2.SetPtEtaPhiM( (*vxcItr)->trackParticle(3)->pt(),
482 (*vxcItr)->trackParticle(3)->eta(),
484 double mass_diTrk1 = (p4_trk1 + p4_trk2).M();
486 }
487
488 double chi2DOF = (*vxcItr)->chiSquared()/(*vxcItr)->numberDoF();
490
491 selectedPsi1Candidates.push_back(*vxcItr);
492 }
493 if(selectedPsi1Candidates.size()==0) return StatusCode::SUCCESS;
494
495 std::vector<std::pair<const xAOD::Vertex*, const xAOD::Vertex*> > candidatePairs;
496 for(auto psi1Itr=selectedPsi1Candidates.cbegin(); psi1Itr!=selectedPsi1Candidates.cend(); ++psi1Itr) {
497 tracksPsi1.clear();
498 for(
size_t i=0;
i<(*psi1Itr)->nTrackParticles();
i++) tracksPsi1.push_back((*psi1Itr)->trackParticle(i));
499 for(auto psi2Itr=selectedPsi2Candidates.cbegin(); psi2Itr!=selectedPsi2Candidates.cend(); ++psi2Itr) {
501 for(
size_t j=0;
j<(*psi2Itr)->nTrackParticles();
j++) {
502 if(std::find(tracksPsi1.cbegin(), tracksPsi1.cend(), (*psi2Itr)->trackParticle(j)) != tracksPsi1.cend()) {
skip =
true;
break; }
503 }
506 for(
size_t ic=0;
ic<candidatePairs.size();
ic++) {
509 if((psi1Vertex == *psi1Itr && psi2Vertex == *psi2Itr) || (psi1Vertex == *psi2Itr && psi2Vertex == *psi1Itr)) {
skip =
true;
break; }
510 }
511 }
513 candidatePairs.push_back(std::pair<const xAOD::Vertex*, const xAOD::Vertex*>(*psi1Itr,*psi2Itr));
514 }
515 }
516
517 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(); } );
519 candidatePairs.erase(candidatePairs.begin()+
m_maxCandidates, candidatePairs.end());
520 }
521
522 for(
size_t ic=0;
ic<candidatePairs.size();
ic++) {
525
526 tracksPsi1.clear();
528 if (tracksPsi1.size() != massesPsi1.size()) {
529 ATH_MSG_ERROR(
"Problems with Psi1 input: number of tracks or track mass inputs is not correct!");
530 }
531 tracksPsi2.clear();
533 if (tracksPsi2.size() != massesPsi2.size()) {
534 ATH_MSG_ERROR(
"Problems with Psi2 input: number of tracks or track mass inputs is not correct!");
535 }
536
537 tracksJpsi1.clear();
540 tracksDiTrk1.clear();
544 }
545 tracksJpsi2.clear();
548 tracksDiTrk2.clear();
552 }
553
554 TLorentzVector p4_moth;
559 }
563 }
565
566
567 std::unique_ptr<Trk::IVKalState> state =
m_iVertexFitter->makeState(ctx);
568
569 int robustness = 0;
571
572
573 std::vector<Trk::VertexID> vrtList;
574
576
579 } else {
581 }
582 vrtList.push_back(vID1);
583
587 } else {
589 }
590 vrtList.push_back(vID2);
591
592 std::vector<const xAOD::TrackParticle*> tp; tp.clear();
593 std::vector<double> tp_masses; tp_masses.clear();
596 std::vector<Trk::VertexID> cnstV; cnstV.clear();
599 }
600 }
602 std::vector<Trk::VertexID> cnstV; cnstV.clear();
605 }
606 }
608 std::vector<Trk::VertexID> cnstV; cnstV.clear();
611 }
612 }
614 std::vector<Trk::VertexID> cnstV; cnstV.clear();
617 }
618 }
619
621
622 bool pass = false;
623 if (result != nullptr) {
624 for(
auto v :
result->vertices()) {
625 if(
v->nTrackParticles()==0) {
626 std::vector<ElementLink<xAOD::TrackParticleContainer> > nullLinkVector;
627 v->setTrackParticleLinks(nullLinkVector);
628 }
629 }
630
632
633
634 result->setSVOwnership(
true);
635
636
639
640 if(chi2CutPassed) {
641 cascadeinfoContainer->push_back(
result.release());
642 pass = true;
643 }
644 }
645
646
647 if(pass) {
649 std::unique_ptr<Trk::IVKalState> state (
m_iVertexFitter->makeState(ctx));
651 std::vector<Trk::VertexID> vrtList_nc;
652
654 vrtList_nc.push_back(vID1_nc);
655
657 vrtList_nc.push_back(vID2_nc);
658
659 std::vector<const xAOD::TrackParticle*> tp; tp.clear();
660 std::vector<double> tp_masses; tp_masses.clear();
662
663 std::unique_ptr<Trk::VxCascadeInfo> result_nc(
m_iVertexFitter->fitCascade(*state));
664
665 if (result_nc != nullptr) {
666 for(auto v : result_nc->vertices()) {
667 if(
v->nTrackParticles()==0) {
668 std::vector<ElementLink<xAOD::TrackParticleContainer> > nullLinkVector;
669 v->setTrackParticleLinks(nullLinkVector);
670 }
671 }
672
674
675
676 result_nc->setSVOwnership(true);
677 cascadeinfoContainer_noConstr->push_back(result_nc.release());
678 }
679 else cascadeinfoContainer_noConstr->push_back(0);
680 }
681 else cascadeinfoContainer_noConstr->push_back(0);
682 }
683 }
684
685 return StatusCode::SUCCESS;
686 }
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