* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-Revision 1.13 1999/10/06 10:15:19 fca
-Correct bug in allocation of layer name and add destructor
-
-Revision 1.12 1999/10/05 08:05:09 fca
-Minor corrections for uninitialised variables.
-
-Revision 1.11 1999/09/29 09:24:20 fca
-Introduction of the Copyright and cvs Log
-
-*/
+/* $Id$ */
///////////////////////////////////////////////////////////////////////////////
-// //
-// Inner Traking System version 5 //
-// This class contains the base procedures for the Inner Tracking System //
-// //
+//
+// Inner Traking System version 5
+// This class contains the base procedures for the Inner Tracking System
+//
// Authors: R. Barbera, B. S. Nilsen.
// version 5.
// Created September 17 1999.
-// //
+//
///////////////////////////////////////////////////////////////////////////////
-
+#include <Riostream.h>
#include <stdio.h>
#include <stdlib.h>
+
+#include <TBRIK.h>
+#include <TCanvas.h>
+#include <TClonesArray.h>
+#include <TFile.h> // only required for Tracking function?
+#include <TGeometry.h>
+#include <TLorentzVector.h>
#include <TMath.h>
+#include <TNode.h>
+#include <TObjArray.h>
+#include <TObjString.h>
+#include <TSystem.h>
+#include <TTUBE.h>
+#include <TVirtualMC.h>
#include "AliRun.h"
-#include "TSystem.h"
#include "AliITShit.h"
-#include "AliITS.h"
-#include "AliITSv5.h"
+#include "AliITSGeant3Geometry.h"
#include "AliITSgeom.h"
+#include "AliITSgeomSDD.h"
+#include "AliITSgeomSPD.h"
+#include "AliITSgeomSSD.h"
+#include "AliITShit.h"
+#include "AliITSv5.h"
+#include "AliRun.h"
+#include "AliMC.h"
ClassImp(AliITSv5)
//_____________________________________________________________________________
AliITSv5::AliITSv5() {
- //
- // Standard constructor for the ITS
- //
- fId5N = 6;
- fId5Name = new char*[fId5N];
- fId5Name[0] = "ITS1";
- fId5Name[1] = "ITS2";
- fId5Name[2] = "ITS3";
- fId5Name[3] = "ITS4";
- fId5Name[4] = "ITS5";
- fId5Name[5] = "ITS6";
-}
-//_____________________________________________________________________________
-AliITSv5::~AliITSv5() {
- //
- // Standard destructor for the ITS
- //
- delete [] fId5Name;
+////////////////////////////////////////////////////////////////////////
+// Standard default constructor for the ITS version 5.
+////////////////////////////////////////////////////////////////////////
+ Int_t i;
+
+ fIdN = 0;
+ fIdName = 0;
+ fIdSens = 0;
+ fEuclidOut = kFALSE; // Don't write Euclide file
+ fGeomDetOut = kFALSE; // Don't write .det file
+ fGeomDetIn = kFALSE; // Don't Read .det file
+ fGeomOldDetIn = kFALSE; // Don't Read old formatted .det file
+ fMajorVersion = IsVersion();
+ fMinorVersion = 1;
+ for(i=0;i<60;i++) fRead[i] = '\0';
+ for(i=0;i<60;i++) fWrite[i] = '\0';
+ for(i=0;i<60;i++) fEuclidGeomDet[i] = '\0';
}
//_____________________________________________________________________________
AliITSv5::AliITSv5(const char *name, const char *title) : AliITS(name, title){
- //
- // Standard constructor for the ITS
- //
- fId5N = 6;
- fId5Name = new char*[fId5N];
- fId5Name[0] = "ITS1";
- fId5Name[1] = "ITS2";
- fId5Name[2] = "ITS3";
- fId5Name[3] = "ITS4";
- fId5Name[4] = "ITS5";
- fId5Name[5] = "ITS6";
+////////////////////////////////////////////////////////////////////////
+// Standard constructor for the ITS version 5.
+////////////////////////////////////////////////////////////////////////
+ Int_t i;
+
+ fIdN = 6;
+ fIdName = new TString[fIdN];
+ fIdName[0] = "ITS1";
+ fIdName[1] = "ITS2";
+ fIdName[2] = "ITS3";
+ fIdName[3] = "ITS4";
+ fIdName[4] = "ITS5";
+ fIdName[5] = "ITS6";
+ fIdSens = new Int_t[fIdN];
+ for (i=0;i<fIdN;i++) fIdSens[i] = 0;
+ fEuclidOut = kFALSE; // Don't write Euclide file
+ fGeomDetOut = kFALSE; // Don't write .det file
+ fGeomDetIn = kFALSE; // Don't Read .det file
+ fGeomOldDetIn = kFALSE; // Don't Read old formatted .det file
+ fMajorVersion = IsVersion();
+ fMinorVersion = 1;
+ for(i=0;i<60;i++) fRead[i] = '\0';
+ for(i=0;i<60;i++) fWrite[i] = '\0';
fEuclidMaterial = "$ALICE_ROOT/Euclid/ITSgeometry_5.tme";
fEuclidGeometry = "$ALICE_ROOT/Euclid/ITSgeometry_5.euc";
+ strncpy(fEuclidGeomDet,"$ALICE_ROOT/ITS/ITSgeometry_v5.det",60);
+ strncpy(fRead,fEuclidGeomDet,60);
+ strncpy(fWrite,fEuclidGeomDet,60);
+}
+//____________________________________________________________________________
+AliITSv5::AliITSv5(const AliITSv5 &source) : AliITS(source){
+////////////////////////////////////////////////////////////////////////
+// Copy Constructor for ITS version 5.
+////////////////////////////////////////////////////////////////////////
+ if(&source == this) return;
+ Warning("Copy Constructor","Not allowed to copy AliITSv5");
+ return;
}
+//_____________________________________________________________________________
+AliITSv5& AliITSv5::operator=(const AliITSv5 &source){
+////////////////////////////////////////////////////////////////////////
+// Assignment operator for the ITS version 5.
+////////////////////////////////////////////////////////////////////////
-
+ if(&source == this) return *this;
+ Warning("= operator","Not allowed to copy AliITSv5");
+ return *this;
+
+}
//_____________________________________________________________________________
-void AliITSv5::CreateMaterials(){
+AliITSv5::~AliITSv5() {
+////////////////////////////////////////////////////////////////////////
+// Standard destructor for the ITS version 5.
+////////////////////////////////////////////////////////////////////////
+}
+//______________________________________________________________________
+void AliITSv5::BuildGeometry(){
+////////////////////////////////////////////////////////////////////////
+// Geometry builder for the ITS version 5.
+////////////////////////////////////////////////////////////////////////
//
- // Read materials for the ITS
+ // Build ITS TNODE geometry for event display using detailed geometry.
+ // This function builds a simple ITS geometry used by the ROOT macro
+ // ITSdisplay.C.
+
+ TNode *top;
+ TNode *nd;
+ //const int kColorITSSPD=kRed;
+ //const int kColorITSSDD=kGreen;
+ const int kColorITSSSD=kBlue;
//
+ AliITSgeom *gm = this->GetITSgeom();
+ if(gm==0) return;
+ top=gAlice->GetGeometry()->GetNode("alice");
+
+ Int_t lay,lad,det,i;
+ Text_t name[10];
+ Float_t xg[3];
+ Float_t rt[9];
+ Double_t rtd[9];
+ TBRIK *box;
+ TRotMatrix *rm;
+ //TCanvas *c1 = new TCanvas("c1","ITS");
+
+ for(lay=1;lay<=2;lay++)
+ for(lad=1;lad<=gm->GetNladders(lay);lad++)
+ for(det=1;det<=gm->GetNdetectors(lay);det++){
+ try {
+ box = new TBRIK ("ActiveSPD","Active volume of SPD","SPD SI DET",
+ 0.64,0.0075,4.19);
+ } catch (...) {
+ cout << "EXCEPTION in box = new TBRIK" << endl;
+ return;
+ }
+ gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
+ gm->GetRotMatrix(lay,lad,det,rt);
+ //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
+ for(i=0;i<9;i++) rtd[i] = rt[i];
+ try {
+ rm = new TRotMatrix(name,name,rtd);
+ } catch (...) {
+ cout << "EXCEPTION in new TRotMatrix" << endl;
+ return;
+ }
+ top->cd();
+ //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
+ try {
+ nd = new TNode("SPD"," ",box,xg[0],xg[1],xg[2],rm);
+ } catch (...) {
+ cout << "EXCEPTION in new TNode" << endl;
+ return;
+ }
+ nd->SetLineColor(kColorITSSSD);
+ fNodes->Add(nd);
+ }
+
+ for(lay=3;lay<=3;lay++)
+ for(lad=1;lad<=gm->GetNladders(lay);lad++)
+ for(det=1;det<=gm->GetNdetectors(lay);det++){
+ try {
+ box = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET",
+ 3.5,0.014,3.763);
+ } catch (...) {
+ cout << "EXCEPTION in box = new TBRIK" << endl;
+ return;
+ }
+ gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
+ gm->GetRotMatrix(lay,lad,det,rt);
+ //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
+ for(i=0;i<9;i++) rtd[i] = rt[i];
+ try {
+ rm = new TRotMatrix(name,name,rtd);
+ } catch (...) {
+ cout << "EXCEPTION in new TRotMatrix" << endl;
+ return;
+ }
+ top->cd();
+ //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
+ try {
+ nd = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm);
+ } catch (...) {
+ cout << "EXCEPTION in new TNode" << endl;
+ return;
+ }
+ nd->SetLineColor(kColorITSSSD);
+ fNodes->Add(nd);
+ }
+
+ for(lay=4;lay<=4;lay++)
+ for(lad=1;lad<=gm->GetNladders(lay);lad++)
+ for(det=1;det<=gm->GetNdetectors(lay);det++){
+ try {
+ box = new TBRIK ("ActiveSDD","Active volume of SDD","SDD SI DET",
+ 3.5,0.014,3.763);
+ } catch (...) {
+ cout << "EXCEPTION in box = new TBRIK" << endl;
+ return;
+ }
+ gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
+ gm->GetRotMatrix(lay,lad,det,rt);
+ //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
+ for(i=0;i<9;i++) rtd[i] = rt[i];
+ try {
+ rm = new TRotMatrix(name,name,rtd);
+ } catch (...) {
+ cout << "EXCEPTION in new TRotMatrix" << endl;
+ return;
+ }
+ top->cd();
+ //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
+ try {
+ nd = new TNode("SDD"," ",box,xg[0],xg[1],xg[2],rm);
+ } catch (...) {
+ cout << "EXCEPTION in new TNode" << endl;
+ return;
+ }
+ nd->SetLineColor(kColorITSSSD);
+ fNodes->Add(nd);
+ }
+ for(lay=5;lay<=5;lay++)
+ for(lad=1;lad<=gm->GetNladders(lay);lad++)
+ for(det=1;det<=gm->GetNdetectors(lay);det++){
+ try {
+ box = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET",
+ 3.65,0.015,2.0);
+ } catch (...) {
+ cout << "EXCEPTION in box = new TBRIK" << endl;
+ return;
+ }
+ gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
+ gm->GetRotMatrix(lay,lad,det,rt);
+ //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
+ for(i=0;i<9;i++) rtd[i] = rt[i];
+ try {
+ rm = new TRotMatrix(name,name,rtd);
+ } catch (...) {
+ cout << "EXCEPTION in new TRotMatrix" << endl;
+ return;
+ }
+ top->cd();
+ //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
+ try {
+ nd = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm);
+ } catch (...) {
+ cout << "EXCEPTION in new TNode" << endl;
+ return;
+ }
+ nd->SetLineColor(kColorITSSSD);
+ fNodes->Add(nd);
+ }
+
+ for(lay=6;lay<=6;lay++)
+ for(lad=1;lad<=gm->GetNladders(lay);lad++)
+ for(det=1;det<=gm->GetNdetectors(lay);det++){
+ try {
+ box = new TBRIK ("ActiveSSD","Active volume of SSD","SSD SI DET",
+ 3.65,0.015,2.0);
+ } catch (...) {
+ cout << "EXCEPTION in box = new TBRIK" << endl;
+ return;
+ }
+
+ gm->GetTrans(lay,lad,det,xg[0],xg[1],xg[2]);
+ gm->GetRotMatrix(lay,lad,det,rt);
+ //sprintf(name,"ROT%1.1d2.2d2.2d",lay,lad,det);
+ for(i=0;i<9;i++) rtd[i] = rt[i];
+ try {
+ rm = new TRotMatrix(name,name,rtd);
+ } catch (...) {
+ cout << "EXCEPTION in new TRotMatrix" << endl;
+ return;
+ }
+ top->cd();
+ //sprintf(name,"ND%1.1d2.2d2.2d",lay,lad,det);
+ try {
+ nd = new TNode("SSD"," ",box,xg[0],xg[1],xg[2],rm);
+ } catch (...) {
+ cout << "EXCEPTION in new TNode" << endl;
+ return;
+ }
+ nd->SetLineColor(kColorITSSSD);
+ fNodes->Add(nd);
+ }
+
+
+}
+//_____________________________________________________________________________
+void AliITSv5::CreateMaterials(){
+////////////////////////////////////////////////////////////////////////
+// Read a file containing the materials for the ITS version 5.
+////////////////////////////////////////////////////////////////////////
char *filtmp;
-//
+
filtmp = gSystem->ExpandPathName(fEuclidMaterial.Data());
-// FILE *file = fopen(fEuclidMaterial.Data(),"r");
+
FILE *file = fopen(filtmp,"r");
if(file) {
fclose(file);
-// gAlice->ReadEuclidMedia(fEuclidMaterial.Data(),this);
- gAlice->ReadEuclidMedia(filtmp,this);
+ ReadEuclidMedia(filtmp);
} else {
- Error("CreateMaterials"," THE MEDIA FILE %s DOES NOT EXIST !",
-// fEuclidMaterial.Data());
- filtmp);
+ Error("CreateMaterials"," THE MEDIA FILE %s DOES NOT EXIST !",filtmp);
exit(1);
} // end if(file)
}
-
//_____________________________________________________________________________
void AliITSv5::CreateGeometry(){
//////////////////////////////////////////////////////////////////////
<pre>
*/
//End_Html
-////////////////////////////////////////////////////////////////////////
-// Read geometry for the ITS
//
+// Read a file containing the geometry for the ITS version 5.
+////////////////////////////////////////////////////////////////////////
+
- Int_t size;
char topvol[5];
char *filtmp;
-//
+
filtmp = gSystem->ExpandPathName(fEuclidGeometry.Data());
FILE *file = fopen(filtmp,"r");
delete [] filtmp;
if(file) {
fclose(file);
- printf("Ready to read Euclid geometry file\n");
- gAlice->ReadEuclid(fEuclidGeometry.Data(),this,topvol);
+ cout << "Ready to read Euclid geometry file" << endl;
+ ReadEuclid(fEuclidGeometry.Data(),topvol);
printf("Read in euclid geometries\n");
} else {
Error("CreateGeometry"," THE GEOM FILE %s DOES NOT EXIST !",
exit(1);
} // end if(file)
//
- //---Place the ITS ghost volume ITSV in its mother volume (ALIC) and make it
- // invisible
+ // Place the ITS ghost volume ITSV in its mother volume (ALIC) and make it
+ // invisible
//
gMC->Gspos("ITSV",1,"ALIC",0,0,0,0,"ONLY");
//
- //---Outputs the geometry tree in the EUCLID/CAD format
+ // Outputs the geometry tree in the EUCLID/CAD format if requested to do so
if (fEuclidOut) {
gMC->WriteEuclid("ITSgeometry", "ITSV", 1, 5);
} // end if (fEuclidOut)
- filtmp = gSystem->ExpandPathName(fEuclidGeometry.Data());
+ cout << "finished with euclid geometrys" << endl;
+}
+//______________________________________________________________________
+void AliITSv5::ReadOldGeometry(const char *filename){
+ // read in the file containing the transformations for the active
+ // volumes for the ITS version 5. This is expected to be in a file
+ // ending in .det. This geometry is kept in the AliITSgeom class.
+ Int_t size;
+ char *filtmp;
+ FILE *file;
+
+ if(fITSgeom!=0) delete fITSgeom;
+ filtmp = gSystem->ExpandPathName(filename);
size = strlen(filtmp);
- if(size>4){
- filtmp[size-3] = 'd'; // change from .euc to .det
+ if(size>4 && fGeomDetIn){
+ filtmp[size-3] = 'd'; // change from .euc to .det
filtmp[size-2] = 'e';
filtmp[size-1] = 't';
- file = fopen(filtmp,"r");
- if(file){ // if file exists use it to fill AliITSgeom structure.
- fclose(file);
- printf("ready to read .det file %s\n",filtmp);
- fITSgeom = new AliITSgeom(filtmp);
- }else{
- fITSgeom = 0;
- // fill AliITSgeom structure from geant structure just filled above
- }// end if(file)
+ file = fopen(filtmp,"r");
+ if(file){ // if file exists use it to fill AliITSgeom structure.
+ fclose(file);
+ fITSgeom = new AliITSgeom(filtmp);
+ fITSgeom->DefineShapes(3); // if fShape isn't defined define it.
+ // Now define the detector types/shapes.
+ fITSgeom->ReSetShape(kSPD,new AliITSgeomSPD300());
+ fITSgeom->ReSetShape(kSDD,new AliITSgeomSDD300());
+ fITSgeom->ReSetShape(kSSD,new AliITSgeomSSD175());
+ }else{
+ fITSgeom = 0;
+ // fill AliITSgeom structure from geant structure just filled above
+ }// end if(file)
delete [] filtmp;
}// end if(size>4)
- printf("finished with euclid geometrys\n");
}
+//______________________________________________________________________
+void AliITSv5::InitAliITSgeom(){
+// Based on the geometry tree defined in Geant 3.21, this
+// routine initilizes the Class AliITSgeom from the Geant 3.21 ITS geometry
+// sturture.
+// if(gMC->IsA()!=TGeant3::Class()) {
+ if(strcmp(gMC->GetName(),"TGeant3")) {
+ Error("InitAliITSgeom",
+ "Wrong Monte Carlo. InitAliITSgeom uses TGeant3 calls");
+ return;
+ } // end if
+ cout << "Reading Geometry transformation directly from Geant 3." << endl;
+ const Int_t nlayers = 6;
+ const Int_t ndeep = 7;
+ Int_t itsGeomTreeNames[nlayers][ndeep],lnam[20],lnum[20];
+ Int_t nlad[nlayers],ndet[nlayers];
+ Double_t t[3],r[10];
+ Float_t par[20],att[20];
+ Int_t npar,natt,idshape,imat,imed;
+ AliITSGeant3Geometry *ig = new AliITSGeant3Geometry();
+ Int_t mod,lay,lad,det,i,j,k;
+ char *names[nlayers][ndeep] = {
+ {"ALIC","ITSV","ITSD","IT12","I132","I186","ITS1"}, // lay=1
+ {"ALIC","ITSV","ITSD","IT12","I132","I131","ITS2"}, // lay=2
+ {"ALIC","ITSV","ITSD","IT34","I004","I302","ITS3"}, // lay=3
+ {"ALIC","ITSV","ITSD","IT34","I005","I402","ITS4"}, // lay=4
+ {"ALIC","ITSV","ITSD","IT56","I565","I562","ITS5"}, // lay=5
+ {"ALIC","ITSV","ITSD","IT56","I569","I566","ITS6"}};// lay=6
+ Int_t itsGeomTreeCopys[nlayers][ndeep] = {{1,1,1,1,10, 2,4}, // lay=1
+ {1,1,1,1,10, 4,4}, // lay=2
+ {1,1,1,1,14, 6,1}, // lay=3
+ {1,1,1,1,22, 8,1}, // lay=4
+ {1,1,1,1,34,23,1}, // lay=5
+ {1,1,1,1,38,26,1}};// lay=6
+
+ // Sorry, but this is not very pritty code. It should be replaced
+ // at some point with a version that can search through the geometry
+ // tree its self.
+ for(i=0;i<20;i++) lnam[i] = lnum[i] = 0;
+ for(i=0;i<nlayers;i++)for(j=0;j<ndeep;j++)
+ strncpy((char*) &itsGeomTreeNames[i][j],names[i][j],4);
+ // itsGeomTreeNames[i][j] = ig->StringToInt(names[i][j]);
+ mod = 0;
+ for(i=0;i<nlayers;i++){
+ k = 1;
+ for(j=0;j<ndeep;j++) if(itsGeomTreeCopys[i][j]!=0)
+ k *= TMath::Abs(itsGeomTreeCopys[i][j]);
+ mod += k;
+ } // end for i
+ if(fITSgeom!=0) delete fITSgeom;
+ nlad[0]=20;nlad[1]=40;nlad[2]=14;nlad[3]=22;nlad[4]=34;nlad[5]=38;
+ ndet[0]=4;ndet[1]=4;ndet[2]=6;ndet[3]=8;ndet[4]=22;ndet[5]=25;
+ fITSgeom = new AliITSgeom(0,6,nlad,ndet,mod);
+ mod = -1;
+ for(lay=1;lay<=nlayers;lay++){
+ for(j=0;j<ndeep;j++) lnam[j] = itsGeomTreeNames[lay-1][j];
+ for(j=0;j<ndeep;j++) lnum[j] = itsGeomTreeCopys[lay-1][j];
+ switch (lay){
+ case 1: case 2: // layers 1 and 2 are a bit special
+ lad = 0;
+ for(j=1;j<=itsGeomTreeCopys[lay-1][4];j++){
+ lnum[4] = j;
+ for(k=1;k<=itsGeomTreeCopys[lay-1][5];k++){
+ lad++;
+ lnum[5] = k;
+ for(det=1;det<=itsGeomTreeCopys[lay-1][6];det++){
+ lnum[6] = det;
+ mod++;
+ ig->GetGeometry(ndeep,lnam,lnum,t,r,idshape,npar,natt,
+ par,att,imat,imed);
+ fITSgeom->CreatMatrix(mod,lay,lad,det,kSPD,t,r);
+ if(!(fITSgeom->IsShapeDefined((Int_t)kSPD)))
+ fITSgeom->ReSetShape(kSPD,
+ new AliITSgeomSPD300());
+ } // end for det
+ } // end for k
+ } // end for j
+ break;
+ case 3: case 4: case 5: case 6: // layers 3-6
+ lnum[6] = 1;
+ for(lad=1;lad<=itsGeomTreeCopys[lay-1][4];lad++){
+ lnum[4] = lad;
+ for(det=1;det<=itsGeomTreeCopys[lay-1][5];det++){
+ lnum[5] = det;
+ mod++;
+ ig->GetGeometry(7,lnam,lnum,t,r,idshape,npar,natt,
+ par,att,imat,imed);
+ switch (lay){
+ case 3: case 4:
+ fITSgeom->CreatMatrix(mod,lay,lad,det,kSDD,t,r);
+ if(!(fITSgeom->IsShapeDefined(kSDD)))
+ fITSgeom->ReSetShape(kSDD,new AliITSgeomSDD300());
+ break;
+ case 5: case 6:
+ fITSgeom->CreatMatrix(mod,lay,lad,det,kSSD,t,r);
+ if(!(fITSgeom->IsShapeDefined(kSSD)))
+ fITSgeom->ReSetShape(kSSD,new AliITSgeomSSD175());
+ break;
+ } // end switch
+ } // end for det
+ } // end for lad
+ break;
+ } // end switch
+ } // end for lay
+ return;
+}
//_____________________________________________________________________________
void AliITSv5::Init(){
- //
- // Initialise the ITS after it has been created
- //
- Int_t i,j,l;
-
- fIdN = fId5N;
- fIdName = new char*[fId5N];
- fIdSens = new Int_t[fId5N];
- for(i=0;i<fId5N;i++) {
- l = strlen(fId5Name[i]);
- fIdName[i] = new char[l+1];
- for(j=0;j<l;j++) fIdName[i][j] = fId5Name[i][j];
- fIdName[i][l] = '\0'; // Null terminate this string.
- } // end for i
+////////////////////////////////////////////////////////////////////////
+// Initialise the ITS after it has been created.
+////////////////////////////////////////////////////////////////////////
+ Int_t i;
+ cout << endl;
+ for(i=0;i<30;i++) cout << "*";cout << " ITSv5_Init ";
+ for(i=0;i<30;i++) cout << "*";cout << endl;
+//
+ if(fRead[0]=='\0') strncpy(fRead,fEuclidGeomDet,60);
+ if(fWrite[0]=='\0') strncpy(fWrite,fEuclidGeomDet,60);
+ if(fGeomDetIn && !fGeomOldDetIn){
+ if(fITSgeom!=0) delete fITSgeom;
+ fITSgeom = new AliITSgeom();
+ fITSgeom->ReadNewFile(fRead);
+ } // end if
+ if(fGeomDetIn && fGeomOldDetIn) ReadOldGeometry(fEuclidGeometry.Data());
+
+ if(!fGeomDetIn) this->InitAliITSgeom();
+ if(fGeomDetOut) fITSgeom->WriteNewFile(fWrite);
AliITS::Init();
-}
-
+//
+ for(i=0;i<72;i++) cout << "*";
+ cout << endl;
+}
//_____________________________________________________________________________
void AliITSv5::StepManager(){
- //
- // Called for every step in the ITS
- //
+////////////////////////////////////////////////////////////////////////
+// Called for every step in the ITS, then calles the AliITShit class
+// creator with the information to be recoreded about that hit.
+// The value of the macro ALIITSPRINTGEOM if set to 1 will allow the
+// printing of information to a file which can be used to create a .det
+// file read in by the routine CreateGeometry(). If set to 0 or any other
+// value except 1, the default behavior, then no such file is created nor
+// it the extra variables and the like used in the printing allocated.
+////////////////////////////////////////////////////////////////////////
Int_t copy, id;
Int_t copy1,copy2;
Float_t hits[8];
hits[6]=gMC->Edep();
hits[7]=gMC->TrackTime();
// Fill hit structure with this new hit.
- new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->CurrentTrack(),vol,hits);
+ new(lhits[fNhits++]) AliITShit(fIshunt,gAlice->GetMCApp()->GetCurrentTrackNumber(),vol,hits);
return;
}
-//____________________________________________________________________________
-void AliITSv5::Streamer(TBuffer &R__b)
-{
- // Stream an object of class AliITSv5.
-
- if (R__b.IsReading()) {
- Version_t R__v = R__b.ReadVersion(); if (R__v) { }
- AliITS::Streamer(R__b);
- // This information does not need to be read. It is "hard wired"
- // into this class via its creators.
- //R__b >> fId5N;
- //R__b.ReadArray(fId5Name);
- } else {
- R__b.WriteVersion(AliITSv5::IsA());
- AliITS::Streamer(R__b);
- // This information does not need to be saved. It is "hard wired"
- // into this class via its creators.
- //R__b << fId5N;
- //R__b.WriteArray(fId5Name, __COUNTER__);
- }
-}
-
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff