Methods to obtain the integral of the field components in the TPC from given point

to the closest cathod plane added. (R. Shahoyan)

diff --git a/STEER/AliMagF.cxx b/STEER/AliMagF.cxx
index 606de825928ac92e33d67bf51df068aadb13b59c..41f3c1b5d02d9d1144493d4a41dae4651bcd9a41 100644 (file)
--- a/STEER/AliMagF.cxx
+++ b/STEER/AliMagF.cxx
@@ -68,6 +68,20 @@ AliMagF::AliMagF(const char *name, const char *title, Int_t integ,
     //
 }
 
+//_______________________________________________________________________
+AliMagF::AliMagF(const AliMagF &src):
+  TNamed(src),
+  fMap(src.fMap),
+  fType(src.fType),
+  fInteg(src.fInteg),
+  fPrecInteg(src.fPrecInteg),
+  fFactor(src.fFactor),
+  fMax(src.fMax),
+  fReadField(src.fReadField)
+{
+    // Copy constructor
+}
+
 //_______________________________________________________________________
 void AliMagF::Field(Float_t*, Float_t *b) const
 {
@@ -77,3 +91,39 @@ void AliMagF::Field(Float_t*, Float_t *b) const
   AliWarning("Undefined MagF Field called, returning 0");
   b[0]=b[1]=b[2]=0;
 }
+
+//_______________________________________________________________________
+void AliMagF::GetTPCInt(Float_t *, Float_t *b) const
+{
+//
+// Obtain the integral of the field components in the TPC from given point
+// to the closest cathod plane
+//
+  AliWarning("Undefined MagF TPCIntegral called, returning 0");
+  b[0]=b[1]=b[2]=0;
+}
+
+//_______________________________________________________________________
+void AliMagF::GetTPCIntCyl(Float_t *, Float_t *b) const
+{
+//    
+// Obtain the integral of the field components in the TPC from given point
+// to the closest cathod plane
+//
+  AliWarning("Undefined MagF TPCIntegral called, returning 0");
+  b[0]=b[1]=b[2]=0;
+}
+
+//_______________________________________________________________________
+AliMagF& AliMagF::operator=(const AliMagF& rhs)
+{
+    // Asignment operator
+    fMap       = rhs.fMap;
+    fType      = rhs.fType;
+    fInteg     = rhs.fInteg;
+    fPrecInteg = rhs.fPrecInteg;
+    fFactor    = rhs.fFactor;
+    fMax       = rhs.fMax;
+    fReadField = rhs.fReadField;
+    return *this;
+}

diff --git a/STEER/AliMagF.h b/STEER/AliMagF.h
index 1665243df6e52c79a8d0d17a33148c8151bc54e2..d9e3250fce807ddf7c4e0710d99093c9d9942785 100644 (file)
--- a/STEER/AliMagF.h
+++ b/STEER/AliMagF.h
@@ -21,18 +21,22 @@ public:
   AliMagF();
   AliMagF(const char *name, const char *title, Int_t integ, 
          Float_t factor = 1., Float_t fmax = 10.);
+  AliMagF(const AliMagF& maps);
   virtual ~AliMagF() {}
-  virtual void Field(Float_t *x, Float_t *b) const;
-  virtual Int_t Type() const {return fType;}
+  AliMagF& operator=(const AliMagF& rhs);
+  virtual void    Field(Float_t *x, Float_t *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;}
   virtual Float_t Max() const {return fMax;}
-  virtual Int_t Map() const {return fMap;}
-  virtual Int_t Integ() const {return fInteg;}
-  virtual Int_t PrecInteg() const {return fPrecInteg;}  
+  virtual Int_t   Map() const {return fMap;}
+  virtual Int_t   Integ() const {return fInteg;}
+  virtual Int_t   PrecInteg() const {return fPrecInteg;}  
   virtual Float_t Factor() const {return fFactor;}
-  virtual void ReadField() {}
+  virtual void    ReadField() {}
   virtual Float_t SolenoidField() const {return 2.;}
-  virtual void SetPrecInteg(Int_t integ) {fPrecInteg = integ;}
-  virtual void SetReadField(Bool_t flag = kTRUE) {fReadField = flag;}
+  virtual void    SetPrecInteg(Int_t integ) {fPrecInteg = integ;}
+  virtual void    SetReadField(Bool_t flag = kTRUE) {fReadField = flag;}
  protected:
   Int_t     fMap;       // Field Map identifier
   Int_t     fType;      // Mag Field type