MCMatching.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 <vector>
#include <string>
 
namespace Belle2 {
  class Particle;
  class MCParticle;

  struct MCMatching {
    static const std::string c_extraInfoMCErrors;

    enum MCErrorFlags {
      c_Correct               = 0,  
      c_MissFSR               = 1, 
      c_MissingResonance      = 2, 
      c_DecayInFlight         = 4, 
      c_MissNeutrino          = 8, 
      c_MissGamma             = 16, 
      c_MissMassiveParticle   = 32, 
      c_MissKlong             = 64, 
      c_MisID                 = 128, 
      c_AddedWrongParticle    = 256, 
      c_InternalError         = 512, 
      c_MissPHOTOS            = 1024, 
      c_AddedRecoBremsPhoton  = 2048, 
    };

    static std::string explainFlags(unsigned int flags);

    static bool setMCTruth(const Particle* particle);

    static int getMCErrors(const Particle* particle, const MCParticle* mcParticle = nullptr);
 

    static int setMCErrorsExtraInfo(Particle* particle, const MCParticle* mcParticle);

    static void fillGenMothers(const MCParticle* mcP, std::vector<int>& genMCPMothers);

    static int findCommonMother(const MCParticle* mcP, const std::vector<int>& firstMothers, int lastMother);

    static bool isFSP(int pdg);

    static bool isFSR(const MCParticle* p);

    static bool isFSRLegacy(const MCParticle* p);

    static bool isRadiativePhoton(const MCParticle* p);

    static int getMissingParticleFlags(const Particle* particle, const MCParticle* mcParticle);

    static int getFlagsOfFSP(const Particle* particle, const MCParticle* mcParticle);

    static int getFlagsOfDaughters(const Particle* daughter, const MCParticle* mcParticle);

    static int getFlagsOfBremsPhotonDaughter(const Particle* daughter, const MCParticle* mcParticle,
                                             const std::vector<const MCParticle*>& genParts);

    static int getNumberOfDaughtersWithoutNeutrinos(const MCParticle* mcParticle);

    static int countMissingParticle(const Particle* particle, const MCParticle* mcParticle,
                                    const std::vector<int>& daughterPDG);

    static int getFlagsIgnoredByProperty(const Particle* particle);

    static int makeDaughterAcceptMask(int daughterProperty);
 
  };

}