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

// $Id$
//
// Class AliMUONChamberGeometry
// -----------------------------
// Class for definititon of the MUON chamber positions in ALIC.
//
// Author: Ivana Hrivnacova, IPN Orsay

#include <TVirtualMC.h>
#include <TGeoMatrix.h>
#include <TObjArray.h>
#include <TArrayI.h>
#include <Riostream.h>

#include "AliMUONChamberGeometry.h"
#include "AliMUONGeometryEnvelope.h"

ClassImp(AliMUONChamberGeometry)

//______________________________________________________________________________
AliMUONChamberGeometry::AliMUONChamberGeometry(Int_t chamberId)
 : TObject(),
   fChamberId(chamberId),
   fMotherVolume("ALIC"),
   fTransformation(0),
   fEnvelopes(0),
   fNofSensVolumeIds(0),
   fSensVolumeIds(0),
   fDebug(kFALSE)
{
// Standard constructor

  // Create the chamber transformation
  fTransformation = new TGeoCombiTrans("");
  fEnvelopes = new TObjArray(20);
  fSensVolumeIds = new TArrayI(20);
}


//______________________________________________________________________________
AliMUONChamberGeometry::AliMUONChamberGeometry()
 : TObject(),
   fChamberId(0),
   fMotherVolume(),
   fTransformation(0),
   fEnvelopes(0),
   fNofSensVolumeIds(0),
   fSensVolumeIds(0),
   fDebug(kFALSE)
{
// Default constructor
}


//______________________________________________________________________________
AliMUONChamberGeometry::AliMUONChamberGeometry(const AliMUONChamberGeometry& rhs)
  : TObject(rhs)
{
  Fatal("Copy constructor", 
        "Copy constructor is not implemented.");
}

//______________________________________________________________________________
AliMUONChamberGeometry::~AliMUONChamberGeometry() {
//

  // Add deleting rotation matrices 
  
  delete fTransformation;
  
  if (fEnvelopes) {
    fEnvelopes->Delete();
    delete fEnvelopes;
  }  
}

//______________________________________________________________________________
AliMUONChamberGeometry& 
AliMUONChamberGeometry::operator = (const AliMUONChamberGeometry& rhs) 
{
  // check assignement to self
  if (this == &rhs) return *this;

  Fatal("operator=", 
        "Assignment operator is not implemented.");
    
  return *this;  
}

//
// private methods
//

//______________________________________________________________________________
AliMUONGeometryEnvelope* 
AliMUONChamberGeometry::FindEnvelope(const TString& name) const
{
// Finds the envelope specified by name.
// ---

  for (Int_t i=0; i<fEnvelopes->GetEntriesFast(); i++) {
    AliMUONGeometryEnvelope* envelope 
      = (AliMUONGeometryEnvelope*)fEnvelopes->At(i);
    
    if (envelope->GetName() == name) return envelope;
  }
  
  return 0;    
}  

