CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
CkfDebugger.cc
Go to the documentation of this file.
1 #include "RecoTracker/DebugTools/interface/CkfDebugger.h"
2 #include "SimDataFormats/TrackingHit/interface/PSimHitContainer.h"
3 #include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
4 #include "TrackingTools/PatternTools/interface/TrajectoryMeasurement.h"
5 #include "TrackingTools/PatternTools/interface/Trajectory.h"
6 #include "Geometry/CommonTopologies/interface/Topology.h"
7 #include "DataFormats/GeometryCommonDetAlgo/interface/MeasurementVector.h"
8 #include "RecoTracker/MeasurementDet/interface/MeasurementTracker.h"
9 #include "MagneticField/Records/interface/IdealMagneticFieldRecord.h"
10 #include "RecoTracker/DebugTools/interface/TSOSFromSimHitFactory.h"
11 #include "TrackingTools/MeasurementDet/interface/MeasurementDet.h"
12 #include "TrackingTools/MaterialEffects/interface/PropagatorWithMaterial.h"
13 #include "TrackingTools/KalmanUpdators/interface/Chi2MeasurementEstimator.h"
14 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
15 
16 #include "RecoTracker/DebugTools/interface/GluedDetFromDetUnit.h"
17 #include "DataFormats/TrackerRecHit2D/interface/SiStripMatchedRecHit2D.h"
18 
19 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
20 #include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"
21 #include "RecoTracker/TkNavigation/interface/TkLayerLess.h"
22 #include "RecoTracker/Record/interface/TrackerRecoGeometryRecord.h"
23 #include "TrackingTools/DetLayers/interface/BarrelDetLayer.h"
24 // #include "RecoTracker/TkDetLayers/interface/PixelBarrelLayer.h"
25 
26 #include "TrackingTools/TransientTrackingRecHit/interface/TrackingRecHitProjector.h"
27 #include "RecoTracker/TransientTrackingRecHit/interface/ProjectedRecHit2D.h"
28 
29 #include <iostream>
30 #include <iomanip>
31 #include <sstream>
32 
33 using namespace std;
34 
35 CkfDebugger::CkfDebugger( edm::EventSetup const & es ):totSeeds(0)
36 {
37  file = new TFile("out.root","recreate");
38  hchi2seedAll = new TH1F("hchi2seedAll","hchi2seedAll",2000,0,200);
39  hchi2seedProb = new TH1F("hchi2seedProb","hchi2seedProb",2000,0,200);
40 
41  edm::ESHandle<TrackerGeometry> tracker;
42  es.get<TrackerDigiGeometryRecord>().get(tracker);
43  theTrackerGeom = &(*tracker);
44 
45  edm::ESHandle<MagneticField> theField;
46  es.get<IdealMagneticFieldRecord>().get(theField);
47  theMagField = &(*theField);
48 
49  //Retrieve tracker topology from geometry
50  edm::ESHandle<TrackerTopology> tTopoHand;
51  es.get<IdealGeometryRecord>().get(tTopoHand);
52  theTopo=tTopoHand.product();
53 
54 
55  for (int i=0; i!=17; i++){
56  dump.push_back(0);
57  }
58 
59  std::stringstream title;
60  for (int i=0; i!=6; i++)
61  for (int j=0; j!=9; j++){
62  if (i==0 && j>2) break;
63  if (i==1 && j>1) break;
64  if (i==2 && j>3) break;
65  if (i==3 && j>2) break;
66  if (i==4 && j>5) break;
67  if (i==5 && j>8) break;
68  dump2[pair<int,int>(i,j)]=0;
69  dump3[pair<int,int>(i,j)]=0;
70  dump4[pair<int,int>(i,j)]=0;
71  dump5[pair<int,int>(i,j)]=0;
72  dump6[pair<int,int>(i,j)]=0;
73  title.str("");
74  title << "pullX_" << i+1 << "-" << j+1 << "_sh-rh";
75  hPullX_shrh[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
76  title.str("");
77  title << "pullY_" << i+1 << "-" << j+1 << "_sh-rh";
78  hPullY_shrh[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
79  title.str("");
80  title << "pullX_" << i+1 << "-" << j+1 << "_sh-st";
81  hPullX_shst[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
82  title.str("");
83  title << "pullY_" << i+1 << "-" << j+1 << "_sh-st";
84  hPullY_shst[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
85  title.str("");
86  title << "pullX_" << i+1 << "-" << j+1 << "_st-rh";
87  hPullX_strh[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
88  title.str("");
89  title << "pullY_" << i+1 << "-" << j+1 << "_st-rh";
90  hPullY_strh[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
91  title.str("");
92  title << "PullGP_X_" << i+1 << "-" << j+1 << "_sh-st";
93  hPullGP_X_shst[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
94  title.str("");
95  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_sh-st";
96  hPullGP_Y_shst[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
97  title.str("");
98  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_sh-st";
99  hPullGP_Z_shst[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
100  if ( ((i==2||i==4)&&(j==0||j==1)) || (i==3||i==5) ){
101  title.str("");
102  title << "pullM_" << i+1 << "-" << j+1 << "_sh-rh";
103  hPullM_shrh[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
104  title.str("");
105  title << "pullS_" << i+1 << "-" << j+1 << "_sh-rh";
106  hPullS_shrh[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
107  title.str("");
108  title << "pullM_" << i+1 << "-" << j+1 << "_sh-st";
109  hPullM_shst[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
110  title.str("");
111  title << "pullS_" << i+1 << "-" << j+1 << "_sh-st";
112  hPullS_shst[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
113  title.str("");
114  title << "pullM_" << i+1 << "-" << j+1 << "_st-rh";
115  hPullM_strh[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
116  title.str("");
117  title << "pullS_" << i+1 << "-" << j+1 << "_st-rh";
118  hPullS_strh[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
119  }
120  }
121 
122  hPullGPXvsGPX_shst = new TH2F("PullGPXvsGPX_shst","PullGPXvsGPX_shst",1000,-50,50,100,-50,50);
123  hPullGPXvsGPY_shst = new TH2F("PullGPXvsGPY_shst","PullGPXvsGPY_shst",1000,-50,50,100,-50,50);
124  hPullGPXvsGPZ_shst = new TH2F("PullGPXvsGPZ_shst","PullGPXvsGPZ_shst",1000,-50,50,200,-100,100);
125  hPullGPXvsGPr_shst = new TH2F("PullGPXvsGPr_shst","PullGPXvsGPr_shst",1000,-50,50,300,-150,150);
126  hPullGPXvsGPeta_shst = new TH2F("PullGPXvsGPeta_shst","PullGPXvsGPeta_shst",1000,-50,50,50,-2.5,2.5);
127  hPullGPXvsGPphi_shst = new TH2F("PullGPXvsGPphi_shst","PullGPXvsGPphi_shst",1000,-50,50,63,0,6.3);
128 
129  seedWithDelta=0;
130  problems=0;
131  no_sim_hit=0;
132  no_layer=0;
133  layer_not_found=0;
134  det_not_found=0;
135  chi2gt30=0;
136  chi2gt30delta=0;
137  chi2gt30deltaSeed=0;
138  chi2ls30=0;
139  simple_hit_not_found=0;
140  no_component=0;
141  only_one_component=0;
142  matched_not_found=0;
143  matched_not_associated=0;
144  partner_det_not_fuond=0;
145  glued_det_not_fuond=0;
146  propagation=0;
147  other=0;
148  totchi2gt30=0;
149 }
150 
151 void CkfDebugger::printSimHits( const edm::Event& iEvent)
152 {
153  edm::LogVerbatim("CkfDebugger") << "\nEVENT #" << iEvent.id();
154 
155  hitAssociator = new TrackerHitAssociator(iEvent);//delete deleteHitAssociator() in TrackCandMaker.cc
156 
157  std::map<unsigned int, std::vector<PSimHit> >& theHitsMap = hitAssociator->SimHitMap;
158  idHitsMap.clear();
159 
160  for (std::map<unsigned int, std::vector<PSimHit> >::iterator it=theHitsMap.begin();
161  it!=theHitsMap.end();it++){
162  for (std::vector<PSimHit>::iterator isim = it->second.begin();
163  isim != it->second.end(); ++isim){
164  idHitsMap[isim->trackId()].push_back(&*isim);
165  }
166  }
167 
168  for (std::map<unsigned int,std::vector<PSimHit*> >::iterator it=idHitsMap.begin();
169  it!=idHitsMap.end();it++){
170  sort(it->second.begin(),it->second.end(),less_mag());
171  for (std::vector<PSimHit*>::iterator isim = it->second.begin();
172  isim != it->second.end(); ++isim){
173  const GeomDetUnit* detUnit = theTrackerGeom->idToDetUnit( DetId((*isim)->detUnitId()));
174  dumpSimHit(SimHit( (*isim), detUnit));
175  }
176  }
177 
178 }
179 
180 void CkfDebugger::dumpSimHit(const SimHit& hit) const
181 {
182  GlobalPoint pos = hit.globalPosition();
183  edm::LogVerbatim("CkfDebugger") << "SimHit pos" << pos
184  << " r=" << pos.perp() << " phi=" << pos.phi()
185  << " trackId=" << hit.trackId()
186  << " particleType=" << hit.particleType()
187  << " pabs=" << hit.pabs()
188  << " processType=" << hit. processType();
189 }
190 
191 
192 bool CkfDebugger::analyseCompatibleMeasurements(const Trajectory& traj,
193  const std::vector<TrajectoryMeasurement>& meas,
194  const MeasurementTracker* aMeasurementTracker,
195  const Propagator* propagator,
196  const Chi2MeasurementEstimatorBase* estimator,
197  const TransientTrackingRecHitBuilder* aTTRHBuilder)
198 {
199  LogTrace("CkfDebugger") << "\nnow in analyseCompatibleMeasurements" ;
200  LogTrace("CkfDebugger") << "number of input hits:" << meas.size() ;
201  for(std::vector<TrajectoryMeasurement>::const_iterator tmpIt=meas.begin();tmpIt!=meas.end();tmpIt++){
202  if (tmpIt->recHit()->isValid()) LogTrace("CkfDebugger") << "valid hit at position:" << tmpIt->recHit()->globalPosition() ;
203  }
204  theForwardPropagator = propagator;
205  theChi2 = estimator;
206  theMeasurementTracker = aMeasurementTracker;
207  theGeomSearchTracker = theMeasurementTracker->geometricSearchTracker();
208  theTTRHBuilder = aTTRHBuilder;
209  unsigned int trajId = 0;
210  if ( !correctTrajectory(traj, trajId)) {
211  LogTrace("CkfDebugger") << "trajectory not correct" ;
212  return true;
213  } // only correct trajectories analysed
214  LogTrace("CkfDebugger") << "correct trajectory" ;
215 
216  if (traj.measurements().size() == 2){
217  if ( testSeed(traj.firstMeasurement().recHit(),traj.lastMeasurement().recHit(),traj.lastMeasurement().updatedState()) == -1 ) {
218  LogTrace("CkfDebugger") << "Seed has delta" ;
219  seedWithDelta++;
220  return false;//true;//false?
221  }
222  }
223 
224  //const PSimHit* correctHit = nextCorrectHit(traj, trajId);
225  //if ( correctHit == 0) return true; // no more simhits on this track
226  std::vector<const PSimHit*> correctHits = nextCorrectHits(traj, trajId);
227  if ( correctHits.size() == 0) return true; // no more simhits on this track
228 
229  for (std::vector<const PSimHit*>::iterator corHit=correctHits.begin();corHit!=correctHits.end();corHit++){
230  for (std::vector<TM>::const_iterator i=meas.begin(); i!=meas.end(); i++) {
231  if (correctMeas( *i, *corHit)) {
232  LogTrace("CkfDebugger") << "Correct hit found at position " << i-meas.begin() ;
233  return true;
234  }
235  }
236  }
237 
238  //debug why the first hit in correctHits is not found
239  //FIXME should loop over all hits
240  const PSimHit* correctHit = *(correctHits.begin());
241 
242  // correct hit not found
243  edm::LogVerbatim("CkfDebugger") << std::endl << "CkfDebugger: problem found: correct hit not found by findCompatibleMeasurements" ;
244  edm::LogVerbatim("CkfDebugger") << "The correct hit position is " << position(correctHit) << " lp " << correctHit->localPosition() ;
245  edm::LogVerbatim("CkfDebugger") << "The size of the meas vector is " << meas.size() ;
246  dump[0]++;problems++;
247 
248  for (std::vector<TM>::const_iterator i=meas.begin(); i!=meas.end(); i++) {
249  edm::LogVerbatim("CkfDebugger") << "Is the hit valid? " << i->recHit()->isValid() ;
250  if (i->recHit()->isValid()) {
251  edm::LogVerbatim("CkfDebugger") << "RecHit at " << i->recHit()->globalPosition()
252  << " layer " << ((i->recHit()->det()->geographicalId().rawId() >>16) & 0xF)
253  << " subdet " << i->recHit()->det()->geographicalId().subdetId()
254  << " Chi2 " << i->estimate() ;
255  }
256  else if (i->recHit()->det() == 0) {
257  edm::LogVerbatim("CkfDebugger") << "Invalid RecHit returned with zero Det pointer" ;
258  }
259  else if (i->recHit()->det() == det(correctHit)) {
260  edm::LogVerbatim("CkfDebugger") << "Invalid hit returned in correct Det" ;
261  }
262  else {
263  edm::LogVerbatim("CkfDebugger") << "Invalid hit returned in Det at gpos " << i->recHit()->det()->position()
264  << " correct Det is at " << det(correctHit)->position() ;
265  }
266  }
267 
268  //Look if the correct RecHit exists
269  std::pair<CTTRHp, double> correctRecHit =
270  analyseRecHitExistance( *correctHit, traj.lastMeasurement().updatedState());
271  if (correctRecHit.first==0 ) {
272  //the hit does not exist or is uncorrectly matched
273  if ( fabs(correctRecHit.second-0)<0.01 ) {dump[1]++;}//other
274  if ( fabs(correctRecHit.second+1)<0.01 ) {dump[8]++;}//propagation
275  if ( fabs(correctRecHit.second+2)<0.01 ) {dump[9]++;}//No component is found
276  if ( fabs(correctRecHit.second+3)<0.01 ) {dump[10]++;}//Partner measurementDet not found
277  if ( fabs(correctRecHit.second+4)<0.01 ) {dump[11]++;}//glued MeasurementDet not found
278  if ( fabs(correctRecHit.second+5)<0.01 ) {dump[12]++;}//matched not found
279  if ( fabs(correctRecHit.second+6)<0.01 ) {dump[13]++;}//Matched not associated
280  if ( fabs(correctRecHit.second+7)<0.01 ) {dump[14]++;}//Only one component is found
281  if ( fabs(correctRecHit.second+8)<0.01 ) {dump[15]++;}//not found (is not a glued det)
282  }
283  else {
284  //the hit exists: why wasn't it found?
285  int result = analyseRecHitNotFound(traj,correctRecHit.first);
286  if (result == 5){
287  if (correctRecHit.second>30) {
288  edm::LogVerbatim("CkfDebugger") << "Outling RecHit at pos=" << correctRecHit.first->globalPosition()
289  << " from SimHit at pos="<< position(correctHit)
290  << " det=" << correctHit->detUnitId() << " process=" << correctHit->processType() ;
291  if (hasDelta(correctHit)){
292  edm::LogVerbatim("CkfDebugger") << "there are deltas on this det" ;
293  chi2gt30delta++;
294  dump5[pair<int,int>((correctRecHit.first->det()->geographicalId().subdetId()-1),(layer(correctRecHit.first->det()))-1)]++;
295  }else{
296  edm::LogVerbatim("CkfDebugger") << "no deltas on this det" ;
297  dump[5]++;
298  chi2gt30++;
299  dump3[pair<int,int>((correctRecHit.first->det()->geographicalId().subdetId()-1),(layer(correctRecHit.first->det()))-1)]++;
300  CTTRHp h1 = traj.measurements()[0].recHit();
301  CTTRHp h2 = traj.measurements()[1].recHit();
302  TSOS t = traj.measurements()[1].updatedState();
303  double chi2 = testSeed(h1,h2,t);
304  if (chi2==-1) {
305  edm::LogVerbatim("CkfDebugger") << "there were deltas in the seed" ;
306  chi2gt30deltaSeed++;
307  }
308  else {
309  hchi2seedProb->Fill(chi2);
310  edm::LogVerbatim("CkfDebugger") << "no deltas in the seed. What is wrong?" ;
311 
312  TSOS detState = theForwardPropagator->propagate( traj.lastMeasurement().updatedState(), correctRecHit.first->det()->surface());
313  TSOS simDetState = theForwardPropagator->propagate( traj.lastMeasurement().updatedState(), det(correctHit)->surface());
314 
315  if (true/*detState.globalMomentum().y()>0*/){
316  int subdetId = correctRecHit.first->det()->geographicalId().subdetId();
317  int layerId = layer(correctRecHit.first->det());
318 
319 
320  LogTrace("CkfDebugger") << "position(correctHit)=" << position(correctHit) ;
321  LogTrace("CkfDebugger") << "correctRecHit.first->globalPosition()=" << correctRecHit.first->globalPosition() ;
322  LogTrace("CkfDebugger") << "detState.globalPosition()=" << detState.globalPosition() ;
323  LogTrace("CkfDebugger") << "simDetState.globalPosition()=" << simDetState.globalPosition() ;
324 
325  LogTrace("CkfDebugger") << "correctHit->localPosition()=" << correctHit->localPosition() ;
326  LogTrace("CkfDebugger") << "correctRecHit.first->localPosition()=" << correctRecHit.first->localPosition() ;
327  LogTrace("CkfDebugger") << "correctRecHit.first->localPositionError()=" << correctRecHit.first->localPositionError() ;
328  LogTrace("CkfDebugger") << "detState.localPosition()=" << detState.localPosition() ;
329  LogTrace("CkfDebugger") << "detState.localError().positionError()=" << detState.localError().positionError() ;
330  LogTrace("CkfDebugger") << "simDetState.localPosition()=" << simDetState.localPosition() ;
331  LogTrace("CkfDebugger") << "simDetState.localError().positionError()=" << simDetState.localError().positionError() ;
332  double pullx_shrh = (correctHit->localPosition().x()-correctRecHit.first->localPosition().x())/
333  sqrt(correctRecHit.first->localPositionError().xx());
334  double pully_shrh = 0;
335  if (correctRecHit.first->localPositionError().yy()!=0)
336  pully_shrh = (correctHit->localPosition().y()-correctRecHit.first->localPosition().y())/
337  sqrt(correctRecHit.first->localPositionError().yy());
338  double pullx_shst = (correctHit->localPosition().x()-simDetState.localPosition().x())/
339  sqrt(simDetState.localError().positionError().xx());
340  double pully_shst = (correctHit->localPosition().y()-simDetState.localPosition().y())/
341  sqrt(simDetState.localError().positionError().yy());
342 
343  LogTrace("CkfDebugger") << "pullx(sh-rh)=" << pullx_shrh ;
344  LogTrace("CkfDebugger") << "pully(sh-rh)=" << pully_shrh ;
345  LogTrace("CkfDebugger") << "pullx(sh-st)=" << pullx_shst ;
346  LogTrace("CkfDebugger") << "pully(sh-st)=" << pully_shst ;
347 
348  LogTrace("CkfDebugger") << "pullx(st-rh)=" << (detState.localPosition().x()-correctRecHit.first->localPosition().x())/
349  sqrt(correctRecHit.first->localPositionError().xx()+detState.localError().positionError().xx()) ;
350 
351  std::pair<double,double> pulls = computePulls(correctRecHit.first, detState);
352  if (subdetId>0 &&subdetId<7 && layerId>0 && layerId<10) {
353  stringstream title;
354  title.str("");
355  title << "pullX_" << subdetId << "-" << layerId << "_sh-rh";
356  hPullX_shrh[title.str()]->Fill( pullx_shrh );
357  title.str("");
358  title << "pullY_" << subdetId << "-" << layerId << "_sh-rh";
359  hPullY_shrh[title.str()]->Fill( pully_shrh );
360  title.str("");
361  title << "pullX_" << subdetId << "-" << layerId <<"_sh-st";
362  hPullX_shst[title.str()]->Fill( pullx_shst );
363  title.str("");
364  title << "pullY_" << subdetId << "-" << layerId <<"_sh-st";
365  hPullY_shst[title.str()]->Fill( pully_shst );
366  title.str("");
367  title << "pullX_" << subdetId << "-" << layerId <<"_st-rh";
368  hPullX_strh[title.str()]->Fill(pulls.first);
369  title.str("");
370  title << "pullY_" << subdetId << "-" << layerId <<"_st-rh";
371  hPullY_strh[title.str()]->Fill(pulls.second);
372 
373  GlobalPoint shGPos = position(correctHit);
374  GlobalPoint stGPos = simDetState.globalPosition();
375  GlobalError stGPosErr = simDetState.cartesianError().position();
376  double pullGPx = (shGPos.x()-stGPos.x())/sqrt(stGPosErr.cxx());
377  title.str("");
378  title << "PullGP_X_" << subdetId << "-" << layerId << "_sh-st";
379  hPullGP_X_shst[title.str()]->Fill(pullGPx);
380  title.str("");
381  title << "PullGP_Y_" << subdetId << "-" << layerId << "_sh-st";
382  hPullGP_Y_shst[title.str()]->Fill((shGPos.y()-stGPos.y())/sqrt(stGPosErr.cyy()));
383  title.str("");
384  title << "PullGP_Z_" << subdetId << "-" << layerId << "_sh-st";
385  hPullGP_Z_shst[title.str()]->Fill((shGPos.z()-stGPos.z())/sqrt(stGPosErr.czz()));
386 
387  if (subdetId==3&&layerId==1){
388  hPullGPXvsGPX_shst->Fill(pullGPx,shGPos.x());
389  hPullGPXvsGPY_shst->Fill(pullGPx,shGPos.y());
390  hPullGPXvsGPZ_shst->Fill(pullGPx,shGPos.z());
391  hPullGPXvsGPr_shst->Fill(pullGPx,shGPos.mag());
392  hPullGPXvsGPeta_shst->Fill(pullGPx,shGPos.eta());
393  hPullGPXvsGPphi_shst->Fill(pullGPx,shGPos.phi());
394  }
395  if (dynamic_cast<const SiStripMatchedRecHit2D*>(correctRecHit.first->hit())) {
396  LogTrace("CkfDebugger") << "MONO HIT";
397  auto m = dynamic_cast<const SiStripMatchedRecHit2D*>(correctRecHit.first->hit())->monoHit();
398  CTTRHp tMonoHit = theTTRHBuilder->build(&m);
399  const PSimHit sMonoHit = *(hitAssociator->associateHit(*tMonoHit->hit()).begin());
400  TSOS monoState = theForwardPropagator->propagate( traj.lastMeasurement().updatedState(), tMonoHit->det()->surface());
401  double pullM_shrh = (sMonoHit.localPosition().x()-tMonoHit->localPosition().x())/
402  sqrt(tMonoHit->localPositionError().xx());
403  double pullM_shst = (sMonoHit.localPosition().x()-monoState.localPosition().x())/
404  sqrt(monoState.localError().positionError().xx());
405  std::pair<double,double> pullsMono = computePulls(tMonoHit, monoState);
406  title.str("");
407  title << "pullM_" << subdetId << "-" << layerId << "_sh-rh";
408  hPullM_shrh[title.str()]->Fill(pullM_shrh);
409  title.str("");
410  title << "pullM_" << subdetId << "-" << layerId << "_sh-st";
411  hPullM_shst[title.str()]->Fill(pullM_shst);
412  title.str("");
413  title << "pullM_" << subdetId << "-" << layerId << "_st-rh";
414  hPullM_strh[title.str()]->Fill(pullsMono.first);
415 
416  LogTrace("CkfDebugger") << "STEREO HIT";
417  auto s= dynamic_cast<const SiStripMatchedRecHit2D*>(correctRecHit.first->hit())->stereoHit();
418  CTTRHp tStereoHit = theTTRHBuilder->build(&s);
419  const PSimHit sStereoHit = *(hitAssociator->associateHit(*tStereoHit->hit()).begin());
420  TSOS stereoState = theForwardPropagator->propagate( traj.lastMeasurement().updatedState(), tStereoHit->det()->surface());
421  double pullS_shrh = (sStereoHit.localPosition().x()-tStereoHit->localPosition().x())/
422  sqrt(tStereoHit->localPositionError().xx());
423  double pullS_shst = (sStereoHit.localPosition().x()-stereoState.localPosition().x())/
424  sqrt(stereoState.localError().positionError().xx());
425  std::pair<double,double> pullsStereo = computePulls(tStereoHit, stereoState);
426  title.str("");
427  title << "pullS_" << subdetId << "-" << layerId << "_sh-rh";
428  hPullS_shrh[title.str()]->Fill(pullS_shrh);
429  title.str("");
430  title << "pullS_" << subdetId << "-" << layerId << "_sh-st";
431  hPullS_shst[title.str()]->Fill(pullS_shst);
432  title.str("");
433  title << "pullS_" << subdetId << "-" << layerId << "_st-rh";
434  hPullS_strh[title.str()]->Fill(pullsStereo.first);
435  }
436  } else
437  edm::LogVerbatim("CkfDebugger") << "unexpected result: wrong det or layer id "
438  << subdetId << " " << layerId << " "
439  << correctRecHit.first->det()->geographicalId().rawId();
440  }
441  }
442  }
443  }
444  else {
445  edm::LogVerbatim("CkfDebugger") << "unexpected result " << correctRecHit.second ;
446  dump[6]++;chi2ls30++;
447  }
448  }
449  else dump[result]++;
450  if (result == 3){
451  dump2[pair<int,int>((correctRecHit.first->det()->geographicalId().subdetId()-1),(layer(correctRecHit.first->det()))-1)]++;
452  }
453  if (result == 4){
454  dump4[pair<int,int>((correctRecHit.first->det()->geographicalId().subdetId()-1),(layer(correctRecHit.first->det()))-1)]++;
455  }
456  if (correctRecHit.second>30) {
457  dump[7]++;totchi2gt30++;
458  }
459  }
460  return false;
461 }
462 
463 
464 bool CkfDebugger::correctTrajectory( const Trajectory& traj,unsigned int& trajId) const
465 {
466  LogTrace("CkfDebugger") << "now in correctTrajectory" ;
467  Trajectory::RecHitContainer hits = traj.recHits();
468 
469  std::vector<SimHitIdpr> currentTrackId = hitAssociator->associateHitId(*hits.front()->hit());
470  if (currentTrackId.size() == 0) return false;
471 
472  for (Trajectory::RecHitContainer::const_iterator rh=hits.begin(); rh!=hits.end(); ++rh) {
473 
474  //if invalid hit exit
475  if (!(*rh)->hit()->isValid()) {
476  //LogTrace("CkfDebugger") << "invalid hit" ;
477  return false;
478  }
479 
480  //if hits from deltas exit
481  bool nogoodhit = true;
482  std::vector<PSimHit> assSimHits = hitAssociator->associateHit(*(*rh)->hit());
483  for (std::vector<PSimHit>::iterator shit=assSimHits.begin();shit!=assSimHits.end();shit++){
484  if (goodSimHit(*shit)) nogoodhit=false;
485  }
486  if (nogoodhit) return false;
487 
488  //all hits must be associated to the same sim track
489  bool test = true;
490  std::vector<SimHitIdpr> nextTrackId = hitAssociator->associateHitId(*(*rh)->hit());
491  for (std::vector<SimHitIdpr>::iterator i=currentTrackId.begin();i!=currentTrackId.end();i++){
492  for (std::vector<SimHitIdpr>::iterator j=nextTrackId.begin();j!=nextTrackId.end();j++){
493  if (i->first == j->first) test = false;
494  //LogTrace("CkfDebugger") << "valid " << *i << " " << *j ;
495  trajId = j->first;
496  }
497  }
498  if (test) {/*LogTrace("CkfDebugger") << "returning false" ;*/return false;}
499  // std::vector<PSimHit*> simTrackHits = idHitsMap[trajId];
500  // if (!goodSimHit(simTrackHits.))
501  }
502  //LogTrace("CkfDebugger") << "returning true" ;
503  return true;
504 }
505 
506 int CkfDebugger::assocTrackId(CTTRHp rechit) const
507 {
508  LogTrace("CkfDebugger") << "now in assocTrackId" ;
509 
510  if (!rechit->hit()->isValid()) {
511  return -1;
512  }
513 
514  std::vector<SimHitIdpr> ids = hitAssociator->associateHitId(*rechit->hit());
515  if (ids.size()!=0) {
516  return ids[0].first;//FIXME if size>1!!
517  }
518  else {
519  return -1;
520  }
521 }
522 
523 
524 vector<const PSimHit*> CkfDebugger::nextCorrectHits( const Trajectory& traj, unsigned int& trajId)
525 {
526  std::vector<const PSimHit*> result;
527  // find the component of the RecHit at largest distance from origin (FIXME: should depend on propagation direction)
528  LogTrace("CkfDebugger") << "now in nextCorrectHits" ;
529  TransientTrackingRecHit::ConstRecHitPointer lastRecHit = traj.lastMeasurement().recHit();
530  TransientTrackingRecHit::RecHitContainer comp = lastRecHit->transientHits();
531  if (!comp.empty()) {
532  float maxR = 0;
533  for (TransientTrackingRecHit::RecHitContainer::const_iterator ch=comp.begin();
534  ch!=comp.end(); ++ch) {
535  if ((*ch)->globalPosition().mag() > maxR) lastRecHit = *ch;
536  maxR = (*ch)->globalPosition().mag();
537  }
538  }
539  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: lastRecHit is at gpos " << lastRecHit->globalPosition()
540  << " layer " << layer((lastRecHit->det()))
541  << " subdet " << lastRecHit->det()->geographicalId().subdetId() ;
542 
543  //find the simHits associated to the recHit
544  const std::vector<PSimHit>& pSimHitVec = hitAssociator->associateHit(*lastRecHit->hit());
545  for (std::vector<PSimHit>::const_iterator shit=pSimHitVec.begin();shit!=pSimHitVec.end();shit++){
546  const GeomDetUnit* detUnit = theTrackerGeom->idToDetUnit( DetId(shit->detUnitId()));
547  LogTrace("CkfDebugger") << "from hitAssociator SimHits are at GP=" << detUnit->toGlobal( shit->localPosition())
548  << " traId=" << shit->trackId() << " particleType " << shit->particleType()
549  << " pabs=" << shit->pabs() << " detUnitId=" << shit->detUnitId() << " layer " << layer((det(&*shit)))
550  << " subdet " << det(&*shit)->geographicalId().subdetId() ;
551  }
552 
553  //choose the simHit from the same track that has the highest tof
554  const PSimHit * lastPSH = 0;
555  if (!pSimHitVec.empty()) {
556  float maxTOF = 0;
557  for (std::vector<PSimHit>::const_iterator ch=pSimHitVec.begin(); ch!=pSimHitVec.end(); ++ch) {
558  if ( ( ch->trackId()== trajId) && (ch->timeOfFlight() > maxTOF) && ( goodSimHit(*ch) )) {
559  lastPSH = &*ch;
560  maxTOF = lastPSH->timeOfFlight();
561  }
562  }
563  }
564  else return result;//return empty vector: no more hits on the sim track
565  if (lastPSH == 0) return result; //return empty vector: no more good hits on the sim track
566  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: corresponding SimHit is at gpos " << position(&*lastPSH) ;
567 
568  //take the simHits on the simTrack that are in the nextLayer (could be > 1 if overlap or matched)
569  std::vector<PSimHit*> trackHits = idHitsMap[trajId];
570  if (fabs((double)(trackHits.back()->detUnitId()-lastPSH->detUnitId()))<1 ) return result;//end of sim track
571  std::vector<PSimHit*>::iterator currentIt = trackHits.end();
572  for (std::vector<PSimHit*>::iterator it=trackHits.begin();
573  it!=trackHits.end();it++){
574  if (goodSimHit(**it) && //good hit
575  ( lastPSH->timeOfFlight()<(*it)->timeOfFlight() ) && //greater tof
576  //( fabs((double)((*it)->detUnitId()-(lastPSH->detUnitId()) ))>1) && //not components of the same matched hit
577  ( (det(lastPSH)->geographicalId().subdetId()!=det(*it)->geographicalId().subdetId()) ||
578  (layer(det(lastPSH))!=layer(det(*it)) ) ) //change layer or detector(tib,tob,...)
579  ){
580  edm::LogVerbatim("CkfDebugger") << "Next good PSimHit is at gpos " << position(*it) ;
581  result.push_back(*it);
582  currentIt = it;
583  break;
584  }
585  }
586  bool samelayer = true;
587  if (currentIt!=(trackHits.end()-1) && currentIt!=trackHits.end()) {
588  for (std::vector<PSimHit*>::iterator nextIt = currentIt; (samelayer && nextIt!=trackHits.end()) ;nextIt++){
589  if (goodSimHit(**nextIt)){
590  if ( (det(*nextIt)->geographicalId().subdetId()==det(*currentIt)->geographicalId().subdetId()) &&
591  (layer(det(*nextIt))==layer(det(*currentIt)) ) ) {
592  result.push_back(*nextIt);
593  }
594  else samelayer = false;
595  }
596  }
597  }
598 
599  return result;
600 }
601 
602 bool CkfDebugger::goodSimHit(const PSimHit& sh) const
603 {
604  if (sh.pabs() > 0.9) return true; // GeV, reject delta rays from association
605  else return false;
606 }
607 
608 
609 bool CkfDebugger::associated(CTTRHp rechit, const PSimHit& pSimHit) const
610 {
611  LogTrace("CkfDebugger") << "now in associated" ;
612 
613  if (!rechit->isValid()) return false;
614  // LogTrace("CkfDebugger") << "rec hit valid" ;
615  const std::vector<PSimHit>& pSimHitVec = hitAssociator->associateHit(*rechit->hit());
616  // LogTrace("CkfDebugger") << "size=" << pSimHitVec.size() ;
617  for (std::vector<PSimHit>::const_iterator shit=pSimHitVec.begin();shit!=pSimHitVec.end();shit++){
618  //const GeomDetUnit* detUnit = theTrackerGeom->idToDetUnit( DetId(shit->detUnitId()));
619  // LogTrace("CkfDebugger") << "pSimHit.timeOfFlight()=" << pSimHit.timeOfFlight()
620  // << " pSimHit.pabs()=" << pSimHit.pabs() << " GP=" << position(&pSimHit);
621  // LogTrace("CkfDebugger") << "(*shit).timeOfFlight()=" << (*shit).timeOfFlight()
622  // << " (*shit).pabs()=" << (*shit).pabs() << " GP=" << detUnit->toGlobal( shit->localPosition());
623  if ( ( fabs((*shit).timeOfFlight()-pSimHit.timeOfFlight())<1e-9 ) &&
624  ( fabs((*shit).pabs()-pSimHit.pabs())<1e-9 ) ) return true;
625  }
626  return false;
627 }
628 
629 bool CkfDebugger::correctMeas( const TM& tm, const PSimHit* correctHit) const
630 {
631  LogTrace("CkfDebugger") << "now in correctMeas" ;
632  CTTRHp recHit = tm.recHit();
633  if (recHit->isValid()) LogTrace("CkfDebugger") << "hit at position:" << recHit->globalPosition() ;
634  TransientTrackingRecHit::RecHitContainer comp = recHit->transientHits();
635  if (comp.empty()) {
636  // LogTrace("CkfDebugger") << "comp.empty()==true" ;
637  return associated( recHit, *correctHit);
638  }
639  else {
640  for (TransientTrackingRecHit::RecHitContainer::const_iterator ch=comp.begin();
641  ch!=comp.end(); ++ch) {
642  if ( associated( recHit, *correctHit)) {
643  // check if the other components are associated to the same trackId
644  for (TransientTrackingRecHit::RecHitContainer::const_iterator ch2=comp.begin();
645  ch2!=comp.end(); ++ch2) {
646  if (ch2 == ch) continue;
648  // LogTrace("CkfDebugger") << "correctHit->trackId()=" << correctHit->trackId() ;
649  bool test=true;
650  std::vector<SimHitIdpr> ids = hitAssociator->associateHitId(*(*ch2)->hit());
651  for (std::vector<SimHitIdpr>::iterator j=ids.begin();j!=ids.end();j++){
652  // LogTrace("CkfDebugger") << "j=" <<j->first;
653  if (correctHit->trackId()==j->first) {
654  test=false;
655  // LogTrace("CkfDebugger") << correctHit->trackId()<< " " <<j->first;
656  }
657  }
658  if (assocTrackId( *ch2) != ((int)( correctHit->trackId())) ) {LogTrace("CkfDebugger") << "returning false 1" ;/*return false;*/}//fixme
659  if (test) {
660  // LogTrace("CkfDebugger") << "returning false 2" ;
661  return false; // not all components from same simtrack
662  }
663  // if (assocTrackId( **ch2) != ((int)( correctHit->trackId())) ) {
664  // return false; // not all components from same simtrack
665  // }
666  }
667  return true; // if all components from same simtrack
668  }
669  }
670  return false;
671  }
672 }
673 
674 //this checks only if there is the rechit on the det where the sim hit is
675 pair<CTTRHp, double> CkfDebugger::analyseRecHitExistance( const PSimHit& sh, const TSOS& startingState)
676 {
677  LogTrace("CkfDebugger") << "now in analyseRecHitExistance" ;
678 
679  std::pair<CTTRHp, double> result;
680 
681  const MeasurementDet* simHitDet = theMeasurementTracker->idToDet( DetId( sh.detUnitId()));
682  TSOS simHitState = TSOSFromSimHitFactory()(sh, *det(&sh), *theMagField);
683  MeasurementDet::RecHitContainer recHits = simHitDet->recHits( simHitState);//take all hits from det
684 
685  //check if the hit is not present or is a problem of association
686  TSOS firstDetState = theForwardPropagator->propagate( startingState, det(&sh)->surface());
687  if (!firstDetState.isValid()) {
688  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: propagation failed from state " << startingState << " to first det surface "
689  << position(&sh) ;
690  propagation++;
691  return std::pair<CTTRHp, double>((CTTRHp)(0),-1);
692  }
693 
694  bool found = false;
695  for ( MeasurementDet::RecHitContainer::const_iterator rh = recHits.begin(); rh != recHits.end(); rh++) {
696  if ( associated( *rh, sh)) {
697  found = true;
698  result = std::pair<CTTRHp, double>(*rh,theChi2->estimate( firstDetState, **rh).second);
699  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: A RecHit associated to the correct Simhit exists at lpos "
700  << (**rh).localPosition()
701  << " gpos " << (**rh).globalPosition()
702  << " layer " << layer((**rh).det())
703  << " subdet " << (**rh).det()->geographicalId().subdetId()
704  << " Chi2 " << theChi2->estimate( firstDetState, **rh).second;
705  }
706  }
707  if (!found) {
708  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: there is no RecHit associated to the correct SimHit." ;
709  edm::LogVerbatim("CkfDebugger") << " There are " << recHits.size() << " RecHits in the simHit DetUnit" ;
710  edm::LogVerbatim("CkfDebugger") << "SH GP=" << position(&sh) << " subdet=" << det(&sh)->geographicalId().subdetId()
711  << " layer=" << layer(det(&sh)) ;
712  int y=0;
713  for (MeasurementDet::RecHitContainer::const_iterator rh = recHits.begin(); rh != recHits.end(); rh++)
714  edm::LogVerbatim("CkfDebugger") << "RH#" << y++ << " GP=" << (**rh).globalPosition() << " subdet=" << (**rh).det()->geographicalId().subdetId()
715  << " layer=" << layer((**rh).det()) ;
716  for ( MeasurementDet::RecHitContainer::const_iterator rh = recHits.begin(); rh != recHits.end(); rh++) {
717  edm::LogVerbatim("CkfDebugger") << "Non-associated RecHit at pos " << (**rh).localPosition() ;
718  }
719  }
720 
721  bool found2 = false;
722  const PSimHit* sh2;
723  StripSubdetector subdet( det(&sh)->geographicalId());
724  if (!subdet.glued()) {
725  edm::LogVerbatim("CkfDebugger") << "The DetUnit is not part of a GluedDet" ;
726  if (found) {
727  if (result.second>30){
728  LogTrace("CkfDebugger") << "rh->parameters()=" << result.first->parameters() ;
729  LogTrace("CkfDebugger") << "rh->parametersError()=" << result.first->parametersError() ;
730  MeasurementExtractor me(firstDetState);
731  AlgebraicVector r(result.first->parameters() - me.measuredParameters(*result.first));
732  LogTrace("CkfDebugger") << "me.measuredParameters(**rh)=" << me.measuredParameters(*result.first) ;
733  LogTrace("CkfDebugger") << "me.measuredError(**rh)=" << me.measuredError(*result.first) ;
734  AlgebraicSymMatrix R(result.first->parametersError() + me.measuredError(*result.first));
735  LogTrace("CkfDebugger") << "r=" << r ;
736  LogTrace("CkfDebugger") << "R=" << R ;
737  int ierr;
738  R.invert(ierr);
739  LogTrace("CkfDebugger") << "R(-1)=" << R ;
740  LogTrace("CkfDebugger") << "chi2=" << R.similarity(r) ;
741  }
742  return result;
743  }
744  else {
745  simple_hit_not_found++;
746  return std::pair<CTTRHp, double>((CTTRHp)(0),-8);//not found (is not a glued det)
747  }
748  } else {
749  edm::LogVerbatim("CkfDebugger") << "The DetUnit is part of a GluedDet" ;
750  DetId partnerDetId = DetId( subdet.partnerDetId());
751 
752  sh2 = pSimHit( sh.trackId(), partnerDetId);
753  if (sh2 == 0) {
754  edm::LogVerbatim("CkfDebugger") << "Partner DetUnit does not have a SimHit from the same track" ;
755  if (found) {
756  //projected rec hit
757  TrackingRecHitProjector<ProjectedRecHit2D> proj;
758  DetId gid = gluedId( subdet);
759  const MeasurementDet* gluedDet = theMeasurementTracker->idToDet( gid);
760  TSOS gluedTSOS = theForwardPropagator->propagate(startingState, gluedDet->geomDet().surface());
761  CTTRHp projHit = proj.project( *result.first,gluedDet->geomDet(),gluedTSOS).get();
762  //LogTrace("CkfDebugger") << proj.project( *result.first,gluedDet->geomDet(),gluedTSOS)->parameters() ;
763  //LogTrace("CkfDebugger") << projHit->parametersError() ;
764  double chi2 = theChi2->estimate(gluedTSOS, *proj.project( *result.first,gluedDet->geomDet(),gluedTSOS)).second;
765  return std::pair<CTTRHp, double>(projHit,chi2);
766  }
767  }
768  else {
769  edm::LogVerbatim("CkfDebugger") << "Partner DetUnit has a good SimHit at gpos " << position(sh2)
770  << " lpos " << sh2->localPosition() ;
771  //}
772 
773  const MeasurementDet* partnerDet = theMeasurementTracker->idToDet( partnerDetId);
774  if (partnerDet == 0) {
775  edm::LogVerbatim("CkfDebugger") << "Partner measurementDet not found!!!" ;
776  partner_det_not_fuond++;
777  return std::pair<CTTRHp, double>((CTTRHp)(0),-3);
778  }
779  TSOS simHitState2 = TSOSFromSimHitFactory()(*sh2, *det(sh2), *theMagField);
780  MeasurementDet::RecHitContainer recHits2 = partnerDet->recHits( simHitState2);
781 
782  TSOS secondDetState = theForwardPropagator->propagate( startingState, det(sh2)->surface());
783  if (!secondDetState.isValid()) {
784  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: propagation failed from state " << startingState << " to second det surface "
785  << position(sh2) ;
786  propagation++;
787  return std::pair<CTTRHp, double>((CTTRHp)(0),-1);
788  }
789 
790  for ( MeasurementDet::RecHitContainer::const_iterator rh = recHits2.begin(); rh != recHits2.end(); rh++) {
791  if ( associated( *rh, *sh2)) {
792  found2 = true;
793  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: A RecHit associated to the correct Simhit exists at lpos "
794  << (**rh).localPosition()
795  << " gpos " << (**rh).globalPosition()
796  << " Chi2 " << theChi2->estimate( secondDetState, **rh).second
797  ;
798  }
799  }
800  if (!found2) {
801  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: there is no RecHit associated to the correct SimHit." ;
802  LogTrace("CkfDebugger") << " There are " << recHits.size() << " RecHits in the simHit DetUnit" ;
803  for ( MeasurementDet::RecHitContainer::const_iterator rh = recHits.begin(); rh != recHits.end(); rh++) {
804  LogTrace("CkfDebugger") << "Non-associated RecHit at pos " << (**rh).localPosition() ;
805  }
806  }
807  }
808  }
809 
810  MeasurementDet::RecHitContainer gluedHits;
811  if (found && found2) {
812  // look in the glued det
813  DetId gid = gluedId( subdet);
814  const MeasurementDet* gluedDet = theMeasurementTracker->idToDet( gid);
815  if ( gluedDet == 0) {
816  edm::LogVerbatim("CkfDebugger") << "CkfDebugger ERROR: glued MeasurementDet not found!" ;
817  glued_det_not_fuond++;
818  return std::pair<CTTRHp, double>((CTTRHp)(0),-4);
819  }
820 
821  TSOS gluedDetState = theForwardPropagator->propagate( startingState, gluedDet->surface());
822  if (!gluedDetState.isValid()) {
823  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: propagation failed from state " << startingState << " to det surface "
824  << gluedDet->position() ;
825  propagation++;
826  return std::pair<CTTRHp, double>((CTTRHp)(0),-1);
827  }
828 
829  gluedHits = gluedDet->recHits( gluedDetState);
830  edm::LogVerbatim("CkfDebugger") << "CkfDebugger: the GluedDet returned " << gluedHits.size() << " hits" ;
831  if (gluedHits.size()==0){
832  edm::LogVerbatim("CkfDebugger") << "Found and associated mono and stereo recHits but not matched!!!" ;
833  matched_not_found++;
834  return std::pair<CTTRHp, double>((CTTRHp)(0),-5);
835  }
836  bool found3 = false;
837  for ( MeasurementDet::RecHitContainer::const_iterator rh = gluedHits.begin(); rh != gluedHits.end(); rh++) {
838  if ( associated( *rh, sh) && associated( *rh, *sh2)) {
839  double chi2 = theChi2->estimate(gluedDetState, **rh).second;
840  edm::LogVerbatim("CkfDebugger") << "Matched hit at lpos " << (**rh).localPosition()
841  << " gpos " << (**rh).globalPosition()
842  << " has Chi2 " << chi2
843  ;
844  result = std::pair<CTTRHp, double>(&**rh,chi2);
845  found3 = true;
846  if (chi2>30){
847  LogTrace("CkfDebugger") << "rh->parameters()=" << (*rh)->parameters() ;
848  LogTrace("CkfDebugger") << "rh->parametersError()=" << (*rh)->parametersError() ;
849  MeasurementExtractor me(gluedDetState);
850  AlgebraicVector r((*rh)->parameters() - me.measuredParameters(**rh));
851  LogTrace("CkfDebugger") << "me.measuredParameters(**rh)=" << me.measuredParameters(**rh) ;
852  LogTrace("CkfDebugger") << "me.measuredError(**rh)=" << me.measuredError(**rh) ;
853  AlgebraicSymMatrix R((*rh)->parametersError() + me.measuredError(**rh));
854  LogTrace("CkfDebugger") << "r=" << r ;
855  LogTrace("CkfDebugger") << "R=" << R ;
856  int ierr;
857  R.invert(ierr);
858  LogTrace("CkfDebugger") << "R(-1)=" << R ;
859  LogTrace("CkfDebugger") << "chi2=" << R.similarity(r) ;
860  }
861  break;
862  }
863  }
864  if (found3) return result;
865  else {
866  edm::LogVerbatim("CkfDebugger") << "Found and associated mono and stereo recHits. Matched found but not associated!!!" ;
867  matched_not_associated++;
868  return std::pair<CTTRHp, double>((CTTRHp)(0),-6);
869  }
870  }
871  else if ( (found && !found2) || (!found && found2) ) {
872  edm::LogVerbatim("CkfDebugger") << "Only one component is found" ;
873  only_one_component++;
874  return std::pair<CTTRHp, double>((CTTRHp)(0),-7);
875  }
876  else {
877  edm::LogVerbatim("CkfDebugger") << "No component is found" ;
878  no_component++;
879  return std::pair<CTTRHp, double>((CTTRHp)(0),-2);
880  }
881  other++;
882  return std::pair<CTTRHp, double>((CTTRHp)(0),0);//other
883 }
884 
885 const PSimHit* CkfDebugger::pSimHit(unsigned int tkId, DetId detId)
886 {
887  for (std::vector<PSimHit*>::iterator shi=idHitsMap[tkId].begin(); shi!=idHitsMap[tkId].end(); ++shi) {
888  if ( (*shi)->detUnitId() == detId.rawId() &&
889  //(shi)->trackId() == tkId &&
890  goodSimHit(**shi) ) {
891  return (*shi);
892  }
893  }
894  return 0;
895 }
896 
897 int CkfDebugger::analyseRecHitNotFound(const Trajectory& traj, CTTRHp correctRecHit)
898 {
899  unsigned int correctDetId = correctRecHit->det()->geographicalId().rawId();
900  int correctLayId = layer(correctRecHit->det());
901  LogTrace("CkfDebugger") << "correct layer id=" << correctLayId ;
902 
903  TSOS currentState( traj.lastMeasurement().updatedState() );
904  std::vector<const DetLayer*> nl = traj.lastLayer()->nextLayers( *currentState.freeState(),traj.direction() );
905  if (nl.empty()) {
906  edm::LogVerbatim("CkfDebugger") << "no compatible layers" ;
907  no_layer++;return 2;
908  }
909 
910  TkLayerLess lless;//FIXME - was lless(traj.direction())
911  const DetLayer* detLayer = 0;
912  bool navLayerAfter = false;
913  bool test = false;
914  for (std::vector<const DetLayer*>::iterator il = nl.begin(); il != nl.end(); il++) {
915  if ( dynamic_cast<const BarrelDetLayer*>(*il) ){
916  const BarrelDetLayer* pbl = dynamic_cast<const BarrelDetLayer*>(*il);
917  LogTrace("CkfDebugger") << "pbl->specificSurface().bounds().length()=" << pbl->specificSurface().bounds().length() ;
918  LogTrace("CkfDebugger") << "pbl->specificSurface().bounds().width()=" << pbl->specificSurface().bounds().width() ;
919  }
920  int layId = layer(((*(*il)->basicComponents().begin())));
921  LogTrace("CkfDebugger") << " subdet=" << (*(*il)->basicComponents().begin())->geographicalId().subdetId() << "layer id=" << layId ;
922  if (layId==correctLayId) {
923  test = true;
924  detLayer = &**il;
925  break;
926  }
927  if ( lless( *il, theGeomSearchTracker->detLayer(correctRecHit->det()->geographicalId()) ))
928  navLayerAfter = true; //it is enough that only one layer is after the correct one?
929  }
930 
931  if (test) {
932  edm::LogVerbatim("CkfDebugger") << "correct layer taken into account. layer id: " << correctLayId ;
933  } else if (navLayerAfter){
934  edm::LogVerbatim("CkfDebugger")<< "SimHit layer after the layers returned by Navigation.";
935  edm::LogVerbatim("CkfDebugger")<< "Probably a missing SimHit." ;
936  edm::LogVerbatim("CkfDebugger")<< "check: " << (correctRecHit->det()->geographicalId().subdetId()) << " " << (layer(correctRecHit->det()));
937  dump6[pair<int,int>((correctRecHit->det()->geographicalId().subdetId()-1),(layer(correctRecHit->det()))-1)]++;
938  no_sim_hit++;return 16;
939  }
940  else {
941  edm::LogVerbatim("CkfDebugger") << "correct layer NOT taken into account. correct layer id: " << correctLayId ;
942  layer_not_found++;
943  return 3;
944  }
945 
946  typedef DetLayer::DetWithState DetWithState;
947  std::vector<DetWithState> compatDets = detLayer->compatibleDets(currentState,*theForwardPropagator,*theChi2);
948  // LogTrace("CkfDebugger") << "DEBUGGER" ;
949  // LogTrace("CkfDebugger") << "runned compatDets." ;
950  // LogTrace("CkfDebugger") << "started from the following TSOS:" ;
951  // LogTrace("CkfDebugger") << currentState ;
952  // LogTrace("CkfDebugger") << "number of dets found=" << compatDets.size() ;
953  // for (std::vector<DetWithState>::iterator det=compatDets.begin();det!=compatDets.end();det++){
954  // unsigned int detId = det->first->geographicalId().rawId();
955  // LogTrace("CkfDebugger") << "detId=" << detId ;
956  // }
957  bool test2 = false;
958  for (std::vector<DetWithState>::iterator det=compatDets.begin();det!=compatDets.end();det++){
959  unsigned int detId = det->first->geographicalId().rawId();
960  // LogTrace("CkfDebugger") << "detId=" << detId
961  // << "\ncorrectRecHit->det()->geographicalId().rawId()=" << correctRecHit->det()->geographicalId().rawId()
962  // << "\ngluedId(correctRecHit->det()->geographicalId()).rawId()=" << gluedId(correctRecHit->det()->geographicalId()).rawId()
963  // ;
964  if (detId==gluedId(correctRecHit->det()->geographicalId()).rawId()) {
965  test2=true;
966  break;
967  }
968  }
969 
970  if (test2){
971  edm::LogVerbatim("CkfDebugger") << "correct det taken into account. correctDetId is: " << correctDetId
972  << ". please check chi2." ;
973  return 5;
974  }
975  else {
976  edm::LogVerbatim("CkfDebugger") << "correct det NOT taken into account. correctDetId: " << correctDetId ;
977  det_not_found++;return 4;
978  }
979 
980 }
981 
982 double CkfDebugger::testSeed(CTTRHp recHit1, CTTRHp recHit2, TSOS state){
983  //edm::LogVerbatim("CkfDebugger") << "CkfDebugger::testSeed";
984  //test Deltas
985  const std::vector<PSimHit>& pSimHitVec1 = hitAssociator->associateHit(*recHit1->hit());
986  const std::vector<PSimHit>& pSimHitVec2 = hitAssociator->associateHit(*recHit2->hit());
987 
988  if ( pSimHitVec1.size()==0 || pSimHitVec2.size()==0 || hasDelta(&(*pSimHitVec1.begin())) || hasDelta(&(*pSimHitVec2.begin())) ) {
989  edm::LogVerbatim("CkfDebugger") << "Seed has delta or problems" ;
990  return -1;
991  }
992 
993  // LogTrace("CkfDebugger") << "state=\n" << state ;
994  // double stlp1 = state.localParameters().vector()[0];
995  // double stlp2 = state.localParameters().vector()[1];
996  // double stlp3 = state.localParameters().vector()[2];
997  // double stlp4 = state.localParameters().vector()[3];
998  // double stlp5 = state.localParameters().vector()[4];
999 
1000  if (pSimHitVec2.size()!=0) {
1001  const PSimHit& simHit = *pSimHitVec2.begin();
1002 
1003  double shlp1 = -1/simHit.momentumAtEntry().mag();
1004  double shlp2 = simHit.momentumAtEntry().x()/simHit.momentumAtEntry().z();
1005  double shlp3 = simHit.momentumAtEntry().y()/simHit.momentumAtEntry().z();
1006  double shlp4 = simHit.localPosition().x();
1007  double shlp5 = simHit.localPosition().y();
1008  AlgebraicVector5 v;
1009  v[0] = shlp1;
1010  v[1] = shlp2;
1011  v[2] = shlp3;
1012  v[3] = shlp4;
1013  v[4] = shlp5;
1014 
1015  // LogTrace("CkfDebugger") << "simHit.localPosition()=" << simHit.localPosition() ;
1016  // LogTrace("CkfDebugger") << "simHit.momentumAtEntry()=" << simHit.momentumAtEntry() ;
1017  // LogTrace("CkfDebugger") << "recHit2->localPosition()=" << recHit2->localPosition() ;
1018  // LogTrace("CkfDebugger") << "recHit2->localPositionError()=" << recHit2->localPositionError() ;
1019  // LogTrace("CkfDebugger") << "state.localPosition()=" << state.localPosition() ;
1020  // LogTrace("CkfDebugger") << "state.localError().positionError()=" << state.localError().positionError() ;
1021 
1022  // LogTrace("CkfDebugger") << "pullx(sh-rh)=" << (simHit.localPosition().x()-recHit2->localPosition().x())/sqrt(recHit2->localPositionError().xx()) ;
1023  // LogTrace("CkfDebugger") << "pullx(sh-st)=" << (simHit.localPosition().x()-state.localPosition().x())/sqrt(state.localError().positionError().xx()) ;
1024  // LogTrace("CkfDebugger") << "pullx(st-rh)=" << (state.localPosition().x()-recHit2->localPosition().x())/
1025  // sqrt(recHit2->localPositionError().xx()+state.localError().positionError().xx()) ;
1026 
1027  // LogTrace("CkfDebugger") << "local parameters" ;
1028  // LogTrace("CkfDebugger") << left;
1029  // LogTrace("CkfDebugger") << setw(15) << stlp1 << setw(15) << shlp1 << setw(15) << sqrt(state.localError().matrix()[0][0])
1030  // << setw(15) << (stlp1-shlp1)/stlp1 << setw(15) << (stlp1-shlp1)/sqrt(state.localError().matrix()[0][0]) ;
1031  // LogTrace("CkfDebugger") << setw(15) << stlp2 << setw(15) << shlp2 << setw(15) << sqrt(state.localError().matrix()[1][1])
1032  // << setw(15) << (stlp2-shlp2)/stlp2 << setw(15) << (stlp2-shlp2)/sqrt(state.localError().matrix()[1][1]) ;
1033  // LogTrace("CkfDebugger") << setw(15) << stlp3 << setw(15) << shlp3 << setw(15) << sqrt(state.localError().matrix()[2][2])
1034  // << setw(15) << (stlp3-shlp3)/stlp3 << setw(15) << (stlp3-shlp3)/sqrt(state.localError().matrix()[2][2]) ;
1035  // LogTrace("CkfDebugger") << setw(15) << stlp4 << setw(15) << shlp4 << setw(15) << sqrt(state.localError().matrix()[3][3])
1036  // << setw(15) << (stlp4-shlp4)/stlp4 << setw(15) << (stlp4-shlp4)/sqrt(state.localError().matrix()[3][3]) ;
1037  // LogTrace("CkfDebugger") << setw(15) << stlp5 << setw(15) << shlp5 << setw(15) << sqrt(state.localError().matrix()[4][4]) <<
1038  // setw(15) << (stlp5-shlp5)/stlp5 << setw(15) << (stlp5-shlp5)/sqrt(state.localError().matrix()[4][4]) ;
1039 
1040  AlgebraicSymMatrix55 R = state.localError().matrix();
1041  R.Invert();
1042  double chi2 = ROOT::Math::Similarity(v-state.localParameters().vector(), R);
1043  LogTrace("CkfDebugger") << "chi2=" << chi2 ;
1044  return chi2;
1045  }
1046 
1047  return 0;//fixme
1048 
1049 }
1050 
1051 
1052 CkfDebugger::~CkfDebugger(){
1053  for (int it=0; it!=((int)(dump.size())); it++)
1054  edm::LogVerbatim("CkfDebugger") << "dump " << it << " " << dump[it] ;
1055 
1056  edm::LogVerbatim("CkfDebugger") ;
1057  edm::LogVerbatim("CkfDebugger") << "seedWithDelta=" << ((double)seedWithDelta/totSeeds) ;
1058  edm::LogVerbatim("CkfDebugger") << "problems=" << ((double)problems/totSeeds) ;
1059  edm::LogVerbatim("CkfDebugger") << "no_sim_hit=" << ((double)no_sim_hit/totSeeds) ;
1060  edm::LogVerbatim("CkfDebugger") << "no_layer=" << ((double)no_layer/totSeeds) ;
1061  edm::LogVerbatim("CkfDebugger") << "layer_not_found=" << ((double)layer_not_found/totSeeds) ;
1062  edm::LogVerbatim("CkfDebugger") << "det_not_found=" << ((double)det_not_found/totSeeds) ;
1063  edm::LogVerbatim("CkfDebugger") << "chi2gt30=" << ((double)chi2gt30/totSeeds) ;
1064  edm::LogVerbatim("CkfDebugger") << "chi2gt30deltaSeed=" << ((double)chi2gt30deltaSeed/totSeeds) ;
1065  edm::LogVerbatim("CkfDebugger") << "chi2gt30delta=" << ((double)chi2gt30delta/totSeeds) ;
1066  edm::LogVerbatim("CkfDebugger") << "chi2ls30=" << ((double)chi2ls30/totSeeds) ;
1067  edm::LogVerbatim("CkfDebugger") << "simple_hit_not_found=" << ((double)simple_hit_not_found/totSeeds) ;
1068  edm::LogVerbatim("CkfDebugger") << "no_component=" << ((double)no_component/totSeeds) ;
1069  edm::LogVerbatim("CkfDebugger") << "only_one_component=" << ((double)only_one_component/totSeeds) ;
1070  edm::LogVerbatim("CkfDebugger") << "matched_not_found=" << ((double)matched_not_found/totSeeds) ;
1071  edm::LogVerbatim("CkfDebugger") << "matched_not_associated=" << ((double)matched_not_associated/totSeeds) ;
1072  edm::LogVerbatim("CkfDebugger") << "partner_det_not_fuond=" << ((double)partner_det_not_fuond/totSeeds) ;
1073  edm::LogVerbatim("CkfDebugger") << "glued_det_not_fuond=" << ((double)glued_det_not_fuond/totSeeds) ;
1074  edm::LogVerbatim("CkfDebugger") << "propagation=" << ((double)propagation/totSeeds) ;
1075  edm::LogVerbatim("CkfDebugger") << "other=" << ((double)other/totSeeds) ;
1076  edm::LogVerbatim("CkfDebugger") << "totchi2gt30=" << ((double)totchi2gt30/totSeeds) ;
1077  edm::LogVerbatim("CkfDebugger") << "totSeeds=" << totSeeds ;
1078  edm::LogVerbatim("CkfDebugger") ;
1079 
1080  edm::LogVerbatim("CkfDebugger") << "layer navigation problems:" ;
1081  for (int i=0; i!=6; i++)
1082  for (int j=0; j!=9; j++){
1083  if (i==0 && j>2) break;
1084  if (i==1 && j>1) break;
1085  if (i==2 && j>3) break;
1086  if (i==3 && j>2) break;
1087  if (i==4 && j>5) break;
1088  if (i==5 && j>8) break;
1089  edm::LogVerbatim("CkfDebugger") << "det=" << i+1 << " lay=" << j+1 << " " << dump2[pair<int,int>(i,j)] ;
1090  }
1091  edm::LogVerbatim("CkfDebugger") << "\nlayer with hit having chi2>30:" ;
1092  for (int i=0; i!=6; i++)
1093  for (int j=0; j!=9; j++){
1094  if (i==0 && j>2) break;
1095  if (i==1 && j>1) break;
1096  if (i==2 && j>3) break;
1097  if (i==3 && j>2) break;
1098  if (i==4 && j>5) break;
1099  if (i==5 && j>8) break;
1100  edm::LogVerbatim("CkfDebugger") << "det=" << i+1 << " lay=" << j+1 << " " << dump3[pair<int,int>(i,j)] ;
1101  }
1102  edm::LogVerbatim("CkfDebugger") << "\nlayer with hit having chi2>30 for delta rays:" ;
1103  for (int i=0; i!=6; i++)
1104  for (int j=0; j!=9; j++){
1105  if (i==0 && j>2) break;
1106  if (i==1 && j>1) break;
1107  if (i==2 && j>3) break;
1108  if (i==3 && j>2) break;
1109  if (i==4 && j>5) break;
1110  if (i==5 && j>8) break;
1111  edm::LogVerbatim("CkfDebugger") << "det=" << i+1 << " lay=" << j+1 << " " << dump5[pair<int,int>(i,j)] ;
1112  }
1113  edm::LogVerbatim("CkfDebugger") << "\nlayer with det not found:" ;
1114  for (int i=0; i!=6; i++)
1115  for (int j=0; j!=9; j++){
1116  if (i==0 && j>2) break;
1117  if (i==1 && j>1) break;
1118  if (i==2 && j>3) break;
1119  if (i==3 && j>2) break;
1120  if (i==4 && j>5) break;
1121  if (i==5 && j>8) break;
1122  edm::LogVerbatim("CkfDebugger") << "det=" << i+1 << " lay=" << j+1 << " " << dump4[pair<int,int>(i,j)] ;
1123  }
1124  edm::LogVerbatim("CkfDebugger") << "\nlayer with correct RecHit after missing Sim Hit:" ;
1125  for (int i=0; i!=6; i++)
1126  for (int j=0; j!=9; j++){
1127  if (i==0 && j>2) break;
1128  if (i==1 && j>1) break;
1129  if (i==2 && j>3) break;
1130  if (i==3 && j>2) break;
1131  if (i==4 && j>5) break;
1132  if (i==5 && j>8) break;
1133  edm::LogVerbatim("CkfDebugger") << "det=" << i+1 << " lay=" << j+1 << " " << dump6[pair<int,int>(i,j)] ;
1134  }
1135  hchi2seedAll->Write();
1136  hchi2seedProb->Write();
1137  std::stringstream title;
1138  for (int i=0; i!=6; i++)
1139  for (int j=0; j!=9; j++){
1140  if (i==0 && j>2) break;
1141  if (i==1 && j>1) break;
1142  if (i==2 && j>3) break;
1143  if (i==3 && j>2) break;
1144  if (i==4 && j>5) break;
1145  if (i==5 && j>8) break;
1146  title.str("");
1147  title << "pullX_" << i+1 << "-" << j+1 << "_sh-rh";
1148  hPullX_shrh[title.str()]->Write();
1149  title.str("");
1150  title << "pullY_" << i+1 << "-" << j+1 << "_sh-rh";
1151  hPullY_shrh[title.str()]->Write();
1152  title.str("");
1153  title << "pullX_" << i+1 << "-" << j+1 << "_sh-st";
1154  hPullX_shst[title.str()]->Write();
1155  title.str("");
1156  title << "pullY_" << i+1 << "-" << j+1 << "_sh-st";
1157  hPullY_shst[title.str()]->Write();
1158  title.str("");
1159  title << "pullX_" << i+1 << "-" << j+1 << "_st-rh";
1160  hPullX_strh[title.str()]->Write();
1161  title.str("");
1162  title << "pullY_" << i+1 << "-" << j+1 << "_st-rh";
1163  hPullY_strh[title.str()]->Write();
1164  title.str("");
1165  title << "PullGP_X_" << i+1 << "-" << j+1 << "_sh-st";
1166  hPullGP_X_shst[title.str()]->Write();
1167  title.str("");
1168  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_sh-st";
1169  hPullGP_Y_shst[title.str()]->Write();
1170  title.str("");
1171  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_sh-st";
1172  hPullGP_Z_shst[title.str()]->Write();
1173  if ( ((i==2||i==4)&&(j==0||j==1)) || (i==3||i==5) ){
1174  title.str("");
1175  title << "pullM_" << i+1 << "-" << j+1 << "_sh-rh";
1176  hPullM_shrh[title.str()]->Write();
1177  title.str("");
1178  title << "pullS_" << i+1 << "-" << j+1 << "_sh-rh";
1179  hPullS_shrh[title.str()]->Write();
1180  title.str("");
1181  title << "pullM_" << i+1 << "-" << j+1 << "_sh-st";
1182  hPullM_shst[title.str()]->Write();
1183  title.str("");
1184  title << "pullS_" << i+1 << "-" << j+1 << "_sh-st";
1185  hPullS_shst[title.str()]->Write();
1186  title.str("");
1187  title << "pullM_" << i+1 << "-" << j+1 << "_st-rh";
1188  hPullM_strh[title.str()]->Write();
1189  title.str("");
1190  title << "pullS_" << i+1 << "-" << j+1 << "_st-rh";
1191  hPullS_strh[title.str()]->Write();
1192  }
1193  }
1194  hPullGPXvsGPX_shst->Write();
1195  hPullGPXvsGPY_shst->Write();
1196  hPullGPXvsGPZ_shst->Write();
1197  hPullGPXvsGPr_shst->Write();
1198  hPullGPXvsGPeta_shst->Write();
1199  hPullGPXvsGPphi_shst->Write();
1200 
1201  //file->Write();
1202  file->Close();
1203 }
MeasurementExtractor::measuredError
AlgebraicSymMatrix measuredError(const TransientTrackingRecHit &)
Definition: MeasurementExtractor.cc:9
CkfDebugger::SimHit::globalPosition
GlobalPoint globalPosition() const
Definition: CkfDebugger.h:78
MeasurementTrackerESProducer_cfi.MeasurementTracker
tuple MeasurementTracker
Definition: MeasurementTrackerESProducer_cfi.py:3
i
int i
Definition: DBlmapReader.cc:9
CkfDebugger::hPullS_shrh
std::map< std::string, TH1F * > hPullS_shrh
Definition: CkfDebugger.h:216
LocalError::xx
float xx() const
Definition: LocalError.h:24
CkfDebugger::theTopo
const TrackerTopology * theTopo
Definition: CkfDebugger.h:107
CkfDebugger::hasDelta
bool hasDelta(const PSimHit *correctHit)
Definition: CkfDebugger.h:134
CkfDebugger::analyseRecHitExistance
std::pair< CTTRHp, double > analyseRecHitExistance(const PSimHit &sh, const TSOS &startingState)
Definition: CkfDebugger.cc:675
CkfDebugger::hPullY_shst
std::map< std::string, TH1F * > hPullY_shst
Definition: CkfDebugger.h:211
CkfDebugger::correctMeas
bool correctMeas(const TM &tm, const PSimHit *correctHit) const
Definition: CkfDebugger.cc:629
CkfDebugger::printSimHits
void printSimHits(const edm::Event &iEvent)
Definition: CkfDebugger.cc:151
CkfDebugger::dump6
std::map< std::pair< int, int >, int > dump6
Definition: CkfDebugger.h:203
PV3DBase::perp
T perp() const
Definition: PV3DBase.h:72
TrajectoryMeasurement::recHit
ConstRecHitPointer const & recHit() const
Definition: TrajectoryMeasurement.h:196
CkfDebugger::theMeasurementTracker
const MeasurementTracker * theMeasurementTracker
Definition: CkfDebugger.h:105
TrajectoryStateOnSurface::localParameters
const LocalTrajectoryParameters & localParameters() const
Definition: TrajectoryStateOnSurface.h:168
CkfDebugger::computePulls
std::pair< double, double > computePulls(CTTRHp recHit, TSOS startingState)
Definition: CkfDebugger.h:160
CkfDebugger::hPullGPXvsGPeta_shst
TH2F * hPullGPXvsGPeta_shst
Definition: CkfDebugger.h:230
CkfDebugger::dumpSimHit
void dumpSimHit(const SimHit &hit) const
Definition: CkfDebugger.cc:180
PSimHit::momentumAtEntry
LocalVector momentumAtEntry() const
The momentum of the track that produced the hit, at entry point.
Definition: PSimHit.h:47
CkfDebugger::theMagField
const MagneticField * theMagField
Definition: CkfDebugger.h:100
MeasurementEstimator::estimate
virtual HitReturnType estimate(const TrajectoryStateOnSurface &ts, const TransientTrackingRecHit &hit) const =0
CkfDebugger::theTTRHBuilder
const TransientTrackingRecHitBuilder * theTTRHBuilder
Definition: CkfDebugger.h:106
GeomDet::toGlobal
GlobalPoint toGlobal(const Local2DPoint &lp) const
Conversion to the global R.F. from the R.F. of the GeomDet.
Definition: GeomDet.h:47
CkfDebugger::hPullGP_Z_shst
std::map< std::string, TH1F * > hPullGP_Z_shst
Definition: CkfDebugger.h:224
CkfDebugger::theGeomSearchTracker
const GeometricSearchTracker * theGeomSearchTracker
Definition: CkfDebugger.h:101
MeasurementDet::RecHitContainer
TransientTrackingRecHit::ConstRecHitContainer RecHitContainer
Definition: MeasurementDet.h:19
MeasurementDet::geomDet
virtual const GeomDet & geomDet() const
Definition: MeasurementDet.h:64
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
TrajectoryStateOnSurface::cartesianError
const CartesianTrajectoryError cartesianError() const
Definition: TrajectoryStateOnSurface.h:162
PV3DBase::y
T y() const
Definition: PV3DBase.h:63
CkfDebugger::hPullY_strh
std::map< std::string, TH1F * > hPullY_strh
Definition: CkfDebugger.h:213
TrajectoryStateOnSurface::globalPosition
GlobalPoint globalPosition() const
Definition: TrajectoryStateOnSurface.h:144
LargeD0_PixelPairStep_cff.propagator
tuple propagator
Definition: LargeD0_PixelPairStep_cff.py:84
CkfDebugger::hPullM_shrh
std::map< std::string, TH1F * > hPullM_shrh
Definition: CkfDebugger.h:215
AlgebraicSymMatrix55
ROOT::Math::SMatrix< double, 5, 5, ROOT::Math::MatRepSym< double, 5 > > AlgebraicSymMatrix55
Definition: AlgebraicROOTObjects.h:22
MeasurementDet::surface
const Surface & surface() const
Definition: MeasurementDet.h:66
CkfDebugger::hPullX_shst
std::map< std::string, TH1F * > hPullX_shst
Definition: CkfDebugger.h:210
CkfDebugger::partner_det_not_fuond
int partner_det_not_fuond
Definition: CkfDebugger.h:248
CkfDebugger::nextCorrectHits
std::vector< const PSimHit * > nextCorrectHits(const Trajectory &, unsigned int &)
Definition: CkfDebugger.cc:524
AlCaHLTBitMon_QueryRunRegistry.comp
string comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:248
CkfDebugger::det_not_found
int det_not_found
Definition: CkfDebugger.h:238
less_mag
Definition: HybridClusterAlgo.h:29
StripSubdetector::partnerDetId
unsigned int partnerDetId() const
Definition: StripSubdetector.h:57
TrackerHitAssociator::associateHitId
std::vector< SimHitIdpr > associateHitId(const TrackingRecHit &thit)
Definition: TrackerHitAssociator.cc:280
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:35
CkfDebugger::seedWithDelta
int seedWithDelta
Definition: CkfDebugger.h:233
CkfDebugger::hPullM_shst
std::map< std::string, TH1F * > hPullM_shst
Definition: CkfDebugger.h:217
LocalTrajectoryError::positionError
LocalError positionError() const
Definition: LocalTrajectoryError.h:82
TransientTrackingRecHit::RecHitContainer
std::vector< ConstRecHitPointer > RecHitContainer
Definition: TransientTrackingRecHit.h:28
GeometricSearchDet::compatibleDets
virtual std::vector< DetWithState > compatibleDets(const TrajectoryStateOnSurface &startingState, const Propagator &prop, const MeasurementEstimator &est) const
Definition: GeometricSearchDet.cc:26
CkfDebugger::theChi2
const MeasurementEstimator * theChi2
Definition: CkfDebugger.h:102
CkfDebugger::matched_not_found
int matched_not_found
Definition: CkfDebugger.h:246
Association.map
dictionary map
Definition: Association.py:205
CkfDebugger::SimHit::trackId
unsigned int trackId() const
Definition: CkfDebugger.h:80
Trajectory::recHits
ConstRecHitContainer recHits(bool splitting=false) const
Definition: Trajectory.cc:67
TrackingRecHitProjector::project
RecHitPointer project(const TransientTrackingRecHit &hit, const GeomDet &det, const TrajectoryStateOnSurface &ts) const
Definition: TrackingRecHitProjector.h:14
Trajectory::direction
PropagationDirection const & direction() const
Definition: Trajectory.cc:196
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:45
MeasurementDet::recHits
virtual RecHitContainer recHits(const TrajectoryStateOnSurface &) const =0
Trajectory::measurements
DataContainer const & measurements() const
Definition: Trajectory.h:215
LocalTrajectoryParameters::vector
AlgebraicVector5 vector() const
Definition: LocalTrajectoryParameters.h:120
CkfDebugger::CkfDebugger
CkfDebugger(edm::EventSetup const &es)
Definition: CkfDebugger.cc:35
CkfDebugger::analyseRecHitNotFound
int analyseRecHitNotFound(const Trajectory &, CTTRHp)
Definition: CkfDebugger.cc:897
iEvent
int iEvent
Definition: GenABIO.cc:243
TrajectoryStateOnSurface::surface
const SurfaceType & surface() const
Definition: TrajectoryStateOnSurface.h:183
StripSubdetector::glued
unsigned int glued() const
glued
Definition: StripSubdetector.h:27
PV3DBase::mag
T mag() const
Definition: PV3DBase.h:67
GeomDet::position
const Surface::PositionType & position() const
The position (origin of the R.F.)
Definition: GeomDet.h:41
CkfDebugger::hPullGP_Y_shst
std::map< std::string, TH1F * > hPullGP_Y_shst
Definition: CkfDebugger.h:223
CkfDebugger::layer
int layer(const GeomDetUnit *det)
Definition: CkfDebugger.h:150
PSimHit::timeOfFlight
float timeOfFlight() const
Definition: PSimHit.h:69
CkfDebugger::testSeed
double testSeed(CTTRHp, CTTRHp, TrajectoryStateOnSurface)
Definition: CkfDebugger.cc:982
CkfDebugger::SimHit::particleType
int particleType() const
Definition: CkfDebugger.h:87
CkfDebugger::matched_not_associated
int matched_not_associated
Definition: CkfDebugger.h:247
LocalError::yy
float yy() const
Definition: LocalError.h:26
PSimHit::localPosition
Local3DPoint localPosition() const
Definition: PSimHit.h:44
CkfDebugger::associated
bool associated(CTTRHp rechit, const PSimHit &sh) const
Definition: CkfDebugger.cc:609
CkfDebugger::only_one_component
int only_one_component
Definition: CkfDebugger.h:245
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
TransientTrackingRecHitBuilder::build
virtual RecHitPointer build(const TrackingRecHit *p) const =0
build a tracking rechit from an existing rechit
Trajectory::lastMeasurement
TrajectoryMeasurement const & lastMeasurement() const
Definition: Trajectory.h:193
PV3DBase::z
T z() const
Definition: PV3DBase.h:64
query.result
tuple result
Definition: query.py:137
CkfDebugger::assocTrackId
int assocTrackId(CTTRHp rechit) const
Definition: CkfDebugger.cc:506
CkfDebugger::hPullGPXvsGPr_shst
TH2F * hPullGPXvsGPr_shst
Definition: CkfDebugger.h:229
CkfDebugger::glued_det_not_fuond
int glued_det_not_fuond
Definition: CkfDebugger.h:249
CkfDebugger::hchi2seedAll
TH1F * hchi2seedAll
Definition: CkfDebugger.h:206
j
int j
Definition: DBlmapReader.cc:9
GeometricSearchTracker::detLayer
const DetLayer * detLayer(const DetId &id) const
obsolete method. Use idToLayer() instead.
Definition: GeometricSearchTracker.h:53
CkfDebugger::hitAssociator
TrackerHitAssociator * hitAssociator
Definition: CkfDebugger.h:104
CkfDebugger::hPullGPXvsGPY_shst
TH2F * hPullGPXvsGPY_shst
Definition: CkfDebugger.h:227
GeomDet::geographicalId
DetId geographicalId() const
The label of this GeomDet.
Definition: GeomDet.h:72
CTTRHp
TransientTrackingRecHit::ConstRecHitPointer CTTRHp
Definition: CkfDebugger.h:39
patCandidatesForDimuonsSequences_cff.tracker
tuple tracker
Definition: patCandidatesForDimuonsSequences_cff.py:69
LocalTrajectoryError::matrix
const AlgebraicSymMatrix55 & matrix() const
Definition: LocalTrajectoryError.h:60
CkfDebugger::hPullX_shrh
std::map< std::string, TH1F * > hPullX_shrh
Definition: CkfDebugger.h:208
TrajectoryStateOnSurface::localError
const LocalTrajectoryError & localError() const
Definition: TrajectoryStateOnSurface.h:180
DetWithState
std::pair< const GeomDet *, TrajectoryStateOnSurface > DetWithState
Definition: RPCRecHitFilter.h:58
PSimHit::pabs
float pabs() const
fast and more accurate access to momentumAtEntry().mag()
Definition: PSimHit.h:63
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:39
CkfDebugger::hPullGP_X_shst
std::map< std::string, TH1F * > hPullGP_X_shst
Definition: CkfDebugger.h:222
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
MeasurementDet::position
const Surface::PositionType & position() const
Definition: MeasurementDet.h:68
Propagator::propagate
virtual TrajectoryStateOnSurface propagate(const FreeTrajectoryState &, const Surface &) const
Definition: Propagator.cc:12
CkfDebugger::theForwardPropagator
const Propagator * theForwardPropagator
Definition: CkfDebugger.h:103
CkfDebugger::correctTrajectory
bool correctTrajectory(const Trajectory &, unsigned int &) const
Definition: CkfDebugger.cc:464
CkfDebugger::chi2gt30deltaSeed
int chi2gt30deltaSeed
Definition: CkfDebugger.h:241
DetId
Definition: DetId.h:20
AlgebraicVector
CLHEP::HepVector AlgebraicVector
Definition: AlgebraicObjects.h:12
CkfDebugger::hPullM_strh
std::map< std::string, TH1F * > hPullM_strh
Definition: CkfDebugger.h:219
Trajectory::firstMeasurement
TrajectoryMeasurement const & firstMeasurement() const
Definition: Trajectory.h:206
CkfDebugger::SimHit::pabs
float pabs() const
Definition: CkfDebugger.h:84
CkfDebugger::hPullS_shst
std::map< std::string, TH1F * > hPullS_shst
Definition: CkfDebugger.h:218
AlgebraicVector5
ROOT::Math::SVector< double, 5 > AlgebraicVector5
Definition: AlgebraicROOTObjects.h:14
CkfDebugger::simple_hit_not_found
int simple_hit_not_found
Definition: CkfDebugger.h:243
CartesianTrajectoryError::position
const GlobalError position() const
Position error submatrix.
Definition: CartesianTrajectoryError.h:46
MeasurementExtractor::measuredParameters
AlgebraicVector measuredParameters(const TransientTrackingRecHit &)
Definition: MeasurementExtractor.cc:3
Trajectory::RecHitContainer
ConstRecHitContainer RecHitContainer
Definition: Trajectory.h:44
CkfDebugger::dump2
std::map< std::pair< int, int >, int > dump2
Definition: CkfDebugger.h:199
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
run_regression.test
tuple test
Definition: run_regression.py:386
CkfDebugger::hPullS_strh
std::map< std::string, TH1F * > hPullS_strh
Definition: CkfDebugger.h:220
PSimHit::processType
unsigned short processType() const
Definition: PSimHit.h:118
CkfDebugger::file
TFile * file
Definition: CkfDebugger.h:205
evf::utils::state
char state
Definition: procUtils.cc:75
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:76
TrackerGeometry::idToDetUnit
virtual const GeomDetUnit * idToDetUnit(DetId) const
Return the pointer to the GeomDetUnit corresponding to a given DetId.
Definition: TrackerGeometry.cc:126
CkfDebugger::theTrackerGeom
const TrackerGeometry * theTrackerGeom
Definition: CkfDebugger.h:99
CkfDebugger::det
const GeomDetUnit * det(const PSimHit *sh) const
Definition: CkfDebugger.h:148
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
begin
#define begin
Definition: vmac.h:31
CkfDebugger::goodSimHit
bool goodSimHit(const PSimHit &sh) const
Definition: CkfDebugger.cc:602
CkfDebugger::hPullY_shrh
std::map< std::string, TH1F * > hPullY_shrh
Definition: CkfDebugger.h:209
CkfDebugger::hPullGPXvsGPX_shst
TH2F * hPullGPXvsGPX_shst
Definition: CkfDebugger.h:226
AlgebraicSymMatrix
CLHEP::HepSymMatrix AlgebraicSymMatrix
Definition: AlgebraicObjects.h:14
TrackerHitAssociator::associateHit
std::vector< PSimHit > associateHit(const TrackingRecHit &thit)
Definition: TrackerHitAssociator.cc:120
PSimHit::trackId
unsigned int trackId() const
Definition: PSimHit.h:102
Trajectory::lastLayer
const DetLayer * lastLayer() const
Redundant method, returns the layer of lastMeasurement() .
Definition: Trajectory.h:296
CkfDebugger::hPullX_strh
std::map< std::string, TH1F * > hPullX_strh
Definition: CkfDebugger.h:212
GeometricSearchDet::DetWithState
std::pair< const GeomDet *, TrajectoryStateOnSurface > DetWithState
Definition: GeometricSearchDet.h:21
CkfDebugger::layer_not_found
int layer_not_found
Definition: CkfDebugger.h:237
CkfDebugger::no_component
int no_component
Definition: CkfDebugger.h:244
TrajectoryMeasurement::updatedState
TrajectoryStateOnSurface const & updatedState() const
Definition: TrajectoryMeasurement.h:184
CkfDebugger::hPullGPXvsGPZ_shst
TH2F * hPullGPXvsGPZ_shst
Definition: CkfDebugger.h:228
CkfDebugger::dump
std::vector< int > dump
Definition: CkfDebugger.h:198
CkfDebugger::pSimHit
const PSimHit * pSimHit(unsigned int tkId, DetId detId)
Definition: CkfDebugger.cc:885
CkfDebugger::dump5
std::map< std::pair< int, int >, int > dump5
Definition: CkfDebugger.h:202
CkfDebugger::position
Global3DPoint position(const PSimHit *sh) const
Definition: CkfDebugger.h:145
CkfDebugger::idHitsMap
std::map< unsigned int, std::vector< PSimHit * > > idHitsMap
Definition: CkfDebugger.h:109
PV3DBase::x
T x() const
Definition: PV3DBase.h:62
BarrelDetLayer::specificSurface
virtual const BoundCylinder & specificSurface() const GCC11_FINAL
Extension of the interface.
Definition: BarrelDetLayer.h:45
gluedId
DetId gluedId(const DetId &du)
Definition: GluedDetFromDetUnit.h:3
CkfDebugger::analyseCompatibleMeasurements
bool analyseCompatibleMeasurements(const Trajectory &, const std::vector< TrajectoryMeasurement > &, const MeasurementTracker *, const Propagator *, const Chi2MeasurementEstimatorBase *, const TransientTrackingRecHitBuilder *)
Definition: CkfDebugger.cc:192
DetLayer::nextLayers
std::vector< const DetLayer * > nextLayers(Args &&...args) const
Definition: DetLayer.h:60
CkfDebugger::hchi2seedProb
TH1F * hchi2seedProb
Definition: CkfDebugger.h:206
PSimHit::detUnitId
unsigned int detUnitId() const
Definition: PSimHit.h:93
CkfDebugger::dump4
std::map< std::pair< int, int >, int > dump4
Definition: CkfDebugger.h:201
CkfDebugger::hPullGPXvsGPphi_shst
TH2F * hPullGPXvsGPphi_shst
Definition: CkfDebugger.h:231
CkfDebugger::chi2gt30delta
int chi2gt30delta
Definition: CkfDebugger.h:240
CkfDebugger::dump3
std::map< std::pair< int, int >, int > dump3
Definition: CkfDebugger.h:200
Generated for CMSSW Reference Manual by doxygen 1.8.5