BOSS 7.0.5
BESIII Offline Software System
Loading...
Searching...
No Matches
Public Member Functions | Static Public Member Functions | Static Public Attributes | List of all members
Hough3D Class Reference

#include <Hough3D.h>

Public Member Functions

 Hough3D ()
 
 Hough3D (const Hough3D &other)
 
 Hough3D (const Hough2D &hough2D, recHitCol hitCol, double bunchtime, double tanl, double z0)
 
 Hough3D (const Hough2D &hough2D, recHitCol hitCol, double bunchtime, double tanl, double z0, const vector< MdcHit * > *mdchit)
 
TrkRecoTrk * getTrk () const
 
double getCirX () const
 
double getCirY () const
 
double getCirR () const
 
double getD0 () const
 
double getPhi0 () const
 
double getOmega () const
 
double getZ0 () const
 
double getTanl () const
 
int getCharge () const
 
int getNfit () const
 
vector< MdcHit * > & getVecForClean ()
 
void setD0 (double d0)
 
void setPhi0 (double phi0)
 
void setOmega (double omega)
 
void setZ0 (double z0)
 
void setTanl (double tanl)
 
void setBunchT0 (double bunchT0)
 
void setRecHit (const recHitCol &rechit)
 
void setCharge (int charge)
 
int fit ()
 
void outerHit ()
 
void printRecHit () const
 
double getPt () const
 
double getPz () const
 
double getP () const
 
double getChi2 () const
 
void print ()
 
 Hough3D ()
 
 Hough3D (const Hough3D &other)
 
 Hough3D (const Hough2D &hough2D, recHitCol hitCol, double bunchtime, double tanl, double z0)
 
 Hough3D (const Hough2D &hough2D, recHitCol hitCol, double bunchtime, double tanl, double z0, const vector< MdcHit * > *mdchit)
 
TrkRecoTrk * getTrk () const
 
double getCirX () const
 
double getCirY () const
 
double getCirR () const
 
double getD0 () const
 
double getPhi0 () const
 
double getOmega () const
 
double getZ0 () const
 
double getTanl () const
 
int getCharge () const
 
int getNfit () const
 
vector< MdcHit * > & getVecForClean ()
 
void setD0 (double d0)
 
void setPhi0 (double phi0)
 
void setOmega (double omega)
 
void setZ0 (double z0)
 
void setTanl (double tanl)
 
void setBunchT0 (double bunchT0)
 
void setRecHit (const recHitCol &rechit)
 
void setCharge (int charge)
 
int fit ()
 
void outerHit ()
 
void printRecHit () const
 
double getPt () const
 
double getPz () const
 
double getP () const
 
double getChi2 () const
 
void print ()
 

Static Public Member Functions

static void setContext (TrkContextEv *context)
 
static void setCalib (const MdcCalibFunSvc *mdcCalibFunSvc)
 
static void setContext (TrkContextEv *context)
 
static void setCalib (const MdcCalibFunSvc *mdcCalibFunSvc)
 

Static Public Attributes

static int m_debug =0
 
static TrkContextEv * _context =NULL
 
static const MdcCalibFunSvc * _mdcCalibFunSvc =NULL
 
static int m_qualityFactor =3
 
static double m_dropTrkChi2Cut =99999
 
static double m_dropTrkDrCut =1
 
static double m_dropTrkDzCut =10
 
static double m_dropTrkNhitCut =9
 
static double m_dropTrkChi2NdfCut =99999
 

Detailed Description

Definition at line 25 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Constructor & Destructor Documentation

◆ Hough3D() [1/8]

Hough3D::Hough3D ( )

Definition at line 16 of file Hough3D.cxx.

16 {
17 //_m_gm = MdcDetector::instance(0);
18}

◆ Hough3D() [2/8]

Hough3D::Hough3D ( const Hough3D &  other)

Definition at line 19 of file Hough3D.cxx.

19 {
20 _circleR=other._circleR;
21 _circleX=other._circleX;
22 _circleY=other._circleY;
23 _charge=other._charge;
24 _d0=other._d0;
25 _phi0=other._phi0;
26 _omega=other._omega;
27 _tanl=other._tanl;
28 _z0=other._z0;
29
30 _pt=other._pt;
31 _pz=other._pz;
32 _p=other._p;
33 _nfit=other._nfit;
34 _bunchT0=other._bunchT0;
35 _chi2_aver=other._chi2_aver;
36 _m_gm = other._m_gm;
37 _recHitVec=other._recHitVec;
38 vec_mdcHit=other.vec_mdcHit;
39}
EvtCyclic3::other
Index other(Index i, Index j)
Definition: EvtCyclic3.cc:118

◆ Hough3D() [3/8]

Hough3D::Hough3D ( const Hough2D &  hough2D,
recHitCol  hitCol,
double  bunchtime,
double  tanl,
double  z0 
)

Definition at line 41 of file Hough3D.cxx.

41 :_tanl(tanl),_z0(z0){
42 //_m_gm = MdcDetector::instance(0);
43 _circleR=hough2D.getCirR();
44 _circleX=hough2D.getCirX();
45 _circleY=hough2D.getCirY();
46 _d0 = hough2D.getD0();
47 _phi0 = hough2D.getPhi0();
48 _omega= hough2D.getOmega();
49 _charge=hough2D.getCharge();
50 _bunchT0=bunchtime;
51 _recHitVec=hitCol;
52 _pt = 0;
53 _p = 0;
54 _pz = 0;
55 _nfit=0;
56 newTrack=NULL;
57}

◆ Hough3D() [4/8]

Hough3D::Hough3D ( const Hough2D &  hough2D,
recHitCol  hitCol,
double  bunchtime,
double  tanl,
double  z0,
const vector< MdcHit * > *  mdchit 
)

Definition at line 58 of file Hough3D.cxx.

58 :_tanl(tanl),_z0(z0),vec_mdcHit(mdchit){
59 //_m_gm = MdcDetector::instance(0);
60 _circleR=hough2D.getCirR();
61 _circleX=hough2D.getCirX();
62 _circleY=hough2D.getCirY();
63 _d0 = hough2D.getD0();
64 _phi0 = hough2D.getPhi0();
65 _omega= hough2D.getOmega();
66 _charge=hough2D.getCharge();
67 _bunchT0=bunchtime;
68 _recHitVec=hitCol;
69 _pt = 0;
70 _p = 0;
71 _pz = 0;
72 _nfit=0;
73 newTrack=NULL;
74}

◆ Hough3D() [5/8]

Hough3D::Hough3D ( )

◆ Hough3D() [6/8]

Hough3D::Hough3D ( const Hough3D &  other)

◆ Hough3D() [7/8]

Hough3D::Hough3D ( const Hough2D &  hough2D,
recHitCol  hitCol,
double  bunchtime,
double  tanl,
double  z0 
)

◆ Hough3D() [8/8]

Hough3D::Hough3D ( const Hough2D &  hough2D,
recHitCol  hitCol,
double  bunchtime,
double  tanl,
double  z0,
const vector< MdcHit * > *  mdchit 
)

Member Function Documentation

◆ fit() [1/2]

int Hough3D::fit ( )

Definition at line 76 of file Hough3D.cxx.

