Removing warnings (icc)
[u/mrichter/AliRoot.git] / STEER / AliMagFMaps.cxx
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$ */
17
18 //------------------------------------------------------------------------
19 // Magnetic field composed by 3 maps: the L3 magnet, extended region, and
20 // dipole magnet
21 // Used in the configuration macros (macros/Config.C, etc.)
22 // Author: Andreas Morsch <andreas.morsch@cern.ch>
23 //------------------------------------------------------------------------
24
25 #include <TFile.h>
26 #include <TSystem.h>
27 #include <TVector.h>
28
29 #include "AliFieldMap.h"
30 #include "AliMagFMaps.h"
31
32 ClassImp(AliMagFMaps)
33
34 //_______________________________________________________________________
35 AliMagFMaps::AliMagFMaps():
36   fSolenoid(0),
37   fSolenoidUser(0.),
38   fL3Option(0),
39   fFieldRead(0)
40 {
41   //
42   // Default constructor
43   //
44   fFieldMap[0] = fFieldMap[1] = fFieldMap[2] = 0;
45 }
46
47 //_______________________________________________________________________
48 AliMagFMaps::AliMagFMaps(const char *name, const char *title, Int_t integ, 
49                          Float_t factor, Float_t fmax, Int_t map, 
50                          Int_t l3):
51   AliMagFC(name,title,integ,factor,fmax),
52   fSolenoid(0),
53   fSolenoidUser(0),
54   fL3Option(l3),
55   fFieldRead(0)
56 {
57   //
58   // Standard constructor
59   //
60   fType         = kConMesh;
61   fFieldMap[0]  = 0;
62   fMap          = map;
63   fL3Option     = l3;
64   ReadField();
65   fFieldRead = 1;
66   //
67   // Don't replicate field information in gAlice
68   for (Int_t i = 0; i < 3; i++)  fFieldMap[i]->SetWriteEnable(0);
69   //
70 }
71
72 //_______________________________________________________________________
73 AliMagFMaps::AliMagFMaps(const AliMagFMaps &magf):
74   AliMagFC(magf),
75   fSolenoid(0),
76   fL3Option(0),
77   fFieldRead(0)
78 {
79   //
80   // Copy constructor
81   //
82   magf.Copy(*this);
83 }
84
85 //_______________________________________________________________________
86 AliMagFMaps::~AliMagFMaps()
87 {
88   //
89   //  Destructor
90   //
91   delete fFieldMap[0];
92   delete fFieldMap[1];
93   delete fFieldMap[2];    
94 }
95
96 //_______________________________________________________________________
97 void AliMagFMaps::ReadField()
98 {
99   //  Read Field Map from file
100   //
101   //  don't read twice
102   //
103   if (fFieldRead) return;
104   fFieldRead = 1;
105   //    
106   char* fname;
107   TFile* file = 0;
108   if (fMap == k2kG) {
109       fSolenoid = 2.;
110       fname = gSystem->ExpandPathName("$(ALICE_ROOT)/data/maps/L3B02.root");
111       file = new TFile(fname);
112       fFieldMap[0] = dynamic_cast<AliFieldMap*>(file->Get("L3B02"));
113       file->Close();
114       delete file;
115       fname = gSystem->ExpandPathName("$(ALICE_ROOT)/data/maps/DipB02.root");
116       file = new TFile(fname);
117       fFieldMap[1] = dynamic_cast<AliFieldMap*>(file->Get("DipB02"));
118       file->Close();
119       delete file;;
120       
121       fname = gSystem->ExpandPathName("$(ALICE_ROOT)/data/maps/ExtB02.root");
122       file = new TFile(fname);
123       fFieldMap[2] = dynamic_cast<AliFieldMap*>(file->Get("ExtB02"));
124       file->Close();
125       delete file;
126   } else if (fMap == k4kG) {
127       fSolenoid = 4.;
128       fname = gSystem->ExpandPathName("$(ALICE_ROOT)/data/maps/L3B04.root");
129       file = new TFile(fname);
130       fFieldMap[0] = dynamic_cast<AliFieldMap*>(file->Get("L3B04"));
131       file->Close();
132       delete file;
133       fname = gSystem->ExpandPathName("$(ALICE_ROOT)/data/maps/DipB04.root");
134       file = new TFile(fname);
135       fFieldMap[1] = dynamic_cast<AliFieldMap*>(file->Get("DipB04"));
136       file->Close();
137       delete file;
138       
139       fname = gSystem->ExpandPathName("$(ALICE_ROOT)/data/maps/ExtB04.root");
140       file = new TFile(fname);
141       fFieldMap[2] = dynamic_cast<AliFieldMap*>(file->Get("ExtB04"));
142       file->Close();
143       delete file;
144   } else if (fMap == k5kG) {
145       fSolenoid = 5.;
146       fname = gSystem->ExpandPathName("$(ALICE_ROOT)/data/maps/L3B05.root");
147       file = new TFile(fname);
148       fFieldMap[0] = dynamic_cast<AliFieldMap*>(file->Get("L3B05"));
149       file->Close();
150       delete file;
151       fname = gSystem->ExpandPathName("$(ALICE_ROOT)/data/maps/DipB05.root");
152       file = new TFile(fname);
153       fFieldMap[1] = dynamic_cast<AliFieldMap*>(file->Get("DipB05"));
154       file->Close();
155       delete file;
156       
157       fname = gSystem->ExpandPathName("$(ALICE_ROOT)/data/maps/ExtB05.root");
158       file = new TFile(fname);
159       fFieldMap[2] = dynamic_cast<AliFieldMap*>(file->Get("ExtB05"));
160       file->Close();
161       delete file;
162   }
163
164   if (!fL3Option) fSolenoidUser = fSolenoid;
165
166 }
167
168 //_______________________________________________________________________
169 Float_t AliMagFMaps::SolenoidField() const
170 {
171   //
172   // Returns max. L3 (solenoid) field strength 
173   // according to field map setting 
174   //
175   return fSolenoid;
176 }
177
178 //_______________________________________________________________________
179 void AliMagFMaps::Field(Float_t *x, Float_t *b)
180 {
181   //
182   // Method to calculate the magnetic field
183   //
184   // --- find the position in the grid ---
185   
186
187     
188     if (!fFieldRead) ReadField();
189
190   //
191   // Field Maps have been calculated for the coordinate system in which 
192   // the Muon Spectrometer is placed at z > 0
193   // Transform coordinates corresponingly 
194   //
195   
196   b[0]=b[1]=b[2]=0;
197   Float_t xm[3];
198   xm[0] = - x[0];
199   xm[1] =   x[1];
200   xm[2] = - x[2];
201   
202   AliFieldMap* map = 0;
203   if (fFieldMap[0]->Inside(xm[0], xm[1], xm[2])) {
204       map = fFieldMap[0];
205       if (!fL3Option && TMath::Abs(xm[2]) < 370.) {
206       //
207       //     Constant L3 field , if this option was selected
208       //
209           b[2] = - fSolenoid;
210           return;
211       } 
212   } else if (fFieldMap[1]->Inside(xm[0], xm[1], xm[2])) {
213     map = fFieldMap[1];
214   } else if (fFieldMap[2]->Inside(xm[0], xm[1], xm[2])) {
215     map = fFieldMap[2];
216   }
217   
218   if(map){
219     map->Field(xm,b);
220     b[0] = - b[0];
221     b[2] = - b[2];
222
223   } else {
224       //This is the ZDC part
225       ZDCField(x, b);
226   }
227
228   
229   if(fFactor!=1) {
230       b[0]*=fFactor;
231       b[1]*=fFactor;
232       b[2]*=fFactor;
233   }
234 }
235
236 //_______________________________________________________________________
237 void AliMagFMaps::Copy(TObject & /* magf */) const
238 {
239   //
240   // Copy *this onto magf -- Not implemented
241   //
242   Fatal("Copy","Not implemented!\n");
243 }
244
245 //_______________________________________________________________________
246 void AliMagFMaps::Streamer(TBuffer &R__b)
247 {
248   // Stream an object of class AliMagFMaps.
249   if (R__b.IsReading()) {
250     AliMagFMaps::Class()->ReadBuffer(R__b, this);
251     fFieldRead = 0;
252   } else {
253     AliMagFMaps::Class()->WriteBuffer(R__b, this);
254   }
255 }