[u/mrichter/AliRoot.git] / ITS / AliITSv5.cxx

///////////////////////////////////////////////////////////////////////////////
//                                                                           //
//  Inner Traking System version 5                                           //
//  This class contains the base procedures for the Inner Tracking System    //
//                                                                           //
//Begin_Html
/*
<img src="gif/AliITSv5Class.gif">
</pre>
<br clear=left>
<font size=+2 color=red>
<p>The responsible person for this module is
<a href="mailto:roberto.barbera@ct.infn.it">Roberto Barbera</a>.
</font>
<pre>
*/
//End_Html
//                                                                           //
///////////////////////////////////////////////////////////////////////////////

#include <stdio.h> 

#include <TMath.h>
#include "AliITSv5.h"
#include "AliRun.h"
#include "stdlib.h"
#include "TSystem.h"
#include "TGeant3.h"

ClassImp(AliITSv5)
 
//_____________________________________________________________________________
AliITSv5::AliITSv5() : AliITS() 
{
  //
  // Default constructor for the ITS
  //
}
 
//_____________________________________________________________________________
AliITSv5::AliITSv5(const char *name, const char *title)
  : AliITS(name, title)
{
  //
  // Standard constructor for the ITS
  //
  fEuclidMaterial="$(ALICE_ROOT)/ITSgeometry_60.tme";
  fEuclidGeometry="$(ALICE_ROOT)/ITSgeometry_60.euc";
}

 
//_____________________________________________________________________________
void AliITSv5::CreateMaterials()
{
  //
  // Read materials for the ITS
  //
  FILE *file = fopen(fEuclidMaterial.Data(),"r");
  if(file) {
    fclose(file);
    gAlice->ReadEuclidMedia(fEuclidMaterial.Data(),2);
  } else {
    Error("CreateMaterials"," THE MEDIA FILE %s DOES NOT EXIST !",fEuclidMaterial.Data());
    exit(1);
  }
}

//_____________________________________________________________________________
void AliITSv5::CreateGeometry()
{
  //
  // Read geometry for the ITS
  //

  AliMC* pMC = AliMC::GetMC();
  
  char topvol[5];
  char *filtmp;
//
  filtmp=gSystem->ExpandPathName(fEuclidGeometry.Data());
  FILE *file = fopen(filtmp,"r");
  delete [] filtmp;
  if(file) {
    fclose(file);
    gAlice->ReadEuclid(fEuclidGeometry.Data(),2,topvol);
  } else {
    Error("CreateGeometry"," THE GEOM FILE %s DOES NOT EXIST !",fEuclidGeometry.Data());
    exit(1);
  }
  //
  // --- Place the ITS ghost volume ITSV in its mother volume (ALIC) and make it
  //     invisible
  //
  pMC->Gspos("ITSV",1,"ALIC",0,0,0,0,"ONLY");
  //
  // --- Outputs the geometry tree in the EUCLID/CAD format
  
    if (fEuclidOut) {
      pMC->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
    }
}

//_____________________________________________________________________________
void AliITSv5::Init()
{
  //
  // Initialise the ITS after it has been created
  //
  AliITS::Init();
} 
 
