1#ifndef EVTDTOKKPIPI0_HH
2#define EVTDTOKKPIPI0_HH
16 void getName(std::string& name);
24 double calPDF(
double P1[],
double P2[],
double P3[],
double P4[],
double & Result);
25 double Spin_factor(
double P1[],
double P2[],
double P3[],
double P4[], Int_t Ang,
double mass1,
double mass2,
double rRes,
double rD, Int_t mode);
27 void Com_Divide(
double a1[2],
double a2[2],
double res[2]);
28 double CalRho4pi(
double s);
30 TComplex propogator(
double mass,
double width,
double sx)
const;
31 double wid(
double mass,
double sa,
double sb,
double sc,
double r, Int_t l)
const;
32 double h(
double m,
double q)
const;
33 double dh(
double mass,
double q0)
const;
34 double f(
double mass,
double sx,
double q0,
double q)
const;
35 double d(
double mass,
double q0)
const;
36 TComplex Flatte_rhoab(
double sa,
double sb,
double sc)
const;
37 double rhoab_L(
double sa,
double sb,
double sc, Int_t L)
const;
38 TComplex propagatorRBW(
double mass,
double width,
double sa,
double sb,
double sc,
double r, Int_t l, Int_t isvaried)
const;
39 TComplex propagatorGS(
double mass,
double width,
double sa,
double sb,
double sc,
double r, Int_t l)
const;
40 void propagatorFlatte(
double mass,
double width,
double sa,
double sb,
double sc,
int r,
double prop[2]);
41 TComplex propagatorK1270(
double mass,
double sx,
double *sb,
double *sc,
double *g)
const;
42 TComplex propagator980(
double mass,
double sx,
double *sb,
double *sc,
double *g)
const;
43 TComplex propagatorf600(
double mass,
double width,
double sx)
const;
44 TComplex propagatorsigma500(
double mass,
double width,
double sx,
double *g)
const;
45 double dot(
double *a1,
double *a2)
const;
46 double Qabcs(
double sa,
double sb,
double sc)
const;
48 double barrier(
double l,
double sa,
double sb,
double sc,
double r,
double mass)
const;
50 void calt1(
double daug1[],
double daug2[],
double t1[])
const;
51 void calt2(
double daug1[],
double daug2[],
double t2[][4])
const;
52 const double sck = 0.493677*0.493677;
53 const double scpi = 0.13957*0.13957;
54 const double snpi = 0.1349768*0.1349768;
55 const double seta = 0.547862*0.547862;
57 double G[4][4],E[4][4][4][4];
1313 double width1270[1206]={
2522 double width1400[1206]={
3731 double gpieta2, gKK2;
TFile f("ana_bhabha660a_dqa_mcPat_zy_old.root")
complex< double > TComplex
****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
void getName(std::string &name)
void decay(EvtParticle *p)