[u/mrichter/AliRoot.git] / PWGGA / EMCALTasks / AliEmcalJet.cxx

// $Id$

#include "AliEmcalJet.h"

ClassImp(AliEmcalJet)

//__________________________________________________________________________________________________
AliEmcalJet::AliEmcalJet(Double_t px, Double_t py, Double_t pz) 
  : AliVParticle(), 
    fPt(TMath::Sqrt(px*px+py*py)), 
    fEta(TMath::ASinH(pz/fPt)),
    fPhi(0), fM(0), fNEF(0), 
    fArea(0), fNch(0), fNn(0),
  fMaxCPt(0), fMaxNPt(0), fClusterIDs(new TArrayI()), fTrackIDs(new TArrayI())
{    
  // Constructor.

  if (fPt != 0) {
    fPhi = TMath::ATan2(py, px);
    if (fPhi<0.) 
      fPhi += 2. * TMath::Pi();
  }
}

//_________________________________________________________________________________________________
AliEmcalJet::AliEmcalJet(Double_t pt, Double_t eta, Double_t phi, Double_t m) :
  AliVParticle(), 
  fPt(pt), fEta(eta), fPhi(phi), 
  fM(m), fNEF(0), fArea(0), 
  fNch(0), fNn(0), 
  fMaxCPt(0), fMaxNPt(0), fClusterIDs(new TArrayI()), fTrackIDs(new TArrayI())
{
  // Constructor.

 if (fPhi<0.) 
   fPhi += 2. * TMath::Pi();
}

//_________________________________________________________________________________________________
AliEmcalJet::AliEmcalJet(const AliEmcalJet &jet) :
  AliVParticle(jet),
  fPt(jet.fPt), fEta(jet.fEta), fPhi(jet.fPhi), 
  fM(jet.fM), fNEF(jet.fNEF), fArea(jet.fArea), 
  fNch(jet.fNch), fNn(jet.fNn),
  fMaxCPt(jet.fMaxCPt), fMaxNPt(jet.fMaxNPt), fClusterIDs(new TArrayI(*(jet.fClusterIDs))), fTrackIDs(new TArrayI(*(jet.fTrackIDs)))
{
  // Constructor.
}

//_________________________________________________________________________________________________
AliEmcalJet &AliEmcalJet::operator=(const AliEmcalJet &jet)
{
  // Assignment operator.

  if (this!=&jet) {
    AliVParticle::operator=(jet);
    fPt     = jet.fPt;
    fEta    = jet.fEta;
    fPhi    = jet.fPhi;
    fM      = jet.fM; 
    fNEF    = jet.fNEF;
    fArea   = jet.fArea; 
    fNch    = jet.fNch; 
    fNn     = jet.fNn;
    fMaxCPt = jet.fMaxCPt; 
    fMaxNPt = jet.fMaxNPt;
  }

  return *this;
}

//__________________________________________________________________________________________________
void AliEmcalJet::GetMom(TLorentzVector &vec) const
{
  // Return momentum as four vector.

  Double_t p = fPt *TMath::CosH(fEta);
  vec.SetPtEtaPhiE(fPt,fEta,fPhi,TMath::Sqrt(p*p+fM*fM));
}

//__________________________________________________________________________________________________
void AliEmcalJet::Print(Option_t* /*option*/) const
{
  // Print jet information.

  printf("Jet pt=%.2f, eta=%.2f, phi=%.2f, area=%.2f, NEF=%.2f\n", fPt, fEta, fPhi, fArea, fNEF);
}
Commit	Line	Data
7df864a3	1	// $Id$
7efbea04	2
914d486c	3	#include "AliEmcalJet.h"
7df864a3	4
914d486c	5	ClassImp(AliEmcalJet)
7efbea04	6
7efbea04	7	//__________________________________________________________________________________________________
914d486c	8	AliEmcalJet::AliEmcalJet(Double_t px, Double_t py, Double_t pz)
7efbea04	9	: AliVParticle(),
	10	fPt(TMath::Sqrt(pxpx+pypy)),
	11	fEta(TMath::ASinH(pz/fPt)),
