/*
$Log$
+Revision 1.12 2002/03/08 16:05:05 nilsen
+Standeard io streamers added to make debugging et al. easier.
+
+Revision 1.11 2001/01/30 09:23:13 hristov
+Streamers removed (R.Brun)
+
+Revision 1.10 2001/01/26 20:01:19 hristov
+Major upgrade of AliRoot code
+
+Revision 1.9 2000/10/02 16:32:51 barbera
+Automatic streamer used and forward declarations added
+
+Revision 1.3.4.7 2000/10/02 15:54:49 barbera
+Automatic streamer used and forward declarations added
+
+Revision 1.8 2000/09/22 12:35:21 nilsen
+Traps placed incase it is used without a properly initilized AliITSgeom class.
+
+Revision 1.7 2000/07/10 16:07:18 fca
+Release version of ITS code
+
+Revision 1.3.4.2 2000/03/04 23:43:57 nilsen
+Fixed up the comments/documentation.
+
+Revision 1.3.4.1 2000/01/12 19:03:32 nilsen
+This is the version of the files after the merging done in December 1999.
+See the ReadMe110100.txt file for details
+
+Revision 1.3 1999/09/29 09:24:20 fca
+Introduction of the Copyright and cvs Log
+
*/
+#include <iostream.h>
#include <TMath.h>
#include <TRandom.h>
#include <TGeometry.h>
#include <TNode.h>
#include <TTUBE.h>
+#include "TParticle.h"
-#include "AliITSgeom.h"
+#include "AliRun.h"
#include "AliITS.h"
+#include "AliITSgeom.h"
#include "AliITShit.h"
-#include "AliRun.h"
ClassImp(AliITShit)
<pre>
*/
//End_Html
+////////////////////////////////////////////////////////////////////////
+// Inline Member functions:
+//
+// AliITShit()
+// The default creator of the AliITShit class.
+//
+// ~AliITShit()
+// The default destructor of the AliITShit class.
+//
+// int GetTrack()
+// See AliHit for a full description. Returns the track number fTrack
+// for this hit.
+//
+// SetTrack(int track)
+// See AliHit for a full description. Sets the track number fTrack
+// for this hit.
+//
+// Int_t GetTrackStatus()
+// Returns the value of the track status flag fStatus. This flag
+// indicates the track status at the time of creating this hit. It is
+// made up of the following 8 status bits from highest order to lowest
+// order bits
+// 0 : IsTrackAlive(): IsTrackStop():IsTrackDisappeared():
+// IsTrackOut():IsTrackExiting():IsTrackEntering():IsTrackInside() .
+// See AliMC for a description of these functions. If the function is
+// true then the bit is set to one, otherwise it is zero.
+//
+// Bool_t StatusInside()
+// Returns kTRUE if the particle producing this hit is still inside
+// the present volume. Returns kFalse if this particle will be in another
+// volume. {bit IsTrackInside is set or not}
+//
+// Bool_t StatusEntering()
+// Returns kTRUE if the particle producing this hit is has just enterd
+// the present volume. Returns kFalse otherwise. {bit IsTrackEntering is
+// set or not}
+//
+// Bool_t StatusExiting()
+// Returns kTRUE if the particle producing this hit is will exit
+// the present volume. Returns kFalse otherwise. {bit IsTrackExiting is set
+// or not}
+//
+// Bool_t StatusOut()
+// Returns kTRUE if the particle producing this hit is goint exit the
+// simulation. Returns kFalse otherwise. {bit IsTrackOut is set or not}
+//
+// Bool_t StatusDisappeared()
+// Returns kTRUE if the particle producing this hit is going to "disappear"
+// for example it has interacted producing some other particles. Returns
+// kFalse otherwise. {bit IsTrackOut is set or not}
+//
+// Bool_t StatusStop()
+// Returns kTRUE if the particle producing this hit is has dropped below
+// its energy cut off producing some other particles. Returns kFalse otherwise.
+// {bit IsTrackOut is set or not}
+//
+// Bool_t StatuAlives()
+// Returns kTRUE if the particle producing this hit is going to continue
+// to be transported. Returns kFalse otherwise. {bit IsTrackOut is set or not}
+//
+// Int_t GetLayer()
+// Returns the layer number, fLayer, for this hit.
+//
+// Int_t GetLadder()
+// Returns the ladder number, fLadder, for this hit.
+//
+// Int_t GetDetector()
+// Returns the detector number, fDet, for this hit.
+//
+// GetDetectorID(Int_t &layer, Int_t &ladder, Int_t &detector)
+// Returns the layer, ladder, and detector numbers, fLayer fLadder fDet,
+// in one call.
+//
+// Float_t GetIonization()
+// Returns the energy lost, fDestep, by the particle creating this hit,
+// in the units defined by the Monte Carlo.
+//
+// GetPositionG(Float_t &x, Float_t &y, Float_t &z)
+// Returns the global position, fX fY fZ, of this hit, in the units
+// define by the Monte Carlo.
+//
+// GetPositionG(Double_t &x, Double_t &y, Double_t &z)
+// Returns the global position, fX fY fZ, of this hit, in the units
+// define by the Monte Carlo.
+//
+// GetPositionG(Float_t &x, Float_t &y, Float_t &z, Float_t &tof)
+// Returns the global position and time of flight, fX fY fZ fTof, of
+// this hit, in the units define by the Monte Carlo.
+//
+// GetPositionG(Double_t &x,Double_t &y,Double_t &z,Double_t &tof)
+// Returns the global position and time of flight, fX fY fZ fTof, of
+// this hit, in the units define by the Monte Carlo.
+//
+// GetPositionL(Double_t &x,Double_t &y,Double_t &z)
+// Returns the local position, fX fY fZ, of this hit in the coordiante
+// of this module, in the units define by the Monte Carlo.
+//
+// GetPositionG(Double_t &x,Double_t &y,Double_t &z,Double_t &tof)
+// Returns the local position and time of flight, fX fY fZ fTof, of
+// this hit in the coordinates of this module, in the units define by the
+// Monte Carlo.
+//
+// Float_t GetXG()
+// Returns the global x position in the units defined by the Monte Carlo.
+//
+// Float_t GetYG()
+// Returns the global y position in the units defined by the Monte Carlo.
+//
+// Float_t GetYG()
+// Returns the global z position in the units defined by the Monte Carlo.
+//
+// Float_t GetTOF()
+// Returns the time of flight, fTof, of this hit, in the units defined
+// by the Monte Carlo.
+//
+// GetMomentumG(Float_t &px, Float_t &py, Float_t &pz)
+// Returns the global momentum, fPx fPy fPz, of the particle that made
+// this hit, in the units define by the Monte Carlo.
+//
+// GetMomentumG(Double_t &px,Double_t &py,Double_t &pz)
+// Returns the global momentum, fPx fPy fPz, of the particle that made
+// this hit, in the units define by the Monte Carlo.
+//
+// GetMomentumL(Double_t &px,Double_t &py,Double_t &pz)
+// Returns the momentum, fPx fPy fPz in coordinate appropreate for this
+// specific module, in the units define by the Monte Carlo.
+//
+// Float_t GetPXG()
+// Returns the global X momentum in the units defined by the Monte Carlo.
+//
+// Float_t GetPYG()
+// Returns the global Y momentum in the units defined by the Monte Carlo.
+//
+// Float_t GetPZG()
+// Returns the global Z momentum in the units defined by the Monte Carlo.
+//
+////////////////////////////////////////////////////////////////////////
//_____________________________________________________________________________
AliITShit::AliITShit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits):
- AliHit(shunt, track){
- //
- // Create ITS hit
- // The creator of the AliITShit class. The variables shunt and
- // track are passed to the creator of the AliHit class. See the AliHit
- // class for a full description. the integer array *vol contains, in order,
- // fLayer = vol[0], fDet = vol[1], fLadder = vol[2], fStatus = vol[3].
- // The array *hits contains, in order, fX = hits[0], fY = hits[1],
- // fZ = hits[2], fPx = hits[3], fPy = hits[4], fPz = hits[5],
- // fDestep = hits[6], and fTof = hits[7].
- //
+ AliHit(shunt, track){
+////////////////////////////////////////////////////////////////////////
+// Create ITS hit
+// The creator of the AliITShit class. The variables shunt and
+// track are passed to the creator of the AliHit class. See the AliHit
+// class for a full description. the integer array *vol contains, in order,
+// fLayer = vol[0], fDet = vol[1], fLadder = vol[2], fStatus = vol[3].
+// The array *hits contains, in order, fX = hits[0], fY = hits[1],
+// fZ = hits[2], fPx = hits[3], fPy = hits[4], fPz = hits[5],
+// fDestep = hits[6], and fTof = hits[7]. In the units of the Monte Carlo
+////////////////////////////////////////////////////////////////////////
fLayer = vol[0]; // Layer number
fLadder = vol[2]; // Ladder number
fDet = vol[1]; // Detector number
fDestep = hits[6]; // Track dE/dx for this step
fTof = hits[7]; // Track Time of Flight for this step
}
-
+//______________________________________________________________________
void AliITShit::GetPositionL(Float_t &x,Float_t &y,Float_t &z){
+////////////////////////////////////////////////////////////////////////
+// Returns the position of this hit in the local coordinates of this
+// module, and in the units of the Monte Carlo.
+////////////////////////////////////////////////////////////////////////
AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
Float_t g[3],l[3];
g[0] = fX;
g[1] = fY;
g[2] = fZ;
- gm->GtoL(fLayer,fLadder,fDet,g,l);
- x = l[0];
- y = l[1];
- z = l[2];
+ if(gm) {
+ gm->GtoL(fLayer,fLadder,fDet,g,l);
+ x = l[0];
+ y = l[1];
+ z = l[2];
+ } else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ // AliITSv7 - SDD case
+ x=fX;
+ y=fZ;
+ z=fY;
+ }
return;
}
-
+//______________________________________________________________________
void AliITShit::GetPositionL(Float_t &x,Float_t &y,Float_t &z,Float_t &tof){
+////////////////////////////////////////////////////////////////////////
+// Returns the position and time of flight of this hit in the local
+// coordinates of this module, and in the units of the Monte Carlo.
+////////////////////////////////////////////////////////////////////////
AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
Float_t g[3],l[3];
g[0] = fX;
g[1] = fY;
g[2] = fZ;
- gm->GtoL(fLayer,fLadder,fDet,g,l);
- x = l[0];
- y = l[1];
- z = l[2];
+ if(gm) {
+ gm->GtoL(fLayer,fLadder,fDet,g,l);
+ x = l[0];
+ y = l[1];
+ z = l[2];
+ } else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ // AliITSv7 - SDD case
+ x=fX;
+ y=fZ;
+ z=fY;
+ }
tof = fTof;
return;
}
-
+//______________________________________________________________________
Float_t AliITShit::GetXL(){
+////////////////////////////////////////////////////////////////////////
+// Returns the x position of this hit in the local coordinates of this
+// module, and in the units of the Monte Carlo.
+////////////////////////////////////////////////////////////////////////
AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
Float_t g[3],l[3];
g[0] = fX;
g[1] = fY;
g[2] = fZ;
- gm->GtoL(fLayer,fLadder,fDet,g,l);
- return l[0];
+ if(gm) {
+ gm->GtoL(fLayer,fLadder,fDet,g,l);
+ return l[0];
+ } else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ return fX;
+ }
}
-
+//______________________________________________________________________
Float_t AliITShit::GetYL(){
+////////////////////////////////////////////////////////////////////////
+// Returns the y position of this hit in the local coordinates of this
+// module, and in the units of the Monte Carlo.
+////////////////////////////////////////////////////////////////////////
AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
Float_t g[3],l[3];
g[0] = fX;
g[1] = fY;
g[2] = fZ;
- gm->GtoL(fLayer,fLadder,fDet,g,l);
- return l[1];
+ if (gm) {
+ gm->GtoL(fLayer,fLadder,fDet,g,l);
+ return l[1];
+ } else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ return fZ;
+ }
}
-
+//______________________________________________________________________
Float_t AliITShit::GetZL(){
+////////////////////////////////////////////////////////////////////////
+// Returns the z position of this hit in the local coordinates of this
+// module, and in the units of the Monte Carlo.
+////////////////////////////////////////////////////////////////////////
AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
Float_t g[3],l[3];
g[0] = fX;
g[1] = fY;
g[2] = fZ;
- gm->GtoL(fLayer,fLadder,fDet,g,l);
- return l[2];
+ if(gm) {
+ gm->GtoL(fLayer,fLadder,fDet,g,l);
+ return l[2];
+ } else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ return fY;
+ }
}
-
+//______________________________________________________________________
void AliITShit::GetMomentumL(Float_t &px,Float_t &py,Float_t &pz){
+////////////////////////////////////////////////////////////////////////
+// Returns the momentum of this hit in the local coordinates of this
+// module, and in the units of the Monte Carlo.
+////////////////////////////////////////////////////////////////////////
+ AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
+ Float_t g[3],l[3];
+
+ g[0] = fPx;
+ g[1] = fPy;
+ g[2] = fPz;
+ if (gm) {
+ gm->GtoLMomentum(fLayer,fLadder,fDet,g,l);
+ px = l[0];
+ py = l[1];
+ pz = l[2];
+ } else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ px=fPx;
+ py=fPy;
+ pz=fPz;
+ }
+ return;
+}
+//______________________________________________________________________
+Float_t AliITShit::GetPXL(){
+////////////////////////////////////////////////////////////////////////
+// Returns the X momentum of this hit in the local coordinates of this
+// module, and in the units of the Monte Carlo.
+////////////////////////////////////////////////////////////////////////
+ AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
+ Float_t g[3],l[3];
+
+ g[0] = fPx;
+ g[1] = fPy;
+ g[2] = fPz;
+ if (gm) {
+ gm->GtoLMomentum(fLayer,fLadder,fDet,g,l);
+ return l[0];
+ } else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ return fPx;
+ }
+}
+//______________________________________________________________________
+Float_t AliITShit::GetPYL(){
+////////////////////////////////////////////////////////////////////////
+// Returns the Y momentum of this hit in the local coordinates of this
+// module, and in the units of the Monte Carlo.
+////////////////////////////////////////////////////////////////////////
+ AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
+ Float_t g[3],l[3];
+
+ g[0] = fPx;
+ g[1] = fPy;
+ g[2] = fPz;
+ if (gm) {
+ gm->GtoLMomentum(fLayer,fLadder,fDet,g,l);
+ return l[1];
+ } else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ return fPy;
+ }
+
+}
+//______________________________________________________________________
+Float_t AliITShit::GetPZL(){
+////////////////////////////////////////////////////////////////////////
+// Returns the Z momentum of this hit in the local coordinates of this
+// module, and in the units of the Monte Carlo.
+////////////////////////////////////////////////////////////////////////
AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
Float_t g[3],l[3];
g[0] = fPx;
g[1] = fPy;
g[2] = fPz;
- gm->GtoLMomentum(fLayer,fLadder,fDet,g,l);
- px = l[0];
- py = l[1];
- pz = l[2];
+ if (gm) {
+ gm->GtoLMomentum(fLayer,fLadder,fDet,g,l);
+ return l[2];
+ } else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ return fPz;
+ }
+
+}
+//___________________________________________________________________________;
+Int_t AliITShit::GetModule(){
+////////////////////////////////////////////////////////////////////////
+// Returns the module index number of the module where this hit was in.
+////////////////////////////////////////////////////////////////////////
+ AliITSgeom *gm = ((AliITS*)gAlice->GetDetector("ITS"))->GetITSgeom();
+
+ if (gm) return gm->GetModuleIndex(fLayer,fLadder,fDet);
+ else {
+ Error("AliITShit","NULL pointer to the geometry! return smth else",gm);
+ return 0;
+ }
+}
+//______________________________________________________________________
+TParticle * AliITShit::GetParticle(){
+////////////////////////////////////////////////////////////////////////
+// Returns the pointer to the TParticle for the particle that created
+// this hit. From the TParticle all kinds of information about this
+// particle can be found. See the TParticle class.
+////////////////////////////////////////////////////////////////////////
+ return gAlice->Particle(GetTrack());
+}
+//----------------------------------------------------------------------
+void AliITShit::Print(ostream *os){
+////////////////////////////////////////////////////////////////////////
+// Standard output format for this class.
+////////////////////////////////////////////////////////////////////////
+#if defined __GNUC__
+#if __GNUC__ > 2
+ ios::fmtflags fmt;
+#else
+ Int_t fmt;
+#endif
+#else
+ Int_t fmt;
+#endif
+
+ fmt = os->setf(ios::scientific); // set scientific floating point output
+ *os << fTrack << " " << fX << " " << fY << " " << fZ << " ";
+ fmt = os->setf(ios::hex); // set hex for fStatus only.
+ *os << fStatus << " ";
+ fmt = os->setf(ios::dec); // every thing else decimel.
+ *os << fLayer << " " << fLadder << " " << fDet << " ";;
+ *os << fPx << " " << fPy << " " << fPz << " ";
+ *os << fDestep << " " << fTof;
+ *os << endl;
+ os->flags(fmt); // reset back to old formating.
return;
}
+//----------------------------------------------------------------------
+void AliITShit::Read(istream *is){
+////////////////////////////////////////////////////////////////////////
+// Standard input format for this class.
+////////////////////////////////////////////////////////////////////////
+
+
+ *is >> fTrack >> fX >> fY >> fZ;
+ *is >> fStatus >> fLayer >> fLadder >> fDet >> fPx >> fPy >> fPz >>
+ fDestep >> fTof;
+ return;
+}
+//----------------------------------------------------------------------
+ostream &operator<<(ostream &os,AliITShit &p){
+////////////////////////////////////////////////////////////////////////
+// Standard output streaming function.
+////////////////////////////////////////////////////////////////////////
+
+ p.Print(&os);
+ return os;
+}
+//----------------------------------------------------------------------
+istream &operator>>(istream &is,AliITShit &r){
+////////////////////////////////////////////////////////////////////////
+// Standard input streaming function.
+////////////////////////////////////////////////////////////////////////
+
+ r.Read(&is);
+ return is;
+}
+//----------------------------------------------------------------------
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff