00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021 #ifndef EVTAMP_HH
00022 #define EVTAMP_HH
00023
00024 #include "EvtGenBase/EvtComplex.hh"
00025 #include "EvtGenBase/EvtPatches.hh"
00026
00027 class EvtSpinDensity;
00028 class EvtId;
00029
00030 class EvtAmp {
00031
00032 public:
00033
00034 EvtAmp();
00035 virtual ~EvtAmp(){}
00036
00037 EvtAmp(const EvtAmp& amp);
00038
00039 void init(EvtId p,int ndaug,EvtId *daug);
00040
00041 void setAmp(int *ind,const EvtComplex& amp);
00042
00043 const EvtComplex& getAmp(int *ind)const;
00044
00045 EvtSpinDensity getSpinDensity();
00046
00047 EvtSpinDensity contract(int i,const EvtAmp& a);
00048 EvtAmp contract(int i, const EvtSpinDensity& rho);
00049
00050
00051 EvtAmp contract(int i, const EvtAmp& a1,const EvtAmp& a2);
00052
00053
00054 EvtSpinDensity getForwardSpinDensity(EvtSpinDensity *rho_list,int k);
00055 EvtSpinDensity getBackwardSpinDensity(EvtSpinDensity *rho_list);
00056
00057 EvtAmp& operator=(const EvtAmp& amp);
00058
00059
00063 void vertex(const EvtComplex& amp);
00064
00068 void vertex(int i1, const EvtComplex& amp);
00069
00073 void vertex(int i1, int i2, const EvtComplex& amp);
00074
00078 void vertex(int i1, int i2, int i3, const EvtComplex& amp);
00079
00083 void vertex(int *i1, const EvtComplex& amp);
00084
00085
00086 void dump();
00087
00088 private:
00089
00090 friend class EvtDecayAmp;
00091 friend class EvtKstarstargamma;
00092
00093 void setNDaug(int n);
00094 void setNState(int parent_states,int *daug_states);
00095
00096
00097 EvtComplex _amp[450];
00098
00099
00100 int _ndaug;
00101
00102
00103 int _pstates;
00104
00105
00106 int dstates[10];
00107
00108
00109 int _dnontrivial[10];
00110
00111
00112 int _nontrivial;
00113
00114
00115 int _nstate[5];
00116
00117 };
00118
00119
00120 #endif
00121