76 {
77 if(_charge==-1){
78// cout<<" charge -1 "<<endl;
79// _d0=-_d0;
80// _omega=-1.*fabs(_omega);
81 }
82 if(_charge==1){
83// cout<<" charge +1 "<<endl;
84// _d0=_d0;
85// _omega=1.*fabs(_omega);
86 //_phi0=_phi0-M_PI;
87 }
88
89 //outerHit(); // delete 4 hit in a wire for kalman
90
91 TrkExchangePar tt(_d0,_phi0,_omega,_z0,_tanl);
92 if (m_debug>0) cout<<"q d0 phi0 omega: "<<_charge<<" "<<_d0<<" "<<_phi0<<" "<<_omega<<" "<<_tanl<<" "<<_z0<<endl;
93 TrkHelixMaker helixfactory;
94 float chisum =0.;
95 newTrack = helixfactory.makeTrack(tt, chisum, *_context, _bunchT0);
96 //vectrk_for_clean.push_back(newTrack);
97 bool permitFlips = true;
98 bool lPickHits = true;
99 helixfactory.setFlipAndDrop(*newTrack, permitFlips, lPickHits);
100 TrkHitList* trkHitList = newTrack->hits();
101 int digiId=0;
102 std::vector<HoughRecHit>::iterator iter = _recHitVec.begin();
103 for (;iter != _recHitVec.end(); iter++, digiId++) {
104 if( (*iter).getflag()!=0 ) continue;
105 //if( (*iter)->calDriftDist(_bunchT0,0)>0.8) continue;
106 //if( (*iter)->driftTime(_bunchT0,0)>800) continue;
107 //cout<<"rechit dist time "<<(*iter).calDriftDist(_bunchT0,0)<<" "<<(*iter).driftTime(_bunchT0,0)<<endl;
108
109 const MdcDigi* aDigi = (*iter).digi();
110 MdcHit* hit;
111 // compare
112 vector<MdcHit*>::const_iterator iter_digi = (*vec_mdcHit).begin();
113 for (;iter_digi != (*vec_mdcHit).end(); iter_digi++) {
114 if( (*iter_digi)->digi() == (*iter).digi() ) {
115 hit = (*iter_digi);
116 }
117 }
118 Identifier id = aDigi->identify();
119 int layer = MdcID::layer(id);
120 int wire = MdcID::wire(id);
121 hit->setCalibSvc(_mdcCalibFunSvc);
122 hit->setCountPropTime(true);
123 // new fltLen and ambig
124 int newAmbig = 0;
125 // calc. position of this point
126 MdcRecoHitOnTrack temp( *hit, newAmbig, _bunchT0*1.e9);//m_bunchT0 nano second here
127 MdcHitOnTrack* newhot = &temp;
128 double flight;
129 double z_flight;
130 //if( layer<8 ) flight = sqrt( ((*iter).getsPos())*((*iter).getsPos())+((*iter).getzPos())*((*iter).getzPos())) ;
131 // if( layer<8 ) flight =(*iter).getsPos();
132 flight = (*iter).getRet().second;
133 double distoTrack= (*iter).getDisToTrack();
134 // if( layer<8 ) newhot->setFltLen( (*iter).getsPos()); //caculate by mc
135 // else newhot->setFltLen( (*iter).getRet().second); //caculate by mc
136 newhot->setFltLen( flight );
137 // if(layer==18) continue;
138
139 int use_in3d=1;
140 // if(hit->driftDist(_bunchT0,0)>1.0) use_in3d=0;
141 if(hit->driftTime(_bunchT0,0)>1000) use_in3d=0;
142 //if(m_debug>0) cout<<"flt : "<<(*iter).getsPos()<<endl;
143 if(m_debug>0) cout<<"flt : "<<flight<<endl;
144 if(m_debug>0) std::cout<<" ("<<layer<<","<<wire<<",rt "<<hit->rawTime()<<",dt "<<hit->driftTime(_bunchT0,0)<<") T0 " << _mdcCalibFunSvc->getT0(layer,wire) <<" timewalk "<< _mdcCalibFunSvc->getTimeWalk(layer, aDigi->getChargeChannel())<< "dist: "<<hit->driftDist(_bunchT0,0)<<" disToTrk "<< distoTrack<<" use?: "<<use_in3d<<endl;
145 // if(hit->driftTime(_bunchT0,0)>800) continue;
146 if(use_in3d==0) continue;
147 trkHitList->appendHot(newhot);
148 }
149 if( m_debug>0) newTrack->printAll(std::cout);
150
151 TrkHitList* qhits = newTrack->hits();
152 TrkErrCode err=qhits->fit();
153 const TrkFit* newFit = newTrack->fitResult();
154 int nActiveHit = newTrack->hots()->nActive();
155 int fit_stat=false;
156 double chi2=-999.;
157 if (!newFit || (newFit->nActive()<5)||err.failure()!=0){
158 //if (!newFit )
159 if(m_debug>0){
160 cout << " global 3d fit failed ";
161 if(newFit) cout <<" nAct "<<newFit->nActive();
162 cout<<"ERR1 failure ="<<err.failure()<<endl;
163 fit_stat=0;
164 }
165 }else{
166 TrkExchangePar par = newFit->helix(0.);
167 if( abs( 1/(par.omega()) )>0.03) {
168 fit_stat = 1;
169 chi2=newFit->chisq();
170 }
171 else {
172 fit_stat = 0;
173 chi2=-999;
174 }
175
176 bool badQuality = false;
177 if(fabs(chi2)>m_qualityFactor*m_dropTrkChi2Cut ){
178 if(m_debug>0){
179 std::cout<<__FILE__<<" "<<__LINE__<<" drop track by chi2 "<<chi2<<" > "<<m_qualityFactor<<" * "<<m_dropTrkChi2Cut <<std::endl;
180 }
181 badQuality=1;
182 }
183 if( fabs(par.d0())>m_qualityFactor*m_dropTrkDrCut) {
184 if(m_debug>0){
185 std::cout<<__FILE__<<" "<<__LINE__<<" drop track by d0 "<<par.d0()<<" > "<<m_qualityFactor<<" * "<<m_dropTrkDrCut <<std::endl;
186 }
187 badQuality=1;
188 }
189 if( fabs(par.z0())>m_qualityFactor*m_dropTrkDzCut) {
190 if(m_debug>0){
191 std::cout<<__FILE__<<" "<<__LINE__<<" drop track by z0 "<<par.z0()<<" > "<<m_qualityFactor<<" * "<<m_dropTrkDzCut <<std::endl;
192 }
193 badQuality=1;
194 }
195 if( (fabs(chi2)/qhits->nHit()) > m_qualityFactor*m_dropTrkChi2NdfCut) {
196 if(m_debug>0){
197 std::cout<<__FILE__<<" "<<__LINE__<<" drop track by chi2/ndf"<<(fabs(chi2)/qhits->nHit()) <<" > "<<m_qualityFactor<<" * "<<m_dropTrkChi2NdfCut<<std::endl;
198 }
199 badQuality=1;
200 }
201 if( nActiveHit <= m_dropTrkNhitCut) {
202 if(m_debug>0){
203 std::cout<<__FILE__<<" "<<__LINE__<<" drop track by nhit"<<nActiveHit <<" <= "<<m_dropTrkNhitCut<<std::endl;
204 }
205 badQuality=1;
206 }
207 //badQuality = false; //not delete track in this stage
208 if( badQuality) fit_stat=0;
209 }
210 if( fit_stat!=1 ){
211 delete newTrack;
212 vectrk_for_clean.pop_back();
213 }
214 if( fit_stat==1 ){
215 if(m_debug>0){
216 cout << " global 3d fit success"<<endl;
217 cout<<__FILE__<<__LINE__<<"AFTER fit 3d "<<endl;
218 newTrack->print(std::cout);
219 }
220 TrkExchangePar par = newFit->helix(0.);
221 _d0=par.d0();
222 _phi0=par.phi0();
223 _omega=par.omega();
224 _tanl=par.tanDip();
225 _z0=par.z0();
226
227 _circleR=_charge/par.omega();
228 _circleX= sin(par.phi0()) *(par.d0()+1/par.omega()) * -1.; //z axis verse,x_babar = -x_belle
229 _circleY= -1.*cos(par.phi0()) *(par.d0()+1/par.omega()) * -1;//???
230 double pxy=fabs(_circleR/333.567);
231 double pz=pxy*par.tanDip();
232 double pxyz=_charge*sqrt(pxy*pxy+pz*pz);
233 _pt=pxy;
234 _pz=pz;
235 _p=pxyz;
236 if(m_debug>0){
237 cout<<" circle after globle 3d: "<<"("<<_circleX<<","<<_circleY<<","<<_circleR<<")"<<endl;
238 cout<<"pt_rec: "<<_pt<<endl;
239 cout<<"pz_rec: "<<_pz<<endl;
240 cout<<"p_rec: "<<_p<<endl;
241 }
242
243 int nfit3d=0;
244 if(m_debug>0) cout<<" hot list:"<<endl;
245 TrkHotList::hot_iterator hotIter= qhits->hotList().begin();
246 int lay=((MdcHit*)(hotIter->hit()))->layernumber();
247 int wir=((MdcHit*)(hotIter->hit()))->wirenumber();
248 while (hotIter!=qhits->hotList().end()) {
249 if(m_debug>0){ cout <<"(" <<((MdcHit*)(hotIter->hit()))->layernumber()
250 <<","<<((MdcHit*)(hotIter->hit()))->wirenumber()
251 <<":"<<hotIter->isActive()<<") ";
252 }
253 // cout << "nuse "<<hotIter->hit()->nUsedHits()<<endl;
254 if (hotIter->isActive()==1) nfit3d++;
255 hotIter++;
256 }
257 _nfit=nfit3d;
258 }
259 _chi2_aver=chi2/_nfit;
260 // for(int i=0;i<t_hitCol.size();i++){
261 // delete t_hitCol[i];
262 // }
263 if(m_debug>0) cout<<" in 3D fit number of Active hits : "<<_nfit<<endl;
264 return fit_stat;
265}
iter
EvtStreamInputIterator< typename Generator::result_type > iter(Generator gen, int N=0)
Definition: EvtStreamInputIterator.hh:100
sin
double sin(const BesAngle a)
Definition: InstallArea/include/MdcGeom/MdcGeom/BesAngle.h:210
cos
double cos(const BesAngle a)
Definition: InstallArea/include/MdcGeom/MdcGeom/BesAngle.h:213
vectrk_for_clean
vector< TrkRecoTrk * > vectrk_for_clean
Definition: Hough3D.cxx:4
MdcCalibFunSvc::getT0
double getT0(int layid, int cellid) const
Definition: MdcCalibFunSvc.cxx:507
MdcCalibFunSvc::getTimeWalk
double getTimeWalk(int layid, double Q) const
Definition: MdcCalibFunSvc.cxx:514
MdcHit::setCalibSvc
void setCalibSvc(const MdcCalibFunSvc *calibSvc)
Definition: MdcHit.cxx:136
MdcHit::driftTime
double driftTime(double tof, double z) const
Definition: MdcHit.cxx:142
MdcHit::setCountPropTime
void setCountPropTime(const bool count)
Definition: InstallArea/include/MdcData/MdcData/MdcHit.h:86
MdcHit::driftDist
double driftDist(double, int, double, double, double) const
Definition: MdcHit.cxx:156
MdcID::layer
static int layer(const Identifier &id)
Values of different levels (failure returns 0)
Definition: MdcID.cxx:49
MdcID::wire
static int wire(const Identifier &id)
Definition: MdcID.cxx:54
RawData::identify
virtual Identifier identify() const
Definition: RawData.cxx:15
RawData::getChargeChannel
unsigned int getChargeChannel() const
Definition: RawData.cxx:45
TrkAbsFit::chisq
virtual double chisq() const =0
TrkFit::nActive
virtual int nActive() const =0
TrkFit::helix
virtual TrkExchangePar helix(double fltL) const =0
TrkFundHit::begin
TrkFundHit::hot_iterator begin() const
Definition: InstallArea/include/TrkBase/TrkBase/TrkFundHit.h:113
TrkHitList::fit
TrkErrCode fit()
Definition: TrkHitList.cxx:59
TrkHitList::appendHot
TrkHitOnTrk * appendHot(const TrkHitOnTrk *theHot)
Definition: TrkHitList.cxx:91
TrkHitList::hotList
const TrkHotList & hotList() const
Definition: TrkHitList.cxx:51
TrkHotList::nActive
virtual int nActive(TrkEnums::TrkViewInfo view=TrkEnums::bothView) const =0
TrkRecoTrk::print
virtual void print(std::ostream &) const
TrkRecoTrk::fitResult
const TrkFit * fitResult() const
Definition: TrkRecoTrk.cxx:387
TrkRecoTrk::printAll
virtual void printAll(std::ostream &) const
TrkSimpleMaker::setFlipAndDrop
bool setFlipAndDrop(TrkRecoTrk &, bool allowFlips, bool allowDrops) const
TrkSimpleMaker::makeTrack
TrkRecoTrk * makeTrack(const TrkExchangePar &helix, const double chi2, const TrkContext &, double trackT0) const

