[u/mrichter/AliRoot.git] / FIT / AliFITv0.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: AliT0v1.cxx 50646 2011-07-18 14:40:16Z alla $ */

/////////////////////////////////////////////////////////////////////
//                                                                 //
// T0 ( T-zero) detector  version 1                             //
//
//Begin Html       
/*
<img src="gif/AliT0v1Class.gif">
*/
//End Html
//                                                                  //
//                                                                  //
//////////////////////////////////////////////////////////////////////

#include <Riostream.h>
#include <stdlib.h>

#include "TGeoCompositeShape.h"
#include "TGeoManager.h"
#include "TGeoMatrix.h"
#include "TGeoVolume.h"
#include "TGeoTube.h"
#include "TGeoBBox.h"
#include "TGeoNode.h"


#include <TGeoGlobalMagField.h>
#include <TGraph.h>
#include <TLorentzVector.h>
#include <TMath.h>
#include <TVirtualMC.h>
#include <TString.h>

#include "AliLog.h"
#include "AliMagF.h"
#include "AliRun.h"

#include "AliFITHits.h"
#include "AliFITv0.h"

#include "AliMC.h"
#include "AliCDBLocal.h"
#include "AliCDBStorage.h"
#include "AliCDBManager.h"
#include "AliCDBEntry.h"
#include "AliTrackReference.h"

ClassImp(AliFITv0)


//--------------------------------------------------------------------
AliFITv0::AliFITv0():  AliFIT(),
		     fIdSens1(0)

{
  //
  // Standart constructor for T0 Detector version 0
}
//--------------------------------------------------------------------
AliFITv0::AliFITv0(const char *name, const char *title):
  AliFIT(name,title),
  fIdSens1(0)
{
  //
  // Standart constructor for T0 Detector version 0
  //
  printf("@@@@@@@@@@@@@@@ AliFITv0::AliFITv0 \n");
  fIshunt = 2; 
}
//_____________________________________________________________________________

AliFITv0::~AliFITv0() 
{
  // desctructor  
}

//-------------------------------------------------------------------------
void AliFITv0::CreateGeometry()
{
  //
  // Create the geometry of FIT Detector version 1 full geometry
  //
  // begin Html
  //

  printf("@@@@@@@@@@@ AliFITv0::CreateGeometry\n");
  Int_t *idtmed = fIdtmed->GetArray();
  /*
  TGeoMedium* kMedAir = gGeoManager->GetMedium("T0_AIR");
  TGeoMedium* kMedMCPGlass = gGeoManager->GetMedium("T0_glass");
  TGeoMedium* kMedOptGlass = gGeoManager->GetMedium("T0_OpticalGlass");
  TGeoMedium* kMedOptGlassCathode = gGeoManager->GetMedium("T0_OpticalGlassCathode");
  */
  Float_t zdetC = 80;
  Float_t zdetA = 373.;
  Int_t idrotm[999];
  Double_t x,y,z;
  Float_t pstart[3] = {6, 20 ,2.6};
  Float_t pinstart[3] = {3,3,2.55};
  Float_t  pmcp[3] = {2.95, 2.95, 1.5}; //MCP
  Float_t ptop[3] = {1.324, 1.324, 1.};//cherenkov radiator
  Float_t preg[3] = {1.324, 1.324, 0.05};//photcathode 
 
  AliMatrix(idrotm[901], 90., 0., 90., 90., 180., 0.);
  
  //-------------------------------------------------------------------
  //  T0 volume 
  //-------------------------------------------------------------------
  
  Float_t x1[20] = {9,  9, 15 ,15 , 9,  
		    3, -3,  3, -3, -9, 
		    -9, -9, -15, -15, -9, 
		    -3, 3, -3, 3, 9}; 
  
  Float_t y1[20] = {3.2, -3.2, 3.2, -3.2, -9.2, 
		    -9, -9, -15, -15, -9.2,
		    -3.2, 3.2, -3.2, 3.2, 9.2,
		    9, 9, 15, 15, 9.2};
  
  
  //mother tube
   TVirtualMC::GetMC()->Gsvolu("0STR","TUBE",idtmed[kAir],pstart,18);
   TVirtualMC::GetMC()->Gspos("0STR",1,"ALIC",0.,0.,-zdetC-pstart[2],idrotm[901],"ONLY");
   TVirtualMC::GetMC()->Gspos("0STR",2,"ALIC",0.,0.,zdetA+pstart[2],0,"ONLY");
   
   //T0 interior
   TVirtualMC::GetMC()->Gsvolu("0INS","BOX",idtmed[kAir],pinstart,3);
    z=-pstart[2]+pinstart[2];
   for (Int_t is=0; is<20; is++) {
     TVirtualMC::GetMC()->Gspos ("0INS", is + 1, "0STR", x1[is], y1[is], z, 0, "ONLY");
     TVirtualMC::GetMC()->Gspos ("0INS", is + 21, "0STR", x1[is], y1[is], z, 0, "ONLY");
     printf(" 0INS is %i x %f y %f z %f \n",is, x1[is],y1[is], z);
   }
   // 
  x=y=0;
  // Entry window (glass)
  TVirtualMC::GetMC()->Gsvolu("0TOP","BOX",idtmed[kAir],ptop,3); //glass
  TVirtualMC::GetMC()->Gsvolu ("0REG", "BOX", idtmed[kSensAir], preg, 3); 
  TVirtualMC::GetMC()->Gsvolu("0MCP","BOX",idtmed[kGlass],pmcp,3); //glass
  Int_t ntops=0;
  Float_t xin=0, yin=0;
  for (Int_t ix=0; ix<4; ix++) {
    xin = - pinstart[0] + 0.35 + (ix+0.5)*2*ptop[0] ;
    for (Int_t iy=0; iy<4; iy++) {
      z = - pinstart[2]+ptop[2];
      yin = - pinstart[1] + 0.35 + (iy+0.5)*2*ptop[1];
      ntops++;
      TVirtualMC::GetMC()->Gspos("0TOP",ntops,"0INS",xin,yin,z,0,"ONLY");
     //   printf(" 0TOP  full x %f y %f z %f \n", xin, yin, z);
      z = -pinstart[2] + 2 * ptop[2] + preg[2];
      TVirtualMC::GetMC()->Gspos ("0REG",ntops, "0INS", xin, yin, z, 0, "ONLY");
      printf(" GEOGEO  %i %i %i %f %f %f %f %f %f", ntops, ix, iy,
	     xin,yin,x1[ntops],y1[ntops],x1[ntops]+xin,y1[ntops]+yin);
    }
  }
  
 // MCP
 //  TGeoVolume* mcp =  gGeoManager->MakeBox("0MCP",kMedMCPGlass, 2.95, 2.95, 1.5);  
   z=-pinstart[2] + 2*ptop[2] + 2*preg[2] + pmcp[2];
   TVirtualMC::GetMC()->Gspos("0MCP",1,"0INS",0,0,z,0,"ONLY");
    
}    
//------------------------------------------------------------------------
void AliFITv0::AddAlignableVolumes() const
{
  //
  // Create entries for alignable volumes associating the symbolic volume
  // name with the corresponding volume path. Needs to be syncronized with
  // eventual changes in the geometry.
  //
  printf("@@@@@@@@@@@@@@@AliFITv0::AddAlignableVolumes()\n");
  TString volPath;
  TString symName, sn;
  TString vpAalign = "/ALIC_1/0STR_1";
  TString vpCalign = "/ALIC_1/0STR_2";
  for (Int_t imod=0; imod<2; imod++)  {
    if (imod==0) {volPath  = vpCalign; symName="/ALIC_1/0STR_1"; }
    if (imod==1) {volPath  = vpAalign; symName="/ALIC_1/0STR_2"; }
    
    AliDebug(2,"--------------------------------------------");
    AliDebug(2,Form("volPath=%s\n",volPath.Data()));
    AliDebug(2,Form("symName=%s\n",symName.Data()));
    AliDebug(2,"--------------------------------------------");
    if(!gGeoManager->SetAlignableEntry(symName.Data(),volPath.Data()))
      AliFatal(Form("Alignable entry %s not created. Volume path %s not valid",
symName.Data(),volPath.Data()));
  }
}   
//------------------------------------------------------------------------
void AliFITv0::CreateMaterials()
{

  printf("@@@@@@@@@@@@AliFITv0::CreateMaterials\n"); 
   Int_t isxfld   = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
   Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
   //   Float_t a,z,d,radl,absl,buf[1];
   // Int_t nbuf;
// AIR
                                                                                
   Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
   Float_t zAir[4]={6.,7.,8.,18.};
   Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
   Float_t dAir = 1.20479E-3;
   Float_t dAir1 = 1.20479E-11;
                                                                               
 // PMT glass SiO2
   Float_t aglass[2]={28.0855,15.9994};
   Float_t zglass[2]={14.,8.};
   Float_t wglass[2]={1.,2.};
   Float_t dglass=2.65;

    
//*** Definition Of avaible T0 materials ***
   AliMixture(1, "Vacuum$", aAir, zAir, dAir1,4,wAir);
   AliMixture(2, "Air$", aAir, zAir, dAir,4,wAir);
   AliMixture( 4, "PMT glass   $",aglass,zglass,dglass,-2,wglass);
   
   AliMedium(1, "FIT_Air$", 2, 0, isxfld, sxmgmx, 10., .1, 1., .003, .003);
   AliMedium(22, "FIT_AirSens$", 2, 1, isxfld, sxmgmx, 10., .1, 1., .003, .003);
   AliMedium(3, "FIT_Vacuum$", 1, 0, isxfld, sxmgmx, 10., .01, .1, .003, .003);
   AliMedium(6, "Glass$", 4, 1, isxfld, sxmgmx, 10., .01, .1, .003, .003);
    
   AliDebugClass(1,": ++++++++++++++Medium set++++++++++");
   
   
}


