Added 2 methos to retrieve the primary particles
[u/mrichter/AliRoot.git] / PHOS / AliPHOSRecParticle.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
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  **************************************************************************/
15 /* $Id$ */
16 //_________________________________________________________________________
17 //  A Reconstructed Particle in PHOS    
18 //  To become a general class of AliRoot ?       
19 //  Why should I put meaningless comments
20 //  just to satisfy
21 //  the code checker                 
22 //       
23 //*-- Author: Yves Schutz (SUBATECH)
24
25
26 // --- ROOT system ---
27
28 // --- Standard library ---
29
30
31 // --- AliRoot header files ---
32
33 #include "AliPHOSRecParticle.h"
34 #include "AliPHOSGetter.h" 
35 #include "TParticle.h"
36
37 ClassImp(AliPHOSRecParticle)
38
39
40 //____________________________________________________________________________
41  AliPHOSRecParticle::AliPHOSRecParticle(const AliPHOSRecParticle & rp)
42 {
43   // copy ctor
44
45   fPHOSTrackSegment = rp.fPHOSTrackSegment ; 
46   fType             = rp.fType ; 
47   fIndexInList      = rp.fIndexInList ;
48
49   fPdgCode     = rp.fPdgCode;
50   fStatusCode  = rp.fStatusCode;
51   fMother[0]   = rp.fMother[0];
52   fMother[1]   = rp.fMother[1];
53   fDaughter[0] = rp.fDaughter[0];
54   fDaughter[1] = rp.fDaughter[1];
55   fWeight      = rp.fWeight;
56   fCalcMass    = rp.fCalcMass;
57   fPx          = rp.fPx;
58   fPy          = rp.fPy;
59   fPz          = rp.fPz;
60   fE           = rp.fE;
61   fVx          = rp.fVx;
62   fVy          = rp.fVy;
63   fVz          = rp.fVz;
64   fVt          = rp.fVt;
65   fPolarTheta  = rp.fPolarTheta;
66   fPolarPhi    = rp.fPolarPhi;
67   fParticlePDG = rp.fParticlePDG; 
68   
69 }
70
71 //____________________________________________________________________________
72 const Int_t AliPHOSRecParticle::GetNPrimaries() const  
73
74
75   Int_t rv = 0 ;
76   AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ; 
77   gime->EmcRecPoint(gime->TrackSegment(GetPHOSTSIndex())->GetEmcIndex())->GetPrimaries(rv) ; 
78   return rv ; 
79 }
80
81 //____________________________________________________________________________
82 const TParticle * AliPHOSRecParticle::GetPrimary(Int_t index) const  
83 {
84   if ( index > GetNPrimaries() ) 
85     cout << "WARNING : AliPHOSRecParticle::GetPrimary -> " << index << " is larger that the number of primaries " 
86          <<  GetNPrimaries() << endl ; 
87   else { 
88     Int_t dummy ; 
89     AliPHOSGetter * gime = AliPHOSGetter::GetInstance() ; 
90     Int_t primaryindex = gime->EmcRecPoint(gime->TrackSegment(GetPHOSTSIndex())->GetEmcIndex())->GetPrimaries(dummy)[index] ; 
91     if (primaryindex >= 10000000) 
92       cout << "WARNING : AliPHOSRecParticle::GetPrimary -> not a signal primary" << endl ; 
93     else 
94       return gime->Primary(primaryindex) ; 
95   }
96   return 0 ; 
97 }