◆ fit() [2/2]

int Hough3D::fit ( )

◆ getCharge() [1/2]

int Hough3D::getCharge ( ) const
inline

Definition at line 42 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

42{return _charge;}

◆ getCharge() [2/2]

int Hough3D::getCharge ( ) const
inline

Definition at line 42 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

42{return _charge;}

◆ getChi2() [1/2]

double Hough3D::getChi2 ( ) const
inline

Definition at line 65 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

65{return _chi2_aver;}

◆ getChi2() [2/2]

double Hough3D::getChi2 ( ) const
inline

Definition at line 65 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

65{return _chi2_aver;}

◆ getCirR() [1/2]

double Hough3D::getCirR ( ) const
inline

Definition at line 36 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

36{return _circleR;}

◆ getCirR() [2/2]

double Hough3D::getCirR ( ) const
inline

Definition at line 36 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

36{return _circleR;}

◆ getCirX() [1/2]

double Hough3D::getCirX ( ) const
inline

Definition at line 34 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

34{return _circleX;}

◆ getCirX() [2/2]

double Hough3D::getCirX ( ) const
inline

Definition at line 34 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

34{return _circleX;}

◆ getCirY() [1/2]

double Hough3D::getCirY ( ) const
inline

Definition at line 35 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

35{return _circleY;}

◆ getCirY() [2/2]

double Hough3D::getCirY ( ) const
inline

Definition at line 35 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

35{return _circleY;}

◆ getD0() [1/2]

double Hough3D::getD0 ( ) const
inline

Definition at line 37 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

37{return _d0;}

◆ getD0() [2/2]

double Hough3D::getD0 ( ) const
inline

Definition at line 37 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

37{return _d0;}

◆ getNfit() [1/2]

int Hough3D::getNfit ( ) const
inline

Definition at line 43 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

43{return _nfit;}

Referenced by HoughTrack::getNfit3D().

◆ getNfit() [2/2]

int Hough3D::getNfit ( ) const
inline

Definition at line 43 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

43{return _nfit;}

◆ getOmega() [1/2]

double Hough3D::getOmega ( ) const
inline

Definition at line 39 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

39{return _omega;}

◆ getOmega() [2/2]

double Hough3D::getOmega ( ) const
inline

Definition at line 39 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

39{return _omega;}

◆ getP() [1/2]

double Hough3D::getP ( ) const
inline

Definition at line 64 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

64{return _p;}

◆ getP() [2/2]

double Hough3D::getP ( ) const
inline

Definition at line 64 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

64{return _p;}

◆ getPhi0() [1/2]

double Hough3D::getPhi0 ( ) const
inline

