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

#include <Kpi.h>

Public Member Functions

 Kpi ()
 ~Kpi ()
bool Getkpimd ()
double Gettagmd ()
double Getmass_bc ()
double GetCQtm ()
double GetdelE_tag ()
Vint Gettagtrk1 ()
HepLorentzVector Gettagp1 ()
Vint GettagGam1 ()
void MTotal (double event, SmartDataPtr< EvtRecTrackCol > evtRecTrkCol, Vint iGood, Vint iGam, double Ebeam, int PID_flag, int Charge_candidate_D)

Detailed Description

Definition at line 20 of file Kpi.h.

Constructor & Destructor Documentation

◆ Kpi()

Kpi::Kpi ( )

Definition at line 28 of file Kpi.cxx.

28{}

◆ ~Kpi()

Kpi::~Kpi ( )

Definition at line 30 of file Kpi.cxx.

30{}

Member Function Documentation

◆ GetCQtm()

double Kpi::GetCQtm ( )
inline

Definition at line 29 of file Kpi.h.

29{ return cqtm; }

Referenced by Sing::Mdset().

◆ GetdelE_tag()

double Kpi::GetdelE_tag ( )
inline

Definition at line 30 of file Kpi.h.

30{ return delE_tag; }

Referenced by Sing::Mdset().

◆ Getkpimd()

bool Kpi::Getkpimd ( )
inline

Definition at line 26 of file Kpi.h.

26{ return kpimd; }

Referenced by Sing::Mdset().

◆ Getmass_bc()

double Kpi::Getmass_bc ( )
inline

Definition at line 28 of file Kpi.h.

28{ return mass_bc; }

Referenced by Sing::Mdset().

◆ GettagGam1()

Vint Kpi::GettagGam1 ( )
inline

Definition at line 33 of file Kpi.h.

33{ return iGamtag; }

Referenced by Sing::Mdset().

◆ Gettagmd()

double Kpi::Gettagmd ( )
inline

Definition at line 27 of file Kpi.h.

27{ return tagmd; }

Referenced by Sing::Mdset().

◆ Gettagp1()

HepLorentzVector Kpi::Gettagp1 ( )
inline

Definition at line 32 of file Kpi.h.

32{ return ptag; }

Referenced by Sing::Mdset().

◆ Gettagtrk1()

Vint Kpi::Gettagtrk1 ( )
inline

Definition at line 31 of file Kpi.h.

31{ return iGoodtag; }

Referenced by Sing::Mdset().

◆ MTotal()

void Kpi::MTotal ( double event,
SmartDataPtr< EvtRecTrackCol > evtRecTrkCol,
Vint iGood,
Vint iGam,
double Ebeam,
int PID_flag,
int Charge_candidate_D )

Definition at line 32 of file Kpi.cxx.

33 {
34
35 int nGood = iGood.size();
36 int nGam = iGam.size();
37
38 iGoodtag.clear();
39 iGamtag.clear();
40
41 double mass_bcgg, delE_tag_temp;
42 int m_chargetag, m_chargek, m_chargepi;
43 int ika_temp, ipi_temp, ipi2_temp, ipi3_temp;
44 HepLorentzVector kmfit1, kmfit2, kmfit3, kmfit4, pddd;
45
46 int cqtm_temp;
47 HepLorentzVector pddd_temp;
48 IDataProviderSvc* eventSvc = NULL;
49 Gaudi::svcLocator()->service( "EventDataSvc", eventSvc );
50 SmartDataPtr<EvtRecEvent> evtRecEvent( eventSvc, EventModel::EvtRec::EvtRecEvent );
51 SmartDataPtr<Event::EventHeader> eventHeader( eventSvc, "/Event/EventHeader" );
52
53 int runNo = eventHeader->runNumber();
54 int rec = eventHeader->eventNumber();
55
56 double xecm = 2 * Ebeam;
57
58 kpimd = false;
59 double tagmode = 0;
60
61 if ( ( evtRecEvent->totalCharged() < 2 ) ) { return; }
62
63 double ecms = xecm;
64 ;
65
66 ISimplePIDSvc* simple_pid;
67 Gaudi::svcLocator()->service( "SimplePIDSvc", simple_pid );
68
69 double deltaE_tem = 0.20;
70 int ncount1 = 0;
71
72 Hep3Vector xorigin( 0, 0, 0 );
73 IVertexDbSvc* vtxsvc;
74 Gaudi::svcLocator()->service( "VertexDbSvc", vtxsvc );
75 if ( vtxsvc->isVertexValid() )
76 {
77 double* dbv = vtxsvc->PrimaryVertex();
78 double* vv = vtxsvc->SigmaPrimaryVertex();
79 xorigin.setX( dbv[0] );
80 xorigin.setY( dbv[1] );
81 xorigin.setZ( dbv[2] );
82 }
83
84 double xv = xorigin.x();
85 double yv = xorigin.y();
86 double zv = xorigin.z();
87
88 HepPoint3D point0( 0., 0., 0. );
89 HepPoint3D IP( xorigin[0], xorigin[1], xorigin[2] );
90
91 HepLorentzVector ptrk1_temp, ptrk2_temp;
92 //////////////////////////////////////////////////////////////////
93 for ( int i = 0; i < evtRecEvent->totalCharged(); i++ )
94 {
95 EvtRecTrackIterator itTrk = evtRecTrkCol->begin() + i;
96
97 int ika = ( *itTrk )->trackId();
98
99 if ( !( *itTrk )->isMdcKalTrackValid() ) continue;
100 RecMdcKalTrack* mdcKalTrk1 = ( *itTrk )->mdcKalTrack();
101 RecMdcKalTrack::setPidType( RecMdcKalTrack::kaon );
102 /////////////////////////////////////////
103 m_chargek = mdcKalTrk1->charge();
104 if ( Charge_candidate_D != 0 )
105 {
106 if ( m_chargek != -Charge_candidate_D ) continue;
107 }
108 if ( Charge_candidate_D == 0 )
109 {
110 if ( abs( m_chargek ) != 1 ) continue;
111 }
112 /////////////////////////////////////////
113 HepVector a1 = mdcKalTrk1->getZHelixK();
114 HepSymMatrix Ea1 = mdcKalTrk1->getZErrorK();
115 VFHelix helixip3_1( point0, a1, Ea1 );
116 helixip3_1.pivot( IP );
117 HepVector vecipa1 = helixip3_1.a();
118
119 double dr1 = fabs( vecipa1[0] );
120 double dz1 = fabs( vecipa1[3] );
121 double costheta1 = cos( mdcKalTrk1->theta() );
122 if ( dr1 >= 1.0 ) continue;
123 if ( dz1 >= 10.0 ) continue;
124 if ( fabs( costheta1 ) >= 0.93 ) continue;
125 /////////////////////////////////////////
126 if ( PID_flag == 5 )
127 {
128 simple_pid->preparePID( *itTrk );
129 if ( simple_pid->probKaon() < 0.0 || simple_pid->probKaon() < simple_pid->probPion() )
130 continue;
131 }
132 /////////////////////////////////////////
133
134 WTrackParameter kam( xmass[3], mdcKalTrk1->getZHelixK(), mdcKalTrk1->getZErrorK() );
135
136 //
137 // select pi
138 //
139 for ( int j = 0; j < evtRecEvent->totalCharged(); j++ )
140 {
141 EvtRecTrackIterator itTrk = evtRecTrkCol->begin() + j;
142
143 int ipi = ( *itTrk )->trackId();
144 if ( ipi == ika ) continue;
145
146 if ( !( *itTrk )->isMdcKalTrackValid() ) continue;
147 RecMdcKalTrack* mdcKalTrk2 = ( *itTrk )->mdcKalTrack();
148 RecMdcKalTrack::setPidType( RecMdcKalTrack::pion );
149
150 /////////////////////////////////////////
151 m_chargepi = mdcKalTrk2->charge();
152 if ( ( m_chargek + m_chargepi ) != 0 ) continue;
153 /////////////////////////////////////////
154 HepVector a2 = mdcKalTrk2->getZHelix();
155 HepSymMatrix Ea2 = mdcKalTrk2->getZError();
156 VFHelix helixip3_2( point0, a2, Ea2 );
157 helixip3_2.pivot( IP );
158 HepVector vecipa2 = helixip3_2.a();
159
160 double dr2 = fabs( vecipa2[0] );
161 double dz2 = fabs( vecipa2[3] );
162 double costheta2 = cos( mdcKalTrk2->theta() );
163 if ( dr2 >= 1.0 ) continue;
164 if ( dz2 >= 10.0 ) continue;
165 if ( fabs( costheta2 ) >= 0.93 ) continue;
166 /////////////////////////////////////////
167 if ( PID_flag == 5 )
168 {
169 simple_pid->preparePID( *itTrk );
170 if ( simple_pid->probPion() < 0.0 || simple_pid->probPion() < simple_pid->probKaon() )
171 continue;
172 }
173 /////////////////////////////////////////
174
175 WTrackParameter pip( xmass[2], mdcKalTrk2->getZHelix(), mdcKalTrk2->getZError() );
176
177 //////////////////////////////////////////////////////////////
178 HepPoint3D vx( xorigin.x(), xorigin.y(), xorigin.z() );
179 HepSymMatrix Evx( 3, 0 );
180 double bx = 1E+6;
181 Evx[0][0] = bx * bx;
182 double by = 1E+6;
183 Evx[1][1] = by * by;
184 double bz = 1E+6;
185 Evx[2][2] = bz * bz;
186 VertexParameter vxpar;
187 vxpar.setVx( vx );
188 vxpar.setEvx( Evx );
189 //////////////////////////////////////////////////////////////
190
191 VertexFit* vtxfit = VertexFit::instance();
192 vtxfit->init();
193 vtxfit->AddTrack( 0, kam );
194 vtxfit->AddTrack( 1, pip );
195 vtxfit->AddVertex( 0, vxpar, 0, 1 );
196 if ( !vtxfit->Fit( 0 ) ) continue;
197 vtxfit->Swim( 0 );
198
199 WTrackParameter wkam = vtxfit->wtrk( 0 );
200 WTrackParameter wpip = vtxfit->wtrk( 1 );
201
202 HepVector kam_val = HepVector( 7, 0 );
203 kam_val = wkam.w();
204 HepVector pip_val = HepVector( 7, 0 );
205 pip_val = wpip.w();
206
207 HepLorentzVector P_KAM( kam_val[0], kam_val[1], kam_val[2], kam_val[3] );
208 HepLorentzVector P_PIP( pip_val[0], pip_val[1], pip_val[2], pip_val[3] );
209
210 P_KAM.boost( -0.011, 0, 0 );
211 P_PIP.boost( -0.011, 0, 0 );
212 pddd = P_KAM + P_PIP;
213
214 double pkpi = pddd.rho();
215
216 double temp1 = ( ecms / 2 ) * ( ecms / 2 ) - pkpi * pkpi;
217 if ( temp1 < 0 ) temp1 = 0;
218 double mass_bc_tem = sqrt( temp1 );
219
220 if ( mass_bc_tem < 1.82 || mass_bc_tem > 1.89 ) continue;
221
222 double delE_tag_tag = ecms / 2 - pddd.e();
223
224 if ( fabs( delE_tag_tag ) < deltaE_tem )
225 {
226 deltaE_tem = fabs( delE_tag_tag );
227 delE_tag_temp = delE_tag_tag;
228 mass_bcgg = mass_bc_tem;
229
230 pddd_temp = pddd;
231 cqtm_temp = m_chargek;
232 ika_temp = ika;
233 ipi_temp = ipi;
234 ncount1 = 1;
235 }
236 }
237 }
238
239 if ( ncount1 == 1 )
240 {
241 tagmode = 11;
242 if ( cqtm_temp < 0 ) tagmode = -11;
243 tagmd = tagmode;
244 mass_bc = mass_bcgg;
245 delE_tag = delE_tag_temp;
246 cqtm = -1.0 * cqtm_temp;
247
248 // cqtm = +1 for D0-->K-Pi+; cqtm = -1 for D0-bar-->K+Pi-;
249 // cqtm = 0 for D0-->K+K- ect. ...; cqtm = 0 means that we do not
250 // know whether it s D0 or D0-bar decays;
251
252 iGoodtag.push_back( ika_temp );
253 iGoodtag.push_back( ipi_temp );
254
255 iGamtag.push_back( 9999 );
256 iGamtag.push_back( 9999 );
257 iGamtag.push_back( 9999 );
258 iGamtag.push_back( 9999 );
259
260 ptag = pddd_temp;
261
262 kpimd = true;
263 }
264}
HepPoint3D
HepGeom::Point3D< double > HepPoint3D
Definition CosmicGenerator.h:36
runNo
int runNo
Definition DQA_TO_DB.cxx:13
EvtRecTrackIterator
EvtRecTrackCol::iterator EvtRecTrackIterator
Definition Event/EvtRecEvent/include/EvtRecEvent/EvtRecTrack.h:106
xmass
const double xmass[5]
Definition Gam4pikp.cxx:35
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187
ISimplePIDSvc::probKaon
virtual double probKaon()=0
ISimplePIDSvc::preparePID
virtual void preparePID(EvtRecTrack *track)=0
ISimplePIDSvc::probPion
virtual double probPion()=0
IVertexDbSvc::isVertexValid
virtual bool isVertexValid()=0
IVertexDbSvc::SigmaPrimaryVertex
virtual double * SigmaPrimaryVertex()=0
IVertexDbSvc::PrimaryVertex
virtual double * PrimaryVertex()=0
TrackPool::AddTrack
void AddTrack(const int number, const double mass, const RecMdcTrack *trk)
Definition TrackPool.cxx:21
VertexFit::wtrk
WTrackParameter wtrk(int n) const
Definition Analysis/VertexFit/include/VertexFit/VertexFit.h:91
VertexFit::init
void init()
Definition VertexFit.cxx:27
VertexFit::AddVertex
void AddVertex(int number, VertexParameter vpar, std::vector< int > lis)
Definition VertexFit.cxx:85
VertexFit::instance
static VertexFit * instance()
Definition VertexFit.cxx:15
VertexFit::Fit
bool Fit()
Definition VertexFit.cxx:288
ecms
const double ecms
Definition inclkstar.cxx:26

Referenced by Sing::Mdset().

The documentation for this class was generated from the following files: