[u/mrichter/AliRoot.git] / MUON / AliMUONSt2GeometryBuilderV2.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$
//
// Class AliMUONSt2GeometryBuilderV2
// -------------------------------
// MUON Station2 coarse geometry construction class.
// Extracted from AliMUONv1
// Dummy version of station 2 with the right DE id (Ch. Finck)


#include <TVirtualMC.h>
#include <TGeoMatrix.h>

#include "AliLog.h"

#include "AliMUONSt2GeometryBuilderV2.h"
#include "AliMUON.h"
#include "AliMUONChamber.h"
#include "AliMUONGeometryModule.h"
#include "AliMUONGeometryEnvelopeStore.h"

ClassImp(AliMUONSt2GeometryBuilderV2)

//______________________________________________________________________________
AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2(AliMUON* muon)
 : AliMUONVGeometryBuilder("st2V2.dat",
                           muon->Chamber(2).GetGeometry(), 
			   muon->Chamber(3).GetGeometry()),
   fMUON(muon)
{
// Standard constructor

}

//______________________________________________________________________________
AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2()
 : AliMUONVGeometryBuilder(),
   fMUON(0)
{
// Default constructor
}


//______________________________________________________________________________
AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2(const AliMUONSt2GeometryBuilderV2& rhs)
  : AliMUONVGeometryBuilder(rhs)
{
// Protected copy constructor

  AliFatal("Copy constructor is not implemented.");
}

//______________________________________________________________________________
AliMUONSt2GeometryBuilderV2::~AliMUONSt2GeometryBuilderV2() {
//
}

//______________________________________________________________________________
AliMUONSt2GeometryBuilderV2& 
AliMUONSt2GeometryBuilderV2::operator = (const AliMUONSt2GeometryBuilderV2& rhs) 
{
// Protected assignement operator

  // check assignement to self
  if (this == &rhs) return *this;

  AliFatal("Assignment operator is not implemented.");
    
  return *this;  
}

//
// public methods
//