Definition at line 38 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

38{return _phi0;}

◆ getPhi0() [2/2]

double Hough3D::getPhi0 ( ) const
inline

Definition at line 38 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

38{return _phi0;}

◆ getPt() [1/2]

double Hough3D::getPt ( ) const
inline

Definition at line 62 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

62{return _pt;}

◆ getPt() [2/2]

double Hough3D::getPt ( ) const
inline

Definition at line 62 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

62{return _pt;}

◆ getPz() [1/2]

double Hough3D::getPz ( ) const
inline

Definition at line 63 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

63{return _pz;}

◆ getPz() [2/2]

double Hough3D::getPz ( ) const
inline

Definition at line 63 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

63{return _pz;}

◆ getTanl() [1/2]

double Hough3D::getTanl ( ) const
inline

Definition at line 41 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

41{return _tanl;}

◆ getTanl() [2/2]

double Hough3D::getTanl ( ) const
inline

Definition at line 41 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

41{return _tanl;}

◆ getTrk() [1/2]

TrkRecoTrk * Hough3D::getTrk ( ) const
inline

Definition at line 32 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

32{return newTrack;}

◆ getTrk() [2/2]

TrkRecoTrk * Hough3D::getTrk ( ) const
inline

Definition at line 32 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

32{return newTrack;}

◆ getVecForClean() [1/2]

vector< MdcHit * > & Hough3D::getVecForClean ( )
inline

Definition at line 44 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

44{return vec_for_clean;}
vec_for_clean
vector< MdcHit * > vec_for_clean
Definition: Hough3D.cxx:3

Referenced by HoughTrack::getVecForClean3D().

◆ getVecForClean() [2/2]

vector< MdcHit * > & Hough3D::getVecForClean ( )
inline

Definition at line 44 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

44{return vec_for_clean;}

◆ getZ0() [1/2]

double Hough3D::getZ0 ( ) const
inline

Definition at line 40 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

40{return _z0;}

◆ getZ0() [2/2]

double Hough3D::getZ0 ( ) const
inline

Definition at line 40 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

40{return _z0;}

◆ outerHit() [1/2]

void Hough3D::outerHit ( )

Definition at line 276 of file Hough3D.cxx.

276 {
277 for(int ihit=0;ihit<_recHitVec.size();ihit++){
278 if( _recHitVec[ihit].calDriftDist(_bunchT0,0)>0.8 || _recHitVec[ihit].driftTime(_bunchT0,0)>800 ) _recHitVec[ihit].setflag(-5); //flag -5 time or drift
279 }
280
281 vector<int> vec_layer_num[43];
282 for(int ilay=0;ilay<43;ilay++){
283 for(int ihit=0;ihit<_recHitVec.size();ihit++){
284 if( _recHitVec[ihit].getLayerId()==ilay && _recHitVec[ihit].getflag()==0 ) {
285 vec_layer_num[ilay].push_back( _recHitVec[ihit].getWireId() );
286 }
287 }
288 std::sort( vec_layer_num[ilay].begin(),vec_layer_num[ilay].end(),less_layer_3D);
289 if( vec_layer_num[ilay].size() < 4 ) continue;
290 for(int j=0;j<vec_layer_num[ilay].size();j++) {
291 // if( (vec_layer_num[ilay][j]+1 == vec_layer_num[ilay][j+1]) && (vec_layer_num[ilay][j+1]+1 == vec_layer_num[ilay][j+2]) && (vec_layer_num[ilay][j+2]+1 == vec_layer_num[ilay][j+3]) ) {
292 for(int jhit=0;jhit<_recHitVec.size();jhit++) {
293 if( (ilay==_recHitVec[jhit].getLayerId()) && (vec_layer_num[ilay][j]==_recHitVec[jhit].getWireId() ) ) {_recHitVec[jhit].setflag(-1);}
294 }
295 //}
296 }
297 }
298}
less_layer_3D
bool less_layer_3D(const int &a, const int &b)
Definition: Hough3D.cxx:272

◆ outerHit() [2/2]

void Hough3D::outerHit ( )

◆ print() [1/2]

void Hough3D::print ( )

Definition at line 267 of file Hough3D.cxx.

267 {
268 cout<<" print hough3d "<<endl;
269 cout<<"par: "<<_d0<<" "<<_phi0<<" "<<_omega<<" "<<_tanl<<" "<<_z0<<endl;
270}

◆ print() [2/2]

void Hough3D::print ( )

◆ printRecHit() [1/2]

void Hough3D::printRecHit ( ) const

◆ printRecHit() [2/2]

void Hough3D::printRecHit ( ) const

◆ setBunchT0() [1/2]

void Hough3D::setBunchT0 ( double  bunchT0)
inline

Definition at line 52 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

52{_bunchT0=bunchT0;}

◆ setBunchT0() [2/2]

void Hough3D::setBunchT0 ( double  bunchT0)
inline

Definition at line 52 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

52{_bunchT0=bunchT0;}

◆ setCalib() [1/2]

static void Hough3D::setCalib ( const MdcCalibFunSvc *  mdcCalibFunSvc)
inlinestatic

Definition at line 56 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

56{_mdcCalibFunSvc=mdcCalibFunSvc;}

Referenced by MdcHoughFinder::initialize().

◆ setCalib() [2/2]

static void Hough3D::setCalib ( const MdcCalibFunSvc *  mdcCalibFunSvc)
inlinestatic

Definition at line 56 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

56{_mdcCalibFunSvc=mdcCalibFunSvc;}

◆ setCharge() [1/2]

void Hough3D::setCharge ( int  charge)
inline

Definition at line 54 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

54{_charge=charge;}
charge
float charge
Definition: prepare_data.cxx:28

◆ setCharge() [2/2]

void Hough3D::setCharge ( int  charge)
inline

Definition at line 54 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

54{_charge=charge;}

◆ setContext() [1/2]

static void Hough3D::setContext ( TrkContextEv *  context)
inlinestatic

Definition at line 55 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

55{_context=context;}

Referenced by MdcHoughFinder::initialize().

◆ setContext() [2/2]

static void Hough3D::setContext ( TrkContextEv *  context)
inlinestatic

Definition at line 55 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

55{_context=context;}

◆ setD0() [1/2]

void Hough3D::setD0 ( double  d0)
inline

Definition at line 46 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

46{_d0=d0;}

◆ setD0() [2/2]

void Hough3D::setD0 ( double  d0)
inline

Definition at line 46 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

46{_d0=d0;}

◆ setOmega() [1/2]

void Hough3D::setOmega ( double  omega)
inline

Definition at line 48 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

48{_omega=omega;}

◆ setOmega() [2/2]

void Hough3D::setOmega ( double  omega)
inline

Definition at line 48 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

48{_omega=omega;}

◆ setPhi0() [1/2]

void Hough3D::setPhi0 ( double  phi0)
inline

Definition at line 47 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

47{_phi0=phi0;}

◆ setPhi0() [2/2]

void Hough3D::setPhi0 ( double  phi0)
inline

Definition at line 47 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

47{_phi0=phi0;}

◆ setRecHit() [1/2]

void Hough3D::setRecHit ( const recHitCol &  rechit)
inline

Definition at line 53 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

53{_recHitVec=rechit;}

◆ setRecHit() [2/2]

void Hough3D::setRecHit ( const recHitCol &  rechit)
inline

Definition at line 53 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

53{_recHitVec=rechit;}

◆ setTanl() [1/2]

void Hough3D::setTanl ( double  tanl)
inline

Definition at line 50 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

50{_tanl=tanl;}

◆ setTanl() [2/2]

void Hough3D::setTanl ( double  tanl)
inline

Definition at line 50 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

50{_tanl=tanl;}

◆ setZ0() [1/2]

void Hough3D::setZ0 ( double  z0)
inline

Definition at line 49 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

49{_z0=z0;}

◆ setZ0() [2/2]

void Hough3D::setZ0 ( double  z0)
inline

Definition at line 49 of file Reconstruction/MdcHoughFinder/MdcHoughFinder-00-00-14/MdcHoughFinder/Hough3D.h.

49{_z0=z0;}

Member Data Documentation

◆ _context

TrkContextEv * Hough3D::_context =NULL
static

Definition at line 68 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Referenced by fit(), and setContext().

◆ _mdcCalibFunSvc

const MdcCalibFunSvc * Hough3D::_mdcCalibFunSvc =NULL
static

Definition at line 69 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Referenced by fit(), and setCalib().

◆ m_debug

int Hough3D::m_debug =0
static

Definition at line 67 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Referenced by fit(), and MdcHoughFinder::initialize().

◆ m_dropTrkChi2Cut

double Hough3D::m_dropTrkChi2Cut =99999
static

Definition at line 72 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Referenced by fit().

◆ m_dropTrkChi2NdfCut

double Hough3D::m_dropTrkChi2NdfCut =99999
static

Definition at line 76 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Referenced by fit().

◆ m_dropTrkDrCut

double Hough3D::m_dropTrkDrCut =1
static

Definition at line 73 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Referenced by fit().

◆ m_dropTrkDzCut

double Hough3D::m_dropTrkDzCut =10
static

Definition at line 74 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Referenced by fit().

◆ m_dropTrkNhitCut

double Hough3D::m_dropTrkNhitCut =9
static

Definition at line 75 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Referenced by fit().

◆ m_qualityFactor

int Hough3D::m_qualityFactor =3
static

Definition at line 71 of file InstallArea/include/MdcHoughFinder/MdcHoughFinder/Hough3D.h.

Referenced by fit().

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