BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
MdcDedxTrk Class Reference

#include <MdcDedxTrk.h>

Public Member Functions

 MdcDedxTrk (RecMdcTrack &trk)
 MdcDedxTrk (RecMdcKalTrack &trk_kal, int pid)
 MdcDedxTrk ()
 ~MdcDedxTrk ()
void set_ExTrk (RecMdcTrack &trk)
void set_ExTrk_Kal (RecMdcKalTrack &trk_kal, int pid)
double cal_dedx (float)
double cal_dedx_bitrunc (float, int, int &)
double cal_dedx_median (float)
double cal_dedx_geometric (float)
double cal_dedx_geometric_trunc (float)
double cal_dedx_harmonic (float)
double cal_dedx_harmonic_trunc (float)
double cal_dedx_transform (int)
double cal_dedx_log (float, int)
double SpaceChargeCorrec (double, double, int, double)
void set_dEdx (int l, float dEdx_meas, int trkalg, int runflag, int vflag[3], double t0, vector< double > &DedxCurve_Parameter, vector< double > &DedxSigma_Parameter, MdcDedxCorrection *)
void set_phlist (const vector< double > &phlist)
void set_phlist_hit (const vector< double > &phlist)
void setVecDedxHits (const DedxHitRefVec &vecdedxhit)
vector< double > get_phlist (void)
vector< double > get_phlist_hit (void)
DedxHitRefVec getVecDedxHits () const
RecMdcTrack * trk_ptr () const
RecMdcKalTrack * trk_ptr_kal () const
int stat () const
int trk_id () const
int nsample () const
int quality () const
double dedx () const
int charge () const
double P () const
double theta () const
double phi () const
double Pt () const
double expect (int pid) const
double exp_sigma (int pid) const
double prob (int pid) const
double chi_dedx (int pid) const
double * pexpect ()
double * pexp_sigma ()
double * pprob ()
double * pchi_dedx ()

Detailed Description

Definition at line 21 of file MdcDedxTrk.h.

Constructor & Destructor Documentation

◆ MdcDedxTrk() [1/3]

MdcDedxTrk::MdcDedxTrk ( RecMdcTrack & trk)

Definition at line 42 of file MdcDedxTrk.cxx.

42 {
43 m_trk = 0;
44 m_trk_kal = 0;
45 m_stat = -1;
46 m_trk_id = 0;
47 m_quality = -99;
48 m_charge = 0;
49 m_P = 0;
50 m_theta = 0;
51 m_phi = 0;
52 m_pl_rp = 0;
53 m_nsample = 0;
54 m_dEdx = 0;
55 m_truncate = 0.7;
56
57 set_ExTrk( trk );
58#ifdef DEBUG
59 std::cout << "MdcDedxTrk(2) constructed!" << std::endl;
60#endif
61}
MdcDedxTrk::set_ExTrk
void set_ExTrk(RecMdcTrack &trk)
Definition MdcDedxTrk.cxx:90

◆ MdcDedxTrk() [2/3]

MdcDedxTrk::MdcDedxTrk ( RecMdcKalTrack & trk_kal,
int pid )

Definition at line 63 of file MdcDedxTrk.cxx.

63 {
64 m_trk = 0;
65 m_trk_kal = 0;
66 m_stat = -1;
67 m_trk_id = 0;
68 m_quality = -99;
69 m_charge = 0;
70 m_P = 0;
71 m_theta = 0;
72 m_phi = 0;
73 m_pl_rp = 0;
74 m_nsample = 0;
75 m_dEdx = 0;
76 m_truncate = 0.7;
77
78 set_ExTrk_Kal( trk_kal, pid );
79#ifdef DEBUG
80 std::cout << "MdcDedxTrk(2) kal constructed!" << std::endl;
81#endif
82}
MdcDedxTrk::set_ExTrk_Kal
void set_ExTrk_Kal(RecMdcKalTrack &trk_kal, int pid)
Definition MdcDedxTrk.cxx:119

◆ MdcDedxTrk() [3/3]

MdcDedxTrk::MdcDedxTrk ( )

Definition at line 14 of file MdcDedxTrk.cxx.

14 {
15 m_trk = 0;
16 m_trk_kal = 0;
17 m_stat = -1;
18 m_trk_id = 0;
19 m_quality = -99;
20 m_charge = 0;
21 m_P = 0;
22 m_theta = 0;
23 m_phi = 0;
24 m_pl_rp = 0;
25 m_nsample = 0;
26
27 m_dEdx = 0;
28 m_truncate = 1;
29 for ( unsigned a = 0; a < 5; a++ )
30 {
31 dedx_exp[a] = 0.0;
32 ex_sigma[a] = 0.0;
33 pid_prob[a] = 0.0;
34 chi_ex[a] = 999.0;
35 }
36
37#ifdef DEBUG
38 std::cout << "MdcDedxTrk(1) constructed!" << std::endl;
39#endif
40}

◆ ~MdcDedxTrk()

MdcDedxTrk::~MdcDedxTrk ( )

Definition at line 84 of file MdcDedxTrk.cxx.

84 {
85#ifdef DEBUG
86 std::cout << "MdcDedxTrk destructed!!!" << std::endl;
87#endif
88}

Member Function Documentation

◆ cal_dedx()

double MdcDedxTrk::cal_dedx ( float truncate)

Definition at line 149 of file MdcDedxTrk.cxx.

149 {
150 m_truncate = truncate;
151 sort( m_phlist.begin(), m_phlist.end() );
152 int nsampl = (int)( m_phlist.size() * truncate );
153 double qSum = 0;
154 unsigned int i = 0;
155 for ( vector<double>::iterator ql = m_phlist.begin(); ql != m_phlist.end(); ql++ )
156 {
157 i++;
158 if ( i <= nsampl ) qSum += ( *ql );
159 }
160
161 float trunc = qSum / nsampl;
162 return trunc;
163 std::cout << "MdcDedxTrk::cal_dedx()!!!" << std::endl;
164}

◆ cal_dedx_bitrunc()

double MdcDedxTrk::cal_dedx_bitrunc ( float truncate,
int alg,
int & usedhit )

Definition at line 166 of file MdcDedxTrk.cxx.

166 {
167 m_truncate = truncate;
168 sort( m_phlist.begin(), m_phlist.end() );
169 int nsampl = (int)( m_phlist.size() * truncate );
170 int smpl = (int)( m_phlist.size() * ( truncate + 0.05 ) );
171 int min_cut = (int)( m_phlist.size() * 0.05 + 0.5 );
172 double qSum = 0;
173 unsigned i = 0;
174 for ( vector<double>::iterator ql = m_phlist.begin(); ql != m_phlist.end(); ql++ )
175 {
176 i++;
177 if ( i <= smpl && i >= min_cut ) qSum += ( *ql );
178 }
179 double trunc = -99;
180 usedhit = smpl - min_cut + 1;
181 if ( usedhit > 0 ) trunc = qSum / usedhit;
182
183 return trunc;
184 std::cout << "MdcDedxTrk::cal_dedx()!!!" << std::endl;
185}

◆ cal_dedx_geometric()

double MdcDedxTrk::cal_dedx_geometric ( float truncate)

Definition at line 197 of file MdcDedxTrk.cxx.

197 {
198 sort( m_phlist.begin(), m_phlist.end() );
199
200 int nsampl = m_phlist.size();
201 double qSum = 1.0;
202 for ( vector<double>::iterator ql = m_phlist.begin(); ql != m_phlist.end(); ql++ )
203 { qSum *= ( *ql ); }
204
205 double trunc = pow( qSum, 1 / double( nsampl ) );
206 return trunc;
207 std::cout << "MdcDedxTrk::cal_dedx()!!!" << std::endl;
208}

◆ cal_dedx_geometric_trunc()

double MdcDedxTrk::cal_dedx_geometric_trunc ( float truncate)

Definition at line 210 of file MdcDedxTrk.cxx.

