[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
// 23/01/2004

#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,
                                          const char* only) 
{
// 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, only);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual,
                                          const TGeoTranslation& translation,
					  const char* only)
{
// 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, only);
  envelope->SetTranslation(translation);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual, 
                                          const TGeoTranslation& translation,
		                          const TGeoRotation& rotation,
					  const char* only)
{
// 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, only);
  envelope->SetRotation(rotation);
  envelope->SetTranslation(translation);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo,
                                          const char* only) 
{
// 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, only);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo,
                                          const TGeoTranslation& translation,
					  const char* only)
{
// 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, only);
  envelope->SetTranslation(translation);

  fEnvelopes->Add(envelope);
}

//______________________________________________________________________________
void  AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo, 
                                          const TGeoTranslation& translation,
		                          const TGeoRotation& rotation,
					  const char* only)
{
// 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, only);
  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
d1cd2474	1	// $Id$
	2	//
	3	// Class AliMUONChamberGeometry
	4	// -----------------------------
5f1df83a	5	// Class for definititon of the MUON chamber positions in ALIC
d1cd2474	6	// Author: Ivana Hrivnacova, IPN Orsay
5f1df83a	7	// 23/01/2004
d1cd2474	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	//______________________________________________________________________________
6b82c1f0	116	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual,
6b82c1f0	117	const char* only)
d1cd2474	118	{
	119	// Adds the volume with the specified name and transformation
	120	// to the list of envelopes.
	121	// ---
	122
	123	if (fDebug) {
	124	cout << "... Adding ";
	125	if (!isVirtual) cout << " non-";
	126	cout << "virtual envelope " << name << endl;
	127	}
	128
	129	AliMUONGeometryEnvelope* envelope
6b82c1f0	130	= new AliMUONGeometryEnvelope(name, isVirtual, only);
d1cd2474	131
	132	fEnvelopes->Add(envelope);
	133	}
	134
	135	//______________________________________________________________________________
	136	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual,
6b82c1f0	137	const TGeoTranslation& translation,
6b82c1f0	138	const char* only)
d1cd2474	139	{
	140	// Adds the volume with the specified name and transformation
	141	// to the list of envelopes.
	142	// ---
	143
	144	if (fDebug) {
	145	cout << "... Adding ";
	146	if (!isVirtual) cout << " non-";
	147	cout << "virtual envelope " << name
	148	<< " with translation" << endl;
	149	}
	150
	151	AliMUONGeometryEnvelope* envelope
6b82c1f0	152	= new AliMUONGeometryEnvelope(name, isVirtual, only);
d1cd2474	153	envelope->SetTranslation(translation);
	154
	155	fEnvelopes->Add(envelope);
	156	}
	157
	158	//______________________________________________________________________________
	159	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Bool_t isVirtual,
	160	const TGeoTranslation& translation,
6b82c1f0	161	const TGeoRotation& rotation,
6b82c1f0	162	const char* only)
d1cd2474	163	{
	164	// Adds the volume with the specified name and transformation
	165	// to the list of envelopes.
	166	// ---
	167
	168	if (fDebug) {
	169	cout << "... Adding ";
	170	if (!isVirtual) cout << " non-";
	171	cout << "virtual envelope " << name
	172	<< " with translation and rotation" << endl;
	173	}
	174
	175	// fEnvelopes->Add(new TGeoCombiTrans(name, translation, rotation));
	176	// would be nice to be so simple
	177
	178	AliMUONGeometryEnvelope* envelope
6b82c1f0	179	= new AliMUONGeometryEnvelope(name, isVirtual, only);
d1cd2474	180	envelope->SetRotation(rotation);
	181	envelope->SetTranslation(translation);
	182
	183	fEnvelopes->Add(envelope);
	184	}
	185
	186	//______________________________________________________________________________
6b82c1f0	187	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo,
6b82c1f0	188	const char* only)
d1cd2474	189	{
	190	// Adds the volume with the specified name and transformation
	191	// to the list of envelopes.
	192	// ---
	193
	194	if (fDebug) {
	195	cout << "... Adding "
	196	<< " non-virtual envelope " << name
	197	<< " with copyNo " << copyNo << endl;
	198	}
	199
	200	AliMUONGeometryEnvelope* envelope
6b82c1f0	201	= new AliMUONGeometryEnvelope(name, copyNo, only);
d1cd2474	202
	203	fEnvelopes->Add(envelope);
	204	}
	205
	206	//______________________________________________________________________________
	207	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo,
