| 13 |
#define ZSPEBP 2.97 //cm |
#define ZSPEBP 2.97 //cm |
| 14 |
#define HORMG 22.57 //cm |
#define HORMG 22.57 //cm |
| 15 |
|
|
| 16 |
#ifdef PAMFIELD |
//#ifdef PAMFIELD |
| 17 |
#include "TrkParams.h" |
#include "TrkParams.h" |
| 18 |
#endif |
//#endif |
|
|
|
| 19 |
|
|
| 20 |
#if ROOT_VERSION_CODE >= 333572 |
#if ROOT_VERSION_CODE >= 333572 |
| 21 |
#include <TVirtualMagField.h> |
#include <TVirtualMagField.h> |
| 74 |
void Field(const Float_t* x, Float_t* b) { |
void Field(const Float_t* x, Float_t* b) { |
| 75 |
|
|
| 76 |
|
|
| 77 |
#ifdef PAMFIELD |
#ifdef PAMFIELD |
| 78 |
float xm[3]={x[0],x[1],x[2]-ZCAL-CALS3-ZSPEBP-HORMG}; |
float xm[3]={x[0],x[1],x[2]-ZCAL-CALS3-ZSPEBP-HORMG}; |
| 79 |
|
|
| 80 |
b[0] = TrkParams::GetBX((float *)xm)*10; |
b[0] = TrkParams::GetBX((float *)xm)*10; |
| 81 |
b[1] = TrkParams::GetBY((float *)xm)*10; |
b[1] = TrkParams::GetBY((float *)xm)*10; |
| 82 |
b[2] = TrkParams::GetBZ((float *)xm)*10; |
b[2] = TrkParams::GetBZ((float *)xm)*10; |
| 83 |
|
|
| 84 |
#else |
#else |
| 85 |
if(x) { |
if(x) { |
| 86 |
for (Int_t i=0; i<3; i++) b[i] = 0.0; |
for (Int_t i=0; i<3; i++) b[i] = 0.0; |
| 87 |
} else { |
} else { |
| 88 |
b[0]=b[1]=b[2]=0.0; |
b[0]=b[1]=b[2]=0.0; |
| 89 |
} |
} |
| 90 |
#endif |
#endif |
| 91 |
|
|
| 92 |
} |
} |
| 93 |
|
|