[u/mrichter/AliRoot.git] / STRUCT / AliABSOv0.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.                  *
 **************************************************************************/

/*
$Log$
Revision 1.6  2000/06/15 09:40:31  morsch
Obsolete typedef keyword removed

Revision 1.5  2000/06/12 19:39:01  morsch
New structure of beam pipe and heating jacket.

Revision 1.4  2000/04/03 08:13:40  fca
Introduce extra scope for non ANSI compliant C++ compilers

Revision 1.3  2000/01/18 17:49:56  morsch
Serious overlap of ABSM with shield corrected
Small error in ARPB parameters corrected

Revision 1.2  2000/01/13 11:23:59  morsch
Last layer of Pb outer angle corrected

Revision 1.1  2000/01/12 15:39:30  morsch
Standard version of ABSO

*/

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  Muon ABSOrber                                                            //
//  This class contains the description of the muon absorber geometry        //
//                                                                           //
//Begin_Html
/*
<img src="picts/AliABSOClass.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 "AliABSOv0.h"
#include "AliRun.h"
#include "AliMC.h"
#include "AliConst.h"

ClassImp(AliABSOv0)
 
//_____________________________________________________________________________
AliABSOv0::AliABSOv0()
{
  //
  // Default constructor
  //
}
 
//_____________________________________________________________________________
AliABSOv0::AliABSOv0(const char *name, const char *title)
       : AliABSO(name,title)
{
  //
  // Standard constructor
  //
  SetMarkerColor(7);
  SetMarkerStyle(2);
  SetMarkerSize(0.4);
}
 
//_____________________________________________________________________________
void AliABSOv0::CreateGeometry()
{
    //
    // Creation of the geometry of the muon absorber
    //
    //Begin_Html
    /*
      <img src="picts/AliABSOv0Tree.gif">
    */
    //End_Html
    //Begin_Html
    /*
      <img src="picts/AliABSOv0.gif">
    */
    //End_Html
    
    //
    //

    enum {kC=1605, kAl=1608, kFe=1609, kCu=1610, kW=1611, kPb=1612,
		  kNiCuW=1620, kVacuum=1615, kAir=1614, kConcrete=1616,
		  kPolyCH2=1617, kSteel=1609, kInsulation=1613};	  
    
    Int_t *idtmed = fIdtmed->GetArray()-1599;
    
    Float_t par[24], cpar[5], cpar0[5], pcpar[12], tpar[3], tpar0[3]; 
    Float_t dz;

#include "ABSOSHILConst.h"
#include "ABSOConst.h"
    Float_t dTube=0.1;                     // tube thickness
    Float_t dInsu=0.5;                     // insulation thickness
    Float_t dEnve=0.1;                     // protective envelope thickness
    Float_t dFree=0.5;                     // clearance thickness


// Mother volume and outer shielding: Pb
  par[0]  = 0.;
  par[1]  = 360.;
  par[2]  = 7.;
    
  par[3]  = -(zRear-zAbsStart)/2.;
  par[4]  = rAbs;
  par[5]  = zAbsStart * TMath::Tan(theta1);

  par[6]  = par[3]+(zNose-zAbsStart);
  par[7]  = rAbs;
  par[8]  = zNose * TMath::Tan(theta1);

  par[9]  = par[3]+(zConeTPC-zAbsStart);
  par[10] = rAbs;
  par[11] = par[8] + (par[9] - par[6]) * TMath::Tan(theta2);

  par[12]  = par[3]+(zOpen-zAbsStart);
  par[13] = rAbs;
  par[14] = par[11] + (par[12] - par[9]) * TMath::Tan(accMax);

  par[15] = par[3]+(zRear-dRear-zAbsStart);
  par[16] = rAbs   + (par[15] - par[12]) * TMath::Tan(thetaOpen1) ;
  par[17] = par[14] + (par[15] - par[12]) * TMath::Tan(accMax);

  par[18] = par[3]+(zRear-dRear-zAbsStart);
  par[19] = (zRear-dRear) * TMath::Tan(accMin);
  par[20] = par[14] + (par[18] - par[12]) * TMath::Tan(accMax);

  par[21] = -par[3];
  par[22] =  zRear* TMath::Tan(accMin);
  par[23] = par[20] + (par[21] - par[18]) * TMath::Tan(accMax);
  gMC->Gsvolu("ABSS", "PCON", idtmed[kPb], par, 24);
  { // Begin local scope for i
      for (Int_t i=4; i<18; i+=3) par[i]  = 0;
  } // End local scope for i
  gMC->Gsvolu("ABSM", "PCON", idtmed[kVacuum+40], par, 24);
  gMC->Gspos("ABSS", 1, "ABSM", 0., 0., 0., 0, "ONLY");