210 {
211 m_truncate = truncate;
212 sort( m_phlist.begin(), m_phlist.end() );
213 int nsampl = (int)( m_phlist.size() * truncate );
214 double qSum = 1.0;
215 unsigned i = 0;
216 for ( vector<double>::iterator ql = m_phlist.begin(); ql != m_phlist.end(); ql++ )
217 {
218 i++;
219 if ( i <= nsampl ) qSum *= ( *ql );
220 }
221
222 double trunc = pow( qSum, 1 / double( nsampl ) );
223 return trunc;
224
225 std::cout << "MdcDedxTrk::cal_dedx()!!!" << std::endl;
226}

◆ cal_dedx_harmonic()

double MdcDedxTrk::cal_dedx_harmonic ( float truncate)

Definition at line 228 of file MdcDedxTrk.cxx.

228 {
229 sort( m_phlist.begin(), m_phlist.end() );
230
231 int nsampl = m_phlist.size();
232 double qSum = 0;
233 for ( vector<double>::iterator ql = m_phlist.begin(); ql != m_phlist.end(); ql++ )
234 { qSum += 1 / ( *ql ); }
235
236 float trunc = nsampl / qSum;
237 return trunc;
238 std::cout << "MdcDedxTrk::cal_dedx()!!!" << std::endl;
239}

◆ cal_dedx_harmonic_trunc()

double MdcDedxTrk::cal_dedx_harmonic_trunc ( float truncate)

Definition at line 241 of file MdcDedxTrk.cxx.

241 {
242 m_truncate = truncate;
243 sort( m_phlist.begin(), m_phlist.end() );
244 int nsampl = (int)( m_phlist.size() * truncate );
245 double qSum = 0;
246 unsigned i = 0;
247 for ( vector<double>::iterator ql = m_phlist.begin(); ql != m_phlist.end(); ql++ )
248 {
249 i++;
250 if ( i <= nsampl ) qSum += 1 / ( *ql );
251 }
252
253 float trunc = nsampl / qSum;
254 return trunc;
255 std::cout << "MdcDedxTrk::cal_dedx()!!!" << std::endl;
256}

◆ cal_dedx_log()

double MdcDedxTrk::cal_dedx_log ( float truncate,
int alg )

Definition at line 272 of file MdcDedxTrk.cxx.

272 {
273 double qSum = 0;
274 unsigned i = 0;
275 for ( vector<double>::iterator ql = m_phlist.begin(); ql != m_phlist.end(); ql++ )
276 {
277 i++;
278 qSum += log( *ql );
279 }
280
281 float trunc = qSum / m_phlist.size();
282 return trunc;
283 std::cout << "MdcDedxTrk::cal_dedx_log()!!!" << std::endl;
284}

◆ cal_dedx_median()

double MdcDedxTrk::cal_dedx_median ( float truncate)

Definition at line 187 of file MdcDedxTrk.cxx.

187 {
188 sort( m_phlist.begin(), m_phlist.end() );
189
190 int nsample = m_phlist.size();
191 double median;
192 if ( fmod( double( nsample ), 2.0 ) ) { median = m_phlist[( nsample - 1 ) / 2]; }
193 else { median = 0.5 * ( m_phlist[nsample / 2] + m_phlist[nsample / 2 - 1] ); }
194 return median;
195}
MdcDedxTrk::nsample
int nsample() const
Definition MdcDedxTrk.h:60

◆ cal_dedx_transform()

double MdcDedxTrk::cal_dedx_transform ( int alg)

Definition at line 258 of file MdcDedxTrk.cxx.

258 {
259 sort( m_phlist.begin(), m_phlist.end() );
260
261 int nsampl = m_phlist.size();
262 double qSum = 0;
263 for ( vector<double>::iterator ql = m_phlist.begin(); ql != m_phlist.end(); ql++ )
264 { qSum += 1 / sqrt( *ql ); }
265
266 float trunc = 1 / qSum;
267
268 return trunc;
269 std::cout << "MdcDedxTrk::cal_dedx()!!!" << std::endl;
270}

◆ charge()

int MdcDedxTrk::charge ( ) const
inline

Definition at line 63 of file MdcDedxTrk.h.

63{ return m_charge; }

◆ chi_dedx()

double MdcDedxTrk::chi_dedx ( int pid) const
inline

Definition at line 119 of file MdcDedxTrk.h.

119 {
120 if ( pid >= 0 && pid < 5 ) return chi_ex[pid];
121 else return -1;
122}

