BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
D0Topipipi0.h
Go to the documentation of this file.
1//--------------------------------------------------------------------------
2//
3// Environment:
4// This software is part of models developed at BES collaboration
5// based on the EvtGen framework. If you use all or part
6// of it, please give an appropriate acknowledgement.
7//
8// Copyright Information: See EvtGen/BesCopyright
9// Copyright (A) 2006 Ping Rong-Gang @IHEP
10//
11// Module: EvtD0Topipipi0.cc
12//
13// Description: Model provided by user, see the mannual
14//
15// Modification history:
16//
17// Liaoyuan Dong Aug. 8, 2022 Module created
18//
19//------------------------------------------------------------------------
20
21#ifndef D0TOPIPIPI0_H
22#define D0TOPIPIPI0_H
23
24#include <complex>
25#include <vector>
26
27using namespace std;
28
29class D0Topipipi0 {
30
31public:
32 D0Topipipi0() {}
33 virtual ~D0Topipipi0();
34
35 void init();
36
37 complex<double> Amp( vector<double> Pip, vector<double> Pim, vector<double> Pi0 );
38
39private:
40 int _nd;
41 int charm;
42 int tagmode;
43
44 void setInput( double* pip, double* pim, double* pi0 );
45 int CalAmp();
46
47 complex<double> m_AmpD0;
48 complex<double> m_AmpDb;
49 complex<double> m_RAmp;
50 complex<double> Get_AmpD0() { return m_AmpD0; }
51 complex<double> Get_AmpDb() { return m_AmpDb; }
52
53 double Get_strongPhase();
54 double arg( complex<double> x );
55 double mag2( complex<double> x );
56 double mag2_itf( complex<double> x, complex<double> y );
57
58 complex<double> CalD0Amp();
59 complex<double> CalDbAmp();
60
61 vector<double> sum_tensor( vector<double> pa, vector<double> pb );
62 double contract_11_0( vector<double> pa, vector<double> pb );
63 double contract_22_0( vector<double> pa, vector<double> pb );
64 vector<double> contract_21_1( vector<double> pa, vector<double> pb );
65 vector<double> contract_22_2( vector<double> pa, vector<double> pb );
66 vector<double> contract_31_2( vector<double> pa, vector<double> pb );
67 vector<double> contract_41_3( vector<double> pa, vector<double> pb );
68 vector<double> contract_42_2( vector<double> pa, vector<double> pb );
69
70 vector<double> OrbitalTensors( vector<double> pa, vector<double> pb, vector<double> pc,
71 double r, int rank );
72 vector<double> ProjectionTensors( vector<double> pa, int rank );
73
74 double fundecaymomentum( double mr2, double m1_2, double m2_2 );
75 double fundecaymomentum2( double mr2, double m1_2, double m2_2 );
76 complex<double> breitwigner( double mx2, double mr, double wr );
77 double h( double m, double q );
78 double dh( double m0, double q0 );
79 double f( double m0, double sx, double q0, double q );
80 double d( double m0, double q0 );
81 double wid( double mass, double sa, double sb, double sc, double r, int l );
82
83 complex<double> GS( double mx2, double mr, double wr, double m1_2, double m2_2, double r,
84 int l );
85 complex<double> irho( double mr2, double m1_2, double m2_2 );
86 complex<double> Flatte( double mx2, double mr, double g1, double g2, double m1a, double m1b,
87 double m2a, double m2b );
88 complex<double> RBW( double mx2, double mr, double wr, double m1_2, double m2_2, double r,
89 int l );
90 double rho22( double sc );
91 complex<double> rhoMTX( int i, int j, double s );
92 complex<double> KMTX( int i, int j, double s );
93 complex<double> IMTX( int i, int j );
94 complex<double> FMTX( double Kijx, double Kijy, double rhojjx, double rhojjy, int i, int j );
95 double FINVMTX( double s, double* FINVx, double* FINVy );
96 complex<double> PVTR( int ID, double s );
97 complex<double> Fvector( double sa, double s0, int l );
98
99 vector<double> p4_Pip;
100 vector<double> p4_Pim;
101 vector<double> p4_Pi0;
102
103 vector<double> g_uv;
104 vector<double> epsilon_uvmn;
105 vector<complex<double>> fitpara;
106
107 double math_pi;
108 double mass_Pion;
109
110 double rRes;
111 double rD;
112 double m_Pi;
113 double m2_Pi;
114 double m_Pi0;
115 double m2_Pi0;
116
117 double m0_rho7700;
118 double w0_rho7700;
119
120 double m0_rho770p;
121 double w0_rho770p;
122
123 double m0_f21270;
124 double w0_f21270;
125
126 double s0_prod;
127
128 double delta;
129};
130
131#endif
mass
double mass
Definition CosmicGenerator.cxx:128
g1
TF1 * g1
Definition DataBase/tau_mode.c:58
x
Double_t x[10]
Definition DataBase/tau_mode.c:49
s
XmlRpcServer s
Definition HelloServer.cpp:11
q
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
Definition KKsem.h:33
D0Topipipi0::~D0Topipipi0
virtual ~D0Topipipi0()
Definition D0Topipipi0.cxx:33
D0Topipipi0::Amp
complex< double > Amp(vector< double > Pip, vector< double > Pim, vector< double > Pi0)
Definition D0Topipipi0.cxx:1335