FairRoot/PandaRoot
Functions
PndCAGBHit.cxx File Reference
#include "PndCAGBHit.h"
#include "PndCAParameters.h"
#include "PndCADef.h"

Go to the source code of this file.

Functions

std::ostream & operator<< (std::ostream &out, const PndCAGBHit &a)
 
std::istream & operator>> (std::istream &in, PndCAGBHit &a)
 

Function Documentation

std::ostream& operator<< ( std::ostream &  out,
const PndCAGBHit &  a 
)

Definition at line 36 of file PndCAGBHit.cxx.

References PndCAGBHit::fC, PndCAGBHit::fErr2R, PndCAGBHit::fHalfLength, PndCAGBHit::fID, PndCAGBHit::fIRow, PndCAGBHit::fPhi, PndCAGBHit::fR, PndCAGBHit::fTubeR, PndCAGBHit::fX, PndCAGBHit::fY, PndCAGBHit::fZ, and out.

37 {
38  double tmp=0;
39  out << a.fX << " " << a.fY << " " << a.fZ << " " << a.fR << std::endl;
40  out << a.fC[0][0] << " " << a.fC[0][1] << " " << a.fC[0][2] << std::endl;// in global c.s.!
41  out << a.fC[1][0] << " " << a.fC[1][1] << " " << a.fC[1][2] << std::endl;
42  out << a.fC[2][0] << " " << a.fC[2][1] << " " << a.fC[2][2] << std::endl;
43  out << a.fErr2R << std::endl;
44  out << a.fTubeR << " " << a.fHalfLength << std::endl;
45  out << tmp << " "<< tmp <<" " << tmp <<std::endl
46  //out << a.fEX << " "<< a.fEY <<" " << a.fEZ <<std::endl
47  << a.fIRow << " " << a.fID << " " << -a.fPhi << std::endl;
48  return out;
49 }
PndCAGBHit::fX
double fX
Definition: PndCAGBHit.h:137
PndCAGBHit::fC
double fC[3][3]
Definition: PndCAGBHit.h:146
PndCAGBHit::fY
double fY
Definition: PndCAGBHit.h:138
PndCAGBHit::fTubeR
double fTubeR
Definition: PndCAGBHit.h:159
PndCAGBHit::fPhi
double fPhi
Definition: PndCAGBHit.h:151
PndCAGBHit::fErr2R
double fErr2R
Definition: PndCAGBHit.h:153
out
TFile * out
Definition: reco_muo.C:20
PndCAGBHit::fZ
double fZ
Definition: PndCAGBHit.h:139
PndCAGBHit::fHalfLength
double fHalfLength
Definition: PndCAGBHit.h:159
PndCAGBHit::fR
double fR
Definition: PndCAGBHit.h:153
std::istream& operator>> ( std::istream &  in,
PndCAGBHit &  a 
)

Definition at line 51 of file PndCAGBHit.cxx.

References c, cos(), PndCAGBHit::fC, PndCAGBHit::fErr2R, PndCAGBHit::fErr2X0, PndCAGBHit::fErr2X1, PndCAGBHit::fErr2X2, PndCAGBHit::fErrX12, PndCAGBHit::fHalfLength, PndCAGBHit::fID, PndCAGBHit::fIRow, PndCAGBHit::fPhi, PndCAGBHit::fR, PndCAGBHit::fTubeR, PndCAGBHit::fX, PndCAGBHit::fXW, PndCAGBHit::fY, PndCAGBHit::fYW, PndCAGBHit::fZ, PndCAGBHit::fZW, s, and sin().