688670de	12	fPhi(0), fM(0), fNEF(0),
688670de	13	fArea(0), fNch(0), fNn(0),
35789a2d	14	fMaxCPt(0), fMaxNPt(0), fClusterIDs(new TArrayI()), fTrackIDs(new TArrayI())
7efbea04	15	{
	16	// Constructor.
	17
	18	if (fPt != 0) {
	19	fPhi = TMath::ATan2(py, px);
	20	if (fPhi<0.)
	21	fPhi += 2. * TMath::Pi();
	22	}
	23	}
	24
9cffd2d0	25	//_________________________________________________________________________________________________
914d486c	26	AliEmcalJet::AliEmcalJet(Double_t pt, Double_t eta, Double_t phi, Double_t m) :
9cffd2d0	27	AliVParticle(),
	28	fPt(pt), fEta(eta), fPhi(phi),
	29	fM(m), fNEF(0), fArea(0),
	30	fNch(0), fNn(0),
35789a2d	31	fMaxCPt(0), fMaxNPt(0), fClusterIDs(new TArrayI()), fTrackIDs(new TArrayI())
9cffd2d0	32	{
	33	// Constructor.
	34
	35	if (fPhi<0.)
	36	fPhi += 2. * TMath::Pi();
	37	}
	38
7efbea04	39	//_________________________________________________________________________________________________
914d486c	40	AliEmcalJet::AliEmcalJet(const AliEmcalJet &jet) :
7efbea04	41	AliVParticle(jet),
	42	fPt(jet.fPt), fEta(jet.fEta), fPhi(jet.fPhi),
	43	fM(jet.fM), fNEF(jet.fNEF), fArea(jet.fArea),
	44	fNch(jet.fNch), fNn(jet.fNn),
35789a2d	45	fMaxCPt(jet.fMaxCPt), fMaxNPt(jet.fMaxNPt), fClusterIDs(new TArrayI((jet.fClusterIDs))), fTrackIDs(new TArrayI((jet.fTrackIDs)))
7df864a3	46	{
7efbea04	47	// Constructor.
	48	}
	49
	50	//_________________________________________________________________________________________________
914d486c	51	AliEmcalJet &AliEmcalJet::operator=(const AliEmcalJet &jet)
7efbea04	52	{
	53	// Assignment operator.
	54
	55	if (this!=&jet) {
	56	AliVParticle::operator=(jet);
	57	fPt = jet.fPt;
	58	fEta = jet.fEta;
	59	fPhi = jet.fPhi;
	60	fM = jet.fM;
	61	fNEF = jet.fNEF;
	62	fArea = jet.fArea;
	63	fNch = jet.fNch;
	64	fNn = jet.fNn;
	65	fMaxCPt = jet.fMaxCPt;
	66	fMaxNPt = jet.fMaxNPt;
	67	}
	68
	69	return *this;
7df864a3	70	}
7efbea04	71
7efbea04	72	//__________________________________________________________________________________________________
914d486c	73	void AliEmcalJet::GetMom(TLorentzVector &vec) const
7efbea04	74	{
	75	// Return momentum as four vector.
	76
	77	Double_t p = fPt *TMath::CosH(fEta);
	78	vec.SetPtEtaPhiE(fPt,fEta,fPhi,TMath::Sqrt(pp+fMfM));
	79	}
	80
	81	//__________________________________________________________________________________________________
914d486c	82	void AliEmcalJet::Print(Option_t* /option/) const
7efbea04	83	{
	84	// Print jet information.
	85
	86	printf("Jet pt=%.2f, eta=%.2f, phi=%.2f, area=%.2f, NEF=%.2f\n", fPt, fEta, fPhi, fArea, fNEF);
	87	}