[u/mrichter/AliRoot.git] / STRUCT / AliMAG.cxx

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  L3 Magnet                                                                //
//                                                                           //
//Begin_Html
/*
<img src="gif/AliMAGClass.gif">
</pre>
<br clear=left>
<font size=+2 color=red>
<p>The responsible person for this module is
<a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
</font>
<pre>
*/
//End_Html
//                                                                           //
///////////////////////////////////////////////////////////////////////////////
 
#include "AliMAG.h"
#include "AliRun.h"
#include "TGeant3.h"
 
ClassImp(AliMAG)
 
//_____________________________________________________________________________
AliMAG::AliMAG()
{
  //
  // Default constructor for L3 magnet
  //
}
 
//_____________________________________________________________________________
AliMAG::AliMAG(const char *name, const char *title)
  : AliDetector(name,title)
{
  //
  // Standard constructor for L3 magnet
  //
  //Begin_Html
  /*
    <img src="gif/aliMAG.gif">
  */
  //End_Html
  
  SetMarkerColor(7);
  SetMarkerStyle(2);
  SetMarkerSize(0.4);
}

//_____________________________________________________________________________
void AliMAG::BuildGeometry()
{
  //
  // Build ROOT TNode geometry for L3 magnet
  // only for detectors
  //
}
 
//_____________________________________________________________________________
void AliMAG::CreateGeometry()
{
  //
  // Create geometry for L3 magnet
  //
  //Begin_Html
  /*
    <img src="gif/mag.gif">
  */
  //End_Html
  //Begin_Html
  /*
    <img src="gif/tree_mag.gif">
  */
  //End_Html
  
  AliMC* pMC = AliMC::GetMC();
  
  Int_t *idtmed = gAlice->Idtmed();
  
  Float_t dpar[13];
  Int_t idrotm[399];
  Float_t par[10];
  
  // ANGLE POLAIRE MAXIMUM 
  
  //  Define Mother 
  
  par[0] = 22.5;
  par[1] = 360.;
  par[2] = 8.;
  par[3] = 2.;
  par[4] = -600.;
  par[5] = 580.;
  par[6] = 790.;
  par[7] = 600.;
  par[8] = 580.;
  par[9] = 790.;
  pMC->Gsvolu("L3MO", "PGON", idtmed[334], par, 10);
  pMC->Gspos("L3MO", 1, "ALIC", 0., -30., 0., 0, "ONLY");
  
  // Define coils 
  
  par[5] = 585.;
  par[6] = 690.;
  par[8] = 585.;
  par[9] = 690.;
  pMC->Gsvolu("L3CO", "PGON", idtmed[328], par, 10);
  pMC->Gspos("L3CO", 1, "L3MO", 0., 0., 0., 0, "ONLY");
  
  par[5] = 580.;
  par[6] = 585.;
  par[8] = 580.;
  par[9] = 585.;
  pMC->Gsvolu("L3C1", "PGON", idtmed[308], par, 10);
  pMC->Gspos("L3C1", 1, "L3MO", 0., 0., 0., 0, "ONLY");
  
  // Define yoke 
  
  par[5] = 690.;
  par[6] = 790.;
  par[8] = 690.;
  par[9] = 790.;
  pMC->Gsvolu("L3YO", "PGON", idtmed[329], par, 10);
  pMC->Gspos("L3YO", 1, "L3MO", 0., 0., 0., 0, "ONLY");
  
  // Define the return yoke of L3 (DOOR) 
  
  par[4] = 600.;
  par[5] = 0.;
  par[6] = 790.;
  par[7] = 700.;
  par[8] = par[5];
  par[9] = par[6];
  pMC->Gsvolu("L3DO", "PGON", idtmed[334], par, 10);
  
  par[4] = 610.;
  par[5] = 0.;
  par[6] = 790.;
  par[7] = 700.;
  par[8] = par[5];
  par[9] = par[6];
  pMC->Gsvolu("L3FR", "PGON", idtmed[329], par, 10);
  
  // INNER LAYER 
  
  par[4] = 600.;
  par[7] = 610.;
  pMC->Gsvolu("L3IR", "PGON", idtmed[309], par, 10);
  
  //     DOOR OPENING 
  
  dpar[0] = 22.5;
  dpar[1] = 360.;
  dpar[2] = 8.;
  dpar[3] = 3.;
  dpar[4] = 610.;
  dpar[5] = 0.;
  dpar[6] = 163.5;
  dpar[7] = 670.;
  dpar[8] = dpar[5];
  dpar[9] = dpar[6];
  dpar[10] = 700.;
  dpar[11] = dpar[5];
  dpar[12] = dpar[6] + 50.;
  pMC->Gsvolu("L3O1", "PGON", idtmed[314], dpar, 13);
  par[4] = 600.;
  par[5] = 0.;
  par[6] = 163.5;
  par[7] = 610.;
  par[8] = 0.;
  par[9] = 163.5;
  pMC->Gsvolu("L3O2", "PGON", idtmed[314], par, 10);
  
  //     THE DOOR OPENING HAS TO BE PLACED WITH 'MANY' SINCE THE REGION 
  //     WILL CONTAIN A MUON CHAMBER, BEAM PIPE AND BEAM SHIELD 
  //     PLACED WITH 'ONLY'. 
  
  pMC->Gspos("L3O1", 1, "L3FR", 0., 30., 0., 0, "MANY");
  pMC->Gspos("L3O2", 1, "L3IR", 0., 30., 0., 0, "MANY");
  
  pMC->Gspos("L3FR", 1, "L3DO", 0., 0., 0., 0, "MANY");
  pMC->Gspos("L3IR", 1, "L3DO", 0., 0., 0., 0, "MANY");
  
  pMC->Gspos("L3DO", 1, "ALIC", 0., -30., 0., 0, "MANY");
  AliMatrix(idrotm[300], 90., 0., 90., 90., 180., 0.);
  pMC->Gspos("L3DO", 2, "ALIC", 0., -30., 0., idrotm[300], "MANY");
}

//_____________________________________________________________________________
void AliMAG::CreateMaterials()
{
  //
  // Create materials for L3 magnet
  //
  
  Int_t   ISXFLD = gAlice->Field()->Integ();
  Float_t SXMGMX = gAlice->Field()->Max();
  
  Float_t epsil, stmin, deemax, tmaxfd, stemax;
  // --- Define the various materials for GEANT --- 
  
  //     Aluminum 
  AliMaterial(9, "Al$", 26.98, 13., 2.7, 8.9, 37.2);
  AliMaterial(29, "Al$", 26.98, 13., 2.7, 8.9, 37.2);
  
  //     Iron 
  AliMaterial(10, "Fe$", 55.85, 26., 7.87, 1.76, 17.1);
  AliMaterial(30, "Fe$", 55.85, 26., 7.87, 1.76, 17.1);
  
  //     Air 
  AliMaterial(15, "Air$", 14.61, 7.3, .001205, 30420., 67500);
  AliMaterial(35, "Air$", 14.61, 7.3, .001205, 30420., 67500);
  
  // **************** 
  //     Defines tracking media parameters. 
  //     Les valeurs sont commentees pour laisser le defaut 
  //     a GEANT (version 3-21, page CONS200), f.m. 
  epsil  = .001;  // Tracking precision, 
  stemax = -1.;   // Maximum displacement for multiple scat 
  tmaxfd = -20.;  // Maximum angle due to field deflection 
  deemax = -.3;   // Maximum fractional energy loss, DLS 
  stmin  = -.8;
  // *************** 
  
  //    IRON 
  
  AliMedium(310, "FE_C0             ", 10, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
  AliMedium(330, "FE_C1             ", 30, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
  
  //     ALUMINUM 

  AliMedium(309, "ALU_C0            ",  9, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
  AliMedium(329, "ALU_C1            ", 29, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
  
  //     AIR 
  
  AliMedium(315, "AIR_C0            ", 15, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
  AliMedium(335, "AIR_C1            ", 35, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
}

//_____________________________________________________________________________
void AliMAG::DrawDetector()
{
  //
  // Draw a shaded view of the L3 magnet
  //
}

//_____________________________________________________________________________
void AliMAG::Init()
{
  //
  // Initialise L3 magnet after it has been built
  Int_t i;
  //
  printf("\n");
  for(i=0;i<35;i++) printf("*");
  printf(" MAG_INIT ");
  for(i=0;i<35;i++) printf("*");
  printf("\n");
  //
  // Here the MAG initialisation code (if any!)
  for(i=0;i<80;i++) printf("*");
  printf("\n");
}

//_____________________________________________________________________________
void AliMAG::StepManager()
{
  //
  // Called at every step in the L3 magnet
  //
}
Commit	Line	Data
fe4da5cc	1	///////////////////////////////////////////////////////////////////////////////
	2	// //
	3	// L3 Magnet //
	4	// //
	5	//Begin_Html
	6	/*
	7	<img src="gif/AliMAGClass.gif">
	8	</pre>
	9	<br clear=left>
	10	<font size=+2 color=red>
	11	<p>The responsible person for this module is
	12	<a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
	13	</font>
	14	<pre>
	15	*/
	16	//End_Html
	17	// //
	18	///////////////////////////////////////////////////////////////////////////////
	19
	20	#include "AliMAG.h"
	21	#include "AliRun.h"
	22	#include "TGeant3.h"
	23
	24	ClassImp(AliMAG)
	25
	26	//_____________________________________________________________________________
	27	AliMAG::AliMAG()
	28	{
	29	//
	30	// Default constructor for L3 magnet
	31	//
	32	}
	33
	34	//_____________________________________________________________________________
	35	AliMAG::AliMAG(const char name, const char title)
	36	: AliDetector(name,title)
	37	{
	38	//
	39	// Standard constructor for L3 magnet
	40	//
	41	//Begin_Html
	42	/*
	43	<img src="gif/aliMAG.gif">
	44	*/
	45	//End_Html
	46
	47	SetMarkerColor(7);
	48	SetMarkerStyle(2);
	49	SetMarkerSize(0.4);
	50	}
	51
	52	//_____________________________________________________________________________
	53	void AliMAG::BuildGeometry()
	54	{
	55	//
	56	// Build ROOT TNode geometry for L3 magnet
	57	// only for detectors
	58	//
	59	}
	60
	61	//_____________________________________________________________________________
	62	void AliMAG::CreateGeometry()
	63	{
	64	//
65	// Create geometry for L3 magnet
66	//
67	//Begin_Html
68	/*
69	<img src="gif/mag.gif">
70	*/
71	//End_Html
72	//Begin_Html
73	/*
74	<img src="gif/tree_mag.gif">
75	*/
76	//End_Html
77
78	AliMC* pMC = AliMC::GetMC();
79
80	Int_t *idtmed = gAlice->Idtmed();
81
82	Float_t dpar[13];
83	Int_t idrotm[399];
84	Float_t par[10];
85
86	// ANGLE POLAIRE MAXIMUM
87
88	// Define Mother
89
90	par[0] = 22.5;
91	par[1] = 360.;
92	par[2] = 8.;
93	par[3] = 2.;
94	par[4] = -600.;
95	par[5] = 580.;
96	par[6] = 790.;
97	par[7] = 600.;
98	par[8] = 580.;
99	par[9] = 790.;
100	pMC->Gsvolu("L3MO", "PGON", idtmed[334], par, 10);
101	pMC->Gspos("L3MO", 1, "ALIC", 0., -30., 0., 0, "ONLY");
102
103	// Define coils
104
105	par[5] = 585.;
106	par[6] = 690.;
107	par[8] = 585.;
108	par[9] = 690.;
109	pMC->Gsvolu("L3CO", "PGON", idtmed[328], par, 10);
110	pMC->Gspos("L3CO", 1, "L3MO", 0., 0., 0., 0, "ONLY");
111
112	par[5] = 580.;
113	par[6] = 585.;
114	par[8] = 580.;
115	par[9] = 585.;
116	pMC->Gsvolu("L3C1", "PGON", idtmed[308], par, 10);
117	pMC->Gspos("L3C1", 1, "L3MO", 0., 0., 0., 0, "ONLY");
118
119	// Define yoke
120
121	par[5] = 690.;
122	par[6] = 790.;
123	par[8] = 690.;
124	par[9] = 790.;
125	pMC->Gsvolu("L3YO", "PGON", idtmed[329], par, 10);
126	pMC->Gspos("L3YO", 1, "L3MO", 0., 0., 0., 0, "ONLY");
127
128	// Define the return yoke of L3 (DOOR)
129
130	par[4] = 600.;
131	par[5] = 0.;
132	par[6] = 790.;
133	par[7] = 700.;
134	par[8] = par[5];
135	par[9] = par[6];
136	pMC->Gsvolu("L3DO", "PGON", idtmed[334], par, 10);
137
138	par[4] = 610.;
139	par[5] = 0.;
140	par[6] = 790.;
141	par[7] = 700.;
142	par[8] = par[5];
143	par[9] = par[6];
144	pMC->Gsvolu("L3FR", "PGON", idtmed[329], par, 10);
145
146	// INNER LAYER
147
148	par[4] = 600.;
149	par[7] = 610.;
150	pMC->Gsvolu("L3IR", "PGON", idtmed[309], par, 10);
151
152	// DOOR OPENING
153
154	dpar[0] = 22.5;
155	dpar[1] = 360.;
156	dpar[2] = 8.;
157	dpar[3] = 3.;
158	dpar[4] = 610.;
159	dpar[5] = 0.;
160	dpar[6] = 163.5;
161	dpar[7] = 670.;
162	dpar[8] = dpar[5];
163	dpar[9] = dpar[6];
164	dpar[10] = 700.;
165	dpar[11] = dpar[5];
166	dpar[12] = dpar[6] + 50.;
167	pMC->Gsvolu("L3O1", "PGON", idtmed[314], dpar, 13);
168	par[4] = 600.;
169	par[5] = 0.;
170	par[6] = 163.5;
171	par[7] = 610.;
172	par[8] = 0.;
173	par[9] = 163.5;
174	pMC->Gsvolu("L3O2", "PGON", idtmed[314], par, 10);
175
176	// THE DOOR OPENING HAS TO BE PLACED WITH 'MANY' SINCE THE REGION
177	// WILL CONTAIN A MUON CHAMBER, BEAM PIPE AND BEAM SHIELD
178	// PLACED WITH 'ONLY'.
179
180	pMC->Gspos("L3O1", 1, "L3FR", 0., 30., 0., 0, "MANY");
181	pMC->Gspos("L3O2", 1, "L3IR", 0., 30., 0., 0, "MANY");
182
183	pMC->Gspos("L3FR", 1, "L3DO", 0., 0., 0., 0, "MANY");
184	pMC->Gspos("L3IR", 1, "L3DO", 0., 0., 0., 0, "MANY");
185
186	pMC->Gspos("L3DO", 1, "ALIC", 0., -30., 0., 0, "MANY");
187	AliMatrix(idrotm[300], 90., 0., 90., 90., 180., 0.);
188	pMC->Gspos("L3DO", 2, "ALIC", 0., -30., 0., idrotm[300], "MANY");
189	}
190
191	//_____________________________________________________________________________
192	void AliMAG::CreateMaterials()
193	{
194	//
195	// Create materials for L3 magnet
196	//
197
198	Int_t ISXFLD = gAlice->Field()->Integ();
199	Float_t SXMGMX = gAlice->Field()->Max();
200
201	Float_t epsil, stmin, deemax, tmaxfd, stemax;
202	// --- Define the various materials for GEANT ---
203
204	// Aluminum
205	AliMaterial(9, "Al$", 26.98, 13., 2.7, 8.9, 37.2);
206	AliMaterial(29, "Al$", 26.98, 13., 2.7, 8.9, 37.2);
207
208	// Iron
209	AliMaterial(10, "Fe$", 55.85, 26., 7.87, 1.76, 17.1);
210	AliMaterial(30, "Fe$", 55.85, 26., 7.87, 1.76, 17.1);
211
212	// Air
213	AliMaterial(15, "Air$", 14.61, 7.3, .001205, 30420., 67500);
214	AliMaterial(35, "Air$", 14.61, 7.3, .001205, 30420., 67500);
215
216	// ****************
217	// Defines tracking media parameters.
218	// Les valeurs sont commentees pour laisser le defaut
219	// a GEANT (version 3-21, page CONS200), f.m.
220	epsil = .001; // Tracking precision,
221	stemax = -1.; // Maximum displacement for multiple scat
222	tmaxfd = -20.; // Maximum angle due to field deflection
223	deemax = -.3; // Maximum fractional energy loss, DLS
224	stmin = -.8;
225	// ***************
226
227	// IRON
228
229	AliMedium(310, "FE_C0 ", 10, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
230	AliMedium(330, "FE_C1 ", 30, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
231
232	// ALUMINUM
233
234	AliMedium(309, "ALU_C0 ", 9, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
235	AliMedium(329, "ALU_C1 ", 29, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
236
237	// AIR
238
239	AliMedium(315, "AIR_C0 ", 15, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
240	AliMedium(335, "AIR_C1 ", 35, 0, ISXFLD, SXMGMX, tmaxfd, stemax, deemax, epsil, stmin);
241	}
242
243	//_____________________________________________________________________________
244	void AliMAG::DrawDetector()
245	{
246	//
247	// Draw a shaded view of the L3 magnet
248	//
249	}
250
251	//_____________________________________________________________________________
252	void AliMAG::Init()
253	{
254	//
255	// Initialise L3 magnet after it has been built
256	Int_t i;
257	//
258	printf("\n");
259	for(i=0;i<35;i++) printf("*");
260	printf(" MAG_INIT ");
261	for(i=0;i<35;i++) printf("*");
262	printf("\n");
263	//
264	// Here the MAG initialisation code (if any!)
265	for(i=0;i<80;i++) printf("*");
266	printf("\n");
267	}
268
269	//_____________________________________________________________________________
270	void AliMAG::StepManager()
271	{
272	//
273	// Called at every step in the L3 magnet
274	//
275	}