◆ dedx()

double MdcDedxTrk::dedx ( ) const
inline

Definition at line 62 of file MdcDedxTrk.h.

62{ return m_dEdx; }

◆ exp_sigma()

double MdcDedxTrk::exp_sigma ( int pid) const
inline

Definition at line 109 of file MdcDedxTrk.h.

109 {
110 if ( pid >= 0 && pid < 5 ) return ex_sigma[pid];
111 else return -1;
112}

◆ expect()

double MdcDedxTrk::expect ( int pid) const
inline

Definition at line 104 of file MdcDedxTrk.h.

104 {
105 if ( pid >= 0 && pid < 5 ) return dedx_exp[pid];
106 else return -1;
107}

◆ get_phlist()

vector< double > MdcDedxTrk::get_phlist ( void )
inline

Definition at line 50 of file MdcDedxTrk.h.

50{ return m_phlist; }

◆ get_phlist_hit()

vector< double > MdcDedxTrk::get_phlist_hit ( void )
inline

Definition at line 51 of file MdcDedxTrk.h.

51{ return m_phlist_hit; }

◆ getVecDedxHits()

DedxHitRefVec MdcDedxTrk::getVecDedxHits ( ) const
inline

Definition at line 52 of file MdcDedxTrk.h.

52{ return m_vecdedxhit; }

◆ nsample()

int MdcDedxTrk::nsample ( ) const
inline

Definition at line 60 of file MdcDedxTrk.h.

60{ return m_nsample; }

Referenced by cal_dedx_median(), and set_dEdx().

◆ P()

double MdcDedxTrk::P ( ) const
inline

Definition at line 64 of file MdcDedxTrk.h.

64{ return m_P; }

◆ pchi_dedx()

double * MdcDedxTrk::pchi_dedx ( )
inline

Definition at line 75 of file MdcDedxTrk.h.

75{ return chi_ex; }

◆ pexp_sigma()

double * MdcDedxTrk::pexp_sigma ( )
inline

Definition at line 73 of file MdcDedxTrk.h.

73{ return ex_sigma; }

◆ pexpect()

double * MdcDedxTrk::pexpect ( )
inline

Definition at line 72 of file MdcDedxTrk.h.

72{ return dedx_exp; }

◆ phi()

double MdcDedxTrk::phi ( ) const
inline

Definition at line 66 of file MdcDedxTrk.h.

66{ return m_phi; }

◆ pprob()

double * MdcDedxTrk::pprob ( )
inline

Definition at line 74 of file MdcDedxTrk.h.

74{ return pid_prob; }

◆ prob()

double MdcDedxTrk::prob ( int pid) const
inline

Definition at line 114 of file MdcDedxTrk.h.

114 {
115 if ( pid >= 0 && pid < 5 ) return pid_prob[pid];
116 else return -1;
117}

◆ Pt()

double MdcDedxTrk::Pt ( ) const
inline

Definition at line 67 of file MdcDedxTrk.h.

67{ return m_P * std::sin( m_theta ); }

◆ quality()

int MdcDedxTrk::quality ( ) const
inline

Definition at line 61 of file MdcDedxTrk.h.

61{ return m_quality; }

◆ set_dEdx()

void MdcDedxTrk::set_dEdx ( int l,
float dEdx_meas,
int trkalg,
int runflag,
int vflag[3],
double t0,
vector< double > & DedxCurve_Parameter,
vector< double > & DedxSigma_Parameter,
MdcDedxCorrection * ex_calib )

Definition at line 286 of file MdcDedxTrk.cxx.

