1 /**************************************************************************
2 * Copyright(c) 1998-1999, 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 **************************************************************************/
18 //_________________________________________________________________________
19 // Hits class for EMCAL
20 // A hit in EMCAL is the sum of all hits in a single segment
21 // from a single enterring particle
22 //*-- Author: Sahal Yacoob (LBL / UCT)
23 // Based on AliPHOSHit
25 // --- Standard library ---
29 #include <Rstrstream.h>
30 #include <Riostream.h>
32 // --- ROOT system ---
33 #include <TLorentzVector.h>
35 // --- AliRoot header files ---
36 #include "AliEMCALHit.h"
39 #include "AliEMCALGeometry.h"
40 #include "AliEMCALGetter.h"
44 //______________________________________________________________________
45 AliEMCALHit::AliEMCALHit(){
63 //______________________________________________________________________
64 AliEMCALHit::AliEMCALHit(const AliEMCALHit & hit){
69 fPrimary = hit.fPrimary ;
78 fIparent = hit.fIparent;
79 fIenergy = hit.fIenergy;
82 //______________________________________________________________________
83 AliEMCALHit::AliEMCALHit(Int_t shunt, Int_t primary, Int_t track,Int_t iparent, Float_t ienergy, Int_t id,
84 Float_t *hits,Float_t *p):AliHit(shunt, track){
86 // Create an EMCAL hit object
103 //______________________________________________________________________
104 Bool_t AliEMCALHit::operator==(AliEMCALHit const &rValue) const{
105 // Two hits are identical if they have the same Id and originat
106 // from the same enterring Particle
109 if ( (fId == rValue.GetId()) && ( fIparent == rValue.GetIparent()) )
114 //______________________________________________________________________
115 AliEMCALHit AliEMCALHit::operator+(const AliEMCALHit &rValue){
116 // Add the energy of the hit
118 fELOS += rValue.GetEnergy() ;
120 if(rValue.GetTime() < fTime)
121 fTime = rValue.GetTime() ;
126 //______________________________________________________________________
127 ostream& operator << (ostream& out,AliEMCALHit& hit){
128 // Print out Id and energy
130 out << "AliEMCALHit:";
131 out << "id=" << hit.GetId();
132 out << ", Eloss=" << hit.GetEnergy();
133 out << ", Time=" << hit.GetTime();
134 out << "GeV , Track no.=" << hit.GetPrimary();
135 out << ", (xyz)=(" << hit.X()<< ","<< hit.Y()<< ","<<hit.Z()<<") cm";
136 out << ", fTrack=" << hit.GetTrack();
137 out << ", P=(" << hit.GetPx() << "," << hit.GetPy() << "," << hit.GetPz()
138 << "," <<hit.GetPe() << ") GeV" ;
139 out << ", Enterring particle ID" << hit.GetIparent();
140 out << ", Enterring particle initial energy = " << hit.GetIenergy() << " GeV" ;