124 if (recoBeamSpotHandle.
isValid()) {
127 edm::LogError(
"UnusableBeamSpot") <<
"No beam spot available from EventSetup";
132 if ((beamVertexState.error().cxx() <= 0.) || (beamVertexState.error().cyy() <= 0.) ||
133 (beamVertexState.error().czz() <= 0.)) {
135 edm::LogError(
"UnusableBeamSpot") <<
"Beamspot with invalid errors " << beamVertexState.error().matrix();
142 for (
auto const& old : oldVertices) {
143 if (!(old.isFake())) {
146 auto result = std::make_unique<reco::VertexCollection>();
161 std::vector<reco::TransientTrack> t_tks;
170 t_tks = (*theB).build(tks,
beamSpot);
173 std::cout <<
"RecoVertex/PrimaryVertexProducer" 174 <<
"Found: " << t_tks.size() <<
" reconstructed tracks" 185 std::cout <<
" clustering returned " <<
clusters.size() <<
" clusters from " << seltks.size()
186 <<
" selected tracks" << std::endl;
191 auto result = std::make_unique<reco::VertexCollection>();
194 std::vector<TransientVertex>
pvs;
195 for (
std::vector<std::vector<reco::TransientTrack> >::const_iterator iclus =
clusters.begin();
201 double meantime = 0.;
204 for (
const auto& tk : *iclus) {
205 const double time = tk.timeExt();
206 const double err = tk.dtErrorExt();
207 const double inverr =
err > 0. ? 1.0 /
err : 0.;
208 const double w = inverr * inverr;
213 meantime = sumwt / sumw;
214 double sumsq = sumwt2 - sumwt * sumwt / sumw;
215 double chisq = iclus->size() > 1 ? sumsq / double(iclus->size() - 1) : sumsq /
double(iclus->size());
216 vartime = chisq / sumw;
220 if (
algorithm->useBeamConstraint && validBS && (iclus->size() > 1)) {
222 }
else if (!(
algorithm->useBeamConstraint) && (iclus->size() > 1)) {
227 if (
f4D and
v.isValid()) {
228 auto err =
v.positionError().matrix4D();
230 auto trkWeightMap3d =
v.weightMap();
231 v =
TransientVertex(
v.position(), meantime,
err,
v.originalTracks(),
v.totalChiSquared(),
v.degreesOfFreedom());
232 v.weightMap(trkWeightMap3d);
240 std::cout <<
"=" <<
v.position().x() <<
" " <<
v.position().y() <<
" " <<
v.position().z();
243 std::cout <<
" cluster size = " << (*iclus).size() << std::endl;
245 std::cout <<
"Invalid fitted vertex, cluster size=" << (*iclus).size() << std::endl;
249 if (
v.isValid() && (
v.degreesOfFreedom() >=
algorithm->minNdof) &&
250 (!validBS || (*(
algorithm->vertexSelector))(
v, beamVertexState)))
255 std::cout <<
"PrimaryVertexProducerAlgorithm::vertices candidates =" <<
pvs.size() << std::endl;
260 <<
"more than half of candidate vertices lost " <<
pvs.size() <<
' ' <<
clusters.size();
262 if (
pvs.empty() && seltks.size() > 5)
264 <<
"no vertex found with " << seltks.size() <<
" tracks and " <<
clusters.size() <<
" vertex-candidates";
267 if (
pvs.size() > 1) {
272 for (std::vector<TransientVertex>::const_iterator
iv =
pvs.begin();
iv !=
pvs.end();
iv++) {
279 if ((bse.cxx() <= 0.) || (bse.cyy() <= 0.) || (bse.czz() <= 0.)) {
286 std::cout <<
"RecoVertex/PrimaryVertexProducer: " 287 <<
"Beamspot with invalid errors " << bse.matrix() << std::endl;
288 std::cout <<
"Will put Vertex derived from dummy-fake BeamSpot into Event.\n";
293 std::cout <<
"RecoVertex/PrimaryVertexProducer: " 294 <<
" will put Vertex derived from BeamSpot into Event.\n";
301 for (reco::VertexCollection::const_iterator
v = vColl.begin();
v != vColl.end(); ++
v) {
302 std::cout <<
"recvtx " << ivtx++ <<
"#trk " << std::setw(3) <<
v->tracksSize() <<
" chi2 " << std::setw(4)
303 <<
v->chi2() <<
" ndof " << std::setw(3) <<
v->ndof() <<
" x " << std::setw(6) <<
v->position().x()
304 <<
" dx " << std::setw(6) <<
v->xError() <<
" y " << std::setw(6) <<
v->position().y() <<
" dy " 305 << std::setw(6) <<
v->yError() <<
" z " << std::setw(6) <<
v->position().z() <<
" dz " << std::setw(6)
308 std::cout <<
" t " << std::setw(6) <<
v->t() <<
" dt " << std::setw(6) <<
v->tError();
edm::EDGetTokenT< reco::BeamSpot > bsToken
const edm::ESGetToken< TransientTrackBuilder, TransientTrackRecord > theTTBToken
std::vector< algo > algorithms
edm::EDGetTokenT< reco::VertexCollection > recoveryVtxToken
T const * product() const
edm::EDGetTokenT< reco::TrackCollection > trkToken
TrackClusterizerInZ * theTrackClusterizer
std::vector< Vertex > VertexCollection
collection of Vertex objects
Log< level::Error, false > LogError
edm::EDGetTokenT< edm::ValueMap< float > > trkTimesToken
edm::EDGetTokenT< edm::ValueMap< float > > trkTimeResosToken
virtual std::vector< std::vector< reco::TransientTrack > > clusterize(const std::vector< reco::TransientTrack > &tracks) const =0
TrackFilterForPVFindingBase * theTrackFilter
bool getData(T &iHolder) const
virtual std::vector< reco::TransientTrack > select(const std::vector< reco::TransientTrack > &tracks) const =0
ROOT::Math::SMatrix< double, 3, 3, ROOT::Math::MatRepSym< double, 3 > > AlgebraicSymMatrix33