Use const for mag field input
authorfca <fca@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 17 Dec 2008 10:20:31 +0000 (10:20 +0000)
committerfca <fca@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 17 Dec 2008 10:20:31 +0000 (10:20 +0000)
16 files changed:
STEER/AliMagF.cxx
STEER/AliMagF.h
STEER/AliMagFC.cxx
STEER/AliMagFC.h
STEER/AliMagFCM.cxx
STEER/AliMagFCM.h
STEER/AliMagFCheb.cxx
STEER/AliMagFCheb.h
STEER/AliMagFDM.cxx
STEER/AliMagFDM.h
STEER/AliMagFMaps.cxx
STEER/AliMagFMaps.h
STEER/AliMagFMapsV1.cxx
STEER/AliMagFMapsV1.h
STEER/AliMagWrapCheb.cxx
STEER/AliMagWrapCheb.h

index a3cf3e2353525d283721b5928b05b69054b6df63..7b00a359f8cfa85b832a491a6614d603f882ed4f 100644 (file)
@@ -86,7 +86,7 @@ AliMagF::AliMagF(const AliMagF &src):
 }
 
 //_______________________________________________________________________
-void AliMagF::Field(float*, float *b) const
+void AliMagF::Field(const float*, float *b) const
 {
   //
   // Method to return the field in one point -- dummy in this case
@@ -96,7 +96,7 @@ void AliMagF::Field(float*, float *b) const
 }
 
 //_______________________________________________________________________
-void AliMagF::Field(double*, double *b) const
+void AliMagF::Field(const double*, double *b) const
 {
   //
   // Method to return the field in one point -- dummy in this case
index c360d53576e2b9266a4999ee4966848628f658a7..d7690827ad1647e21c42ad41734ffc3c55ce5f15 100644 (file)
@@ -24,8 +24,8 @@ public:
   AliMagF(const AliMagF& maps);
   virtual ~AliMagF() {}
   AliMagF& operator=(const AliMagF& rhs);
-  virtual void    Field(float *x, float *b)                  const;
-  virtual void    Field(double *x, double *b)                const;
+  virtual void    Field(const float *x, float *b)                  const;
+  virtual void    Field(const double *x, double *b)                const;
   virtual void    GetTPCInt(Float_t *xyz, Float_t *b)        const;
   virtual void    GetTPCIntCyl(Float_t *rphiz, Float_t *b)   const;
   virtual Int_t   Type() const {return fType;}
index 7f2897999101e451402958201b7d9ca64c74d26d..0cd875b16c7e91f45cc19982207232447364c43b 100644 (file)
@@ -68,7 +68,7 @@ AliMagFC::AliMagFC(const char *name, const char *title, Int_t integ,
 }
 
 //________________________________________
-void AliMagFC::Field(float *x, float *b) const
+void AliMagFC::Field(const float *x, float *b) const
 {
   //
   // Method to return the field in a point
@@ -99,7 +99,7 @@ void AliMagFC::Field(float *x, float *b) const
   }
 }
 //________________________________________
-void AliMagFC::Field(double *x, double *b) const
+void AliMagFC::Field(const double *x, double *b) const
 {
   //
   // Method to return the field in a point
@@ -131,7 +131,7 @@ void AliMagFC::Field(double *x, double *b) const
 }
 
 //___________________________________________________
-void AliMagFC::ZDCField(float *x, float *b) const
+void AliMagFC::ZDCField(const float *x, float *b) const
 {
   // ---- This is the ZDC part
   
@@ -275,7 +275,7 @@ void AliMagFC::ZDCField(float *x, float *b) const
   }
 }
 