//
// Steel envelope
//
  par[4] = par[5] -dSteel;
  par[7] = par[8] -dSteel;
  par[10]= par[11]-dSteel;  
  par[13]= par[14]-dSteel;  
  par[16]= par[17]-dSteel;  
  par[19]= par[20]-dSteel;  
  par[22]= par[23]-dSteel;  
  gMC->Gsvolu("ABST", "PCON", idtmed[kSteel], par, 24);
  gMC->Gspos("ABST", 1, "ABSS", 0., 0., 0., 0, "ONLY");
//
// Polyethylene shield
// 
  cpar[0] = (zRear - zConeTPC) / 2.;
  cpar[1] = zConeTPC * TMath::Tan(accMax);
  cpar[2] = cpar[1] + dPoly;
  cpar[3] = zRear * TMath::Tan(accMax);
  cpar[4] = cpar[3] + dPoly;
  gMC->Gsvolu("APOL", "CONE", idtmed[kPolyCH2+40], cpar, 5);
  dz = (zRear-zAbsStart)/2.-cpar[0];
  gMC->Gspos("APOL", 1, "ABSS", 0., 0., dz, 0, "ONLY");

//
// Tungsten nose to protect TPC
// 
  cpar[0] = (zNose - zAbsStart) / 2.;
  cpar[1] = zAbsStart * TMath::Tan(accMax);
  cpar[2] = zAbsStart * TMath::Tan(theta1)-dSteel;
  cpar[3] = zNose * TMath::Tan(accMax);
  cpar[4] = zNose * TMath::Tan(theta1)-dSteel;
  gMC->Gsvolu("ANOS", "CONE", idtmed[kW], cpar, 5);
//
  dz = -(zRear-zAbsStart)/2.+cpar[0];
  gMC->Gspos("ANOS", 1, "ABSS", 0., 0., dz, 0, "ONLY");
