1 /**************************************************************************
2 * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
17 //-------------------------------------------------------------------------
18 // Realisation of AliVParticle for MC Particles
19 // Basically a stripped down AliMCParicle / TParticle
20 // with minimum information on MC tracks
21 // Author: Christian Klein-Bösing, CERN
22 //-------------------------------------------------------------------------
25 #include "AliAODMCParticle.h"
26 #include "AliAODEvent.h"
28 #include "TDatabasePDG.h"
29 #include "TParticle.h"
30 #include "TClonesArray.h"
33 ClassImp(AliAODMCParticle)
35 TString AliAODMCParticle::fgkStdBranchName("mcparticles");
37 AliAODMCParticle::AliAODMCParticle():
52 // Default Constructor
53 fDaughter[0] = fDaughter[1] = 0;
57 AliAODMCParticle::AliAODMCParticle(AliMCParticle* mcpart, Int_t label,Int_t flag):
59 fPdgCode(mcpart->Particle()->GetPdgCode()),
62 fMother(mcpart->GetMother()),
63 fPx(mcpart->Particle()->Px()),
64 fPy(mcpart->Particle()->Py()),
65 fPz(mcpart->Particle()->Pz()),
66 fE(mcpart->Particle()->Energy()),
67 fVx(mcpart->Particle()->Vx()),
68 fVy(mcpart->Particle()->Vy()),
69 fVz(mcpart->Particle()->Vz()),
70 fVt(mcpart->Particle()->T())
72 fDaughter[0] = mcpart->GetFirstDaughter();
73 fDaughter[1] = mcpart->GetLastDaughter();
77 AliAODMCParticle::AliAODMCParticle(const AliAODMCParticle& mcPart) :
79 fPdgCode(mcPart.fPdgCode),
81 fLabel(mcPart.fLabel),
82 fMother(mcPart.fMother),
93 fDaughter[0] = mcPart.fDaughter[0];
94 fDaughter[1] = mcPart.fDaughter[1];
97 AliAODMCParticle& AliAODMCParticle::operator=(const AliAODMCParticle& mcPart)
101 AliVParticle::operator=(mcPart);
102 fPdgCode = mcPart.fPdgCode;
103 fFlag = mcPart.fFlag;
104 fLabel = mcPart.fLabel;
105 fMother = mcPart.fMother;
114 fDaughter[0] = mcPart.fDaughter[0];
115 fDaughter[1] = mcPart.fDaughter[1];
122 Double_t AliAODMCParticle::M() const
124 TParticlePDG* pdg = TDatabasePDG::Instance()->GetParticle(fPdgCode);
126 return (pdg->Mass());
128 return GetCalcMass();
133 Short_t AliAODMCParticle::Charge() const
135 TParticlePDG* pdg = TDatabasePDG::Instance()->GetParticle(fPdgCode);
137 return (Short_t (pdg->Charge()));
143 void AliAODMCParticle::Print(const Option_t */*opt*/) const {
144 if(TDatabasePDG::Instance()->GetParticle(fPdgCode)){
145 Printf(">>> PDG (%d) : %s",fPdgCode,TDatabasePDG::Instance()->GetParticle(fPdgCode)->GetName());
148 Printf(">>> PDG (%d) : %s",fPdgCode,"Unknown");
150 Printf(">> P(%3.3f,%3.3f,%3.3f) V((%3.3f,%3.3f,%3.3f,%3.3f)",fPx,fPy,fPz,fVx,fVy,fVz,fVt);
151 Printf("> Mother %d, First Daughter %d Last Daughter %d , Status %d, PhysicalPrimary %d",
152 fMother,fDaughter[0],fDaughter[1],GetStatus(),
153 IsPhysicalPrimary());