//______________________________________________________________________________
void AliMUONSt2GeometryBuilderV2::CreateGeometry() 
{
// From AliMUONv1::CreateGeometry()

//
//********************************************************************
//                            Station 2                             **
//********************************************************************
     // indices 1 and 2 for first and second chambers in the station
     // iChamber (first chamber) kept for other quanties than Z,
     // assumed to be the same in both chambers

     // Get tracking medias Ids     
     Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
     Int_t idAir= idtmed[1100]; // medium 1
  //     Int_t idAlu1=idtmed[1103]; // medium 4
     //     Int_t idAlu2=idtmed[1104]; // medium 5
     Int_t idGas=idtmed[1108];  // medium 9 = Ar-CO2 gas (80%+20%)

     const Float_t kDeltaQuad = 0.01;//2.6; dummy value til we find the good value
     const Float_t kDeltaZ = 6.5/2.;

     // Rotation matrices in the x-y plane  
     // phi= 0 deg
     Int_t irot1;
     fMUON->AliMatrix(irot1,  90.,   0., 90.,  90., 0., 0.);
     // phi= 90 deg
     Int_t irot2;
     fMUON->AliMatrix(irot2,  90.,  90., 90., 180., 0., 0.);

     AliMUONChamber* iChamber = &fMUON->Chamber(2);
     //     AliMUONChamber* iChamber1 = iChamber;
     //     AliMUONChamber* iChamber2 = &fMUON->Chamber(3);
     
     // Half of the total thickness of frame crosses (including DAlu)
     // for each chamber in stations 1 and 2:
     // 3% of X0 of composite material,
     // but taken as Aluminium here, with same thickness in number of X0
     //     Float_t dframez = 3. * 8.9 / 100;
     // DGas and DAlu not changed from standard values
     //     Double_t zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
             // The same parameters are defined in builder for station 1 
     
     // Mother volume
     // Outer excess and inner recess for mother volume radius
     // with respect to ROuter and RInner
     //     Float_t dframep=.001; // Value for station 3 should be 6 ...
     // Width (RdPhi) of the frame crosses for stations 1 and 2 (cm)
     // Float_t dframep1=.001;
     //     Float_t phi=2*TMath::Pi()/12/2;
             // The same parameters are defined in builder for station 1 (30deg)

     Float_t tpar[5];
     //     Double_t dstation = (-iChamber2->Z()) - (-iChamber1->Z());


     Float_t posx, posy, posz;
//   Chamber Material represented by Alu sheet
     tpar[0]= iChamber->RInner();
     tpar[1]= iChamber->ROuter();
     tpar[2] = (iChamber->DGas())/2;
     tpar[3] = 0.;
     tpar[4] = 90.;
  
//     
//   Sensitive volumes
     gMC->Gsvolu("S03G", "TUBS", idGas, tpar, 5);
     gMC->Gsvolu("S04G", "TUBS", idGas, tpar, 5);

     Int_t detElemId;

     posx = kDeltaQuad;
     posy = kDeltaQuad;
     posz = -kDeltaZ;

     detElemId = 351;
     gMC->Gsvolu("LE01", "TUBS", idAir, tpar, 5);
     GetEnvelopes(2)->AddEnvelope("LE01", detElemId, true, TGeoTranslation(posx, posy, posz),
				   TGeoRotation("rot1",90,0,90,90,0,0) );
     detElemId = 451;
     gMC->Gsvolu("LF01", "TUBS", idAir, tpar, 5);
     GetEnvelopes(3)->AddEnvelope("LF01", detElemId, true, TGeoTranslation(posx, posy, posz),
				   TGeoRotation("rot1",90,0,90,90,0,0) );
     detElemId = 300;
     gMC->Gsvolu("LE02", "TUBS", idAir, tpar, 5);
     GetEnvelopes(2)->AddEnvelope("LE02", detElemId, true, TGeoTranslation(-posx, posy,-posz),
				  TGeoRotation("rot2",90,180,90,90,180,0) );
     detElemId = 400;
     gMC->Gsvolu("LF02", "TUBS", idAir, tpar, 5);
     GetEnvelopes(3)->AddEnvelope("LF02", detElemId, true, TGeoTranslation(-posx, posy,-posz),
				  TGeoRotation("rot2",90,180,90,90,180,0) );
     detElemId = 350;
     gMC->Gsvolu("LE03", "TUBS", idAir, tpar, 5);
     GetEnvelopes(2)->AddEnvelope("LE03", detElemId, true, TGeoTranslation(posx, -posy, -posz),
				  TGeoRotation("rot3",90,0,90,270,180,0) );
     detElemId = 450;
     gMC->Gsvolu("LF03", "TUBS", idAir, tpar, 5);
     GetEnvelopes(3)->AddEnvelope("LF03", detElemId, true, TGeoTranslation(posx, -posy, -posz),
				  TGeoRotation("rot3",90,0,90,270,180,0) );
     detElemId = 301;
     gMC->Gsvolu("LE04", "TUBS", idAir, tpar, 5);
     GetEnvelopes(2)->AddEnvelope("LE04", detElemId, true, TGeoTranslation(-posx, -posy, posz),
				  TGeoRotation("rot4",90,180,90,270,0,0) );
     detElemId = 401;
     gMC->Gsvolu("LF04", "TUBS", idAir, tpar, 5);
     GetEnvelopes(3)->AddEnvelope("LF04", detElemId, true, TGeoTranslation(-posx, -posy, posz),
				  TGeoRotation("rot4",90,180,90,270,0,0) );

     GetEnvelopes(2)->AddEnvelopeConstituent("S03G", "LE01", 1);
     GetEnvelopes(3)->AddEnvelopeConstituent("S04G", "LF01", 1);

     GetEnvelopes(2)->AddEnvelopeConstituent("S03G", "LE02", 2);
     GetEnvelopes(3)->AddEnvelopeConstituent("S04G", "LF02", 2);

     GetEnvelopes(2)->AddEnvelopeConstituent("S03G", "LE03", 3);
     GetEnvelopes(3)->AddEnvelopeConstituent("S04G", "LF03", 3);

     GetEnvelopes(2)->AddEnvelopeConstituent("S03G", "LE04", 4);
     GetEnvelopes(3)->AddEnvelopeConstituent("S04G", "LF04", 4);

}

