diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_ATLxk.h b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_ATLxk.h index ff071845c6d5b51a298e0458a3005eb6a444962b..f603d96b4ed292e01cc9af0247f9caa88236898d 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_ATLxk.h +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_ATLxk.h @@ -160,8 +160,8 @@ namespace InDet { float m_zmaxB ; float m_ftrig ; float m_ftrigW ; - float r_rmax ; - float r_rstep ; + float m_r_rmax ; + float m_r_rstep ; float m_dzver ; float m_dzdrver ; float m_diver ; @@ -178,26 +178,26 @@ namespace InDet { float m_ipt2C ; float m_COFK ; float m_umax ; - int r_size ; - int r_first ; - int rf_size ; - int rfz_size ; - std::list<InDet::SiSpacePointForSeed*>* r_Sorted ; - std::list<InDet::SiSpacePointForSeed*> rfz_Sorted [ 583] ; - std::list<InDet::SiSpacePointForSeed*> rfzv_Sorted[ 300] ; - std::list<InDet::SiSpacePointForSeed*> l_spforseed ; - std::list<InDet::SiSpacePointForSeed*>::iterator i_spforseed; + int m_r_size ; + int m_r_first ; + int m_rf_size ; + int m_rfz_size ; + std::list<InDet::SiSpacePointForSeed*>* m_r_Sorted ; + std::list<InDet::SiSpacePointForSeed*> m_rfz_Sorted [ 583] ; + std::list<InDet::SiSpacePointForSeed*> m_rfzv_Sorted[ 300] ; + std::list<InDet::SiSpacePointForSeed*> m_l_spforseed ; + std::list<InDet::SiSpacePointForSeed*>::iterator m_i_spforseed; std::list<InDet::SiSpacePointForSeed*>::iterator m_rMin ; int m_ns,m_nsaz,m_nsazv ; int m_fNmax,m_fvNmax ; int m_fNmin,m_fvNmin ; int m_zMin ; - int m_nr ; int* r_index ; int* r_map ; - int m_nrfz , rfz_index [583], rfz_map [583] ; - int m_nrfzv , rfzv_index [300], rfzv_map [300] ; - int rfz_b[583],rfz_t[593],rfz_ib[583][9],rfz_it[583][9] ; - int rfzv_n[300],rfzv_i[300][6] ; + int m_nr ; int* m_r_index ; int* m_r_map ; + int m_nrfz , m_rfz_index [583], m_rfz_map [583] ; + int m_nrfzv , m_rfzv_index [300], m_rfzv_map [300] ; + int m_rfz_b[583],m_rfz_t[593],m_rfz_ib[583][9],m_rfz_it[583][9] ; + int m_rfzv_n[300],m_rfzv_i[300][6] ; float m_sF ; float m_sFv ; @@ -216,9 +216,9 @@ namespace InDet { InDet::SiSpacePointsSeed* m_seedOutput ; - std::list<InDet::SiSpacePointsProSeed*> l_seeds ; - std::list<InDet::SiSpacePointsProSeed*>::iterator i_seed ; - std::list<InDet::SiSpacePointsProSeed*>::iterator i_seede ; + std::list<InDet::SiSpacePointsProSeed*> m_l_seeds ; + std::list<InDet::SiSpacePointsProSeed*>::iterator m_i_seed ; + std::list<InDet::SiSpacePointsProSeed*>::iterator m_i_seede ; std::multimap<float,InDet::SiSpacePointsProSeed*> m_seeds ; std::multimap<float,InDet::SiSpacePointsProSeed*>::iterator m_seed ; @@ -228,7 +228,7 @@ namespace InDet { int m_maxOneSize ; int m_nOneSeeds ; int m_fillOneSeeds; - std::set<float> l_vertex ; + std::set<float> m_l_vertex ; std::vector<std::pair<float,InDet::SiSpacePointForSeed*>> m_CmSp; /////////////////////////////////////////////////////////////////// @@ -311,14 +311,14 @@ namespace InDet { { if(m_nspoint==3) { do { - if(i_seed==i_seede) {findNext(); if(i_seed==i_seede) return 0;} ++i_seed; + if(m_i_seed==m_i_seede) {findNext(); if(m_i_seed==m_i_seede) return 0;} ++m_i_seed; } while(!(*m_seed++).second->set3(*m_seedOutput)); return(m_seedOutput); } else { - if(i_seed==i_seede) {findNext(); if(i_seed==i_seede) return 0;} - (*i_seed++)->set2(*m_seedOutput); return(m_seedOutput); + if(m_i_seed==m_i_seede) {findNext(); if(m_i_seed==m_i_seede) return 0;} + (*m_i_seed++)->set2(*m_seedOutput); return(m_seedOutput); } return 0; } @@ -330,7 +330,7 @@ namespace InDet { if(Zv < m_zminU || Zv > m_zmaxU) return false; if(!m_isvertex) return true; - std::set<float>::iterator v=l_vertex.begin(),ve=l_vertex.end(); + std::set<float>::iterator v=m_l_vertex.begin(),ve=m_l_vertex.end(); float dZmin = fabs((*v)-Zv); for(++v; v!=ve; ++v) { @@ -358,12 +358,12 @@ namespace InDet { if((z*z )<(x*x+y*y)) return 0; } - if(i_spforseed!=l_spforseed.end()) { - sps = (*i_spforseed++); sps->set(sp,r); + if(m_i_spforseed!=m_l_spforseed.end()) { + sps = (*m_i_spforseed++); sps->set(sp,r); } else { - l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); - i_spforseed = l_spforseed.end(); + m_l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); + m_i_spforseed = m_l_spforseed.end(); } return sps; @@ -378,13 +378,13 @@ namespace InDet { { SiSpacePointForSeed* p3 = 0; - if(i_seede!=l_seeds.end()) { - SiSpacePointsProSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsProSeed* s = (*m_i_seede++); s->set(p1,p2,p3,z); } else { - l_seeds.push_back(new SiSpacePointsProSeed(p1,p2,p3,z)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsProSeed(p1,p2,p3,z)); + m_i_seede = m_l_seeds.end(); } } diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_BeamGas.h b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_BeamGas.h index d075566e7090cd9a33e415db0225e2e1e18eca30..104141d8a10a847c030ff6e99999907e4e3cde99 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_BeamGas.h +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_BeamGas.h @@ -147,8 +147,8 @@ namespace InDet { float m_dzdrmax ; float m_zmin ; float m_zmax ; - float r_rmax ; - float r_rstep ; + float m_r_rmax ; + float m_r_rstep ; float m_dzver ; float m_dzdrver ; float m_diver ; @@ -160,24 +160,24 @@ namespace InDet { float m_COF ; float m_K ; - int r_size ; - int rf_size ; - int rfz_size ; - std::list<InDet::SiSpacePointForSeed*>* r_Sorted ; - std::list<InDet::SiSpacePointForSeed*> rf_Sorted [ 53] ; - std::list<InDet::SiSpacePointForSeed*> rfz_Sorted [ 583] ; - std::list<InDet::SiSpacePointForSeed*> l_spforseed ; - std::list<InDet::SiSpacePointForSeed*>::iterator i_spforseed; + int m_r_size ; + int m_rf_size ; + int m_rfz_size ; + std::list<InDet::SiSpacePointForSeed*>* m_r_Sorted ; + std::list<InDet::SiSpacePointForSeed*> m_rf_Sorted [ 53] ; + std::list<InDet::SiSpacePointForSeed*> m_rfz_Sorted [ 583] ; + std::list<InDet::SiSpacePointForSeed*> m_l_spforseed ; + std::list<InDet::SiSpacePointForSeed*>::iterator m_i_spforseed; std::list<InDet::SiSpacePointForSeed*>::iterator m_rMin ; int m_ns,m_nsaz ; int m_fNmax ; int m_fNmin ; int m_zMin ; - int m_nr ; int* r_index ; int* r_map ; - int m_nrf , rf_index [ 53], rf_map [ 53] ; - int m_nrfz , rfz_index [583], rfz_map [583] ; - int rfz_b[583],rfz_t[593],rfz_ib[583][12],rfz_it[583][12] ; + int m_nr ; int* m_r_index ; int* m_r_map ; + int m_nrf , m_rf_index [ 53], m_rf_map [ 53] ; + int m_nrfz , m_rfz_index [583], m_rfz_map [583] ; + int m_rfz_b[583],m_rfz_t[593],m_rfz_ib[583][12],m_rfz_it[583][12] ; float m_sF ; /////////////////////////////////////////////////////////////////// @@ -193,9 +193,9 @@ namespace InDet { float * m_V ; float * m_Zo ; - std::list<InDet::SiSpacePointsSeed*> l_seeds ; - std::list<InDet::SiSpacePointsSeed*>::iterator i_seed ; - std::list<InDet::SiSpacePointsSeed*>::iterator i_seede ; + std::list<InDet::SiSpacePointsSeed*> m_l_seeds ; + std::list<InDet::SiSpacePointsSeed*>::iterator m_i_seed ; + std::list<InDet::SiSpacePointsSeed*>::iterator m_i_seede ; std::multimap<float,InDet::SiSpacePointsSeed*> m_mapOneSeeds; InDet::SiSpacePointsSeed* m_OneSeeds ; @@ -272,10 +272,10 @@ namespace InDet { inline const SiSpacePointsSeed* SiSpacePointsSeedMaker_BeamGas::next() { - if(i_seed==i_seede) { - findNext(); if(i_seed==i_seede) return 0; + if(m_i_seed==m_i_seede) { + findNext(); if(m_i_seed==m_i_seede) return 0; } - return(*i_seed++); + return(*m_i_seed++); } inline bool SiSpacePointsSeedMaker_BeamGas::isZCompatible @@ -295,12 +295,12 @@ namespace InDet { float r[3]; convertToBeamFrameWork(sp,r); - if(i_spforseed!=l_spforseed.end()) { - sps = (*i_spforseed++); sps->set(sp,r); + if(m_i_spforseed!=m_l_spforseed.end()) { + sps = (*m_i_spforseed++); sps->set(sp,r); } else { - l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); - i_spforseed = l_spforseed.end(); + m_l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); + m_i_spforseed = m_l_spforseed.end(); } return sps; @@ -314,16 +314,16 @@ namespace InDet { (const Trk::SpacePoint*& p1,const Trk::SpacePoint*& p2, const float& z) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); s->erase ( ); s->add (p1); s->add (p2); s->setZVertex(double(z)); } else { - l_seeds.push_back(new SiSpacePointsSeed(p1,p2,z)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(p1,p2,z)); + m_i_seede = m_l_seeds.end(); } } @@ -335,8 +335,8 @@ namespace InDet { (const Trk::SpacePoint*& p1,const Trk::SpacePoint*& p2, const Trk::SpacePoint*& p3,const float& z) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); s->erase ( ); s->add (p1); s->add (p2); @@ -344,8 +344,8 @@ namespace InDet { s->setZVertex(double(z)); } else { - l_seeds.push_back(new SiSpacePointsSeed(p1,p2,p3,z)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(p1,p2,p3,z)); + m_i_seede = m_l_seeds.end(); } } @@ -361,13 +361,13 @@ namespace InDet { for(; l!=le; ++l) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); *s = *(*l).second; } else { - l_seeds.push_back(new SiSpacePointsSeed(*(*l).second)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(*(*l).second)); + m_i_seede = m_l_seeds.end(); } } } diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_Cosmic.h b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_Cosmic.h index 0d6a56ea9c1566093b914f4c4925560eacb351ab..d4427f2e4c4e45295d224efe73dadb89d7b53618 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_Cosmic.h +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_Cosmic.h @@ -142,31 +142,31 @@ namespace InDet { float m_dzdrmax ; float m_zmin ; float m_zmax ; - float r_rmax ; - float r_rstep ; + float m_r_rmax ; + float m_r_rstep ; float m_diver ; float m_dazmax ; float m_diverpps ; float m_diversss ; float m_etamin, m_etamax ; - int r_size ; - int rf_size ; - int rfz_size ; - std::list<InDet::SiSpacePointForSeed*>* r_Sorted ; - std::list<InDet::SiSpacePointForSeed*> rf_Sorted [ 53] ; - std::list<InDet::SiSpacePointForSeed*> rfz_Sorted [ 583] ; - std::list<InDet::SiSpacePointForSeed*> l_spforseed ; - std::list<InDet::SiSpacePointForSeed*>::iterator i_spforseed; + int m_r_size ; + int m_rf_size ; + int m_rfz_size ; + std::list<InDet::SiSpacePointForSeed*>* m_r_Sorted ; + std::list<InDet::SiSpacePointForSeed*> m_rf_Sorted [ 53] ; + std::list<InDet::SiSpacePointForSeed*> m_rfz_Sorted [ 583] ; + std::list<InDet::SiSpacePointForSeed*> m_l_spforseed ; + std::list<InDet::SiSpacePointForSeed*>::iterator m_i_spforseed; std::list<InDet::SiSpacePointForSeed*>::iterator m_rMin ; int m_ns,m_nsaz ; int m_fNmax ; int m_fNmin ; int m_zMin ; - int m_nr ; int* r_index ; int* r_map ; - int m_nrf , rf_index [ 53], rf_map [ 53] ; - int m_nrfz , rfz_index [583], rfz_map [583] ; - int rfz_b[583],rfz_t[593],rfz_ib[583][12],rfz_it[583][12] ; + int m_nr ; int* m_r_index ; int* m_r_map ; + int m_nrf , m_rf_index [ 53], m_rf_map [ 53] ; + int m_nrfz , m_rfz_index [583], m_rfz_map [583] ; + int m_rfz_b[583],m_rfz_t[593],m_rfz_ib[583][12],m_rfz_it[583][12] ; float m_sF ; /////////////////////////////////////////////////////////////////// @@ -181,9 +181,9 @@ namespace InDet { float * m_U ; float * m_V ; - std::multimap<float,InDet::SiSpacePointsSeed*> l_seeds; - std::multimap<float,InDet::SiSpacePointsSeed*>::iterator i_seed ; - std::multimap<float,InDet::SiSpacePointsSeed*>::iterator i_seede; + std::multimap<float,InDet::SiSpacePointsSeed*> m_l_seeds; + std::multimap<float,InDet::SiSpacePointsSeed*>::iterator m_i_seed ; + std::multimap<float,InDet::SiSpacePointsSeed*>::iterator m_i_seede; InDet::SiSpacePointsSeed* m_seeds; /////////////////////////////////////////////////////////////////// @@ -249,8 +249,8 @@ namespace InDet { inline const SiSpacePointsSeed* SiSpacePointsSeedMaker_Cosmic::next() { - if(i_seed==i_seede) return 0; - SiSpacePointsSeed* sp = (*i_seed).second; ++i_seed; + if(m_i_seed==m_i_seede) return 0; + SiSpacePointsSeed* sp = (*m_i_seed).second; ++m_i_seed; return(sp); } @@ -267,12 +267,12 @@ namespace InDet { float(sp->globalPosition().y()), float(sp->globalPosition().z())}; - if(i_spforseed!=l_spforseed.end()) { - sps = (*i_spforseed++); sps->set(sp,r); + if(m_i_spforseed!=m_l_spforseed.end()) { + sps = (*m_i_spforseed++); sps->set(sp,r); } else { - l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); - i_spforseed = l_spforseed.end(); + m_l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); + m_i_spforseed = m_l_spforseed.end(); } return sps; @@ -289,12 +289,12 @@ namespace InDet { float r[3]; r[0]=sp->globalPosition().x(); r[1]=sp->globalPosition().y(); r[2]=sp->globalPosition().z(); - if(i_spforseed!=l_spforseed.end()) { - sps = (*i_spforseed++); sps->set(sp,r,sc); + if(m_i_spforseed!=m_l_spforseed.end()) { + sps = (*m_i_spforseed++); sps->set(sp,r,sc); } else { - l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r,sc))); - i_spforseed = l_spforseed.end(); + m_l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r,sc))); + m_i_spforseed = m_l_spforseed.end(); } return sps; @@ -315,21 +315,21 @@ namespace InDet { m_seeds[m_nseeds].add (p1); m_seeds[m_nseeds].add (p2); m_seeds[m_nseeds].setZVertex(0.); - l_seeds.insert(std::make_pair(z,m_seeds+m_nseeds)); + m_l_seeds.insert(std::make_pair(z,m_seeds+m_nseeds)); ++m_nseeds; } else { - std::multimap<float,InDet::SiSpacePointsSeed*>::reverse_iterator l = l_seeds.rbegin(); + std::multimap<float,InDet::SiSpacePointsSeed*>::reverse_iterator l = m_l_seeds.rbegin(); if((*l).first <= z) return; InDet::SiSpacePointsSeed* s = (*l).second; - l_seeds.erase((*l).first); + m_l_seeds.erase((*l).first); s->erase ( ); s->add (p1); s->add (p2); s->setZVertex(0.); - l_seeds.insert(std::make_pair(z,s)); + m_l_seeds.insert(std::make_pair(z,s)); } } @@ -348,22 +348,22 @@ namespace InDet { m_seeds[m_nseeds].add (p2); m_seeds[m_nseeds].add (p3); m_seeds[m_nseeds].setZVertex(0.); - l_seeds.insert(std::make_pair(z,m_seeds+m_nseeds)); + m_l_seeds.insert(std::make_pair(z,m_seeds+m_nseeds)); ++m_nseeds; } else { - std::multimap<float,InDet::SiSpacePointsSeed*>::reverse_iterator l = l_seeds.rbegin(); + std::multimap<float,InDet::SiSpacePointsSeed*>::reverse_iterator l = m_l_seeds.rbegin(); if((*l).first <= z) return; InDet::SiSpacePointsSeed* s = (*l).second; - l_seeds.erase((*l).first); + m_l_seeds.erase((*l).first); s->erase ( ); s->add (p1); s->add (p2); s->add (p3); s->setZVertex(0.); - l_seeds.insert(std::make_pair(z,s)); + m_l_seeds.insert(std::make_pair(z,s)); } } diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_HeavyIon.h b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_HeavyIon.h index 01b28d07e3005d20dacec20fb9572f0c69e6c860..07874b53764a259d238aa91e3311c8927783bb25 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_HeavyIon.h +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_HeavyIon.h @@ -149,8 +149,8 @@ namespace InDet { float m_zmaxU ; float m_ftrig ; float m_ftrigW ; - float r_rmax ; - float r_rstep ; + float m_r_rmax ; + float m_r_rstep ; float m_dzver ; float m_dzdrver ; float m_diver ; @@ -162,25 +162,25 @@ namespace InDet { float m_ipt2 ; float m_COF ; float m_K ; - int r_size ; - int rf_size ; - int rfz_size ; - std::list<InDet::SiSpacePointForSeed*>* r_Sorted ; - std::list<InDet::SiSpacePointForSeed*> rfz_Sorted [ 583] ; - std::list<InDet::SiSpacePointForSeed*> rfzv_Sorted[ 300] ; - std::list<InDet::SiSpacePointForSeed*> l_spforseed ; - std::list<InDet::SiSpacePointForSeed*>::iterator i_spforseed; + int m_r_size ; + int m_rf_size ; + int m_rfz_size ; + std::list<InDet::SiSpacePointForSeed*>* m_r_Sorted ; + std::list<InDet::SiSpacePointForSeed*> m_rfz_Sorted [ 583] ; + std::list<InDet::SiSpacePointForSeed*> m_rfzv_Sorted[ 300] ; + std::list<InDet::SiSpacePointForSeed*> m_l_spforseed ; + std::list<InDet::SiSpacePointForSeed*>::iterator m_i_spforseed; std::list<InDet::SiSpacePointForSeed*>::iterator m_rMin ; int m_ns,m_nsaz,m_nsazv ; int m_fNmax,m_fvNmax ; int m_fNmin,m_fvNmin ; int m_zMin ; - int m_nr ; int* r_index ; int* r_map ; - int m_nrfz , rfz_index [583], rfz_map [583] ; - int m_nrfzv , rfzv_index [300], rfzv_map [300] ; - int rfz_b[583],rfz_t[593],rfz_ib[583][9],rfz_it[583][9] ; - int rfzv_n[300],rfzv_i[300][6] ; + int m_nr ; int* m_r_index ; int* m_r_map ; + int m_nrfz , m_rfz_index [583], m_rfz_map [583] ; + int m_nrfzv , m_rfzv_index [300], m_rfzv_map [300] ; + int m_rfz_b[583],m_rfz_t[593],m_rfz_ib[583][9],m_rfz_it[583][9] ; + int m_rfzv_n[300],m_rfzv_i[300][6] ; float m_sF ; float m_sFv ; @@ -197,15 +197,15 @@ namespace InDet { float * m_V ; float * m_Er ; - std::list<InDet::SiSpacePointsSeed*> l_seeds ; - std::list<InDet::SiSpacePointsSeed*>::iterator i_seed ; - std::list<InDet::SiSpacePointsSeed*>::iterator i_seede ; + std::list<InDet::SiSpacePointsSeed*> m_l_seeds ; + std::list<InDet::SiSpacePointsSeed*>::iterator m_i_seed ; + std::list<InDet::SiSpacePointsSeed*>::iterator m_i_seede ; std::multimap<float,InDet::SiSpacePointsSeed*> m_mapOneSeeds; InDet::SiSpacePointsSeed* m_OneSeeds ; int m_maxOneSize ; int m_nOneSeeds ; - std::list<float> l_vertex ; + std::list<float> m_l_vertex ; /////////////////////////////////////////////////////////////////// // Beam geometry @@ -290,8 +290,8 @@ namespace InDet { inline const SiSpacePointsSeed* SiSpacePointsSeedMaker_HeavyIon::next() { - if(i_seed==i_seede) {findNext(); if(i_seed==i_seede) return 0;} - return(*i_seed++); + if(m_i_seed==m_i_seede) {findNext(); if(m_i_seed==m_i_seede) return 0;} + return(*m_i_seed++); } inline bool SiSpacePointsSeedMaker_HeavyIon::isZCompatible @@ -300,7 +300,7 @@ namespace InDet { if(Zv < m_zmin || Zv > m_zmax) return false; if(!m_izvertex ) return true ; - std::list<float>::iterator v=l_vertex.begin(),ve=l_vertex.end(); + std::list<float>::iterator v=m_l_vertex.begin(),ve=m_l_vertex.end(); float dZmin = fabs((*v)-Zv); ++v; @@ -313,7 +313,7 @@ namespace InDet { inline float SiSpacePointsSeedMaker_HeavyIon::dZVertexMin(float& Z) { float dZm = 1.E10; - std::list<float>::iterator v=l_vertex.begin(),ve=l_vertex.end(); + std::list<float>::iterator v=m_l_vertex.begin(),ve=m_l_vertex.end(); for(; v!=ve; ++v) {float dZ = fabs((*v)-Z); if(dZ<dZm) dZm = dZ;} return dZm; @@ -330,12 +330,12 @@ namespace InDet { float r[3]; convertToBeamFrameWork(sp,r); - if(i_spforseed!=l_spforseed.end()) { - sps = (*i_spforseed++); sps->set(sp,r); + if(m_i_spforseed!=m_l_spforseed.end()) { + sps = (*m_i_spforseed++); sps->set(sp,r); } else { - l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); - i_spforseed = l_spforseed.end(); + m_l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); + m_i_spforseed = m_l_spforseed.end(); } return sps; @@ -349,16 +349,16 @@ namespace InDet { (const Trk::SpacePoint*& p1,const Trk::SpacePoint*& p2, const float& z) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); s->erase ( ); s->add (p1); s->add (p2); s->setZVertex(double(z)); } else { - l_seeds.push_back(new SiSpacePointsSeed(p1,p2,z)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(p1,p2,z)); + m_i_seede = m_l_seeds.end(); } } @@ -370,8 +370,8 @@ namespace InDet { (const Trk::SpacePoint*& p1,const Trk::SpacePoint*& p2, const Trk::SpacePoint*& p3,const float& z) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); s->erase ( ); s->add (p1); s->add (p2); @@ -379,8 +379,8 @@ namespace InDet { s->setZVertex(double(z)); } else { - l_seeds.push_back(new SiSpacePointsSeed(p1,p2,p3,z)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(p1,p2,p3,z)); + m_i_seede = m_l_seeds.end(); } } @@ -396,13 +396,13 @@ namespace InDet { for(; l!=le; ++l) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); *s = *(*l).second; } else { - l_seeds.push_back(new SiSpacePointsSeed(*(*l).second)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(*(*l).second)); + m_i_seede = m_l_seeds.end(); } } } diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_ITK.h b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_ITK.h index a2414a7ad387eb54c44b1df81b3f68c2dc754811..8a469ff1eefb3047217a7af50949736f99c8a613 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_ITK.h +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_ITK.h @@ -156,9 +156,9 @@ namespace InDet { float m_zmaxB ; float m_ftrig ; float m_ftrigW ; - float r_rmax ; - float r_rmin ; - float r_rstep ; + float m_r_rmax ; + float m_r_rmin ; + float m_r_rstep ; float m_dzver ; float m_dzdrver ; float m_diver ; @@ -176,26 +176,26 @@ namespace InDet { float m_COFK ; float m_umax ; float m_dzmaxPPP ; - int r_size ; - int r_first ; - int rf_size ; - int rfz_size ; - std::list<InDet::SiSpacePointForSeedITK*>* r_Sorted ; - std::list<InDet::SiSpacePointForSeedITK*> rfz_Sorted [ 583] ; - std::list<InDet::SiSpacePointForSeedITK*> rfzv_Sorted[ 300] ; - std::list<InDet::SiSpacePointForSeedITK*> l_spforseed ; - std::list<InDet::SiSpacePointForSeedITK*>::iterator i_spforseed; - std::list<InDet::SiSpacePointForSeedITK*>::iterator m_rMin ; + int m_r_size ; + int m_r_first ; + int m_rf_size ; + int m_rfz_size ; + std::list<InDet::SiSpacePointForSeedITK*>* m_r_Sorted ; + std::list<InDet::SiSpacePointForSeedITK*> m_rfz_Sorted [ 583] ; + std::list<InDet::SiSpacePointForSeedITK*> m_rfzv_Sorted[ 300] ; + std::list<InDet::SiSpacePointForSeedITK*> m_l_spforseed ; + std::list<InDet::SiSpacePointForSeedITK*>::iterator m_i_spforseed; + std::list<InDet::SiSpacePointForSeedITK*>::iterator m_rMin ; int m_ns,m_nsaz,m_nsazv ; int m_fNmax,m_fvNmax ; int m_fNmin,m_fvNmin ; int m_zMin ; - int m_nr ; int* r_index ; int* r_map ; - int m_nrfz , rfz_index [583], rfz_map [583] ; - int m_nrfzv , rfzv_index [300], rfzv_map [300] ; - int rfz_b[583],rfz_t[593],rfz_ib[583][9],rfz_it[583][9] ; - int rfzv_n[300],rfzv_i[300][6] ; + int m_nr ; int* m_r_index ; int* m_r_map ; + int m_nrfz , m_rfz_index [583], m_rfz_map [583] ; + int m_nrfzv , m_rfzv_index [300], m_rfzv_map [300] ; + int m_rfz_b[583],m_rfz_t[593],m_rfz_ib[583][9],m_rfz_it[583][9] ; + int m_rfzv_n[300],m_rfzv_i[300][6] ; float m_sF ; float m_sFv ; @@ -216,9 +216,9 @@ namespace InDet { InDet::SiSpacePointsSeed* m_seedOutput ; - std::list<InDet::SiSpacePointsProSeedITK*> l_seeds; - std::list<InDet::SiSpacePointsProSeedITK*>::iterator i_seed ; - std::list<InDet::SiSpacePointsProSeedITK*>::iterator i_seede ; + std::list<InDet::SiSpacePointsProSeedITK*> m_l_seeds; + std::list<InDet::SiSpacePointsProSeedITK*>::iterator m_i_seed ; + std::list<InDet::SiSpacePointsProSeedITK*>::iterator m_i_seede ; std::multimap<float,InDet::SiSpacePointsProSeedITK*> m_seeds ; std::multimap<float,InDet::SiSpacePointsProSeedITK*>::iterator m_seed ; @@ -228,7 +228,7 @@ namespace InDet { int m_maxOneSize ; int m_nOneSeeds ; int m_fillOneSeeds; - std::set<float> l_vertex ; + std::set<float> m_l_vertex ; std::vector<std::pair<float,InDet::SiSpacePointForSeedITK*>> m_CmSp; /////////////////////////////////////////////////////////////////// @@ -318,14 +318,14 @@ namespace InDet { { if(m_nspoint==3) { do { - if(i_seed==i_seede) {findNext(); if(i_seed==i_seede) return 0;} ++i_seed; + if(m_i_seed==m_i_seede) {findNext(); if(m_i_seed==m_i_seede) return 0;} ++m_i_seed; } while(!(*m_seed++).second->set3(*m_seedOutput)); return(m_seedOutput); } else { - if(i_seed==i_seede) {findNext(); if(i_seed==i_seede) return 0;} - (*i_seed++)->set2(*m_seedOutput); return(m_seedOutput); + if(m_i_seed==m_i_seede) {findNext(); if(m_i_seed==m_i_seede) return 0;} + (*m_i_seed++)->set2(*m_seedOutput); return(m_seedOutput); } return 0; } @@ -337,7 +337,7 @@ namespace InDet { if(Zv < m_zminU || Zv > m_zmaxU) return false; if(!m_isvertex) return true; - std::set<float>::iterator v=l_vertex.begin(),ve=l_vertex.end(); + std::set<float>::iterator v=m_l_vertex.begin(),ve=m_l_vertex.end(); float dZmin = fabs((*v)-Zv); for(++v; v!=ve; ++v) { @@ -365,12 +365,12 @@ namespace InDet { if((z*z )<(x*x+y*y)) return 0; } - if(i_spforseed!=l_spforseed.end()) { - sps = (*i_spforseed++); sps->set(sp,r); + if(m_i_spforseed!=m_l_spforseed.end()) { + sps = (*m_i_spforseed++); sps->set(sp,r); } else { - l_spforseed.push_back((sps=new SiSpacePointForSeedITK(sp,r))); - i_spforseed = l_spforseed.end(); + m_l_spforseed.push_back((sps=new SiSpacePointForSeedITK(sp,r))); + m_i_spforseed = m_l_spforseed.end(); } return sps; @@ -385,13 +385,13 @@ namespace InDet { { SiSpacePointForSeedITK* p3 = 0; - if(i_seede!=l_seeds.end()) { - SiSpacePointsProSeedITK* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsProSeedITK* s = (*m_i_seede++); s->set(p1,p2,p3,z); } else { - l_seeds.push_back(new SiSpacePointsProSeedITK(p1,p2,p3,z)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsProSeedITK(p1,p2,p3,z)); + m_i_seede = m_l_seeds.end(); } } diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_LowMomentum.h b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_LowMomentum.h index 2fc6e00e2c52b2e2e780c6e3dcfa1ab10c0399fa..7e9f1512c02cf0221408699849ead107e7456eae 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_LowMomentum.h +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_LowMomentum.h @@ -150,8 +150,8 @@ namespace InDet { float m_zmax ; float m_zminU ; float m_zmaxU ; - float r_rmax ; - float r_rstep ; + float m_r_rmax ; + float m_r_rstep ; float m_dzver ; float m_dzdrver ; float m_diver ; @@ -162,22 +162,22 @@ namespace InDet { float m_iptmin ; float m_iptmax ; - int r_size ; - int rf_size ; - int rfz_size ; - std::list<InDet::SiSpacePointForSeed*>* r_Sorted ; - std::list<InDet::SiSpacePointForSeed*> rfz_Sorted [ 220] ; - std::list<InDet::SiSpacePointForSeed*> l_spforseed ; - std::list<InDet::SiSpacePointForSeed*>::iterator i_spforseed; + int m_r_size ; + int m_rf_size ; + int m_rfz_size ; + std::list<InDet::SiSpacePointForSeed*>* m_r_Sorted ; + std::list<InDet::SiSpacePointForSeed*> m_rfz_Sorted [ 220] ; + std::list<InDet::SiSpacePointForSeed*> m_l_spforseed ; + std::list<InDet::SiSpacePointForSeed*>::iterator m_i_spforseed; std::list<InDet::SiSpacePointForSeed*>::iterator m_rMin ; int m_ns,m_nsaz ; int m_fNmax ; int m_fNmin ; int m_zMin ; - int m_nr ; int* r_index ; int* r_map ; - int m_nrfz , rfz_index [220], rfz_map [220] ; - int rfz_b[220],rfz_t[220],rfz_ib[220][9],rfz_it[220][9] ; + int m_nr ; int* m_r_index ; int* m_r_map ; + int m_nrfz , m_rfz_index [220], m_rfz_map [220] ; + int m_rfz_b[220],m_rfz_t[220],m_rfz_ib[220][9],m_rfz_it[220][9] ; float m_sF ; /////////////////////////////////////////////////////////////////// @@ -193,15 +193,15 @@ namespace InDet { float * m_V ; float * m_Zo ; - std::list<InDet::SiSpacePointsSeed*> l_seeds ; - std::list<InDet::SiSpacePointsSeed*>::iterator i_seed ; - std::list<InDet::SiSpacePointsSeed*>::iterator i_seede ; + std::list<InDet::SiSpacePointsSeed*> m_l_seeds ; + std::list<InDet::SiSpacePointsSeed*>::iterator m_i_seed ; + std::list<InDet::SiSpacePointsSeed*>::iterator m_i_seede ; std::multimap<float,InDet::SiSpacePointsSeed*> m_mapOneSeeds; InDet::SiSpacePointsSeed* m_OneSeeds ; int m_maxOneSize; int m_nOneSeeds ; - std::list<float> l_vertex ; + std::list<float> m_l_vertex ; /////////////////////////////////////////////////////////////////// // Beam geometry @@ -276,11 +276,11 @@ namespace InDet { inline const SiSpacePointsSeed* SiSpacePointsSeedMaker_LowMomentum::next() { - if(i_seed==i_seede) { + if(m_i_seed==m_i_seede) { findNext(); - if(i_seed==i_seede) return 0; + if(m_i_seed==m_i_seede) return 0; } - return(*i_seed++); + return(*m_i_seed++); } inline bool SiSpacePointsSeedMaker_LowMomentum::isZCompatible @@ -288,7 +288,7 @@ namespace InDet { { if(Zv < m_zmin || Zv > m_zmax) return false; - std::list<float>::iterator v=l_vertex.begin(),ve=l_vertex.end(); + std::list<float>::iterator v=m_l_vertex.begin(),ve=m_l_vertex.end(); if(v==ve) return true; float dZmin = fabs((*v)-Zv); ++v; @@ -310,12 +310,12 @@ namespace InDet { float r[3]; convertToBeamFrameWork(sp,r); - if(i_spforseed!=l_spforseed.end()) { - sps = (*i_spforseed++); sps->set(sp,r); + if(m_i_spforseed!=m_l_spforseed.end()) { + sps = (*m_i_spforseed++); sps->set(sp,r); } else { - l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); - i_spforseed = l_spforseed.end(); + m_l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); + m_i_spforseed = m_l_spforseed.end(); } return sps; @@ -329,16 +329,16 @@ namespace InDet { (const Trk::SpacePoint*& p1,const Trk::SpacePoint*& p2, const float& z) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); s->erase ( ); s->add (p1); s->add (p2); s->setZVertex(double(z)); } else { - l_seeds.push_back(new SiSpacePointsSeed(p1,p2,z)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(p1,p2,z)); + m_i_seede = m_l_seeds.end(); } } @@ -350,8 +350,8 @@ namespace InDet { (const Trk::SpacePoint*& p1,const Trk::SpacePoint*& p2, const Trk::SpacePoint*& p3,const float& z) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); s->erase ( ); s->add (p1); s->add (p2); @@ -359,8 +359,8 @@ namespace InDet { s->setZVertex(double(z)); } else { - l_seeds.push_back(new SiSpacePointsSeed(p1,p2,p3,z)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(p1,p2,p3,z)); + m_i_seede = m_l_seeds.end(); } } @@ -376,13 +376,13 @@ namespace InDet { for(; l!=le; ++l) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); *s = *(*l).second; } else { - l_seeds.push_back(new SiSpacePointsSeed(*(*l).second)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(*(*l).second)); + m_i_seede = m_l_seeds.end(); } } } diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_Trigger.h b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_Trigger.h index 1ccced7147d934031ff032f435a93b8e927ea819..e5dd983dbd51a8edcd57104c74510637a4921f72 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_Trigger.h +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/SiSpacePointsSeedTool_xk/SiSpacePointsSeedMaker_Trigger.h @@ -148,8 +148,8 @@ namespace InDet { float m_zmaxB ; float m_ftrig ; float m_ftrigW ; - float r_rmax ; - float r_rstep ; + float m_r_rmax ; + float m_r_rstep ; float m_dzver ; float m_dzdrver ; float m_diver ; @@ -161,25 +161,25 @@ namespace InDet { float m_ipt2 ; float m_COF ; float m_K ; - int r_size ; - int rf_size ; - int rfz_size ; - std::list<InDet::SiSpacePointForSeed*>* r_Sorted ; - std::list<InDet::SiSpacePointForSeed*> rfz_Sorted [ 583] ; - std::list<InDet::SiSpacePointForSeed*> rfzv_Sorted[ 300] ; - std::list<InDet::SiSpacePointForSeed*> l_spforseed ; - std::list<InDet::SiSpacePointForSeed*>::iterator i_spforseed; + int m_r_size ; + int m_rf_size ; + int m_rfz_size ; + std::list<InDet::SiSpacePointForSeed*>* m_r_Sorted ; + std::list<InDet::SiSpacePointForSeed*> m_rfz_Sorted [ 583] ; + std::list<InDet::SiSpacePointForSeed*> m_rfzv_Sorted[ 300] ; + std::list<InDet::SiSpacePointForSeed*> m_l_spforseed ; + std::list<InDet::SiSpacePointForSeed*>::iterator m_i_spforseed; std::list<InDet::SiSpacePointForSeed*>::iterator m_rMin ; int m_ns,m_nsaz,m_nsazv ; int m_fNmax,m_fvNmax ; int m_fNmin,m_fvNmin ; int m_zMin ; - int m_nr ; int* r_index ; int* r_map ; - int m_nrfz , rfz_index [583], rfz_map [583] ; - int m_nrfzv , rfzv_index [300], rfzv_map [300] ; - int rfz_b[583],rfz_t[593],rfz_ib[583][9],rfz_it[583][9] ; - int rfzv_n[300],rfzv_i[300][6] ; + int m_nr ; int* m_r_index ; int* m_r_map ; + int m_nrfz , m_rfz_index [583], m_rfz_map [583] ; + int m_nrfzv , m_rfzv_index [300], m_rfzv_map [300] ; + int m_rfz_b[583],m_rfz_t[593],m_rfz_ib[583][9],m_rfz_it[583][9] ; + int m_rfzv_n[300],m_rfzv_i[300][6] ; float m_sF ; float m_sFv ; @@ -196,9 +196,9 @@ namespace InDet { float * m_V ; float * m_Er ; - std::list<InDet::SiSpacePointsSeed*> l_seeds ; - std::list<InDet::SiSpacePointsSeed*>::iterator i_seed ; - std::list<InDet::SiSpacePointsSeed*>::iterator i_seede ; + std::list<InDet::SiSpacePointsSeed*> m_l_seeds ; + std::list<InDet::SiSpacePointsSeed*>::iterator m_i_seed ; + std::list<InDet::SiSpacePointsSeed*>::iterator m_i_seede ; std::multimap<float,InDet::SiSpacePointsSeed*> m_mapOneSeeds; std::multimap<float,InDet::SiSpacePointsSeed*> m_mapSeeds ; @@ -208,7 +208,7 @@ namespace InDet { InDet::SiSpacePointsSeed* m_OneSeeds ; int m_maxOneSize ; int m_nOneSeeds ; - std::list<float> l_vertex ; + std::list<float> m_l_vertex ; /////////////////////////////////////////////////////////////////// // Beam geometry @@ -284,9 +284,9 @@ namespace InDet { inline const SiSpacePointsSeed* SiSpacePointsSeedMaker_Trigger::next() { - if(i_seed==i_seede) {findNext(); if(i_seed==i_seede) return 0;} - if(m_mode==0 || m_mode==1) return(*i_seed++); - ++i_seed; + if(m_i_seed==m_i_seede) {findNext(); if(m_i_seed==m_i_seede) return 0;} + if(m_mode==0 || m_mode==1) return(*m_i_seed++); + ++m_i_seed; return (*m_seed++).second; } @@ -295,7 +295,7 @@ namespace InDet { { if(Zv < m_zmin || Zv > m_zmax) return false; - std::list<float>::iterator v=l_vertex.begin(),ve=l_vertex.end(); + std::list<float>::iterator v=m_l_vertex.begin(),ve=m_l_vertex.end(); if(v==ve) return true; float dZmin = fabs((*v)-Zv); ++v; @@ -308,7 +308,7 @@ namespace InDet { inline float SiSpacePointsSeedMaker_Trigger::dZVertexMin(float& Z) { - std::list<float>::iterator v=l_vertex.begin(),ve=l_vertex.end(); + std::list<float>::iterator v=m_l_vertex.begin(),ve=m_l_vertex.end(); if(v==ve) return 0.; float dZm = 1.E10; @@ -327,12 +327,12 @@ namespace InDet { float r[3]; convertToBeamFrameWork(sp,r); - if(i_spforseed!=l_spforseed.end()) { - sps = (*i_spforseed++); sps->set(sp,r); + if(m_i_spforseed!=m_l_spforseed.end()) { + sps = (*m_i_spforseed++); sps->set(sp,r); } else { - l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); - i_spforseed = l_spforseed.end(); + m_l_spforseed.push_back((sps=new SiSpacePointForSeed(sp,r))); + m_i_spforseed = m_l_spforseed.end(); } return sps; @@ -346,16 +346,16 @@ namespace InDet { (const Trk::SpacePoint*& p1,const Trk::SpacePoint*& p2, const float& z) { - if(i_seede!=l_seeds.end()) { - SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + SiSpacePointsSeed* s = (*m_i_seede++); s->erase ( ); s->add (p1); s->add (p2); s->setZVertex(double(z)); } else { - l_seeds.push_back(new SiSpacePointsSeed(p1,p2,z)); - i_seede = l_seeds.end(); + m_l_seeds.push_back(new SiSpacePointsSeed(p1,p2,z)); + m_i_seede = m_l_seeds.end(); } } @@ -378,16 +378,16 @@ namespace InDet { (*s0->spacePoints().begin())->clusterList().second ? q+=1000. : q+=10000.; } - if(i_seede!=l_seeds.end()) { - InDet::SiSpacePointsSeed* s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + InDet::SiSpacePointsSeed* s = (*m_i_seede++); *s = *s0; m_mapSeeds.insert(std::make_pair(q,s)); } else { InDet::SiSpacePointsSeed* s = new InDet::SiSpacePointsSeed(*s0); - l_seeds.push_back(s); - i_seede = l_seeds.end(); + m_l_seeds.push_back(s); + m_i_seede = m_l_seeds.end(); m_mapSeeds.insert(std::make_pair(q,s)); } } diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_ATLxk.cxx b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_ATLxk.cxx index 7b1342da4d493444132fc0d87667e3f4af793b47..3d9fef9d596029883632af5fe18cdb2693784aef 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_ATLxk.cxx +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_ATLxk.cxx @@ -67,11 +67,11 @@ InDet::SiSpacePointsSeedMaker_ATLxk::SiSpacePointsSeedMaker_ATLxk m_diversss = 50 ; m_divermax = 20. ; m_dazmax = .02 ; - r_rmax = 600. ; - r_rstep = 2. ; - r_Sorted = 0 ; - r_index = 0 ; - r_map = 0 ; + m_r_rmax = 600. ; + m_r_rstep = 2. ; + m_r_Sorted = 0 ; + m_r_index = 0 ; + m_r_map = 0 ; m_maxsizeSP = 5000 ; m_maxOneSize= 5 ; m_SP = 0 ; @@ -107,8 +107,8 @@ InDet::SiSpacePointsSeedMaker_ATLxk::SiSpacePointsSeedMaker_ATLxk declareProperty("etaMin" ,m_etamin ); declareProperty("etaMax" ,m_etamax ); declareProperty("pTmin" ,m_ptmin ); - declareProperty("radMax" ,r_rmax ); - declareProperty("radStep" ,r_rstep ); + declareProperty("radMax" ,m_r_rmax ); + declareProperty("radStep" ,m_r_rstep ); declareProperty("maxSize" ,m_maxsize ); declareProperty("maxSizeSP" ,m_maxsizeSP ); declareProperty("minZ" ,m_zmin ); @@ -150,20 +150,20 @@ InDet::SiSpacePointsSeedMaker_ATLxk::SiSpacePointsSeedMaker_ATLxk InDet::SiSpacePointsSeedMaker_ATLxk::~SiSpacePointsSeedMaker_ATLxk() { - if(r_index ) delete [] r_index ; - if(r_map ) delete [] r_map ; - if(r_Sorted) delete [] r_Sorted; + if(m_r_index ) delete [] m_r_index ; + if(m_r_map ) delete [] m_r_map ; + if(m_r_Sorted) delete [] m_r_Sorted; // Delete seeds // - for(i_seed=l_seeds.begin(); i_seed!=l_seeds.end (); ++i_seed) { - delete *i_seed; + for(m_i_seed=m_l_seeds.begin(); m_i_seed!=m_l_seeds.end (); ++m_i_seed) { + delete *m_i_seed; } // Delete space points for reconstruction // - i_spforseed=l_spforseed.begin(); - for(; i_spforseed!=l_spforseed.end(); ++i_spforseed) { - delete *i_spforseed; + m_i_spforseed=m_l_spforseed.begin(); + for(; m_i_spforseed!=m_l_spforseed.end(); ++m_i_spforseed) { + delete *m_i_spforseed; } if(m_seedOutput) delete m_seedOutput; @@ -261,22 +261,22 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::newEvent(int iteration) m_ipt2K = m_ipt2/(m_K*m_K); m_ipt2C = m_ipt2*m_COF ; m_COFK = m_COF*(m_K*m_K) ; - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); } else { - r_first = 0; fillLists(); return; + m_r_first = 0; fillLists(); return; } m_dbm ? m_checketa = false : m_checketa = m_dzdrmin > 1.; - float irstep = 1./r_rstep; - int irmax = r_size-1 ; - for(int i=0; i!=m_nr; ++i) {int n = r_index[i]; r_map[n] = 0; r_Sorted[n].clear();} + float irstep = 1./m_r_rstep; + int irmax = m_r_size-1 ; + for(int i=0; i!=m_nr; ++i) {int n = m_r_index[i]; m_r_map[n] = 0; m_r_Sorted[n].clear();} m_ns = m_nr = 0; // Get pixels space points containers from store gate // - r_first = 0; + m_r_first = 0; if(!m_dbm && m_pixel) { // m_spacepointsPixel = 0; @@ -294,7 +294,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::newEvent(int iteration) for(; sp != spe; ++sp) { - if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > r_rmax) continue; + if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > m_r_rmax) continue; // Remove DBM space points // @@ -305,14 +305,14 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::newEvent(int iteration) InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); if(!sps) continue; int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; - if(ir > r_first) r_first = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; + if(ir > m_r_first) m_r_first = ir; ++m_ns; } } } - ++r_first; + ++m_r_first; } // Get sct space points containers from store gate @@ -334,13 +334,13 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::newEvent(int iteration) for(; sp != spe; ++sp) { - if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > r_rmax) continue; + if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); if(!sps) continue; int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -359,13 +359,13 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::newEvent(int iteration) for (; sp!=spe; ++sp) { - if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > r_rmax) continue; + if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); if(!sps) continue; int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -395,21 +395,21 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::newEvent(int iteration) // const InDetDD::SiDetectorElement* de= static_cast<const InDetDD::SiDetectorElement*>((*sp)->clusterList().first->detectorElement()); - if(!de || !de->isDBM() || (*sp)->r() > r_rmax) continue; + if(!de || !de->isDBM() || (*sp)->r() > m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); if(!sps) continue; int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; - if(ir > r_first) r_first = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; + if(ir > m_r_first) m_r_first = ir; ++m_ns; } } } - ++r_first; + ++m_r_first; } - if(iteration < 0) r_first = 0; + if(iteration < 0) m_r_first = 0; fillLists(); } @@ -442,15 +442,15 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::newRegion m_ipt2C = m_ipt2*m_COF ; m_COFK = m_COF*(m_K*m_K) ; - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; - int irmax = r_size-1 ; + float irstep = 1./m_r_rstep; + int irmax = m_r_size-1 ; - r_first = 0 ; + m_r_first = 0 ; m_checketa = false ; - for(int i=0; i!=m_nr; ++i) {int n = r_index[i]; r_map[n] = 0; r_Sorted[n].clear();} + for(int i=0; i!=m_nr; ++i) {int n = m_r_index[i]; m_r_map[n] = 0; m_r_Sorted[n].clear();} m_ns = m_nr = 0; // Get pixels space points containers from store gate @@ -476,11 +476,11 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r > r_rmax) continue; + float r = (*sp)->r(); if(r > m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -510,11 +510,11 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r > r_rmax) continue; + float r = (*sp)->r(); if(r > m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -562,7 +562,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::find2Sp(const std::list<Trk::Vertex>& if(newv || !m_state || m_nspoint!=2 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); m_state = 1 ; m_nspoint = 2 ; m_nlist = 0 ; @@ -573,7 +573,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::find2Sp(const std::list<Trk::Vertex>& m_zMin = 0 ; production2Sp (); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -594,7 +594,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::find3Sp(const std::list<Trk::Vertex>& bool newv = newVertices(lv); if(newv || !m_state || m_nspoint!=3 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 3 ; m_nlist = 0 ; @@ -605,7 +605,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::find3Sp(const std::list<Trk::Vertex>& m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_seeds.begin(); if(m_outputlevel<=0) { @@ -627,7 +627,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::find3Sp(const std::list<Trk::Vertex>& bool newv = newVertices(lv); if(newv || !m_state || m_nspoint!=3 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 3 ; m_nlist = 0 ; @@ -638,7 +638,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::find3Sp(const std::list<Trk::Vertex>& m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_seeds.begin(); if(m_outputlevel<=0) { @@ -662,7 +662,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::findVSp (const std::list<Trk::Vertex>& if(newv || !m_state || m_nspoint!=4 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 4 ; m_nlist = 0 ; @@ -673,7 +673,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::findVSp (const std::list<Trk::Vertex>& m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_seeds.begin(); if(m_outputlevel<=0) { @@ -735,10 +735,10 @@ MsgStream& InDet::SiSpacePointsSeedMaker_ATLxk::dumpConditions( MsgStream& out ) <<std::setw(12)<<std::setprecision(5)<<m_rapcut <<" |"<<std::endl; out<<"| max radius SP | " - <<std::setw(12)<<std::setprecision(5)<<r_rmax + <<std::setw(12)<<std::setprecision(5)<<m_r_rmax <<" |"<<std::endl; out<<"| radius step | " - <<std::setw(12)<<std::setprecision(5)<<r_rstep + <<std::setw(12)<<std::setprecision(5)<<m_r_rstep <<" |"<<std::endl; out<<"| min Z-vertex position | " <<std::setw(12)<<std::setprecision(5)<<m_zmin @@ -847,7 +847,7 @@ MsgStream& InDet::SiSpacePointsSeedMaker_ATLxk::dumpEvent( MsgStream& out ) cons <<std::setw(12)<<m_nsazv <<" |"<<std::endl; out<<"| seeds | " - <<std::setw(12)<<l_seeds.size() + <<std::setw(12)<<m_l_seeds.size() <<" |"<<std::endl; out<<"|---------------------------------------------------------------------|" <<std::endl; @@ -922,13 +922,13 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::findNext () { if(m_endlist) return; - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); if (m_mode==0 || m_mode==1) production2Sp (); else if(m_mode==2 || m_mode==3) production3Sp (); else if(m_mode==5 || m_mode==6) production3Sp (); - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_seeds.begin(); ++m_nlist; } @@ -939,21 +939,21 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::findNext () bool InDet::SiSpacePointsSeedMaker_ATLxk::newVertices(const std::list<Trk::Vertex>& lV) { - unsigned int s1 = l_vertex.size(); + unsigned int s1 = m_l_vertex.size(); unsigned int s2 = lV .size(); m_isvertex = false; if(s1==0 && s2==0) return false; std::list<Trk::Vertex>::const_iterator v; - l_vertex.clear(); + m_l_vertex.clear(); if(s2 == 0) return false; m_isvertex = true; - for(v=lV.begin(); v!=lV.end(); ++v) {l_vertex.insert(float((*v).position().z()));} + for(v=lV.begin(); v!=lV.end(); ++v) {m_l_vertex.insert(float((*v).position().z()));} - m_zminU = (*l_vertex. begin())-20.; if( m_zminU < m_zmin) m_zminU = m_zmin; - m_zmaxU = (*l_vertex.rbegin())+20.; if( m_zmaxU > m_zmax) m_zmaxU = m_zmax; + m_zminU = (*m_l_vertex. begin())-20.; if( m_zminU < m_zmin) m_zminU = m_zmin; + m_zmaxU = (*m_l_vertex.rbegin())+20.; if( m_zmaxU > m_zmax) m_zmaxU = m_zmax; return false; } @@ -968,7 +968,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::buildFrameWork() // For DBM reconstruction we use new parameters // - if ( m_dbm ) {r_rmax = 150.; r_rstep = .5; m_drmin = 2.; m_ptmin = 10.; m_diver = 150.;} + if ( m_dbm ) {m_r_rmax = 150.; m_r_rstep = .5; m_drmin = 2.; m_ptmin = 10.; m_diver = 150.;} else if(m_ptmin < 100.) m_ptmin = 100.; if(m_diversss < m_diver ) m_diversss = m_diver ; @@ -978,7 +978,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::buildFrameWork() m_dzdrmax0 = 1./tan(2.*atan(exp(-m_etamax))); m_dzdrmin0 = 1./tan(2.*atan(exp(-m_etamin))); - m_r3max = r_rmax ; + m_r3max = m_r_rmax ; m_COF = 134*.05*9. ; m_ipt = 1./fabs(.9*m_ptmin) ; m_ipt2 = m_ipt*m_ipt ; @@ -988,27 +988,27 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::buildFrameWork() // Build radius sorted containers // - r_size = int((r_rmax+.1)/r_rstep); - r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[r_size]; - r_index = new int[r_size]; - r_map = new int[r_size]; - m_nr = 0; for(int i=0; i!=r_size; ++i) {r_index[i]=0; r_map[i]=0;} + m_r_size = int((m_r_rmax+.1)/m_r_rstep); + m_r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[m_r_size]; + m_r_index = new int[m_r_size]; + m_r_map = new int[m_r_size]; + m_nr = 0; for(int i=0; i!=m_r_size; ++i) {m_r_index[i]=0; m_r_map[i]=0;} // Build radius-azimuthal sorted containers // const float pi2 = 2.*M_PI ; const int NFmax = 53 ; const float sFmax = float(NFmax )/pi2; - const float m_sFmin = 100./60. ; + const float sFmin = 100./60. ; float ptm = 400.; if(!m_dbm && m_ptmin < ptm) ptm = m_ptmin; - m_sF = ptm /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < m_sFmin) m_sF = m_sFmin; + m_sF = ptm /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < sFmin) m_sF = sFmin; m_fNmax = int(pi2*m_sF); if(m_fNmax >=NFmax) m_fNmax = NFmax-1; // Build radius-azimuthal-Z sorted containers // - m_nrfz = 0; for(int i=0; i!=583; ++i) {rfz_index [i]=0; rfz_map [i]=0;} + m_nrfz = 0; for(int i=0; i!=583; ++i) {m_rfz_index [i]=0; m_rfz_map [i]=0;} // Build radius-azimuthal-Z sorted containers for Z-vertices // @@ -1016,7 +1016,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::buildFrameWork() const float sFvmax = float(NFtmax)/pi2; m_sFv = m_ptmin/120. ; if(m_sFv >sFvmax) m_sFv = sFvmax; m_fvNmax = int(pi2*m_sFv); if(m_fvNmax>=NFtmax) m_fvNmax = NFtmax-1; - m_nrfzv = 0; for(int i=0; i!=300; ++i) {rfzv_index[i]=0; rfzv_map[i]=0;} + m_nrfzv = 0; for(int i=0; i!=300; ++i) {m_rfzv_index[i]=0; m_rfzv_map[i]=0;} // Build maps for radius-azimuthal-Z sorted collections // @@ -1032,62 +1032,62 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::buildFrameWork() int a = f *11+z; int b = fb*11+z; int c = ft*11+z; - rfz_b [a] = 3; rfz_t [a] = 3; - rfz_ib[a][0] = a; rfz_it[a][0] = a; - rfz_ib[a][1] = b; rfz_it[a][1] = b; - rfz_ib[a][2] = c; rfz_it[a][2] = c; + m_rfz_b [a] = 3; m_rfz_t [a] = 3; + m_rfz_ib[a][0] = a; m_rfz_it[a][0] = a; + m_rfz_ib[a][1] = b; m_rfz_it[a][1] = b; + m_rfz_ib[a][2] = c; m_rfz_it[a][2] = c; if (z==5) { - rfz_t [a] = 9 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; - rfz_it[a][6] = a-1; - rfz_it[a][7] = b-1; - rfz_it[a][8] = c-1; + m_rfz_t [a] = 9 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; + m_rfz_it[a][6] = a-1; + m_rfz_it[a][7] = b-1; + m_rfz_it[a][8] = c-1; } else if(z> 5) { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a-1; - rfz_ib[a][4] = b-1; - rfz_ib[a][5] = c-1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a-1; + m_rfz_ib[a][4] = b-1; + m_rfz_ib[a][5] = c-1; if(z<10) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; } } else { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a+1; - rfz_ib[a][4] = b+1; - rfz_ib[a][5] = c+1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a+1; + m_rfz_ib[a][4] = b+1; + m_rfz_ib[a][5] = c+1; if(z>0) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a-1; - rfz_it[a][4] = b-1; - rfz_it[a][5] = c-1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a-1; + m_rfz_it[a][4] = b-1; + m_rfz_it[a][5] = c-1; } } if (z==3) { - rfz_b[a] = 9; - rfz_ib[a][6] = a+2; - rfz_ib[a][7] = b+2; - rfz_ib[a][8] = c+2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a+2; + m_rfz_ib[a][7] = b+2; + m_rfz_ib[a][8] = c+2; } else if(z==7) { - rfz_b[a] = 9; - rfz_ib[a][6] = a-2; - rfz_ib[a][7] = b-2; - rfz_ib[a][8] = c-2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a-2; + m_rfz_ib[a][7] = b-2; + m_rfz_ib[a][8] = c-2; } } } @@ -1106,21 +1106,21 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::buildFrameWork() int a = f *3+z; int b = fb*3+z; int c = ft*3+z; - rfzv_n[a] = 3; - rfzv_i[a][0] = a; - rfzv_i[a][1] = b; - rfzv_i[a][2] = c; + m_rfzv_n[a] = 3; + m_rfzv_i[a][0] = a; + m_rfzv_i[a][1] = b; + m_rfzv_i[a][2] = c; if (z>1) { - rfzv_n[a] = 6; - rfzv_i[a][3] = a-1; - rfzv_i[a][4] = b-1; - rfzv_i[a][5] = c-1; + m_rfzv_n[a] = 6; + m_rfzv_i[a][3] = a-1; + m_rfzv_i[a][4] = b-1; + m_rfzv_i[a][5] = c-1; } else if(z<1) { - rfzv_n[a] = 6; - rfzv_i[a][3] = a+1; - rfzv_i[a][4] = b+1; - rfzv_i[a][5] = c+1; + m_rfzv_n[a] = 6; + m_rfzv_i[a][3] = a+1; + m_rfzv_i[a][4] = b+1; + m_rfzv_i[a][5] = c+1; } } } @@ -1136,8 +1136,8 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::buildFrameWork() if(!m_seedOutput) m_seedOutput = new InDet::SiSpacePointsSeed(); - i_seed = l_seeds.begin(); - i_seede = l_seeds.end (); + m_i_seed = m_l_seeds.begin(); + m_i_seede = m_l_seeds.end (); } /////////////////////////////////////////////////////////////////// @@ -1200,9 +1200,9 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::fillLists() int ir0 = 0; bool ibl = false; - for(int i=r_first; i!=r_size; ++i) { + for(int i=m_r_first; i!=m_r_size; ++i) { - if(!r_map[i]) continue; r = r_Sorted[i].begin(); re = r_Sorted[i].end(); + if(!m_r_map[i]) continue; r = m_r_Sorted[i].begin(); re = m_r_Sorted[i].end(); if(!ir0) ir0 = i; if( m_iteration) { @@ -1232,7 +1232,7 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::fillLists() } int n = f*11+z; ++m_nsaz; - rfz_Sorted[n].push_back(*r); if(!rfz_map[n]++) rfz_index[m_nrfz++] = n; + m_rfz_Sorted[n].push_back(*r); if(!m_rfz_map[n]++) m_rfz_index[m_nrfz++] = n; if(!m_iteration && (*r)->spacepoint->clusterList().second == 0 && z>=3 && z<=7) { z<=4 ? z=0 : z>=6 ? z=2 : z=1; @@ -1241,11 +1241,11 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::fillLists() // f = int(F*m_sFv); f<0 ? f+=m_fvNmax : f> m_fvNmax ? f-=m_fvNmax : f=f; n = f*3+z; ++m_nsazv; - rfzv_Sorted[n].push_back(*r); if(!rfzv_map[n]++) rfzv_index[m_nrfzv++] = n; + m_rfzv_Sorted[n].push_back(*r); if(!m_rfzv_map[n]++) m_rfzv_index[m_nrfzv++] = n; } } } - if(!m_sct && ir0 && float(ir0)*r_rstep < 43.) m_umax = -200.; + if(!m_sct && ir0 && float(ir0)*m_r_rstep < 43.) m_umax = -200.; m_state = 0; } @@ -1256,13 +1256,13 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::fillLists() void InDet::SiSpacePointsSeedMaker_ATLxk::erase() { for(int i=0; i!=m_nrfz; ++i) { - int n = rfz_index[i]; rfz_map[n] = 0; - rfz_Sorted[n].clear(); + int n = m_rfz_index[i]; m_rfz_map[n] = 0; + m_rfz_Sorted[n].clear(); } for(int i=0; i!=m_nrfzv; ++i) { - int n = rfzv_index[i]; rfzv_map[n] = 0; - rfzv_Sorted[n].clear(); + int n = m_rfzv_index[i]; m_rfzv_map[n] = 0; + m_rfzv_Sorted[n].clear(); } m_state = 0; m_nsaz = 0; @@ -1307,9 +1307,9 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::production2Sp() int z = 0; if(!m_endlist) z = m_zMin; for(; z!=3; ++z) { - int a = f*3+z; if(!rfzv_map[a]) continue; - r0 = rfzv_Sorted[a].begin(); - r0e = rfzv_Sorted[a].end (); + int a = f*3+z; if(!m_rfzv_map[a]) continue; + r0 = m_rfzv_Sorted[a].begin(); + r0e = m_rfzv_Sorted[a].end (); if(!m_endlist) {r0 = m_rMin; m_endlist = true;} @@ -1327,14 +1327,14 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::production2Sp() // Bottom links production // - int NB = rfzv_n[a]; + int NB = m_rfzv_n[a]; for(int i=0; i!=NB; ++i) { - int an = rfzv_i[a][i]; - if(!rfzv_map[an]) continue; + int an = m_rfzv_i[a][i]; + if(!m_rfzv_map[an]) continue; - r = rfzv_Sorted[an].begin(); - re = rfzv_Sorted[an].end (); + r = m_rfzv_Sorted[an].begin(); + re = m_rfzv_Sorted[an].end (); for(; r!=re; ++r) { @@ -1399,19 +1399,19 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::production3Sp() for(; z!=11; ++z) { - int a = f *11+ZI[z]; if(!rfz_map[a]) continue; + int a = f *11+ZI[z]; if(!m_rfz_map[a]) continue; int NB = 0, NT = 0; - for(int i=0; i!=rfz_b[a]; ++i) { + for(int i=0; i!=m_rfz_b[a]; ++i) { - int an = rfz_ib[a][i]; - if(!rfz_map[an]) continue; - rb [NB] = rfz_Sorted[an].begin(); rbe[NB++] = rfz_Sorted[an].end(); + int an = m_rfz_ib[a][i]; + if(!m_rfz_map[an]) continue; + rb [NB] = m_rfz_Sorted[an].begin(); rbe[NB++] = m_rfz_Sorted[an].end(); } - for(int i=0; i!=rfz_t[a]; ++i) { + for(int i=0; i!=m_rfz_t[a]; ++i) { - int an = rfz_it[a][i]; - if(!rfz_map[an]) continue; - rt [NT] = rfz_Sorted[an].begin(); rte[NT++] = rfz_Sorted[an].end(); + int an = m_rfz_it[a][i]; + if(!m_rfz_map[an]) continue; + rt [NT] = m_rfz_Sorted[an].begin(); rte[NT++] = m_rfz_Sorted[an].end(); } if(!m_trigger) production3Sp (rb,rbe,rt,rte,NB,NT,nseed); else production3SpTrigger(rb,rbe,rt,rte,NB,NT,nseed); @@ -1871,14 +1871,14 @@ void InDet::SiSpacePointsSeedMaker_ATLxk::fillSeeds () if(l!=lf && s->spacepoint0()->radius() < 43. && w > -200.) continue; if(!s->setQuality(w)) continue; - if(i_seede!=l_seeds.end()) { - s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + s = (*m_i_seede++); *s = *(*l).second; } else { s = new SiSpacePointsProSeed(*(*l).second); - l_seeds.push_back(s); - i_seede = l_seeds.end(); + m_l_seeds.push_back(s); + m_i_seede = m_l_seeds.end(); } if (s->spacepoint0()->spacepoint->clusterList().second) w-=3000.; diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_BeamGas.cxx b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_BeamGas.cxx index 2011508e0315194bb3d429a19bbf6e2aa3b44c5a..b6205781b3d8a95e3917d4d11267f7035f45c836 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_BeamGas.cxx +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_BeamGas.cxx @@ -59,11 +59,11 @@ InDet::SiSpacePointsSeedMaker_BeamGas::SiSpacePointsSeedMaker_BeamGas m_diver = 10. ; m_diverpps = 1.7 ; m_dazmax = .02 ; - r_rmax = 600. ; - r_rstep = 2. ; - r_Sorted = 0 ; - r_index = 0 ; - r_map = 0 ; + m_r_rmax = 600. ; + m_r_rstep = 2. ; + m_r_Sorted = 0 ; + m_r_index = 0 ; + m_r_map = 0 ; m_maxsizeSP = 1500 ; m_maxOneSize= 5 ; m_SP = 0 ; @@ -91,8 +91,8 @@ InDet::SiSpacePointsSeedMaker_BeamGas::SiSpacePointsSeedMaker_BeamGas declareProperty("usePixel" ,m_pixel ); declareProperty("useSCT" ,m_sct ); declareProperty("pTmin" ,m_ptmin ); - declareProperty("radMax" ,r_rmax ); - declareProperty("radStep" ,r_rstep ); + declareProperty("radMax" ,m_r_rmax ); + declareProperty("radStep" ,m_r_rstep ); declareProperty("maxSize" ,m_maxsize ); declareProperty("maxSizeSP" ,m_maxsizeSP ); declareProperty("minZ" ,m_zmin ); @@ -126,20 +126,20 @@ InDet::SiSpacePointsSeedMaker_BeamGas::SiSpacePointsSeedMaker_BeamGas InDet::SiSpacePointsSeedMaker_BeamGas::~SiSpacePointsSeedMaker_BeamGas() { - if(r_index ) delete [] r_index ; - if(r_map ) delete [] r_map ; - if(r_Sorted) delete [] r_Sorted; + if(m_r_index ) delete [] m_r_index ; + if(m_r_map ) delete [] m_r_map ; + if(m_r_Sorted) delete [] m_r_Sorted; // Delete seeds // - for(i_seed=l_seeds.begin(); i_seed!=l_seeds.end (); ++i_seed) { - delete *i_seed; + for(m_i_seed=m_l_seeds.begin(); m_i_seed!=m_l_seeds.end (); ++m_i_seed) { + delete *m_i_seed; } // Delete space points for reconstruction // - i_spforseed=l_spforseed.begin(); - for(; i_spforseed!=l_spforseed.end(); ++i_spforseed) { - delete *i_spforseed; + m_i_spforseed=m_l_spforseed.begin(); + for(; m_i_spforseed!=m_l_spforseed.end(); ++m_i_spforseed) { + delete *m_i_spforseed; } if(m_SP) delete [] m_SP; @@ -226,9 +226,9 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::newEvent (int) } else m_K = 2./(300.* 5. ); - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; + float irstep = 1./m_r_rstep; // Get pixels space points containers from store gate // @@ -249,13 +249,13 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; int ir = int(r*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -280,13 +280,13 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; int ir = int(r*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -305,13 +305,13 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::newEvent (int) for (; sp!=spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; int ir = int(r*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -338,9 +338,9 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::newRegion } else m_K = 2./(300.* 5. ); - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; + float irstep = 1./m_r_rstep; // Get pixels space points containers from store gate // @@ -365,13 +365,13 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; int ir = int(r*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -401,13 +401,13 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; int ir = int(r*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -438,7 +438,7 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::find2Sp(const std::list<Trk::Vertex> if(!m_state || m_nspoint!=2 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); m_state = 1 ; m_nspoint = 2 ; m_nlist = 0 ; @@ -448,7 +448,7 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::find2Sp(const std::list<Trk::Vertex> m_zMin = 0 ; production2Sp (); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -466,7 +466,7 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::find3Sp(const std::list<Trk::Vertex> if(!m_state || m_nspoint!=3 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 3 ; m_nlist = 0 ; @@ -476,7 +476,7 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::find3Sp(const std::list<Trk::Vertex> m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -499,7 +499,7 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::findVSp (const std::list<Trk::Vertex if(!m_state || m_nspoint!=4 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 4 ; m_nlist = 0 ; @@ -509,7 +509,7 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::findVSp (const std::list<Trk::Vertex m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -570,10 +570,10 @@ MsgStream& InDet::SiSpacePointsSeedMaker_BeamGas::dumpConditions( MsgStream& out <<std::setw(12)<<std::setprecision(5)<<m_rapcut <<" |"<<std::endl; out<<"| max radius SP | " - <<std::setw(12)<<std::setprecision(5)<<r_rmax + <<std::setw(12)<<std::setprecision(5)<<m_r_rmax <<" |"<<std::endl; out<<"| radius step | " - <<std::setw(12)<<std::setprecision(5)<<r_rstep + <<std::setw(12)<<std::setprecision(5)<<m_r_rstep <<" |"<<std::endl; out<<"| min Z-vertex position | " <<std::setw(12)<<std::setprecision(5)<<m_zmin @@ -664,7 +664,7 @@ MsgStream& InDet::SiSpacePointsSeedMaker_BeamGas::dumpEvent( MsgStream& out ) co <<std::setw(12)<<m_nsaz <<" |"<<std::endl; out<<"| seeds | " - <<std::setw(12)<<l_seeds.size() + <<std::setw(12)<<m_l_seeds.size() <<" |"<<std::endl; out<<"|---------------------------------------------------------------------|" <<std::endl; @@ -739,12 +739,12 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::findNext () { if(m_endlist) return; - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); if (m_mode==0 || m_mode==1) production2Sp (); else if(m_mode==2 || m_mode==3) production3Sp (); else if(m_mode==5 || m_mode==6) production3Sp (); - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); ++m_nlist; } @@ -767,27 +767,27 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::buildFrameWork() // Build radius sorted containers // - r_size = int((r_rmax+.1)/r_rstep); - r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[r_size]; - r_index = new int[r_size]; - r_map = new int[r_size]; - m_nr = 0; for(int i=0; i!=r_size; ++i) {r_index[i]=0; r_map[i]=0;} + m_r_size = int((m_r_rmax+.1)/m_r_rstep); + m_r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[m_r_size]; + m_r_index = new int[m_r_size]; + m_r_map = new int[m_r_size]; + m_nr = 0; for(int i=0; i!=m_r_size; ++i) {m_r_index[i]=0; m_r_map[i]=0;} // Build radius-azimuthal sorted containers // const float pi2 = 2.*M_PI ; const int NFmax = 53 ; const float sFmax = float(NFmax )/pi2; - const float m_sFmin = 100./60. ; + const float sFmin = 100./60. ; - m_sF = m_ptmin /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < m_sFmin) m_sF = m_sFmin; + m_sF = m_ptmin /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < sFmin) m_sF = sFmin; m_fNmax = int(pi2*m_sF); if(m_fNmax >=NFmax) m_fNmax = NFmax-1; - m_nrf = 0; for(int i=0; i!= 53; ++i) {rf_index [i]=0; rf_map [i]=0;} + m_nrf = 0; for(int i=0; i!= 53; ++i) {m_rf_index [i]=0; m_rf_map [i]=0;} // Build radius-azimuthal-Z sorted containers // - m_nrfz = 0; for(int i=0; i!=583; ++i) {rfz_index [i]=0; rfz_map [i]=0;} + m_nrfz = 0; for(int i=0; i!=583; ++i) {m_rfz_index [i]=0; m_rfz_map [i]=0;} // Build maps for radius-azimuthal-Z sorted collections // @@ -803,48 +803,48 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::buildFrameWork() int a = f *11+z; int b = fb*11+z; int c = ft*11+z; - rfz_b [a] = 3; rfz_t [a] = 3; - rfz_ib[a][0] = a; rfz_it[a][0] = a; - rfz_ib[a][1] = b; rfz_it[a][1] = b; - rfz_ib[a][2] = c; rfz_it[a][2] = c; + m_rfz_b [a] = 3; m_rfz_t [a] = 3; + m_rfz_ib[a][0] = a; m_rfz_it[a][0] = a; + m_rfz_ib[a][1] = b; m_rfz_it[a][1] = b; + m_rfz_ib[a][2] = c; m_rfz_it[a][2] = c; if (z==5) { - rfz_t [a] = 9 ; rfz_b [a] = 9 ; - rfz_it[a][3] = a+1; rfz_ib[a][3] = a+1; - rfz_it[a][4] = b+1; rfz_ib[a][4] = b+1; - rfz_it[a][5] = c+1; rfz_ib[a][5] = c+1; - rfz_it[a][6] = a-1; rfz_ib[a][6] = a-1; - rfz_it[a][7] = b-1; rfz_ib[a][7] = b-1; - rfz_it[a][8] = c-1; rfz_ib[a][8] = c-1; + m_rfz_t [a] = 9 ; m_rfz_b [a] = 9 ; + m_rfz_it[a][3] = a+1; m_rfz_ib[a][3] = a+1; + m_rfz_it[a][4] = b+1; m_rfz_ib[a][4] = b+1; + m_rfz_it[a][5] = c+1; m_rfz_ib[a][5] = c+1; + m_rfz_it[a][6] = a-1; m_rfz_ib[a][6] = a-1; + m_rfz_it[a][7] = b-1; m_rfz_ib[a][7] = b-1; + m_rfz_it[a][8] = c-1; m_rfz_ib[a][8] = c-1; } else if(z> 5) { - rfz_b [a] = 6 ; rfz_t [a] = 6 ; - rfz_ib[a][3] = a-1; rfz_it[a][3] = a-1; - rfz_ib[a][4] = b-1; rfz_it[a][4] = b-1; - rfz_ib[a][5] = c-1; rfz_it[a][5] = c-1; + m_rfz_b [a] = 6 ; m_rfz_t [a] = 6 ; + m_rfz_ib[a][3] = a-1; m_rfz_it[a][3] = a-1; + m_rfz_ib[a][4] = b-1; m_rfz_it[a][4] = b-1; + m_rfz_ib[a][5] = c-1; m_rfz_it[a][5] = c-1; if(z<10) { - rfz_t [a] = 9 ; rfz_b [a] = 9 ; - rfz_it[a][6] = a+1; rfz_ib[a][6] = a+1; - rfz_it[a][7] = b+1; rfz_ib[a][7] = b+1; - rfz_it[a][8] = c+1; rfz_ib[a][8] = c+1; + m_rfz_t [a] = 9 ; m_rfz_b [a] = 9 ; + m_rfz_it[a][6] = a+1; m_rfz_ib[a][6] = a+1; + m_rfz_it[a][7] = b+1; m_rfz_ib[a][7] = b+1; + m_rfz_it[a][8] = c+1; m_rfz_ib[a][8] = c+1; } } else { - rfz_b [a] = 6 ; rfz_t [a] = 6 ; - rfz_ib[a][3] = a+1; rfz_it[a][3] = a+1; - rfz_ib[a][4] = b+1; rfz_it[a][4] = b+1; - rfz_ib[a][5] = c+1; rfz_it[a][5] = c+1; + m_rfz_b [a] = 6 ; m_rfz_t [a] = 6 ; + m_rfz_ib[a][3] = a+1; m_rfz_it[a][3] = a+1; + m_rfz_ib[a][4] = b+1; m_rfz_it[a][4] = b+1; + m_rfz_ib[a][5] = c+1; m_rfz_it[a][5] = c+1; if(z>0) { - rfz_t [a] = 9 ; rfz_b [a] = 9 ; - rfz_it[a][6] = a-1; rfz_ib[a][6] = a-1; - rfz_it[a][7] = b-1; rfz_ib[a][7] = b-1; - rfz_it[a][8] = c-1; rfz_ib[a][8] = c-1; + m_rfz_t [a] = 9 ; m_rfz_b [a] = 9 ; + m_rfz_it[a][6] = a-1; m_rfz_ib[a][6] = a-1; + m_rfz_it[a][7] = b-1; m_rfz_ib[a][7] = b-1; + m_rfz_it[a][8] = c-1; m_rfz_ib[a][8] = c-1; } } } @@ -859,8 +859,8 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::buildFrameWork() if(!m_Zo) m_Zo = new float[m_maxsizeSP]; if(!m_OneSeeds) m_OneSeeds = new InDet::SiSpacePointsSeed [m_maxOneSize]; - i_seed = l_seeds.begin(); - i_seede = l_seeds.end (); + m_i_seed = m_l_seeds.begin(); + m_i_seede = m_l_seeds.end (); } /////////////////////////////////////////////////////////////////// @@ -923,18 +923,18 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::fillLists() const float pi2 = 2.*M_PI; std::list<InDet::SiSpacePointForSeed*>::iterator r; - for(int i=0; i!= r_size; ++i) { + for(int i=0; i!= m_r_size; ++i) { - if(!r_map[i]) continue; r = r_Sorted[i].begin(); + if(!m_r_map[i]) continue; r = m_r_Sorted[i].begin(); - while(r!=r_Sorted[i].end()) { + while(r!=m_r_Sorted[i].end()) { // Azimuthal angle sort // float F = (*r)->phi(); if(F<0.) F+=pi2; int f = int(F*m_sF); f<0 ? f = m_fNmax : f>m_fNmax ? f = 0 : f=f; - rf_Sorted[f].push_back(*r); if(!rf_map[f]++) rf_index[m_nrf++] = f; + m_rf_Sorted[f].push_back(*r); if(!m_rf_map[f]++) m_rf_index[m_nrf++] = f; int z; float Z = (*r)->z(); @@ -947,10 +947,10 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::fillLists() Z>-250.?z=5:Z>-450.?z=4:Z>-925.?z=3:Z>-1400.?z=2:Z>-2500.?z=1:z= 0; } int n = f*11+z; ++m_nsaz; - rfz_Sorted[n].push_back(*r); if(!rfz_map[n]++) rfz_index[m_nrfz++] = n; - r_Sorted[i].erase(r++); + m_rfz_Sorted[n].push_back(*r); if(!m_rfz_map[n]++) m_rfz_index[m_nrfz++] = n; + m_r_Sorted[i].erase(r++); } - r_map[i] = 0; + m_r_map[i] = 0; } m_nr = 0; m_state = 0; @@ -963,18 +963,18 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::fillLists() void InDet::SiSpacePointsSeedMaker_BeamGas::erase() { for(int i=0; i!=m_nr; ++i) { - int n = r_index[i]; r_map[n] = 0; - r_Sorted[n].erase(r_Sorted[n].begin(),r_Sorted[n].end()); + int n = m_r_index[i]; m_r_map[n] = 0; + m_r_Sorted[n].erase(m_r_Sorted[n].begin(),m_r_Sorted[n].end()); } for(int i=0; i!=m_nrf; ++i) { - int n = rf_index[i]; rf_map[n] = 0; - rf_Sorted[n].erase(rf_Sorted[n].begin(),rf_Sorted[n].end()); + int n = m_rf_index[i]; m_rf_map[n] = 0; + m_rf_Sorted[n].erase(m_rf_Sorted[n].begin(),m_rf_Sorted[n].end()); } for(int i=0; i!=m_nrfz; ++i) { - int n = rfz_index[i]; rfz_map[n] = 0; - rfz_Sorted[n].erase(rfz_Sorted[n].begin(),rfz_Sorted[n].end()); + int n = m_rfz_index[i]; m_rfz_map[n] = 0; + m_rfz_Sorted[n].erase(m_rfz_Sorted[n].begin(),m_rfz_Sorted[n].end()); } m_state = 0; @@ -1016,19 +1016,19 @@ void InDet::SiSpacePointsSeedMaker_BeamGas::production3Sp() for(; z!=11; ++z) { - int a = f *11+ZI[z]; if(!rfz_map[a]) continue; + int a = f *11+ZI[z]; if(!m_rfz_map[a]) continue; int NB = 0, NT = 0; - for(int i=0; i!=rfz_b[a]; ++i) { + for(int i=0; i!=m_rfz_b[a]; ++i) { - int an = rfz_ib[a][i]; - if(!rfz_map[an]) continue; - rb [NB] = rfz_Sorted[an].begin(); rbe[NB++] = rfz_Sorted[an].end(); + int an = m_rfz_ib[a][i]; + if(!m_rfz_map[an]) continue; + rb [NB] = m_rfz_Sorted[an].begin(); rbe[NB++] = m_rfz_Sorted[an].end(); } - for(int i=0; i!=rfz_t[a]; ++i) { + for(int i=0; i!=m_rfz_t[a]; ++i) { - int an = rfz_it[a][i]; - if(!rfz_map[an]) continue; - rt [NT] = rfz_Sorted[an].begin(); rte[NT++] = rfz_Sorted[an].end(); + int an = m_rfz_it[a][i]; + if(!m_rfz_map[an]) continue; + rt [NT] = m_rfz_Sorted[an].begin(); rte[NT++] = m_rfz_Sorted[an].end(); } production3Sp(rb,rbe,rt,rte,NB,NT,nseed); if(!m_endlist) {m_fNmin=f; m_zMin = z; return;} diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_Cosmic.cxx b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_Cosmic.cxx index fbf6800ab1cafb1c93e00451f6199d0f75b87a43..896a2844f6b32e91e5e59120773a536be6949c7d 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_Cosmic.cxx +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_Cosmic.cxx @@ -49,11 +49,11 @@ InDet::SiSpacePointsSeedMaker_Cosmic::SiSpacePointsSeedMaker_Cosmic m_zmin = -10000. ; m_zmax = +10000. ; m_diver = 10. ; - r_rmax = 600. ; - r_rstep = 5. ; - r_Sorted = 0 ; - r_index = 0 ; - r_map = 0 ; + m_r_rmax = 600. ; + m_r_rstep = 5. ; + m_r_Sorted = 0 ; + m_r_index = 0 ; + m_r_map = 0 ; m_maxsizeSP = 1500 ; m_SP = 0 ; m_R = 0 ; @@ -74,8 +74,8 @@ InDet::SiSpacePointsSeedMaker_Cosmic::SiSpacePointsSeedMaker_Cosmic declareProperty("usePixel" ,m_pixel ); declareProperty("useSCT" ,m_sct ); declareProperty("pTmin" ,m_ptmin ); - declareProperty("radMax" ,r_rmax ); - declareProperty("radStep" ,r_rstep ); + declareProperty("radMax" ,m_r_rmax ); + declareProperty("radStep" ,m_r_rstep ); declareProperty("maxSize" ,m_maxsize ); declareProperty("maxSizeSP" ,m_maxsizeSP ); declareProperty("minZ" ,m_zmin ); @@ -103,15 +103,15 @@ InDet::SiSpacePointsSeedMaker_Cosmic::SiSpacePointsSeedMaker_Cosmic InDet::SiSpacePointsSeedMaker_Cosmic::~SiSpacePointsSeedMaker_Cosmic() { - if(r_index ) delete [] r_index ; - if(r_map ) delete [] r_map ; - if(r_Sorted) delete [] r_Sorted; + if(m_r_index ) delete [] m_r_index ; + if(m_r_map ) delete [] m_r_map ; + if(m_r_Sorted) delete [] m_r_Sorted; // Delete space points for reconstruction // - i_spforseed=l_spforseed.begin(); - for(; i_spforseed!=l_spforseed.end(); ++i_spforseed) { - delete *i_spforseed; + m_i_spforseed=m_l_spforseed.begin(); + for(; m_i_spforseed!=m_l_spforseed.end(); ++m_i_spforseed) { + delete *m_i_spforseed; } if(m_SP) delete [] m_SP; @@ -181,10 +181,11 @@ StatusCode InDet::SiSpacePointsSeedMaker_Cosmic::finalize() void InDet::SiSpacePointsSeedMaker_Cosmic::newEvent (int) { - if(!m_pixel && !m_sct) return; erase(); - i_spforseed = l_spforseed.begin(); + if(!m_pixel && !m_sct) return; + erase(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; + float irstep = 1./m_r_rstep; float errorsc[4] = {16.,16.,100.,16.}; // Get pixels space points containers from store gate @@ -206,14 +207,14 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; - int ir = int(((*sp)->globalPosition().y()+r_rmax)*irstep); + int ir = int(((*sp)->globalPosition().y()+m_r_rmax)*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp),errorsc); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -238,14 +239,14 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; - int ir = int(((*sp)->globalPosition().y()+r_rmax)*irstep); + int ir = int(((*sp)->globalPosition().y()+m_r_rmax)*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp),errorsc); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -264,14 +265,14 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::newEvent (int) for (; sp!=spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; - int ir = int(((*sp)->globalPosition().y()+r_rmax)*irstep); + int ir = int(((*sp)->globalPosition().y()+m_r_rmax)*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp),errorsc); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -289,9 +290,9 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::newRegion { if(!m_pixel && !m_sct) return; erase(); - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; + float irstep = 1./m_r_rstep; float errorsc[4] = {16.,16.,100.,16.}; // Get pixels space points containers from store gate @@ -317,13 +318,13 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; - int ir = int(((*sp)->globalPosition().y()+r_rmax)*irstep); + int ir = int(((*sp)->globalPosition().y()+m_r_rmax)*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp),errorsc); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -353,13 +354,13 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; - int ir = int(((*sp)->globalPosition().y()+r_rmax)*irstep); + int ir = int(((*sp)->globalPosition().y()+m_r_rmax)*irstep); InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp),errorsc); - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -387,7 +388,7 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::find2Sp(const std::list<Trk::Vertex>& { int mode; lv.begin()!=lv.end() ? mode = 1 : mode = 0; - m_nseeds = 0; l_seeds.erase(l_seeds.begin(),l_seeds.end()); + m_nseeds = 0; m_l_seeds.erase(m_l_seeds.begin(),m_l_seeds.end()); if( !m_state || m_nspoint!=2 || m_mode!=mode || m_nlist) { @@ -401,8 +402,8 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::find2Sp(const std::list<Trk::Vertex>& production2Sp (); } - i_seed = l_seeds.begin(); - i_seede = l_seeds.end (); + m_i_seed = m_l_seeds.begin(); + m_i_seede = m_l_seeds.end (); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -418,7 +419,7 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::find3Sp(const std::list<Trk::Vertex>& { int mode; lv.begin()!=lv.end() ? mode = 3 : mode = 2; - m_nseeds = 0; l_seeds.erase(l_seeds.begin(),l_seeds.end()); + m_nseeds = 0; m_l_seeds.erase(m_l_seeds.begin(),m_l_seeds.end()); if(!m_state || m_nspoint!=3 || m_mode!=mode || m_nlist) { @@ -432,8 +433,8 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::find3Sp(const std::list<Trk::Vertex>& production3Sp(); } - i_seed = l_seeds.begin(); - i_seede = l_seeds.end (); + m_i_seed = m_l_seeds.begin(); + m_i_seede = m_l_seeds.end (); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -467,8 +468,8 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::findVSp (const std::list<Trk::Vertex> production3Sp(); } - i_seed = l_seeds.begin(); - i_seede = l_seeds.end (); + m_i_seed = m_l_seeds.begin(); + m_i_seede = m_l_seeds.end (); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -481,7 +482,8 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::findVSp (const std::list<Trk::Vertex> MsgStream& InDet::SiSpacePointsSeedMaker_Cosmic::dump( MsgStream& out ) const { - if(m_nprint) return dumpEvent(out); return dumpConditions(out); + if(m_nprint) return dumpEvent(out); + return dumpConditions(out); } /////////////////////////////////////////////////////////////////// @@ -525,10 +527,10 @@ MsgStream& InDet::SiSpacePointsSeedMaker_Cosmic::dumpConditions( MsgStream& out <<std::setw(12)<<std::setprecision(5)<<m_rapcut <<" |"<<std::endl; out<<"| max radius SP | " - <<std::setw(12)<<std::setprecision(5)<<r_rmax + <<std::setw(12)<<std::setprecision(5)<<m_r_rmax <<" |"<<std::endl; out<<"| radius step | " - <<std::setw(12)<<std::setprecision(5)<<r_rstep + <<std::setw(12)<<std::setprecision(5)<<m_r_rstep <<" |"<<std::endl; out<<"| min Z-vertex position | " <<std::setw(12)<<std::setprecision(5)<<m_zmin @@ -566,7 +568,7 @@ MsgStream& InDet::SiSpacePointsSeedMaker_Cosmic::dumpEvent( MsgStream& out ) con <<std::setw(12)<<m_nsaz <<" |"<<std::endl; out<<"| seeds | " - <<std::setw(12)<<l_seeds.size() + <<std::setw(12)<<m_l_seeds.size() <<" |"<<std::endl; out<<"|---------------------------------------------------------------------|" <<std::endl; @@ -647,27 +649,27 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::buildFrameWork() // Build radius sorted containers // - r_size = int((r_rmax+.1)/r_rstep)*2; - r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[r_size]; - r_index = new int[r_size]; - r_map = new int[r_size]; - m_nr = 0; for(int i=0; i!=r_size; ++i) {r_index[i]=0; r_map[i]=0;} + m_r_size = int((m_r_rmax+.1)/m_r_rstep)*2; + m_r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[m_r_size]; + m_r_index = new int[m_r_size]; + m_r_map = new int[m_r_size]; + m_nr = 0; for(int i=0; i!=m_r_size; ++i) {m_r_index[i]=0; m_r_map[i]=0;} // Build radius-azimuthal sorted containers // const float pi2 = 2.*M_PI ; const int NFmax = 53 ; const float sFmax = float(NFmax )/pi2; - const float m_sFmin = 50./60. ; + const float sFmin = 50./60. ; - m_sF = 50./60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < m_sFmin) m_sF = m_sFmin; + m_sF = 50./60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < sFmin) m_sF = sFmin; m_fNmax = int(pi2*m_sF); if(m_fNmax >=NFmax) m_fNmax = NFmax-1; - m_nrf = 0; for(int i=0; i!= 53; ++i) {rf_index [i]=0; rf_map [i]=0;} + m_nrf = 0; for(int i=0; i!= 53; ++i) {m_rf_index [i]=0; m_rf_map [i]=0;} // Build radius-azimuthal-Z sorted containers // - m_nrfz = 0; for(int i=0; i!=583; ++i) {rfz_index [i]=0; rfz_map [i]=0;} + m_nrfz = 0; for(int i=0; i!=583; ++i) {m_rfz_index [i]=0; m_rfz_map [i]=0;} // Build maps for radius-azimuthal-Z sorted collections // @@ -683,62 +685,62 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::buildFrameWork() int a = f *11+z; int b = fb*11+z; int c = ft*11+z; - rfz_b [a] = 3; rfz_t [a] = 3; - rfz_ib[a][0] = a; rfz_it[a][0] = a; - rfz_ib[a][1] = b; rfz_it[a][1] = b; - rfz_ib[a][2] = c; rfz_it[a][2] = c; + m_rfz_b [a] = 3; m_rfz_t [a] = 3; + m_rfz_ib[a][0] = a; m_rfz_it[a][0] = a; + m_rfz_ib[a][1] = b; m_rfz_it[a][1] = b; + m_rfz_ib[a][2] = c; m_rfz_it[a][2] = c; if (z==5) { - rfz_t [a] = 9 ; rfz_b [a] = 9 ; - rfz_it[a][3] = a+1; rfz_ib[a][3] = a+1; - rfz_it[a][4] = b+1; rfz_ib[a][4] = b+1; - rfz_it[a][5] = c+1; rfz_ib[a][5] = c+1; - rfz_it[a][6] = a-1; rfz_ib[a][6] = a-1; - rfz_it[a][7] = b-1; rfz_ib[a][7] = b-1; - rfz_it[a][8] = c-1; rfz_ib[a][8] = c-1; + m_rfz_t [a] = 9 ; m_rfz_b [a] = 9 ; + m_rfz_it[a][3] = a+1; m_rfz_ib[a][3] = a+1; + m_rfz_it[a][4] = b+1; m_rfz_ib[a][4] = b+1; + m_rfz_it[a][5] = c+1; m_rfz_ib[a][5] = c+1; + m_rfz_it[a][6] = a-1; m_rfz_ib[a][6] = a-1; + m_rfz_it[a][7] = b-1; m_rfz_ib[a][7] = b-1; + m_rfz_it[a][8] = c-1; m_rfz_ib[a][8] = c-1; } else if(z> 5) { - rfz_b [a] = 6 ; rfz_t [a] = 6 ; - rfz_ib[a][3] = a-1; rfz_it[a][3] = a-1; - rfz_ib[a][4] = b-1; rfz_it[a][4] = b-1; - rfz_ib[a][5] = c-1; rfz_it[a][5] = c-1; + m_rfz_b [a] = 6 ; m_rfz_t [a] = 6 ; + m_rfz_ib[a][3] = a-1; m_rfz_it[a][3] = a-1; + m_rfz_ib[a][4] = b-1; m_rfz_it[a][4] = b-1; + m_rfz_ib[a][5] = c-1; m_rfz_it[a][5] = c-1; if(z<10) { - rfz_t [a] = 9 ; rfz_b [a] = 9 ; - rfz_it[a][6] = a+1; rfz_ib[a][6] = a+1; - rfz_it[a][7] = b+1; rfz_ib[a][7] = b+1; - rfz_it[a][8] = c+1; rfz_ib[a][8] = c+1; + m_rfz_t [a] = 9 ; m_rfz_b [a] = 9 ; + m_rfz_it[a][6] = a+1; m_rfz_ib[a][6] = a+1; + m_rfz_it[a][7] = b+1; m_rfz_ib[a][7] = b+1; + m_rfz_it[a][8] = c+1; m_rfz_ib[a][8] = c+1; } } else { - rfz_b [a] = 6 ; rfz_t [a] = 6 ; - rfz_ib[a][3] = a+1; rfz_it[a][3] = a+1; - rfz_ib[a][4] = b+1; rfz_it[a][4] = b+1; - rfz_ib[a][5] = c+1; rfz_it[a][5] = c+1; + m_rfz_b [a] = 6 ; m_rfz_t [a] = 6 ; + m_rfz_ib[a][3] = a+1; m_rfz_it[a][3] = a+1; + m_rfz_ib[a][4] = b+1; m_rfz_it[a][4] = b+1; + m_rfz_ib[a][5] = c+1; m_rfz_it[a][5] = c+1; if(z>0) { - rfz_t [a] = 9 ; rfz_b [a] = 9 ; - rfz_it[a][6] = a-1; rfz_ib[a][6] = a-1; - rfz_it[a][7] = b-1; rfz_ib[a][7] = b-1; - rfz_it[a][8] = c-1; rfz_ib[a][8] = c-1; + m_rfz_t [a] = 9 ; m_rfz_b [a] = 9 ; + m_rfz_it[a][6] = a-1; m_rfz_ib[a][6] = a-1; + m_rfz_it[a][7] = b-1; m_rfz_ib[a][7] = b-1; + m_rfz_it[a][8] = c-1; m_rfz_ib[a][8] = c-1; } } if (z==3) { - rfz_b[a] = 12; rfz_t [a] = 12; - rfz_ib[a][ 9] = a+2; rfz_it[a][ 9] = a+2; - rfz_ib[a][10] = b+2; rfz_it[a][10] = b+2; - rfz_ib[a][11] = c+2; rfz_it[a][11] = c+2; + m_rfz_b[a] = 12; m_rfz_t [a] = 12; + m_rfz_ib[a][ 9] = a+2; m_rfz_it[a][ 9] = a+2; + m_rfz_ib[a][10] = b+2; m_rfz_it[a][10] = b+2; + m_rfz_ib[a][11] = c+2; m_rfz_it[a][11] = c+2; } else if(z==7) { - rfz_b[a] = 12; rfz_t [a] = 12; - rfz_ib[a][ 9] = a-2; rfz_it[a][ 9] = a-2; - rfz_ib[a][10] = b-2; rfz_it[a][10] = b-2; - rfz_ib[a][11] = c-2; rfz_it[a][11] = c-2; + m_rfz_b[a] = 12; m_rfz_t [a] = 12; + m_rfz_ib[a][ 9] = a-2; m_rfz_it[a][ 9] = a-2; + m_rfz_ib[a][10] = b-2; m_rfz_it[a][10] = b-2; + m_rfz_ib[a][11] = c-2; m_rfz_it[a][11] = c-2; } } } @@ -761,11 +763,11 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::fillLists() const float pi2 = 2.*M_PI; std::list<InDet::SiSpacePointForSeed*>::iterator r; - for(int i=0; i!= r_size; ++i) { + for(int i=0; i!= m_r_size; ++i) { - if(!r_map[i]) continue; r = r_Sorted[i].begin(); + if(!m_r_map[i]) continue; r = m_r_Sorted[i].begin(); - while(r!=r_Sorted[i].end()) { + while(r!=m_r_Sorted[i].end()) { // Azimuthal angle sort // @@ -773,7 +775,7 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::fillLists() int f = int(F*m_sF); f<0 ? f = m_fNmax : f>m_fNmax ? f = 0 : f=f; f = 1; - rf_Sorted[f].push_back(*r); if(!rf_map[f]++) rf_index[m_nrf++] = f; + m_rf_Sorted[f].push_back(*r); if(!m_rf_map[f]++) m_rf_index[m_nrf++] = f; int z; float Z = (*r)->z(); @@ -786,10 +788,10 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::fillLists() Z>-250.?z=5:Z>-450.?z=4:Z>-925.?z=3:Z>-1400.?z=2:Z>-2500.?z=1:z= 0; } int n = f*11+z; ++m_nsaz; - rfz_Sorted[n].push_back(*r); if(!rfz_map[n]++) rfz_index[m_nrfz++] = n; - r_Sorted[i].erase(r++); + m_rfz_Sorted[n].push_back(*r); if(!m_rfz_map[n]++) m_rfz_index[m_nrfz++] = n; + m_r_Sorted[i].erase(r++); } - r_map[i] = 0; + m_r_map[i] = 0; } m_nr = 0; m_state = 0; @@ -802,18 +804,18 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::fillLists() void InDet::SiSpacePointsSeedMaker_Cosmic::erase() { for(int i=0; i!=m_nr; ++i) { - int n = r_index[i]; r_map[n] = 0; - r_Sorted[n].erase(r_Sorted[n].begin(),r_Sorted[n].end()); + int n = m_r_index[i]; m_r_map[n] = 0; + m_r_Sorted[n].erase(m_r_Sorted[n].begin(),m_r_Sorted[n].end()); } for(int i=0; i!=m_nrf; ++i) { - int n = rf_index[i]; rf_map[n] = 0; - rf_Sorted[n].erase(rf_Sorted[n].begin(),rf_Sorted[n].end()); + int n = m_rf_index[i]; m_rf_map[n] = 0; + m_rf_Sorted[n].erase(m_rf_Sorted[n].begin(),m_rf_Sorted[n].end()); } for(int i=0; i!=m_nrfz; ++i) { - int n = rfz_index[i]; rfz_map[n] = 0; - rfz_Sorted[n].erase(rfz_Sorted[n].begin(),rfz_Sorted[n].end()); + int n = m_rfz_index[i]; m_rfz_map[n] = 0; + m_rfz_Sorted[n].erase(m_rfz_Sorted[n].begin(),m_rfz_Sorted[n].end()); } m_state = 0; @@ -865,19 +867,19 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::production3Sp() for(; z!=11; ++z) { - int a = f *11+ZI[z]; if(!rfz_map[a]) continue; + int a = f *11+ZI[z]; if(!m_rfz_map[a]) continue; int NB = 0, NT = 0; - for(int i=0; i!=rfz_b[a]; ++i) { + for(int i=0; i!=m_rfz_b[a]; ++i) { - int an = rfz_ib[a][i]; - if(!rfz_map[an]) continue; - rb [NB] = rfz_Sorted[an].begin(); rbe[NB++] = rfz_Sorted[an].end(); + int an = m_rfz_ib[a][i]; + if(!m_rfz_map[an]) continue; + rb [NB] = m_rfz_Sorted[an].begin(); rbe[NB++] = m_rfz_Sorted[an].end(); } - for(int i=0; i!=rfz_t[a]; ++i) { + for(int i=0; i!=m_rfz_t[a]; ++i) { - int an = rfz_it[a][i]; - if(!rfz_map[an]) continue; - rt [NT] = rfz_Sorted[an].begin(); rte[NT++] = rfz_Sorted[an].end(); + int an = m_rfz_it[a][i]; + if(!m_rfz_map[an]) continue; + rt [NT] = m_rfz_Sorted[an].begin(); rte[NT++] = m_rfz_Sorted[an].end(); } production3Sp(rb,rbe,rt,rte,NB,NT,K,ipt); if(!m_endlist) {m_fNmin=f; m_zMin = z; return;} @@ -907,19 +909,19 @@ void InDet::SiSpacePointsSeedMaker_Cosmic::production3SpWithoutField() for(; z!=11; ++z) { - int a = f *11+ZI[z]; if(!rfz_map[a]) continue; + int a = f *11+ZI[z]; if(!m_rfz_map[a]) continue; int NB = 0, NT = 0; - for(int i=0; i!=rfz_b[a]; ++i) { + for(int i=0; i!=m_rfz_b[a]; ++i) { - int an = rfz_ib[a][i]; - if(!rfz_map[an]) continue; - rb [NB] = rfz_Sorted[an].begin(); rbe[NB++] = rfz_Sorted[an].end(); + int an = m_rfz_ib[a][i]; + if(!m_rfz_map[an]) continue; + rb [NB] = m_rfz_Sorted[an].begin(); rbe[NB++] = m_rfz_Sorted[an].end(); } - for(int i=0; i!=rfz_t[a]; ++i) { + for(int i=0; i!=m_rfz_t[a]; ++i) { - int an = rfz_it[a][i]; - if(!rfz_map[an]) continue; - rt [NT] = rfz_Sorted[an].begin(); rte[NT++] = rfz_Sorted[an].end(); + int an = m_rfz_it[a][i]; + if(!m_rfz_map[an]) continue; + rt [NT] = m_rfz_Sorted[an].begin(); rte[NT++] = m_rfz_Sorted[an].end(); } production3SpWithoutField(rb,rbe,rt,rte,NB,NT,ipt); if(!m_endlist) {m_fNmin=f; m_zMin = z; return;} diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_HeavyIon.cxx b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_HeavyIon.cxx index 056fb3772bd49a909003ca2a792e7dc4824b7062..78db556995b6039289d3a7794eb8e2aef2463a40 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_HeavyIon.cxx +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_HeavyIon.cxx @@ -59,11 +59,11 @@ InDet::SiSpacePointsSeedMaker_HeavyIon::SiSpacePointsSeedMaker_HeavyIon m_diverpps = 1.7 ; m_diversss = 1000 ; m_dazmax = .02 ; - r_rmax = 600. ; - r_rstep = 2. ; - r_Sorted = 0 ; - r_index = 0 ; - r_map = 0 ; + m_r_rmax = 600. ; + m_r_rstep = 2. ; + m_r_Sorted = 0 ; + m_r_index = 0 ; + m_r_map = 0 ; m_maxsizeSP = 4000 ; m_maxOneSize= 6 ; m_SP = 0 ; @@ -91,8 +91,8 @@ InDet::SiSpacePointsSeedMaker_HeavyIon::SiSpacePointsSeedMaker_HeavyIon declareProperty("usePixel" ,m_pixel ); declareProperty("useSCT" ,m_sct ); declareProperty("pTmin" ,m_ptmin ); - declareProperty("radMax" ,r_rmax ); - declareProperty("radStep" ,r_rstep ); + declareProperty("radMax" ,m_r_rmax ); + declareProperty("radStep" ,m_r_rstep ); declareProperty("maxSize" ,m_maxsize ); declareProperty("maxSizeSP" ,m_maxsizeSP ); declareProperty("minZ" ,m_zmin ); @@ -131,20 +131,20 @@ InDet::SiSpacePointsSeedMaker_HeavyIon::SiSpacePointsSeedMaker_HeavyIon InDet::SiSpacePointsSeedMaker_HeavyIon::~SiSpacePointsSeedMaker_HeavyIon() { - if(r_index ) delete [] r_index ; - if(r_map ) delete [] r_map ; - if(r_Sorted) delete [] r_Sorted; + if(m_r_index ) delete [] m_r_index ; + if(m_r_map ) delete [] m_r_map ; + if(m_r_Sorted) delete [] m_r_Sorted; // Delete seeds // - for(i_seed=l_seeds.begin(); i_seed!=l_seeds.end (); ++i_seed) { - delete *i_seed; + for(m_i_seed=m_l_seeds.begin(); m_i_seed!=m_l_seeds.end (); ++m_i_seed) { + delete *m_i_seed; } // Delete space points for reconstruction // - i_spforseed=l_spforseed.begin(); - for(; i_spforseed!=l_spforseed.end(); ++i_spforseed) { - delete *i_spforseed; + m_i_spforseed=m_l_spforseed.begin(); + for(; m_i_spforseed!=m_l_spforseed.end(); ++m_i_spforseed) { + delete *m_i_spforseed; } if(m_SP) delete [] m_SP; @@ -220,10 +220,10 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::newEvent (int) } else m_K = 2./(300.* 5. ); - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; - int irmax = r_size-1 ; + float irstep = 1./m_r_rstep; + int irmax = m_r_size-1 ; // Get pixels space points containers from store gate // @@ -244,12 +244,12 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r < 43. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r < 43. || r>=m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -274,12 +274,12 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -307,10 +307,10 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::newRegion } else m_K = 2./(300.* 5. ); - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; - int irmax = r_size-1 ; + float irstep = 1./m_r_rstep; + int irmax = m_r_size-1 ; // Get pixels space points containers from store gate // @@ -335,12 +335,12 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -370,12 +370,12 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -408,7 +408,7 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::find2Sp(const std::list<Trk::Vertex if(newv || !m_state || m_nspoint!=2 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); m_state = 1 ; m_nspoint = 2 ; m_nlist = 0 ; @@ -419,7 +419,7 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::find2Sp(const std::list<Trk::Vertex m_zMin = 0 ; production2Sp (); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -440,7 +440,7 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::find3Sp(const std::list<Trk::Vertex if(newv || !m_state || m_nspoint!=3 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 3 ; m_nlist = 0 ; @@ -451,7 +451,7 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::find3Sp(const std::list<Trk::Vertex m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -477,7 +477,7 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::findVSp (const std::list<Trk::Verte if(newv || !m_state || m_nspoint!=4 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 4 ; m_nlist = 0 ; @@ -488,7 +488,7 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::findVSp (const std::list<Trk::Verte m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -501,7 +501,8 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::findVSp (const std::list<Trk::Verte MsgStream& InDet::SiSpacePointsSeedMaker_HeavyIon::dump( MsgStream& out ) const { - if(m_nprint) return dumpEvent(out); return dumpConditions(out); + if(m_nprint) return dumpEvent(out); + return dumpConditions(out); } /////////////////////////////////////////////////////////////////// @@ -549,10 +550,10 @@ MsgStream& InDet::SiSpacePointsSeedMaker_HeavyIon::dumpConditions( MsgStream& ou <<std::setw(12)<<std::setprecision(5)<<m_rapcut <<" |"<<std::endl; out<<"| max radius SP | " - <<std::setw(12)<<std::setprecision(5)<<r_rmax + <<std::setw(12)<<std::setprecision(5)<<m_r_rmax <<" |"<<std::endl; out<<"| radius step | " - <<std::setw(12)<<std::setprecision(5)<<r_rstep + <<std::setw(12)<<std::setprecision(5)<<m_r_rstep <<" |"<<std::endl; out<<"| min Z-vertex position | " <<std::setw(12)<<std::setprecision(5)<<m_zmin @@ -661,7 +662,7 @@ MsgStream& InDet::SiSpacePointsSeedMaker_HeavyIon::dumpEvent( MsgStream& out ) c <<std::setw(12)<<m_nsazv <<" |"<<std::endl; out<<"| seeds | " - <<std::setw(12)<<l_seeds.size() + <<std::setw(12)<<m_l_seeds.size() <<" |"<<std::endl; out<<"|---------------------------------------------------------------------|" <<std::endl; @@ -736,12 +737,12 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::findNext () { if(m_endlist) return; - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); if (m_mode==0 || m_mode==1) production2Sp (); else if(m_mode==2 || m_mode==3) production3Sp (); else if(m_mode==5 || m_mode==6) production3Sp (); - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); ++m_nlist; } @@ -751,17 +752,17 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::findNext () bool InDet::SiSpacePointsSeedMaker_HeavyIon::newVertices(const std::list<Trk::Vertex>& lV) { - unsigned int s1 = l_vertex.size(); + unsigned int s1 = m_l_vertex.size(); unsigned int s2 = lV .size(); if(s1==0 && s2==0) return false; std::list<Trk::Vertex>::const_iterator v; - l_vertex.erase(l_vertex.begin(),l_vertex.end()); + m_l_vertex.erase(m_l_vertex.begin(),m_l_vertex.end()); for(v=lV.begin(); v!=lV.end(); ++v) { - l_vertex.push_back(float((*v).position().z())); - if(l_vertex.size() >= m_maxNumberVertices) break; + m_l_vertex.push_back(float((*v).position().z())); + if(m_l_vertex.size() >= m_maxNumberVertices) break; } return false; } @@ -776,7 +777,7 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::buildFrameWork() m_rapcut = fabs(m_rapcut) ; m_dzdrmax = 1./tan(2.*atan(exp(-m_rapcut))); m_dzdrmin =-m_dzdrmax ; - m_r3max = r_rmax ; + m_r3max = m_r_rmax ; m_COF = 134*.05*9. ; m_ipt = 1./fabs(.9*m_ptmin) ; m_ipt2 = m_ipt*m_ipt ; @@ -786,25 +787,25 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::buildFrameWork() // Build radius sorted containers // - r_size = int((r_rmax+.1)/r_rstep); - r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[r_size]; - r_index = new int[r_size]; - r_map = new int[r_size]; - m_nr = 0; for(int i=0; i!=r_size; ++i) {r_index[i]=0; r_map[i]=0;} + m_r_size = int((m_r_rmax+.1)/m_r_rstep); + m_r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[m_r_size]; + m_r_index = new int[m_r_size]; + m_r_map = new int[m_r_size]; + m_nr = 0; for(int i=0; i!=m_r_size; ++i) {m_r_index[i]=0; m_r_map[i]=0;} // Build radius-azimuthal sorted containers // const float pi2 = 2.*M_PI ; const int NFmax = 53 ; const float sFmax = float(NFmax )/pi2; - const float m_sFmin = 100./60. ; + const float sFmin = 100./60. ; - m_sF = m_ptmin /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < m_sFmin) m_sF = m_sFmin; + m_sF = m_ptmin /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < sFmin) m_sF = sFmin; m_fNmax = int(pi2*m_sF); if(m_fNmax >=NFmax) m_fNmax = NFmax-1; // Build radius-azimuthal-Z sorted containers // - m_nrfz = 0; for(int i=0; i!=583; ++i) {rfz_index [i]=0; rfz_map [i]=0;} + m_nrfz = 0; for(int i=0; i!=583; ++i) {m_rfz_index [i]=0; m_rfz_map [i]=0;} // Build radius-azimuthal-Z sorted containers for Z-vertices // @@ -812,7 +813,7 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::buildFrameWork() const float sFvmax = float(NFtmax)/pi2; m_sFv = m_ptmin/120. ; if(m_sFv >sFvmax) m_sFv = sFvmax; m_fvNmax = int(pi2*m_sFv); if(m_fvNmax>=NFtmax) m_fvNmax = NFtmax-1; - m_nrfzv = 0; for(int i=0; i!=300; ++i) {rfzv_index[i]=0; rfzv_map[i]=0;} + m_nrfzv = 0; for(int i=0; i!=300; ++i) {m_rfzv_index[i]=0; m_rfzv_map[i]=0;} // Build maps for radius-azimuthal-Z sorted collections // @@ -828,62 +829,62 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::buildFrameWork() int a = f *11+z; int b = fb*11+z; int c = ft*11+z; - rfz_b [a] = 3; rfz_t [a] = 3; - rfz_ib[a][0] = a; rfz_it[a][0] = a; - rfz_ib[a][1] = b; rfz_it[a][1] = b; - rfz_ib[a][2] = c; rfz_it[a][2] = c; + m_rfz_b [a] = 3; m_rfz_t [a] = 3; + m_rfz_ib[a][0] = a; m_rfz_it[a][0] = a; + m_rfz_ib[a][1] = b; m_rfz_it[a][1] = b; + m_rfz_ib[a][2] = c; m_rfz_it[a][2] = c; if (z==5) { - rfz_t [a] = 9 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; - rfz_it[a][6] = a-1; - rfz_it[a][7] = b-1; - rfz_it[a][8] = c-1; + m_rfz_t [a] = 9 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; + m_rfz_it[a][6] = a-1; + m_rfz_it[a][7] = b-1; + m_rfz_it[a][8] = c-1; } else if(z> 5) { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a-1; - rfz_ib[a][4] = b-1; - rfz_ib[a][5] = c-1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a-1; + m_rfz_ib[a][4] = b-1; + m_rfz_ib[a][5] = c-1; if(z<10) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; } } else { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a+1; - rfz_ib[a][4] = b+1; - rfz_ib[a][5] = c+1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a+1; + m_rfz_ib[a][4] = b+1; + m_rfz_ib[a][5] = c+1; if(z>0) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a-1; - rfz_it[a][4] = b-1; - rfz_it[a][5] = c-1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a-1; + m_rfz_it[a][4] = b-1; + m_rfz_it[a][5] = c-1; } } if (z==3) { - rfz_b[a] = 9; - rfz_ib[a][6] = a+2; - rfz_ib[a][7] = b+2; - rfz_ib[a][8] = c+2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a+2; + m_rfz_ib[a][7] = b+2; + m_rfz_ib[a][8] = c+2; } else if(z==7) { - rfz_b[a] = 9; - rfz_ib[a][6] = a-2; - rfz_ib[a][7] = b-2; - rfz_ib[a][8] = c-2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a-2; + m_rfz_ib[a][7] = b-2; + m_rfz_ib[a][8] = c-2; } } } @@ -902,21 +903,21 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::buildFrameWork() int a = f *3+z; int b = fb*3+z; int c = ft*3+z; - rfzv_n[a] = 3; - rfzv_i[a][0] = a; - rfzv_i[a][1] = b; - rfzv_i[a][2] = c; + m_rfzv_n[a] = 3; + m_rfzv_i[a][0] = a; + m_rfzv_i[a][1] = b; + m_rfzv_i[a][2] = c; if (z>1) { - rfzv_n[a] = 6; - rfzv_i[a][3] = a-1; - rfzv_i[a][4] = b-1; - rfzv_i[a][5] = c-1; + m_rfzv_n[a] = 6; + m_rfzv_i[a][3] = a-1; + m_rfzv_i[a][4] = b-1; + m_rfzv_i[a][5] = c-1; } else if(z<1) { - rfzv_n[a] = 6; - rfzv_i[a][3] = a+1; - rfzv_i[a][4] = b+1; - rfzv_i[a][5] = c+1; + m_rfzv_n[a] = 6; + m_rfzv_i[a][3] = a+1; + m_rfzv_i[a][4] = b+1; + m_rfzv_i[a][5] = c+1; } } } @@ -930,8 +931,8 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::buildFrameWork() if(!m_Zo) m_Zo = new float[m_maxsizeSP]; if(!m_OneSeeds) m_OneSeeds = new InDet::SiSpacePointsSeed [m_maxOneSize]; - i_seed = l_seeds.begin(); - i_seede = l_seeds.end (); + m_i_seed = m_l_seeds.begin(); + m_i_seede = m_l_seeds.end (); } /////////////////////////////////////////////////////////////////// @@ -990,11 +991,11 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::fillLists() const float pi2 = 2.*M_PI; std::list<InDet::SiSpacePointForSeed*>::iterator r; - for(int i=0; i!= r_size; ++i) { + for(int i=0; i!= m_r_size; ++i) { - if(!r_map[i]) continue; r = r_Sorted[i].begin(); + if(!m_r_map[i]) continue; r = m_r_Sorted[i].begin(); - while(r!=r_Sorted[i].end()) { + while(r!=m_r_Sorted[i].end()) { // Azimuthal angle sort // @@ -1013,7 +1014,7 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::fillLists() Z>-250.?z=5:Z>-450.?z=4:Z>-925.?z=3:Z>-1400.?z=2:Z>-2500.?z=1:z= 0; } int n = f*11+z; ++m_nsaz; - rfz_Sorted[n].push_back(*r); if(!rfz_map[n]++) rfz_index[m_nrfz++] = n; + m_rfz_Sorted[n].push_back(*r); if(!m_rfz_map[n]++) m_rfz_index[m_nrfz++] = n; if((*r)->spacepoint->clusterList().second == 0 && z>=3 && z<=7) { z<=4 ? z=0 : z>=6 ? z=2 : z=1; @@ -1022,11 +1023,11 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::fillLists() // f = int(F*m_sFv); f<0 ? f+=m_fvNmax : f> m_fvNmax ? f-=m_fvNmax : f=f; n = f*3+z; ++m_nsazv; - rfzv_Sorted[n].push_back(*r); if(!rfzv_map[n]++) rfzv_index[m_nrfzv++] = n; + m_rfzv_Sorted[n].push_back(*r); if(!m_rfzv_map[n]++) m_rfzv_index[m_nrfzv++] = n; } - r_Sorted[i].erase(r++); + m_r_Sorted[i].erase(r++); } - r_map[i] = 0; + m_r_map[i] = 0; } m_nr = 0; m_state = 0; @@ -1039,18 +1040,18 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::fillLists() void InDet::SiSpacePointsSeedMaker_HeavyIon::erase() { for(int i=0; i!=m_nr; ++i) { - int n = r_index[i]; r_map[n] = 0; - r_Sorted[n].clear(); + int n = m_r_index[i]; m_r_map[n] = 0; + m_r_Sorted[n].clear(); } for(int i=0; i!=m_nrfz; ++i) { - int n = rfz_index[i]; rfz_map[n] = 0; - rfz_Sorted[n].clear(); + int n = m_rfz_index[i]; m_rfz_map[n] = 0; + m_rfz_Sorted[n].clear(); } for(int i=0; i!=m_nrfzv; ++i) { - int n = rfzv_index[i]; rfzv_map[n] = 0; - rfzv_Sorted[n].clear(); + int n = m_rfzv_index[i]; m_rfzv_map[n] = 0; + m_rfzv_Sorted[n].clear(); } m_state = 0; m_ns = 0; @@ -1081,9 +1082,9 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::production2Sp() int z = 0; if(!m_endlist) z = m_zMin; for(; z!=3; ++z) { - int a = f*3+z; if(!rfzv_map[a]) continue; - r0 = rfzv_Sorted[a].begin(); - r0e = rfzv_Sorted[a].end (); + int a = f*3+z; if(!m_rfzv_map[a]) continue; + r0 = m_rfzv_Sorted[a].begin(); + r0e = m_rfzv_Sorted[a].end (); if(!m_endlist) {r0 = m_rMin; m_endlist = true;} @@ -1101,14 +1102,14 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::production2Sp() // Bottom links production // - int NB = rfzv_n[a]; + int NB = m_rfzv_n[a]; for(int i=0; i!=NB; ++i) { - int an = rfzv_i[a][i]; - if(!rfzv_map[an]) continue; + int an = m_rfzv_i[a][i]; + if(!m_rfzv_map[an]) continue; - r = rfzv_Sorted[an].begin(); - re = rfzv_Sorted[an].end (); + r = m_rfzv_Sorted[an].begin(); + re = m_rfzv_Sorted[an].end (); for(; r!=re; ++r) { @@ -1172,19 +1173,19 @@ void InDet::SiSpacePointsSeedMaker_HeavyIon::production3Sp() for(; z!=11; ++z) { - int a = f *11+ZI[z]; if(!rfz_map[a]) continue; + int a = f *11+ZI[z]; if(!m_rfz_map[a]) continue; int NB = 0, NT = 0; - for(int i=0; i!=rfz_b[a]; ++i) { + for(int i=0; i!=m_rfz_b[a]; ++i) { - int an = rfz_ib[a][i]; - if(!rfz_map[an]) continue; - rb [NB] = rfz_Sorted[an].begin(); rbe[NB++] = rfz_Sorted[an].end(); + int an = m_rfz_ib[a][i]; + if(!m_rfz_map[an]) continue; + rb [NB] = m_rfz_Sorted[an].begin(); rbe[NB++] = m_rfz_Sorted[an].end(); } - for(int i=0; i!=rfz_t[a]; ++i) { + for(int i=0; i!=m_rfz_t[a]; ++i) { - int an = rfz_it[a][i]; - if(!rfz_map[an]) continue; - rt [NT] = rfz_Sorted[an].begin(); rte[NT++] = rfz_Sorted[an].end(); + int an = m_rfz_it[a][i]; + if(!m_rfz_map[an]) continue; + rt [NT] = m_rfz_Sorted[an].begin(); rte[NT++] = m_rfz_Sorted[an].end(); } if(m_izvertex) { if(!m_trigger) production3Sp (rb,rbe,rt,rte,NB,NT,nseed); diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_ITK.cxx b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_ITK.cxx index 767120bf27441a48eb8bd1d0d804aa686ca5a058..d13678b3e362797b0b3388cbaec41f4591f7859a 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_ITK.cxx +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_ITK.cxx @@ -67,14 +67,14 @@ InDet::SiSpacePointsSeedMaker_ITK::SiSpacePointsSeedMaker_ITK m_diversss = 20 ; m_divermax = 20. ; m_dazmax = .02 ; - r_rmax = 1100. ; - r_rmin = 0. ; + m_r_rmax = 1100. ; + m_r_rmin = 0. ; m_umax = 0. ; - r_rstep = 2. ; + m_r_rstep = 2. ; m_dzmaxPPP = 600. ; - r_Sorted = 0 ; - r_index = 0 ; - r_map = 0 ; + m_r_Sorted = 0 ; + m_r_index = 0 ; + m_r_map = 0 ; m_maxsizeSP = 5000 ; m_maxOneSize= 5 ; m_SP = 0 ; @@ -112,9 +112,9 @@ InDet::SiSpacePointsSeedMaker_ITK::SiSpacePointsSeedMaker_ITK declareProperty("etaMin" ,m_etamin ); declareProperty("etaMax" ,m_etamax ); declareProperty("pTmin" ,m_ptmin ); - declareProperty("radMax" ,r_rmax ); - declareProperty("radMin" ,r_rmin ); - declareProperty("radStep" ,r_rstep ); + declareProperty("radMax" ,m_r_rmax ); + declareProperty("radMin" ,m_r_rmin ); + declareProperty("radStep" ,m_r_rstep ); declareProperty("maxSize" ,m_maxsize ); declareProperty("maxSizeSP" ,m_maxsizeSP ); declareProperty("minZ" ,m_zmin ); @@ -157,20 +157,20 @@ InDet::SiSpacePointsSeedMaker_ITK::SiSpacePointsSeedMaker_ITK InDet::SiSpacePointsSeedMaker_ITK::~SiSpacePointsSeedMaker_ITK() { - if(r_index ) delete [] r_index ; - if(r_map ) delete [] r_map ; - if(r_Sorted) delete [] r_Sorted; + if(m_r_index ) delete [] m_r_index ; + if(m_r_map ) delete [] m_r_map ; + if(m_r_Sorted) delete [] m_r_Sorted; // Delete seeds // - for(i_seed=l_seeds.begin(); i_seed!=l_seeds.end (); ++i_seed) { - delete *i_seed; + for(m_i_seed=m_l_seeds.begin(); m_i_seed!=m_l_seeds.end (); ++m_i_seed) { + delete *m_i_seed; } // Delete space points for reconstruction // - i_spforseed=l_spforseed.begin(); - for(; i_spforseed!=l_spforseed.end(); ++i_spforseed) { - delete *i_spforseed; + m_i_spforseed=m_l_spforseed.begin(); + for(; m_i_spforseed!=m_l_spforseed.end(); ++m_i_spforseed) { + delete *m_i_spforseed; } if(m_seedOutput) delete m_seedOutput; @@ -221,7 +221,7 @@ StatusCode InDet::SiSpacePointsSeedMaker_ITK::initialize() } } - if(r_rmax < 1100.) r_rmax = 1100.; + if(m_r_rmax < 1100.) m_r_rmax = 1100.; // Build framework // @@ -274,22 +274,22 @@ void InDet::SiSpacePointsSeedMaker_ITK::newEvent(int iteration) m_ipt2K = m_ipt2/(m_K*m_K); m_ipt2C = m_ipt2*m_COF ; m_COFK = m_COF*(m_K*m_K) ; - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); } else { - r_first = 0; fillLists(); return; + m_r_first = 0; fillLists(); return; } m_checketa = m_dzdrmin > 1.; - float irstep = 1./r_rstep; - int irmax = r_size-1 ; - for(int i=0; i!=m_nr; ++i) {int n = r_index[i]; r_map[n] = 0; r_Sorted[n].clear();} + float irstep = 1./m_r_rstep; + int irmax = m_r_size-1 ; + for(int i=0; i!=m_nr; ++i) {int n = m_r_index[i]; m_r_map[n] = 0; m_r_Sorted[n].clear();} m_ns = m_nr = 0; // Get pixels space points containers from store gate // - r_first = 0; + m_r_first = 0; if(m_pixel) { // m_spacepointsPixel = 0; @@ -307,19 +307,19 @@ void InDet::SiSpacePointsSeedMaker_ITK::newEvent(int iteration) for(; sp != spe; ++sp) { - if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > r_rmax || (*sp)->r() < r_rmin ) continue; + if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > m_r_rmax || (*sp)->r() < m_r_rmin ) continue; InDet::SiSpacePointForSeedITK* sps = newSpacePoint((*sp)); if(!sps) continue; int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; - if(ir > r_first) r_first = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; + if(ir > m_r_first) m_r_first = ir; ++m_ns; } } } - ++r_first; + ++m_r_first; } // Get sct space points containers from store gate @@ -340,13 +340,13 @@ void InDet::SiSpacePointsSeedMaker_ITK::newEvent(int iteration) for(; sp != spe; ++sp) { - if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > r_rmax || (*sp)->r() < r_rmin ) continue; + if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > m_r_rmax || (*sp)->r() < m_r_rmin ) continue; InDet::SiSpacePointForSeedITK* sps = newSpacePoint((*sp)); if(!sps) continue; int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -365,20 +365,20 @@ void InDet::SiSpacePointsSeedMaker_ITK::newEvent(int iteration) for (; sp!=spe; ++sp) { - if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > r_rmax || (*sp)->r() < r_rmin) continue; + if ((m_useassoTool && isUsed(*sp)) || (*sp)->r() > m_r_rmax || (*sp)->r() < m_r_rmin) continue; InDet::SiSpacePointForSeedITK* sps = newSpacePoint((*sp)); if(!sps) continue; int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } } } - if(iteration < 0) r_first = 0; + if(iteration < 0) m_r_first = 0; fillLists(); } @@ -410,15 +410,15 @@ void InDet::SiSpacePointsSeedMaker_ITK::newRegion m_ipt2C = m_ipt2*m_COF ; m_COFK = m_COF*(m_K*m_K) ; - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; - int irmax = r_size-1 ; + float irstep = 1./m_r_rstep; + int irmax = m_r_size-1 ; - r_first = 0 ; + m_r_first = 0 ; m_checketa = false ; - for(int i=0; i!=m_nr; ++i) {int n = r_index[i]; r_map[n] = 0; r_Sorted[n].clear();} + for(int i=0; i!=m_nr; ++i) {int n = m_r_index[i]; m_r_map[n] = 0; m_r_Sorted[n].clear();} m_ns = m_nr = 0; // Get pixels space points containers from store gate @@ -444,11 +444,11 @@ void InDet::SiSpacePointsSeedMaker_ITK::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r > r_rmax || r < r_rmin) continue; + float r = (*sp)->r(); if(r > m_r_rmax || r < m_r_rmin) continue; InDet::SiSpacePointForSeedITK* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -478,11 +478,11 @@ void InDet::SiSpacePointsSeedMaker_ITK::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r > r_rmax || r < r_rmin) continue; + float r = (*sp)->r(); if(r > m_r_rmax || r < m_r_rmin) continue; InDet::SiSpacePointForSeedITK* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -530,7 +530,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::find2Sp(const std::list<Trk::Vertex>& lv if(newv || !m_state || m_nspoint!=2 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); m_state = 1 ; m_nspoint = 2 ; m_nlist = 0 ; @@ -541,7 +541,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::find2Sp(const std::list<Trk::Vertex>& lv m_zMin = 0 ; production2Sp (); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -562,7 +562,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::find3Sp(const std::list<Trk::Vertex>& lv bool newv = newVertices(lv); if(newv || !m_state || m_nspoint!=3 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 3 ; m_nlist = 0 ; @@ -573,7 +573,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::find3Sp(const std::list<Trk::Vertex>& lv m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_seeds.begin(); if(m_outputlevel<=0) { @@ -595,7 +595,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::find3Sp(const std::list<Trk::Vertex>& lv bool newv = newVertices(lv); if(newv || !m_state || m_nspoint!=3 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 3 ; m_nlist = 0 ; @@ -606,7 +606,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::find3Sp(const std::list<Trk::Vertex>& lv m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_seeds.begin(); if(m_outputlevel<=0) { @@ -630,7 +630,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::findVSp (const std::list<Trk::Vertex>& l if(newv || !m_state || m_nspoint!=4 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 4 ; m_nlist = 0 ; @@ -641,7 +641,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::findVSp (const std::list<Trk::Vertex>& l m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_seeds.begin(); if(m_outputlevel<=0) { @@ -703,13 +703,13 @@ MsgStream& InDet::SiSpacePointsSeedMaker_ITK::dumpConditions( MsgStream& out ) c <<std::setw(12)<<std::setprecision(5)<<m_rapcut <<" |"<<std::endl; out<<"| max radius SP | " - <<std::setw(12)<<std::setprecision(5)<<r_rmax + <<std::setw(12)<<std::setprecision(5)<<m_r_rmax <<" |"<<std::endl; out<<"| min radius SP | " - <<std::setw(12)<<std::setprecision(5)<<r_rmin + <<std::setw(12)<<std::setprecision(5)<<m_r_rmin <<" |"<<std::endl; out<<"| radius step | " - <<std::setw(12)<<std::setprecision(5)<<r_rstep + <<std::setw(12)<<std::setprecision(5)<<m_r_rstep <<" |"<<std::endl; out<<"| min Z-vertex position | " <<std::setw(12)<<std::setprecision(5)<<m_zmin @@ -818,7 +818,7 @@ MsgStream& InDet::SiSpacePointsSeedMaker_ITK::dumpEvent( MsgStream& out ) const <<std::setw(12)<<m_nsazv <<" |"<<std::endl; out<<"| seeds | " - <<std::setw(12)<<l_seeds.size() + <<std::setw(12)<<m_l_seeds.size() <<" |"<<std::endl; out<<"|---------------------------------------------------------------------|" <<std::endl; @@ -893,13 +893,13 @@ void InDet::SiSpacePointsSeedMaker_ITK::findNext () { if(m_endlist) return; - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); if (m_mode==0 || m_mode==1) production2Sp (); else if(m_mode==2 || m_mode==3) production3Sp (); else if(m_mode==5 || m_mode==6) production3Sp (); - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_seeds.begin(); ++m_nlist; } @@ -910,21 +910,21 @@ void InDet::SiSpacePointsSeedMaker_ITK::findNext () bool InDet::SiSpacePointsSeedMaker_ITK::newVertices(const std::list<Trk::Vertex>& lV) { - unsigned int s1 = l_vertex.size(); + unsigned int s1 = m_l_vertex.size(); unsigned int s2 = lV .size(); m_isvertex = false; if(s1==0 && s2==0) return false; std::list<Trk::Vertex>::const_iterator v; - l_vertex.clear(); + m_l_vertex.clear(); if(s2 == 0) return false; m_isvertex = true; - for(v=lV.begin(); v!=lV.end(); ++v) {l_vertex.insert(float((*v).position().z()));} + for(v=lV.begin(); v!=lV.end(); ++v) {m_l_vertex.insert(float((*v).position().z()));} - m_zminU = (*l_vertex. begin())-20.; if( m_zminU < m_zmin) m_zminU = m_zmin; - m_zmaxU = (*l_vertex.rbegin())+20.; if( m_zmaxU > m_zmax) m_zmaxU = m_zmax; + m_zminU = (*m_l_vertex. begin())-20.; if( m_zminU < m_zmin) m_zminU = m_zmin; + m_zmaxU = (*m_l_vertex.rbegin())+20.; if( m_zmaxU > m_zmax) m_zmaxU = m_zmax; return false; } @@ -955,28 +955,28 @@ void InDet::SiSpacePointsSeedMaker_ITK::buildFrameWork() // Build radius sorted containers // - r_size = int((r_rmax+.1)/r_rstep); - r_Sorted = new std::list<InDet::SiSpacePointForSeedITK*>[r_size]; - r_index = new int[r_size]; - r_map = new int[r_size]; - m_nr = 0; for(int i=0; i!=r_size; ++i) {r_index[i]=0; r_map[i]=0;} + m_r_size = int((m_r_rmax+.1)/m_r_rstep); + m_r_Sorted = new std::list<InDet::SiSpacePointForSeedITK*>[m_r_size]; + m_r_index = new int[m_r_size]; + m_r_map = new int[m_r_size]; + m_nr = 0; for(int i=0; i!=m_r_size; ++i) {m_r_index[i]=0; m_r_map[i]=0;} // Build radius-azimuthal sorted containers // const float pi2 = 2.*M_PI ; const int NFmax = 53 ; const float sFmax = float(NFmax )/pi2; - const float m_sFmin = 100./60. ; + const float sFmin = 100./60. ; float ptm = 400.; if(m_ptmin < ptm) ptm = m_ptmin; - m_sF = ptm /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < m_sFmin) m_sF = m_sFmin; + m_sF = ptm /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < sFmin) m_sF = sFmin; m_fNmax = int(pi2*m_sF); if(m_fNmax >=NFmax) m_fNmax = NFmax-1; // Build radius-azimuthal-Z sorted containers // - m_nrfz = 0; for(int i=0; i!=583; ++i) {rfz_index [i]=0; rfz_map [i]=0;} + m_nrfz = 0; for(int i=0; i!=583; ++i) {m_rfz_index [i]=0; m_rfz_map [i]=0;} // Build radius-azimuthal-Z sorted containers for Z-vertices // @@ -984,7 +984,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::buildFrameWork() const float sFvmax = float(NFtmax)/pi2; m_sFv = m_ptmin/120. ; if(m_sFv >sFvmax) m_sFv = sFvmax; m_fvNmax = int(pi2*m_sFv); if(m_fvNmax>=NFtmax) m_fvNmax = NFtmax-1; - m_nrfzv = 0; for(int i=0; i!=300; ++i) {rfzv_index[i]=0; rfzv_map[i]=0;} + m_nrfzv = 0; for(int i=0; i!=300; ++i) {m_rfzv_index[i]=0; m_rfzv_map[i]=0;} // Build maps for radius-azimuthal-Z sorted collections // @@ -1000,62 +1000,62 @@ void InDet::SiSpacePointsSeedMaker_ITK::buildFrameWork() int a = f *11+z; int b = fb*11+z; int c = ft*11+z; - rfz_b [a] = 3; rfz_t [a] = 3; - rfz_ib[a][0] = a; rfz_it[a][0] = a; - rfz_ib[a][1] = b; rfz_it[a][1] = b; - rfz_ib[a][2] = c; rfz_it[a][2] = c; + m_rfz_b [a] = 3; m_rfz_t [a] = 3; + m_rfz_ib[a][0] = a; m_rfz_it[a][0] = a; + m_rfz_ib[a][1] = b; m_rfz_it[a][1] = b; + m_rfz_ib[a][2] = c; m_rfz_it[a][2] = c; if (z==5) { - rfz_t [a] = 9 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; - rfz_it[a][6] = a-1; - rfz_it[a][7] = b-1; - rfz_it[a][8] = c-1; + m_rfz_t [a] = 9 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; + m_rfz_it[a][6] = a-1; + m_rfz_it[a][7] = b-1; + m_rfz_it[a][8] = c-1; } else if(z> 5) { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a-1; - rfz_ib[a][4] = b-1; - rfz_ib[a][5] = c-1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a-1; + m_rfz_ib[a][4] = b-1; + m_rfz_ib[a][5] = c-1; if(z<10) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; } } else { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a+1; - rfz_ib[a][4] = b+1; - rfz_ib[a][5] = c+1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a+1; + m_rfz_ib[a][4] = b+1; + m_rfz_ib[a][5] = c+1; if(z>0) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a-1; - rfz_it[a][4] = b-1; - rfz_it[a][5] = c-1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a-1; + m_rfz_it[a][4] = b-1; + m_rfz_it[a][5] = c-1; } } if (z==3) { - rfz_b[a] = 9; - rfz_ib[a][6] = a+2; - rfz_ib[a][7] = b+2; - rfz_ib[a][8] = c+2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a+2; + m_rfz_ib[a][7] = b+2; + m_rfz_ib[a][8] = c+2; } else if(z==7) { - rfz_b[a] = 9; - rfz_ib[a][6] = a-2; - rfz_ib[a][7] = b-2; - rfz_ib[a][8] = c-2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a-2; + m_rfz_ib[a][7] = b-2; + m_rfz_ib[a][8] = c-2; } } } @@ -1074,21 +1074,21 @@ void InDet::SiSpacePointsSeedMaker_ITK::buildFrameWork() int a = f *3+z; int b = fb*3+z; int c = ft*3+z; - rfzv_n[a] = 3; - rfzv_i[a][0] = a; - rfzv_i[a][1] = b; - rfzv_i[a][2] = c; + m_rfzv_n[a] = 3; + m_rfzv_i[a][0] = a; + m_rfzv_i[a][1] = b; + m_rfzv_i[a][2] = c; if (z>1) { - rfzv_n[a] = 6; - rfzv_i[a][3] = a-1; - rfzv_i[a][4] = b-1; - rfzv_i[a][5] = c-1; + m_rfzv_n[a] = 6; + m_rfzv_i[a][3] = a-1; + m_rfzv_i[a][4] = b-1; + m_rfzv_i[a][5] = c-1; } else if(z<1) { - rfzv_n[a] = 6; - rfzv_i[a][3] = a+1; - rfzv_i[a][4] = b+1; - rfzv_i[a][5] = c+1; + m_rfzv_n[a] = 6; + m_rfzv_i[a][3] = a+1; + m_rfzv_i[a][4] = b+1; + m_rfzv_i[a][5] = c+1; } } } @@ -1106,8 +1106,8 @@ void InDet::SiSpacePointsSeedMaker_ITK::buildFrameWork() if(!m_seedOutput) m_seedOutput = new InDet::SiSpacePointsSeed(); - i_seed = l_seeds.begin(); - i_seede = l_seeds.end (); + m_i_seed = m_l_seeds.begin(); + m_i_seede = m_l_seeds.end (); } /////////////////////////////////////////////////////////////////// @@ -1207,9 +1207,9 @@ void InDet::SiSpacePointsSeedMaker_ITK::fillLists() int ir0 =0; - for(int i=r_first; i!=r_size; ++i) { + for(int i=m_r_first; i!=m_r_size; ++i) { - if(!r_map[i]) continue; r = r_Sorted[i].begin(); re = r_Sorted[i].end(); + if(!m_r_map[i]) continue; r = m_r_Sorted[i].begin(); re = m_r_Sorted[i].end(); if(!ir0) ir0 = i; if(m_iteration && (*r)->spacepoint->clusterList().second) break; @@ -1234,7 +1234,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::fillLists() } int n = f*11+z; ++m_nsaz; - rfz_Sorted[n].push_back(*r); if(!rfz_map[n]++) rfz_index[m_nrfz++] = n; + m_rfz_Sorted[n].push_back(*r); if(!m_rfz_map[n]++) m_rfz_index[m_nrfz++] = n; if(!m_iteration && (*r)->spacepoint->clusterList().second == 0 && z>=3 && z<=7) { z<=4 ? z=0 : z>=6 ? z=2 : z=1; @@ -1243,7 +1243,7 @@ void InDet::SiSpacePointsSeedMaker_ITK::fillLists() // f = int(F*m_sFv); f<0 ? f+=m_fvNmax : f> m_fvNmax ? f-=m_fvNmax : f=f; n = f*3+z; ++m_nsazv; - rfzv_Sorted[n].push_back(*r); if(!rfzv_map[n]++) rfzv_index[m_nrfzv++] = n; + m_rfzv_Sorted[n].push_back(*r); if(!m_rfzv_map[n]++) m_rfzv_index[m_nrfzv++] = n; } } } @@ -1257,13 +1257,13 @@ void InDet::SiSpacePointsSeedMaker_ITK::fillLists() void InDet::SiSpacePointsSeedMaker_ITK::erase() { for(int i=0; i!=m_nrfz; ++i) { - int n = rfz_index[i]; rfz_map[n] = 0; - rfz_Sorted[n].clear(); + int n = m_rfz_index[i]; m_rfz_map[n] = 0; + m_rfz_Sorted[n].clear(); } for(int i=0; i!=m_nrfzv; ++i) { - int n = rfzv_index[i]; rfzv_map[n] = 0; - rfzv_Sorted[n].clear(); + int n = m_rfzv_index[i]; m_rfzv_map[n] = 0; + m_rfzv_Sorted[n].clear(); } m_state = 0; m_nsaz = 0; @@ -1307,9 +1307,9 @@ void InDet::SiSpacePointsSeedMaker_ITK::production2Sp() int z = 0; if(!m_endlist) z = m_zMin; for(; z!=3; ++z) { - int a = f*3+z; if(!rfzv_map[a]) continue; - r0 = rfzv_Sorted[a].begin(); - r0e = rfzv_Sorted[a].end (); + int a = f*3+z; if(!m_rfzv_map[a]) continue; + r0 = m_rfzv_Sorted[a].begin(); + r0e = m_rfzv_Sorted[a].end (); if(!m_endlist) {r0 = m_rMin; m_endlist = true;} @@ -1327,14 +1327,14 @@ void InDet::SiSpacePointsSeedMaker_ITK::production2Sp() // Bottom links production // - int NB = rfzv_n[a]; + int NB = m_rfzv_n[a]; for(int i=0; i!=NB; ++i) { - int an = rfzv_i[a][i]; - if(!rfzv_map[an]) continue; + int an = m_rfzv_i[a][i]; + if(!m_rfzv_map[an]) continue; - r = rfzv_Sorted[an].begin(); - re = rfzv_Sorted[an].end (); + r = m_rfzv_Sorted[an].begin(); + re = m_rfzv_Sorted[an].end (); for(; r!=re; ++r) { @@ -1399,19 +1399,19 @@ void InDet::SiSpacePointsSeedMaker_ITK::production3Sp() for(; z!=11; ++z) { - int a = f *11+ZI[z]; if(!rfz_map[a]) continue; + int a = f *11+ZI[z]; if(!m_rfz_map[a]) continue; int NB = 0, NT = 0; - for(int i=0; i!=rfz_b[a]; ++i) { + for(int i=0; i!=m_rfz_b[a]; ++i) { - int an = rfz_ib[a][i]; - if(!rfz_map[an]) continue; - rb [NB] = rfz_Sorted[an].begin(); rbe[NB++] = rfz_Sorted[an].end(); + int an = m_rfz_ib[a][i]; + if(!m_rfz_map[an]) continue; + rb [NB] = m_rfz_Sorted[an].begin(); rbe[NB++] = m_rfz_Sorted[an].end(); } - for(int i=0; i!=rfz_t[a]; ++i) { + for(int i=0; i!=m_rfz_t[a]; ++i) { - int an = rfz_it[a][i]; - if(!rfz_map[an]) continue; - rt [NT] = rfz_Sorted[an].begin(); rte[NT++] = rfz_Sorted[an].end(); + int an = m_rfz_it[a][i]; + if(!m_rfz_map[an]) continue; + rt [NT] = m_rfz_Sorted[an].begin(); rte[NT++] = m_rfz_Sorted[an].end(); } if(m_iteration == 0 && m_iteration0 ==0) production3SpSSS(rb,rbe,rt,rte,NB,NT,nseed); @@ -2080,14 +2080,14 @@ void InDet::SiSpacePointsSeedMaker_ITK::fillSeeds () if(l!=lf && s->spacepoint0()->radius() < 43. && w > -200.) continue; if(!s->setQuality(w)) continue; - if(i_seede!=l_seeds.end()) { - s = (*i_seede++); + if(m_i_seede!=m_l_seeds.end()) { + s = (*m_i_seede++); *s = *(*l).second; } else { s = new SiSpacePointsProSeedITK(*(*l).second); - l_seeds.push_back(s); - i_seede = l_seeds.end(); + m_l_seeds.push_back(s); + m_i_seede = m_l_seeds.end(); } if (s->spacepoint0()->spacepoint->clusterList().second) w-=3000.; diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_LowMomentum.cxx b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_LowMomentum.cxx index fe858ae13e6a2c4d34c8f13e78c984fa27519ab0..2e2c7cc39b07afa65982f994df2fae2f5c39bbcb 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_LowMomentum.cxx +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_LowMomentum.cxx @@ -60,11 +60,11 @@ InDet::SiSpacePointsSeedMaker_LowMomentum::SiSpacePointsSeedMaker_LowMomentum m_diverpps = 1.2 ; m_dazmax = .02 ; m_iptmax = 1./400. ; - r_rmax = 200. ; - r_rstep = 2. ; - r_Sorted = 0 ; - r_index = 0 ; - r_map = 0 ; + m_r_rmax = 200. ; + m_r_rstep = 2. ; + m_r_Sorted = 0 ; + m_r_index = 0 ; + m_r_map = 0 ; m_maxsizeSP = 1500 ; m_maxOneSize = 5 ; m_SP = 0 ; @@ -95,8 +95,8 @@ InDet::SiSpacePointsSeedMaker_LowMomentum::SiSpacePointsSeedMaker_LowMomentum declareProperty("useSCT" ,m_sct ); declareProperty("pTmin" ,m_ptmin ); declareProperty("pTmax" ,m_ptmax ); - declareProperty("radMax" ,r_rmax ); - declareProperty("radStep" ,r_rstep ); + declareProperty("radMax" ,m_r_rmax ); + declareProperty("radStep" ,m_r_rstep ); declareProperty("maxSize" ,m_maxsize ); declareProperty("maxSizeSP" ,m_maxsizeSP ); declareProperty("minZ" ,m_zmin ); @@ -130,20 +130,20 @@ InDet::SiSpacePointsSeedMaker_LowMomentum::SiSpacePointsSeedMaker_LowMomentum InDet::SiSpacePointsSeedMaker_LowMomentum::~SiSpacePointsSeedMaker_LowMomentum() { - if(r_index ) delete [] r_index ; - if(r_map ) delete [] r_map ; - if(r_Sorted) delete [] r_Sorted; + if(m_r_index ) delete [] m_r_index ; + if(m_r_map ) delete [] m_r_map ; + if(m_r_Sorted) delete [] m_r_Sorted; // Delete seeds // - for(i_seed=l_seeds.begin(); i_seed!=l_seeds.end (); ++i_seed) { - delete *i_seed; + for(m_i_seed=m_l_seeds.begin(); m_i_seed!=m_l_seeds.end (); ++m_i_seed) { + delete *m_i_seed; } // Delete space points for reconstruction // - i_spforseed=l_spforseed.begin(); - for(; i_spforseed!=l_spforseed.end(); ++i_spforseed) { - delete *i_spforseed; + m_i_spforseed=m_l_spforseed.begin(); + for(; m_i_spforseed!=m_l_spforseed.end(); ++m_i_spforseed) { + delete *m_i_spforseed; } if(m_SP) delete [] m_SP; @@ -221,11 +221,11 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::newEvent (int) { m_trigger = false; if(!m_pixel && !m_sct) return; erase(); - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); buildBeamFrameWork(); - float irstep = 1./r_rstep; - int irmax = r_size-1 ; + float irstep = 1./m_r_rstep; + int irmax = m_r_size-1 ; // Get pixels space points containers from store gate // @@ -246,14 +246,14 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -278,14 +278,14 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -303,11 +303,11 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::newRegion { m_trigger = false; if(!m_pixel && !m_sct) return; erase(); - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); buildBeamFrameWork(); - int irmax = r_size-1 ; - float irstep = 1./r_rstep; + int irmax = m_r_size-1 ; + float irstep = 1./m_r_rstep; // Get pixels space points containers from store gate // @@ -332,14 +332,14 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -369,14 +369,14 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; if(m_useassoTool && isUsed(*sp) ) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -407,7 +407,7 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::find2Sp(const std::list<Trk::Ver if(newv || !m_state || m_nspoint!=2 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); m_state = 1 ; m_nspoint = 2 ; m_nlist = 0 ; @@ -417,7 +417,7 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::find2Sp(const std::list<Trk::Ver m_zMin = 0 ; production2Sp (); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -436,7 +436,7 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::find3Sp(const std::list<Trk::Ver if(newv || !m_state || m_nspoint!=3 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 3 ; m_nlist = 0 ; @@ -446,7 +446,7 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::find3Sp(const std::list<Trk::Ver m_zMin = 0 ; production3Sp (); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -471,7 +471,7 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::findVSp (const std::list<Trk::Ve if(newv || !m_state || m_nspoint!=4 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 4 ; m_nlist = 0 ; @@ -481,7 +481,7 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::findVSp (const std::list<Trk::Ve m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -547,10 +547,10 @@ MsgStream& InDet::SiSpacePointsSeedMaker_LowMomentum::dumpConditions( MsgStream& <<std::setw(12)<<std::setprecision(5)<<m_rapcut <<" |"<<std::endl; out<<"| max radius SP | " - <<std::setw(12)<<std::setprecision(5)<<r_rmax + <<std::setw(12)<<std::setprecision(5)<<m_r_rmax <<" |"<<std::endl; out<<"| radius step | " - <<std::setw(12)<<std::setprecision(5)<<r_rstep + <<std::setw(12)<<std::setprecision(5)<<m_r_rstep <<" |"<<std::endl; out<<"| min Z-vertex position | " <<std::setw(12)<<std::setprecision(5)<<m_zmin @@ -641,7 +641,7 @@ MsgStream& InDet::SiSpacePointsSeedMaker_LowMomentum::dumpEvent( MsgStream& out <<std::setw(12)<<m_nsaz <<" |"<<std::endl; out<<"| seeds | " - <<std::setw(12)<<l_seeds.size() + <<std::setw(12)<<m_l_seeds.size() <<" |"<<std::endl; out<<"|---------------------------------------------------------------------|" <<std::endl; @@ -716,7 +716,7 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::findNext () { if(m_endlist) return; - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); if (m_mode==0 || m_mode==1) production2Sp (); else if(m_mode==2 || m_mode==3) { @@ -729,7 +729,7 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::findNext () } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); ++m_nlist; } @@ -739,16 +739,16 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::findNext () bool InDet::SiSpacePointsSeedMaker_LowMomentum::newVertices(const std::list<Trk::Vertex>& lV) { - unsigned int s1 = l_vertex.size(); + unsigned int s1 = m_l_vertex.size(); unsigned int s2 = lV .size(); if(s1==0 && s2==0) return false; std::list<Trk::Vertex>::const_iterator v; - l_vertex.erase(l_vertex.begin(),l_vertex.end()); + m_l_vertex.erase(m_l_vertex.begin(),m_l_vertex.end()); for(v=lV.begin(); v!=lV.end(); ++v) { - l_vertex.push_back(float((*v).position().z())); + m_l_vertex.push_back(float((*v).position().z())); } return false; } @@ -765,17 +765,17 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::buildFrameWork() m_rapcut = fabs(m_rapcut) ; m_dzdrmax = 1./tan(2.*atan(exp(-m_rapcut))); m_dzdrmin =-m_dzdrmax ; - m_r3max = r_rmax ; + m_r3max = m_r_rmax ; m_ns = m_nsaz = m_nr = m_nrfz = 0; // Build radius sorted containers // - r_size = int((r_rmax+.1)/r_rstep); - r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[r_size]; - r_index = new int[r_size]; - r_map = new int[r_size]; - m_nr = 0; for(int i=0; i!=r_size; ++i) {r_index[i]=0; r_map[i]=0;} + m_r_size = int((m_r_rmax+.1)/m_r_rstep); + m_r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[m_r_size]; + m_r_index = new int[m_r_size]; + m_r_map = new int[m_r_size]; + m_nr = 0; for(int i=0; i!=m_r_size; ++i) {m_r_index[i]=0; m_r_map[i]=0;} // Build radius-azimuthal sorted containers // @@ -790,7 +790,7 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::buildFrameWork() // Build radius-azimuthal-Z sorted containers // - m_nrfz = 0; for(int i=0; i!=220; ++i) {rfz_index [i]=0; rfz_map [i]=0;} + m_nrfz = 0; for(int i=0; i!=220; ++i) {m_rfz_index [i]=0; m_rfz_map [i]=0;} // Build maps for radius-azimuthal-Z sorted collections @@ -807,62 +807,62 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::buildFrameWork() int a = f *11+z; int b = fb*11+z; int c = ft*11+z; - rfz_b [a] = 3; rfz_t [a] = 3; - rfz_ib[a][0] = a; rfz_it[a][0] = a; - rfz_ib[a][1] = b; rfz_it[a][1] = b; - rfz_ib[a][2] = c; rfz_it[a][2] = c; + m_rfz_b [a] = 3; m_rfz_t [a] = 3; + m_rfz_ib[a][0] = a; m_rfz_it[a][0] = a; + m_rfz_ib[a][1] = b; m_rfz_it[a][1] = b; + m_rfz_ib[a][2] = c; m_rfz_it[a][2] = c; if (z==5) { - rfz_t [a] = 9 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; - rfz_it[a][6] = a-1; - rfz_it[a][7] = b-1; - rfz_it[a][8] = c-1; + m_rfz_t [a] = 9 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; + m_rfz_it[a][6] = a-1; + m_rfz_it[a][7] = b-1; + m_rfz_it[a][8] = c-1; } else if(z> 5) { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a-1; - rfz_ib[a][4] = b-1; - rfz_ib[a][5] = c-1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a-1; + m_rfz_ib[a][4] = b-1; + m_rfz_ib[a][5] = c-1; if(z<10) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; } } else { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a+1; - rfz_ib[a][4] = b+1; - rfz_ib[a][5] = c+1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a+1; + m_rfz_ib[a][4] = b+1; + m_rfz_ib[a][5] = c+1; if(z>0) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a-1; - rfz_it[a][4] = b-1; - rfz_it[a][5] = c-1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a-1; + m_rfz_it[a][4] = b-1; + m_rfz_it[a][5] = c-1; } } if (z==3) { - rfz_b[a] = 9; - rfz_ib[a][6] = a+2; - rfz_ib[a][7] = b+2; - rfz_ib[a][8] = c+2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a+2; + m_rfz_ib[a][7] = b+2; + m_rfz_ib[a][8] = c+2; } else if(z==7) { - rfz_b[a] = 9; - rfz_ib[a][6] = a-2; - rfz_ib[a][7] = b-2; - rfz_ib[a][8] = c-2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a-2; + m_rfz_ib[a][7] = b-2; + m_rfz_ib[a][8] = c-2; } } } @@ -876,8 +876,8 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::buildFrameWork() if(!m_Zo) m_Zo = new float[m_maxsizeSP]; if(!m_OneSeeds) m_OneSeeds = new InDet::SiSpacePointsSeed [m_maxOneSize]; - i_seed = l_seeds.begin(); - i_seede = l_seeds.end (); + m_i_seed = m_l_seeds.begin(); + m_i_seede = m_l_seeds.end (); } /////////////////////////////////////////////////////////////////// @@ -937,11 +937,11 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::fillLists() const float pi2 = 2.*M_PI; std::list<InDet::SiSpacePointForSeed*>::iterator r; - for(int i=0; i!= r_size; ++i) { + for(int i=0; i!= m_r_size; ++i) { - if(!r_map[i]) continue; r = r_Sorted[i].begin(); + if(!m_r_map[i]) continue; r = m_r_Sorted[i].begin(); - while(r!=r_Sorted[i].end()) { + while(r!=m_r_Sorted[i].end()) { // Azimuthal angle sort // @@ -960,10 +960,10 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::fillLists() Z>-250.?z=5:Z>-450.?z=4:Z>-925.?z=3:Z>-1400.?z=2:Z>-2500.?z=1:z= 0; } int n = f*11+z; ++m_nsaz; - rfz_Sorted[n].push_back(*r); if(!rfz_map[n]++) rfz_index[m_nrfz++] = n; - r_Sorted[i].erase(r++); + m_rfz_Sorted[n].push_back(*r); if(!m_rfz_map[n]++) m_rfz_index[m_nrfz++] = n; + m_r_Sorted[i].erase(r++); } - r_map[i] = 0; + m_r_map[i] = 0; } m_nr = 0; m_state = 0; @@ -976,13 +976,13 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::fillLists() void InDet::SiSpacePointsSeedMaker_LowMomentum::erase() { for(int i=0; i!=m_nr; ++i) { - int n = r_index[i]; r_map[n] = 0; - r_Sorted[n].clear(); + int n = m_r_index[i]; m_r_map[n] = 0; + m_r_Sorted[n].clear(); } for(int i=0; i!=m_nrfz; ++i) { - int n = rfz_index[i]; rfz_map[n] = 0; - rfz_Sorted[n].clear(); + int n = m_rfz_index[i]; m_rfz_map[n] = 0; + m_rfz_Sorted[n].clear(); } m_state = 0; @@ -1032,19 +1032,19 @@ void InDet::SiSpacePointsSeedMaker_LowMomentum::production3Sp() for(; z!=11; ++z) { - int a = f *11+ZI[z]; if(!rfz_map[a]) continue; + int a = f *11+ZI[z]; if(!m_rfz_map[a]) continue; int NB = 0, NT = 0; - for(int i=0; i!=rfz_b[a]; ++i) { + for(int i=0; i!=m_rfz_b[a]; ++i) { - int an = rfz_ib[a][i]; - if(!rfz_map[an]) continue; - rb [NB] = rfz_Sorted[an].begin(); rbe[NB++] = rfz_Sorted[an].end(); + int an = m_rfz_ib[a][i]; + if(!m_rfz_map[an]) continue; + rb [NB] = m_rfz_Sorted[an].begin(); rbe[NB++] = m_rfz_Sorted[an].end(); } - for(int i=0; i!=rfz_t[a]; ++i) { + for(int i=0; i!=m_rfz_t[a]; ++i) { - int an = rfz_it[a][i]; - if(!rfz_map[an]) continue; - rt [NT] = rfz_Sorted[an].begin(); rte[NT++] = rfz_Sorted[an].end(); + int an = m_rfz_it[a][i]; + if(!m_rfz_map[an]) continue; + rt [NT] = m_rfz_Sorted[an].begin(); rte[NT++] = m_rfz_Sorted[an].end(); } production3Sp(rb,rbe,rt,rte,NB,NT,nseed,K); if(!m_endlist) {m_fNmin=f; m_zMin = z; return;} diff --git a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_Trigger.cxx b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_Trigger.cxx index 3bc840198c2ce51d25ed8402f01a7f3a77715da8..aefda1caf4b73f16b646e8bb4a5a98eef6cb9d12 100755 --- a/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_Trigger.cxx +++ b/InnerDetector/InDetRecTools/SiSpacePointsSeedTool_xk/src/SiSpacePointsSeedMaker_Trigger.cxx @@ -59,11 +59,11 @@ InDet::SiSpacePointsSeedMaker_Trigger::SiSpacePointsSeedMaker_Trigger m_diverpps = 1.7 ; m_diversss = 1000. ; m_dazmax = .02 ; - r_rmax = 600. ; - r_rstep = 2. ; - r_Sorted = 0 ; - r_index = 0 ; - r_map = 0 ; + m_r_rmax = 600. ; + m_r_rstep = 2. ; + m_r_Sorted = 0 ; + m_r_index = 0 ; + m_r_map = 0 ; m_maxsizeSP = 1500 ; m_maxOneSize= 3 ; m_SP = 0 ; @@ -93,8 +93,8 @@ InDet::SiSpacePointsSeedMaker_Trigger::SiSpacePointsSeedMaker_Trigger declareProperty("usePixel" ,m_pixel ); declareProperty("useSCT" ,m_sct ); declareProperty("pTmin" ,m_ptmin ); - declareProperty("radMax" ,r_rmax ); - declareProperty("radStep" ,r_rstep ); + declareProperty("radMax" ,m_r_rmax ); + declareProperty("radStep" ,m_r_rstep ); declareProperty("maxSize" ,m_maxsize ); declareProperty("maxSizeSP" ,m_maxsizeSP ); declareProperty("minZ" ,m_zmin ); @@ -133,20 +133,20 @@ InDet::SiSpacePointsSeedMaker_Trigger::SiSpacePointsSeedMaker_Trigger InDet::SiSpacePointsSeedMaker_Trigger::~SiSpacePointsSeedMaker_Trigger() { - if(r_index ) delete [] r_index ; - if(r_map ) delete [] r_map ; - if(r_Sorted) delete [] r_Sorted; + if(m_r_index ) delete [] m_r_index ; + if(m_r_map ) delete [] m_r_map ; + if(m_r_Sorted) delete [] m_r_Sorted; // Delete seeds // - for(i_seed=l_seeds.begin(); i_seed!=l_seeds.end (); ++i_seed) { - delete *i_seed; + for(m_i_seed=m_l_seeds.begin(); m_i_seed!=m_l_seeds.end (); ++m_i_seed) { + delete *m_i_seed; } // Delete space points for reconstruction // - i_spforseed=l_spforseed.begin(); - for(; i_spforseed!=l_spforseed.end(); ++i_spforseed) { - delete *i_spforseed; + m_i_spforseed=m_l_spforseed.begin(); + for(; m_i_spforseed!=m_l_spforseed.end(); ++m_i_spforseed) { + delete *m_i_spforseed; } if(m_SP) delete [] m_SP; @@ -222,10 +222,10 @@ void InDet::SiSpacePointsSeedMaker_Trigger::newEvent (int) } else m_K = 2./(300.* 5. ); - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; - int irmax = r_size-1 ; + float irstep = 1./m_r_rstep; + int irmax = m_r_size-1 ; // Get pixels space points containers from store gate // @@ -246,13 +246,13 @@ void InDet::SiSpacePointsSeedMaker_Trigger::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -278,13 +278,13 @@ void InDet::SiSpacePointsSeedMaker_Trigger::newEvent (int) for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -304,13 +304,13 @@ void InDet::SiSpacePointsSeedMaker_Trigger::newEvent (int) for (; sp!=spe; ++sp) { float r = (*sp)->r(); - if(r<0. || r>=r_rmax) continue; + if(r<0. || r>=m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -337,10 +337,10 @@ void InDet::SiSpacePointsSeedMaker_Trigger::newRegion } else m_K = 2./(300.* 5. ); - i_spforseed = l_spforseed.begin(); + m_i_spforseed = m_l_spforseed.begin(); - float irstep = 1./r_rstep; - int irmax = r_size-1 ; + float irstep = 1./m_r_rstep; + int irmax = m_r_size-1 ; // Get pixels space points containers from store gate // @@ -365,13 +365,13 @@ void InDet::SiSpacePointsSeedMaker_Trigger::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -401,13 +401,13 @@ void InDet::SiSpacePointsSeedMaker_Trigger::newRegion for(; sp != spe; ++sp) { - float r = (*sp)->r(); if(r<0. || r>=r_rmax) continue; + float r = (*sp)->r(); if(r<0. || r>=m_r_rmax) continue; InDet::SiSpacePointForSeed* sps = newSpacePoint((*sp)); int ir = int(sps->radius()*irstep); if(ir>irmax) ir = irmax; - r_Sorted[ir].push_back(sps); ++r_map[ir]; - if(r_map[ir]==1) r_index[m_nr++] = ir; + m_r_Sorted[ir].push_back(sps); ++m_r_map[ir]; + if(m_r_map[ir]==1) m_r_index[m_nr++] = ir; ++m_ns; } } @@ -453,7 +453,7 @@ void InDet::SiSpacePointsSeedMaker_Trigger::find2Sp(const std::list<Trk::Vertex> if(newv || !m_state || m_nspoint!=2 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); m_state = 1 ; m_nspoint = 2 ; m_nlist = 0 ; @@ -464,7 +464,7 @@ void InDet::SiSpacePointsSeedMaker_Trigger::find2Sp(const std::list<Trk::Vertex> m_zMin = 0 ; production2Sp (); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); if(m_outputlevel<=0) { m_nprint=1; msg(MSG::DEBUG)<<(*this)<<endmsg; @@ -484,7 +484,7 @@ void InDet::SiSpacePointsSeedMaker_Trigger::find3Sp(const std::list<Trk::Vertex> if(newv || !m_state || m_nspoint!=3 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 3 ; m_nlist = 0 ; @@ -495,7 +495,7 @@ void InDet::SiSpacePointsSeedMaker_Trigger::find3Sp(const std::list<Trk::Vertex> m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_mapSeeds.begin(); m_seede = m_mapSeeds.end (); @@ -523,7 +523,7 @@ void InDet::SiSpacePointsSeedMaker_Trigger::findVSp (const std::list<Trk::Vertex if(newv || !m_state || m_nspoint!=4 || m_mode!=mode || m_nlist) { - i_seede = l_seeds.begin() ; + m_i_seede = m_l_seeds.begin() ; m_state = 1 ; m_nspoint = 4 ; m_nlist = 0 ; @@ -534,7 +534,7 @@ void InDet::SiSpacePointsSeedMaker_Trigger::findVSp (const std::list<Trk::Vertex m_zMin = 0 ; production3Sp(); } - i_seed = l_seeds.begin (); + m_i_seed = m_l_seeds.begin (); m_seed = m_mapSeeds.begin(); m_seede = m_mapSeeds.end (); @@ -597,10 +597,10 @@ MsgStream& InDet::SiSpacePointsSeedMaker_Trigger::dumpConditions( MsgStream& out <<std::setw(12)<<std::setprecision(5)<<m_rapcut <<" |"<<std::endl; out<<"| max radius SP | " - <<std::setw(12)<<std::setprecision(5)<<r_rmax + <<std::setw(12)<<std::setprecision(5)<<m_r_rmax <<" |"<<std::endl; out<<"| radius step | " - <<std::setw(12)<<std::setprecision(5)<<r_rstep + <<std::setw(12)<<std::setprecision(5)<<m_r_rstep <<" |"<<std::endl; out<<"| min Z-vertex position | " <<std::setw(12)<<std::setprecision(5)<<m_zmin @@ -709,7 +709,7 @@ MsgStream& InDet::SiSpacePointsSeedMaker_Trigger::dumpEvent( MsgStream& out ) co <<std::setw(12)<<m_nsazv <<" |"<<std::endl; out<<"| seeds | " - <<std::setw(12)<<l_seeds.size() + <<std::setw(12)<<m_l_seeds.size() <<" |"<<std::endl; out<<"|---------------------------------------------------------------------|" <<std::endl; @@ -784,13 +784,13 @@ void InDet::SiSpacePointsSeedMaker_Trigger::findNext () { if(m_endlist) return; - i_seede = l_seeds.begin(); + m_i_seede = m_l_seeds.begin(); if (m_mode==0 || m_mode==1) production2Sp (); else if(m_mode==2 || m_mode==3) production3Sp (); else if(m_mode==5 || m_mode==6) production3Sp (); - i_seed = l_seeds.begin(); + m_i_seed = m_l_seeds.begin(); m_seed = m_mapSeeds.begin(); m_seede = m_mapSeeds.end (); ++m_nlist; @@ -802,17 +802,17 @@ void InDet::SiSpacePointsSeedMaker_Trigger::findNext () bool InDet::SiSpacePointsSeedMaker_Trigger::newVertices(const std::list<Trk::Vertex>& lV) { - unsigned int s1 = l_vertex.size(); + unsigned int s1 = m_l_vertex.size(); unsigned int s2 = lV .size(); if(s1==0 && s2==0) return false; std::list<Trk::Vertex>::const_iterator v; - l_vertex.erase(l_vertex.begin(),l_vertex.end()); + m_l_vertex.erase(m_l_vertex.begin(),m_l_vertex.end()); for(v=lV.begin(); v!=lV.end(); ++v) { - l_vertex.push_back(float((*v).position().z())); - if(l_vertex.size() >= m_maxNumberVertices) break; + m_l_vertex.push_back(float((*v).position().z())); + if(m_l_vertex.size() >= m_maxNumberVertices) break; } return false; } @@ -827,7 +827,7 @@ void InDet::SiSpacePointsSeedMaker_Trigger::buildFrameWork() m_rapcut = fabs(m_rapcut) ; m_dzdrmax = 1./tan(2.*atan(exp(-m_rapcut))); m_dzdrmin =-m_dzdrmax ; - m_r3max = r_rmax ; + m_r3max = m_r_rmax ; m_COF = 134*.05*9. ; m_ipt = 1./fabs(.9*m_ptmin) ; m_ipt2 = m_ipt*m_ipt ; @@ -837,26 +837,26 @@ void InDet::SiSpacePointsSeedMaker_Trigger::buildFrameWork() // Build radius sorted containers // - r_size = int((r_rmax+.1)/r_rstep); - r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[r_size]; - r_index = new int[r_size]; - r_map = new int[r_size]; - m_nr = 0; for(int i=0; i!=r_size; ++i) {r_index[i]=0; r_map[i]=0;} + m_r_size = int((m_r_rmax+.1)/m_r_rstep); + m_r_Sorted = new std::list<InDet::SiSpacePointForSeed*>[m_r_size]; + m_r_index = new int[m_r_size]; + m_r_map = new int[m_r_size]; + m_nr = 0; for(int i=0; i!=m_r_size; ++i) {m_r_index[i]=0; m_r_map[i]=0;} // Build radius-azimuthal sorted containers // const float pi2 = 2.*M_PI ; const int NFmax = 53 ; const float sFmax = float(NFmax )/pi2; - const float m_sFmin = 100./60. ; + const float sFmin = 100./60. ; - m_sF = m_ptmin /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < m_sFmin) m_sF = m_sFmin; + m_sF = m_ptmin /60. ; if(m_sF >sFmax ) m_sF = sFmax ; else if(m_sF < sFmin) m_sF = sFmin; m_fNmax = int(pi2*m_sF); if(m_fNmax >=NFmax) m_fNmax = NFmax-1; // Build radius-azimuthal-Z sorted containers // - m_nrfz = 0; for(int i=0; i!=583; ++i) {rfz_index [i]=0; rfz_map [i]=0;} + m_nrfz = 0; for(int i=0; i!=583; ++i) {m_rfz_index [i]=0; m_rfz_map [i]=0;} // Build radius-azimuthal-Z sorted containers for Z-vertices // @@ -864,7 +864,7 @@ void InDet::SiSpacePointsSeedMaker_Trigger::buildFrameWork() const float sFvmax = float(NFtmax)/pi2; m_sFv = m_ptmin/120. ; if(m_sFv >sFvmax) m_sFv = sFvmax; m_fvNmax = int(pi2*m_sFv); if(m_fvNmax>=NFtmax) m_fvNmax = NFtmax-1; - m_nrfzv = 0; for(int i=0; i!=300; ++i) {rfzv_index[i]=0; rfzv_map[i]=0;} + m_nrfzv = 0; for(int i=0; i!=300; ++i) {m_rfzv_index[i]=0; m_rfzv_map[i]=0;} // Build maps for radius-azimuthal-Z sorted collections // @@ -880,62 +880,62 @@ void InDet::SiSpacePointsSeedMaker_Trigger::buildFrameWork() int a = f *11+z; int b = fb*11+z; int c = ft*11+z; - rfz_b [a] = 3; rfz_t [a] = 3; - rfz_ib[a][0] = a; rfz_it[a][0] = a; - rfz_ib[a][1] = b; rfz_it[a][1] = b; - rfz_ib[a][2] = c; rfz_it[a][2] = c; + m_rfz_b [a] = 3; m_rfz_t [a] = 3; + m_rfz_ib[a][0] = a; m_rfz_it[a][0] = a; + m_rfz_ib[a][1] = b; m_rfz_it[a][1] = b; + m_rfz_ib[a][2] = c; m_rfz_it[a][2] = c; if (z==5) { - rfz_t [a] = 9 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; - rfz_it[a][6] = a-1; - rfz_it[a][7] = b-1; - rfz_it[a][8] = c-1; + m_rfz_t [a] = 9 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; + m_rfz_it[a][6] = a-1; + m_rfz_it[a][7] = b-1; + m_rfz_it[a][8] = c-1; } else if(z> 5) { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a-1; - rfz_ib[a][4] = b-1; - rfz_ib[a][5] = c-1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a-1; + m_rfz_ib[a][4] = b-1; + m_rfz_ib[a][5] = c-1; if(z<10) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a+1; - rfz_it[a][4] = b+1; - rfz_it[a][5] = c+1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a+1; + m_rfz_it[a][4] = b+1; + m_rfz_it[a][5] = c+1; } } else { - rfz_b [a] = 6 ; - rfz_ib[a][3] = a+1; - rfz_ib[a][4] = b+1; - rfz_ib[a][5] = c+1; + m_rfz_b [a] = 6 ; + m_rfz_ib[a][3] = a+1; + m_rfz_ib[a][4] = b+1; + m_rfz_ib[a][5] = c+1; if(z>0) { - rfz_t [a] = 6 ; - rfz_it[a][3] = a-1; - rfz_it[a][4] = b-1; - rfz_it[a][5] = c-1; + m_rfz_t [a] = 6 ; + m_rfz_it[a][3] = a-1; + m_rfz_it[a][4] = b-1; + m_rfz_it[a][5] = c-1; } } if (z==3) { - rfz_b[a] = 9; - rfz_ib[a][6] = a+2; - rfz_ib[a][7] = b+2; - rfz_ib[a][8] = c+2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a+2; + m_rfz_ib[a][7] = b+2; + m_rfz_ib[a][8] = c+2; } else if(z==7) { - rfz_b[a] = 9; - rfz_ib[a][6] = a-2; - rfz_ib[a][7] = b-2; - rfz_ib[a][8] = c-2; + m_rfz_b[a] = 9; + m_rfz_ib[a][6] = a-2; + m_rfz_ib[a][7] = b-2; + m_rfz_ib[a][8] = c-2; } } } @@ -954,21 +954,21 @@ void InDet::SiSpacePointsSeedMaker_Trigger::buildFrameWork() int a = f *3+z; int b = fb*3+z; int c = ft*3+z; - rfzv_n[a] = 3; - rfzv_i[a][0] = a; - rfzv_i[a][1] = b; - rfzv_i[a][2] = c; + m_rfzv_n[a] = 3; + m_rfzv_i[a][0] = a; + m_rfzv_i[a][1] = b; + m_rfzv_i[a][2] = c; if (z>1) { - rfzv_n[a] = 6; - rfzv_i[a][3] = a-1; - rfzv_i[a][4] = b-1; - rfzv_i[a][5] = c-1; + m_rfzv_n[a] = 6; + m_rfzv_i[a][3] = a-1; + m_rfzv_i[a][4] = b-1; + m_rfzv_i[a][5] = c-1; } else if(z<1) { - rfzv_n[a] = 6; - rfzv_i[a][3] = a+1; - rfzv_i[a][4] = b+1; - rfzv_i[a][5] = c+1; + m_rfzv_n[a] = 6; + m_rfzv_i[a][3] = a+1; + m_rfzv_i[a][4] = b+1; + m_rfzv_i[a][5] = c+1; } } } @@ -982,8 +982,8 @@ void InDet::SiSpacePointsSeedMaker_Trigger::buildFrameWork() if(!m_Zo) m_Zo = new float[m_maxsizeSP]; if(!m_OneSeeds) m_OneSeeds = new InDet::SiSpacePointsSeed [m_maxOneSize]; - i_seed = l_seeds.begin(); - i_seede = l_seeds.end (); + m_i_seed = m_l_seeds.begin(); + m_i_seede = m_l_seeds.end (); m_seed = m_mapSeeds.begin(); m_seede = m_mapSeeds.end (); } @@ -1044,11 +1044,11 @@ void InDet::SiSpacePointsSeedMaker_Trigger::fillLists() const float pi2 = 2.*M_PI; std::list<InDet::SiSpacePointForSeed*>::iterator r; - for(int i=0; i!= r_size; ++i) { + for(int i=0; i!= m_r_size; ++i) { - if(!r_map[i]) continue; r = r_Sorted[i].begin(); + if(!m_r_map[i]) continue; r = m_r_Sorted[i].begin(); - while(r!=r_Sorted[i].end()) { + while(r!=m_r_Sorted[i].end()) { // Azimuthal angle sort // @@ -1063,7 +1063,7 @@ void InDet::SiSpacePointsSeedMaker_Trigger::fillLists() if(Z>0.) {Z < 250.? z=5 : z=6;} else {Z >-250.? z=5 : z=4;} int n = f*11+z; ++m_nsaz; - rfz_Sorted[n].push_back(*r); if(!rfz_map[n]++) rfz_index[m_nrfz++] = n; + m_rfz_Sorted[n].push_back(*r); if(!m_rfz_map[n]++) m_rfz_index[m_nrfz++] = n; if((*r)->spacepoint->clusterList().second == 0 && z>=3 && z<=7) { z<=4 ? z=0 : z>=6 ? z=2 : z=1; @@ -1072,11 +1072,11 @@ void InDet::SiSpacePointsSeedMaker_Trigger::fillLists() // f = int(F*m_sFv); f<0 ? f+=m_fvNmax : f> m_fvNmax ? f-=m_fvNmax : f=f; n = f*3+z; ++m_nsazv; - rfzv_Sorted[n].push_back(*r); if(!rfzv_map[n]++) rfzv_index[m_nrfzv++] = n; + m_rfzv_Sorted[n].push_back(*r); if(!m_rfzv_map[n]++) m_rfzv_index[m_nrfzv++] = n; } - r_Sorted[i].erase(r++); + m_r_Sorted[i].erase(r++); } - r_map[i] = 0; + m_r_map[i] = 0; } m_nr = 0; m_state = 0; @@ -1089,18 +1089,18 @@ void InDet::SiSpacePointsSeedMaker_Trigger::fillLists() void InDet::SiSpacePointsSeedMaker_Trigger::erase() { for(int i=0; i!=m_nr; ++i) { - int n = r_index[i]; r_map[n] = 0; - r_Sorted[n].clear(); + int n = m_r_index[i]; m_r_map[n] = 0; + m_r_Sorted[n].clear(); } for(int i=0; i!=m_nrfz; ++i) { - int n = rfz_index[i]; rfz_map[n] = 0; - rfz_Sorted[n].clear(); + int n = m_rfz_index[i]; m_rfz_map[n] = 0; + m_rfz_Sorted[n].clear(); } for(int i=0; i!=m_nrfzv; ++i) { - int n = rfzv_index[i]; rfzv_map[n] = 0; - rfzv_Sorted[n].clear(); + int n = m_rfzv_index[i]; m_rfzv_map[n] = 0; + m_rfzv_Sorted[n].clear(); } m_state = 0; m_ns = 0; @@ -1131,9 +1131,9 @@ void InDet::SiSpacePointsSeedMaker_Trigger::production2Sp() int z = 0; if(!m_endlist) z = m_zMin; for(; z!=3; ++z) { - int a = f*3+z; if(!rfzv_map[a]) continue; - r0 = rfzv_Sorted[a].begin(); - r0e = rfzv_Sorted[a].end (); + int a = f*3+z; if(!m_rfzv_map[a]) continue; + r0 = m_rfzv_Sorted[a].begin(); + r0e = m_rfzv_Sorted[a].end (); if(!m_endlist) {r0 = m_rMin; m_endlist = true;} @@ -1151,14 +1151,14 @@ void InDet::SiSpacePointsSeedMaker_Trigger::production2Sp() // Bottom links production // - int NB = rfzv_n[a]; + int NB = m_rfzv_n[a]; for(int i=0; i!=NB; ++i) { - int an = rfzv_i[a][i]; - if(!rfzv_map[an]) continue; + int an = m_rfzv_i[a][i]; + if(!m_rfzv_map[an]) continue; - r = rfzv_Sorted[an].begin(); - re = rfzv_Sorted[an].end (); + r = m_rfzv_Sorted[an].begin(); + re = m_rfzv_Sorted[an].end (); for(; r!=re; ++r) { @@ -1222,19 +1222,19 @@ void InDet::SiSpacePointsSeedMaker_Trigger::production3Sp() for(; z!=11; ++z) { - int a = f *11+ZI[z]; if(!rfz_map[a]) continue; + int a = f *11+ZI[z]; if(!m_rfz_map[a]) continue; int NB = 0, NT = 0; - for(int i=0; i!=rfz_b[a]; ++i) { + for(int i=0; i!=m_rfz_b[a]; ++i) { - int an = rfz_ib[a][i]; - if(!rfz_map[an]) continue; - rb [NB] = rfz_Sorted[an].begin(); rbe[NB++] = rfz_Sorted[an].end(); + int an = m_rfz_ib[a][i]; + if(!m_rfz_map[an]) continue; + rb [NB] = m_rfz_Sorted[an].begin(); rbe[NB++] = m_rfz_Sorted[an].end(); } - for(int i=0; i!=rfz_t[a]; ++i) { + for(int i=0; i!=m_rfz_t[a]; ++i) { - int an = rfz_it[a][i]; - if(!rfz_map[an]) continue; - rt [NT] = rfz_Sorted[an].begin(); rte[NT++] = rfz_Sorted[an].end(); + int an = m_rfz_it[a][i]; + if(!m_rfz_map[an]) continue; + rt [NT] = m_rfz_Sorted[an].begin(); rte[NT++] = m_rfz_Sorted[an].end(); } if(!m_trigger) production3Sp (rb,rbe,rt,rte,NB,NT,nseed); else production3SpTrigger(rb,rbe,rt,rte,NB,NT,nseed);