Tool to calculate the momentum balance significance from the detailed information (parameters, scattering angles) of a track fit. More...

#include <MuonMomentumBalanceSignificanceTool.h>

Inheritance diagram for Rec::MuonMomentumBalanceSignificanceTool:

Collaboration diagram for Rec::MuonMomentumBalanceSignificanceTool:

Public Member Functions
	~MuonMomentumBalanceSignificanceTool ()=default
double	momentumBalanceSignificance (const xAOD::Muon &muon) const
	Calculate momentum (im)balance significance of a muon (method will simply step down to the relevant track.
double	momentumBalanceSignificance (const Trk::Track &track) const
	Calculate momentum (im)balance significance of a track, which needs to be a full ATLAS track.

Detailed Description

Tool to calculate the momentum balance significance from the detailed information (parameters, scattering angles) of a track fit.

Definition at line 20 of file MuonMomentumBalanceSignificanceTool.h.

Constructor & Destructor Documentation

◆ ~MuonMomentumBalanceSignificanceTool()

Rec::MuonMomentumBalanceSignificanceTool::~MuonMomentumBalanceSignificanceTool ( )

default

Member Function Documentation

◆ momentumBalanceSignificance() [1/2]

double Rec::MuonMomentumBalanceSignificanceTool::momentumBalanceSignificance ( const Trk::Track & track ) const

Calculate momentum (im)balance significance of a track, which needs to be a full ATLAS track.

Definition at line 35 of file MuonMomentumBalanceSignificanceTool.cxx.

                                                                                                       {
        // find the TSOS carrying caloEnergy (TSOS type CaloDeposit)
        // compare parameters with those from previous TSOS
        double energyBalance = 0.;
        const Trk::EnergyLoss* energyLoss = nullptr;
        const Trk::TrackParameters* previousParameters = nullptr;
        for (Trk::TrackStates::const_iterator s = track.trackStateOnSurfaces()->begin();
             s != track.trackStateOnSurfaces()->end(); ++s) {
            if (!(**s).trackParameters()) continue;
            if ((**s).materialEffectsOnTrack()) {
                const Trk::MaterialEffectsOnTrack* meot = dynamic_cast<const Trk::MaterialEffectsOnTrack*>((**s).materialEffectsOnTrack());
                if (!meot) continue;
                energyLoss = meot->energyLoss();
                if ((**s).type(Trk::TrackStateOnSurface::CaloDeposit) && energyLoss && previousParameters) {
                    energyBalance = previousParameters->momentum().mag() - energyLoss->deltaE() - (**s).trackParameters()->momentum().mag();
 
                    if (fabs(energyBalance) < energyLoss->sigmaDeltaE() || energyLoss->sigmaMinusDeltaE() <= 0.0 ||
                        energyLoss->sigmaPlusDeltaE() <= 0.0)
 
                    {
                        ATH_MSG_VERBOSE(std::setiosflags(std::ios::fixed)
                                        << " momentum balance  " << std::setw(6) << std::setprecision(2) << energyBalance / GeV
                                        << "   significance " << std::setw(6) << std::setprecision(1)
                                        << energyBalance / energyLoss->sigmaDeltaE() << "  p before/after calo" << std::setw(7)
                                        << std::setprecision(2) << previousParameters->momentum().mag() / GeV << " /" << std::setw(7)
                                        << std::setprecision(2) << (**s).trackParameters()->momentum().mag() / GeV
                                        << "   energy deposit sigma  " << energyLoss->sigmaDeltaE() / GeV << " GeV");
                        energyBalance /= energyLoss->sigmaDeltaE();
                        if (fabs(energyBalance) > 1.0)
                            ATH_MSG_DEBUG(" momentum balance using symmetric error "
                                          << "as asymmetric error undefined ");
                    } else if (energyBalance < 0.) {
                        ATH_MSG_VERBOSE(std::setiosflags(std::ios::fixed)
                                        << " momentum balance  " << std::setw(6) << std::setprecision(2) << energyBalance / GeV
                                        << "   significance " << std::setw(6) << std::setprecision(1)
                                        << energyBalance / energyLoss->sigmaDeltaE() << "  p before/after calo" << std::setw(7)
                                        << std::setprecision(2) << previousParameters->momentum().mag() / GeV << " /" << std::setw(7)
                                        << std::setprecision(2) << (**s).trackParameters()->momentum().mag() / GeV
                                        << "   energy deposit sigma- " << energyLoss->sigmaMinusDeltaE() / GeV << " GeV");
                        energyBalance /= energyLoss->sigmaMinusDeltaE();
                    } else {
                        ATH_MSG_VERBOSE(std::setiosflags(std::ios::fixed)
                                        << " momentum balance  " << std::setw(6) << std::setprecision(2) << energyBalance / GeV
                                        << "   significance " << std::setw(6) << std::setprecision(1)
                                        << energyBalance / energyLoss->sigmaDeltaE() << "  p before/after calo" << std::setw(7)
                                        << std::setprecision(2) << previousParameters->momentum().mag() / GeV << " /" << std::setw(7)
                                        << std::setprecision(2) << (**s).trackParameters()->momentum().mag() / GeV
                                        << "   energy deposit sigma+ " << energyLoss->sigmaPlusDeltaE() / GeV << " GeV");
                        energyBalance /= energyLoss->sigmaPlusDeltaE();
                    }
                    break;
                }
            }
 
            // update previous parameters
            previousParameters = (**s).trackParameters();
        }
 
        return energyBalance;
    }

◆ momentumBalanceSignificance() [2/2]

double Rec::MuonMomentumBalanceSignificanceTool::momentumBalanceSignificance ( const xAOD::Muon & muon ) const

Calculate momentum (im)balance significance of a muon (method will simply step down to the relevant track.

Definition at line 20 of file MuonMomentumBalanceSignificanceTool.cxx.

                                                                                                      {
        // exclude statistical combination
        if (muon.author() == xAOD::Muon::Author::STACO) return 0.;
        using enum xAOD::Muon::TrackParticleType;
        // get combined track
        const Trk::Track* theTrack =
            muon.trackParticle(CombinedTrackParticle) ? muon.trackParticle(CombinedTrackParticle)->track() : nullptr;
 
        if (theTrack == nullptr) {
            ATH_MSG_DEBUG("No primary author original track for refitted muon, stop calculation...");
            return 0.0;
        } else
            return momentumBalanceSignificance(*theTrack);
    }