//______________________________________________________________________________
void AliMUONSt2GeometryBuilderV2::SetTransformations()
{
// Defines the transformations for the station2 chambers.
// ---

  AliMUONChamber* iChamber1 = &fMUON->Chamber(2);
  Double_t zpos1 = - iChamber1->Z(); 
  iChamber1->GetGeometry()
    ->SetTranslation(TGeoTranslation(0., 0., zpos1));

  AliMUONChamber* iChamber2 = &fMUON->Chamber(3);
  Double_t zpos2 = - iChamber2->Z(); 
  iChamber2->GetGeometry()
    ->SetTranslation(TGeoTranslation(0., 0., zpos2));
}

//______________________________________________________________________________
void AliMUONSt2GeometryBuilderV2::SetSensitiveVolumes()
{
// Defines the sensitive volumes for station2 chambers.
// ---

  GetGeometry(2)->SetSensitiveVolume("S03G");
  GetGeometry(3)->SetSensitiveVolume("S04G");
}
Commit	Line	Data
a2b41b89	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	// Class AliMUONSt2GeometryBuilderV2
	19	// -------------------------------
	20	// MUON Station2 coarse geometry construction class.
	21	// Extracted from AliMUONv1
e6b1b13c	22	// Dummy version of station 2 with the right DE id (Ch. Finck)
a2b41b89	23
	24
	25	#include <TVirtualMC.h>
	26	#include <TGeoMatrix.h>
	27
	28	#include "AliLog.h"
	29
	30	#include "AliMUONSt2GeometryBuilderV2.h"
	31	#include "AliMUON.h"
	32	#include "AliMUONChamber.h"
	33	#include "AliMUONGeometryModule.h"
	34	#include "AliMUONGeometryEnvelopeStore.h"
	35
	36	ClassImp(AliMUONSt2GeometryBuilderV2)
	37
	38	//______________________________________________________________________________
	39	AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2(AliMUON* muon)