6b82c1f0	208	const TGeoTranslation& translation,
6b82c1f0	209	const char* only)
d1cd2474	210	{
	211	// Adds the volume with the specified name and transformation
	212	// to the list of envelopes.
	213	// ---
	214
	215	if (fDebug) {
	216	cout << "... Adding "
	217	<< " non-virtual envelope " << name
	218	<< " with copyNo " << copyNo
	219	<< " with translation " << endl;
	220	}
	221
	222	AliMUONGeometryEnvelope* envelope
6b82c1f0	223	= new AliMUONGeometryEnvelope(name, copyNo, only);
d1cd2474	224	envelope->SetTranslation(translation);
	225
	226	fEnvelopes->Add(envelope);
	227	}
	228
	229	//______________________________________________________________________________
	230	void AliMUONChamberGeometry::AddEnvelope(const TString& name, Int_t copyNo,
	231	const TGeoTranslation& translation,
6b82c1f0	232	const TGeoRotation& rotation,
6b82c1f0	233	const char* only)
d1cd2474	234	{
	235	// Adds the volume with the specified name and transformation
	236	// to the list of envelopes.
	237	// ---
	238
	239	if (fDebug) {
	240	cout << "... Adding "
	241	<< " non-virtual envelope " << name
	242	<< " with copyNo " << copyNo
	243	<< " with translation and rotation" << endl;
	244	}
	245
	246	// fEnvelopes->Add(new TGeoCombiTrans(name, translation, rotation));
	247	// would be nice to be so simple
	248
	249	AliMUONGeometryEnvelope* envelope
6b82c1f0	250	= new AliMUONGeometryEnvelope(name, copyNo, only);
d1cd2474	251	envelope->SetRotation(rotation);
	252	envelope->SetTranslation(translation);
	253
	254	fEnvelopes->Add(envelope);
	255	}
	256
	257	//______________________________________________________________________________
	258	void AliMUONChamberGeometry::AddEnvelopeConstituent(const TString& name,
	259	const TString& envName, Int_t copyNo)
	260	{
	261	// Adds the volume with the specified name and transformation
	262	// to the list of envelopes.
	263	// ---
	264
	265	if (fDebug) {
	266	cout << "... Adding constituent " << name
	267	<< " to envelope " << envName
	268	<< " with copyNo " << copyNo << endl;
	269	}
	270
	271	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
	272
	273	if (!envelope) {
	274	// add warning
	275	return;
	276	}
	277
	278	envelope->AddConstituent(name, copyNo);
	279	}
	280
	281	//______________________________________________________________________________
	282	void AliMUONChamberGeometry::AddEnvelopeConstituent(const TString& name,
	283	const TString& envName, Int_t copyNo,
	284	const TGeoTranslation& translation)
	285	{
	286	// Adds the volume with the specified name and transformation
	287	// to the list of envelopes.
	288	// ---
	289
	290	if (fDebug) {
	291	cout << "... Adding constituent " << name
	292	<< " to envelope " << envName
	293	<< " with copyNo " << copyNo
	294	<< " with translation" << endl;
	295	}
	296
	297	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
	298
	299	if (!envelope) {
	300	// add warning
	301	return;
	302	}
	303
	304	envelope->AddConstituent(name, copyNo, translation);
	305	}
	306
	307	//______________________________________________________________________________
	308	void AliMUONChamberGeometry::AddEnvelopeConstituent(const TString& name,
	309	const TString& envName, Int_t copyNo,
	310	const TGeoTranslation& translation,
	311	const TGeoRotation& rotation)
	312	{
	313	// Adds the volume with the specified name and transformation
	314	// to the list of envelopes.
315	// ---
316
317	if (fDebug) {
318	cout << "... Adding constituent " << name
319	<< " to envelope " << envName
320	<< " with copyNo " << copyNo
321	<< " with translation and rotation" << endl;
322	}
323
324	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
325
326	if (!envelope) {
327	// add warning
328	return;
329	}
330
331	envelope->AddConstituent(name, copyNo, translation, rotation);
332	}
333
334	//______________________________________________________________________________
335	void AliMUONChamberGeometry::AddEnvelopeConstituentParam(const TString& name,
336	const TString& envName, Int_t copyNo,
337	Int_t npar, Double_t* param)
338	{
339	// Adds the volume with the specified name and transformation
340	// to the list of envelopes.
341	// ---
342
343	if (fDebug) {
344	cout << "... Adding parameterised constituent " << name
345	<< " to envelope " << envName
346	<< " with copyNo " << copyNo << endl;
347	}
348
349	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
350
351	if (!envelope) {
352	// add warning
353	return;
354	}
355
356	envelope->AddConstituentParam(name, copyNo, npar, param);
357	}
358
359	//______________________________________________________________________________
360	void AliMUONChamberGeometry::AddEnvelopeConstituentParam(const TString& name,
361	const TString& envName, Int_t copyNo,
362	const TGeoTranslation& translation,
363	Int_t npar, Double_t* param)
364	{
365	// Adds the volume with the specified name and transformation
366	// to the list of envelopes.
367	// ---
368
369	if (fDebug) {
370	cout << "... Adding parameterised constituent " << name
371	<< " to envelope " << envName
372	<< " with copyNo " << copyNo
373	<< " with translation" << endl;
374	}
375
376	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
377
378	if (!envelope) {
379	// add warning
380	return;
381	}
382
383	envelope->AddConstituentParam(name, copyNo, translation, npar, param);
384	}
385
386	//______________________________________________________________________________
387	void AliMUONChamberGeometry::AddEnvelopeConstituentParam(const TString& name,
388	const TString& envName, Int_t copyNo,
389	const TGeoTranslation& translation,
390	const TGeoRotation& rotation,
391	Int_t npar, Double_t* param)
392	{
393	// Adds the volume with the specified name and transformation
394	// to the list of envelopes.
395	// ---
396
397	if (fDebug) {
398	cout << "... Adding parameterised constituent " << name
399	<< " to envelope " << envName
400	<< " with copyNo " << copyNo
401	<< " with translation and rotation" << endl;
402	}
403
404	AliMUONGeometryEnvelope* envelope = FindEnvelope(envName);
405
406	if (!envelope) {
407	// add warning
408	return;
409	}
410
411	envelope->AddConstituentParam(name, copyNo, translation, rotation, npar, param);
412	}
413
414	//______________________________________________________________________________
415	void AliMUONChamberGeometry::SetTranslation(const TGeoTranslation& translation)
416	{
417	// Sets the chamber position wrt ALIC.
418	// ---
419
420	fTransformation
421	->SetTranslation(const_cast<Double_t*>(translation.GetTranslation()));
422	}
423
424	//______________________________________________________________________________
425	void AliMUONChamberGeometry::SetRotation(const TGeoRotation& rotation)
426	{
427	// Sets the chamber rotation wrt ALIC.
428	// ---
429
430	TGeoRotation* rot = new TGeoRotation();
431	rot->SetMatrix(const_cast<Double_t*>(rotation.GetRotationMatrix()));
432
433	fTransformation->SetRotation(rot);
434	}
435
436	//______________________________________________________________________________
437	void AliMUONChamberGeometry::SetSensitiveVolume(Int_t volId)
438	{
439	// Adds the volume specified by volId to the list of sensitive
440	// volumes
441
442	fSensVolumeIds->AddAt(volId,fNofSensVolumeIds++);
443	}
444
445	//______________________________________________________________________________
446	void AliMUONChamberGeometry::SetSensitiveVolume(const TString& volName)
447	{
448	// Adds the volume specified by volId to the list of sensitive
449	// volumes
450
451	fSensVolumeIds->AddAt(gMC->VolId(volName),fNofSensVolumeIds++);
452	}
453
454	//______________________________________________________________________________
455	Bool_t AliMUONChamberGeometry::IsSensitiveVolume(Int_t volId) const
456	{
457	// Checks if the volume specified by volId is present in the list
458	// of sensitive volumes.
459
460	for (Int_t i=0; i<fNofSensVolumeIds; i++) {
461	if (fSensVolumeIds->At(i) == volId) return kTRUE;
462	}
463	return kFALSE;
464	}