STEER/AliMagF.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 // Basic magnetic field class
  20 // Used in all the detectors, and also in the traking classes
  21 // Author:
  22 //----------------------------------------------------------------------
  23
  24 #include "AliLog.h"
  25 #include "AliMagF.h"
  26
  27 ClassImp(AliMagF)
  28
  29 //_______________________________________________________________________
  30 AliMagF::AliMagF():
  31   fMap(0),
  32   fType(0),
  33   fInteg(0),
  34   fPrecInteg(1),
  35   fFactor(0),
  36   fMax(0),
  37   fReadField(1)
  38 {
  39   //
  40   // Default constructor
  41   //
  42 }
  43
  44 //_______________________________________________________________________
  45 AliMagF::AliMagF(const char *name, const char *title, Int_t integ,
  46                  Float_t factor, Float_t fmax):
  47   TNamed(name,title),
  48   fMap(0),
  49   fType(0),
  50   fInteg(0),
  51   fPrecInteg(1),
  52   fFactor(factor),
  53   fMax(fmax),
  54   fReadField(1)
  55 {
  56   //
  57   // Standard constructor
  58   //
  59     if(integ<0 || integ > 2) {
  60       AliWarning(Form(
  61               "Invalid magnetic field flag: %5d; Helix tracking chosen instead"
  62               ,integ));
  63       fInteg = 2;
  64     } else {
  65       fInteg = integ;
  66     }
  67
  68     if (fInteg == 0) fPrecInteg = 0;
  69
  70     fType = kUndef;
  71     //
  72 }
  73
  74 //_______________________________________________________________________
  75 AliMagF::AliMagF(const AliMagF &src):
  76   TNamed(src),
  77   fMap(src.fMap),
  78   fType(src.fType),
  79   fInteg(src.fInteg),
  80   fPrecInteg(src.fPrecInteg),
  81   fFactor(src.fFactor),
  82   fMax(src.fMax),
  83   fReadField(src.fReadField)
  84 {
  85     // Copy constructor
  86 }
  87
  88 //_______________________________________________________________________
  89 void AliMagF::Field(Float_t*, Float_t *b) const
  90 {
  91   //
  92   // Method to return the field in one point -- dummy in this case
  93   //
  94   AliWarning("Undefined MagF Field called, returning 0");
  95   b[0]=b[1]=b[2]=0;
  96 }
  97
  98 //_______________________________________________________________________
  99 void AliMagF::GetTPCInt(Float_t *, Float_t *b) const
 100 {
 101 //
 102 // Obtain the integral of the field components in the TPC from given point
 103 // to the closest cathod plane
 104 //
 105   AliWarning("Undefined MagF TPCIntegral called, returning 0");
 106   b[0]=b[1]=b[2]=0;
 107 }
 108
 109 //_______________________________________________________________________
 110 void AliMagF::GetTPCIntCyl(Float_t *, Float_t *b) const
 111 {
 112 //
 113 // Obtain the integral of the field components in the TPC from given point
 114 // to the closest cathod plane
 115 //
 116   AliWarning("Undefined MagF TPCIntegral called, returning 0");
 117   b[0]=b[1]=b[2]=0;
 118 }
 119
 120 //_______________________________________________________________________
 121 AliMagF& AliMagF::operator=(const AliMagF& rhs)
 122 {
 123     // Asignment operator
 124     fMap       = rhs.fMap;
 125     fType      = rhs.fType;
 126     fInteg     = rhs.fInteg;
 127     fPrecInteg = rhs.fPrecInteg;
 128     fFactor    = rhs.fFactor;
 129     fMax       = rhs.fMax;
 130     fReadField = rhs.fReadField;
 131     return *this;
 132 }
 133
 134 void AliMagF::SetPrecInteg(Int_t integ)
 135 {
 136     if (fInteg > 0) {
 137         fPrecInteg = integ;
 138     }
 139     else if (integ != 0)
 140     {
 141         AliWarning("Precision integration flag set to 0 \n");
 142         fPrecInteg = 0;
 143     }
 144 }