AliReducedJetParticle.cxx

Go to the documentation of this file.

/**************************************************************************
 * Copyright(c) 1998-2015, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/
#include <TLorentzVector.h>

#include "AliReducedJetParticle.h"

ClassImp(HighPtTracks::AliReducedJetParticle)

namespace HighPtTracks {

AliReducedJetParticle::AliReducedJetParticle() :
  TObject(),
  fPx(0),
  fPy(0),
  fPz(0),
  fE(0),
  fDr(0),
  fPdgCode(0),
  fNTPCTrackReferences(0),
  fMatchedTrackContainer(NULL)
{
  fMatchedTrackContainer = new TObjArray();
  fMatchedTrackContainer->SetOwner();
}

AliReducedJetParticle::AliReducedJetParticle(double px, double py, double pz, double e, int pdgcode) :
  TObject(),
  fPx(px),
  fPy(py),
  fPz(pz),
  fE(e),
  fDr(0),
  fPdgCode(pdgcode),
  fNTPCTrackReferences(0),
  fMatchedTrackContainer(NULL)
{
  fMatchedTrackContainer = new TObjArray();
  fMatchedTrackContainer->SetOwner();
}

AliReducedJetParticle::AliReducedJetParticle(const AliReducedJetParticle& ref):
  TObject(ref),
  fPx(ref.fPx),
  fPy(ref.fPy),
  fPz(ref.fPz),
  fE(ref.fE),
  fDr(ref.fDr),
  fPdgCode(ref.fPdgCode),
  fNTPCTrackReferences(ref.fNTPCTrackReferences),
  fMatchedTrackContainer(NULL)
{
  fMatchedTrackContainer = new TObjArray;
  fMatchedTrackContainer->SetOwner();

  for(TIter trackIter = TIter(ref.fMatchedTrackContainer).Begin(); trackIter != TIter::End(); ++trackIter){
    AliReducedMatchedTrack *othermatched = static_cast<AliReducedMatchedTrack *>(*trackIter);
    fMatchedTrackContainer->Add(new AliReducedMatchedTrack(*othermatched));
  }
}

AliReducedJetParticle& AliReducedJetParticle::operator=(const AliReducedJetParticle& ref) {
  if(this != &ref){
    delete fMatchedTrackContainer;

    fMatchedTrackContainer = new TObjArray;
    fMatchedTrackContainer->SetOwner();

    for(TIter trackIter = TIter(ref.fMatchedTrackContainer).Begin(); trackIter != TIter::End(); ++trackIter){
      AliReducedMatchedTrack *othermatched = static_cast<AliReducedMatchedTrack *>(*trackIter);
      fMatchedTrackContainer->Add(new AliReducedMatchedTrack(*othermatched));
    }

  }
  return *this;
}

AliReducedJetParticle::~AliReducedJetParticle() {
  delete fMatchedTrackContainer;
}

void AliReducedJetParticle::FillLorentzVector(TLorentzVector& ref) const {
  ref.SetPxPyPzE(fPx, fPy, fPz, fE);
}

void AliReducedJetParticle::AddMatchedTrack(AliReducedMatchedTrack *trk) {
  fMatchedTrackContainer->Add(trk);
}

double AliReducedJetParticle::GetDeltaPt(int itrk) const {
  if(itrk >=  fMatchedTrackContainer->GetEntries()) return -1000.;
  TLorentzVector partvec;
  partvec.SetPxPyPzE(fPx, fPy, fPz, fE);
  double particlePt = partvec.Pt();
  AliReducedMatchedTrack *mytrk = static_cast<AliReducedMatchedTrack *>(fMatchedTrackContainer->At(itrk));
  return (mytrk->Pt()-particlePt)/particlePt;
}

} /* namespace HighPtTracks */