//
// Tungsten inner shield
//
  Float_t zW=zTwoDeg+.1;
  Float_t dZ = zW+(zRear-dRear-zW)/2.;
  //
  pcpar[0]  = 0.;
  pcpar[1]  = 360.;
  pcpar[2]  = 3.;
  pcpar[3]  = zW-dZ;
  pcpar[4]  = rAbs;
  pcpar[5]  = zW * TMath::Tan(accMin);
  pcpar[6]  = zOpen-dZ;
  pcpar[7]  = rAbs;
  pcpar[8]  = zOpen * TMath::Tan(accMin);
  pcpar[9]  = zRear-dRear-dZ;
  pcpar[10] = rAbs+(zRear-dRear-zOpen) * TMath::Tan(thetaOpen1);
  pcpar[11] = (zRear-dRear) * TMath::Tan(accMin);
  
  gMC->Gsvolu("AWIN", "PCON", idtmed[kNiCuW+40], pcpar, 12);
  //
  dz=(zW+zRear-dRear)/2-(zAbsStart+zRear)/2.;
  gMC->Gspos("AWIN", 1, "ABSS", 0., 0., dz, 0, "ONLY");

  //     Inner tracking region
  //
  //     mother volume: Pb
  //
  pcpar[0]  = 0.;
  pcpar[1]  = 360.;
  pcpar[2]  = 3.;
  pcpar[3]  = -(zRear-zAbsStart)/2.;
  pcpar[4]  = rAbs;
  pcpar[5]  = zAbsStart * TMath::Tan(accMax);
  pcpar[6]  = pcpar[3]+(zTwoDeg-zAbsStart);
  pcpar[7]  = rAbs;
  pcpar[8]  = zTwoDeg * TMath::Tan(accMax);
  pcpar[9]  = -pcpar[3];
  pcpar[10] = zRear * TMath::Tan(accMin);
  pcpar[11] = zRear * TMath::Tan(accMax);
  gMC->Gsvolu("AITR", "PCON", idtmed[kPb], pcpar, 12);
  //
  // special Pb medium for last 5 cm of Pb
  Float_t zr=zRear-2.-0.001;
  cpar[0] = 1.0;
  cpar[1] = zr * TMath::Tan(thetaR);
  cpar[2] = zr * TMath::Tan(accMax);
  cpar[3] = cpar[1] + TMath::Tan(thetaR) * 2;
  cpar[4] = cpar[2] + TMath::Tan(accMax) * 2;
  gMC->Gsvolu("ARPB", "CONE", idtmed[kPb], cpar, 5);
  dz=(zRear-zAbsStart)/2.-cpar[0]-0.001;
  gMC->Gspos("ARPB", 1, "AITR", 0., 0., dz, 0, "ONLY");
  //
  //     concrete cone: concrete 
  //
  pcpar[9]  = pcpar[3]+(zRear-dRear-zAbsStart);
  pcpar[10] = (zRear-dRear) * TMath::Tan(accMin);
  pcpar[11] = (zRear-dRear) * TMath::Tan(accMax);
  gMC->Gsvolu("ACON", "PCON", idtmed[kConcrete+40], pcpar, 12);
  gMC->Gspos("ACON", 1, "AITR", 0., 0., 0., 0, "ONLY");
  //
  //     carbon cone: carbon
  //
  pcpar[9]  = pcpar[3]+(zAbsCc-zAbsStart);
  pcpar[10]  = zAbsCc * TMath::Tan(accMin);
  pcpar[11]  = zAbsCc * TMath::Tan(accMax);
  gMC->Gsvolu("ACAR", "PCON", idtmed[kC+40], pcpar, 12);
  gMC->Gspos("ACAR", 1, "ACON", 0., 0., 0., 0, "ONLY");
 //
 //     carbon cone outer region
 //
  cpar[0]  = 10.;
  cpar[1]  = rAbs;
  cpar[2]  = zAbsStart* TMath::Tan(accMax);
  cpar[3]  = rAbs;
  cpar[4]  = cpar[2]+2. * cpar[0] * TMath::Tan(accMax);

  gMC->Gsvolu("ACAO", "CONE", idtmed[kC], cpar, 5);
  dz=-(zRear-zAbsStart)/2.+cpar[0];
  gMC->Gspos("ACAO", 1, "ACAR", 0., 0., dz, 0, "ONLY");
  //
  //     inner W shield
  Float_t epsi=0.;
  Float_t repsi=1.;
  
  zr=zRear-(dRear-epsi);
  cpar[0] = (dRear-epsi)/2.;
  cpar[1] = zr * TMath::Tan(accMin);
  cpar[2] = zr * TMath::Tan(thetaR*repsi);
  cpar[3] = cpar[1] + TMath::Tan(accMin) * (dRear-epsi);
  cpar[4] = cpar[2] + TMath::Tan(thetaR*repsi) * (dRear-epsi);
  gMC->Gsvolu("ARW0", "CONE", idtmed[kNiCuW+40], cpar, 5);
  dz=(zRear-zAbsStart)/2.-cpar[0];
  gMC->Gspos("ARW0", 1, "AITR", 0., 0., dz, 0, "ONLY");
  //
  // special W medium for last 5 cm of W
  zr=zRear-5;
  cpar[0] = 2.5;
  cpar[1] = zr * TMath::Tan(accMin);
  cpar[2] = zr * TMath::Tan(thetaR*repsi);
  cpar[3] = cpar[1] + TMath::Tan(accMin) * 5.;
  cpar[4] = cpar[2] + TMath::Tan(thetaR*repsi) * 5.;
  gMC->Gsvolu("ARW1", "CONE", idtmed[kNiCuW+20], cpar, 5);
  dz=(dRear-epsi)/2.-cpar[0];
  gMC->Gspos("ARW1", 1, "ARW0", 0., 0., dz, 0, "ONLY");
  //
  // PolyEthylene Layers
  Float_t drMin=TMath::Tan(thetaR) * 5;
  Float_t drMax=TMath::Tan(accMax) * 5;
  gMC->Gsvolu("ARPE", "CONE", idtmed[kPolyCH2], cpar, 0);
  cpar[0]=2.5;
  { // Begin local scope for i
      for (Int_t i=0; i<3; i++) {
	  zr=zRear-dRear+5+i*10.;
	  cpar[1] = zr * TMath::Tan(thetaR);
	  cpar[2] = zr * TMath::Tan(accMax);
	  cpar[3] = cpar[1] + drMin;
	  cpar[4] = cpar[2] + drMax;
	  dz=(zRear-zAbsStart)/2.-cpar[0]-5.-(2-i)*10;
	  gMC->Gsposp("ARPE", i+1, "AITR", 0., 0., dz, 0, "ONLY",cpar,5);
      }
  } // End local scope for i
  gMC->Gspos("AITR", 1, "ABSS", 0., 0., 0., 0, "ONLY");	
  dz = (zRear-zAbsStart)/2.+zAbsStart;
  gMC->Gspos("ABSM", 1, "ALIC", 0., 0., dz, 0, "ONLY");	
//
//
// vacuum system
//
// pipe and heating jackets
//
//
// cylindrical piece
  tpar0[2]=(zOpen-zAbsStart)/2;
  tpar0[0]=rVacu;
  tpar0[1]=rAbs;
  gMC->Gsvolu("AV11", "TUBE", idtmed[kSteel+40], tpar0, 3);
//
// insulation

  tpar[2]=tpar0[2];
  tpar[0]=rVacu+dTube;
  tpar[1]=tpar[0]+dInsu;
  gMC->Gsvolu("AI11", "TUBE", idtmed[kInsulation+40], tpar, 3);
  gMC->Gspos("AI11", 1, "AV11", 0., 0., 0., 0, "ONLY"); 
//
// clearance 
  tpar[0]=tpar[1]+dEnve;
  tpar[1]=tpar[0]+dFree;
  gMC->Gsvolu("AP11", "TUBE", idtmed[kAir+40], tpar, 3);
  gMC->Gspos("AP11", 1, "AV11", 0., 0., 0., 0, "ONLY"); 
//
  dz=-(zRear-zAbsStart)/2.+tpar0[2];
  gMC->Gspos("AV11", 1, "ABSM", 0., 0., dz, 0, "ONLY"); 
//
// conical piece

  cpar0[0]=(zRear-dRear-zOpen)/2;
  cpar0[1]=rVacu-0.05;
  cpar0[2]=rAbs;
  Float_t dR=2.*cpar0[0]*TMath::Tan(thetaOpen1);
  cpar0[3]=cpar0[1]+dR;
  cpar0[4]=cpar0[2]+dR;
  gMC->Gsvolu("AV21", "CONE", idtmed[kSteel+40], cpar0, 5);
  dTube+=0.05;

//
// insulation
  cpar[0]=cpar0[0];
  cpar[1]=cpar0[1]+dTube;
  cpar[2]=cpar0[1]+dTube+dInsu;
  cpar[3]=cpar0[3]+dTube;
  cpar[4]=cpar0[3]+dTube+dInsu;
  gMC->Gsvolu("AI21", "CONE", idtmed[kInsulation+40], cpar, 5);
  gMC->Gspos("AI21", 1, "AV21", 0., 0., 0., 0, "ONLY"); 
//
// clearance
  cpar[1]=cpar0[1]+dTube+dInsu+dEnve;
  cpar[2]=rAbs;
  cpar[3]=cpar0[1]+dTube+dInsu+dEnve+dR;
  cpar[4]=rAbs+dR;

  gMC->Gsvolu("AP21", "CONE", idtmed[kAir+40], cpar, 5);
  gMC->Gspos("AP21", 1, "AV21", 0., 0., 0., 0, "ONLY"); 
  
  dz=(zRear-zAbsStart)/2.-cpar0[0]-dRear;
  gMC->Gspos("AV21", 1, "ABSM", 0., 0., dz, 0, "ONLY"); 
}

//_____________________________________________________________________________

void AliABSOv0::Init()
{
  //
  // Initialisation of the muon absorber after it has been built
  Int_t i;
  //
  printf("\n");
  for(i=0;i<35;i++) printf("*");
  printf(" ABSOv0_INIT ");
  for(i=0;i<35;i++) printf("*");
  printf("\n");
  //
  for(i=0;i<80;i++) printf("*");
  printf("\n");
}
Commit	Line	Data
43d014f0	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	/*
	17	$Log$
94de3818	18	Revision 1.6 2000/06/15 09:40:31 morsch
	19	Obsolete typedef keyword removed
	20
56f14ac2	21	Revision 1.5 2000/06/12 19:39:01 morsch
	22	New structure of beam pipe and heating jacket.
	23
b2ef5f06	24	Revision 1.4 2000/04/03 08:13:40 fca
	25	Introduce extra scope for non ANSI compliant C++ compilers
	26
74c0d076	27	Revision 1.3 2000/01/18 17:49:56 morsch
	28	Serious overlap of ABSM with shield corrected
	29	Small error in ARPB parameters corrected
	30
85d164ab	31	Revision 1.2 2000/01/13 11:23:59 morsch
	32	Last layer of Pb outer angle corrected
	33
54fba644	34	Revision 1.1 2000/01/12 15:39:30 morsch
b2ef5f06	35	Standard version of ABSO
54fba644	36
43d014f0	37	*/
	38
	39	///////////////////////////////////////////////////////////////////////////////
	40	// //
	41	// Muon ABSOrber //
	42	// This class contains the description of the muon absorber geometry //
	43	// //
	44	//Begin_Html
	45	/*
	46	<img src="picts/AliABSOClass.gif">
	47	</pre>
	48	<br clear=left>
	49	<font size=+2 color=red>
	50	<p>The responsible person for this module is
	51	<a href="mailto:andreas.morsch@cern.ch">Andreas Morsch</a>.
	52	</font>
	53	<pre>
	54	*/
	55	//End_Html
	56	// //
	57	// //
	58	///////////////////////////////////////////////////////////////////////////////
	59
	60	#include "AliABSOv0.h"
	61	#include "AliRun.h"
