#include <MucMagneticField.h>
Public Member Functions | |
MucMagneticField () | |
~MucMagneticField () | |
void | getMucField (int part, int layer, int mat, HepPoint3D &r, HepVector3D &b) |
void | readPar () |
std::string | getPath () |
Private Attributes | |
double | bipx0 [npar] |
double | bipy0 [npar] |
double | bipz0 [npar] |
double | bipx1 [npar] |
double | bipy1 [npar] |
double | bipz1 [npar] |
double | bipx2 [npar] |
double | bipy2 [npar] |
double | bipz2 [npar] |
double | bipx3 [npar] |
double | bipy3 [npar] |
double | bipz3 [npar] |
double | bipx4 [npar] |
double | bipy4 [npar] |
double | bipz4 [npar] |
double | bipx5 [npar] |
double | bipy5 [npar] |
double | bipz5 [npar] |
double | bipx6 [npar] |
double | bipy6 [npar] |
double | bipz6 [npar] |
double | bipx7 [npar] |
double | bipy7 [npar] |
double | bipz7 [npar] |
double | bipx8 [npar] |
double | bipy8 [npar] |
double | bipz8 [npar] |
double | bapx0 [npar] |
double | bapy0 [npar] |
double | bapz0 [npar] |
double | bapx1 [npar] |
double | bapy1 [npar] |
double | bapz1 [npar] |
double | bapx2 [npar] |
double | bapy2 [npar] |
double | bapz2 [npar] |
double | bapx3 [npar] |
double | bapy3 [npar] |
double | bapz3 [npar] |
double | bapx4 [npar] |
double | bapy4 [npar] |
double | bapz4 [npar] |
double | bapx5 [npar] |
double | bapy5 [npar] |
double | bapz5 [npar] |
double | bapx6 [npar] |
double | bapy6 [npar] |
double | bapz6 [npar] |
double | bapx7 [npar] |
double | bapy7 [npar] |
double | bapz7 [npar] |
double | aipx0 [npar1] |
double | aipy0 [npar1] |
double | aipz0 [npar1] |
double | aipx1 [npar1] |
double | aipy1 [npar1] |
double | aipz1 [npar1] |
double | aipx2 [npar1] |
double | aipy2 [npar1] |
double | aipz2 [npar1] |
double | aipx3 [npar1] |
double | aipy3 [npar1] |
double | aipz3 [npar1] |
double | aipx4 [npar1] |
double | aipy4 [npar1] |
double | aipz4 [npar1] |
double | aipx5 [npar1] |
double | aipy5 [npar1] |
double | aipz5 [npar1] |
double | aipx6 [npar1] |
double | aipy6 [npar1] |
double | aipz6 [npar1] |
double | aipx7 [npar1] |
double | aipy7 [npar1] |
double | aipz7 [npar1] |
double | aipx8 [npar1] |
double | aipy8 [npar1] |
double | aipz8 [npar1] |
double | aapx0 [npar1] |
double | aapy0 [npar1] |
double | aapz0 [npar1] |
double | aapx1 [npar1] |
double | aapy1 [npar1] |
double | aapz1 [npar1] |
double | aapx2 [npar1] |
double | aapy2 [npar1] |
double | aapz2 [npar1] |
double | aapx3 [npar1] |
double | aapy3 [npar1] |
double | aapz3 [npar1] |
double | aapx4 [npar1] |
double | aapy4 [npar1] |
double | aapz4 [npar1] |
double | aapx5 [npar1] |
double | aapy5 [npar1] |
double | aapz5 [npar1] |
double | aapx6 [npar1] |
double | aapy6 [npar1] |
double | aapz6 [npar1] |
double | aapx7 [npar1] |
double | aapy7 [npar1] |
double | aapz7 [npar1] |
double | bp [npar] |
double | ep [npar1] |
std::string | filename |
std::string | path |
Static Private Attributes | |
static const int | npar = 18 |
static const int | npar1 = 17 |
Definition at line 15 of file MucMagneticField.h.
MucMagneticField::MucMagneticField | ( | ) |
Definition at line 10 of file MucMagneticField.cxx.
References path, readPar(), and deljobs::string.
00011 { 00012 path = std::string(getenv( "MAGNETICFIELDROOT" )); 00013 #else 00014 MucMagneticField::MucMagneticField(const std::string Path) 00015 { 00016 path = Path; 00017 #endif 00018 readPar(); 00019 }
MucMagneticField::~MucMagneticField | ( | ) |
void MucMagneticField::getMucField | ( | int | part, | |
int | layer, | |||
int | mat, | |||
HepPoint3D & | r, | |||
HepVector3D & | b | |||
) |
Definition at line 23 of file MucMagneticField.cxx.
References aapx0, aapx1, aapx2, aapx3, aapx4, aapx5, aapx6, aapx7, aapy0, aapy1, aapy2, aapy3, aapy4, aapy5, aapy6, aapy7, aapz0, aapz1, aapz2, aapz3, aapz4, aapz5, aapz6, aapz7, aipx0, aipx1, aipx2, aipx3, aipx4, aipx5, aipx6, aipx7, aipx8, aipy0, aipy1, aipy2, aipy3, aipy4, aipy5, aipy6, aipy7, aipy8, aipz0, aipz1, aipz2, aipz3, aipz4, aipz5, aipz6, aipz7, aipz8, bapx0, bapx1, bapx2, bapx3, bapx4, bapx5, bapx6, bapx7, bapy0, bapy1, bapy2, bapy3, bapy4, bapy5, bapy6, bapy7, bapz0, bapz1, bapz2, bapz3, bapz4, bapz5, bapz6, bapz7, bipx0, bipx1, bipx2, bipx3, bipx4, bipx5, bipx6, bipx7, bipx8, bipy0, bipy1, bipy2, bipy3, bipy4, bipy5, bipy6, bipy7, bipy8, bipz0, bipz1, bipz2, bipz3, bipz4, bipz5, bipz6, bipz7, bipz8, bp, ep, genRecEmupikp::i, and x.
Referenced by MagneticFieldSvc::fieldVector().
00025 { 00026 double x,y,z,bmf,emf; 00027 x = r.x()/m; 00028 y = r.y()/m; 00029 z = r.z()/m; 00030 //barrel, layer is the number of iron layer, mat=0 iron, mat=1 air 00031 if(part==1) 00032 { 00033 if(layer==0&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx0[i]; } 00034 if(layer==0&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx0[i]; } 00035 if(layer==1&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx1[i]; } 00036 if(layer==1&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx1[i]; } 00037 if(layer==2&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx2[i]; } 00038 if(layer==2&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx2[i]; } 00039 if(layer==3&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx3[i]; } 00040 if(layer==3&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx3[i]; } 00041 if(layer==4&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx4[i]; } 00042 if(layer==4&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx4[i]; } 00043 if(layer==5&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx5[i]; } 00044 if(layer==5&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx5[i]; } 00045 if(layer==6&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx6[i]; } 00046 if(layer==6&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx6[i]; } 00047 if(layer==7&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx7[i]; } 00048 if(layer==7&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapx7[i]; } 00049 if(layer==8&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipx8[i]; } 00050 bmf = bp[0]+bp[1]*x+bp[2]*x*x+bp[3]*z+bp[4]*z*z+bp[5]*z*z*z+bp[6]*z*z*z*z 00051 +bp[7]*z*z*z*z*z+bp[8]*z*z*z*z*z*z+bp[9]*z*z*z*z*z*z*z+bp[10]*y+bp[11]*y*y 00052 +bp[12]*y*y*y+bp[13]*x*z+bp[14]*z*y+bp[15]*x*y+bp[16]*z*y*y+bp[17]*z*z*y; 00053 b[0] = bmf*tesla; 00054 00055 if(layer==0&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy0[i]; } 00056 if(layer==0&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy0[i]; } 00057 if(layer==1&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy1[i]; } 00058 if(layer==1&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy1[i]; } 00059 if(layer==2&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy2[i]; } 00060 if(layer==2&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy2[i]; } 00061 if(layer==3&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy3[i]; } 00062 if(layer==3&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy3[i]; } 00063 if(layer==4&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy4[i]; } 00064 if(layer==4&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy4[i]; } 00065 if(layer==5&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy5[i]; } 00066 if(layer==5&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy5[i]; } 00067 if(layer==6&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy6[i]; } 00068 if(layer==6&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy6[i]; } 00069 if(layer==7&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy7[i]; } 00070 if(layer==7&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapy7[i]; } 00071 if(layer==8&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipy8[i]; } 00072 bmf = bp[0]+bp[1]*x+bp[2]*x*x+bp[3]*z+bp[4]*z*z+bp[5]*z*z*z+bp[6]*z*z*z*z 00073 +bp[7]*z*z*z*z*z+bp[8]*z*z*z*z*z*z+bp[9]*z*z*z*z*z*z*z+bp[10]*y+bp[11]*y*y 00074 +bp[12]*y*y*y+bp[13]*x*z+bp[14]*z*y+bp[15]*x*y+bp[16]*z*y*y+bp[17]*z*z*y; 00075 b[1] = bmf*tesla; 00076 00077 if(layer==0&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz0[i]; } 00078 if(layer==0&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz0[i]; } 00079 if(layer==1&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz1[i]; } 00080 if(layer==1&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz1[i]; } 00081 if(layer==2&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz2[i]; } 00082 if(layer==2&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz2[i]; } 00083 if(layer==3&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz3[i]; } 00084 if(layer==3&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz3[i]; } 00085 if(layer==4&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz4[i]; } 00086 if(layer==4&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz4[i]; } 00087 if(layer==5&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz5[i]; } 00088 if(layer==5&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz5[i]; } 00089 if(layer==6&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz6[i]; } 00090 if(layer==6&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz6[i]; } 00091 if(layer==7&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz7[i]; } 00092 if(layer==7&&mat==1) { for(int i=0; i<18; i++) bp[i] = bapz7[i]; } 00093 if(layer==8&&mat==0) { for(int i=0; i<18; i++) bp[i] = bipz8[i]; } 00094 bmf = bp[0]+bp[1]*x+bp[2]*x*x+bp[3]*z+bp[4]*z*z+bp[5]*z*z*z+bp[6]*z*z*z*z 00095 +bp[7]*z*z*z*z*z+bp[8]*z*z*z*z*z*z+bp[9]*z*z*z*z*z*z*z+bp[10]*y+bp[11]*y*y 00096 +bp[12]*y*y*y+bp[13]*x*z+bp[14]*z*y+bp[15]*x*y+bp[16]*z*y*y+bp[17]*z*z*y; 00097 b[2] = bmf*tesla; 00098 } 00099 //endcaps, layer is the number of iron layer, total 9. m=0 iron, m=1 air 00100 if(part==0) 00101 { 00102 if(layer==0&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx0[i]; } 00103 if(layer==0&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx0[i]; } 00104 if(layer==1&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx1[i]; } 00105 if(layer==1&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx1[i]; } 00106 if(layer==2&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx2[i]; } 00107 if(layer==2&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx2[i]; } 00108 if(layer==3&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx3[i]; } 00109 if(layer==3&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx3[i]; } 00110 if(layer==4&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx4[i]; } 00111 if(layer==4&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx4[i]; } 00112 if(layer==5&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx5[i]; } 00113 if(layer==5&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx5[i]; } 00114 if(layer==6&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx6[i]; } 00115 if(layer==6&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx6[i]; } 00116 if(layer==7&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx7[i]; } 00117 if(layer==7&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapx7[i]; } 00118 if(layer==8&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipx8[i]; } 00119 emf = ep[0]+ep[1]*x+ep[2]*x*x+ep[3]*x*x*x+ep[4]*x*x*x*x+ep[5]*x*x*x*x*x+ep[6]*x*x*x*x*x*x 00120 +ep[7]*x*x*x*x*x*x*x+ep[8]*z+ep[9]*y+ep[10]*y*y+ep[11]*y*y*y+ep[12]*x*z+ep[13]*z*y 00121 +ep[14]*y*x+ep[15]*y*x*x+ep[16]*y*y*x; 00122 b[0] = emf*tesla; 00123 00124 if(layer==0&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy0[i]; } 00125 if(layer==0&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy0[i]; } 00126 if(layer==1&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy1[i]; } 00127 if(layer==1&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy1[i]; } 00128 if(layer==2&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy2[i]; } 00129 if(layer==2&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy2[i]; } 00130 if(layer==3&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy3[i]; } 00131 if(layer==3&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy3[i]; } 00132 if(layer==4&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy4[i]; } 00133 if(layer==4&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy4[i]; } 00134 if(layer==5&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy5[i]; } 00135 if(layer==5&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy5[i]; } 00136 if(layer==6&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy6[i]; } 00137 if(layer==6&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy6[i]; } 00138 if(layer==7&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy7[i]; } 00139 if(layer==7&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapy7[i]; } 00140 if(layer==8&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipy8[i]; } 00141 emf = ep[0]+ep[1]*x+ep[2]*x*x+ep[3]*x*x*x+ep[4]*x*x*x*x+ep[5]*x*x*x*x*x+ep[6]*x*x*x*x*x*x 00142 +ep[7]*x*x*x*x*x*x*x+ep[8]*z+ep[9]*y+ep[10]*y*y+ep[11]*y*y*y+ep[12]*x*z+ep[13]*z*y 00143 +ep[14]*y*x+ep[15]*y*x*x+ep[16]*y*y*x; 00144 b[1] = emf*tesla; 00145 00146 if(layer==0&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz0[i]; } 00147 if(layer==0&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz0[i]; } 00148 if(layer==1&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz1[i]; } 00149 if(layer==1&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz1[i]; } 00150 if(layer==2&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz2[i]; } 00151 if(layer==2&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz2[i]; } 00152 if(layer==3&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz3[i]; } 00153 if(layer==3&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz3[i]; } 00154 if(layer==4&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz4[i]; } 00155 if(layer==4&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz4[i]; } 00156 if(layer==5&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz5[i]; } 00157 if(layer==5&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz5[i]; } 00158 if(layer==6&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz6[i]; } 00159 if(layer==6&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz6[i]; } 00160 if(layer==7&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz7[i]; } 00161 if(layer==7&&mat==1) { for(int i=0; i<17; i++) ep[i] = aapz7[i]; } 00162 if(layer==8&&mat==0) { for(int i=0; i<17; i++) ep[i] = aipz8[i]; } 00163 emf = ep[0]+ep[1]*x+ep[2]*x*x+ep[3]*x*x*x+ep[4]*x*x*x*x+ep[5]*x*x*x*x*x+ep[6]*x*x*x*x*x*x 00164 +ep[7]*x*x*x*x*x*x*x+ep[8]*z+ep[9]*y+ep[10]*y*y+ep[11]*y*y*y+ep[12]*x*z+ep[13]*z*y 00165 +ep[14]*y*x+ep[15]*y*x*x+ep[16]*y*y*x; 00166 b[2] = emf*tesla; 00167 } 00168 }
std::string MucMagneticField::getPath | ( | ) | [inline] |
void MucMagneticField::readPar | ( | ) |
Definition at line 169 of file MucMagneticField.cxx.
References aapx0, aapx1, aapx2, aapx3, aapx4, aapx5, aapx6, aapx7, aapy0, aapy1, aapy2, aapy3, aapy4, aapy5, aapy6, aapy7, aapz0, aapz1, aapz2, aapz3, aapz4, aapz5, aapz6, aapz7, aipx0, aipx1, aipx2, aipx3, aipx4, aipx5, aipx6, aipx7, aipx8, aipy0, aipy1, aipy2, aipy3, aipy4, aipy5, aipy6, aipy7, aipy8, aipz0, aipz1, aipz2, aipz3, aipz4, aipz5, aipz6, aipz7, aipz8, bapx0, bapx1, bapx2, bapx3, bapx4, bapx5, bapx6, bapx7, bapy0, bapy1, bapy2, bapy3, bapy4, bapy5, bapy6, bapy7, bapz0, bapz1, bapz2, bapz3, bapz4, bapz5, bapz6, bapz7, bipx0, bipx1, bipx2, bipx3, bipx4, bipx5, bipx6, bipx7, bipx8, bipy0, bipy1, bipy2, bipy3, bipy4, bipy5, bipy6, bipy7, bipy8, bipz0, bipz1, bipz2, bipz3, bipz4, bipz5, bipz6, bipz7, bipz8, filename, genRecEmupikp::i, ganga-rec::j, path, and deljobs::string.
Referenced by MucMagneticField().
00170 { 00171 filename = path + 00172 std::string( "/dat/barIronpar.txt"); 00173 ifstream infile1(filename.c_str(),ios_base::in); 00174 if(infile1) cout<<"*** Read field map: "<<filename<<endl; 00175 for(int i=0; i<18; i++) 00176 { 00177 infile1>>bipx0[i]>>bipy0[i]>>bipz0[i]>>bipx1[i]>>bipy1[i]>>bipz1[i] 00178 >>bipx2[i]>>bipy2[i]>>bipz2[i]>>bipx3[i]>>bipy3[i]>>bipz3[i] 00179 >>bipx4[i]>>bipy4[i]>>bipz4[i]>>bipx5[i]>>bipy5[i]>>bipz5[i] 00180 >>bipx6[i]>>bipy6[i]>>bipz6[i]>>bipx7[i]>>bipy7[i]>>bipz7[i] 00181 >>bipx8[i]>>bipy8[i]>>bipz8[i]; 00182 } 00183 00184 filename = path + 00185 std::string( "/dat/barAirpar.txt"); 00186 00187 ifstream infile2(filename.c_str(),ios_base::in); 00188 if(infile2) cout<<"*** Read field map: "<<filename<<endl; 00189 for(int j=0; j<18; j++) 00190 { 00191 infile2>>bapx0[j]>>bapy0[j]>>bapz0[j]>>bapx1[j]>>bapy1[j]>>bapz1[j] 00192 >>bapx2[j]>>bapy2[j]>>bapz2[j]>>bapx3[j]>>bapy3[j]>>bapz3[j] 00193 >>bapx4[j]>>bapy4[j]>>bapz4[j]>>bapx5[j]>>bapy5[j]>>bapz5[j] 00194 >>bapx6[j]>>bapy6[j]>>bapz6[j]>>bapx7[j]>>bapy7[j]>>bapz7[j]; 00195 } 00196 00197 filename = path + 00198 std::string( "/dat/endIronpar.txt"); 00199 00200 ifstream infile3(filename.c_str(),ios_base::in); 00201 if(infile3) cout<<"*** Read field map: "<<filename<<endl; 00202 for(int i=0; i<17; i++) 00203 { 00204 infile3>>aipx0[i]>>aipy0[i]>>aipz0[i]>>aipx1[i]>>aipy1[i]>>aipz1[i] 00205 >>aipx2[i]>>aipy2[i]>>aipz2[i]>>aipx3[i]>>aipy3[i]>>aipz3[i] 00206 >>aipx4[i]>>aipy4[i]>>aipz4[i]>>aipx5[i]>>aipy5[i]>>aipz5[i] 00207 >>aipx6[i]>>aipy6[i]>>aipz6[i]>>aipx7[i]>>aipy7[i]>>aipz7[i] 00208 >>aipx8[i]>>aipy8[i]>>aipz8[i]; 00209 } 00210 00211 filename = path + 00212 std::string( "/dat/endAirpar.txt"); 00213 00214 ifstream infile4(filename.c_str(),ios_base::in); 00215 if(infile4) cout<<"*** Read field map: "<<filename<<endl; 00216 for(int j=0; j<17; j++) 00217 { 00218 infile4>>aapx0[j]>>aapy0[j]>>aapz0[j]>>aapx1[j]>>aapy1[j]>>aapz1[j] 00219 >>aapx2[j]>>aapy2[j]>>aapz2[j]>>aapx3[j]>>aapy3[j]>>aapz3[j] 00220 >>aapx4[j]>>aapy4[j]>>aapz4[j]>>aapx5[j]>>aapy5[j]>>aapz5[j] 00221 >>aapx6[j]>>aapy6[j]>>aapz6[j]>>aapx7[j]>>aapy7[j]>>aapz7[j]; 00222 } 00223 }
double MucMagneticField::aapx0[npar1] [private] |
double MucMagneticField::aapx1[npar1] [private] |
double MucMagneticField::aapx2[npar1] [private] |
double MucMagneticField::aapx3[npar1] [private] |
double MucMagneticField::aapx4[npar1] [private] |
double MucMagneticField::aapx5[npar1] [private] |
double MucMagneticField::aapx6[npar1] [private] |
double MucMagneticField::aapx7[npar1] [private] |
double MucMagneticField::aapy0[npar1] [private] |
double MucMagneticField::aapy1[npar1] [private] |
double MucMagneticField::aapy2[npar1] [private] |
double MucMagneticField::aapy3[npar1] [private] |
double MucMagneticField::aapy4[npar1] [private] |
double MucMagneticField::aapy5[npar1] [private] |
double MucMagneticField::aapy6[npar1] [private] |
double MucMagneticField::aapy7[npar1] [private] |
double MucMagneticField::aapz0[npar1] [private] |
double MucMagneticField::aapz1[npar1] [private] |
double MucMagneticField::aapz2[npar1] [private] |
double MucMagneticField::aapz3[npar1] [private] |
double MucMagneticField::aapz4[npar1] [private] |
double MucMagneticField::aapz5[npar1] [private] |
double MucMagneticField::aapz6[npar1] [private] |
double MucMagneticField::aapz7[npar1] [private] |
double MucMagneticField::aipx0[npar1] [private] |
double MucMagneticField::aipx1[npar1] [private] |
double MucMagneticField::aipx2[npar1] [private] |
double MucMagneticField::aipx3[npar1] [private] |
double MucMagneticField::aipx4[npar1] [private] |
double MucMagneticField::aipx5[npar1] [private] |
double MucMagneticField::aipx6[npar1] [private] |
double MucMagneticField::aipx7[npar1] [private] |
double MucMagneticField::aipx8[npar1] [private] |
double MucMagneticField::aipy0[npar1] [private] |
double MucMagneticField::aipy1[npar1] [private] |
double MucMagneticField::aipy2[npar1] [private] |
double MucMagneticField::aipy3[npar1] [private] |
double MucMagneticField::aipy4[npar1] [private] |
double MucMagneticField::aipy5[npar1] [private] |
double MucMagneticField::aipy6[npar1] [private] |
double MucMagneticField::aipy7[npar1] [private] |
double MucMagneticField::aipy8[npar1] [private] |
double MucMagneticField::aipz0[npar1] [private] |
double MucMagneticField::aipz1[npar1] [private] |
double MucMagneticField::aipz2[npar1] [private] |
double MucMagneticField::aipz3[npar1] [private] |
double MucMagneticField::aipz4[npar1] [private] |
double MucMagneticField::aipz5[npar1] [private] |
double MucMagneticField::aipz6[npar1] [private] |
double MucMagneticField::aipz7[npar1] [private] |
double MucMagneticField::aipz8[npar1] [private] |
double MucMagneticField::bapx0[npar] [private] |
double MucMagneticField::bapx1[npar] [private] |
double MucMagneticField::bapx2[npar] [private] |
double MucMagneticField::bapx3[npar] [private] |
double MucMagneticField::bapx4[npar] [private] |
double MucMagneticField::bapx5[npar] [private] |
double MucMagneticField::bapx6[npar] [private] |
double MucMagneticField::bapx7[npar] [private] |
double MucMagneticField::bapy0[npar] [private] |
double MucMagneticField::bapy1[npar] [private] |
double MucMagneticField::bapy2[npar] [private] |
double MucMagneticField::bapy3[npar] [private] |
double MucMagneticField::bapy4[npar] [private] |
double MucMagneticField::bapy5[npar] [private] |
double MucMagneticField::bapy6[npar] [private] |
double MucMagneticField::bapy7[npar] [private] |
double MucMagneticField::bapz0[npar] [private] |
double MucMagneticField::bapz1[npar] [private] |
double MucMagneticField::bapz2[npar] [private] |
double MucMagneticField::bapz3[npar] [private] |
double MucMagneticField::bapz4[npar] [private] |
double MucMagneticField::bapz5[npar] [private] |
double MucMagneticField::bapz6[npar] [private] |
double MucMagneticField::bapz7[npar] [private] |
double MucMagneticField::bipx0[npar] [private] |
double MucMagneticField::bipx1[npar] [private] |
double MucMagneticField::bipx2[npar] [private] |
double MucMagneticField::bipx3[npar] [private] |
double MucMagneticField::bipx4[npar] [private] |
double MucMagneticField::bipx5[npar] [private] |
double MucMagneticField::bipx6[npar] [private] |
double MucMagneticField::bipx7[npar] [private] |
double MucMagneticField::bipx8[npar] [private] |
double MucMagneticField::bipy0[npar] [private] |
double MucMagneticField::bipy1[npar] [private] |
double MucMagneticField::bipy2[npar] [private] |
double MucMagneticField::bipy3[npar] [private] |
double MucMagneticField::bipy4[npar] [private] |
double MucMagneticField::bipy5[npar] [private] |
double MucMagneticField::bipy6[npar] [private] |
double MucMagneticField::bipy7[npar] [private] |
double MucMagneticField::bipy8[npar] [private] |
double MucMagneticField::bipz0[npar] [private] |
double MucMagneticField::bipz1[npar] [private] |
double MucMagneticField::bipz2[npar] [private] |
double MucMagneticField::bipz3[npar] [private] |
double MucMagneticField::bipz4[npar] [private] |
double MucMagneticField::bipz5[npar] [private] |
double MucMagneticField::bipz6[npar] [private] |
double MucMagneticField::bipz7[npar] [private] |
double MucMagneticField::bipz8[npar] [private] |
double MucMagneticField::bp[npar] [private] |
double MucMagneticField::ep[npar1] [private] |
std::string MucMagneticField::filename [private] |
const int MucMagneticField::npar = 18 [static, private] |
Definition at line 31 of file MucMagneticField.h.
const int MucMagneticField::npar1 = 17 [static, private] |
Definition at line 51 of file MucMagneticField.h.
std::string MucMagneticField::path [private] |
Definition at line 75 of file MucMagneticField.h.
Referenced by getPath(), MucMagneticField(), and readPar().