288 {
289#ifdef DEBUG
290 cout << "in MdcDedxTrk::set_dEdx() landau: " << landau << " dedx: " << dEdx
291 << " nsample(): " << nsample() << " ph size: " << m_phlist.size() << " m_P: " << m_P
292 << " theta: " << m_theta << " pl-rp: " << m_pl_rp << endl;
293#endif
294
295 m_dEdx = dEdx;
296 int dedxhit_use = (int)( m_phlist.size() * m_truncate );
297
298 // some old data with od methods
299 if ( runflag == 1 || runflag == 2 )
300 ex_calib->dedx_pid_exp_old( landau, runflag, (float)dEdx, (int)dedxhit_use, (float)m_P,
301 (float)m_theta, (float)t0, (float)m_pl_rp, dedx_exp, ex_sigma,
302 pid_prob, chi_ex );
303 // for 2009 psip data and after
304 else
305 ex_calib->dedx_pid_exp( vflag, (float)dEdx, trkalg, (int)dedxhit_use, (float)m_P,
306 (float)m_theta, (float)t0, (float)m_pl_rp, dedx_exp, ex_sigma,
307 pid_prob, chi_ex, DedxCurve_Parameter, DedxSigma_Parameter );
308
309#ifdef DEBUG
310 std::cout << "MdcDedxTrk::set_dEdx()!!!" << std::endl;
311#endif
312}
MdcDedxCorrection::dedx_pid_exp_old
void dedx_pid_exp_old(int landau, int runflag, float dedx, int Nhit, float mom, float theta, float t0, float lsamp, double dedx_ex[5], double ex_sigma[5], double pid_prob[5], double chi_dedx[5]) const
Definition MdcDedxCorrection.cxx:33
MdcDedxCorrection::dedx_pid_exp
void dedx_pid_exp(int vflag[3], float dedx, int trkalg, int Nhit, float mom, float theta, float t0, float lsamp, double dedx_ex[5], double ex_sigma[5], double pid_prob[5], double chi_dedx[5], std::vector< double > &DedxCurve_Parameter, std::vector< double > &DedxSigma_Parameter) const
Definition MdcDedxCorrection.cxx:143

◆ set_ExTrk()

void MdcDedxTrk::set_ExTrk ( RecMdcTrack & trk)

Definition at line 90 of file MdcDedxTrk.cxx.

90 {
91 for ( unsigned a = 0; a < 5; a++ )
92 {
93 dedx_exp[a] = 0.0;
94 ex_sigma[a] = 0.0;
95 pid_prob[a] = 0.0;
96 chi_ex[a] = 0.0;
97 }
98 m_stat = 1;
99 m_trk = &trk;
100 m_trk_id = trk.trackId();
101 m_quality = trk.stat();
102
103 m_charge = ( trk.helix( 2 ) > 0 ) ? 1 : -1;
104 float m_Pt = 1.0 / fabs( trk.helix( 2 ) );
105 m_P = m_Pt * sqrt( 1 + trk.helix( 4 ) * trk.helix( 4 ) );
106 m_theta = M_PI_2 - atan( trk.helix( 4 ) );
107 m_phi =
108 ( trk.helix( 1 ) < 3 * M_PI_2 ) ? trk.helix( 1 ) + M_PI_2 : trk.helix( 1 ) - 3 * M_PI_2;
109 m_nsample = trk.getNhits();
110 m_pl_rp = 1.5;
111 // cout<<"set_ExTrk: "<<trk.helix(0)<<" "<<trk.helix(1)<<" "<<trk.helix(2)<<"
112 // "<<trk.helix(3)<<" "<<trk.helix(4)<<endl;
113
114#ifdef DEBUG
115 std::cout << "MdcDedxTrk::set_ExTrk(&trk)!!!" << std::endl;
116#endif
117}
DstMdcTrack::helix
const HepVector helix() const
......
Definition DstMdcTrack.cxx:147

Referenced by MdcDedxTrk().

◆ set_ExTrk_Kal()

void MdcDedxTrk::set_ExTrk_Kal ( RecMdcKalTrack & trk_kal,
int pid )

Definition at line 119 of file MdcDedxTrk.cxx.

119 {
120 DstMdcKalTrack* dstTrk = &trk_kal;
121 for ( unsigned a = 0; a < 5; a++ )
122 {
123 dedx_exp[a] = 0.0;
124 ex_sigma[a] = 0.0;
125 pid_prob[a] = 0.0;
126 chi_ex[a] = 0.0;
127 }
128 if ( pid < 0 || pid > 4 ) pid = 2;
129 m_stat = 1;
130 m_trk_kal = &trk_kal;
131 m_trk_id = trk_kal.trackId();
132 m_quality = dstTrk->getStat( pid );
133
134 HepVector kal_helix = dstTrk->getFHelix( pid );
135 m_charge = ( kal_helix[2] > 0 ) ? 1 : -1;
136 float m_Pt = 1.0 / fabs( kal_helix[2] );
137 m_P = m_Pt * sqrt( 1 + kal_helix[4] * kal_helix[4] );
138 m_theta = M_PI_2 - atan( kal_helix[4] );
139 m_phi = ( kal_helix[1] < 3 * M_PI_2 ) ? kal_helix[1] + M_PI_2 : kal_helix[1] - 3 * M_PI_2;
140 m_nsample = trk_kal.getNhits( pid );
141 m_pl_rp = 1.5;
142 // cout<<"set_ExTrk_Kal: "<<kal_helix[0]<<" "<<kal_helix[1]<<" "<<kal_helix[2]<<"
143 // "<<kal_helix[3]<<" "<<kal_helix[4]<<endl;
144#ifdef DEBUG
145 std::cout << "MdcDedxTrk::set_ExTrk_Kal(&trk_kal)!!!" << std::endl;
146#endif
147}
DstMdcKalTrack::getFHelix
const HepVector & getFHelix(const int pid) const
Definition Event/DstEvent/include/DstEvent/DstMdcKalTrack.h:178

Referenced by MdcDedxTrk().

◆ set_phlist()

void MdcDedxTrk::set_phlist ( const vector< double > & phlist)
inline

Definition at line 46 of file MdcDedxTrk.h.

46{ m_phlist = phlist; }

Referenced by MdcDedxRecon::execute(), MdcDedxRecon::kaltrackrec(), and MdcDedxRecon::mdctrackrec().

◆ set_phlist_hit()

void MdcDedxTrk::set_phlist_hit ( const vector< double > & phlist)
inline

Definition at line 47 of file MdcDedxTrk.h.

47{ m_phlist_hit = phlist; }

Referenced by MdcDedxRecon::execute(), MdcDedxRecon::kaltrackrec(), and MdcDedxRecon::mdctrackrec().

◆ setVecDedxHits()

void MdcDedxTrk::setVecDedxHits ( const DedxHitRefVec & vecdedxhit)
inline

Definition at line 48 of file MdcDedxTrk.h.

48{ m_vecdedxhit = vecdedxhit; }

Referenced by MdcDedxRecon::execute(), MdcDedxRecon::kaltrackrec(), and MdcDedxRecon::mdctrackrec().

◆ SpaceChargeCorrec()

double MdcDedxTrk::SpaceChargeCorrec ( double m_theta,
double mom,
int Particle,
double dEdx )

Definition at line 315 of file MdcDedxTrk.cxx.

