[u/mrichter/AliRoot.git] / MUON / AliMUONHit.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

// MUON classe for MonteCarlo Hits, inherited from AliHit for the 
// In addition to the ALiHit data member fX, fY, fZ and fTrack, AliMUONHit contains some info about the particle crossing the chamber:
// Impulsion: fPtot, fPx, fPy and fPz
// Reference position at the center of the chamber (wire plane) fXref, fYref and fZref
// Cumulated path along the active volume fTlength for spliting of hits for very inclined tracks 
// Energy loss of the particle inside the gas active volume.
// Incident fTheta and fPhi angle with respect of the wire plane of the chamber.
//

#include <TMath.h>

#include "AliMUONHit.h"
#include "AliMUONGeometryStore.h"
#include "AliLog.h"

ClassImp(AliMUONHit)
 
//___________________________________________
AliMUONHit::AliMUONHit()
  : AliHit() 
{
// Default constructor
}

//___________________________________________
AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits):
	AliHit(shunt, track)
{
// Constructor
// TBR
   
    fIsDetElemId = kFALSE;
    fDetElemId = vol[0];
    fParticle  = hits[0];
    fX         = hits[1];
    fY         = hits[2];
    fZ         = hits[3];
    fTheta     = hits[4];
    fPhi       = hits[5];
    fTlength   = hits[6];
    fEloss     = hits[7];
    fPHfirst   = (Int_t) hits[8];
    fPHlast    = (Int_t) hits[9];
    fPTot      = hits[10];
    fPx        = hits[11];
    fPy        = hits[12];
    fPz        = hits[13];
    fAge       = hits[14];
    fXref      = 0.;
    fYref      = 0.;
    fZref      = 0.;
}

//___________________________________________
AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits, 
                       Bool_t /*isDetElemId*/) :
	AliHit(shunt, track)
{
// Constructor
   
    fIsDetElemId = kTRUE;
    fDetElemId = vol[0];
    fParticle  = hits[0];
    fX         = hits[1];
    fY         = hits[2];
    fZ         = hits[3];
    fTheta     = hits[4];
    fPhi       = hits[5];
    fTlength   = hits[6];
    fEloss     = hits[7];
    fPHfirst   = (Int_t) hits[8];
    fPHlast    = (Int_t) hits[9];
    fPTot      = hits[10];
    fPx        = hits[11];
    fPy        = hits[12];
    fPz        = hits[13];
    fAge       = hits[14];
    fXref      = 0.;
    fYref      = 0.;
    fZref      = 0.;
}

//___________________________________________
AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t iChamber, Int_t idpart, 
		       Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum, 
		       Float_t theta, Float_t phi, Float_t length, Float_t destep):
	AliHit(shunt, track)
{
// Constructor
// TBR

    fIsDetElemId = kFALSE;
    fDetElemId = iChamber;
    fParticle  = idpart;
    fX         = X;
    fY         = Y;
    fZ         = Z;
    fTheta     = theta;
    fPhi       = phi;
    fTlength   = length;
    fEloss     = destep;
    fPHfirst   = 0;
    fPHlast    = 0;
    fPTot      = momentum;
    fPx        = momentum * TMath::Sin(theta) * TMath::Cos(phi);
    fPy        = momentum * TMath::Sin(theta) * TMath::Sin(phi);
    fPz        = momentum * TMath::Cos(theta) ;
    fAge       = tof;
    fXref      = 0.;
    fYref      = 0.;
    fZref      = 0.;
}

//___________________________________________
AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t detElemId, Int_t idpart, 
		       Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum, 
		       Float_t theta, Float_t phi, Float_t length, Float_t destep,
		       Bool_t /*isDetElemId*/):
	AliHit(shunt, track)
{
// Constructor
    fIsDetElemId = kTRUE;
    fDetElemId = detElemId;
    fParticle  = idpart;
    fX         = X;
    fY         = Y;
    fZ         = Z;
    fTheta     = theta;
    fPhi       = phi;
    fTlength   = length;
    fEloss     = destep;
    fPHfirst   = 0;
    fPHlast    = 0;
    fPTot      = momentum;
    fPx        = momentum * TMath::Sin(theta) * TMath::Cos(phi);
    fPy        = momentum * TMath::Sin(theta) * TMath::Sin(phi);
    fPz        = momentum * TMath::Cos(theta) ;
    fAge       = tof;
    fXref      = 0.;
    fYref      = 0.;
    fZref      = 0.;
}