52 {
53  double tmp;
54  UNUSED_PARAM1(tmp);
55  in >> a.fX >> a.fY >> a.fZ >> a.fR
56  >> a.fC[0][0] >> a.fC[0][1] >> a.fC[0][2] // in global c.s.!
57  >> a.fC[1][0] >> a.fC[1][1] >> a.fC[1][2]
58  >> a.fC[2][0] >> a.fC[2][1] >> a.fC[2][2]
59  >> a.fErr2R
60  >> a.fTubeR >> a.fHalfLength
61  >> tmp >> tmp >> tmp
62  //>> a.fEX >> a.fEY >> a.fEZ
63  >> a.fIRow >> a.fID >> a.fPhi;
64  a.fPhi = - a.fPhi; // TODO: fix input files
65  a.fXW = a.fX;
66  a.fYW = a.fY;
67  a.fZW = a.fZ;
68 
69  // cout << " T0 " << endl;
70  // cout << a.fC[0][0] << " " << a.fC[0][1] << " " << a.fC[0][2] << endl;
71  // cout << a.fC[1][0] << " " << a.fC[1][1] << " " << a.fC[1][2] << endl;
72  // cout << a.fC[2][0] << " " << a.fC[2][1] << " " << a.fC[2][2] << endl;
73 
74  const double s = sin(a.fPhi);
75  const double c = cos(a.fPhi);
76  // double T[3][3];
77  // T[0][0] = -s; T[0][1] = c; T[0][2] = 0;
78  // T[1][0] = c; T[1][1] = s; T[1][2] = 0;
79  // T[2][0] = 0; T[2][1] = 0; T[2][2] = -1;
80 
81  // double TC[3][3];
82  // TC[0][0] = T[0][0]*a.fC[0][0] + T[0][1]*a.fC[1][0] + T[0][2]*a.fC[2][0];
83  // TC[0][1] = T[0][0]*a.fC[0][1] + T[0][1]*a.fC[1][1] + T[0][2]*a.fC[2][1];
84  // TC[0][2] = T[0][0]*a.fC[0][2] + T[0][1]*a.fC[1][2] + T[0][2]*a.fC[2][2];
85  // TC[1][0] = T[1][0]*a.fC[0][0] + T[1][1]*a.fC[1][0] + T[1][2]*a.fC[2][0];
86  // TC[1][1] = T[1][0]*a.fC[0][1] + T[1][1]*a.fC[1][1] + T[1][2]*a.fC[2][1];
87  // TC[1][2] = T[1][0]*a.fC[0][2] + T[1][1]*a.fC[1][2] + T[1][2]*a.fC[2][2];
88  // TC[2][0] = T[2][0]*a.fC[0][0] + T[2][1]*a.fC[1][0] + T[2][2]*a.fC[2][0];
89  // TC[2][1] = T[2][0]*a.fC[0][1] + T[2][1]*a.fC[1][1] + T[2][2]*a.fC[2][1];
90  // TC[2][2] = T[2][0]*a.fC[0][2] + T[2][1]*a.fC[1][2] + T[2][2]*a.fC[2][2];
91 
92  // double TCT[3][3]; // local C
93  // // TCT'
94  // TCT[0][0] = TC[0][0]*T[0][0] + TC[0][1]*T[0][1] + TC[0][2]*T[0][2];
95  // TCT[0][1] = TC[0][0]*T[1][0] + TC[0][1]*T[1][1] + TC[0][2]*T[1][2];
96  // TCT[0][2] = TC[0][0]*T[2][0] + TC[0][1]*T[2][1] + TC[0][2]*T[2][2];
97  // TCT[1][0] = TC[1][0]*T[0][0] + TC[1][1]*T[0][1] + TC[1][2]*T[0][2];
98  // TCT[1][1] = TC[1][0]*T[1][0] + TC[1][1]*T[1][1] + TC[1][2]*T[1][2];
99  // TCT[1][2] = TC[1][0]*T[2][0] + TC[1][1]*T[2][1] + TC[1][2]*T[2][2];
100  // TCT[2][0] = TC[2][0]*T[0][0] + TC[2][1]*T[0][1] + TC[2][2]*T[0][2];
101  // TCT[2][1] = TC[2][0]*T[1][0] + TC[2][1]*T[1][1] + TC[2][2]*T[1][2];
102  // TCT[2][2] = TC[2][0]*T[2][0] + TC[2][1]*T[2][1] + TC[2][2]*T[2][2];
103 
104  // // cout << " T1 " << endl;
105  // // cout << TCT[0][0] << " " << TCT[0][1] << " " << TCT[0][2] << endl;
106  // // cout << TCT[1][0] << " " << TCT[1][1] << " " << TCT[1][2] << endl;
107  // // cout << TCT[2][0] << " " << TCT[2][1] << " " << TCT[2][2] << endl;
108 
109  // double x,y,z;
110  // PndCAParameters::GlobalToCALocal( a.fEX, a.fEY, a.fEZ, a.fPhi, x, y, z );
111  // // cout << " A " << x << " " << y << " " << z << endl;
112  // double T2[3][3];
113  // T2[0][0] = 1; T2[0][1] = 0; T2[0][2] = 0;
114  // T2[1][0] = 0; T2[1][1] = z; T2[1][2] = -y;
115  // T2[2][0] = 0; T2[2][1] = y; T2[2][2] = z;
116 
117  // TC[0][0] = T2[0][0]*TCT[0][0] + T2[0][1]*TCT[1][0] + T2[0][2]*TCT[2][0];
118  // TC[0][1] = T2[0][0]*TCT[0][1] + T2[0][1]*TCT[1][1] + T2[0][2]*TCT[2][1];
119  // TC[0][2] = T2[0][0]*TCT[0][2] + T2[0][1]*TCT[1][2] + T2[0][2]*TCT[2][2];
120  // TC[1][0] = T2[1][0]*TCT[0][0] + T2[1][1]*TCT[1][0] + T2[1][2]*TCT[2][0];
121  // TC[1][1] = T2[1][0]*TCT[0][1] + T2[1][1]*TCT[1][1] + T2[1][2]*TCT[2][1];
122  // TC[1][2] = T2[1][0]*TCT[0][2] + T2[1][1]*TCT[1][2] + T2[1][2]*TCT[2][2];
123  // TC[2][0] = T2[2][0]*TCT[0][0] + T2[2][1]*TCT[1][0] + T2[2][2]*TCT[2][0];
124  // TC[2][1] = T2[2][0]*TCT[0][1] + T2[2][1]*TCT[1][1] + T2[2][2]*TCT[2][1];
125  // TC[2][2] = T2[2][0]*TCT[0][2] + T2[2][1]*TCT[1][2] + T2[2][2]*TCT[2][2];
126 
127  // TCT[0][0] = TC[0][0]*T2[0][0] + TC[0][1]*T2[0][1] + TC[0][2]*T2[0][2];
128  // TCT[0][1] = TC[0][0]*T2[1][0] + TC[0][1]*T2[1][1] + TC[0][2]*T2[1][2];
129  // TCT[0][2] = TC[0][0]*T2[2][0] + TC[0][1]*T2[2][1] + TC[0][2]*T2[2][2];
130  // TCT[1][0] = TC[1][0]*T2[0][0] + TC[1][1]*T2[0][1] + TC[1][2]*T2[0][2];
131  // TCT[1][1] = TC[1][0]*T2[1][0] + TC[1][1]*T2[1][1] + TC[1][2]*T2[1][2];
132  // TCT[1][2] = TC[1][0]*T2[2][0] + TC[1][1]*T2[2][1] + TC[1][2]*T2[2][2];
133  // TCT[2][0] = TC[2][0]*T2[0][0] + TC[2][1]*T2[0][1] + TC[2][2]*T2[0][2];
134  // TCT[2][1] = TC[2][0]*T2[1][0] + TC[2][1]*T2[1][1] + TC[2][2]*T2[1][2];
135  // TCT[2][2] = TC[2][0]*T2[2][0] + TC[2][1]*T2[2][1] + TC[2][2]*T2[2][2];
136 
137  // cout << " T2 " << endl;
138  // cout << TCT[0][0] << " " << TCT[0][1] << " " << TCT[0][2] << endl;
139  // cout << TCT[1][0] << " " << TCT[1][1] << " " << TCT[1][2] << endl;
140  // cout << TCT[2][0] << " " << TCT[2][1] << " " << TCT[2][2] << endl;
141 
142  // const float s = y;
143  // const float c = z;
144 
145  // const float s = sin( a.fPhi );
146  // const float c = cos( a.fPhi );
147  // T = sA cA
148  // cA -sA
149  a.fErr2X0 = s*s*a.fC[0][0] + 2*s*c*a.fC[0][1] + c*c*a.fC[1][1];
150  //a.fErrX01 = c*s*(a.fC[1][1] - a.fC[0][0]) + (c*c - s*s)*a.fC[0][1];
151  a.fErr2X1 = c*c*a.fC[0][0] - 2*s*c*a.fC[0][1] + s*s*a.fC[1][1];
152  a.fErrX12 = -(c*a.fC[0][2] - s*a.fC[1][2]);
153 
154  a.fErr2X2 = a.fC[2][2];
155  return in;
156 }
cos
friend F32vec4 cos(const F32vec4 &a)
Definition: P4_F32vec4.h:112
PndCAGBHit::fX
double fX
Definition: PndCAGBHit.h:137
PndCAGBHit::fErrX12
float fErrX12
Definition: PndCAGBHit.h:143
sin
friend F32vec4 sin(const F32vec4 &a)
Definition: P4_F32vec4.h:111
s
TLorentzVector s
Definition: Pnd2DStar.C:50
PndCAGBHit::fC
double fC[3][3]
Definition: PndCAGBHit.h:146
PndCAGBHit::fY
double fY
Definition: PndCAGBHit.h:138
PndCAGBHit::fTubeR
double fTubeR
Definition: PndCAGBHit.h:159
PndCAGBHit::fZW
double fZW
Definition: PndCAGBHit.h:157
PndCAGBHit::fXW
double fXW
Definition: PndCAGBHit.h:155
PndCAGBHit::fPhi
double fPhi
Definition: PndCAGBHit.h:151
PndCAGBHit::fErr2R
double fErr2R
Definition: PndCAGBHit.h:153
PndCAGBHit::fErr2X0
float fErr2X0
Definition: PndCAGBHit.h:141
c
GeV c
Definition: ReadMesDecCorrRes_PiT.C:28
PndCAGBHit::fZ
double fZ
Definition: PndCAGBHit.h:139
PndCAGBHit::fYW
double fYW
Definition: PndCAGBHit.h:156
PndCAGBHit::fHalfLength
double fHalfLength
Definition: PndCAGBHit.h:159
PndCAGBHit::fR
double fR
Definition: PndCAGBHit.h:153
PndCAGBHit::fErr2X1
float fErr2X1
Definition: PndCAGBHit.h:142
PndCAGBHit::fErr2X2
float fErr2X2
Definition: PndCAGBHit.h:144
Generated on Sun Aug 18 2019 03:02:35 for FairRoot/PandaRoot by   doxygen 1.8.8