MCMatching.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 *                                                                        *
 * Copyright(C) 2013-2018 - Belle II Collaboration                        *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Anze Zupanc, Christian Pulvermacher, Yo Sato             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#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 Belle2::Particle* particle);
 
    static int getMCErrors(const Belle2::Particle* particle, const Belle2::MCParticle* mcParticle = nullptr);
 
 
    static int setMCErrorsExtraInfo(Belle2::Particle* particle, const Belle2::MCParticle* mcParticle);
 
    static void fillGenMothers(const Belle2::MCParticle* mcP, std::vector<int>& genMCPMothers);
 
    static int findCommonMother(const Belle2::MCParticle* mcP, const std::vector<int>& firstMothers, int lastMother);
 
    static bool isFSP(int pdg);
 
    static bool isFSR(const Belle2::MCParticle* p);
 
    static bool isFSRLegacy(const Belle2::MCParticle* p);
 
    static bool isRadiativePhoton(const Belle2::MCParticle* p);
 
    static int getMissingParticleFlags(const Belle2::Particle* particle, const Belle2::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 Belle2::Particle* particle, const Belle2::MCParticle* mcParticle,
                                    const std::vector<int>& daughterPDG);
 
    static int getFlagsIgnoredByProperty(const Belle2::Particle* particle);
 
    static int makeDaughterAcceptMask(int daughterProperty);
 
  };
}