From 36549c881e1bcdb14c1bdb54b3af55b897dc6040 Mon Sep 17 00:00:00 2001
From: barbera <barbera@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Fri, 2 Mar 2001 19:36:40 +0000
Subject: [PATCH] new class to take into account ITS material distribution in 
 tracking v1

---
 ITS/AliITSRad.cxx | 127 ++++++++++++++++++++++++++++++++++++++++++++++
 ITS/AliITSRad.h   |  50 ++++++++++++++++++
 2 files changed, 177 insertions(+)
 create mode 100644 ITS/AliITSRad.cxx
 create mode 100644 ITS/AliITSRad.h

diff --git a/ITS/AliITSRad.cxx b/ITS/AliITSRad.cxx
new file mode 100644
index 00000000000..54b0d4c892c
--- /dev/null
+++ b/ITS/AliITSRad.cxx
@@ -0,0 +1,127 @@
+#include <TMatrix.h>
+#include <TObjArray.h>
+#include <iostream.h>
+#include <fstream.h>
+
+
+#include "AliITSRad.h"
+
+
+ClassImp(AliITSRad)
+
+AliITSRad::AliITSRad(Int_t iimax, Int_t jjmax) {
+
+  imax=iimax;
+  jmax=jjmax;
+  
+  fmrad1 = new TMatrix(imax,jmax);
+  fmrad2 = new TMatrix(imax,jmax);
+  fmrad3 = new TMatrix(imax,jmax);
+  fmrad4 = new TMatrix(imax,jmax);
+  fmrad5 = new TMatrix(imax,jmax);
+  fmrad6 = new TMatrix(imax,jmax);
+
+  ifstream in("ITSlegov5.map");
+  Int_t i,j;
+  
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   in>>(*fmrad1)(i,j);
+	 }
+  }
+	 
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   in>>(*fmrad2)(i,j);
+	 }
+  } 
+	 
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   in>>(*fmrad3)(i,j);
+	 }
+  } 
+	 
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   in>>(*fmrad4)(i,j);
+	 }
+  } 
+	 
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   in>>(*fmrad5)(i,j);
+	 }
+  }
+  
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   in>>(*fmrad6)(i,j);
+	 }
+  }
+	 	 
+  in.close();
+
+/*
+/////////////////////////////////////////////////////////////////////////////////////////////////////////// 
+//                       Stampe provvisorie delle matrici su rad.out         
+//	 
+  ofstream out("rad.out");
+  	 
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   out<<(*fmrad1)(i,j)<<" ";
+	 }
+	 out<<"\n";
+  }
+  out<<"\n";
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   out<<(*fmrad2)(i,j)<<" ";
+	 }
+	 out<<"\n";
+  }
+  out<<"\n";
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   out<<(*fmrad3)(i,j)<<" ";
+	 }
+	 out<<"\n";
+  }
+  out<<"\n";
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   out<<(*fmrad4)(i,j)<<" ";
+	 }
+	 out<<"\n";
+  }
+  out<<"\n"; 
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   out<<(*fmrad5)(i,j)<<" ";
+	 }
+	 out<<"\n";
+  }
+  out<<"\n";
+  for(i=0; i<imax; i++) {
+    for(j=0; j<jmax; j++) {
+	   out<<(*fmrad6)(i,j)<<" ";
+	 }
+	 out<<"\n";
+  }
+  
+  out.close(); 
+///////////////////////////////////////////////////////////////////////////////////////////////////////////	 
+*/
+}
+
+AliITSRad::~AliITSRad() {
+
+  delete fmrad1;
+  delete fmrad2;
+  delete fmrad3;
+  delete fmrad4;
+  delete fmrad5;
+  delete fmrad6;
+  
+}
diff --git a/ITS/AliITSRad.h b/ITS/AliITSRad.h
new file mode 100644
index 00000000000..0bdeb1e1952
--- /dev/null
+++ b/ITS/AliITSRad.h
@@ -0,0 +1,50 @@
+#ifndef ALIITSRAD_H
+#define ALIITSRAD_H
+
+#include <TObject.h>
+#include <TMatrix.h>
+
+
+class TObjArray;
+class TMatrix;
+
+
+//                  ITS Class to calculate the radiation lenght matrix
+//Origin  A. Badala' and G.S. Pappalardo:  e-mail Angela.Badala@ct.infn.it, Giuseppe.S.Pappalardo@ct.infn.it
+//
+
+
+class AliITSRad : public TObject { 
+
+public:
+
+  AliITSRad(Int_t iimax, Int_t jjmax);         // class constructor
+  ~AliITSRad();                                // class destructor
+  
+  Int_t Getimax() {return imax;}               // return the first dimension of the matrices
+  Int_t Getjmax() {return jmax;}               // return the second dimension of the matrices
+  
+  TMatrix &GetRadMatrix1() {return *fmrad1;}   // return the radiation lengh matrix for layer 1
+  TMatrix &GetRadMatrix2() {return *fmrad2;}   // return the radiation lengh matrix for layer 2
+  TMatrix &GetRadMatrix3() {return *fmrad3;}   // return the radiation lengh matrix for layer 3
+  TMatrix &GetRadMatrix4() {return *fmrad4;}   // return the radiation lengh matrix for layer 4
+  TMatrix &GetRadMatrix5() {return *fmrad5;}   // return the radiation lengh matrix for layer 5
+  TMatrix &GetRadMatrix6() {return *fmrad6;}   // return the radiation lengh matrix for layer 6
+  
+private:
+
+  Int_t           imax;        // first dimension of the matrices
+  Int_t           jmax;        // second dimension of the matrices
+  
+  TMatrix         *fmrad1;     // matrix of the radiation lenghts for layer 1
+  TMatrix         *fmrad2;     // matrix of the radiation lenghts for layer 2
+  TMatrix         *fmrad3;     // matrix of the radiation lenghts for layer 3
+  TMatrix         *fmrad4;     // matrix of the radiation lenghts for layer 4
+  TMatrix         *fmrad5;     // matrix of the radiation lenghts for layer 5
+  TMatrix         *fmrad6;     // matrix of the radiation lenghts for layer 6
+
+  ClassDef(AliITSRad, 1)
+};
+
+#endif
+
-- 
2.39.3