-void AliMagFC::ZDCField(double *x, double *b) const
+void AliMagFC::ZDCField(const double *x, double *b) const
 {
   // ---- This is the ZDC part
   
index 55cc14bcf51645100198b82c177e98bcb0209447..d17e0de1931a993e214de06558b1295bc89901a6 100644 (file)
@@ -23,11 +23,11 @@ public:
   AliMagFC(const char *name, const char *title, Int_t integ, 
           Float_t factor, Float_t fmax);
   virtual ~AliMagFC(){}
-  virtual void Field(float *x, float *b)      const;
-  virtual void Field(double *x, double *b)    const;
+  virtual void Field(const float *x, float *b)      const;
+  virtual void Field(const double *x, double *b)    const;
   virtual void ReadField() {}
-  virtual void ZDCField(float *x, float *b)   const;
-  virtual void ZDCField(double *x, double *b) const;
+  virtual void ZDCField(const float *x, float *b)   const;
+  virtual void ZDCField(const double *x, double *b) const;
   virtual void SetBeamType(BeamType_t type)      {fBeamType    = type;}
   virtual void SetBeamEnergy(Float_t energy)     {fBeamEnergy  = energy;}
   virtual void SetCompensatorMagnet(Bool_t flag) {fCompensator = flag;}
index cd11c36f873efc0beb24363eb1e9d274af51609e..e3b8759a071bf119c99e93c9e7fc899bbb581507 100644 (file)
@@ -111,7 +111,7 @@ AliMagFCM::AliMagFCM(const AliMagFCM &magf):
 }
 
 //_______________________________________________________________________
-void AliMagFCM::Field(float *x, float *b) const
+void AliMagFCM::Field(const float *x, float *b) const
 {
   //
   // Method to calculate the magnetic field
@@ -211,7 +211,7 @@ void AliMagFCM::Field(float *x, float *b) const
 }
 
 //_______________________________________________________________________
-void AliMagFCM::Field(double *x, double *b) const
+void AliMagFCM::Field(const double *x, double *b) const
 {
   //
   // Method to calculate the magnetic field
index 77e353155626e2a9f8317ca62d1c1ddc71bbdb03..fb774f3cc0737388ee6f5215f1e321e37ab4ccb3 100644 (file)
@@ -24,8 +24,8 @@ public:
           Float_t factor, Float_t fmax);
   AliMagFCM(const AliMagFCM &mag);
   virtual ~AliMagFCM() {delete fB;}
-  virtual void Field(float *x, float *b) const;
-  virtual void Field(double *x, double *b) const;
+  virtual void Field(const float *x, float *b) const;
+  virtual void Field(const double *x, double *b) const;
   virtual void ReadField();
   virtual void    SetSolenoidField(Float_t field = 2.) {fSolenoid = field;}
   virtual Float_t SolenoidField() const {
index 058a8777bfd6adc71aa0af81c3b36dac436fc57c..a98b83ac3308e358dc474aef18cab35df72a6397 100644 (file)
@@ -237,7 +237,7 @@ void AliMagFCheb::Clear(const Option_t *)
 }
 
 //__________________________________________________________________________________________
-void AliMagFCheb::Field(float *xyz, float *b) const
+void AliMagFCheb::Field(const float *xyz, float *b) const
 {
   // compute field in cartesian coordinates. If point is outside of the parameterized region
   // get it at closest valid point
@@ -272,7 +272,7 @@ void AliMagFCheb::Field(float *xyz, float *b) const
 }
 
 //__________________________________________________________________________________________
-void AliMagFCheb::Field(double *xyz, double *b) const
+void AliMagFCheb::Field(const double *xyz, double *b) const
 {
   // compute field in cartesian coordinates. If point is outside of the parameterized region
   // get it at closest valid point
@@ -307,7 +307,7 @@ void AliMagFCheb::Field(double *xyz, double *b) const
 }
 
 //__________________________________________________________________________________________
-void AliMagFCheb::GetTPCInt(Float_t *xyz, Float_t *b) const
+void AliMagFCheb::GetTPCInt(const Float_t *xyz, Float_t *b) const
 {
   // compute TPC region field integral in cartesian coordinates.
   // If point is outside of the parameterized region get it at closeset valid point
index d5f34debe45d05ab883e2c971d848f73d384cda8..724f0e5585a80deaab1839c65ae0d2372d208eb4 100644 (file)
@@ -56,7 +56,7 @@ class AliMagFCheb: public TNamed
   //
   Int_t      GetNParamsSol()                              const {return fNParamsSol;}
   Int_t      GetNSegZSol()                                const {return fNSegZSol;}
-  float*     GetSegZSol() const {return fSegZSol;}
+  Float_t*     GetSegZSol() const {return fSegZSol;}
   //
   Int_t      GetNParamsTPCInt()                           const {return fNParamsTPCInt;}
   Int_t      GetNSegZTPCInt()                             const {return fNSegZTPCInt;}
@@ -82,14 +82,14 @@ class AliMagFCheb: public TNamed
   //
   virtual void Print(Option_t * = "")                     const;
   //
-  virtual void Field(float *xyz, float *b)                const;
-  virtual void Field(double *xyz, double *b)              const;
+  virtual void Field(const Float_t *xyz, Float_t *b)                const;
+  virtual void Field(const Double_t *xyz, Double_t *b)              const;
   //
-  virtual void FieldCyl(const float *rphiz, float *b)     const;
-  virtual void FieldCyl(const double *rphiz, double *b)   const;
+  virtual void FieldCyl(const Float_t *rphiz, Float_t *b)     const;
+  virtual void FieldCyl(const Double_t *rphiz, Double_t *b)   const;
   //
-  virtual void GetTPCInt(Float_t *xyz, Float_t *b)        const;
-  virtual void GetTPCIntCyl(Float_t *rphiz, Float_t *b)   const;
+  virtual void GetTPCInt(const Float_t *xyz, Float_t *b)        const;
+  virtual void GetTPCIntCyl(const Float_t *rphiz, Float_t *b)   const;
   //
   template <class T>
     Int_t      FindDipSegment(const T *xyz)               const; 
@@ -112,8 +112,8 @@ class AliMagFCheb: public TNamed
   //
   AliMagFCheb(const char* inputFile);
   void       SaveData(const char* outfile)              const;
-  Int_t      SegmentDipDimension(float** seg,const TObjArray* par,int npar, int dim, 
-                                float xmn,float xmx,float ymn,float ymx,float zmn,float zmx);
+  Int_t      SegmentDipDimension(Float_t** seg,const TObjArray* par,int npar, int dim, 
+                                Float_t xmn,Float_t xmx,Float_t ymn,Float_t ymx,Float_t zmn,Float_t zmx);
   //
   void       AddParamSol(const AliCheb3D* param);
   void       AddParamTPCInt(const AliCheb3D* param);
@@ -127,8 +127,8 @@ class AliMagFCheb: public TNamed
 #endif
   //
  protected:
-    virtual void FieldCylSol(const float *rphiz, float *b)      const;
-    virtual void FieldCylSol(const double *rphiz, double *b)    const;
+    virtual void FieldCylSol(const Float_t *rphiz, Float_t *b)      const;
+    virtual void FieldCylSol(const Double_t *rphiz, Double_t *b)    const;
   //
  protected:
   //
@@ -186,7 +186,7 @@ class AliMagFCheb: public TNamed
 
 
 //__________________________________________________________________________________________
-inline void AliMagFCheb::FieldCyl(const float *rphiz, float *b) const
+inline void AliMagFCheb::FieldCyl(const Float_t *rphiz, Float_t *b) const
 {
   // compute field in Cylindircal coordinates
   //  if (rphiz[2]<GetMinZSol() || rphiz[2]>GetMaxZSol() || rphiz[0]>GetMaxRSol()) {for (int i=3;i--;) b[i]=0; return;}
@@ -195,7 +195,7 @@ inline void AliMagFCheb::FieldCyl(const float *rphiz, float *b) const
 
 
 //__________________________________________________________________________________________
-inline void AliMagFCheb::FieldCyl(const double *rphiz, double *b) const
+inline void AliMagFCheb::FieldCyl(const Double_t *rphiz, Double_t *b) const
 {
   // compute field in Cylindircal coordinates
   //  if (rphiz[2]<GetMinZSol() || rphiz[2]>GetMaxZSol() || rphiz[0]>GetMaxRSol()) {for (int i=3;i--;) b[i]=0; return;}
@@ -278,7 +278,7 @@ template <class T>
 Int_t    AliMagFCheb::FindDipSegment(const T *xyz) const 
 {
   // find the segment containing point xyz. If it is outside find the closest segment 
-  int xid,yid,zid = TMath::BinarySearch(fNZSegDip,fSegZDip,(float)xyz[2]); // find zsegment
+  int xid,yid,zid = TMath::BinarySearch(fNZSegDip,fSegZDip,(Float_t)xyz[2]); // find zsegment
   int ysegBeg = fBegSegYDip[zid];
   //
   for (yid=0;yid<fNSegYDip[zid];yid++) if (xyz[1]<fSegYDip[ysegBeg+yid]) break;
index 23ff5af41ca167d92c0b9f29a94d98ea1246901d..28c509a86d8d5b8a70d459beb96315c9281b9c8f 100644 (file)
@@ -108,7 +108,7 @@ AliMagFDM::AliMagFDM(const char *name, const char *title, Int_t integ,
 }
 
 //_______________________________________________________________________
-void AliMagFDM::Field(float *xfi, float *b) const
+void AliMagFDM::Field(const float *xfi, float *b) const
 {
   //
   // Main routine to compute the field in a point
@@ -336,7 +336,7 @@ if ((kfZbg/100<xL3[2] && xL3[2]<=zCmin && r0<=rPmax) || ((zCmin<xL3[2] && xL3[2]
 }
 
 //_______________________________________________________________________
-void AliMagFDM::Field(double *xfi, double *b) const
+void AliMagFDM::Field(const double *xfi, double *b) const
 {
   //
   // Main routine to compute the field in a point
index 031e1dd45651d8eb23463387c96d8f20ca6847cc..ae97c34a6d2f5cf93054bc6b794fb843c49ba8b8 100644 (file)
@@ -22,8 +22,8 @@ public:
   AliMagFDM(const char *name, const char *title, Int_t integ,
            Float_t factor, Float_t fmax);
   virtual ~AliMagFDM(){} 
-  virtual void Field(float *x, float *b) const;
-  virtual void Field(double *x, double *b) const;
+  virtual void Field(const float *x, float *b) const;
+  virtual void Field(const double *x, double *b) const;
   virtual void ReadField(); 
   virtual void SetSolenoidField(Float_t field = 2.) {fSolenoid = field;}
   virtual Float_t SolenoidField() const {
index 9051b9210a3060130b7fc535c88dcda32bcd87fc..7db65add35693b42ac9607d3d04a120f07433351 100644 (file)
@@ -187,7 +187,7 @@ Float_t AliMagFMaps::SolenoidField() const
 }
 
 //_______________________________________________________________________
-void AliMagFMaps::Field(float *x, float *b) const
+void AliMagFMaps::Field(const float *x, float *b) const
 {
   //
   // Method to calculate the magnetic field
@@ -245,7 +245,7 @@ void AliMagFMaps::Field(float *x, float *b) const
 }
 
 //_______________________________________________________________________
-void AliMagFMaps::Field(double *x, double *b) const
+void AliMagFMaps::Field(const double *x, double *b) const
 {
   //
   // Method to calculate the magnetic field
index 2d5ced5b18ac6514ee3d7b0b597eb3a2eada9c0c..3256e24cf0b198e0e01c447bd0c261c8e4fff403 100644 (file)
@@ -26,8 +26,8 @@ public:
                Int_t l3 = 1);
     AliMagFMaps(const AliMagFMaps &mag);
     virtual ~AliMagFMaps();
-    virtual void    Field(float  *x, float  *b) const;
-    virtual void    Field(double *x, double *b) const;
+    virtual void    Field(const float  *x, float  *b) const;
+    virtual void    Field(const double *x, double *b) const;
     AliFieldMap* FieldMap(Int_t i) {return fFieldMap[i];}
     virtual void ReadField();
     virtual Float_t SolenoidField() const;
index 0eae1a93a37a3b3f68f93b8e441351f78d678486..6fcc4a1348bb3c21519e4b79e59b6fb52d04b2a8 100644 (file)
@@ -94,7 +94,7 @@ AliMagFMapsV1::~AliMagFMapsV1()
 }
 
 //_______________________________________________________________________
-void AliMagFMapsV1::Field(float *x, float *b) const
+void AliMagFMapsV1::Field(const float *x, float *b) const
 {
   //
   // Method to calculate the magnetic field at position x
@@ -125,7 +125,7 @@ void AliMagFMapsV1::Field(float *x, float *b) const
 }
 
 //_______________________________________________________________________
-void AliMagFMapsV1::Field(double *x, double *b) const
+void AliMagFMapsV1::Field(const double *x, double *b) const
 {
   //
   // Method to calculate the magnetic field at position x
index b28be52001d34d2fbfd4369a2784f64f4449f0ca..b94d6422283f08739b93477e4fb5f5fa18c604bc 100644 (file)
@@ -25,8 +25,8 @@ public:
     AliMagFMapsV1(const AliMagFMapsV1& maps);             
     AliMagFMapsV1& operator=(const AliMagFMapsV1& maps) {maps.Copy(*this); return *this;}
     virtual ~AliMagFMapsV1();
-    virtual void    Field(float *x, float *b) const;
-    virtual void    Field(double *x, double *b) const;
+    virtual void    Field(const float *x, float *b) const;
+    virtual void    Field(const double *x, double *b) const;
     virtual Float_t SolenoidField() const;
     AliMagFCheb* GetMeasuredMap()                   const {return fMeasuredMap;}
     void SetMeasuredMap(AliMagFCheb* parm)               {if (parm) delete parm; fMeasuredMap = parm;}
index d0e397b73c5694c84e31069120cd0d8ceaa9a506..d2937a3fdde41274d3e884c180dd194ea784b43f 100644 (file)
@@ -88,7 +88,7 @@ AliMagWrapCheb::~AliMagWrapCheb()
 }
 
 //_______________________________________________________________________
-void AliMagWrapCheb::GetTPCInt(Float_t *xyz, Float_t *b) const
+void AliMagWrapCheb::GetTPCInt(const Float_t *xyz, Float_t *b) const
 {
   // Method to calculate the integral of magnetic integral from xyz to nearest cathode plane
   //
@@ -98,7 +98,7 @@ void AliMagWrapCheb::GetTPCInt(Float_t *xyz, Float_t *b) const
 }
 
 //_______________________________________________________________________
-void AliMagWrapCheb::GetTPCIntCyl(Float_t *rphiz, Float_t *b) const
+void AliMagWrapCheb::GetTPCIntCyl(const Float_t *rphiz, Float_t *b) const
 {
   // Method to calculate the integral of magnetic integral from point to nearest cathode plane
   // in cylindrical coordiates ( -pi<phi<pi convention )
@@ -108,7 +108,7 @@ void AliMagWrapCheb::GetTPCIntCyl(Float_t *rphiz, Float_t *b) const
 }
 
 //_______________________________________________________________________
-void AliMagWrapCheb::Field(float *xyz, float *b) const
+void AliMagWrapCheb::Field(const Float_t *xyz, Float_t *b) const
 {
   // Method to calculate the field at point  xyz
   //
@@ -124,7 +124,7 @@ void AliMagWrapCheb::Field(float *xyz, float *b) const
 }
 
 //_______________________________________________________________________
-void AliMagWrapCheb::Field(double *xyz, double *b) const
+void AliMagWrapCheb::Field(const Double_t *xyz, Double_t *b) const
 {
   // Method to calculate the field at point  xyz
   //
index d623cf9e78d604296feed5e4d7d547c8a5a46843..189be0e735b7473a2c4be29e2f355370c64cbe3a 100644 (file)
@@ -26,10 +26,10 @@ public:
   AliMagWrapCheb& operator=(const AliMagWrapCheb& maps);
   virtual ~AliMagWrapCheb();
   //
-  virtual void Field(float *x, float *b)                        const;
-  virtual void Field(double *x, double *b)                      const;
-  virtual void GetTPCInt(Float_t *xyz, Float_t *b)              const;
-  virtual void GetTPCIntCyl(Float_t *rphiz, Float_t *b)         const;
+  virtual void Field(const Float_t *x, Float_t *b)                        const;
+  virtual void Field(const Double_t *x, Double_t *b)                      const;
+  virtual void GetTPCInt(const Float_t *xyz, Float_t *b)              const;
+  virtual void GetTPCIntCyl(const Float_t *rphiz, Float_t *b)         const;
   //
   AliMagFCheb* GetMeasuredMap()                                 const {return fMeasuredMap;}
   void SetMeasuredMap(AliMagFCheb* parm)                        {if (fMeasuredMap) delete fMeasuredMap; fMeasuredMap = parm;}