//-----------------------------------------------------------------------------------------------
AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t iChamber, Int_t idpart, 
		       Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum, 
		       Float_t theta, Float_t phi, Float_t length, Float_t destep,
		       Float_t Xref,Float_t Yref,Float_t Zref):
	AliHit(shunt, track)
{
// Constructor
// TBR

    fIsDetElemId = kFALSE;
    fDetElemId = iChamber;
    fParticle  = idpart;
    fX         = X;
    fY         = Y;
    fZ         = Z;
    fTheta     = theta;
    fPhi       = phi;
    fTlength   = length;
    fEloss     = destep;
    fPHfirst   = 0;
    fPHlast    = 0;
    fPTot      = momentum;
    fPx        = momentum * TMath::Sin(theta) * TMath::Cos(phi);
    fPy        = momentum * TMath::Sin(theta) * TMath::Sin(phi);
    fPz        = momentum * TMath::Cos(theta) ;
    fAge       = tof;
    fXref      = Xref;
    fYref      = Yref;
    fZref      = Zref;
}
//-----------------------------------------------------------------------------------------------
AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t detElemId, Int_t idpart, 
		       Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum, 
		       Float_t theta, Float_t phi, Float_t length, Float_t destep,
		       Float_t Xref,Float_t Yref,Float_t Zref,
		       Bool_t /*isDetElemId*/):
	AliHit(shunt, track)
{
// Constructor
    fIsDetElemId = kTRUE;
    fDetElemId = detElemId;
    fParticle  = idpart;
    fX         = X;
    fY         = Y;
    fZ         = Z;
    fTheta     = theta;
    fPhi       = phi;
    fTlength   = length;
    fEloss     = destep;
    fPHfirst   = 0;
    fPHlast    = 0;
    fPTot      = momentum;
    fPx        = momentum * TMath::Sin(theta) * TMath::Cos(phi);
    fPy        = momentum * TMath::Sin(theta) * TMath::Sin(phi);
    fPz        = momentum * TMath::Cos(theta) ;
    fAge       = tof;
    fXref      = Xref;
    fYref      = Yref;
    fZref      = Zref;
}

//-----------------------------------------------------------------------------------------------
Int_t AliMUONHit::DetElemId()const
{
// Return detection element ID

  if (!fIsDetElemId) {
    AliWarning("Detection element Id is not defined.");
    return 0;
  }  
  // end of TBR
  
  return fDetElemId;
}

//-----------------------------------------------------------------------------------------------
Int_t  AliMUONHit::Chamber()  const
{  
// Return chamber ID

  if (!fIsDetElemId) 
    return fDetElemId;
  else  
    return AliMUONGeometryStore::GetModuleId(fDetElemId)+1;  
}
Commit	Line	Data
a9e2aefa	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
88cb7938	16	/* $Id$ */
0e83b872	17
34c512f6	18	// MUON classe for MonteCarlo Hits, inherited from AliHit for the
	19	// In addition to the ALiHit data member fX, fY, fZ and fTrack, AliMUONHit contains some info about the particle crossing the chamber:
	20	// Impulsion: fPtot, fPx, fPy and fPz
	21	// Reference position at the center of the chamber (wire plane) fXref, fYref and fZref
	22	// Cumulated path along the active volume fTlength for spliting of hits for very inclined tracks
	23	// Energy loss of the particle inside the gas active volume.
	24	// Incident fTheta and fPhi angle with respect of the wire plane of the chamber.
	25	//
	26
30178c30	27	#include <TMath.h>
30178c30	28
a9e2aefa	29	#include "AliMUONHit.h"
e5df8ea1	30	#include "AliMUONGeometryStore.h"
d59b2ec3	31	#include "AliLog.h"
a9e2aefa	32
	33	ClassImp(AliMUONHit)
	34
30178c30	35	//___________________________________________
	36	AliMUONHit::AliMUONHit()
	37	: AliHit()
	38	{
	39	// Default constructor
	40	}
	41
a9e2aefa	42	//___________________________________________
0e83b872	43	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t vol, Float_t hits):
a9e2aefa	44	AliHit(shunt, track)
	45	{
	46	// Constructor
d59b2ec3	47	// TBR
	48
	49	fIsDetElemId = kFALSE;
	50	fDetElemId = vol[0];
2eb55fab	51	fParticle = hits[0];
	52	fX = hits[1];
	53	fY = hits[2];
	54	fZ = hits[3];
	55	fTheta = hits[4];
	56	fPhi = hits[5];
	57	fTlength = hits[6];
	58	fEloss = hits[7];
	59	fPHfirst = (Int_t) hits[8];
	60	fPHlast = (Int_t) hits[9];
	61	fPTot = hits[10];
	62	fPx = hits[11];
	63	fPy = hits[12];
	64	fPz = hits[13];
	65	fAge = hits[14];
1391e633	66	fXref = 0.;
	67	fYref = 0.;
	68	fZref = 0.;
a9e2aefa	69	}
d59b2ec3	70
	71	//___________________________________________
	72	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t vol, Float_t hits,
	73	Bool_t /isDetElemId/) :
	74	AliHit(shunt, track)
	75	{
	76	// Constructor
	77
	78	fIsDetElemId = kTRUE;
	79	fDetElemId = vol[0];
	80	fParticle = hits[0];
	81	fX = hits[1];
	82	fY = hits[2];
	83	fZ = hits[3];
	84	fTheta = hits[4];
	85	fPhi = hits[5];
	86	fTlength = hits[6];
	87	fEloss = hits[7];
	88	fPHfirst = (Int_t) hits[8];
	89	fPHlast = (Int_t) hits[9];
	90	fPTot = hits[10];
	91	fPx = hits[11];
	92	fPy = hits[12];
	93	fPz = hits[13];
	94	fAge = hits[14];
	95	fXref = 0.;
	96	fYref = 0.;
	97	fZref = 0.;
	98	}
	99
0e83b872	100	//___________________________________________
1391e633	101	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t iChamber, Int_t idpart,
	102	Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
	103	Float_t theta, Float_t phi, Float_t length, Float_t destep):
0e83b872	104	AliHit(shunt, track)
	105	{
	106	// Constructor
d59b2ec3	107	// TBR
	108
	109	fIsDetElemId = kFALSE;
	110	fDetElemId = iChamber;
0e83b872	111	fParticle = idpart;
	112	fX = X;
	113	fY = Y;
	114	fZ = Z;
	115	fTheta = theta;
	116	fPhi = phi;
	117	fTlength = length;
	118	fEloss = destep;
	119	fPHfirst = 0;
	120	fPHlast = 0;
	121	fPTot = momentum;
c55ca9e5	122	fPx = momentum * TMath::Sin(theta) * TMath::Cos(phi);
c55ca9e5	123	fPy = momentum * TMath::Sin(theta) * TMath::Sin(phi);
e516b01d	124	fPz = momentum * TMath::Cos(theta) ;
0e83b872	125	fAge = tof;
1391e633	126	fXref = 0.;
	127	fYref = 0.;
	128	fZref = 0.;
0e83b872	129	}
d59b2ec3	130
	131	//___________________________________________
	132	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t detElemId, Int_t idpart,
	133	Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
	134	Float_t theta, Float_t phi, Float_t length, Float_t destep,
	135	Bool_t /isDetElemId/):
	136	AliHit(shunt, track)
	137	{
	138	// Constructor
	139	fIsDetElemId = kTRUE;
	140	fDetElemId = detElemId;
	141	fParticle = idpart;
	142	fX = X;
	143	fY = Y;
	144	fZ = Z;
	145	fTheta = theta;
	146	fPhi = phi;
	147	fTlength = length;
	148	fEloss = destep;
	149	fPHfirst = 0;
	150	fPHlast = 0;
	151	fPTot = momentum;
	152	fPx = momentum * TMath::Sin(theta) * TMath::Cos(phi);
	153	fPy = momentum * TMath::Sin(theta) * TMath::Sin(phi);
	154	fPz = momentum * TMath::Cos(theta) ;
	155	fAge = tof;
	156	fXref = 0.;
	157	fYref = 0.;
	158	fZref = 0.;
	159	}
	160
