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