]> git.uio.no Git - u/mrichter/AliRoot.git/blob - ITS/AliITSRad.cxx
Pyhtia comparison extended
[u/mrichter/AliRoot.git] / ITS / AliITSRad.cxx
1 #include <TMatrix.h>
2 #include <TObjArray.h>
3 #include <Riostream.h>
4 #include <Riostream.h>
5
6
7 #include "AliITSRad.h"
8
9
10 ClassImp(AliITSRad)
11 //______________________________________________________________________
12 AliITSRad::AliITSRad() {
13     // Default constructor
14
15
16   fmrad1 = 0;
17   fmrad2 = 0;
18   fmrad3 = 0;
19   fmrad4 = 0;
20   fmrad5 = 0;
21   fmrad6 = 0;
22 }
23 AliITSRad::AliITSRad(Int_t iimax, Int_t jjmax) {
24
25   imax=iimax;
26   jmax=jjmax;
27   
28   fmrad1 = new TMatrix(imax,jmax);
29   fmrad2 = new TMatrix(imax,jmax);
30   fmrad3 = new TMatrix(imax,jmax);
31   fmrad4 = new TMatrix(imax,jmax);
32   fmrad5 = new TMatrix(imax,jmax);
33   fmrad6 = new TMatrix(imax,jmax);
34
35   ifstream in("ITSlegov5.map");
36   Int_t i,j;
37   
38   for(i=0; i<imax; i++) {
39     for(j=0; j<jmax; j++) {
40            in>>(*fmrad1)(i,j);
41          }
42   }
43          
44   for(i=0; i<imax; i++) {
45     for(j=0; j<jmax; j++) {
46            in>>(*fmrad2)(i,j);
47          }
48   } 
49          
50   for(i=0; i<imax; i++) {
51     for(j=0; j<jmax; j++) {
52            in>>(*fmrad3)(i,j);
53          }
54   } 
55          
56   for(i=0; i<imax; i++) {
57     for(j=0; j<jmax; j++) {
58            in>>(*fmrad4)(i,j);
59          }
60   } 
61          
62   for(i=0; i<imax; i++) {
63     for(j=0; j<jmax; j++) {
64            in>>(*fmrad5)(i,j);
65          }
66   }
67   
68   for(i=0; i<imax; i++) {
69     for(j=0; j<jmax; j++) {
70            in>>(*fmrad6)(i,j);
71          }
72   }
73                  
74   in.close();
75
76 /*
77 /////////////////////////////////////////////////////////////////////////////////////////////////////////// 
78 //                       Stampe provvisorie delle matrici su rad.out         
79 //       
80   ofstream out("rad.out");
81          
82   for(i=0; i<imax; i++) {
83     for(j=0; j<jmax; j++) {
84            out<<(*fmrad1)(i,j)<<" ";
85          }
86          out<<"\n";
87   }
88   out<<"\n";
89   for(i=0; i<imax; i++) {
90     for(j=0; j<jmax; j++) {
91            out<<(*fmrad2)(i,j)<<" ";
92          }
93          out<<"\n";
94   }
95   out<<"\n";
96   for(i=0; i<imax; i++) {
97     for(j=0; j<jmax; j++) {
98            out<<(*fmrad3)(i,j)<<" ";
99          }
100          out<<"\n";
101   }
102   out<<"\n";
103   for(i=0; i<imax; i++) {
104     for(j=0; j<jmax; j++) {
105            out<<(*fmrad4)(i,j)<<" ";
106          }
107          out<<"\n";
108   }
109   out<<"\n"; 
110   for(i=0; i<imax; i++) {
111     for(j=0; j<jmax; j++) {
112            out<<(*fmrad5)(i,j)<<" ";
113          }
114          out<<"\n";
115   }
116   out<<"\n";
117   for(i=0; i<imax; i++) {
118     for(j=0; j<jmax; j++) {
119            out<<(*fmrad6)(i,j)<<" ";
120          }
121          out<<"\n";
122   }
123   
124   out.close(); 
125 ///////////////////////////////////////////////////////////////////////////////////////////////////////////      
126 */
127 }
128
129 AliITSRad::~AliITSRad() {
130
131   delete fmrad1;
132   delete fmrad2;
133   delete fmrad3;
134   delete fmrad4;
135   delete fmrad5;
136   delete fmrad6;
137   
138 }