Transition to NewIO

[u/mrichter/AliRoot.git] / STEER / AliMagFCM.cxx

diff --git a/STEER/AliMagFCM.cxx b/STEER/AliMagFCM.cxx
index 3810a0b080898a02be26a7caae824dadbfc71c83..b11a19e95146f6e07a7d8fcdb1033ba295f6d15b 100644 (file)
--- a/STEER/AliMagFCM.cxx
+++ b/STEER/AliMagFCM.cxx
@@ -13,48 +13,94 @@
  * provided "as is" without express or implied warranty.                  *
  **************************************************************************/
 
-/*
-$Log$
-Revision 1.5  2000/12/01 11:20:27  alibrary
-Corrector dipole removed from ZDC
+/* $Id$ */
 
-Revision 1.4  2000/11/30 07:12:49  alibrary
-Introducing new Rndm and QA classes
-
-Revision 1.3  2000/11/10 18:09:55  fca
-New field map for the ZDC
-
-Revision 1.2  2000/07/12 08:56:25  fca
-Coding convention correction and warning removal
-
-Revision 1.1  2000/07/11 18:24:59  fca
-Coding convention corrections + few minor bug fixes
+//-----------------------------------------------------------------------
+//  Class for Alice magnetic field with constant mesh
+//  Used in the configuration macros (macros/Config.C, etc.)
+//  Author:
+//-----------------------------------------------------------------------
+#include "TSystem.h"
 
-*/
 #include "TVector.h"
 
 #include "AliMagFCM.h"
-#include "TSystem.h"
 
 ClassImp(AliMagFCM)
 
-//________________________________________
+//_______________________________________________________________________
+AliMagFCM::AliMagFCM():
+  fXbeg(0),
+  fYbeg(0),
+  fZbeg(0),
+  fXdel(0),
+  fYdel(0),
+  fZdel(0),
+  fSolenoid(0),
+  fXdeli(0),
+  fYdeli(0),
+  fZdeli(0),
+  fXn(0),
+  fYn(0),
+  fZn(0),
+  fB(0)
+{
+  //
+  // Standard constructor
+  //
+  fType = kConMesh;
+  fMap  = 2;
+  SetSolenoidField();
+}
+
+//_______________________________________________________________________
 AliMagFCM::AliMagFCM(const char *name, const char *title, const Int_t integ, 
-                    const Float_t factor, const Float_t fmax)
-  : AliMagF(name,title,integ,factor,fmax)
+                     const Float_t factor, const Float_t fmax):
+  AliMagF(name,title,integ,factor,fmax),
+  fXbeg(0),
+  fYbeg(0),
+  fZbeg(0),
+  fXdel(0),
+  fYdel(0),
+  fZdel(0),
+  fSolenoid(0),
+  fXdeli(0),
+  fYdeli(0),
+  fZdeli(0),
+  fXn(0),
+  fYn(0),
+  fZn(0),
+  fB(0)
 {
   //
   // Standard constructor
   //
   fType = kConMesh;
   fMap  = 2;
-  
-  printf("Constant Mesh Field %s created: map= %d, factor= %f, file= %s\n",
-        fName.Data(), fMap, factor,fTitle.Data());
+  SetSolenoidField();
+
+  if(fDebug>-1) Info("ctor",
+     "%s: Constant Mesh Field %s created: map= %d, factor= %f, file= %s\n",
+        ClassName(),fName.Data(), fMap, factor,fTitle.Data());
 }
 
-//________________________________________
-AliMagFCM::AliMagFCM(const AliMagFCM &magf)
+//_______________________________________________________________________
+AliMagFCM::AliMagFCM(const AliMagFCM &magf):
+  AliMagF(magf),
+  fXbeg(0),
+  fYbeg(0),
+  fZbeg(0),
+  fXdel(0),
+  fYdel(0),
+  fZdel(0),
+  fSolenoid(0),
+  fXdeli(0),
+  fYdeli(0),
+  fZdeli(0),
+  fXn(0),
+  fYn(0),
+  fZn(0),
+  fB(0)
 {
   //
   // Copy constructor
@@ -62,7 +108,7 @@ AliMagFCM::AliMagFCM(const AliMagFCM &magf)
   magf.Copy(*this);
 }
 
-//________________________________________
+//_______________________________________________________________________
 void AliMagFCM::Field(Float_t *x, Float_t *b)
 {
   //
@@ -77,7 +123,7 @@ void AliMagFCM::Field(Float_t *x, Float_t *b)
   
   b[0]=b[1]=b[2]=0;
   if(-700<x[2] && x[2]<fZbeg && x[0]*x[0]+(x[1]+30)*(x[1]+30) < 560*560) {
-    b[2]=2;
+    b[2]= fSolenoid;
   } else  {
     Bool_t infield=(fZbeg<=x[2] && x[2]<fZbeg+fZdel*(fZn-1)
                    &&  ( fXbeg <= TMath::Abs(x[0]) && TMath::Abs(x[0]) < fXbeg+fXdel*(fXn-1) )
@@ -191,7 +237,7 @@ void AliMagFCM::Field(Float_t *x, Float_t *b)
   }
 }
 
-//________________________________________
+//_______________________________________________________________________
 void AliMagFCM::ReadField()
 {
   // 
@@ -201,15 +247,15 @@ void AliMagFCM::ReadField()
   Int_t ix, iy, iz, ipx, ipy, ipz;
   Float_t bx, by, bz;
   char *fname;
-  printf("Reading Magnetic Field %s from file %s\n",fName.Data(),fTitle.Data());
+  if(fDebug) printf("%s: Reading Magnetic Field %s from file %s\n",ClassName(),fName.Data(),fTitle.Data());
   fname = gSystem->ExpandPathName(fTitle.Data());
   magfile=fopen(fname,"r");
   delete [] fname;
   if (magfile) {
     fscanf(magfile,"%d %d %d %f %f %f %f %f %f",
           &fXn, &fYn, &fZn, &fXdel, &fYdel, &fZdel, &fXbeg, &fYbeg, &fZbeg);
-    printf("fXn %d, fYn %d, fZn %d, fXdel %f, fYdel %f, fZdel %f, fXbeg %f, fYbeg %f, fZbeg %f\n",
-          fXn, fYn, fZn, fXdel, fYdel, fZdel, fXbeg, fYbeg, fZbeg);
+    if(fDebug>1) printf("%s: fXn %d, fYn %d, fZn %d, fXdel %f, fYdel %f, fZdel %f, fXbeg %f, fYbeg %f, fZbeg %f\n",
+                       ClassName(),fXn, fYn, fZn, fXdel, fYdel, fZdel, fXbeg, fYbeg, fZbeg);
     fXdeli=1./fXdel;
     fYdeli=1./fYdel;
     fZdeli=1./fZdel;
@@ -228,12 +274,12 @@ void AliMagFCM::ReadField()
       }
     }
   } else { 
-    printf("File %s not found !\n",fTitle.Data());
+    printf("%s: File %s not found !\n",ClassName(),fTitle.Data());
     exit(1);
   }
 }
 
-//________________________________________
+//_______________________________________________________________________
 void AliMagFCM::Copy(AliMagFCM & /* magf */) const
 {
   //
@@ -242,9 +288,3 @@ void AliMagFCM::Copy(AliMagFCM & /* magf */) const
   Fatal("Copy","Not implemented!\n");
 }
 
-//________________________________________
-AliMagFCM & AliMagFCM::operator =(const AliMagFCM &magf)
-{
-  magf.Copy(*this);
-  return *this;
-}