[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 class 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 "Riostream.h"

#include <TMath.h>
#include <TString.h>

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


/// \cond CLASSIMP
ClassImp(AliMUONHit)
/// \endcond
 
//___________________________________________
AliMUONHit::AliMUONHit()
  : AliHit(), 
    fIsDetElemId(kTRUE),
    fDetElemId(0),
    fParticle(0),
    fTheta(0),
    fPhi(0),
    fTlength(0),
    fEloss(0),
    fAge(0),
    fPHfirst(0),
    fPHlast(0),
    fPTot(0),
    fPx(0),
    fPy(0),
    fPz(0),
    fXref(0),
    fYref(0),
    fZref(0)
{
/// Default constructor
}

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

    fX         = hits[1];
    fY         = hits[2];
    fZ         = hits[3];
}

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

//___________________________________________
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),
    fIsDetElemId(kFALSE),
    fDetElemId(iChamber),
    fParticle(idpart),
    fTheta(theta),
    fPhi(phi),
    fTlength(length),
    fEloss(destep),
    fAge(tof),
    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)),
    fXref(0),
    fYref(0),
    fZref(0)
{
/// Constructor
/// \deprecated TBR

    fX         = X;
    fY         = Y;
    fZ         = Z;
}

//___________________________________________
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),
    fIsDetElemId(kTRUE),
    fDetElemId(detElemId),
    fParticle(idpart),
    fTheta(theta),
    fPhi(phi),
    fTlength(length),
    fEloss(destep),
    fAge(tof),
    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)),
    fXref(0),
    fYref(0),
    fZref(0)
{
/// Constructor
    fX         = X;
    fY         = Y;
    fZ         = Z;
}

//-----------------------------------------------------------------------------------------------
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),
    fIsDetElemId(kFALSE),
    fDetElemId(iChamber),
    fParticle(idpart),
    fTheta(theta),
    fPhi(phi),
    fTlength(length),
    fEloss(destep),
    fAge(tof),
    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)),
    fXref(Xref),
    fYref(Yref),
    fZref(Zref)
{
/// Constructor
/// \deprecated TBR

    fX         = X;
    fY         = Y;
    fZ         = Z;
}
//-----------------------------------------------------------------------------------------------
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),
    fIsDetElemId(kTRUE),
    fDetElemId(detElemId),
    fParticle(idpart),
    fTheta(theta),
    fPhi(phi),
    fTlength(length),
    fEloss(destep),
    fAge(tof),
    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)),
    fXref(Xref),
    fYref(Yref),
    fZref(Zref)
{
/// Constructor

    fX         = X;
    fY         = Y;
    fZ         = Z;
}

//-----------------------------------------------------------------------------------------------
AliMUONHit::~AliMUONHit()
{
/// Dectructor
}

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

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

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

  if (!fIsDetElemId) 
    return fDetElemId;
  else  
    return AliMUONGeometryStore::GetModuleId(fDetElemId)+1;  
}