ff707ad4	40	: AliMUONVGeometryBuilder("st2V2.dat",
a2b41b89	41	muon->Chamber(2).GetGeometry(),
	42	muon->Chamber(3).GetGeometry()),
	43	fMUON(muon)
	44	{
	45	// Standard constructor
	46
	47	}
	48
	49	//______________________________________________________________________________
	50	AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2()
	51	: AliMUONVGeometryBuilder(),
	52	fMUON(0)
	53	{
	54	// Default constructor
	55	}
	56
	57
	58	//______________________________________________________________________________
	59	AliMUONSt2GeometryBuilderV2::AliMUONSt2GeometryBuilderV2(const AliMUONSt2GeometryBuilderV2& rhs)
	60	: AliMUONVGeometryBuilder(rhs)
	61	{
	62	// Protected copy constructor
	63
	64	AliFatal("Copy constructor is not implemented.");
	65	}
	66
	67	//______________________________________________________________________________
	68	AliMUONSt2GeometryBuilderV2::~AliMUONSt2GeometryBuilderV2() {
	69	//
	70	}
	71
	72	//______________________________________________________________________________
	73	AliMUONSt2GeometryBuilderV2&
	74	AliMUONSt2GeometryBuilderV2::operator = (const AliMUONSt2GeometryBuilderV2& rhs)
	75	{
	76	// Protected assignement operator
	77
	78	// check assignement to self
	79	if (this == &rhs) return *this;
	80
	81	AliFatal("Assignment operator is not implemented.");
	82
	83	return *this;
	84	}
	85
	86	//
	87	// public methods
	88	//
	89
	90	//______________________________________________________________________________
	91	void AliMUONSt2GeometryBuilderV2::CreateGeometry()
	92	{
	93	// From AliMUONv1::CreateGeometry()
	94
	95	//
	96	//********************************************************************
	97	// Station 2 **
	98	//********************************************************************
	99	// indices 1 and 2 for first and second chambers in the station
	100	// iChamber (first chamber) kept for other quanties than Z,
	101	// assumed to be the same in both chambers
	102
	103	// Get tracking medias Ids
	104	Int_t *idtmed = fMUON->GetIdtmed()->GetArray()-1099;
105	Int_t idAir= idtmed[1100]; // medium 1
106	// Int_t idAlu1=idtmed[1103]; // medium 4
107	// Int_t idAlu2=idtmed[1104]; // medium 5
108	Int_t idGas=idtmed[1108]; // medium 9 = Ar-CO2 gas (80%+20%)
109
e6b1b13c	110	const Float_t kDeltaQuad = 0.01;//2.6; dummy value til we find the good value
a2b41b89	111	const Float_t kDeltaZ = 6.5/2.;
	112
	113	// Rotation matrices in the x-y plane
	114	// phi= 0 deg
	115	Int_t irot1;
	116	fMUON->AliMatrix(irot1, 90., 0., 90., 90., 0., 0.);
	117	// phi= 90 deg
	118	Int_t irot2;
	119	fMUON->AliMatrix(irot2, 90., 90., 90., 180., 0., 0.);
	120
	121	AliMUONChamber* iChamber = &fMUON->Chamber(2);
	122	// AliMUONChamber* iChamber1 = iChamber;
	123	// AliMUONChamber* iChamber2 = &fMUON->Chamber(3);
	124
	125	// Half of the total thickness of frame crosses (including DAlu)
	126	// for each chamber in stations 1 and 2:
	127	// 3% of X0 of composite material,
	128	// but taken as Aluminium here, with same thickness in number of X0
	129	// Float_t dframez = 3. * 8.9 / 100;
	130	// DGas and DAlu not changed from standard values
	131	// Double_t zfpos=-(iChamber->DGas()+dframez+iChamber->DAlu())/2;
	132	// The same parameters are defined in builder for station 1
	133
	134	// Mother volume
	135	// Outer excess and inner recess for mother volume radius
	136	// with respect to ROuter and RInner
	137	// Float_t dframep=.001; // Value for station 3 should be 6 ...
	138	// Width (RdPhi) of the frame crosses for stations 1 and 2 (cm)
	139	// Float_t dframep1=.001;
	140	// Float_t phi=2*TMath::Pi()/12/2;
	141	// The same parameters are defined in builder for station 1 (30deg)
	142
	143	Float_t tpar[5];
	144	// Double_t dstation = (-iChamber2->Z()) - (-iChamber1->Z());
	145
	146
	147	Float_t posx, posy, posz;
	148	// Chamber Material represented by Alu sheet
	149	tpar[0]= iChamber->RInner();
	150	tpar[1]= iChamber->ROuter();
	151	tpar[2] = (iChamber->DGas())/2;
	152	tpar[3] = 0.;
	153	tpar[4] = 90.;
	154
	155	//
	156	// Sensitive volumes
	157	gMC->Gsvolu("S03G", "TUBS", idGas, tpar, 5);
	158	gMC->Gsvolu("S04G", "TUBS", idGas, tpar, 5);
	159
	160	Int_t detElemId;
	161
	162	posx = kDeltaQuad;
	163	posy = kDeltaQuad;
e6b1b13c	164	posz = -kDeltaZ;
a2b41b89	165
e7addd77	166	detElemId = 351;
a2b41b89	167	gMC->Gsvolu("LE01", "TUBS", idAir, tpar, 5);
	168	GetEnvelopes(2)->AddEnvelope("LE01", detElemId, true, TGeoTranslation(posx, posy, posz),
	169	TGeoRotation("rot1",90,0,90,90,0,0) );
e7addd77	170	detElemId = 451;
a2b41b89	171	gMC->Gsvolu("LF01", "TUBS", idAir, tpar, 5);
	172	GetEnvelopes(3)->AddEnvelope("LF01", detElemId, true, TGeoTranslation(posx, posy, posz),
	173	TGeoRotation("rot1",90,0,90,90,0,0) );
e7addd77	174	detElemId = 300;
a2b41b89	175	gMC->Gsvolu("LE02", "TUBS", idAir, tpar, 5);
	176	GetEnvelopes(2)->AddEnvelope("LE02", detElemId, true, TGeoTranslation(-posx, posy,-posz),
	177	TGeoRotation("rot2",90,180,90,90,180,0) );
e7addd77	178	detElemId = 400;
a2b41b89	179	gMC->Gsvolu("LF02", "TUBS", idAir, tpar, 5);
	180	GetEnvelopes(3)->AddEnvelope("LF02", detElemId, true, TGeoTranslation(-posx, posy,-posz),
	181	TGeoRotation("rot2",90,180,90,90,180,0) );
e7addd77	182	detElemId = 350;
a2b41b89	183	gMC->Gsvolu("LE03", "TUBS", idAir, tpar, 5);
	184	GetEnvelopes(2)->AddEnvelope("LE03", detElemId, true, TGeoTranslation(posx, -posy, -posz),
	185	TGeoRotation("rot3",90,0,90,270,180,0) );
e7addd77	186	detElemId = 450;
a2b41b89	187	gMC->Gsvolu("LF03", "TUBS", idAir, tpar, 5);
	188	GetEnvelopes(3)->AddEnvelope("LF03", detElemId, true, TGeoTranslation(posx, -posy, -posz),
	189	TGeoRotation("rot3",90,0,90,270,180,0) );
e7addd77	190	detElemId = 301;
a2b41b89	191	gMC->Gsvolu("LE04", "TUBS", idAir, tpar, 5);
	192	GetEnvelopes(2)->AddEnvelope("LE04", detElemId, true, TGeoTranslation(-posx, -posy, posz),
	193	TGeoRotation("rot4",90,180,90,270,0,0) );
e7addd77	194	detElemId = 401;
a2b41b89	195	gMC->Gsvolu("LF04", "TUBS", idAir, tpar, 5);
	196	GetEnvelopes(3)->AddEnvelope("LF04", detElemId, true, TGeoTranslation(-posx, -posy, posz),
	197	TGeoRotation("rot4",90,180,90,270,0,0) );
	198
	199	GetEnvelopes(2)->AddEnvelopeConstituent("S03G", "LE01", 1);
	200	GetEnvelopes(3)->AddEnvelopeConstituent("S04G", "LF01", 1);
	201
	202	GetEnvelopes(2)->AddEnvelopeConstituent("S03G", "LE02", 2);
	203	GetEnvelopes(3)->AddEnvelopeConstituent("S04G", "LF02", 2);
	204
	205	GetEnvelopes(2)->AddEnvelopeConstituent("S03G", "LE03", 3);
	206	GetEnvelopes(3)->AddEnvelopeConstituent("S04G", "LF03", 3);
	207
	208	GetEnvelopes(2)->AddEnvelopeConstituent("S03G", "LE04", 4);
	209	GetEnvelopes(3)->AddEnvelopeConstituent("S04G", "LF04", 4);
	210
	211	}
	212
	213	//______________________________________________________________________________
	214	void AliMUONSt2GeometryBuilderV2::SetTransformations()
	215	{
	216	// Defines the transformations for the station2 chambers.
	217	// ---
	218
	219	AliMUONChamber* iChamber1 = &fMUON->Chamber(2);
	220	Double_t zpos1 = - iChamber1->Z();
	221	iChamber1->GetGeometry()
	222	->SetTranslation(TGeoTranslation(0., 0., zpos1));
	223
	224	AliMUONChamber* iChamber2 = &fMUON->Chamber(3);
	225	Double_t zpos2 = - iChamber2->Z();
	226	iChamber2->GetGeometry()
	227	->SetTranslation(TGeoTranslation(0., 0., zpos2));
	228	}
	229
	230	//______________________________________________________________________________
	231	void AliMUONSt2GeometryBuilderV2::SetSensitiveVolumes()
	232	{
	233	// Defines the sensitive volumes for station2 chambers.
	234	// ---
	235
	236	GetGeometry(2)->SetSensitiveVolume("S03G");
	237	GetGeometry(3)->SetSensitiveVolume("S04G");
	238	}