//_____________________________________________________________________________
void AliITSv5::StepManager()
{
  //
  // Called for every step in the ITS
  //
  Int_t         copy, id;
  Int_t         copy1,copy2;
  Float_t       hits[7];
  Int_t         vol[3];
  Float_t       position[3];
  Float_t       momentum[4];
  TClonesArray &lhits = *fHits;
  AliMC* pMC = AliMC::GetMC();
  //
  if(pMC->TrackCharge() && pMC->Edep()) {
    //
    // Only entering charged tracks
    if((id=pMC->CurrentVol(0,copy))==fIdSens1) {  
      vol[0]=1;
      id=pMC->CurrentVolOff(0,0,copy);        //detector copy in the ladder = 1<->4  (ITS1)
      vol[1]=copy;
      pMC->CurrentVolOff(1,0,copy1);      //ladder copy in the module   = 1<->2  (I186)
      pMC->CurrentVolOff(2,0,copy2);      //module copy in the layer    = 1<->10 (I132)
      vol[2]=copy1+(copy2-1)*2;               //# of ladders in one module  = 2          
    } else if(id==fIdSens2) {
      vol[0]=2;
      id=pMC->CurrentVolOff(0,0,copy);        //detector copy in the ladder = 1<->4  (ITS2)        
      vol[1]=copy;
      pMC->CurrentVolOff(1,0,copy1);      //ladder copy in the module   = 1<->4  (I131)
      pMC->CurrentVolOff(2,0,copy2);      //module copy in the layer    = 1<->10 (I132)
      vol[2]=copy1+(copy2-1)*4;   	          //# of ladders in one module  = 4   
    } else if(id==fIdSens3) {
      vol[0]=3;
      id=pMC->CurrentVolOff(1,0,copy);        //detector copy in the ladder = 1<->5  (ITS3 is inside I314)
      vol[1]=copy;
      id=pMC->CurrentVolOff(2,0,copy);        //ladder copy in the layer    = 1<->12 (I316)
      vol[2]=copy;             
    } else if(id==fIdSens4) {
      vol[0]=4;
      id=pMC->CurrentVolOff(1,0,copy);	      //detector copy in the ladder = 1<->8  (ITS4 is inside I414)
      vol[1]=copy;
      id=pMC->CurrentVolOff(2,0,copy);        //ladder copy in the layer    = 1<->22 (I417)
      vol[2]=copy;               
    } else if(id==fIdSens5) {
      vol[0]=5;
      id=pMC->CurrentVolOff(1,0,copy);	      //detector copy in the ladder = 1<->23  (ITS5 is inside I562)	  
      vol[1]=copy;
      id=pMC->CurrentVolOff(2,0,copy);        //ladder copy in the layer    = 1<->34 (I565)
      vol[2]=copy;               
    } else if(id==fIdSens6) {
      vol[0]=6;
      id=pMC->CurrentVolOff(1,0,copy);	      //detector copy in the ladder = 1<->26  (ITS6 is inside I566)	  
      vol[1]=copy;
      id=pMC->CurrentVolOff(2,0,copy);        //ladder copy in the layer    = 1<->38 (I569)
      vol[2]=copy;                      
    } else return;
    pMC->TrackPosition(position);
    pMC->TrackMomentum(momentum);
    hits[0]=position[0];
    hits[1]=position[1];
    hits[2]=position[2];          
    hits[3]=momentum[0]*momentum[3];
    hits[4]=momentum[1]*momentum[3];
    hits[5]=momentum[2]*momentum[3];        
    hits[6]=pMC->Edep();
    new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
  }      
}
Commit	Line	Data
fe4da5cc	1	///////////////////////////////////////////////////////////////////////////////
	2	// //
	3	// Inner Traking System version 5 //
	4	// This class contains the base procedures for the Inner Tracking System //
	5	// //
	6	//Begin_Html
	7	/*
	8	<img src="gif/AliITSv5Class.gif">
	9	</pre>
	10	<br clear=left>
	11	<font size=+2 color=red>
	12	<p>The responsible person for this module is
	13	<a href="mailto:roberto.barbera@ct.infn.it">Roberto Barbera</a>.
	14	</font>
	15	<pre>
	16	*/
	17	//End_Html
	18	// //
	19	///////////////////////////////////////////////////////////////////////////////
	20
	21	#include <stdio.h>
	22
	23	#include <TMath.h>
	24	#include "AliITSv5.h"
	25	#include "AliRun.h"
	26	#include "stdlib.h"
	27	#include "TSystem.h"
	28	#include "TGeant3.h"
	29
	30	ClassImp(AliITSv5)
	31
	32	//_____________________________________________________________________________
	33	AliITSv5::AliITSv5() : AliITS()
	34	{
	35	//
	36	// Default constructor for the ITS
	37	//
	38	}
	39
	40	//_____________________________________________________________________________
	41	AliITSv5::AliITSv5(const char name, const char title)
	42	: AliITS(name, title)
	43	{
	44	//
	45	// Standard constructor for the ITS
	46	//
	47	fEuclidMaterial="$(ALICE_ROOT)/ITSgeometry_60.tme";
	48	fEuclidGeometry="$(ALICE_ROOT)/ITSgeometry_60.euc";
	49	}
	50
	51
	52	//_____________________________________________________________________________
	53	void AliITSv5::CreateMaterials()
	54	{
	55	//
	56	// Read materials for the ITS
	57	//
	58	FILE *file = fopen(fEuclidMaterial.Data(),"r");
	59	if(file) {
	60	fclose(file);
	61	gAlice->ReadEuclidMedia(fEuclidMaterial.Data(),2);
	62	} else {
	63	Error("CreateMaterials"," THE MEDIA FILE %s DOES NOT EXIST !",fEuclidMaterial.Data());
	64	exit(1);
65	}
66	}
67
68	//_____________________________________________________________________________
69	void AliITSv5::CreateGeometry()
70	{
71	//
72	// Read geometry for the ITS
73	//
74
75	AliMC* pMC = AliMC::GetMC();
76
77	char topvol[5];
78	char *filtmp;
79	//
80	filtmp=gSystem->ExpandPathName(fEuclidGeometry.Data());
81	FILE *file = fopen(filtmp,"r");
82	delete [] filtmp;
83	if(file) {
84	fclose(file);
85	gAlice->ReadEuclid(fEuclidGeometry.Data(),2,topvol);
86	} else {
87	Error("CreateGeometry"," THE GEOM FILE %s DOES NOT EXIST !",fEuclidGeometry.Data());
88	exit(1);
89	}
90	//
91	// --- Place the ITS ghost volume ITSV in its mother volume (ALIC) and make it
92	// invisible
93	//
94	pMC->Gspos("ITSV",1,"ALIC",0,0,0,0,"ONLY");
95	//
96	// --- Outputs the geometry tree in the EUCLID/CAD format
97
98	if (fEuclidOut) {
99	pMC->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
100	}
101	}
102
103	//_____________________________________________________________________________
104	void AliITSv5::Init()
105	{
106	//
107	// Initialise the ITS after it has been created
108	//
109	AliITS::Init();
110	}
111
112	//_____________________________________________________________________________
113	void AliITSv5::StepManager()
114	{
115	//
116	// Called for every step in the ITS
117	//
118	Int_t copy, id;
119	Int_t copy1,copy2;
120	Float_t hits[7];
121	Int_t vol[3];
122	Float_t position[3];
123	Float_t momentum[4];
124	TClonesArray &lhits = *fHits;
125	AliMC* pMC = AliMC::GetMC();
126	//
127	if(pMC->TrackCharge() && pMC->Edep()) {
128	//
129	// Only entering charged tracks
130	if((id=pMC->CurrentVol(0,copy))==fIdSens1) {
131	vol[0]=1;
132	id=pMC->CurrentVolOff(0,0,copy); //detector copy in the ladder = 1<->4 (ITS1)
133	vol[1]=copy;
134	pMC->CurrentVolOff(1,0,copy1); //ladder copy in the module = 1<->2 (I186)
135	pMC->CurrentVolOff(2,0,copy2); //module copy in the layer = 1<->10 (I132)
136	vol[2]=copy1+(copy2-1)*2; //# of ladders in one module = 2
137	} else if(id==fIdSens2) {
138	vol[0]=2;
139	id=pMC->CurrentVolOff(0,0,copy); //detector copy in the ladder = 1<->4 (ITS2)
140	vol[1]=copy;
141	pMC->CurrentVolOff(1,0,copy1); //ladder copy in the module = 1<->4 (I131)
142	pMC->CurrentVolOff(2,0,copy2); //module copy in the layer = 1<->10 (I132)
143	vol[2]=copy1+(copy2-1)*4; //# of ladders in one module = 4
144	} else if(id==fIdSens3) {
145	vol[0]=3;
146	id=pMC->CurrentVolOff(1,0,copy); //detector copy in the ladder = 1<->5 (ITS3 is inside I314)
147	vol[1]=copy;
148	id=pMC->CurrentVolOff(2,0,copy); //ladder copy in the layer = 1<->12 (I316)
149	vol[2]=copy;
150	} else if(id==fIdSens4) {
151	vol[0]=4;
152	id=pMC->CurrentVolOff(1,0,copy); //detector copy in the ladder = 1<->8 (ITS4 is inside I414)
153	vol[1]=copy;
154	id=pMC->CurrentVolOff(2,0,copy); //ladder copy in the layer = 1<->22 (I417)
155	vol[2]=copy;
156	} else if(id==fIdSens5) {
157	vol[0]=5;
158	id=pMC->CurrentVolOff(1,0,copy); //detector copy in the ladder = 1<->23 (ITS5 is inside I562)
159	vol[1]=copy;
160	id=pMC->CurrentVolOff(2,0,copy); //ladder copy in the layer = 1<->34 (I565)
161	vol[2]=copy;
162	} else if(id==fIdSens6) {
163	vol[0]=6;
164	id=pMC->CurrentVolOff(1,0,copy); //detector copy in the ladder = 1<->26 (ITS6 is inside I566)
165	vol[1]=copy;
166	id=pMC->CurrentVolOff(2,0,copy); //ladder copy in the layer = 1<->38 (I569)
167	vol[2]=copy;
168	} else return;
169	pMC->TrackPosition(position);
170	pMC->TrackMomentum(momentum);
171	hits[0]=position[0];
172	hits[1]=position[1];
173	hits[2]=position[2];
174	hits[3]=momentum[0]*momentum[3];
175	hits[4]=momentum[1]*momentum[3];
176	hits[5]=momentum[2]*momentum[3];
177	hits[6]=pMC->Edep();
178	new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
179	}
180	}