Fix for definitions for CINT
[u/mrichter/AliRoot.git] / TEvtGen / EvtGenBase / EvtRaritaSchwinger.hh
1 //--------------------------------------------------------------------------
2 //
3 // Environment:
4 //      This software is part of the EvtGen package developed jointly
5 //      for the BaBar and CLEO collaborations.  If you use all or part
6 //      of it, please give an appropriate acknowledgement.
7 //
8 // Copyright Information: See EvtGen/COPYRIGHT
9 //      Copyright (C) 2000      Caltech, UCSB
10 //
11 // Module: EvtGen/EvtRaritaSchwing.hh
12 //
13 // Description:Class to handle spin 3/2
14 //
15 // Modification history:
16 //
17 //    RYD     April 23, 2000         Module created
18 //
19 //------------------------------------------------------------------------
20
21 #ifndef EVTRARITASCHWINGER_HH
22 #define EVTRARITASCHWINGER_HH
23
24 #include "EvtGenBase/EvtComplex.hh"
25 #include "EvtGenBase/EvtVector4C.hh"
26 #include "EvtGenBase/EvtVector4R.hh"
27 #include "EvtGenBase/EvtDiracSpinor.hh"
28
29 class EvtRaritaSchwinger ;
30 EvtRaritaSchwinger rotateEuler(const EvtRaritaSchwinger& rs,
31                                double alpha,double beta,double gamma);
32 EvtRaritaSchwinger boostTo(const EvtRaritaSchwinger& rs,
33                            const EvtVector4R p4);
34 EvtRaritaSchwinger boostTo(const EvtRaritaSchwinger& rs,
35                            const EvtVector3R boost);
36 EvtRaritaSchwinger dirProd(EvtVector4R v,EvtDiracSpinor u);
37 EvtRaritaSchwinger dirProd(EvtVector4C v,EvtDiracSpinor u);
38 EvtRaritaSchwinger operator+(const EvtRaritaSchwinger& u1,
39                              const EvtRaritaSchwinger& u2); 
40 EvtRaritaSchwinger operator-(const EvtRaritaSchwinger& u1,
41                              const EvtRaritaSchwinger& u2); 
42 EvtComplex operator*(const EvtRaritaSchwinger& u1,
43                      const EvtRaritaSchwinger& u2);
44
45 class EvtRaritaSchwinger{
46
47   friend EvtRaritaSchwinger rotateEuler(const EvtRaritaSchwinger& rs,
48                                         double alpha,double beta,double gamma);
49   friend EvtRaritaSchwinger boostTo(const EvtRaritaSchwinger& rs,
50                                     const EvtVector4R p4);
51   friend EvtRaritaSchwinger boostTo(const EvtRaritaSchwinger& rs,
52                                     const EvtVector3R boost);
53
54   friend EvtRaritaSchwinger dirProd(EvtVector4R v,EvtDiracSpinor u);
55   friend EvtRaritaSchwinger dirProd(EvtVector4C v,EvtDiracSpinor u);
56
57   friend EvtRaritaSchwinger operator+(const EvtRaritaSchwinger& u1,
58                                   const EvtRaritaSchwinger& u2); 
59   friend EvtRaritaSchwinger operator-(const EvtRaritaSchwinger& u1,
60                                   const EvtRaritaSchwinger& u2); 
61
62   friend EvtComplex operator*(const EvtRaritaSchwinger& u1,
63                                   const EvtRaritaSchwinger& u2); 
64
65 public:
66
67   inline EvtRaritaSchwinger();
68   virtual ~EvtRaritaSchwinger();
69   inline EvtRaritaSchwinger(const EvtRaritaSchwinger& rs);
70   inline EvtRaritaSchwinger& operator=(const EvtRaritaSchwinger& rs);
71
72   void set(int i,int j,const EvtComplex& sp);
73
74   void applyRotateEuler(double alpha,double beta,double gamma);
75   void applyBoostTo(const EvtVector4R p4);
76   void applyBoostTo(const EvtVector3R boost);
77  
78   EvtRaritaSchwinger& operator+=(const EvtRaritaSchwinger& u2);
79   EvtRaritaSchwinger& operator-=(const EvtRaritaSchwinger& u2);
80
81   EvtComplex get(int i,int j) const; 
82   friend std::ostream& operator<<(std::ostream& s, const EvtRaritaSchwinger& rs); 
83
84   EvtVector4C getVector(int i) const;
85   EvtDiracSpinor getSpinor(int i) const;
86
87   void setVector(int i,const EvtVector4C& v);
88   void setSpinor(int i,const EvtDiracSpinor& sp);
89
90
91   
92 private:
93   
94   //First index in spinor index, second is Lorentz index.
95   EvtComplex _rs[4][4];
96
97 };
98
99 EvtRaritaSchwinger::EvtRaritaSchwinger(){
100
101   int i,j;
102   for(i=0;i<4;i++){
103     for(j=0;j<4;j++){
104       _rs[i][j]=0.0;
105     }
106   }
107
108 }
109
110 EvtRaritaSchwinger::EvtRaritaSchwinger(const EvtRaritaSchwinger& rs){
111
112   int i,j;
113   for(i=0;i<4;i++){
114     for(j=0;j<4;j++){
115       _rs[i][j]=rs._rs[i][j];
116     }
117   }
118
119 }
120
121 EvtRaritaSchwinger& EvtRaritaSchwinger::operator=(const EvtRaritaSchwinger& rs){
122
123   int i,j;
124   for(i=0;i<4;i++){
125     for(j=0;j<4;j++){
126       _rs[i][j]=rs._rs[i][j];
127     }
128   }
129
130   return *this;
131
132 }
133
134 #endif
135
136