315 {
316 const int par_cand( 5 );
317 const float Charge_Mass[par_cand] = { 0.00051100, 0.10566, 0.13957, 0.4937, 0.93827 };
318 double beta_G;
319 double e_Par[5] = { 143.349, 1.7315, 0.192616, 2.90437, 1.08248 };
320 double Beta_Gamma[22] = { 0.373026, 0.479605, 0.586184, 0.692763, 0.799342, 782.779,
321 1565.56, 2348.34, 17.2727, 18.1245, 1.43297, 2.14946,
322 12.1803, 13.6132, 6.62515, 10.4109, 14.1967, 17.9825,
323 21.7683, 26.0274, 30.7596, 35.4919 };
324 double K_par[22] = { 4.64411e-05, 5.86544e-05, 8.05289e-05, 8.46981e-05, 8.92014e-05,
325 4.74517e-05, 4.51684e-05, 5.32732e-05, 6.12803e-05, 6.14592e-05,
326 8.08608e-05, 6.73184e-05, 5.46448e-05, 6.1377e-05, 6.57385e-05,
327 7.03053e-05, 6.61171e-05, 6.86824e-05, 6.246e-05, 7.25988e-05,
328 7.11034e-05, 6.24924e-05 };
329 double D_par[22] = { 0.0871504, 0.0956379, 0.117193, 0.118647, 0.127203, 0.0566449,
330 0.0529198, 0.0642525, 0.0764562, 0.081341, 0.0952263, 0.0987536,
331 0.0639901, 0.0845994, 0.0777062, 0.0823206, 0.0783874, 0.079537,
332 0.0815792, 0.0849875, 0.0824751, 0.0776296 };
333 double DSqr_par[22] = { 0.00759519, 0.0091466, 0.0137341, 0.0140772, 0.0161807,
334 0.00320864, 0.00280051, 0.00412839, 0.00584555, 0.00661636,
335 0.00906805, 0.00975227, 0.00409473, 0.00715706, 0.00603826,
336 0.00677668, 0.00614458, 0.00632613, 0.00665516, 0.00722288,
337 0.00680214, 0.00602635 };
338
339 beta_G = mom / Charge_Mass[Particle];
340 if ( beta_G < 0.3 ) beta_G = 0.3;
341 double bet = beta_G / TMath::Sqrt( beta_G * beta_G + 1 );
342 double fterm = TMath::Log( e_Par[2] + 1 / pow( beta_G, e_Par[4] ) );
343 double fitval =
344 e_Par[0] / pow( bet, e_Par[3] ) * ( e_Par[1] - pow( bet, e_Par[3] ) - fterm );
345 TGraphErrors* gr1 = new TGraphErrors( 22, Beta_Gamma, K_par, 0, 0 );
346 TGraphErrors* gr2 = new TGraphErrors( 22, Beta_Gamma, DSqr_par, 0, 0 );
347
348 double par[3];
349 par[0] = fitval;
350 par[1] = gr1->Eval( m_theta );
351 par[2] = gr2->Eval( m_theta );
352 Double_t y = fabs( cos( m_theta ) );
353 double electron_par[3] = { 334.032, 6.20658e-05, 0.00525673 };
354 double arg = TMath::Sqrt( y * y + par[2] );
355 // double cal_factor =par[0]*TMath::Exp(-(par[1]* par[0])/arg);
356 double cal_factor = TMath::Exp( -( par[1] * par[0] ) / arg );
357 double arg_electron = TMath::Sqrt( y * y + electron_par[2] );
358 // double electron_factor = electron_par[0]*TMath::Exp(-(electron_par[1]*
359 // electron_par[0])/arg);
360 double electron_factor = TMath::Exp( -( electron_par[1] * electron_par[0] ) / arg_electron );
361 // cout<<"cal_factor = "<<cal_factor<<" electron_factor = "<<electron_factor<<endl;
362 double dedx_cal = dEdx / ( cal_factor / electron_factor );
363 // double dedx_cal = dEdx/cal_factor;
364 // cout<<"m_theta= "<<m_theta<<" y ="<<y<<" beta_G = "<<beta_G <<" dEdx =
365 // "<<dEdx<<" cal dedx = "<<dedx_cal<<endl;
366 delete gr1;
367 delete gr2;
368 return dedx_cal;
369}
arg
double arg(const EvtComplex &c)
Definition EvtComplex.hh:195
Particle
HepMC::GenParticle Particle
Definition GenIMCselector.h:16
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187

Referenced by MdcDedxCorrection::dedx_pid_exp_old().

◆ stat()

int MdcDedxTrk::stat ( ) const
inline

Definition at line 57 of file MdcDedxTrk.h.

57{ return m_stat; }

◆ theta()

double MdcDedxTrk::theta ( ) const
inline

Definition at line 65 of file MdcDedxTrk.h.

65{ return m_theta; }

◆ trk_id()

int MdcDedxTrk::trk_id ( ) const
inline

Definition at line 58 of file MdcDedxTrk.h.

58{ return m_trk_id; }

Referenced by MdcDedxRecon::execute(), MdcDedxRecon::kaltrackrec(), and MdcDedxRecon::mdctrackrec().

◆ trk_ptr()

RecMdcTrack * MdcDedxTrk::trk_ptr ( ) const
inline

Definition at line 54 of file MdcDedxTrk.h.

54{ return m_trk; }

◆ trk_ptr_kal()

RecMdcKalTrack * MdcDedxTrk::trk_ptr_kal ( ) const
inline

Definition at line 55 of file MdcDedxTrk.h.

55{ return m_trk_kal; }
The documentation for this class was generated from the following files: