]>
Commit | Line | Data |
---|---|---|
2eedd4ed | 1 | #include "AliKFConversionMother.h" |
2 | #include "AliKFParticle.h" | |
3 | #include "TMath.h" | |
4 | #include "TVector3.h" | |
5 | ||
6 | using namespace std; | |
7 | ||
8 | ClassImp(AliKFConversionMother) | |
9 | ||
10 | AliKFConversionMother::AliKFConversionMother() : | |
11 | AliKFParticle(), | |
92efd725 | 12 | fMCLabel(-1), |
13 | fOpeningAngle(-1), | |
14 | fAlpha(-1) | |
2eedd4ed | 15 | |
16 | { | |
17 | //Default constructor | |
18 | fLabel[0] = -1; | |
19 | fLabel[1] = -1; | |
20 | } | |
21 | ||
22 | /*AliKFConversionMother::AliKFConversionMother(const AliKFParticle& d1, const AliKFParticle& d2) : | |
23 | AliKFParticle(d1,d2), | |
24 | fOpeningAngle(-1), | |
25 | fAlpha(-1) | |
26 | ||
27 | { | |
28 | //Default constructor | |
29 | fLabel[0] = -1; | |
30 | fLabel[1] = -1; | |
31 | ||
32 | // Calculate Opening Angle | |
33 | TVector3 v1(d1.GetPx(),d1.GetPy(),d1.GetPz()); | |
34 | TVector3 v2(d2.GetPx(),d2.GetPy(),d2.GetPz()); | |
35 | fOpeningAngle=v1.Angle(v2); | |
36 | // Calculate Alpha | |
37 | if((d1.GetE()+d2.GetE()) != 0){ | |
38 | fAlpha=TMath::Abs((d1.GetE()-d2.GetE())/(d1.GetE()+d2.GetE())); | |
39 | } | |
40 | }*/ | |
41 | ||
42 | AliKFConversionMother::AliKFConversionMother(const AliKFConversionPhoton& d1, const AliKFConversionPhoton& d2) : | |
43 | AliKFParticle(d1,d2), | |
92efd725 | 44 | fMCLabel(-1), |
2eedd4ed | 45 | fOpeningAngle(-1), |
46 | fAlpha(-1) | |
47 | ||
48 | { | |
49 | //Default constructor | |
50 | fLabel[0] = -1; | |
51 | fLabel[1] = -1; | |
52 | ||
53 | // Calculate Opening Angle | |
54 | TVector3 v1(d1.GetPx(),d1.GetPy(),d1.GetPz()); | |
55 | TVector3 v2(d2.GetPx(),d2.GetPy(),d2.GetPz()); | |
56 | fOpeningAngle=v1.Angle(v2); | |
57 | // Calculate Alpha | |
58 | if((d1.GetE()+d2.GetE()) != 0){ | |
59 | fAlpha=TMath::Abs((d1.GetE()-d2.GetE())/(d1.GetE()+d2.GetE())); | |
60 | } | |
61 | } | |
62 | ||
63 | AliKFConversionMother::AliKFConversionMother(const AliKFConversionMother & original) : | |
64 | AliKFParticle(original), | |
92efd725 | 65 | fMCLabel(original.fMCLabel), |
2eedd4ed | 66 | fOpeningAngle(original.fOpeningAngle), |
67 | fAlpha(original.fAlpha) | |
68 | { | |
69 | //Copy constructor | |
70 | fLabel[0] = original.fLabel[0]; | |
71 | fLabel[1] = original.fLabel[1]; | |
72 | } | |
73 | ||
74 | ||
75 | AliKFConversionMother & AliKFConversionMother::operator = (const AliKFConversionMother & /*source*/) | |
76 | { | |
77 | // assignment operator | |
78 | return *this; | |
79 | } | |
80 | ||
81 | Double_t AliKFConversionMother::GetRapidity() | |
82 | { | |
83 | Double_t rapidity; | |
84 | if(GetE() - GetPz() <= 0 || GetE() + GetPz() <= 0){ | |
92efd725 | 85 | AliWarning("Error: |Pz| > E !!!! "); |
2eedd4ed | 86 | rapidity=8.; |
87 | } else { | |
88 | rapidity = 0.5*(TMath::Log((GetE() +GetPz()) / (GetE()-GetPz()))); | |
89 | } | |
90 | ||
91 | return rapidity; | |
92 | } | |
93 | ||
94 | Double_t AliKFConversionMother::Phi() const | |
95 | { | |
96 | Double_t phi = AliKFParticle::GetPhi(); | |
97 | if (phi < 0.) phi += 2. * TMath::Pi(); | |
98 | return phi; | |
99 | } |