94de3818	62	#include "AliMC.h"
43d014f0	63	#include "AliConst.h"
b2ef5f06	64
43d014f0	65	ClassImp(AliABSOv0)
	66
	67	//_____________________________________________________________________________
	68	AliABSOv0::AliABSOv0()
	69	{
	70	//
	71	// Default constructor
	72	//
	73	}
	74
	75	//_____________________________________________________________________________
	76	AliABSOv0::AliABSOv0(const char name, const char title)
	77	: AliABSO(name,title)
	78	{
	79	//
	80	// Standard constructor
	81	//
	82	SetMarkerColor(7);
	83	SetMarkerStyle(2);
	84	SetMarkerSize(0.4);
	85	}
	86
	87	//_____________________________________________________________________________
	88	void AliABSOv0::CreateGeometry()
	89	{
b2ef5f06	90	//
	91	// Creation of the geometry of the muon absorber
	92	//
	93	//Begin_Html
	94	/*
	95	<img src="picts/AliABSOv0Tree.gif">
	96	*/
	97	//End_Html
	98	//Begin_Html
	99	/*
	100	<img src="picts/AliABSOv0.gif">
	101	*/
	102	//End_Html
	103
	104	//
	105	//
56f14ac2	106
56f14ac2	107	enum {kC=1605, kAl=1608, kFe=1609, kCu=1610, kW=1611, kPb=1612,
b2ef5f06	108	kNiCuW=1620, kVacuum=1615, kAir=1614, kConcrete=1616,
	109	kPolyCH2=1617, kSteel=1609, kInsulation=1613};
	110
	111	Int_t *idtmed = fIdtmed->GetArray()-1599;
43d014f0	112
b2ef5f06	113	Float_t par[24], cpar[5], cpar0[5], pcpar[12], tpar[3], tpar0[3];
	114	Float_t dz;
	115
	116	#include "ABSOSHILConst.h"
	117	#include "ABSOConst.h"
	118	Float_t dTube=0.1; // tube thickness
	119	Float_t dInsu=0.5; // insulation thickness
	120	Float_t dEnve=0.1; // protective envelope thickness
	121	Float_t dFree=0.5; // clearance thickness
	122
	123
	124	// Mother volume and outer shielding: Pb
43d014f0	125	par[0] = 0.;
	126	par[1] = 360.;
	127	par[2] = 7.;
b2ef5f06	128
	129	par[3] = -(zRear-zAbsStart)/2.;
	130	par[4] = rAbs;
	131	par[5] = zAbsStart * TMath::Tan(theta1);
43d014f0	132
b2ef5f06	133	par[6] = par[3]+(zNose-zAbsStart);
	134	par[7] = rAbs;
	135	par[8] = zNose * TMath::Tan(theta1);
43d014f0	136
b2ef5f06	137	par[9] = par[3]+(zConeTPC-zAbsStart);
b2ef5f06	138	par[10] = rAbs;
43d014f0	139	par[11] = par[8] + (par[9] - par[6]) * TMath::Tan(theta2);
43d014f0	140
b2ef5f06	141	par[12] = par[3]+(zOpen-zAbsStart);
	142	par[13] = rAbs;
	143	par[14] = par[11] + (par[12] - par[9]) * TMath::Tan(accMax);
43d014f0	144
b2ef5f06	145	par[15] = par[3]+(zRear-dRear-zAbsStart);
	146	par[16] = rAbs + (par[15] - par[12]) * TMath::Tan(thetaOpen1) ;
	147	par[17] = par[14] + (par[15] - par[12]) * TMath::Tan(accMax);
43d014f0	148
b2ef5f06	149	par[18] = par[3]+(zRear-dRear-zAbsStart);
	150	par[19] = (zRear-dRear) * TMath::Tan(accMin);
	151	par[20] = par[14] + (par[18] - par[12]) * TMath::Tan(accMax);
43d014f0	152
43d014f0	153	par[21] = -par[3];
b2ef5f06	154	par[22] = zRear* TMath::Tan(accMin);
	155	par[23] = par[20] + (par[21] - par[18]) * TMath::Tan(accMax);
	156	gMC->Gsvolu("ABSS", "PCON", idtmed[kPb], par, 24);
74c0d076	157	{ // Begin local scope for i
b2ef5f06	158	for (Int_t i=4; i<18; i+=3) par[i] = 0;
74c0d076	159	} // End local scope for i
b2ef5f06	160	gMC->Gsvolu("ABSM", "PCON", idtmed[kVacuum+40], par, 24);
43d014f0	161	gMC->Gspos("ABSS", 1, "ABSM", 0., 0., 0., 0, "ONLY");
	162
	163	//
	164	// Steel envelope
	165	//
b2ef5f06	166	par[4] = par[5] -dSteel;
	167	par[7] = par[8] -dSteel;
	168	par[10]= par[11]-dSteel;
	169	par[13]= par[14]-dSteel;
	170	par[16]= par[17]-dSteel;
	171	par[19]= par[20]-dSteel;
	172	par[22]= par[23]-dSteel;
	173	gMC->Gsvolu("ABST", "PCON", idtmed[kSteel], par, 24);
43d014f0	174	gMC->Gspos("ABST", 1, "ABSS", 0., 0., 0., 0, "ONLY");
	175	//
	176	// Polyethylene shield
	177	//
b2ef5f06	178	cpar[0] = (zRear - zConeTPC) / 2.;
	179	cpar[1] = zConeTPC * TMath::Tan(accMax);
	180	cpar[2] = cpar[1] + dPoly;
	181	cpar[3] = zRear * TMath::Tan(accMax);
	182	cpar[4] = cpar[3] + dPoly;
	183	gMC->Gsvolu("APOL", "CONE", idtmed[kPolyCH2+40], cpar, 5);
	184	dz = (zRear-zAbsStart)/2.-cpar[0];
43d014f0	185	gMC->Gspos("APOL", 1, "ABSS", 0., 0., dz, 0, "ONLY");
	186
	187	//
	188	// Tungsten nose to protect TPC
	189	//
b2ef5f06	190	cpar[0] = (zNose - zAbsStart) / 2.;
	191	cpar[1] = zAbsStart * TMath::Tan(accMax);
	192	cpar[2] = zAbsStart * TMath::Tan(theta1)-dSteel;
	193	cpar[3] = zNose * TMath::Tan(accMax);
	194	cpar[4] = zNose * TMath::Tan(theta1)-dSteel;
	195	gMC->Gsvolu("ANOS", "CONE", idtmed[kW], cpar, 5);
	196	//
	197	dz = -(zRear-zAbsStart)/2.+cpar[0];
43d014f0	198	gMC->Gspos("ANOS", 1, "ABSS", 0., 0., dz, 0, "ONLY");
	199	//
	200	// Tungsten inner shield
	201	//
b2ef5f06	202	Float_t zW=zTwoDeg+.1;
b2ef5f06	203	Float_t dZ = zW+(zRear-dRear-zW)/2.;
43d014f0	204	//
b2ef5f06	205	pcpar[0] = 0.;
	206	pcpar[1] = 360.;
	207	pcpar[2] = 3.;
	208	pcpar[3] = zW-dZ;
	209	pcpar[4] = rAbs;
	210	pcpar[5] = zW * TMath::Tan(accMin);
	211	pcpar[6] = zOpen-dZ;
	212	pcpar[7] = rAbs;
	213	pcpar[8] = zOpen * TMath::Tan(accMin);
	214	pcpar[9] = zRear-dRear-dZ;
	215	pcpar[10] = rAbs+(zRear-dRear-zOpen) * TMath::Tan(thetaOpen1);
	216	pcpar[11] = (zRear-dRear) * TMath::Tan(accMin);
	217
	218	gMC->Gsvolu("AWIN", "PCON", idtmed[kNiCuW+40], pcpar, 12);
	219	//
	220	dz=(zW+zRear-dRear)/2-(zAbsStart+zRear)/2.;
43d014f0	221	gMC->Gspos("AWIN", 1, "ABSS", 0., 0., dz, 0, "ONLY");
	222
	223	// Inner tracking region
	224	//
	225	// mother volume: Pb
	226	//
	227	pcpar[0] = 0.;
	228	pcpar[1] = 360.;
	229	pcpar[2] = 3.;
b2ef5f06	230	pcpar[3] = -(zRear-zAbsStart)/2.;
	231	pcpar[4] = rAbs;
	232	pcpar[5] = zAbsStart * TMath::Tan(accMax);
	233	pcpar[6] = pcpar[3]+(zTwoDeg-zAbsStart);
	234	pcpar[7] = rAbs;
	235	pcpar[8] = zTwoDeg * TMath::Tan(accMax);
85d164ab	236	pcpar[9] = -pcpar[3];
b2ef5f06	237	pcpar[10] = zRear * TMath::Tan(accMin);
	238	pcpar[11] = zRear * TMath::Tan(accMax);
	239	gMC->Gsvolu("AITR", "PCON", idtmed[kPb], pcpar, 12);
43d014f0	240	//
43d014f0	241	// special Pb medium for last 5 cm of Pb
b2ef5f06	242	Float_t zr=zRear-2.-0.001;
	243	cpar[0] = 1.0;
	244	cpar[1] = zr * TMath::Tan(thetaR);
	245	cpar[2] = zr * TMath::Tan(accMax);
	246	cpar[3] = cpar[1] + TMath::Tan(thetaR) * 2;
	247	cpar[4] = cpar[2] + TMath::Tan(accMax) * 2;
	248	gMC->Gsvolu("ARPB", "CONE", idtmed[kPb], cpar, 5);
	249	dz=(zRear-zAbsStart)/2.-cpar[0]-0.001;
43d014f0	250	gMC->Gspos("ARPB", 1, "AITR", 0., 0., dz, 0, "ONLY");
43d014f0	251	//
	252	// concrete cone: concrete
	253	//
b2ef5f06	254	pcpar[9] = pcpar[3]+(zRear-dRear-zAbsStart);
	255	pcpar[10] = (zRear-dRear) * TMath::Tan(accMin);
	256	pcpar[11] = (zRear-dRear) * TMath::Tan(accMax);
	257	gMC->Gsvolu("ACON", "PCON", idtmed[kConcrete+40], pcpar, 12);
43d014f0	258	gMC->Gspos("ACON", 1, "AITR", 0., 0., 0., 0, "ONLY");
	259	//
	260	// carbon cone: carbon
	261	//
b2ef5f06	262	pcpar[9] = pcpar[3]+(zAbsCc-zAbsStart);
	263	pcpar[10] = zAbsCc * TMath::Tan(accMin);
	264	pcpar[11] = zAbsCc * TMath::Tan(accMax);
	265	gMC->Gsvolu("ACAR", "PCON", idtmed[kC+40], pcpar, 12);
43d014f0	266	gMC->Gspos("ACAR", 1, "ACON", 0., 0., 0., 0, "ONLY");
b2ef5f06	267	//
	268	// carbon cone outer region
	269	//
	270	cpar[0] = 10.;
	271	cpar[1] = rAbs;
	272	cpar[2] = zAbsStart* TMath::Tan(accMax);
	273	cpar[3] = rAbs;
	274	cpar[4] = cpar[2]+2. * cpar[0] * TMath::Tan(accMax);
	275
	276	gMC->Gsvolu("ACAO", "CONE", idtmed[kC], cpar, 5);
	277	dz=-(zRear-zAbsStart)/2.+cpar[0];
	278	gMC->Gspos("ACAO", 1, "ACAR", 0., 0., dz, 0, "ONLY");
43d014f0	279	//
43d014f0	280	// inner W shield
b2ef5f06	281	Float_t epsi=0.;
	282	Float_t repsi=1.;
	283
	284	zr=zRear-(dRear-epsi);
	285	cpar[0] = (dRear-epsi)/2.;
	286	cpar[1] = zr * TMath::Tan(accMin);
	287	cpar[2] = zr * TMath::Tan(thetaR*repsi);
	288	cpar[3] = cpar[1] + TMath::Tan(accMin) * (dRear-epsi);
	289	cpar[4] = cpar[2] + TMath::Tan(thetaRrepsi) (dRear-epsi);
	290	gMC->Gsvolu("ARW0", "CONE", idtmed[kNiCuW+40], cpar, 5);
	291	dz=(zRear-zAbsStart)/2.-cpar[0];
43d014f0	292	gMC->Gspos("ARW0", 1, "AITR", 0., 0., dz, 0, "ONLY");
	293	//
	294	// special W medium for last 5 cm of W
b2ef5f06	295	zr=zRear-5;
43d014f0	296	cpar[0] = 2.5;
b2ef5f06	297	cpar[1] = zr * TMath::Tan(accMin);
	298	cpar[2] = zr * TMath::Tan(thetaR*repsi);
	299	cpar[3] = cpar[1] + TMath::Tan(accMin) * 5.;
	300	cpar[4] = cpar[2] + TMath::Tan(thetaRrepsi) 5.;
	301	gMC->Gsvolu("ARW1", "CONE", idtmed[kNiCuW+20], cpar, 5);
	302	dz=(dRear-epsi)/2.-cpar[0];
43d014f0	303	gMC->Gspos("ARW1", 1, "ARW0", 0., 0., dz, 0, "ONLY");
	304	//
	305	// PolyEthylene Layers
b2ef5f06	306	Float_t drMin=TMath::Tan(thetaR) * 5;
	307	Float_t drMax=TMath::Tan(accMax) * 5;
	308	gMC->Gsvolu("ARPE", "CONE", idtmed[kPolyCH2], cpar, 0);
43d014f0	309	cpar[0]=2.5;
74c0d076	310	{ // Begin local scope for i
b2ef5f06	311	for (Int_t i=0; i<3; i++) {
	312	zr=zRear-dRear+5+i*10.;
	313	cpar[1] = zr * TMath::Tan(thetaR);
	314	cpar[2] = zr * TMath::Tan(accMax);
	315	cpar[3] = cpar[1] + drMin;
	316	cpar[4] = cpar[2] + drMax;
	317	dz=(zRear-zAbsStart)/2.-cpar[0]-5.-(2-i)*10;
	318	gMC->Gsposp("ARPE", i+1, "AITR", 0., 0., dz, 0, "ONLY",cpar,5);
	319	}
74c0d076	320	} // End local scope for i
43d014f0	321	gMC->Gspos("AITR", 1, "ABSS", 0., 0., 0., 0, "ONLY");
b2ef5f06	322	dz = (zRear-zAbsStart)/2.+zAbsStart;
43d014f0	323	gMC->Gspos("ABSM", 1, "ALIC", 0., 0., dz, 0, "ONLY");
	324	//
	325	//
	326	// vacuum system
	327	//
	328	// pipe and heating jackets
	329	//
	330	//
	331	// cylindrical piece
b2ef5f06	332	tpar0[2]=(zOpen-zAbsStart)/2;
	333	tpar0[0]=rVacu;
	334	tpar0[1]=rAbs;
	335	gMC->Gsvolu("AV11", "TUBE", idtmed[kSteel+40], tpar0, 3);
43d014f0	336	//
43d014f0	337	// insulation
b2ef5f06	338
43d014f0	339	tpar[2]=tpar0[2];
b2ef5f06	340	tpar[0]=rVacu+dTube;
	341	tpar[1]=tpar[0]+dInsu;
	342	gMC->Gsvolu("AI11", "TUBE", idtmed[kInsulation+40], tpar, 3);
43d014f0	343	gMC->Gspos("AI11", 1, "AV11", 0., 0., 0., 0, "ONLY");
43d014f0	344	//
b2ef5f06	345	// clearance
	346	tpar[0]=tpar[1]+dEnve;
	347	tpar[1]=tpar[0]+dFree;
	348	gMC->Gsvolu("AP11", "TUBE", idtmed[kAir+40], tpar, 3);
43d014f0	349	gMC->Gspos("AP11", 1, "AV11", 0., 0., 0., 0, "ONLY");
b2ef5f06	350	//
b2ef5f06	351	dz=-(zRear-zAbsStart)/2.+tpar0[2];
43d014f0	352	gMC->Gspos("AV11", 1, "ABSM", 0., 0., dz, 0, "ONLY");
43d014f0	353	//
43d014f0	354	// conical piece
b2ef5f06	355
	356	cpar0[0]=(zRear-dRear-zOpen)/2;
	357	cpar0[1]=rVacu-0.05;
	358	cpar0[2]=rAbs;
	359	Float_t dR=2.cpar0[0]TMath::Tan(thetaOpen1);
	360	cpar0[3]=cpar0[1]+dR;
	361	cpar0[4]=cpar0[2]+dR;
	362	gMC->Gsvolu("AV21", "CONE", idtmed[kSteel+40], cpar0, 5);
	363	dTube+=0.05;
	364
43d014f0	365	//
	366	// insulation
	367	cpar[0]=cpar0[0];
b2ef5f06	368	cpar[1]=cpar0[1]+dTube;
	369	cpar[2]=cpar0[1]+dTube+dInsu;
	370	cpar[3]=cpar0[3]+dTube;
	371	cpar[4]=cpar0[3]+dTube+dInsu;
	372	gMC->Gsvolu("AI21", "CONE", idtmed[kInsulation+40], cpar, 5);
43d014f0	373	gMC->Gspos("AI21", 1, "AV21", 0., 0., 0., 0, "ONLY");
	374	//
	375	// clearance
b2ef5f06	376	cpar[1]=cpar0[1]+dTube+dInsu+dEnve;
	377	cpar[2]=rAbs;
	378	cpar[3]=cpar0[1]+dTube+dInsu+dEnve+dR;
	379	cpar[4]=rAbs+dR;
	380
	381	gMC->Gsvolu("AP21", "CONE", idtmed[kAir+40], cpar, 5);
43d014f0	382	gMC->Gspos("AP21", 1, "AV21", 0., 0., 0., 0, "ONLY");
43d014f0	383
b2ef5f06	384	dz=(zRear-zAbsStart)/2.-cpar0[0]-dRear;
43d014f0	385	gMC->Gspos("AV21", 1, "ABSM", 0., 0., dz, 0, "ONLY");
43d014f0	386	}
	387
	388	//_____________________________________________________________________________
	389
	390	void AliABSOv0::Init()
	391	{
	392	//
	393	// Initialisation of the muon absorber after it has been built
	394	Int_t i;
	395	//
	396	printf("\n");
	397	for(i=0;i<35;i++) printf("*");
	398	printf(" ABSOv0_INIT ");
	399	for(i=0;i<35;i++) printf("*");
	400	printf("\n");
	401	//
	402	for(i=0;i<80;i++) printf("*");
	403	printf("\n");
	404	}
	405
	406
	407
	408
	409
	410
	411
	412
	413
	414