1391e633	161	//-----------------------------------------------------------------------------------------------
	162	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t iChamber, Int_t idpart,
	163	Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
	164	Float_t theta, Float_t phi, Float_t length, Float_t destep,
	165	Float_t Xref,Float_t Yref,Float_t Zref):
	166	AliHit(shunt, track)
	167	{
	168	// Constructor
d59b2ec3	169	// TBR
	170
	171	fIsDetElemId = kFALSE;
	172	fDetElemId = iChamber;
1391e633	173	fParticle = idpart;
	174	fX = X;
	175	fY = Y;
	176	fZ = Z;
	177	fTheta = theta;
	178	fPhi = phi;
	179	fTlength = length;
	180	fEloss = destep;
	181	fPHfirst = 0;
	182	fPHlast = 0;
	183	fPTot = momentum;
c55ca9e5	184	fPx = momentum * TMath::Sin(theta) * TMath::Cos(phi);
c55ca9e5	185	fPy = momentum * TMath::Sin(theta) * TMath::Sin(phi);
e516b01d	186	fPz = momentum * TMath::Cos(theta) ;
1391e633	187	fAge = tof;
	188	fXref = Xref;
	189	fYref = Yref;
	190	fZref = Zref;
	191	}
	192	//-----------------------------------------------------------------------------------------------
d59b2ec3	193	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t detElemId, Int_t idpart,
	194	Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
	195	Float_t theta, Float_t phi, Float_t length, Float_t destep,
	196	Float_t Xref,Float_t Yref,Float_t Zref,
	197	Bool_t /isDetElemId/):
	198	AliHit(shunt, track)
	199	{
	200	// Constructor
	201	fIsDetElemId = kTRUE;
	202	fDetElemId = detElemId;
	203	fParticle = idpart;
	204	fX = X;
	205	fY = Y;
	206	fZ = Z;
	207	fTheta = theta;
	208	fPhi = phi;
	209	fTlength = length;
	210	fEloss = destep;
	211	fPHfirst = 0;
	212	fPHlast = 0;
	213	fPTot = momentum;
	214	fPx = momentum * TMath::Sin(theta) * TMath::Cos(phi);
	215	fPy = momentum * TMath::Sin(theta) * TMath::Sin(phi);
	216	fPz = momentum * TMath::Cos(theta) ;
	217	fAge = tof;
	218	fXref = Xref;
	219	fYref = Yref;
	220	fZref = Zref;
	221	}
	222
	223	//-----------------------------------------------------------------------------------------------
	224	Int_t AliMUONHit::DetElemId()const
	225	{
	226	// Return detection element ID
	227
	228	if (!fIsDetElemId) {
	229	AliWarning("Detection element Id is not defined.");
	230	return 0;
	231	}
	232	// end of TBR
	233
	234	return fDetElemId;
	235	}
	236
	237	//-----------------------------------------------------------------------------------------------
	238	Int_t AliMUONHit::Chamber() const
	239	{
	240	// Return chamber ID
	241
8690a0a6	242	if (!fIsDetElemId)
d59b2ec3	243	return fDetElemId;
d59b2ec3	244	else
e5df8ea1	245	return AliMUONGeometryStore::GetModuleId(fDetElemId)+1;
d59b2ec3	246	}
a9e2aefa	247