]> git.uio.no Git - u/mrichter/AliRoot.git/blob - STEER/AliMagF.cxx
Additional protection: no update of V0 momentum without mag. field
[u/mrichter/AliRoot.git] / 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     fType = kUndef;
68     //
69 }
70
71 //_______________________________________________________________________
72 AliMagF::AliMagF(const AliMagF &src):
73   TNamed(src),
74   fMap(src.fMap),
75   fType(src.fType),
76   fInteg(src.fInteg),
77   fPrecInteg(src.fPrecInteg),
78   fFactor(src.fFactor),
79   fMax(src.fMax),
80   fReadField(src.fReadField)
81 {
82     // Copy constructor
83 }
84
85 //_______________________________________________________________________
86 void AliMagF::Field(Float_t*, Float_t *b) const
87 {
88   //
89   // Method to return the field in one point -- dummy in this case
90   //
91   AliWarning("Undefined MagF Field called, returning 0");
92   b[0]=b[1]=b[2]=0;
93 }
94
95 //_______________________________________________________________________
96 void AliMagF::GetTPCInt(Float_t *, Float_t *b) const
97 {
98 //
99 // Obtain the integral of the field components in the TPC from given point
100 // to the closest cathod plane
101 //
102   AliWarning("Undefined MagF TPCIntegral called, returning 0");
103   b[0]=b[1]=b[2]=0;
104 }
105
106 //_______________________________________________________________________
107 void AliMagF::GetTPCIntCyl(Float_t *, Float_t *b) const
108 {
109 //    
110 // Obtain the integral of the field components in the TPC from given point
111 // to the closest cathod plane
112 //
113   AliWarning("Undefined MagF TPCIntegral called, returning 0");
114   b[0]=b[1]=b[2]=0;
115 }
116
117 //_______________________________________________________________________
118 AliMagF& AliMagF::operator=(const AliMagF& rhs)
119 {
120     // Asignment operator
121     fMap       = rhs.fMap;
122     fType      = rhs.fType;
123     fInteg     = rhs.fInteg;
124     fPrecInteg = rhs.fPrecInteg;
125     fFactor    = rhs.fFactor;
126     fMax       = rhs.fMax;
127     fReadField = rhs.fReadField;
128     return *this;
129 }