Take absolute value of daughter label in MatchToMC, to keep fakes in PbPb.

[u/mrichter/AliRoot.git] / STEER / AliAODPWG4Particle.cxx

/**************************************************************************\r
 * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *\r
 *                                                                        *\r
 * Author: The ALICE Off-line Project.                                    *\r
 * Contributors are mentioned in the code where appropriate.              *\r
 *                                                                        *\r
 * Permission to use, copy, modify and distribute this software and its   *\r
 * documentation strictly for non-commercial purposes is hereby granted   *\r
 * without fee, provided that the above copyright notice appears in all   *\r
 * copies and that both the copyright notice and this permission notice   *\r
 * appear in the supporting documentation. The authors make no claims     *\r
 * about the suitability of this software for any purpose. It is          *\r
 * provided "as is" without express or implied warranty.                  *\r
 **************************************************************************/\r
\r
/* $Id:   AliAODPWG4Particle.h $ */\r
\r
//-------------------------------------------------------------------------\r
//     AOD class for photon and other particles storage and \r
//     correlation studies\r
//     Author: Yves Schutz, CERN, Gustavo Conesa, INFN\r
//-------------------------------------------------------------------------\r
\r
//-- ROOT system --\r
\r
//-- Analysis system\r
#include "AliAODPWG4Particle.h"\r
\r
ClassImp(AliAODPWG4Particle)\r
\r
\r
//______________________________________________________________________________\r
AliAODPWG4Particle::AliAODPWG4Particle() :\r
AliVParticle(),\r
fMomentum(0),fPdg(-1), fTag(0), fBtag(-1), fLabel(-1), fCaloLabel(), fTrackLabel(),\r
fDetector(""), fDisp(0), fTof(0), fCharged(0), fTagged(0), fBadDist(0), fFidArea(0), fInputFileIndex(0)\r
{\r
  // constructor\r
  fCaloLabel [0] = -1;\r
  fCaloLabel [1] = -1;\r
  fTrackLabel[0] = -1;\r
  fTrackLabel[1] = -1;\r
}\r
\r
//______________________________________________________________________________\r
AliAODPWG4Particle::AliAODPWG4Particle(Double_t px, Double_t py, Double_t pz, Double_t e):\r
  AliVParticle(),\r
  fMomentum(0),fPdg(-1), fTag(0), fBtag(-1), fLabel(-1),fCaloLabel(), fTrackLabel(),\r
  fDetector(""), fDisp(0), fTof(0),fCharged(0), fTagged(0), fBadDist(0), fFidArea(0), fInputFileIndex(0)\r
{\r
  // constructor\r
  fMomentum = new TLorentzVector(px, py, pz, e);\r
  \r
  fCaloLabel [0] = -1;\r
  fCaloLabel [1] = -1;\r
  fTrackLabel[0] = -1;\r
  fTrackLabel[1] = -1;	\r
}\r
\r
//______________________________________________________________________________\r
AliAODPWG4Particle::AliAODPWG4Particle(TLorentzVector & p):\r
  AliVParticle(),\r
  fMomentum(0),fPdg(-1), fTag(0), fBtag(-1), fLabel(-1),fCaloLabel(), fTrackLabel(),\r
  fDetector(""), fDisp(0), fTof(0), fCharged(0), fTagged(0), fBadDist(0), fFidArea(0), fInputFileIndex(0)\r
{\r
  // constructor\r
  fMomentum = new TLorentzVector(p);\r
  \r
  fCaloLabel [0] = -1;\r
  fCaloLabel [1] = -1;\r
  fTrackLabel[0] = -1;\r
  fTrackLabel[1] = -1;\r
}\r
\r
\r
//______________________________________________________________________________\r
AliAODPWG4Particle::~AliAODPWG4Particle() \r
{\r
  // destructor\r
    delete fMomentum;\r
}\r
\r
//______________________________________________________________________________\r
AliAODPWG4Particle::AliAODPWG4Particle(const AliAODPWG4Particle& part) :\r
  AliVParticle(part),\r
  fMomentum(0), fPdg(part.fPdg), fTag(part.fTag), fBtag(part.fBtag), fLabel(part.fLabel), \r
  fCaloLabel(), fTrackLabel(), fDetector(part.fDetector),fDisp(part.fDisp), \r
  fTof(part.fTof), fCharged(part.fCharged), fTagged(part.fTagged), fBadDist(part.fBadDist), \r
  fFidArea(part.fFidArea), fInputFileIndex(part.fInputFileIndex)\r
{\r
  // Copy constructor\r
  fMomentum = new TLorentzVector(*part.fMomentum);\r
  \r
  fCaloLabel [0] = part.fCaloLabel[0];\r
  fCaloLabel [1] = part.fCaloLabel[1];\r
  fTrackLabel[0] = part.fTrackLabel[0];\r
  fTrackLabel[1] = part.fTrackLabel[1];\r
}\r
\r
//______________________________________________________________________________\r
AliAODPWG4Particle& AliAODPWG4Particle::operator=(const AliAODPWG4Particle& part)\r
{\r
  // Assignment operator\r
  if(this!=&part) {\r
    \r
    fPdg   = part.fPdg;\r
    fTag   = part.fTag;\r
    fBtag  = part.fBtag;\r
    fLabel = part.fLabel;\r
	\r
    fCaloLabel [0] = part.fCaloLabel[0];\r
    fCaloLabel [1] = part.fCaloLabel[1];\r
    fTrackLabel[0] = part.fTrackLabel[0];\r
    fTrackLabel[1] = part.fTrackLabel[1];\r
	\r
    fDetector = part.fDetector;\r
    fDisp     = part.fDisp; \r
    fTof      = part.fTof; \r
    fCharged  = part.fCharged; \r
	fTagged   = part.fTagged;\r
    fBadDist  = part.fBadDist;\r
	fFidArea  = part.fFidArea;\r
	fInputFileIndex =  part.fInputFileIndex;\r
	  \r
    if (fMomentum ) delete fMomentum;	\r
    fMomentum = new TLorentzVector(*part.fMomentum);\r
  }\r
  \r
  return *this;\r
}\r
\r
\r
//_______________________________________________________________\r
Bool_t AliAODPWG4Particle::IsPIDOK(const Int_t ipid, const Int_t pdgwanted) const{\r
  // returns true if particle satisfies given PID criterium\r
	switch(ipid){\r
	case 0: return kTRUE ; //No PID at all\r
	case 1: \r
	  {\r
	    if (fPdg == pdgwanted) return kTRUE; \r
	    else return kFALSE; //Overall PID calculated with bayesian methods.\r
	  }\r
	case 2: return fDisp ;   //only dispersion cut\r
	case 3: return fTof ;    //Only TOF cut\r
	case 4: return fCharged ;    //Only Charged cut\r
	case 5: return fDisp && fTof ;  //Dispersion and TOF\r
	case 6: return fDisp && fCharged ;  //Dispersion and Charged\r
	case 7: return fTof  && fCharged ;  //TOF and Charged\r
	case 8: return fDisp && fTof && fCharged ; // all 3 cuts\r
	default: return kFALSE ; //Not known combination\r
	}\r
}\r
\r
//______________________________________________________________________________\r
void AliAODPWG4Particle::Print(Option_t* /*option*/) const \r
{\r
  // Print information of all data members\r
  printf("Particle 4-vector:\n");\r
  printf("     E  = %13.3f", E() );\r
  printf("     Px = %13.3f", Px());\r
  printf("     Py = %13.3f", Py());\r
  printf("     Pz = %13.3f\n", Pz());\r
  printf("PID bits :\n");\r
  printf("     TOF        : %d",fTof);\r
  printf("     Charged    : %d",fCharged);\r
  printf("     Dispersion : %d\n",fDisp);\r
  printf("PDG       : %d\n",fPdg);\r
  printf("Tag       : %d\n",fTag); \r
  printf("Btag      : %d\n",fBtag);  \r
  printf("Pi0 Tag   : %d\n",fTagged);  \r
  printf("Dist. to bad channel : %d\n",fBadDist);  \r
  printf("Fid Area  : %d\n",fFidArea);  \r
  printf("Input File Index : %d\n",fInputFileIndex);  \r
  printf("Detector  : %s\n",fDetector.Data());\r
  \r
}\r