MCTruthVariables.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <Math/Vector4D.h>
 
#include <vector>
#include <map>
 
namespace Belle2 {
  class MCParticle;
  class Particle;
 
  namespace Variable {
 
    double particleMCInitialParticle(const Particle* particle);
 
    double particleMCVirtualParticle(const Particle* particle);
 
    double particleMCFSRParticle(const Particle* particle);
 
    double particleMCPhotosParticle(const Particle* particle);
 
    double particleMCISRParticle(const Particle* particle);
 
    double isSignal(const Particle* part);
 
    double isSignalAcceptWrongFSPs(const Particle* part);
 
    double isSignalAcceptMissingNeutrino(const Particle* part);
 
    double isSignalAcceptMissingMassive(const Particle* part);
 
    double isSignalAcceptMissingGamma(const Particle* part);
 
    double isSignalAcceptBremsPhotons(const Particle* part);
 
    double isSignalAcceptMissing(const Particle* part);
 
    double isWrongCharge(const Particle* particle);
 
    double isCloneTrack(const Particle* particle);
 
    double isOrHasCloneTrack(const Particle* particle);
 
    double isMisidentified(const Particle* particle);
 
    double genNthMotherPDG(const Particle* part, const std::vector<double>& daughterIDs);
 
    double genNthMotherIndex(const Particle* part, const std::vector<double>& daughterIDs);
 
    double genQ2PmPd(const Particle* part, const std::vector<double>& daughter_indices);
 
    double genMotherPDG(const Particle* particle);
 
    double genMotherIndex(const Particle* particle);
 
    double genMotherP(const Particle* particle);
 
    double genParticleIndex(const Particle* particle);
 
    double particleMCMatchPDGCode(const Particle* particle);
 
    double particleMCErrors(const Particle* particle);
 
    double particleMCMatchWeight(const Particle* particle);
 
    double particleNumberOfMCMatch(const Particle* particle);
 
    double particleMCPrimaryParticle(const Particle* particle);
 
    /*
     * return the true momentum transfer to lepton pair in a B (semi-) leptonic B meson decay
     */
    double particleMCMomentumTransfer2(const Particle* part);
 
    double particleMCMatchDecayTime(const Particle* particle);
 
 
    double particleMCMatchLifeTime(const Particle* particle);
 
    double particleMCMatchPX(const Particle* particle);
 
    double particleMCMatchPY(const Particle* particle);
 
    double particleMCMatchPZ(const Particle* particle);
 
    double particleMCMatchPT(const Particle* particle);
 
    double particleMCMatchE(const Particle* particle);
 
    double particleMCMatchP(const Particle* particle);
 
    ROOT::Math::PxPyPzEVector MCInvisibleP4(const MCParticle* mcparticle);
 
    double particleMCCosThetaBetweenParticleAndNominalB(const Particle* particle);
 
    double particleMCRecoilMass(const Particle* particle);
 
    double generatorEventWeight(const Particle*);
 
    int tauPlusMcMode(const Particle*);
 
    int tauMinusMcMode(const Particle*);
 
    int tauPlusMcProng(const Particle*);
 
    int tauMinusMcProng(const Particle*);
 
    double isReconstructible(const Particle*);
 
    double isTrackFound(const Particle*);
 
    double seenInPXD(const Particle*);
 
    double seenInSVD(const Particle*);
 
    double seenInCDC(const Particle*);
 
    double seenInTOP(const Particle*);
 
    double seenInECL(const Particle*);
 
    double seenInARICH(const Particle*);
 
    double seenInKLM(const Particle*);
 
    int genNStepsToDaughter(const Particle* particle, const std::vector<double>& arguments);
 
    int genNMissingDaughter(const Particle* particle, const std::vector<double>& arguments);
 
    double getHEREnergy(const Particle*);
 
    double getLEREnergy(const Particle*);
 
    double getCrossingAngleX(const Particle*);
 
    double getCrossingAngleY(const Particle*);
 
    double particleClusterMatchWeight(const Particle* particle);
 
    double particleClusterBestMCMatchWeight(const Particle*);
 
    double particleClusterBestMCPDGCode(const Particle*);
 
    double particleClusterTotalMCMatchWeight(const Particle*);
 
    void getKlongWeightMap(const Particle* particle, std::map<int, double>& mapMCParticleIndxAndWeight);
 
    double particleClusterTotalMCMatchWeightForKlong(const Particle*);
 
    double particleClusterTotalMCMatchWeightForBestKlong(const Particle*);
 
    double isBBCrossfeed(const Particle*);
 
    int ancestorBIndex(const Particle*);
  }
}