EvtB0toKsKK.h

Go to the documentation of this file.

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2012-2018 - Belle II Collaboration                        *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Jeremy Dalseno                                           *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include "EvtGenBase/EvtDecayAmp.hh"
#include "EvtGenBase/EvtParticle.hh"
 
#include <fstream>
 
namespace Belle2 {
  class EvtB0toKsKK : public  EvtDecayAmp {
 
  public:
 
    /* Constructor. */
    //EvtB0toKsKK() {}
    EvtB0toKsKK() : m_alpha_kpkmnr(0.0), m_alpha_kskpnr(0.0), m_alpha_kskmnr(0.0) {};
 
    /* Destructor. */
    virtual ~EvtB0toKsKK();
 
    std::string getName();  
 
    EvtDecayBase* clone();  
 
    void init();  
 
    void initProbMax();  
 
    void decay(EvtParticle* p); 
 
    EvtVector4R umu(const EvtVector4R& p4a, const EvtVector4R& p4b,
                    const EvtVector4R& p4c);  
    EvtVector4R Smu(const EvtVector4R& p4a, const EvtVector4R& p4b,
                    const EvtVector4R& p4c);  
    EvtVector4R Lmu(const EvtVector4R& p4a, const EvtVector4R& p4b,
                    const EvtVector4R& p4c);  
    EvtTensor4C gmunu_tilde(const EvtVector4R& p4a,
                            const EvtVector4R& p4b,
                            const EvtVector4R& p4c); 
    EvtTensor4C Tmunu(const EvtVector4R& p4a, const EvtVector4R& p4b,
                      const EvtVector4R& p4c);     
    EvtTensor4C Multiply(const EvtTensor4C& t1,
                         const EvtTensor4C& t2);   
    EvtTensor4C RaiseIndices(const EvtTensor4C& t);  
    void RaiseIndex(EvtVector4R& vector);            
    EvtTensor4C Mmunu(const EvtVector4R& p4a, const EvtVector4R& p4b,
                      const EvtVector4R& p4c);       
    double BWBF(const double& q, const unsigned int& L); 
    double BWBF(const double& q, const double& q0,
                const unsigned int& L);              
    EvtComplex BreitWigner(const double& m, const double& m0,
                           const double& Gamma0,
                           const double& q, const double& q0,
                           const unsigned int& L);   
    EvtVector4R Boost(const EvtVector4R& p4,
                      const EvtVector4R& boost);     
    double p(const double& mab, const double& M, const double& mc);  
    double q(const double& mab, const double& ma, const double& mb); 
    EvtComplex Flatte_k(const double& s, const double& m_h); 
    EvtComplex Flatte(const double& m, const double& m0);    
 
    EvtComplex A_f0ks(const EvtVector4R& p4ks,
                      const EvtVector4R& p4kp, const EvtVector4R& p4km); 
    EvtComplex A_phiks(const EvtVector4R& p4ks,
                       const EvtVector4R& p4kp, const EvtVector4R& p4km); 
    EvtComplex A_fxks(const EvtVector4R& p4ks,
                      const EvtVector4R& p4kp, const EvtVector4R& p4km); 
    EvtComplex A_chic0ks(const EvtVector4R& p4ks,
                         const EvtVector4R& p4kp, const EvtVector4R& p4km); 
    EvtComplex A_kknr(const EvtVector4R& p4k1, const EvtVector4R& p4k2,
                      const double& m_alpha_kk);    
 
  private:

    EvtComplex m_a_f0ks_;      
    EvtComplex m_a_phiks_;     
    EvtComplex m_a_fxks_;      
    EvtComplex m_a_chic0ks_;   
    EvtComplex m_a_kpkmnr_;    
    EvtComplex m_a_kskpnr_;    
    EvtComplex m_a_kskmnr_;    
 
    EvtComplex m_abar_f0ks_;    
    EvtComplex m_abar_phiks_;   
    EvtComplex m_abar_fxks_;    
    EvtComplex m_abar_chic0ks_; 
    EvtComplex m_abar_kpkmnr_;  
    EvtComplex m_abar_kskpnr_;  
    EvtComplex m_abar_kskmnr_;  
 
    double m_alpha_kpkmnr;      
    double m_alpha_kskpnr;      
    double m_alpha_kskmnr;      
 
    std::ofstream m_debugfile_; 
  }; 
 
} // Belle 2 Namespace