TauDecayModeModule.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 <framework/core/Module.h>
#include <mdst/dataobjects/MCParticle.h>
#include <framework/datastore/StoreArray.h>
#include <framework/datastore/StoreObjPtr.h>
#include <framework/dataobjects/EventMetaData.h>
#include <string>
#include <analysis/dataobjects/TauPairDecay.h>
namespace Belle2 {
  class TauDecayModeModule : public Module {
 
  public:
 
    TauDecayModeModule();
 
    virtual void initialize() override;
 
    virtual void event() override;
 
  protected:
 
    std::string m_printmode; 
  private:
 
    StoreObjPtr<TauPairDecay> m_tauDecay;
 
    StoreObjPtr<EventMetaData> m_event_metadata; 
    int m_taum_no;     
    int m_taup_no;     
    StoreArray<MCParticle> MCParticles;
    std::map<std::string, int> mode_decay_minus;
    std::map<std::string, int> mode_decay_plus;
    std::string m_tauminusdecaymode;
    std::string m_tauplusdecaymode;
 
    Int_t m_mmode;
    Int_t m_pmode;
 
    Int_t m_mprong;
    Int_t m_pprong;
 
    double m_megstar;
 
    double m_pegstar;
 
    bool tauPair;
    int numOfTauMinus;
    int numOfTauPlus;
    int idOfTauMinus;
    int idOfTauPlus;
    std::string m_file_minus;
    std::string m_file_plus;
 
    static constexpr int finalStatePDGs[5] = { 11, 13, 211, 321, 2212 };
 
    static constexpr int Neutrinos[3] = { 12, 14, 16 };
 
    static constexpr int OrderedList[46] = {
      16, -16, 14, -14, 12, -12, // neutrinos
      11, -11, 13, -13, -211, 211, -321, 321, -2212, 2212, // charged final state particles
      111, 310, 130, 221, 223, 331, 333, //neutral mesons
      -213, 213, 113, //rho
      -323, 323, 313, -313, //Kstar
      -20213, 20213, //a1
      -9000211, 9000211, 9000111, -10211, 10211, 10111, //a0
      -10213, 10213, //b1
      20223, 9010221, //f1,f0
      3122, -3122, //lambda
      94144,//alpha
      22//gamma
    };
 
    bool m_isEtaPizPizPizFromTauMinus;
    bool m_isEtaPizPizPizFromTauPlus;
 
    bool m_isOmegaPimPipFromTauMinus;
    bool m_isOmegaPimPipFromTauPlus;
 
    //
    std::vector<int> vec_nut;
    std::vector<int> vec_anut;
    std::vector<int> vec_numu;
    std::vector<int> vec_anumu;
    std::vector<int> vec_nue; 
    std::vector<int> vec_anue; 
    std::vector<int> vec_em;
    std::vector<int> vec_ep;
    std::vector<int> vec_mum;
    std::vector<int> vec_mup;
    std::vector<int> vec_pim;
    std::vector<int> vec_pip;
    std::vector<int> vec_km;
    std::vector<int> vec_kp;
    std::vector<int> vec_apro;
    std::vector<int> vec_pro;
    std::vector<int> vec_pi0;
    std::vector<int> vec_k0s;
    std::vector<int> vec_k0l;
    std::vector<int> vec_eta;
    std::vector<int> vec_omega;
    std::vector<int> vec_etapr;
    std::vector<int> vec_phi;
    std::vector<int> vec_rhom;
    std::vector<int> vec_rhop;
    std::vector<int> vec_rho0;
    std::vector<int> vec_kstarm;
    std::vector<int> vec_kstarp;
    std::vector<int> vec_kstar0; 
    std::vector<int> vec_kstar0_br; 
    std::vector<int> vec_a1m;
    std::vector<int> vec_a1p;
    std::vector<int> vec_a00_980;
    std::vector<int> vec_a0m_980;
    std::vector<int> vec_a0p_980;
    std::vector<int> vec_a00_1450;
    std::vector<int> vec_a0m_1450;
    std::vector<int> vec_a0p_1450;
    std::vector<int> vec_b1m;
    std::vector<int> vec_b1p; 
    std::vector<int> vec_f1;
    std::vector<int> vec_f0;
    std::vector<int> vec_lambda;
    std::vector<int> vec_lmb_br;
    std::vector<int> vec_alpha;
    std::vector<int> vec_gam;
    std::vector<int> vec_radgam_taum;
    std::vector<int> vec_radgam_taup;
    std::vector<int> vec_dau_tauminus;
    std::vector<int> vec_dau_tauplus;
    void AnalyzeTauPairEvent();
    int TauolaBelle2DecayMode(const std::string& s, int chg);
    int getRecursiveMotherCharge(const MCParticle* mc);
    void IdentifyTauPair();
    int getProngOfDecay(const MCParticle& mc);
    double getEgstar(const std::vector<int>& vec_radgam, const MCParticle& mc);
  };
 
}