@@ -405,10 +405,6 @@ namespace mkfit {
405405 fread (&beamSpot_, sizeof (BeamSpot), 1 , fp);
406406 }
407407
408- if (Config::kludgeCmsHitErrors) {
409- kludge_cms_hit_errors ();
410- }
411-
412408 if (!Config::silent)
413409 printf (" Read complete, %d simtracks on file.\n " , nt);
414410 }
@@ -472,41 +468,6 @@ namespace mkfit {
472468
473469 // ------------------------------------------------------------------------------
474470
475- void Event::kludge_cms_hit_errors () {
476- // Enforce Vxy on all layers, Vz on pixb only.
477-
478- const float Exy = 15 * 1e-4 , Vxy = Exy * Exy;
479- const float Ez = 30 * 1e-4 , Vz = Ez * Ez;
480-
481- int nl = layerHits_.size ();
482-
483- int cnt = 0 ;
484-
485- for (int il = 0 ; il < nl; il++) {
486- if (layerHits_[il].empty ())
487- continue ;
488-
489- for (Hit &h : layerHits_[il]) {
490- SVector6 &c = h.error_nc ();
491-
492- float vxy = c[0 ] + c[2 ];
493- if (vxy < Vxy) {
494- c[0 ] *= Vxy / vxy;
495- c[2 ] *= Vxy / vxy;
496- ++cnt;
497- }
498- if (il < 4 && c[5 ] < Vz) {
499- c[5 ] = Vz;
500- ++cnt;
501- }
502- }
503- }
504-
505- printf (" Event::kludge_cms_hit_errors processed %d layers, kludged %d entries.\n " , nl, cnt);
506- }
507-
508- // ------------------------------------------------------------------------------
509-
510471 int Event::clean_cms_simtracks () {
511472 // Sim tracks from cmssw have the following issues:
512473 // - hits are not sorted by layer;
0 commit comments