//-------------------------------------------------------------------
void AliFITv0::Init()
{
  // Initialises version 0 of the Forward Multiplicity Detector
  //
  AliFIT::Init();
  fIdSens1=TVirtualMC::GetMC()->VolId("0REG");
  
  AliDebug(1,Form("%s: *** FIT version 0 initialized ***\n",ClassName()));
}

//-------------------------------------------------------------------

void AliFITv0::StepManager()
{
  //
  // Called for every step in the T0 Detector
  //
  Int_t id,copy,copy1;
  static Float_t hits[6];
  static Int_t vol[3];
  TLorentzVector pos;
  TLorentzVector mom;
  
  //   TClonesArray &lhits = *fHits;
  
  if(!TVirtualMC::GetMC()->IsTrackAlive()) return; // particle has disappeared
  
  id=TVirtualMC::GetMC()->CurrentVolID(copy);  
  // Check the sensetive volume
  if(id==fIdSens1 ) { 
    if(TVirtualMC::GetMC()->IsTrackEntering()) {
      TVirtualMC::GetMC()->CurrentVolOffID(1,copy1);
      vol[1] = copy1;
      vol[0]=copy;
      TVirtualMC::GetMC()->TrackPosition(pos);
      hits[0] = pos[0];
      hits[1] = pos[1];
      hits[2] = pos[2];
      if(pos[2]<0) vol[2] = 0;
      else vol[2] = 1 ;
      printf(" volumes pmt %i mcp %i side %i x %f y %f z %f\n",  vol[0], vol[1],  vol[2], hits[0], hits[1], hits[2] );
      
      Float_t etot=TVirtualMC::GetMC()->Etot();
      hits[3]=etot;
      Int_t iPart= TVirtualMC::GetMC()->TrackPid();
      Int_t partID=TVirtualMC::GetMC()->IdFromPDG(iPart);
      hits[4]=partID;
      Float_t ttime=TVirtualMC::GetMC()->TrackTime();
      hits[5]=ttime*1e12;
      AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(),vol,hits);
      // Create a track reference at the exit of photocatode
    }
    
    //charge particle 
    if ( TVirtualMC::GetMC()->TrackCharge() )
      AddTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber(), AliTrackReference::kFIT);
  } //sensitive

}
Commit	Line	Data
c1c44db3	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: AliT0v1.cxx 50646 2011-07-18 14:40:16Z alla $ */
	17
	18	/////////////////////////////////////////////////////////////////////
	19	// //
	20	// T0 ( T-zero) detector version 1 //
	21	//
	22	//Begin Html
	23	/*
	24	<img src="gif/AliT0v1Class.gif">
	25	*/
	26	//End Html
	27	// //
	28	// //
	29	//////////////////////////////////////////////////////////////////////
	30
	31	#include <Riostream.h>
	32	#include <stdlib.h>
	33
	34	#include "TGeoCompositeShape.h"
	35	#include "TGeoManager.h"
	36	#include "TGeoMatrix.h"
	37	#include "TGeoVolume.h"
	38	#include "TGeoTube.h"
	39	#include "TGeoBBox.h"
	40	#include "TGeoNode.h"
	41
	42
	43	#include <TGeoGlobalMagField.h>
	44	#include <TGraph.h>
	45	#include <TLorentzVector.h>
	46	#include <TMath.h>
	47	#include <TVirtualMC.h>
	48	#include <TString.h>
	49
	50	#include "AliLog.h"
	51	#include "AliMagF.h"
	52	#include "AliRun.h"
	53
	54	#include "AliFITHits.h"
	55	#include "AliFITv0.h"
	56
	57	#include "AliMC.h"
	58	#include "AliCDBLocal.h"
	59	#include "AliCDBStorage.h"
	60	#include "AliCDBManager.h"
	61	#include "AliCDBEntry.h"
	62	#include "AliTrackReference.h"
	63
	64	ClassImp(AliFITv0)
65
66
67	//--------------------------------------------------------------------
68	AliFITv0::AliFITv0(): AliFIT(),
69	fIdSens1(0)
70
71	{
72	//
73	// Standart constructor for T0 Detector version 0
74	}
75	//--------------------------------------------------------------------
76	AliFITv0::AliFITv0(const char name, const char title):
77	AliFIT(name,title),
78	fIdSens1(0)
79	{
80	//
81	// Standart constructor for T0 Detector version 0
82	//
83	printf("@@@@@@@@@@@@@@@ AliFITv0::AliFITv0 \n");
84	fIshunt = 2;
85	}
86	//_____________________________________________________________________________
87
88	AliFITv0::~AliFITv0()
89	{
90	// desctructor
91	}
92
93	//-------------------------------------------------------------------------
94	void AliFITv0::CreateGeometry()
95	{
96	//
97	// Create the geometry of FIT Detector version 1 full geometry
98	//
99	// begin Html
100	//
101
102	printf("@@@@@@@@@@@ AliFITv0::CreateGeometry\n");
103	Int_t *idtmed = fIdtmed->GetArray();
104	/*
105	TGeoMedium* kMedAir = gGeoManager->GetMedium("T0_AIR");
106	TGeoMedium* kMedMCPGlass = gGeoManager->GetMedium("T0_glass");
107	TGeoMedium* kMedOptGlass = gGeoManager->GetMedium("T0_OpticalGlass");
108	TGeoMedium* kMedOptGlassCathode = gGeoManager->GetMedium("T0_OpticalGlassCathode");
109	*/
110	Float_t zdetC = 80;
111	Float_t zdetA = 373.;
112	Int_t idrotm[999];
113	Double_t x,y,z;
114	Float_t pstart[3] = {6, 20 ,2.6};
115	Float_t pinstart[3] = {3,3,2.55};
116	Float_t pmcp[3] = {2.95, 2.95, 1.5}; //MCP
117	Float_t ptop[3] = {1.324, 1.324, 1.};//cherenkov radiator
118	Float_t preg[3] = {1.324, 1.324, 0.05};//photcathode
119
120	AliMatrix(idrotm[901], 90., 0., 90., 90., 180., 0.);
121
122	//-------------------------------------------------------------------
123	// T0 volume
124	//-------------------------------------------------------------------
125
126	Float_t x1[20] = {9, 9, 15 ,15 , 9,
127	3, -3, 3, -3, -9,
128	-9, -9, -15, -15, -9,
129	-3, 3, -3, 3, 9};
130
131	Float_t y1[20] = {3.2, -3.2, 3.2, -3.2, -9.2,
132	-9, -9, -15, -15, -9.2,
133	-3.2, 3.2, -3.2, 3.2, 9.2,
134	9, 9, 15, 15, 9.2};
135
136
137	//mother tube
2942f542	138	TVirtualMC::GetMC()->Gsvolu("0STR","TUBE",idtmed[kAir],pstart,18);
	139	TVirtualMC::GetMC()->Gspos("0STR",1,"ALIC",0.,0.,-zdetC-pstart[2],idrotm[901],"ONLY");
	140	TVirtualMC::GetMC()->Gspos("0STR",2,"ALIC",0.,0.,zdetA+pstart[2],0,"ONLY");
c1c44db3	141
c1c44db3	142	//T0 interior
2942f542	143	TVirtualMC::GetMC()->Gsvolu("0INS","BOX",idtmed[kAir],pinstart,3);
c1c44db3	144	z=-pstart[2]+pinstart[2];
c1c44db3	145	for (Int_t is=0; is<20; is++) {
2942f542	146	TVirtualMC::GetMC()->Gspos ("0INS", is + 1, "0STR", x1[is], y1[is], z, 0, "ONLY");
2942f542	147	TVirtualMC::GetMC()->Gspos ("0INS", is + 21, "0STR", x1[is], y1[is], z, 0, "ONLY");
c1c44db3	148	printf(" 0INS is %i x %f y %f z %f \n",is, x1[is],y1[is], z);
	149	}
	150	//
	151	x=y=0;
	152	// Entry window (glass)
2942f542	153	TVirtualMC::GetMC()->Gsvolu("0TOP","BOX",idtmed[kAir],ptop,3); //glass
	154	TVirtualMC::GetMC()->Gsvolu ("0REG", "BOX", idtmed[kSensAir], preg, 3);
	155	TVirtualMC::GetMC()->Gsvolu("0MCP","BOX",idtmed[kGlass],pmcp,3); //glass
c1c44db3	156	Int_t ntops=0;
	157	Float_t xin=0, yin=0;
	158	for (Int_t ix=0; ix<4; ix++) {
	159	xin = - pinstart[0] + 0.35 + (ix+0.5)2ptop[0] ;
	160	for (Int_t iy=0; iy<4; iy++) {
	161	z = - pinstart[2]+ptop[2];
	162	yin = - pinstart[1] + 0.35 + (iy+0.5)2ptop[1];
	163	ntops++;
2942f542	164	TVirtualMC::GetMC()->Gspos("0TOP",ntops,"0INS",xin,yin,z,0,"ONLY");
c1c44db3	165	// printf(" 0TOP full x %f y %f z %f \n", xin, yin, z);
c1c44db3	166	z = -pinstart[2] + 2 * ptop[2] + preg[2];
2942f542	167	TVirtualMC::GetMC()->Gspos ("0REG",ntops, "0INS", xin, yin, z, 0, "ONLY");
c1c44db3	168	printf(" GEOGEO %i %i %i %f %f %f %f %f %f", ntops, ix, iy,
	169	xin,yin,x1[ntops],y1[ntops],x1[ntops]+xin,y1[ntops]+yin);
	170	}
	171	}
	172
	173	// MCP
	174	// TGeoVolume* mcp = gGeoManager->MakeBox("0MCP",kMedMCPGlass, 2.95, 2.95, 1.5);
	175	z=-pinstart[2] + 2ptop[2] + 2preg[2] + pmcp[2];
2942f542	176	TVirtualMC::GetMC()->Gspos("0MCP",1,"0INS",0,0,z,0,"ONLY");
c1c44db3	177
	178	}
	179	//------------------------------------------------------------------------
	180	void AliFITv0::AddAlignableVolumes() const
	181	{
	182	//
	183	// Create entries for alignable volumes associating the symbolic volume
	184	// name with the corresponding volume path. Needs to be syncronized with
	185	// eventual changes in the geometry.
	186	//
	187	printf("@@@@@@@@@@@@@@@AliFITv0::AddAlignableVolumes()\n");
	188	TString volPath;
	189	TString symName, sn;
	190	TString vpAalign = "/ALIC_1/0STR_1";
	191	TString vpCalign = "/ALIC_1/0STR_2";
	192	for (Int_t imod=0; imod<2; imod++) {
	193	if (imod==0) {volPath = vpCalign; symName="/ALIC_1/0STR_1"; }
	194	if (imod==1) {volPath = vpAalign; symName="/ALIC_1/0STR_2"; }
	195
	196	AliDebug(2,"--------------------------------------------");
	197	AliDebug(2,Form("volPath=%s\n",volPath.Data()));
	198	AliDebug(2,Form("symName=%s\n",symName.Data()));
	199	AliDebug(2,"--------------------------------------------");
	200	if(!gGeoManager->SetAlignableEntry(symName.Data(),volPath.Data()))
	201	AliFatal(Form("Alignable entry %s not created. Volume path %s not valid",
	202	symName.Data(),volPath.Data()));
	203	}
	204	}
	205	//------------------------------------------------------------------------
	206	void AliFITv0::CreateMaterials()
	207	{
	208
	209	printf("@@@@@@@@@@@@AliFITv0::CreateMaterials\n");
	210	Int_t isxfld = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Integ();
	211	Float_t sxmgmx = ((AliMagF*)TGeoGlobalMagField::Instance()->GetField())->Max();
	212	// Float_t a,z,d,radl,absl,buf[1];
	213	// Int_t nbuf;
	214	// AIR
	215
	216	Float_t aAir[4]={12.0107,14.0067,15.9994,39.948};
	217	Float_t zAir[4]={6.,7.,8.,18.};
	218	Float_t wAir[4]={0.000124,0.755267,0.231781,0.012827};
	219	Float_t dAir = 1.20479E-3;
	220	Float_t dAir1 = 1.20479E-11;
	221
	222	// PMT glass SiO2
	223	Float_t aglass[2]={28.0855,15.9994};
	224	Float_t zglass[2]={14.,8.};
	225	Float_t wglass[2]={1.,2.};
	226	Float_t dglass=2.65;
	227
	228
	229	//* Definition Of avaible T0 materials *
	230	AliMixture(1, "Vacuum$", aAir, zAir, dAir1,4,wAir);
	231	AliMixture(2, "Air$", aAir, zAir, dAir,4,wAir);
	232	AliMixture( 4, "PMT glass $",aglass,zglass,dglass,-2,wglass);
	233
	234	AliMedium(1, "FIT_Air$", 2, 0, isxfld, sxmgmx, 10., .1, 1., .003, .003);
	235	AliMedium(22, "FIT_AirSens$", 2, 1, isxfld, sxmgmx, 10., .1, 1., .003, .003);
	236	AliMedium(3, "FIT_Vacuum$", 1, 0, isxfld, sxmgmx, 10., .01, .1, .003, .003);
	237	AliMedium(6, "Glass$", 4, 1, isxfld, sxmgmx, 10., .01, .1, .003, .003);
	238
	239	AliDebugClass(1,": ++++++++++++++Medium set++++++++++");
	240
241
242	}
243
244
245	//-------------------------------------------------------------------
246	void AliFITv0::Init()
247	{
248	// Initialises version 0 of the Forward Multiplicity Detector
249	//
250	AliFIT::Init();
2942f542	251	fIdSens1=TVirtualMC::GetMC()->VolId("0REG");
c1c44db3	252
	253	AliDebug(1,Form("%s: * FIT version 0 initialized *\n",ClassName()));
	254	}
	255
	256	//-------------------------------------------------------------------
	257
	258	void AliFITv0::StepManager()
	259	{
	260	//
	261	// Called for every step in the T0 Detector
	262	//
	263	Int_t id,copy,copy1;
	264	static Float_t hits[6];
	265	static Int_t vol[3];
	266	TLorentzVector pos;
	267	TLorentzVector mom;
	268
	269	// TClonesArray &lhits = *fHits;
	270
2942f542	271	if(!TVirtualMC::GetMC()->IsTrackAlive()) return; // particle has disappeared
c1c44db3	272
2942f542	273	id=TVirtualMC::GetMC()->CurrentVolID(copy);
c1c44db3	274	// Check the sensetive volume
c1c44db3	275	if(id==fIdSens1 ) {
2942f542	276	if(TVirtualMC::GetMC()->IsTrackEntering()) {
2942f542	277	TVirtualMC::GetMC()->CurrentVolOffID(1,copy1);
c1c44db3	278	vol[1] = copy1;
c1c44db3	279	vol[0]=copy;
2942f542	280	TVirtualMC::GetMC()->TrackPosition(pos);
c1c44db3	281	hits[0] = pos[0];
	282	hits[1] = pos[1];
	283	hits[2] = pos[2];
	284	if(pos[2]<0) vol[2] = 0;
	285	else vol[2] = 1 ;
	286	printf(" volumes pmt %i mcp %i side %i x %f y %f z %f\n", vol[0], vol[1], vol[2], hits[0], hits[1], hits[2] );
	287
2942f542	288	Float_t etot=TVirtualMC::GetMC()->Etot();
c1c44db3	289	hits[3]=etot;
2942f542	290	Int_t iPart= TVirtualMC::GetMC()->TrackPid();
2942f542	291	Int_t partID=TVirtualMC::GetMC()->IdFromPDG(iPart);
c1c44db3	292	hits[4]=partID;
2942f542	293	Float_t ttime=TVirtualMC::GetMC()->TrackTime();
c1c44db3	294	hits[5]=ttime*1e12;
	295	AddHit(gAlice->GetMCApp()->GetCurrentTrackNumber(),vol,hits);
	296	// Create a track reference at the exit of photocatode
	297	}
	298
	299	//charge particle
2942f542	300	if ( TVirtualMC::GetMC()->TrackCharge() )
c1c44db3	301	AddTrackReference(gAlice->GetMCApp()->GetCurrentTrackNumber(), AliTrackReference::kFIT);
	302	} //sensitive
	303
	304	}
	305