[u/mrichter/AliRoot.git] / PHOS / AliPHOSv0.cxx

/////////////////////////////////////////////////////////
//  Manager and hits classes for set:PHOS version 0    //
/////////////////////////////////////////////////////////
 
// --- ROOT system ---
#include "TH1.h"
#include "TRandom.h"
#include "TFile.h"
#include "TTree.h"
#include "TBRIK.h"
#include "TNode.h"

// --- galice header files ---
#include "AliPHOSv0.h"
#include "AliRun.h"
#include "AliConst.h"
#include "AliMC.h" 

ClassImp(AliPHOSv0)

//______________________________________________________________________________


AliPHOSv0::AliPHOSv0() : AliPHOS()
{
  fIdSens=0;
}
 
//______________________________________________________________________________

AliPHOSv0::AliPHOSv0(const char *name, const char *title)
          : AliPHOS(name, title)
{
  fIdSens=0;
}
 
//___________________________________________
void AliPHOSv0::Init()
{
  AliMC* pMC = AliMC::GetMC();
  
  fIdSens=pMC->VolId("PXTL");
}

//___________________________________________
void AliPHOSv0::CreateGeometry()
{
// *** DEFINITION OF THE -0.25<y<0.25 TILTED GEOMETRY OF THE PHOS *** 
// ORIGIN    : NICK VAN EIJNDHOVEN 

  AliMC* pMC = AliMC::GetMC();
  
    Float_t pphi;
    Float_t r, dptcb[3], dpair[3], dphos[3], dpucp[3], dpasp[3], dpcpv[3];
    Float_t dpxtl[3];
    Float_t yo;
    Int_t idrotm[99];
    Float_t xp1, yp1, xp2, yp2;
    
    Int_t *idtmed = gAlice->Idtmed();

// --- Dimensions of PbWO4 crystal --- 
      const Float_t XTL_X=2.2;
      const Float_t XTL_Y=18.;
      const Float_t XTL_Z=2.2;
// --- Tyvek wrapper thickness 
      const Float_t PAP_THICK=0.01;
// --- CPV thickness --- 
      const Float_t CPV_Y=0.5;
// --- Polystyrene Foam Outer Cover dimensions --- 
      const Float_t FOC_X=214.6;
      const Float_t FOC_Y=80.;
      const Float_t FOC_Z=260.;
// --- Inner AIR volume dimensions --- 
      const Float_t AIR_X=206.;
      const Float_t AIR_Y=66.;
      const Float_t AIR_Z=244.;
// --- Tyvek Crystal Block dimensions --- 
      const Float_t TCB_X=198.;
      const Float_t TCB_Y=25.;
      const Float_t TCB_Z=234.;
// --- Upper Cooling Plate thickness --- 
      const Float_t UCP_Y=0.06;
// --- Al Support Plate thickness --- 
      const Float_t ASP_Y=10.;
//--- Distance from IP to Foam Outer Cover top plate (needs to be 447.) ---
      const Float_t FOC_R=467.;
//--- Distance from IP to Crystal Block top Surface (needs to be 460.) ---
      const Float_t CBS_R=480.;

// --- Dimensions of volumes --- 


// --- Define PHOS box volume, fill with Polystyrene foam --- 
    dphos[0] = FOC_X/2.;
    dphos[1] = FOC_Y/2.;
    dphos[2] = FOC_Z/2.;
    pMC->Gsvolu("PHOS", "BOX ", idtmed[703], dphos, 3);

// --- Define air-filled box, place inside PHOS --- 
    dpair[0] = AIR_X/2.;
    dpair[1] = AIR_Y/2.;
    dpair[2] = AIR_Z/2.;
    pMC->Gsvolu("PAIR", "BOX ", idtmed[798], dpair, 3);
    pMC->Gspos("PAIR", 1, "PHOS", 0., 0., 0., 0, "ONLY");

// --- Define Upper Cooling Panel --- 
// --- place it right behind upper foam plate --- 
    dpucp[0] = TCB_X/2.;
    dpucp[1] = UCP_Y/2.;
    dpucp[2] = TCB_Z/2.;
    pMC->Gsvolu("PUCP", "BOX ", idtmed[701], dpucp, 3);
    yo = (AIR_Y-UCP_Y)/2.;
    pMC->Gspos("PUCP", 1, "PAIR", 0., yo, 0., 0, "ONLY");

// --- Define Crystal Block, fill with Tyvek, position inside PAIR --- 
    dptcb[0] = TCB_X/2.;
    dptcb[1] = TCB_Y/2.;
    dptcb[2] = TCB_Z/2.;
    pMC->Gsvolu("PTCB", "BOX ", idtmed[702], dptcb, 3);
// --- Divide PTCB in X and Z directions -- 
    pMC->Gsdvn("PSEC", "PTCB", 11, 1);
    pMC->Gsdvn("PMOD", "PSEC", 13, 3);
    pMC->Gsdvn("PSTR", "PMOD", 8, 1);
    pMC->Gsdvn("PCEL", "PSTR", 8, 3);
    yo = (FOC_Y-TCB_Y)/2. -(CBS_R-FOC_R);
    pMC->Gspos("PTCB", 1, "PAIR", 0., yo, 0., 0, "ONLY");

// --- Define PbWO4 crystal volume, place inside PCEL --- 
    dpxtl[0] = XTL_X/2.;
    dpxtl[1] = XTL_Y/2.;
    dpxtl[2] = XTL_Z/2.;
    pMC->Gsvolu("PXTL", "BOX ", idtmed[699], dpxtl, 3);
    yo = (TCB_Y-XTL_Y)/2. - PAP_THICK;
    pMC->Gspos("PXTL", 1, "PCEL", 0., yo, 0., 0, "ONLY");

// --- Define Al Support Plate, position it inside PAIR --- 
// --- right beneath PTCB --- 
    dpasp[0] = AIR_X/2.;
    dpasp[1] = ASP_Y/2.;
    dpasp[2] = AIR_Z/2.;
    pMC->Gsvolu("PASP", "BOX ", idtmed[701], dpasp, 3);
    yo = (FOC_Y-ASP_Y)/2. - (CBS_R-FOC_R+TCB_Y);
    pMC->Gspos("PASP", 1, "PAIR", 0., yo, 0., 0, "ONLY");

// --- Define CPV volume, DON'T PLACE IT YET --- 
    dpcpv[0] = TCB_X/2.;
    dpcpv[1] = CPV_Y/2.;
    dpcpv[2] = TCB_Z/2.;
    pMC->Gsvolu("PCPV", "BOX ", idtmed[700], dpcpv, 3);
// --- Divide in X and Z direction (same way as PTCB) --- 
    pMC->Gsdvn("PCSE", "PCPV", 11, 1);
    pMC->Gsdvn("PCMO", "PCSE", 13, 3);
    pMC->Gsdvn("PCST", "PCMO", 8, 1);
    pMC->Gsdvn("PCCE", "PCST", 8, 3);

// --- Position various PHOS units in ALICE setup --- 
// --- PHOS itself first --- 
    r     = FOC_R+FOC_Y/2.;
    pphi  = TMath::ATan(FOC_X/(2.*FOC_R));
    xp1   = -r * TMath::Sin(pphi * 3.);
    yp1   = -r * TMath::Cos(pphi * 3.);
    xp2   = -r * TMath::Sin(pphi);
    yp2   = -r * TMath::Cos(pphi);
    pphi *= kRaddeg;
    AliMatrix(idrotm[0], 90.,-3*pphi, 90., 90-3*pphi, 0., 0.);
    AliMatrix(idrotm[1], 90.,  -pphi, 90., 90-pphi,   0., 0.);
    AliMatrix(idrotm[2], 90.,   pphi, 90., 90+pphi,   0., 0.);
    AliMatrix(idrotm[3], 90., 3*pphi, 90., 90+3*pphi, 0., 0.);
    pMC->Gspos("PHOS", 1, "ALIC", xp1, yp1, 0., idrotm[0], "ONLY");
    pMC->Gspos("PHOS", 2, "ALIC", xp2, yp2, 0., idrotm[1], "ONLY");
    pMC->Gspos("PHOS", 3, "ALIC",-xp2, yp2, 0., idrotm[2], "ONLY");
    pMC->Gspos("PHOS", 4, "ALIC",-xp1, yp1, 0., idrotm[3], "ONLY");

// --- Now position PCPV so that its plates are right on top of --- 
// --- corresponding PHOS supermodules (previously called cradles) --- 
    r    = FOC_R-CPV_Y/2.;
    pphi = TMath::ATan(FOC_X/(2.*FOC_R));
    xp1  = -r * TMath::Sin(pphi * 3.);
    yp1  = -r * TMath::Cos(pphi * 3.);
    xp2  = -r * TMath::Sin(pphi);
    yp2  = -r * TMath::Cos(pphi);
    pMC->Gspos("PCPV", 1, "ALIC", xp1, yp1, 0., idrotm[0], "ONLY");
    pMC->Gspos("PCPV", 2, "ALIC", xp2, yp2, 0., idrotm[1], "ONLY");
    pMC->Gspos("PCPV", 3, "ALIC",-xp2, yp2, 0., idrotm[2], "ONLY");
    pMC->Gspos("PCPV", 4, "ALIC",-xp1, yp1, 0., idrotm[3], "ONLY");

// --- Set modules seen without tree for drawings --- 
    pMC->Gsatt("PMOD", "SEEN", -2);
    pMC->Gsatt("PCMO", "SEEN", -2);
}
 
//___________________________________________
void AliPHOSv0::CreateMaterials()
{
// *** DEFINITION OF AVAILABLE PHOS MATERIALS *** 
// ORIGIN    : NICK VAN EIJNDHOVEN 

  AliMC* pMC = AliMC::GetMC();
  
    Int_t   ISXFLD = gAlice->Field()->Integ();
    Float_t SXMGMX = gAlice->Field()->Max();
    
// --- The PbWO4 crystals --- 
    Float_t ax[3] = { 207.19,183.85,16. };
    Float_t zx[3] = { 82.,74.,8. };
    Float_t wx[3] = { 1.,1.,4. };
    Float_t dx    = 8.28;
// --- The polysterene scintillator (CH) --- 
    Float_t ap[2] = { 12.011,1.00794 };
    Float_t zp[2] = { 6.,1. };
    Float_t wp[2] = { 1.,1. };
    Float_t dp    = 1.032;
// --- Tyvek (CnH2n) 
    Float_t at[2] = { 12.011,1.00794 };
    Float_t zt[2] = { 6.,1. };
    Float_t wt[2] = { 1.,2. };
    Float_t dt    = .331;
// --- Polystyrene foam --- 
    Float_t af[2] = { 12.011,1.00794 };
    Float_t zf[2] = { 6.,1. };
    Float_t wf[2] = { 1.,1. };
    Float_t df    = .3;

    Int_t *idtmed = gAlice->Idtmed();
    
    AliMixture( 0, "PbWO4$",       ax, zx, dx, -3, wx);
    AliMixture( 1, "Polystyrene$", ap, zp, dp, -2, wp);
    AliMaterial(2, "Al$",          26.98, 13., 2.7, 8.9, 999.);
// ---                                Absorption length^ is ignored --- 
    AliMixture( 3, "Tyvek$", at, zt, dt, -2, wt);
    AliMixture( 4, "Foam$",  af, zf, df, -2, wf);
    AliMaterial(9, "Air$", 14.61, 7.3, .001205, 30420., 67500);

    AliMedium(700, "PHOS Xtal    $", 0, 1, ISXFLD, SXMGMX, 10., .1, .1, .1, .1);
    AliMedium(701, "CPV scint.   $", 1, 1, ISXFLD, SXMGMX, 10., .1, .1, .1, .1);
    AliMedium(702, "Al parts     $", 2, 0, ISXFLD, SXMGMX, 10., .1, .1, .001, .001);
    AliMedium(703, "Tyvek wrapper$", 3, 0, ISXFLD, SXMGMX, 10., .1, .1, .001, .001);
    AliMedium(704, "Polyst. foam $", 4, 0, ISXFLD, SXMGMX, 10., .1, .1, .1, .1);
    AliMedium(799, "Air          $", 9, 0, ISXFLD, SXMGMX, 10., 1., .1, .1, 10.);

// --- Generate explicitly delta rays in aluminium parts --- 
    pMC->Gstpar(idtmed[701], "LOSS", 3.);
    pMC->Gstpar(idtmed[701], "DRAY", 1.);
}

void AliPHOSv0::StepManager()
{

  AliMC* pMC = AliMC::GetMC();
  
  TClonesArray &lhits = *fHits;
  Int_t copy, i;
  Int_t vol[5];
  Float_t hits[4];
  if(pMC->CurrentVol(0,copy) == fIdSens) {
    //
    //We are in the sensitive volume
    for(i=0;i<4;i++) {
      pMC->CurrentVolOff(i+1,0,copy);
      vol[4-i]=copy;
    }
    pMC->CurrentVolOff(7,0,copy);
    vol[0]=copy;
    pMC->TrackPosition(hits);
    hits[3]=pMC->Edep();
    new(lhits[fNhits++]) AliPHOShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
  }
}
Commit	Line	Data
fe4da5cc	1	/////////////////////////////////////////////////////////
	2	// Manager and hits classes for set:PHOS version 0 //
	3	/////////////////////////////////////////////////////////
	4
	5	// --- ROOT system ---
	6	#include "TH1.h"
	7	#include "TRandom.h"
	8	#include "TFile.h"
	9	#include "TTree.h"
	10	#include "TBRIK.h"
	11	#include "TNode.h"
	12
	13	// --- galice header files ---
	14	#include "AliPHOSv0.h"
	15	#include "AliRun.h"
	16	#include "AliConst.h"
	17	#include "AliMC.h"
	18
	19	ClassImp(AliPHOSv0)
	20
	21	//______________________________________________________________________________
	22
	23
	24	AliPHOSv0::AliPHOSv0() : AliPHOS()
	25	{
	26	fIdSens=0;
	27	}
	28
	29	//______________________________________________________________________________
	30
	31	AliPHOSv0::AliPHOSv0(const char name, const char title)
	32	: AliPHOS(name, title)
	33	{
	34	fIdSens=0;
	35	}
	36
	37	//___________________________________________
	38	void AliPHOSv0::Init()
	39	{
	40	AliMC* pMC = AliMC::GetMC();
	41
	42	fIdSens=pMC->VolId("PXTL");
	43	}
	44
	45	//___________________________________________
	46	void AliPHOSv0::CreateGeometry()
	47	{
	48	// * DEFINITION OF THE -0.25<y<0.25 TILTED GEOMETRY OF THE PHOS *
	49	// ORIGIN : NICK VAN EIJNDHOVEN
	50
	51	AliMC* pMC = AliMC::GetMC();
	52
	53	Float_t pphi;
	54	Float_t r, dptcb[3], dpair[3], dphos[3], dpucp[3], dpasp[3], dpcpv[3];
	55	Float_t dpxtl[3];
	56	Float_t yo;
	57	Int_t idrotm[99];
	58	Float_t xp1, yp1, xp2, yp2;
	59
	60	Int_t *idtmed = gAlice->Idtmed();
	61
	62	// --- Dimensions of PbWO4 crystal ---
	63	const Float_t XTL_X=2.2;
	64	const Float_t XTL_Y=18.;
65	const Float_t XTL_Z=2.2;
66	// --- Tyvek wrapper thickness
67	const Float_t PAP_THICK=0.01;
68	// --- CPV thickness ---
69	const Float_t CPV_Y=0.5;
70	// --- Polystyrene Foam Outer Cover dimensions ---
71	const Float_t FOC_X=214.6;
72	const Float_t FOC_Y=80.;
73	const Float_t FOC_Z=260.;
74	// --- Inner AIR volume dimensions ---
75	const Float_t AIR_X=206.;
76	const Float_t AIR_Y=66.;
77	const Float_t AIR_Z=244.;
78	// --- Tyvek Crystal Block dimensions ---
79	const Float_t TCB_X=198.;
80	const Float_t TCB_Y=25.;
81	const Float_t TCB_Z=234.;
82	// --- Upper Cooling Plate thickness ---
83	const Float_t UCP_Y=0.06;
84	// --- Al Support Plate thickness ---
85	const Float_t ASP_Y=10.;
86	//--- Distance from IP to Foam Outer Cover top plate (needs to be 447.) ---
87	const Float_t FOC_R=467.;
88	//--- Distance from IP to Crystal Block top Surface (needs to be 460.) ---
89	const Float_t CBS_R=480.;
90
91	// --- Dimensions of volumes ---
92
93
94	// --- Define PHOS box volume, fill with Polystyrene foam ---
95	dphos[0] = FOC_X/2.;
96	dphos[1] = FOC_Y/2.;
97	dphos[2] = FOC_Z/2.;
98	pMC->Gsvolu("PHOS", "BOX ", idtmed[703], dphos, 3);
99
100	// --- Define air-filled box, place inside PHOS ---
101	dpair[0] = AIR_X/2.;
102	dpair[1] = AIR_Y/2.;
103	dpair[2] = AIR_Z/2.;
104	pMC->Gsvolu("PAIR", "BOX ", idtmed[798], dpair, 3);
105	pMC->Gspos("PAIR", 1, "PHOS", 0., 0., 0., 0, "ONLY");
106
107	// --- Define Upper Cooling Panel ---
108	// --- place it right behind upper foam plate ---
109	dpucp[0] = TCB_X/2.;
110	dpucp[1] = UCP_Y/2.;
111	dpucp[2] = TCB_Z/2.;
112	pMC->Gsvolu("PUCP", "BOX ", idtmed[701], dpucp, 3);
113	yo = (AIR_Y-UCP_Y)/2.;
114	pMC->Gspos("PUCP", 1, "PAIR", 0., yo, 0., 0, "ONLY");
115
116	// --- Define Crystal Block, fill with Tyvek, position inside PAIR ---
117	dptcb[0] = TCB_X/2.;
118	dptcb[1] = TCB_Y/2.;
119	dptcb[2] = TCB_Z/2.;
120	pMC->Gsvolu("PTCB", "BOX ", idtmed[702], dptcb, 3);
121	// --- Divide PTCB in X and Z directions --
122	pMC->Gsdvn("PSEC", "PTCB", 11, 1);
123	pMC->Gsdvn("PMOD", "PSEC", 13, 3);
124	pMC->Gsdvn("PSTR", "PMOD", 8, 1);
125	pMC->Gsdvn("PCEL", "PSTR", 8, 3);
126	yo = (FOC_Y-TCB_Y)/2. -(CBS_R-FOC_R);
127	pMC->Gspos("PTCB", 1, "PAIR", 0., yo, 0., 0, "ONLY");
128
129	// --- Define PbWO4 crystal volume, place inside PCEL ---
130	dpxtl[0] = XTL_X/2.;
131	dpxtl[1] = XTL_Y/2.;
132	dpxtl[2] = XTL_Z/2.;
133	pMC->Gsvolu("PXTL", "BOX ", idtmed[699], dpxtl, 3);
134	yo = (TCB_Y-XTL_Y)/2. - PAP_THICK;
135	pMC->Gspos("PXTL", 1, "PCEL", 0., yo, 0., 0, "ONLY");
136
137	// --- Define Al Support Plate, position it inside PAIR ---
138	// --- right beneath PTCB ---
139	dpasp[0] = AIR_X/2.;
140	dpasp[1] = ASP_Y/2.;
141	dpasp[2] = AIR_Z/2.;
142	pMC->Gsvolu("PASP", "BOX ", idtmed[701], dpasp, 3);
143	yo = (FOC_Y-ASP_Y)/2. - (CBS_R-FOC_R+TCB_Y);
144	pMC->Gspos("PASP", 1, "PAIR", 0., yo, 0., 0, "ONLY");
145
146	// --- Define CPV volume, DON'T PLACE IT YET ---
147	dpcpv[0] = TCB_X/2.;
148	dpcpv[1] = CPV_Y/2.;
149	dpcpv[2] = TCB_Z/2.;
150	pMC->Gsvolu("PCPV", "BOX ", idtmed[700], dpcpv, 3);
151	// --- Divide in X and Z direction (same way as PTCB) ---
152	pMC->Gsdvn("PCSE", "PCPV", 11, 1);
153	pMC->Gsdvn("PCMO", "PCSE", 13, 3);
154	pMC->Gsdvn("PCST", "PCMO", 8, 1);
155	pMC->Gsdvn("PCCE", "PCST", 8, 3);
156
157	// --- Position various PHOS units in ALICE setup ---
158	// --- PHOS itself first ---
159	r = FOC_R+FOC_Y/2.;
160	pphi = TMath::ATan(FOC_X/(2.*FOC_R));
161	xp1 = -r * TMath::Sin(pphi * 3.);
162	yp1 = -r * TMath::Cos(pphi * 3.);
163	xp2 = -r * TMath::Sin(pphi);
164	yp2 = -r * TMath::Cos(pphi);
165	pphi *= kRaddeg;
166	AliMatrix(idrotm[0], 90.,-3pphi, 90., 90-3pphi, 0., 0.);
167	AliMatrix(idrotm[1], 90., -pphi, 90., 90-pphi, 0., 0.);
168	AliMatrix(idrotm[2], 90., pphi, 90., 90+pphi, 0., 0.);
169	AliMatrix(idrotm[3], 90., 3pphi, 90., 90+3pphi, 0., 0.);
170	pMC->Gspos("PHOS", 1, "ALIC", xp1, yp1, 0., idrotm[0], "ONLY");
171	pMC->Gspos("PHOS", 2, "ALIC", xp2, yp2, 0., idrotm[1], "ONLY");
172	pMC->Gspos("PHOS", 3, "ALIC",-xp2, yp2, 0., idrotm[2], "ONLY");
173	pMC->Gspos("PHOS", 4, "ALIC",-xp1, yp1, 0., idrotm[3], "ONLY");
174
175	// --- Now position PCPV so that its plates are right on top of ---
176	// --- corresponding PHOS supermodules (previously called cradles) ---
177	r = FOC_R-CPV_Y/2.;
178	pphi = TMath::ATan(FOC_X/(2.*FOC_R));
179	xp1 = -r * TMath::Sin(pphi * 3.);
180	yp1 = -r * TMath::Cos(pphi * 3.);
181	xp2 = -r * TMath::Sin(pphi);
182	yp2 = -r * TMath::Cos(pphi);
183	pMC->Gspos("PCPV", 1, "ALIC", xp1, yp1, 0., idrotm[0], "ONLY");
184	pMC->Gspos("PCPV", 2, "ALIC", xp2, yp2, 0., idrotm[1], "ONLY");
185	pMC->Gspos("PCPV", 3, "ALIC",-xp2, yp2, 0., idrotm[2], "ONLY");
186	pMC->Gspos("PCPV", 4, "ALIC",-xp1, yp1, 0., idrotm[3], "ONLY");
187
188	// --- Set modules seen without tree for drawings ---
189	pMC->Gsatt("PMOD", "SEEN", -2);
190	pMC->Gsatt("PCMO", "SEEN", -2);
191	}
192
193	//___________________________________________
194	void AliPHOSv0::CreateMaterials()
195	{
196	// * DEFINITION OF AVAILABLE PHOS MATERIALS *
197	// ORIGIN : NICK VAN EIJNDHOVEN
198
199	AliMC* pMC = AliMC::GetMC();
200
201	Int_t ISXFLD = gAlice->Field()->Integ();
202	Float_t SXMGMX = gAlice->Field()->Max();
203
204	// --- The PbWO4 crystals ---
205	Float_t ax[3] = { 207.19,183.85,16. };
206	Float_t zx[3] = { 82.,74.,8. };
207	Float_t wx[3] = { 1.,1.,4. };
208	Float_t dx = 8.28;
209	// --- The polysterene scintillator (CH) ---
210	Float_t ap[2] = { 12.011,1.00794 };
211	Float_t zp[2] = { 6.,1. };
212	Float_t wp[2] = { 1.,1. };
213	Float_t dp = 1.032;
214	// --- Tyvek (CnH2n)
215	Float_t at[2] = { 12.011,1.00794 };
216	Float_t zt[2] = { 6.,1. };
217	Float_t wt[2] = { 1.,2. };
218	Float_t dt = .331;
219	// --- Polystyrene foam ---
220	Float_t af[2] = { 12.011,1.00794 };
221	Float_t zf[2] = { 6.,1. };
222	Float_t wf[2] = { 1.,1. };
223	Float_t df = .3;
224
225	Int_t *idtmed = gAlice->Idtmed();
226
227	AliMixture( 0, "PbWO4$", ax, zx, dx, -3, wx);
228	AliMixture( 1, "Polystyrene$", ap, zp, dp, -2, wp);
229	AliMaterial(2, "Al$", 26.98, 13., 2.7, 8.9, 999.);
230	// --- Absorption length^ is ignored ---
231	AliMixture( 3, "Tyvek$", at, zt, dt, -2, wt);
232	AliMixture( 4, "Foam$", af, zf, df, -2, wf);
233	AliMaterial(9, "Air$", 14.61, 7.3, .001205, 30420., 67500);
234
235	AliMedium(700, "PHOS Xtal $", 0, 1, ISXFLD, SXMGMX, 10., .1, .1, .1, .1);
236	AliMedium(701, "CPV scint. $", 1, 1, ISXFLD, SXMGMX, 10., .1, .1, .1, .1);
237	AliMedium(702, "Al parts $", 2, 0, ISXFLD, SXMGMX, 10., .1, .1, .001, .001);
238	AliMedium(703, "Tyvek wrapper$", 3, 0, ISXFLD, SXMGMX, 10., .1, .1, .001, .001);
239	AliMedium(704, "Polyst. foam $", 4, 0, ISXFLD, SXMGMX, 10., .1, .1, .1, .1);
240	AliMedium(799, "Air $", 9, 0, ISXFLD, SXMGMX, 10., 1., .1, .1, 10.);
241
242	// --- Generate explicitly delta rays in aluminium parts ---
243	pMC->Gstpar(idtmed[701], "LOSS", 3.);
244	pMC->Gstpar(idtmed[701], "DRAY", 1.);
245	}
246
247	void AliPHOSv0::StepManager()
248	{
249
250	AliMC* pMC = AliMC::GetMC();
251
252	TClonesArray &lhits = *fHits;
253	Int_t copy, i;
254	Int_t vol[5];
255	Float_t hits[4];
256	if(pMC->CurrentVol(0,copy) == fIdSens) {
257	//
258	//We are in the sensitive volume
259	for(i=0;i<4;i++) {
260	pMC->CurrentVolOff(i+1,0,copy);
261	vol[4-i]=copy;
262	}
263	pMC->CurrentVolOff(7,0,copy);
264	vol[0]=copy;
265	pMC->TrackPosition(hits);
266	hits[3]=pMC->Edep();
267	new(lhits[fNhits++]) AliPHOShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
268	}
269	}