87 {
89
90 std::vector<const xAOD::TrackParticle*> tmpp(particles);
93 if( tst != tmpp.end()) {
95 return {nullptr};
96 }
97
98 static const SG::ConstAccessor<float> vyAcc("vy");
99 bool fullxAOD=false; if(vyAcc.isAvailable(*particles[0])) fullxAOD=true;
100
102 BEAM.makePrivateStore();
104 std::vector<float> defaultCovar {0.015f*0.015f,0.f,0.015f*0.015f,0.f,0.f,1.e6};
106 float beamtiltX=0.;
107 float beamtiltY=0.;
109 beamcov.setIdentity();
110
112 BEAM.
setX(eventINFO->beamPosX());
113 BEAM.
setY(eventINFO->beamPosY());
114 BEAM.
setZ(eventINFO->beamPosZ());
115 beamcov(0,0) = eventINFO->beamPosSigmaX() * eventINFO->beamPosSigmaX();
116 beamcov(1,0) = beamcov(0,1) = eventINFO->beamPosSigmaXY();
117 beamcov(1,1) = eventINFO->beamPosSigmaY() * eventINFO->beamPosSigmaY();
118 beamcov(2,2) = eventINFO->beamPosSigmaZ() * eventINFO->beamPosSigmaZ()*1.e6;
120 beamtiltX= eventINFO->beamTiltXZ();
121 beamtiltY= eventINFO->beamTiltYZ();
122 }
123 SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle {
m_beamSpotKey, ctx };
124
125 if(beamSpotHandle.
isValid() && fullxAOD ){
127 BEAM.
setPosition(beamSpotHandle->beamVtx().position());
128 beamcov = beamSpotHandle->beamVtx().covariancePosition();
129 beamcov(2,2) *= 1.e6;
131 beamtiltX= beamSpotHandle->beamTilt(0);
132 beamtiltY= beamSpotHandle->beamTilt(1);
133 }
140
141 if(fullxAOD){
142
143 if( beamtiltX!=0. || beamtiltY!=0.) {
144 std::unique_ptr<xAOD::Vertex> iniVertex =
m_fitterSvc->fit(ctx,particles,BEAM);
145 if(!iniVertex) return {nullptr};
146 BEAM.
setX(BEAM.
x()+beamtiltX*iniVertex->z());
147 BEAM.
setY(BEAM.
y()+beamtiltY*iniVertex->z());
148 }
150 }
151
152
153
154
156 std::unique_ptr< SG::AuxStoreInternal > pAux;
158 std::vector<const xAOD::TrackParticle*> wrkTrkC(NPRT);
159 pAux = std::make_unique< SG::AuxStoreInternal >();
160 TPC.setStore( pAux.get() );
162 static const SG::ConstAccessor<float> vzAcc("vz");
163 for(
int i=0;
i<NPRT;
i++){
165 if(!TPC[i])return {nullptr};
166 const float mvx= (eventINFO) ? eventINFO->beamPosX() : 0.;
167 const float mvy= (eventINFO) ? eventINFO->beamPosY() : 0.;
168 const float mvz= (vzAcc.isAvailable(*particles[i])) ?
particles[
i]->vz() : 0.;
169 TPC[
i]->setParametersOrigin( mvx, mvy, mvz);
171 }
172
173
174 if( beamtiltX!=0. || beamtiltY!=0.) {
175 std::unique_ptr<xAOD::Vertex> iniVertex =
m_fitterSvc->fit(ctx,wrkTrkC,BEAM);
176 if(!iniVertex) return {nullptr};
177 BEAM.
setX(BEAM.
x()+beamtiltX*iniVertex->z());
178 BEAM.
setY(BEAM.
y()+beamtiltY*iniVertex->z());
179 }
181 }
#define AmgSymMatrix(dim)
float z() const
Returns the z position.
void setCovariancePosition(const AmgSymMatrix(3)&covariancePosition)
Sets the vertex covariance matrix.
void setZ(float value)
Sets the z position.
void setX(float value)
Sets the x position.
void setY(float value)
Sets the y position.
float y() const
Returns the y position.
void setCovariance(const std::vector< float > &value)
Sets the covariance matrix as a simple vector of values.
void setPosition(const Amg::Vector3D &position)
Sets the 3-position.
float x() const
Returns the x position.
const std::vector< float > & covariance() const
Returns the covariance matrix as a simple vector of values.
DataModel_detail::iterator< DVL > unique(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of unique for DataVector/List.
void sort(typename DataModel_detail::iterator< DVL > beg, typename DataModel_detail::iterator< DVL > end)
Specialization of sort for DataVector/List.
Vertex_v1 Vertex
Define the latest version of the vertex class.
1.16.1