//
// public methods
//

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual) 
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding ";
    if (!isVirtual) cout << " non-";
    cout << "virtual envelope " << name << endl;
  }  

  AliMUONGeometryEnvelope* envelope 
    = new AliMUONGeometryEnvelope(name, isVirtual);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual,
                                          const TGeoTranslation& translation)
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding ";
    if (!isVirtual) cout << " non-";
    cout << "virtual envelope " << name 
         << " with translation" << endl;
  }  

  AliMUONGeometryEnvelope* envelope 
    = new AliMUONGeometryEnvelope(name, isVirtual);
  envelope->SetTranslation(translation);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual, 
                                          const TGeoTranslation& translation,
		                          const TGeoRotation& rotation)
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding ";
    if (!isVirtual) cout << " non-";
    cout << "virtual envelope " << name 
         << " with translation and rotation" << endl;
  }  

  // fEnvelopes->Add(new TGeoCombiTrans(name, translation, rotation));
           // would be nice to be so simple 

  AliMUONGeometryEnvelope* envelope 
    = new AliMUONGeometryEnvelope(name, isVirtual);
  envelope->SetRotation(rotation);
  envelope->SetTranslation(translation);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo) 
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding "
         << " non-virtual envelope " << name 
         << " with copyNo " << copyNo << endl;
   }  

  AliMUONGeometryEnvelope* envelope 
    = new AliMUONGeometryEnvelope(name, copyNo);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo,
                                          const TGeoTranslation& translation)
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding "
         << " non-virtual envelope " << name 
         << " with copyNo " << copyNo
         << " with translation " << endl;
  }  

  AliMUONGeometryEnvelope* envelope 
    = new AliMUONGeometryEnvelope(name, copyNo);
  envelope->SetTranslation(translation);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo, 
                                          const TGeoTranslation& translation,
		                          const TGeoRotation& rotation)
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding "
         << " non-virtual envelope " << name 
         << " with copyNo " << copyNo
         << " with translation and rotation" << endl;
  }  

  // fEnvelopes->Add(new TGeoCombiTrans(name, translation, rotation));
           // would be nice to be so simple 

  AliMUONGeometryEnvelope* envelope 
    = new AliMUONGeometryEnvelope(name, copyNo);
  envelope->SetRotation(rotation);
  envelope->SetTranslation(translation);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelopeConstituent(const TString& name, 
                                         const TString& envName, Int_t copyNo) 
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding constituent " << name
         << " to envelope " << envName 
         << " with copyNo " << copyNo << endl;
  }  

  AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
  
  if (!envelope) {
    // add warning
    return;
  }  
   
  envelope->AddConstituent(name, copyNo);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelopeConstituent(const TString& name, 
                                          const TString& envName, Int_t copyNo,
                                          const TGeoTranslation& translation)
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding constituent " << name
         << " to envelope " << envName 
         << " with copyNo " << copyNo
         << " with translation" << endl;
  }  

  AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
  
  if (!envelope) {
    // add warning
    return;
  }  
   
  envelope->AddConstituent(name, copyNo, translation);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelopeConstituent(const TString& name, 
                                          const TString& envName, Int_t copyNo, 
                                          const TGeoTranslation& translation,
		                          const TGeoRotation& rotation)
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding constituent " << name
         << " to envelope " << envName 
         << " with copyNo " << copyNo
         << " with translation and rotation" << endl;
  }  

  AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
  
  if (!envelope) {
    // add warning
    return;
  }  
   
  envelope->AddConstituent(name, copyNo, translation, rotation);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelopeConstituentParam(const TString& name, 
                                         const TString& envName, Int_t copyNo,
					 Int_t npar, Double_t* param) 
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding parameterised constituent " << name
         << " to envelope " << envName 
         << " with copyNo " << copyNo << endl;
  }  

  AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
  
  if (!envelope) {
    // add warning
    return;
  }  
   
  envelope->AddConstituentParam(name, copyNo, npar, param);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelopeConstituentParam(const TString& name, 
                                          const TString& envName, Int_t copyNo,
                                          const TGeoTranslation& translation,
					  Int_t npar, Double_t* param)
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding parameterised constituent " << name
         << " to envelope " << envName 
         << " with copyNo " << copyNo
         << " with translation" << endl;
  }  

  AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
  
  if (!envelope) {
    // add warning
    return;
  }  
   
  envelope->AddConstituentParam(name, copyNo, translation, npar, param);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelopeConstituentParam(const TString& name, 
                                          const TString& envName, Int_t copyNo, 
                                          const TGeoTranslation& translation,
		                          const TGeoRotation& rotation,
					  Int_t npar, Double_t* param)
{
// Adds the volume with the specified name and transformation
// to the list of envelopes.
// ---  					   

  if (fDebug) {
    cout << "... Adding parameterised constituent " << name
         << " to envelope " << envName 
         << " with copyNo " << copyNo
         << " with translation and rotation" << endl;
  }  

  AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
  
  if (!envelope) {
    // add warning
    return;
  }  
   
  envelope->AddConstituentParam(name, copyNo, translation, rotation, npar, param);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::SetTranslation(const TGeoTranslation& translation)
{
// Sets the chamber position wrt ALIC.
// ---

  fTransformation
    ->SetTranslation(const_cast<Double_t*>(translation.GetTranslation()));
}  

//______________________________________________________________________________
void  AliMUONChamberGeometry::SetRotation(const TGeoRotation& rotation)
{
// Sets the chamber rotation wrt ALIC.
// ---

  TGeoRotation* rot = new TGeoRotation();
  rot->SetMatrix(const_cast<Double_t*>(rotation.GetRotationMatrix()));

  fTransformation->SetRotation(rot);
}  

//______________________________________________________________________________
void  AliMUONChamberGeometry::SetSensitiveVolume(Int_t volId)
{
// Adds the volume specified by volId to the list of sensitive
// volumes
  
  fSensVolumeIds->AddAt(volId,fNofSensVolumeIds++);
}      

//______________________________________________________________________________
void  AliMUONChamberGeometry::SetSensitiveVolume(const TString& volName)
{
// Adds the volume specified by volId to the list of sensitive
// volumes

  fSensVolumeIds->AddAt(gMC->VolId(volName),fNofSensVolumeIds++);
}      

//______________________________________________________________________________
Bool_t AliMUONChamberGeometry::IsSensitiveVolume(Int_t volId) const
{
// Checks if the volume specified by volId is present in the list
// of sensitive volumes.

  for (Int_t i=0; i<fNofSensVolumeIds; i++) {
      if (fSensVolumeIds->At(i) == volId) return kTRUE;
  }
  return kFALSE;
}
Commit	Line	Data
	1	// $Id$
	2	//
	3	// Class AliMUONChamberGeometry
	4	// -----------------------------
	5	// Class for definititon of the MUON chamber positions in ALIC.
	6	//
	7	// Author: Ivana Hrivnacova, IPN Orsay
	8
	9	#include <TVirtualMC.h>
	10	#include <TGeoMatrix.h>
	11	#include <TObjArray.h>
	12	#include <TArrayI.h>
	13	#include <Riostream.h>
	14
	15	#include "AliMUONChamberGeometry.h"
	16	#include "AliMUONGeometryEnvelope.h"
	17
	18	ClassImp(AliMUONChamberGeometry)
	19
	20	//______________________________________________________________________________
	21	AliMUONChamberGeometry::AliMUONChamberGeometry(Int_t chamberId)
	22	: TObject(),
	23	fChamberId(chamberId),
	24	fMotherVolume("ALIC"),
	25	fTransformation(0),
	26	fEnvelopes(0),
	27	fNofSensVolumeIds(0),
	28	fSensVolumeIds(0),
	29	fDebug(kFALSE)
	30	{
	31	// Standard constructor
	32
	33	// Create the chamber transformation
	34	fTransformation = new TGeoCombiTrans("");
	35	fEnvelopes = new TObjArray(20);
	36	fSensVolumeIds = new TArrayI(20);
	37	}
	38
	39
	40	//______________________________________________________________________________
	41	AliMUONChamberGeometry::AliMUONChamberGeometry()
	42	: TObject(),
	43	fChamberId(0),
	44	fMotherVolume(),
	45	fTransformation(0),
	46	fEnvelopes(0),
	47	fNofSensVolumeIds(0),
	48	fSensVolumeIds(0),
	49	fDebug(kFALSE)
	50	{
	51	// Default constructor
	52	}
	53
	54
	55	//______________________________________________________________________________
	56	AliMUONChamberGeometry::AliMUONChamberGeometry(const AliMUONChamberGeometry& rhs)
	57	: TObject(rhs)
	58	{
	59	Fatal("Copy constructor",
	60	"Copy constructor is not implemented.");
	61	}
	62
	63	//______________________________________________________________________________
	64	AliMUONChamberGeometry::~AliMUONChamberGeometry() {
	65	//
	66
	67	// Add deleting rotation matrices
	68
	69	delete fTransformation;
	70
	71	if (fEnvelopes) {
	72	fEnvelopes->Delete();
	73	delete fEnvelopes;
	74	}
	75	}
	76
	77	//______________________________________________________________________________
	78	AliMUONChamberGeometry&
	79	AliMUONChamberGeometry::operator = (const AliMUONChamberGeometry& rhs)
	80	{
	81	// check assignement to self
	82	if (this == &rhs) return *this;
	83
	84	Fatal("operator=",
	85	"Assignment operator is not implemented.");
	86
	87	return *this;
	88	}
	89
	90	//
	91	// private methods
	92	//
	93
	94	//______________________________________________________________________________
	95	AliMUONGeometryEnvelope*
	96	AliMUONChamberGeometry::FindEnvelope(const TString& name) const
	97	{
	98	// Finds the envelope specified by name.
	99	// ---
	100
	101	for (Int_t i=0; i<fEnvelopes->GetEntriesFast(); i++) {
	102	AliMUONGeometryEnvelope* envelope
	103	= (AliMUONGeometryEnvelope*)fEnvelopes->At(i);
	104
	105	if (envelope->GetName() == name) return envelope;
	106	}
	107
	108	return 0;
	109	}
	110
	111	//
	112	// public methods
	113	//
	114
	115	//______________________________________________________________________________
	116	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual)
	117	{
	118	// Adds the volume with the specified name and transformation
	119	// to the list of envelopes.
	120	// ---
	121
	122	if (fDebug) {
	123	cout << "... Adding ";
	124	if (!isVirtual) cout << " non-";
	125	cout << "virtual envelope " << name << endl;
	126	}
	127
	128	AliMUONGeometryEnvelope* envelope
	129	= new AliMUONGeometryEnvelope(name, isVirtual);
	130
	131	fEnvelopes->Add(envelope);
	132	}
	133
	134	//______________________________________________________________________________
	135	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual,
	136	const TGeoTranslation& translation)
	137	{
	138	// Adds the volume with the specified name and transformation
	139	// to the list of envelopes.
	140	// ---
	141
	142	if (fDebug) {
	143	cout << "... Adding ";
	144	if (!isVirtual) cout << " non-";
	145	cout << "virtual envelope " << name
	146	<< " with translation" << endl;
	147	}
	148
	149	AliMUONGeometryEnvelope* envelope
	150	= new AliMUONGeometryEnvelope(name, isVirtual);
	151	envelope->SetTranslation(translation);
	152
	153	fEnvelopes->Add(envelope);
	154	}
	155
	156	//______________________________________________________________________________
	157	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual,
	158	const TGeoTranslation& translation,
	159	const TGeoRotation& rotation)
	160	{
	161	// Adds the volume with the specified name and transformation
	162	// to the list of envelopes.
	163	// ---
	164
	165	if (fDebug) {
	166	cout << "... Adding ";
	167	if (!isVirtual) cout << " non-";
	168	cout << "virtual envelope " << name
	169	<< " with translation and rotation" << endl;
	170	}
	171
	172	// fEnvelopes->Add(new TGeoCombiTrans(name, translation, rotation));
	173	// would be nice to be so simple
	174
	175	AliMUONGeometryEnvelope* envelope
	176	= new AliMUONGeometryEnvelope(name, isVirtual);
	177	envelope->SetRotation(rotation);
	178	envelope->SetTranslation(translation);
	179
	180	fEnvelopes->Add(envelope);
	181	}
	182
	183	//______________________________________________________________________________
	184	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo)
	185	{
	186	// Adds the volume with the specified name and transformation
	187	// to the list of envelopes.
	188	// ---
	189
	190	if (fDebug) {
	191	cout << "... Adding "
	192	<< " non-virtual envelope " << name
	193	<< " with copyNo " << copyNo << endl;
	194	}
	195
	196	AliMUONGeometryEnvelope* envelope
	197	= new AliMUONGeometryEnvelope(name, copyNo);
	198
	199	fEnvelopes->Add(envelope);
	200	}
	201
	202	//______________________________________________________________________________
	203	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo,
	204	const TGeoTranslation& translation)
	205	{
	206	// Adds the volume with the specified name and transformation
	207	// to the list of envelopes.
	208	// ---
	209
	210	if (fDebug) {
	211	cout << "... Adding "
	212	<< " non-virtual envelope " << name
	213	<< " with copyNo " << copyNo
	214	<< " with translation " << endl;
	215	}
	216
	217	AliMUONGeometryEnvelope* envelope
	218	= new AliMUONGeometryEnvelope(name, copyNo);
	219	envelope->SetTranslation(translation);
	220
	221	fEnvelopes->Add(envelope);
	222	}
	223
	224	//______________________________________________________________________________
	225	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo,
	226	const TGeoTranslation& translation,
	227	const TGeoRotation& rotation)
	228	{
	229	// Adds the volume with the specified name and transformation
	230	// to the list of envelopes.
	231	// ---
	232
	233	if (fDebug) {
	234	cout << "... Adding "
	235	<< " non-virtual envelope " << name
	236	<< " with copyNo " << copyNo
	237	<< " with translation and rotation" << endl;
	238	}
	239
	240	// fEnvelopes->Add(new TGeoCombiTrans(name, translation, rotation));
	241	// would be nice to be so simple
	242
	243	AliMUONGeometryEnvelope* envelope
	244	= new AliMUONGeometryEnvelope(name, copyNo);
	245	envelope->SetRotation(rotation);
	246	envelope->SetTranslation(translation);
	247
	248	fEnvelopes->Add(envelope);
	249	}
	250
	251	//______________________________________________________________________________
	252	void AliMUONChamberGeometry::AddEnvelopeConstituent(const TString& name,
	253	const TString& envName, Int_t copyNo)
	254	{
	255	// Adds the volume with the specified name and transformation
	256	// to the list of envelopes.
	257	// ---
	258
	259	if (fDebug) {
	260	cout << "... Adding constituent " << name
	261	<< " to envelope " << envName
	262	<< " with copyNo " << copyNo << endl;
	263	}
	264
	265	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
	266
	267	if (!envelope) {
	268	// add warning
	269	return;
	270	}
	271
	272	envelope->AddConstituent(name, copyNo);
	273	}
	274
	275	//______________________________________________________________________________
	276	void AliMUONChamberGeometry::AddEnvelopeConstituent(const TString& name,
	277	const TString& envName, Int_t copyNo,
	278	const TGeoTranslation& translation)
	279	{
	280	// Adds the volume with the specified name and transformation
	281	// to the list of envelopes.
	282	// ---
	283
	284	if (fDebug) {
	285	cout << "... Adding constituent " << name
	286	<< " to envelope " << envName
	287	<< " with copyNo " << copyNo
	288	<< " with translation" << endl;
	289	}
	290
	291	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
	292
	293	if (!envelope) {
	294	// add warning
	295	return;
	296	}
	297
	298	envelope->AddConstituent(name, copyNo, translation);
	299	}
	300
	301	//______________________________________________________________________________
	302	void AliMUONChamberGeometry::AddEnvelopeConstituent(const TString& name,
	303	const TString& envName, Int_t copyNo,
	304	const TGeoTranslation& translation,
	305	const TGeoRotation& rotation)
	306	{
	307	// Adds the volume with the specified name and transformation
	308	// to the list of envelopes.
	309	// ---
	310
	311	if (fDebug) {
	312	cout << "... Adding constituent " << name
	313	<< " to envelope " << envName
	314	<< " with copyNo " << copyNo
	315	<< " with translation and rotation" << endl;
	316	}
	317
	318	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
	319
	320	if (!envelope) {
	321	// add warning
	322	return;
	323	}
	324
	325	envelope->AddConstituent(name, copyNo, translation, rotation);
	326	}
	327
	328	//______________________________________________________________________________
	329	void AliMUONChamberGeometry::AddEnvelopeConstituentParam(const TString& name,
	330	const TString& envName, Int_t copyNo,
	331	Int_t npar, Double_t* param)
	332	{
	333	// Adds the volume with the specified name and transformation
	334	// to the list of envelopes.
	335	// ---
	336
	337	if (fDebug) {
	338	cout << "... Adding parameterised constituent " << name
	339	<< " to envelope " << envName
	340	<< " with copyNo " << copyNo << endl;
	341	}
	342
	343	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
	344
	345	if (!envelope) {
	346	// add warning
	347	return;
	348	}
	349
	350	envelope->AddConstituentParam(name, copyNo, npar, param);
	351	}
	352
	353	//______________________________________________________________________________
	354	void AliMUONChamberGeometry::AddEnvelopeConstituentParam(const TString& name,
	355	const TString& envName, Int_t copyNo,
	356	const TGeoTranslation& translation,
	357	Int_t npar, Double_t* param)
	358	{
	359	// Adds the volume with the specified name and transformation
	360	// to the list of envelopes.
	361	// ---
	362
	363	if (fDebug) {
	364	cout << "... Adding parameterised constituent " << name
	365	<< " to envelope " << envName
	366	<< " with copyNo " << copyNo
	367	<< " with translation" << endl;
	368	}
	369
	370	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
	371
	372	if (!envelope) {
	373	// add warning
	374	return;
	375	}
	376
	377	envelope->AddConstituentParam(name, copyNo, translation, npar, param);
	378	}
	379
	380	//______________________________________________________________________________
	381	void AliMUONChamberGeometry::AddEnvelopeConstituentParam(const TString& name,
	382	const TString& envName, Int_t copyNo,
	383	const TGeoTranslation& translation,
	384	const TGeoRotation& rotation,
	385	Int_t npar, Double_t* param)
	386	{
	387	// Adds the volume with the specified name and transformation
	388	// to the list of envelopes.
	389	// ---
	390
	391	if (fDebug) {
	392	cout << "... Adding parameterised constituent " << name
	393	<< " to envelope " << envName
	394	<< " with copyNo " << copyNo
	395	<< " with translation and rotation" << endl;
	396	}
	397
	398	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
	399
	400	if (!envelope) {
	401	// add warning
	402	return;
	403	}
	404
	405	envelope->AddConstituentParam(name, copyNo, translation, rotation, npar, param);
	406	}
	407
	408	//______________________________________________________________________________
	409	void AliMUONChamberGeometry::SetTranslation(const TGeoTranslation& translation)
	410	{
	411	// Sets the chamber position wrt ALIC.
	412	// ---
	413
	414	fTransformation
	415	->SetTranslation(const_cast<Double_t*>(translation.GetTranslation()));
	416	}
	417
	418	//______________________________________________________________________________
	419	void AliMUONChamberGeometry::SetRotation(const TGeoRotation& rotation)
	420	{
	421	// Sets the chamber rotation wrt ALIC.
	422	// ---
	423
	424	TGeoRotation* rot = new TGeoRotation();
	425	rot->SetMatrix(const_cast<Double_t*>(rotation.GetRotationMatrix()));
	426
	427	fTransformation->SetRotation(rot);
	428	}
	429
	430	//______________________________________________________________________________
	431	void AliMUONChamberGeometry::SetSensitiveVolume(Int_t volId)
	432	{
	433	// Adds the volume specified by volId to the list of sensitive
	434	// volumes
	435
	436	fSensVolumeIds->AddAt(volId,fNofSensVolumeIds++);
	437	}
	438
	439	//______________________________________________________________________________
	440	void AliMUONChamberGeometry::SetSensitiveVolume(const TString& volName)
	441	{
	442	// Adds the volume specified by volId to the list of sensitive
	443	// volumes
	444
	445	fSensVolumeIds->AddAt(gMC->VolId(volName),fNofSensVolumeIds++);
	446	}
	447
	448	//______________________________________________________________________________
	449	Bool_t AliMUONChamberGeometry::IsSensitiveVolume(Int_t volId) const
	450	{
	451	// Checks if the volume specified by volId is present in the list
	452	// of sensitive volumes.
	453
	454	for (Int_t i=0; i<fNofSensVolumeIds; i++) {
	455	if (fSensVolumeIds->At(i) == volId) return kTRUE;
	456	}
	457	return kFALSE;
	458	}