1 #include "AliKFConversionMother.h"
2 #include "AliKFParticle.h"
8 ClassImp(AliKFConversionMother)
10 AliKFConversionMother::AliKFConversionMother() :
22 /*AliKFConversionMother::AliKFConversionMother(const AliKFParticle& d1, const AliKFParticle& d2) :
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);
37 if((d1.GetE()+d2.GetE()) != 0){
38 fAlpha=TMath::Abs((d1.GetE()-d2.GetE())/(d1.GetE()+d2.GetE()));
42 AliKFConversionMother::AliKFConversionMother(const AliKFConversionPhoton& d1, const AliKFConversionPhoton& d2) :
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);
58 if((d1.GetE()+d2.GetE()) != 0){
59 fAlpha=TMath::Abs((d1.GetE()-d2.GetE())/(d1.GetE()+d2.GetE()));
63 AliKFConversionMother::AliKFConversionMother(const AliKFConversionMother & original) :
64 AliKFParticle(original),
65 fMCLabel(original.fMCLabel),
66 fOpeningAngle(original.fOpeningAngle),
67 fAlpha(original.fAlpha)
70 fLabel[0] = original.fLabel[0];
71 fLabel[1] = original.fLabel[1];
75 AliKFConversionMother & AliKFConversionMother::operator = (const AliKFConversionMother & /*source*/)
77 // assignment operator
81 Double_t AliKFConversionMother::GetRapidity()
84 if(GetE() - GetPz() <= 0 || GetE() + GetPz() <= 0){
85 AliWarning("Error: |Pz| > E !!!! ");
88 rapidity = 0.5*(TMath::Log((GetE() +GetPz()) / (GetE()-GetPz())));
94 Double_t AliKFConversionMother::Phi() const
96 Double_t phi = AliKFParticle::GetPhi();
97 if (phi < 0.) phi += 2. * TMath::Pi();