PWGGA/GammaConv/AliAODConversionMother.cxx

   1 #include "AliAODConversionMother.h"
   2 #include "AliKFConversionMother.h"
   3 #include "AliKFParticle.h"
   4
   5
   6 // Author D. Lohner (Daniel.Lohner@cern.ch)
   7
   8 using namespace std;
   9
  10 ClassImp(AliAODConversionMother)
  11
  12 AliAODConversionMother::AliAODConversionMother() :
  13 AliAODConversionParticle(),
  14 fChi2(-1),
  15 fOpeningAngle(-1),
  16     fAlpha(-1),
  17     fWeight(1)
  18 {
  19         fLabel[0] = -1;
  20         fLabel[1] = -1;
  21 }
  22
  23 AliAODConversionMother::AliAODConversionMother(AliKFConversionMother *kf):
  24 AliAODConversionParticle(),
  25 fChi2(kf->GetChi2()),
  26 fOpeningAngle(kf->GetOpeningAngle()),
  27 fAlpha(kf->GetAlpha()),
  28 fWeight(1)
  29 {
  30     // Set 4momentu
  31     SetPxPyPzE(kf->GetPx(),kf->GetPy(),kf->GetPz(),kf->GetE());
  32
  33      //Set Decay Photon Labels
  34     fLabel[0]=kf->GetGammaLabel(0);
  35     fLabel[1]=kf->GetGammaLabel(1);
  36 }
  37
  38 AliAODConversionMother::AliAODConversionMother(AliAODConversionPhoton *y1,AliAODConversionPhoton *y2):
  39 AliAODConversionParticle(),
  40 fChi2(-1),
  41 fOpeningAngle(-1),
  42 fAlpha(-1),
  43 fWeight(1)
  44 {
  45     // Set 4momentum
  46     SetPxPyPzE(y1->Px()+y2->Px(),y1->Py()+y2->Py(),y1->Pz()+y2->Pz(),y1->E()+y2->E());
  47
  48     // Calculate Opening Angle
  49     TVector3 v1(y1->Px(),y1->Py(),y1->Pz());
  50     TVector3 v2(y2->Px(),y2->Py(),y2->Pz());
  51     fOpeningAngle=v1.Angle(v2);
  52
  53     // Calculate Alpha
  54     if((y1->E()+y2->E()) != 0){
  55         fAlpha=TMath::Abs((y1->E()-y2->E())/(y1->E()+y2->E()));
  56     }
  57
  58     // Set Chi2 to the mean chi2 of gammas
  59  //   fChi2=0.5*(y1->GetChi2perNDF()+y2->GetChi2perNDF());
  60
  61     //Set Decay Photon Labels
  62     fLabel[0]=-1;
  63     fLabel[1]=-1;
  64
  65 }
  66
  67 AliAODConversionMother::~AliAODConversionMother() {
  68 // empty standard destructor
  69
  70 }
  71