void AliMUONHit::Print(Option_t* opt) const
{
  //
  // Printing hit information 
  // "full" option for printing all the information about the hit
  //
  TString sopt(opt);
  sopt.ToUpper();
 
  if ( sopt.Contains("FULL") ) { 
    cout <<"<AliMUONHit>: Geant track="   << setw(4)  << Track() <<
      ", DetEle="        << setw(4)  << DetElemId() <<  
      ", (x,y,z)=(" << setw(7) << setprecision(5) << X() << "," << setw(7) << setprecision(5) << Y() <<  "," << setw(7) << setprecision(5) << Z() << 
      " )cm, Delta E=" << setw(8) << setprecision(3) << Eloss() << " GeV" << endl;
  }
  else {
    cout << "<AliMUONHit>: DetEle="        << setw(4)  << DetElemId() << 
      ", (x,y,z)=(" << setw(7) << setprecision(5) << X() << "," << setw(7) << setprecision(5) << Y() <<  "," << setw(7) << setprecision(5) << Z() << 
      " ) cm" <<endl;
  }
    
}
Commit	Line	Data
	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
	16	/* $Id$ */
	17
	18	// MUON class 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	#include "Riostream.h"
	27
	28	#include <TMath.h>
	29	#include <TString.h>
	30
	31	#include "AliLog.h"
	32	#include "AliMUONHit.h"
	33	#include "AliMUONGeometryStore.h"
	34
	35
	36	/// \cond CLASSIMP
	37	ClassImp(AliMUONHit)
	38	/// \endcond
	39
	40	//___________________________________________
	41	AliMUONHit::AliMUONHit()
	42	: AliHit(),
	43	fIsDetElemId(kTRUE),
	44	fDetElemId(0),
	45	fParticle(0),
	46	fTheta(0),
	47	fPhi(0),
	48	fTlength(0),
	49	fEloss(0),
	50	fAge(0),
	51	fPHfirst(0),
	52	fPHlast(0),
	53	fPTot(0),
	54	fPx(0),
	55	fPy(0),
	56	fPz(0),
	57	fXref(0),
	58	fYref(0),
	59	fZref(0)
	60	{
	61	/// Default constructor
	62	}
	63
	64	//___________________________________________
	65	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t vol, Float_t hits)
	66	: AliHit(shunt, track),
	67	fIsDetElemId(kFALSE),
	68	fDetElemId(vol[0]),
	69	fParticle(hits[0]),
	70	fTheta(hits[4]),
	71	fPhi(hits[5]),
	72	fTlength(hits[6]),
	73	fEloss(hits[7]),
	74	fAge(hits[14]),
	75	fPHfirst((Int_t)hits[8]),
	76	fPHlast((Int_t)hits[9]),
	77	fPTot(hits[10]),
	78	fPx(hits[11]),
	79	fPy(hits[12]),
	80	fPz(hits[13]),
	81	fXref(0),
	82	fYref(0),
	83	fZref(0)
	84	{
	85	/// Constructor
	86	/// \deprecated TBR
	87
	88	fX = hits[1];
	89	fY = hits[2];
	90	fZ = hits[3];
	91	}
	92
	93	//___________________________________________
	94	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t vol, Float_t hits,
	95	Bool_t /isDetElemId/)
	96	: AliHit(shunt, track),
	97	fIsDetElemId(kTRUE),
	98	fDetElemId(vol[0]),
	99	fParticle(hits[0]),
	100	fTheta(hits[4]),
	101	fPhi(hits[5]),
	102	fTlength(hits[6]),
	103	fEloss(hits[7]),
	104	fAge(hits[14]),
	105	fPHfirst((Int_t)hits[8]),
	106	fPHlast((Int_t)hits[9]),
	107	fPTot(hits[10]),
	108	fPx(hits[11]),
	109	fPy(hits[12]),
	110	fPz(hits[13]),
	111	fXref(0),
	112	fYref(0),
	113	fZref(0)
	114	{
	115	/// Constructor
	116
	117	fX = hits[1];
	118	fY = hits[2];
	119	fZ = hits[3];
	120	}
	121
	122	//___________________________________________
	123	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t iChamber, Int_t idpart,
	124	Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
	125	Float_t theta, Float_t phi, Float_t length, Float_t destep)
	126	: AliHit(shunt, track),
	127	fIsDetElemId(kFALSE),
	128	fDetElemId(iChamber),
	129	fParticle(idpart),
	130	fTheta(theta),
	131	fPhi(phi),
	132	fTlength(length),
	133	fEloss(destep),
	134	fAge(tof),
	135	fPHfirst(0),
	136	fPHlast(0),
	137	fPTot(momentum),
	138	fPx(momentum * TMath::Sin(theta) * TMath::Cos(phi)),
	139	fPy(momentum * TMath::Sin(theta) * TMath::Sin(phi)),
	140	fPz(momentum * TMath::Cos(theta)),
	141	fXref(0),
	142	fYref(0),
	143	fZref(0)
	144	{
	145	/// Constructor
	146	/// \deprecated TBR
	147
	148	fX = X;
	149	fY = Y;
	150	fZ = Z;
	151	}
	152
	153	//___________________________________________
	154	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t detElemId, Int_t idpart,
	155	Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
	156	Float_t theta, Float_t phi, Float_t length, Float_t destep,
	157	Bool_t /isDetElemId/)
	158	: AliHit(shunt, track),
	159	fIsDetElemId(kTRUE),
	160	fDetElemId(detElemId),
	161	fParticle(idpart),
	162	fTheta(theta),
	163	fPhi(phi),
	164	fTlength(length),
	165	fEloss(destep),
	166	fAge(tof),
	167	fPHfirst(0),
	168	fPHlast(0),
	169	fPTot(momentum),
	170	fPx(momentum * TMath::Sin(theta) * TMath::Cos(phi)),
	171	fPy(momentum * TMath::Sin(theta) * TMath::Sin(phi)),
	172	fPz(momentum * TMath::Cos(theta)),
	173	fXref(0),
	174	fYref(0),
	175	fZref(0)
	176	{
	177	/// Constructor
	178	fX = X;
	179	fY = Y;
	180	fZ = Z;
	181	}
	182
	183	//-----------------------------------------------------------------------------------------------
	184	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t iChamber, Int_t idpart,
	185	Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
	186	Float_t theta, Float_t phi, Float_t length, Float_t destep,
	187	Float_t Xref,Float_t Yref,Float_t Zref)
	188	: AliHit(shunt, track),
	189	fIsDetElemId(kFALSE),
	190	fDetElemId(iChamber),
	191	fParticle(idpart),
	192	fTheta(theta),
	193	fPhi(phi),
	194	fTlength(length),
	195	fEloss(destep),
	196	fAge(tof),
	197	fPHfirst(0),
	198	fPHlast(0),
	199	fPTot(momentum),
	200	fPx(momentum * TMath::Sin(theta) * TMath::Cos(phi)),
	201	fPy(momentum * TMath::Sin(theta) * TMath::Sin(phi)),
	202	fPz(momentum * TMath::Cos(theta)),
	203	fXref(Xref),
	204	fYref(Yref),
	205	fZref(Zref)
	206	{
	207	/// Constructor
	208	/// \deprecated TBR
	209
	210	fX = X;
	211	fY = Y;
	212	fZ = Z;
	213	}
	214	//-----------------------------------------------------------------------------------------------
	215	AliMUONHit::AliMUONHit(Int_t shunt, Int_t track, Int_t detElemId, Int_t idpart,
	216	Float_t X, Float_t Y, Float_t Z, Float_t tof, Float_t momentum,
	217	Float_t theta, Float_t phi, Float_t length, Float_t destep,
	218	Float_t Xref,Float_t Yref,Float_t Zref,
	219	Bool_t /isDetElemId/)
	220	: AliHit(shunt, track),
	221	fIsDetElemId(kTRUE),
	222	fDetElemId(detElemId),
	223	fParticle(idpart),
	224	fTheta(theta),
	225	fPhi(phi),
	226	fTlength(length),
	227	fEloss(destep),
	228	fAge(tof),
	229	fPHfirst(0),
	230	fPHlast(0),
	231	fPTot(momentum),
	232	fPx(momentum * TMath::Sin(theta) * TMath::Cos(phi)),
	233	fPy(momentum * TMath::Sin(theta) * TMath::Sin(phi)),
	234	fPz(momentum * TMath::Cos(theta)),
	235	fXref(Xref),
	236	fYref(Yref),
	237	fZref(Zref)
	238	{
	239	/// Constructor
	240
	241	fX = X;
	242	fY = Y;
	243	fZ = Z;
	244	}
	245
	246	//-----------------------------------------------------------------------------------------------
	247	AliMUONHit::~AliMUONHit()
	248	{
	249	/// Dectructor
	250	}
	251
	252	//-----------------------------------------------------------------------------------------------
	253	Int_t AliMUONHit::DetElemId()const
	254	{
	255	/// Return detection element ID
	256
	257	if (!fIsDetElemId) {
	258	AliWarning("Detection element Id is not defined.");
	259	return 0;
	260	}
	261
	262	return fDetElemId;
	263	}
	264
	265	//-----------------------------------------------------------------------------------------------
	266	Int_t AliMUONHit::Chamber() const
	267	{
	268	/// Return chamber ID
	269
	270	if (!fIsDetElemId)
	271	return fDetElemId;
	272	else
	273	return AliMUONGeometryStore::GetModuleId(fDetElemId)+1;
	274	}
	275
	276	void AliMUONHit::Print(Option_t* opt) const
	277	{
	278	//
	279	// Printing hit information
	280	// "full" option for printing all the information about the hit
	281	//
	282	TString sopt(opt);
	283	sopt.ToUpper();
	284
	285	if ( sopt.Contains("FULL") ) {
	286	cout <<"<AliMUONHit>: Geant track=" << setw(4) << Track() <<
	287	", DetEle=" << setw(4) << DetElemId() <<
	288	", (x,y,z)=(" << setw(7) << setprecision(5) << X() << "," << setw(7) << setprecision(5) << Y() << "," << setw(7) << setprecision(5) << Z() <<
	289	" )cm, Delta E=" << setw(8) << setprecision(3) << Eloss() << " GeV" << endl;
	290	}
	291	else {
	292	cout << "<AliMUONHit>: DetEle=" << setw(4) << DetElemId() <<
	293	", (x,y,z)=(" << setw(7) << setprecision(5) << X() << "," << setw(7) << setprecision(5) << Y() << "," << setw(7) << setprecision(5) << Z() <<
	294	" ) cm" <<endl;
	295	}
	296
	297	}