GeneralTHDMunitarity.h

Go to the documentation of this file.

/* 
 * Copyright (C) 2016 HEPfit Collaboration
 * All rights reserved.
 *
 * For the licensing terms see doc/COPYING.
 */
 
#ifndef GENERALTHDMUNITARITY_H
#define GENERALTHDMUNITARITY_H
 
#include "ThObservable.h"
#include <gslpp.h>
 
class GeneralTHDM;
class GeneralTHDMcache;
 
class unitarity_GTHDM {
public:
    
    unitarity_GTHDM(const StandardModel& SM_i);
    
    virtual ~unitarity_GTHDM();
    
    bool CalcSeigen21(gslpp::matrix<gslpp::complex>& Seigvec_i, gslpp::vector<double>& Seigval_i);
    
    bool CalcSeigen01(gslpp::matrix<gslpp::complex>& Seigvec_i, gslpp::vector<double>& Seigval_i);
    
    bool CalcSeigen00(gslpp::matrix<gslpp::complex>& Seigvec_i, gslpp::vector<double>& Seigval_i);
    
    gslpp::vector<double> getSeigen21();
    
    gslpp::vector<double> getSeigen01();
    
    gslpp::vector<double> getSeigen00();
    
    
private:
    const GeneralTHDM& myGTHDM;
    
    gslpp::matrix<gslpp::complex> Smat21, Smat01, Smat00; 
    gslpp::matrix<gslpp::complex> Seigvec21, Seigvec01, Seigvec00;
    gslpp::vector<double> Seigval21, Seigval01, Seigval00;
};
 
 
class unitarity1_GTHDM: public ThObservable {
public:
 
    unitarity1_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
 
class unitarity2_GTHDM: public ThObservable {
public:
 
    unitarity2_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity3_GTHDM: public ThObservable {
public:
 
    unitarity3_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity5_GTHDM: public ThObservable {
public:
 
    unitarity5_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity6_GTHDM: public ThObservable {
public:
 
    unitarity6_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity7_GTHDM: public ThObservable {
public:
 
    unitarity7_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity8_GTHDM: public ThObservable {
public:
 
    unitarity8_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity9_GTHDM: public ThObservable {
public:
 
    unitarity9_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity10_GTHDM: public ThObservable {
public:
 
    unitarity10_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity11_GTHDM: public ThObservable {
public:
 
    unitarity11_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity12_GTHDM: public ThObservable {
public:
 
    unitarity12_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    unitarity_GTHDM myunitarity_GTHDM;
};
 
class unitarity4_GTHDM: public ThObservable {
public:
 
    unitarity4_GTHDM(const StandardModel& SM_i);
 
    double computeThValue();
 
private:
    const GeneralTHDM& myGTHDM;
};
 
#endif /* GENERALTHDMUNITARITY_H */