CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
TestHits.cc
Go to the documentation of this file.
1 #include "RecoTracker/DebugTools/interface/TestHits.h"
2 
3 #include "FWCore/MessageLogger/interface/MessageLogger.h"
4 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"
5 #include "TrackingTools/TrajectoryState/interface/TrajectoryStateTransform.h"
6 #include "TrackingTools/PatternTools/interface/TransverseImpactPointExtrapolator.h"
7 #include "TrackingTools/TransientTrack/interface/TransientTrack.h"
8 #include "RecoTracker/TrackProducer/interface/TrackingRecHitLessFromGlobalPosition.h"
9 #include "TrackingTools/PatternTools/interface/TSCPBuilderNoMaterial.h"
10 #include "TrackingTools/TrackFitters/interface/TrajectoryStateCombiner.h"
11 #include <TDirectory.h>
12 #include "Geometry/TrackerGeometryBuilder/interface/GluedGeomDet.h"
13 #include "TrackingTools/TrackFitters/interface/KFTrajectoryFitter.h"
14 
15 #include "TrackingTools/TrackFitters/interface/TrajectoryFitter.h"
16 #include "DataFormats/SiStripDetId/interface/TIBDetId.h"
17 #include "DataFormats/SiStripDetId/interface/TOBDetId.h"
18 #include "DataFormats/SiStripDetId/interface/TIDDetId.h"
19 #include "DataFormats/SiStripDetId/interface/TECDetId.h"
20 #include "DataFormats/SiPixelDetId/interface/PXBDetId.h"
21 #include "DataFormats/SiPixelDetId/interface/PXFDetId.h"
22 
23 typedef TrajectoryStateOnSurface TSOS;
24 typedef TransientTrackingRecHit::ConstRecHitPointer CTTRHp;
25 using namespace std;
26 using namespace edm;
27 
28 TestHits::TestHits(const edm::ParameterSet& iConfig):
29  conf_(iConfig){
30  LogTrace("TestHits") << conf_<< std::endl;
31  propagatorName = conf_.getParameter<std::string>("Propagator");
32  builderName = conf_.getParameter<std::string>("TTRHBuilder");
33  srcName = conf_.getParameter<std::string>("src");
34  fname = conf_.getParameter<std::string>("Fitter");
35  mineta = conf_.getParameter<double>("mineta");
36  maxeta = conf_.getParameter<double>("maxeta");
37 }
38 
39 TestHits::~TestHits(){}
40 
41 void TestHits::beginRun(edm::Run & run, const edm::EventSetup& iSetup)
42 {
43 
44  iSetup.get<TrackerDigiGeometryRecord>().get(theG);
45  iSetup.get<IdealMagneticFieldRecord>().get(theMF);
46  iSetup.get<TrackingComponentsRecord>().get(propagatorName,thePropagator);
47  iSetup.get<TransientRecHitRecord>().get(builderName,theBuilder);
48  iSetup.get<TrajectoryFitter::Record>().get(fname, fit);
49 
50  file = new TFile("testhits.root","recreate");
51  for (int i=0; i!=6; i++)
52  for (int j=0; j!=9; j++){
53  if (i==0 && j>2) break;
54  if (i==1 && j>1) break;
55  if (i==2 && j>3) break;
56  if (i==3 && j>2) break;
57  if (i==4 && j>5) break;
58  if (i==5 && j>8) break;
59  title.str("");
60  title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts";
61  hPullGP_X_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
62  title.str("");
63  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts";
64  hPullGP_Y_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
65  title.str("");
66  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts";
67  hPullGP_Z_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
68  title.str("");
69  title << "Chi2Increment_" << i+1 << "-" << j+1;
70  hChi2Increment[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,0,100);
71 
72  title.str("");
73  title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts";
74  hPullGM_X_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
75  title.str("");
76  title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts";
77  hPullGM_Y_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
78  title.str("");
79  title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts";
80  hPullGM_Z_ts[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
81 
82  title.str("");
83  title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr";
84  hPullGP_X_tr[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
85  title.str("");
86  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr";
87  hPullGP_Y_tr[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
88  title.str("");
89  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr";
90  hPullGP_Z_tr[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
91 
92  title.str("");
93  title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs";
94  hPullGP_X_rs[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
95  title.str("");
96  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs";
97  hPullGP_Y_rs[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
98  title.str("");
99  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs";
100  hPullGP_Z_rs[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
101 
102  if ( ((i==2||i==4)&&(j==0||j==1)) || (i==3||i==5) ){
103  //mono
104  title.str("");
105  title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts_mono";
106  hPullGP_X_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
107  title.str("");
108  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts_mono";
109  hPullGP_Y_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
110  title.str("");
111  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts_mono";
112  hPullGP_Z_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
113 
114  title.str("");
115  title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts_mono";
116  hPullGM_X_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
117  title.str("");
118  title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts_mono";
119  hPullGM_Y_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
120  title.str("");
121  title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts_mono";
122  hPullGM_Z_ts_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
123 
124  title.str("");
125  title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr_mono";
126  hPullGP_X_tr_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
127  title.str("");
128  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr_mono";
129  hPullGP_Y_tr_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
130  title.str("");
131  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr_mono";
132  hPullGP_Z_tr_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
133 
134  title.str("");
135  title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs_mono";
136  hPullGP_X_rs_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
137  title.str("");
138  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs_mono";
139  hPullGP_Y_rs_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
140  title.str("");
141  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs_mono";
142  hPullGP_Z_rs_mono[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
143 
144  //stereo
145  title.str("");
146  title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts_stereo";
147  hPullGP_X_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
148  title.str("");
149  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts_stereo";
150  hPullGP_Y_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
151  title.str("");
152  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts_stereo";
153  hPullGP_Z_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
154 
155  title.str("");
156  title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts_stereo";
157  hPullGM_X_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
158  title.str("");
159  title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts_stereo";
160  hPullGM_Y_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
161  title.str("");
162  title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts_stereo";
163  hPullGM_Z_ts_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
164 
165  title.str("");
166  title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr_stereo";
167  hPullGP_X_tr_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
168  title.str("");
169  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr_stereo";
170  hPullGP_Y_tr_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
171  title.str("");
172  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr_stereo";
173  hPullGP_Z_tr_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
174 
175  title.str("");
176  title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs_stereo";
177  hPullGP_X_rs_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
178  title.str("");
179  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs_stereo";
180  hPullGP_Y_rs_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
181  title.str("");
182  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs_stereo";
183  hPullGP_Z_rs_stereo[title.str()] = new TH1F(title.str().c_str(),title.str().c_str(),1000,-50,50);
184  }
185  }
186  hTotChi2Increment = new TH1F("TotChi2Increment","TotChi2Increment",1000,0,100);
187  hProcess_vs_Chi2 = new TH2F("Process_vs_Chi2","Process_vs_Chi2",1000,0,100,17,-0.5,16.5);
188  hClsize_vs_Chi2 = new TH2F("Clsize_vs_Chi2","Clsize_vs_Chi2",1000,0,100,17,-0.5,16.5);
189 }
190 
191 
192 void TestHits::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
193 {
194  LogTrace("TestHits") << "\nnew event";
195 
196  iEvent.getByLabel(srcName,theTCCollection );
197  hitAssociator = new TrackerHitAssociator(iEvent);
198 
199  TrajectoryStateCombiner combiner;
200 
201  for (TrackCandidateCollection::const_iterator i=theTCCollection->begin(); i!=theTCCollection->end();i++){
202 
203  LogTrace("TestHits") << "\n*****************new candidate*****************" << std::endl;
204 
205  const TrackCandidate * theTC = &(*i);
206  PTrajectoryStateOnDet state = theTC->trajectoryStateOnDet();
207  const TrackCandidate::range& recHitVec=theTC->recHits();
208 
209  //convert PTrajectoryStateOnDet to TrajectoryStateOnSurface
210 
211 
212  DetId detId(state.detId());
213  TrajectoryStateOnSurface theTSOS=
214  trajectoryStateTransform::transientState(state, &(theG->idToDet(detId)->surface()),theMF.product());
215 
216  if (theTSOS.globalMomentum().eta()>maxeta || theTSOS.globalMomentum().eta()<mineta) continue;
217 
218  //convert the TrackingRecHit vector to a TransientTrackingRecHit vector
219  TransientTrackingRecHit::RecHitContainer hits;
220 
221  for (edm::OwnVector<TrackingRecHit>::const_iterator i=recHitVec.first;
222  i!=recHitVec.second; i++){
223  hits.push_back(theBuilder->build(&(*i) ));
224  }
225 
226  //call the fitter
227  std::vector<Trajectory> result = fit->fit(theTC->seed(), hits, theTSOS);
228  if (result.size()==0) continue;
229  std::vector<TrajectoryMeasurement> vtm = result[0].measurements();
230  double tchi2 = 0;
231 
232  TSOS lastState = theTSOS;
233  for (std::vector<TrajectoryMeasurement>::iterator tm=vtm.begin(); tm!=vtm.end();tm++){
234 
235  TransientTrackingRecHit::ConstRecHitPointer rhit = tm->recHit();
236  if ((rhit)->isValid()==0&&rhit->det()!=0) continue;
237  LogTrace("TestHits") << "*****************new hit*****************" ;
238 
239  int subdetId = rhit->det()->geographicalId().subdetId();
240  int layerId = 0;
241  DetId id = rhit->det()->geographicalId();
242  if (id.subdetId()==3) layerId = ((TIBDetId)(id)).layer();
243  if (id.subdetId()==5) layerId = ((TOBDetId)(id)).layer();
244  if (id.subdetId()==1) layerId = ((PXBDetId)(id)).layer();
245  if (id.subdetId()==4) layerId = ((TIDDetId)(id)).wheel();
246  if (id.subdetId()==6) layerId = ((TECDetId)(id)).wheel();
247  if (id.subdetId()==2) layerId = ((PXFDetId)(id)).disk();
248  LogTrace("TestHits") << "subdetId=" << subdetId << " layerId=" << layerId ;
249 
250  double delta = 99999;
251  LocalPoint rhitLPv = rhit->localPosition();
252 
253  std::vector<PSimHit> assSimHits = hitAssociator->associateHit(*(rhit->hit()));
254  if (assSimHits.size()==0) continue;
255  PSimHit shit;
256  for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){
257  if ((m->localPosition()-rhitLPv).mag()<delta) {
258  shit=*m;
259  delta = (m->localPosition()-rhitLPv).mag();
260  }
261  }
262 
263  TSOS currentState = tm->forwardPredictedState();
264  if (tm->backwardPredictedState().isValid())
265  currentState = combiner(tm->backwardPredictedState(), tm->forwardPredictedState());
266  TSOS updatedState = tm->updatedState();
267  tchi2+=tm->estimate();
268 
269  //plot chi2 increment
270  double chi2increment = tm->estimate();
271  LogTrace("TestHits") << "tm->estimate()=" << tm->estimate();
272  title.str("");
273  title << "Chi2Increment_" << subdetId << "-" << layerId;
274  hChi2Increment[title.str()]->Fill( chi2increment );
275  hTotChi2Increment->Fill( chi2increment );
276  hProcess_vs_Chi2->Fill( chi2increment, shit.processType() );
277  if (dynamic_cast<const SiPixelRecHit*>(rhit->hit()))
278  hClsize_vs_Chi2->Fill( chi2increment, ((const SiPixelRecHit*)(rhit->hit()))->cluster()->size() );
279  if (dynamic_cast<const SiStripRecHit2D*>(rhit->hit()))
280  hClsize_vs_Chi2->Fill( chi2increment, ((const SiStripRecHit2D*)(rhit->hit()))->cluster()->amplitudes().size() );
281 
282  //test hits
283  const Surface * surf = &( (rhit)->det()->surface() );
284  LocalVector shitLMom;
285  LocalPoint shitLPos;
286  if (dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())) {
287  double rechitmatchedx = rhit->localPosition().x();
288  double rechitmatchedy = rhit->localPosition().y();
289  double mindist = 999999;
290  double distx, disty;
291  std::pair<LocalPoint,LocalVector> closestPair;
292  const StripGeomDetUnit* stripDet =(StripGeomDetUnit*) ((const GluedGeomDet *)(rhit)->det())->stereoDet();
293  const BoundPlane& plane = (rhit)->det()->surface();
294  for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++) {
295  //project simhit;
296  std::pair<LocalPoint,LocalVector> hitPair = projectHit((*m),stripDet,plane);
297  distx = fabs(rechitmatchedx - hitPair.first.x());
298  disty = fabs(rechitmatchedy - hitPair.first.y());
299  double dist = distx*distx+disty*disty;
300  if(sqrt(dist)<mindist){
301  mindist = dist;
302  closestPair = hitPair;
303  }
304  }
305  shitLPos = closestPair.first;
306  shitLMom = closestPair.second;
307  } else {
308  shitLPos = shit.localPosition();
309  shitLMom = shit.momentumAtEntry();
310  }
311  GlobalVector shitGMom = surf->toGlobal(shitLMom);
312  GlobalPoint shitGPos = surf->toGlobal(shitLPos);
313 
314  GlobalVector tsosGMom = currentState.globalMomentum();
315  GlobalError tsosGMEr(currentState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3,3));
316  GlobalPoint tsosGPos = currentState.globalPosition();
317  GlobalError tsosGPEr = currentState.cartesianError().position();
318 
319  GlobalPoint rhitGPos = (rhit)->globalPosition();
320  GlobalError rhitGPEr = (rhit)->globalPositionError();
321 
322  double pullGPX_rs = (rhitGPos.x()-shitGPos.x())/sqrt(rhitGPEr.cxx());
323  double pullGPY_rs = (rhitGPos.y()-shitGPos.y())/sqrt(rhitGPEr.cyy());
324  double pullGPZ_rs = (rhitGPos.z()-shitGPos.z())/sqrt(rhitGPEr.czz());
325  //double pullGPX_rs = (rhitGPos.x()-shitGPos.x());
326  //double pullGPY_rs = (rhitGPos.y()-shitGPos.y());
327  //double pullGPZ_rs = (rhitGPos.z()-shitGPos.z());
328 
329  LogTrace("TestHits") << "rs" << std::endl;
330  LogVerbatim("TestHits") << "assSimHits.size()=" << assSimHits.size() ;
331  LogVerbatim("TestHits") << "tsos globalPos =" << tsosGPos ;
332  LogVerbatim("TestHits") << "sim hit globalPos=" << shitGPos ;
333  LogVerbatim("TestHits") << "rec hit globalPos=" << rhitGPos ;
334  LogVerbatim("TestHits") << "geographicalId =" << rhit->det()->geographicalId().rawId() ;
335  LogVerbatim("TestHits") << "surface position =" << surf->position() ;
336 
337  title.str("");
338  title << "PullGP_X_" << subdetId << "-" << layerId << "_rs";
339  hPullGP_X_rs[title.str()]->Fill( pullGPX_rs );
340  title.str("");
341  title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs";
342  hPullGP_Y_rs[title.str()]->Fill( pullGPY_rs );
343  title.str("");
344  title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs";
345  hPullGP_Z_rs[title.str()]->Fill( pullGPZ_rs );
346 
347  double pullGPX_tr = (tsosGPos.x()-rhitGPos.x())/sqrt(tsosGPEr.cxx()+rhitGPEr.cxx());
348  double pullGPY_tr = (tsosGPos.y()-rhitGPos.y())/sqrt(tsosGPEr.cyy()+rhitGPEr.cyy());
349  double pullGPZ_tr = (tsosGPos.z()-rhitGPos.z())/sqrt(tsosGPEr.czz()+rhitGPEr.czz());
350  //double pullGPX_tr = (tsosGPos.x()-rhitGPos.x());
351  //double pullGPY_tr = (tsosGPos.y()-rhitGPos.y());
352  //double pullGPZ_tr = (tsosGPos.z()-rhitGPos.z());
353 
354  LogTrace("TestHits") << "tr" << std::endl;
355 
356  title.str("");
357  title << "PullGP_X_" << subdetId << "-" << layerId << "_tr";
358  hPullGP_X_tr[title.str()]->Fill( pullGPX_tr );
359  title.str("");
360  title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr";
361  hPullGP_Y_tr[title.str()]->Fill( pullGPY_tr );
362  title.str("");
363  title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr";
364  hPullGP_Z_tr[title.str()]->Fill( pullGPZ_tr );
365 
366  double pullGPX_ts = (tsosGPos.x()-shitGPos.x())/sqrt(tsosGPEr.cxx());
367  double pullGPY_ts = (tsosGPos.y()-shitGPos.y())/sqrt(tsosGPEr.cyy());
368  double pullGPZ_ts = (tsosGPos.z()-shitGPos.z())/sqrt(tsosGPEr.czz());
369  //double pullGPX_ts = (tsosGPos.x()-shitGPos.x());
370  //double pullGPY_ts = (tsosGPos.y()-shitGPos.y());
371  //double pullGPZ_ts = (tsosGPos.z()-shitGPos.z());
372 
373  LogTrace("TestHits") << "ts1" << std::endl;
374 
375  title.str("");
376  title << "PullGP_X_" << subdetId << "-" << layerId << "_ts";
377  hPullGP_X_ts[title.str()]->Fill( pullGPX_ts );
378  title.str("");
379  title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts";
380  hPullGP_Y_ts[title.str()]->Fill( pullGPY_ts );
381  title.str("");
382  title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts";
383  hPullGP_Z_ts[title.str()]->Fill( pullGPZ_ts );
384 
385  double pullGMX_ts = (tsosGMom.x()-shitGMom.x())/sqrt(tsosGMEr.cxx());
386  double pullGMY_ts = (tsosGMom.y()-shitGMom.y())/sqrt(tsosGMEr.cyy());
387  double pullGMZ_ts = (tsosGMom.z()-shitGMom.z())/sqrt(tsosGMEr.czz());
388  //double pullGMX_ts = (tsosGMom.x()-shitGMom.x());
389  //double pullGMY_ts = (tsosGMom.y()-shitGMom.y());
390  //double pullGMZ_ts = (tsosGMom.z()-shitGMom.z());
391 
392  LogTrace("TestHits") << "ts2" << std::endl;
393 
394  title.str("");
395  title << "PullGM_X_" << subdetId << "-" << layerId << "_ts";
396  hPullGM_X_ts[title.str()]->Fill( pullGMX_ts );
397  title.str("");
398  title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts";
399  hPullGM_Y_ts[title.str()]->Fill( pullGMY_ts );
400  title.str("");
401  title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts";
402  hPullGM_Z_ts[title.str()]->Fill( pullGMZ_ts );
403 
404  if (dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())) {
405  //mono
406  LogTrace("TestHits") << "MONO HIT" << std::endl;
407  auto m = dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())->monoHit();
408  CTTRHp tMonoHit = theBuilder->build(&m);
409  if (tMonoHit==0) continue;
410  vector<PSimHit> assMonoSimHits = hitAssociator->associateHit(*tMonoHit->hit());
411  if (assMonoSimHits.size()==0) continue;
412  const PSimHit sMonoHit = *(assSimHits.begin());
413  const Surface * monoSurf = &( tMonoHit->det()->surface() );
414  if (monoSurf==0) continue;
415  TSOS monoState = thePropagator->propagate(lastState,*monoSurf);
416  if (monoState.isValid()==0) continue;
417 
418  LocalVector monoShitLMom = sMonoHit.momentumAtEntry();
419  GlobalVector monoShitGMom = monoSurf->toGlobal(monoShitLMom);
420  LocalPoint monoShitLPos = sMonoHit.localPosition();
421  GlobalPoint monoShitGPos = monoSurf->toGlobal(monoShitLPos);
422 
423  GlobalVector monoTsosGMom = monoState.globalMomentum();
424  GlobalError monoTsosGMEr(monoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3,3));
425  GlobalPoint monoTsosGPos = monoState.globalPosition();
426  GlobalError monoTsosGPEr = monoState.cartesianError().position();
427 
428  GlobalPoint monoRhitGPos = tMonoHit->globalPosition();
429  GlobalError monoRhitGPEr = tMonoHit->globalPositionError();
430 
431  double pullGPX_rs_mono = (monoRhitGPos.x()-monoShitGPos.x())/sqrt(monoRhitGPEr.cxx());
432  double pullGPY_rs_mono = (monoRhitGPos.y()-monoShitGPos.y())/sqrt(monoRhitGPEr.cyy());
433  double pullGPZ_rs_mono = (monoRhitGPos.z()-monoShitGPos.z())/sqrt(monoRhitGPEr.czz());
434  //double pullGPX_rs_mono = (monoRhitGPos.x()-monoShitGPos.x());
435  //double pullGPY_rs_mono = (monoRhitGPos.y()-monoShitGPos.y());
436  //double pullGPZ_rs_mono = (monoRhitGPos.z()-monoShitGPos.z());
437 
438  title.str("");
439  title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_mono";
440  hPullGP_X_rs_mono[title.str()]->Fill( pullGPX_rs_mono );
441  title.str("");
442  title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_mono";
443  hPullGP_Y_rs_mono[title.str()]->Fill( pullGPY_rs_mono );
444  title.str("");
445  title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_mono";
446  hPullGP_Z_rs_mono[title.str()]->Fill( pullGPZ_rs_mono );
447 
448  double pullGPX_tr_mono = (monoTsosGPos.x()-monoRhitGPos.x())/sqrt(monoTsosGPEr.cxx()+monoRhitGPEr.cxx());
449  double pullGPY_tr_mono = (monoTsosGPos.y()-monoRhitGPos.y())/sqrt(monoTsosGPEr.cyy()+monoRhitGPEr.cyy());
450  double pullGPZ_tr_mono = (monoTsosGPos.z()-monoRhitGPos.z())/sqrt(monoTsosGPEr.czz()+monoRhitGPEr.czz());
451  //double pullGPX_tr_mono = (monoTsosGPos.x()-monoRhitGPos.x());
452  //double pullGPY_tr_mono = (monoTsosGPos.y()-monoRhitGPos.y());
453  //double pullGPZ_tr_mono = (monoTsosGPos.z()-monoRhitGPos.z());
454 
455  title.str("");
456  title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_mono";
457  hPullGP_X_tr_mono[title.str()]->Fill( pullGPX_tr_mono );
458  title.str("");
459  title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_mono";
460  hPullGP_Y_tr_mono[title.str()]->Fill( pullGPY_tr_mono );
461  title.str("");
462  title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_mono";
463  hPullGP_Z_tr_mono[title.str()]->Fill( pullGPZ_tr_mono );
464 
465  double pullGPX_ts_mono = (monoTsosGPos.x()-monoShitGPos.x())/sqrt(monoTsosGPEr.cxx());
466  double pullGPY_ts_mono = (monoTsosGPos.y()-monoShitGPos.y())/sqrt(monoTsosGPEr.cyy());
467  double pullGPZ_ts_mono = (monoTsosGPos.z()-monoShitGPos.z())/sqrt(monoTsosGPEr.czz());
468  //double pullGPX_ts_mono = (monoTsosGPos.x()-monoShitGPos.x());
469  //double pullGPY_ts_mono = (monoTsosGPos.y()-monoShitGPos.y());
470  //double pullGPZ_ts_mono = (monoTsosGPos.z()-monoShitGPos.z());
471 
472  title.str("");
473  title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_mono";
474  hPullGP_X_ts_mono[title.str()]->Fill( pullGPX_ts_mono );
475  title.str("");
476  title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_mono";
477  hPullGP_Y_ts_mono[title.str()]->Fill( pullGPY_ts_mono );
478  title.str("");
479  title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_mono";
480  hPullGP_Z_ts_mono[title.str()]->Fill( pullGPZ_ts_mono );
481 
482  double pullGMX_ts_mono = (monoTsosGMom.x()-monoShitGMom.x())/sqrt(monoTsosGMEr.cxx());
483  double pullGMY_ts_mono = (monoTsosGMom.y()-monoShitGMom.y())/sqrt(monoTsosGMEr.cyy());
484  double pullGMZ_ts_mono = (monoTsosGMom.z()-monoShitGMom.z())/sqrt(monoTsosGMEr.czz());
485  //double pullGMX_ts_mono = (monoTsosGMom.x()-monoShitGMom.x());
486  //double pullGMY_ts_mono = (monoTsosGMom.y()-monoShitGMom.y());
487  //double pullGMZ_ts_mono = (monoTsosGMom.z()-monoShitGMom.z());
488 
489  title.str("");
490  title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_mono";
491  hPullGM_X_ts_mono[title.str()]->Fill( pullGMX_ts_mono );
492  title.str("");
493  title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_mono";
494  hPullGM_Y_ts_mono[title.str()]->Fill( pullGMY_ts_mono );
495  title.str("");
496  title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_mono";
497  hPullGM_Z_ts_mono[title.str()]->Fill( pullGMZ_ts_mono );
498 
499  //stereo
500  LogTrace("TestHits") << "STEREO HIT" << std::endl;
501  auto s = dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())->stereoHit();
502  CTTRHp tStereoHit =
503  theBuilder->build(&s);
504  if (tStereoHit==0) continue;
505  vector<PSimHit> assStereoSimHits = hitAssociator->associateHit(*tStereoHit->hit());
506  if (assStereoSimHits.size()==0) continue;
507  const PSimHit sStereoHit = *(assSimHits.begin());
508  const Surface * stereoSurf = &( tStereoHit->det()->surface() );
509  if (stereoSurf==0) continue;
510  TSOS stereoState = thePropagator->propagate(lastState,*stereoSurf);
511  if (stereoState.isValid()==0) continue;
512 
513  LocalVector stereoShitLMom = sStereoHit.momentumAtEntry();
514  GlobalVector stereoShitGMom = stereoSurf->toGlobal(stereoShitLMom);
515  LocalPoint stereoShitLPos = sStereoHit.localPosition();
516  GlobalPoint stereoShitGPos = stereoSurf->toGlobal(stereoShitLPos);
517 
518  GlobalVector stereoTsosGMom = stereoState.globalMomentum();
519  GlobalError stereoTsosGMEr(stereoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3,3));
520  GlobalPoint stereoTsosGPos = stereoState.globalPosition();
521  GlobalError stereoTsosGPEr = stereoState.cartesianError().position();
522 
523  GlobalPoint stereoRhitGPos = tStereoHit->globalPosition();
524  GlobalError stereoRhitGPEr = tStereoHit->globalPositionError();
525 
526  double pullGPX_rs_stereo = (stereoRhitGPos.x()-stereoShitGPos.x())/sqrt(stereoRhitGPEr.cxx());
527  double pullGPY_rs_stereo = (stereoRhitGPos.y()-stereoShitGPos.y())/sqrt(stereoRhitGPEr.cyy());
528  double pullGPZ_rs_stereo = (stereoRhitGPos.z()-stereoShitGPos.z())/sqrt(stereoRhitGPEr.czz());
529  //double pullGPX_rs_stereo = (stereoRhitGPos.x()-stereoShitGPos.x());
530  //double pullGPY_rs_stereo = (stereoRhitGPos.y()-stereoShitGPos.y());
531  //double pullGPZ_rs_stereo = (stereoRhitGPos.z()-stereoShitGPos.z());
532 
533  title.str("");
534  title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_stereo";
535  hPullGP_X_rs_stereo[title.str()]->Fill( pullGPX_rs_stereo );
536  title.str("");
537  title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_stereo";
538  hPullGP_Y_rs_stereo[title.str()]->Fill( pullGPY_rs_stereo );
539  title.str("");
540  title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_stereo";
541  hPullGP_Z_rs_stereo[title.str()]->Fill( pullGPZ_rs_stereo );
542 
543  double pullGPX_tr_stereo = (stereoTsosGPos.x()-stereoRhitGPos.x())/sqrt(stereoTsosGPEr.cxx()+stereoRhitGPEr.cxx());
544  double pullGPY_tr_stereo = (stereoTsosGPos.y()-stereoRhitGPos.y())/sqrt(stereoTsosGPEr.cyy()+stereoRhitGPEr.cyy());
545  double pullGPZ_tr_stereo = (stereoTsosGPos.z()-stereoRhitGPos.z())/sqrt(stereoTsosGPEr.czz()+stereoRhitGPEr.czz());
546  //double pullGPX_tr_stereo = (stereoTsosGPos.x()-stereoRhitGPos.x());
547  //double pullGPY_tr_stereo = (stereoTsosGPos.y()-stereoRhitGPos.y());
548  //double pullGPZ_tr_stereo = (stereoTsosGPos.z()-stereoRhitGPos.z());
549 
550  title.str("");
551  title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_stereo";
552  hPullGP_X_tr_stereo[title.str()]->Fill( pullGPX_tr_stereo );
553  title.str("");
554  title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_stereo";
555  hPullGP_Y_tr_stereo[title.str()]->Fill( pullGPY_tr_stereo );
556  title.str("");
557  title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_stereo";
558  hPullGP_Z_tr_stereo[title.str()]->Fill( pullGPZ_tr_stereo );
559 
560  double pullGPX_ts_stereo = (stereoTsosGPos.x()-stereoShitGPos.x())/sqrt(stereoTsosGPEr.cxx());
561  double pullGPY_ts_stereo = (stereoTsosGPos.y()-stereoShitGPos.y())/sqrt(stereoTsosGPEr.cyy());
562  double pullGPZ_ts_stereo = (stereoTsosGPos.z()-stereoShitGPos.z())/sqrt(stereoTsosGPEr.czz());
563  //double pullGPX_ts_stereo = (stereoTsosGPos.x()-stereoShitGPos.x());
564  //double pullGPY_ts_stereo = (stereoTsosGPos.y()-stereoShitGPos.y());
565  //double pullGPZ_ts_stereo = (stereoTsosGPos.z()-stereoShitGPos.z());
566 
567  title.str("");
568  title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_stereo";
569  hPullGP_X_ts_stereo[title.str()]->Fill( pullGPX_ts_stereo );
570  title.str("");
571  title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_stereo";
572  hPullGP_Y_ts_stereo[title.str()]->Fill( pullGPY_ts_stereo );
573  title.str("");
574  title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_stereo";
575  hPullGP_Z_ts_stereo[title.str()]->Fill( pullGPZ_ts_stereo );
576 
577  double pullGMX_ts_stereo = (stereoTsosGMom.x()-stereoShitGMom.x())/sqrt(stereoTsosGMEr.cxx());
578  double pullGMY_ts_stereo = (stereoTsosGMom.y()-stereoShitGMom.y())/sqrt(stereoTsosGMEr.cyy());
579  double pullGMZ_ts_stereo = (stereoTsosGMom.z()-stereoShitGMom.z())/sqrt(stereoTsosGMEr.czz());
580  //double pullGMX_ts_stereo = (stereoTsosGMom.x()-stereoShitGMom.x());
581  //double pullGMY_ts_stereo = (stereoTsosGMom.y()-stereoShitGMom.y());
582  //double pullGMZ_ts_stereo = (stereoTsosGMom.z()-stereoShitGMom.z());
583 
584  title.str("");
585  title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_stereo";
586  hPullGM_X_ts_stereo[title.str()]->Fill( pullGMX_ts_stereo );
587  title.str("");
588  title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_stereo";
589  hPullGM_Y_ts_stereo[title.str()]->Fill( pullGMY_ts_stereo );
590  title.str("");
591  title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_stereo";
592  hPullGM_Z_ts_stereo[title.str()]->Fill( pullGMZ_ts_stereo );
593  }
594  lastState = updatedState;
595  }
596  LogTrace("TestHits") << "traj chi2=" << tchi2 ;
597  LogTrace("TestHits") << "track chi2=" << result[0].chiSquared() ;
598  }
599  delete hitAssociator;
600  LogTrace("TestHits") << "end of event" << std::endl;
601 }
602 // TSOS lastState = theTSOS;
603 // for (std::vector<TrajectoryMeasurement>::iterator tm=vtm.begin(); tm!=vtm.end();tm++){
604 
605 // TransientTrackingRecHit::ConstRecHitPointer rhit = tm->recHit();
606 // if ((rhit)->isValid()==0&&rhit->det()!=0) continue;
607 
608 // // //test hits
609 // // TSOS lastState, currentState, updatedState;
610 // // updatedState=theTSOS;
611 // // for (TransientTrackingRecHit::RecHitContainer::iterator rhit=hits.begin()+1;rhit!=hits.end();++rhit){
612 
613 // lastState=updatedState;
614 
615 // LogTrace("TestHits") << "new hit" << std::endl;
616 
617 // if ((*rhit)->isValid()==0) continue;
618 
619 // int subdetId = (*rhit)->det()->geographicalId().subdetId();
620 // int layerId = 0;
621 // DetId id = (*rhit)->det()->geographicalId();
622 // if (id.subdetId()==3) layerId = ((TIBDetId)(id)).layer();
623 // if (id.subdetId()==5) layerId = ((TOBDetId)(id)).layer();
624 // if (id.subdetId()==1) layerId = ((PXBDetId)(id)).layer();
625 // if (id.subdetId()==4) layerId = ((TIDDetId)(id)).wheel();
626 // if (id.subdetId()==6) layerId = ((TECDetId)(id)).wheel();
627 // if (id.subdetId()==2) layerId = ((PXFDetId)(id)).disk();
628 // const Surface * surf = &( (*rhit)->det()->surface() );
629 // currentState=thePropagator->propagate(lastState,*surf);
630 // if (currentState.isValid()==0) continue;
631 // updatedState=theUpdator->update(currentState,**rhit);
632 
633 // double delta = 99999;
634 // LocalPoint rhitLP = rhit->localPosition();
635 
636 // std::vector<PSimHit> assSimHits = hitAssociator->associateHit(*(*rhit)->hit());
637 // if (assSimHits.size()==0) continue;
638 // PSimHit shit;
639 // for(std::vector<PSimHit>::const_iterator m=assSimHits.begin(); m<assSimHits.end(); m++){
640 // if ((*m-rhitLP).mag()<delta) {
641 // shit=*m;
642 // delta = (*m-rhitLP).mag();
643 // }
644 // }
645 
646 // LocalVector shitLMom = shit.momentumAtEntry();
647 // GlobalVector shitGMom = surf->toGlobal(shitLMom);
648 // LocalPoint shitLPos = shit.localPosition();
649 // GlobalPoint shitGPos = surf->toGlobal(shitLPos);
650 
651 // GlobalVector tsosGMom = currentState.globalMomentum();
652 // GlobalError tsosGMEr(currentState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3,3));
653 // GlobalPoint tsosGPos = currentState.globalPosition();
654 // GlobalError tsosGPEr = currentState.cartesianError().position();
655 
656 // GlobalPoint rhitGPos = (*rhit)->globalPosition();
657 // GlobalError rhitGPEr = (*rhit)->globalPositionError();
658 
659 // double pullGPX_rs = (rhitGPos.x()-shitGPos.x())/sqrt(rhitGPEr.cxx());
660 // double pullGPY_rs = (rhitGPos.y()-shitGPos.y())/sqrt(rhitGPEr.cyy());
661 // double pullGPZ_rs = (rhitGPos.z()-shitGPos.z())/sqrt(rhitGPEr.czz());
662 // // double pullGPX_rs = (rhitGPos.x()-shitGPos.x())/*/sqrt(rhitGPEr.cxx())*/;
663 // // double pullGPY_rs = (rhitGPos.y()-shitGPos.y())/*/sqrt(rhitGPEr.cyy())*/;
664 // // double pullGPZ_rs = (rhitGPos.z()-shitGPos.z())/*/sqrt(rhitGPEr.czz())*/;
665 
666 // //plot chi2 increment
667 // MeasurementExtractor me(currentState);
668 // double chi2increment = computeChi2Increment(me,*rhit);
669 // LogTrace("TestHits") << "chi2increment=" << chi2increment << std::endl;
670 // title.str("");
671 // title << "Chi2Increment_" << subdetId << "-" << layerId;
672 // hChi2Increment[title.str()]->Fill( chi2increment );
673 // hTotChi2Increment->Fill( chi2increment );
674 // hChi2_vs_Process->Fill( chi2increment, shit.processType() );
675 // if (dynamic_cast<const SiPixelRecHit*>((*rhit)->hit()))
676 // hChi2_vs_clsize->Fill( chi2increment, ((const SiPixelRecHit*)(*rhit)->hit())->cluster()->size() );
677 // if (dynamic_cast<const SiStripRecHit2D*>((*rhit)->hit()))
678 // hChi2_vs_clsize->Fill( chi2increment, ((const SiStripRecHit2D*)(*rhit)->hit())->cluster()->amplitudes().size() );
679 
680 // LogTrace("TestHits") << "rs" << std::endl;
681 
682 // title.str("");
683 // title << "PullGP_X_" << subdetId << "-" << layerId << "_rs";
684 // hPullGP_X_rs[title.str()]->Fill( pullGPX_rs );
685 // title.str("");
686 // title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs";
687 // hPullGP_Y_rs[title.str()]->Fill( pullGPY_rs );
688 // title.str("");
689 // title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs";
690 // hPullGP_Z_rs[title.str()]->Fill( pullGPZ_rs );
691 
692 // double pullGPX_tr = (tsosGPos.x()-rhitGPos.x())/sqrt(tsosGPEr.cxx()+rhitGPEr.cxx());
693 // double pullGPY_tr = (tsosGPos.y()-rhitGPos.y())/sqrt(tsosGPEr.cyy()+rhitGPEr.cyy());
694 // double pullGPZ_tr = (tsosGPos.z()-rhitGPos.z())/sqrt(tsosGPEr.czz()+rhitGPEr.czz());
695 // // double pullGPX_tr = (tsosGPos.x()-rhitGPos.x())/*/sqrt(tsosGPEr.cxx()+rhitGPEr.cxx())*/;
696 // // double pullGPY_tr = (tsosGPos.y()-rhitGPos.y())/*/sqrt(tsosGPEr.cyy()+rhitGPEr.cyy())*/;
697 // // double pullGPZ_tr = (tsosGPos.z()-rhitGPos.z())/*/sqrt(tsosGPEr.czz()+rhitGPEr.czz())*/;
698 
699 // LogTrace("TestHits") << "tr" << std::endl;
700 
701 // title.str("");
702 // title << "PullGP_X_" << subdetId << "-" << layerId << "_tr";
703 // hPullGP_X_tr[title.str()]->Fill( pullGPX_tr );
704 // title.str("");
705 // title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr";
706 // hPullGP_Y_tr[title.str()]->Fill( pullGPY_tr );
707 // title.str("");
708 // title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr";
709 // hPullGP_Z_tr[title.str()]->Fill( pullGPZ_tr );
710 
711 // double pullGPX_ts = (tsosGPos.x()-shitGPos.x())/sqrt(tsosGPEr.cxx());
712 // double pullGPY_ts = (tsosGPos.y()-shitGPos.y())/sqrt(tsosGPEr.cyy());
713 // double pullGPZ_ts = (tsosGPos.z()-shitGPos.z())/sqrt(tsosGPEr.czz());
714 // // double pullGPX_ts = (tsosGPos.x()-shitGPos.x())/*/sqrt(tsosGPEr.cxx())*/;
715 // // double pullGPY_ts = (tsosGPos.y()-shitGPos.y())/*/sqrt(tsosGPEr.cyy())*/;
716 // // double pullGPZ_ts = (tsosGPos.z()-shitGPos.z())/*/sqrt(tsosGPEr.czz())*/;
717 
718 // LogTrace("TestHits") << "ts1" << std::endl;
719 
720 // title.str("");
721 // title << "PullGP_X_" << subdetId << "-" << layerId << "_ts";
722 // hPullGP_X_ts[title.str()]->Fill( pullGPX_ts );
723 // title.str("");
724 // title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts";
725 // hPullGP_Y_ts[title.str()]->Fill( pullGPY_ts );
726 // title.str("");
727 // title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts";
728 // hPullGP_Z_ts[title.str()]->Fill( pullGPZ_ts );
729 
730 // double pullGMX_ts = (tsosGMom.x()-shitGMom.x())/sqrt(tsosGMEr.cxx());
731 // double pullGMY_ts = (tsosGMom.y()-shitGMom.y())/sqrt(tsosGMEr.cyy());
732 // double pullGMZ_ts = (tsosGMom.z()-shitGMom.z())/sqrt(tsosGMEr.czz());
733 // // double pullGMX_ts = (tsosGMom.x()-shitGMom.x())/*/sqrt(tsosGMEr.cxx())*/;
734 // // double pullGMY_ts = (tsosGMom.y()-shitGMom.y())/*/sqrt(tsosGMEr.cyy())*/;
735 // // double pullGMZ_ts = (tsosGMom.z()-shitGMom.z())/*/sqrt(tsosGMEr.czz())*/;
736 
737 // LogTrace("TestHits") << "ts2" << std::endl;
738 
739 // title.str("");
740 // title << "PullGM_X_" << subdetId << "-" << layerId << "_ts";
741 // hPullGM_X_ts[title.str()]->Fill( pullGMX_ts );
742 // title.str("");
743 // title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts";
744 // hPullGM_Y_ts[title.str()]->Fill( pullGMY_ts );
745 // title.str("");
746 // title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts";
747 // hPullGM_Z_ts[title.str()]->Fill( pullGMZ_ts );
748 
749 // if (dynamic_cast<const SiStripMatchedRecHit2D*>((*rhit)->hit())) {
750 // //mono
751 // LogTrace("TestHits") << "MONO HIT" << std::endl;
752 // CTTRHp tMonoHit =
753 // theBuilder->build(dynamic_cast<const SiStripMatchedRecHit2D*>((*rhit)->hit())->monoHit());
754 // if (tMonoHit==0) continue;
755 // vector<PSimHit> assMonoSimHits = hitAssociator->associateHit(*tMonoHit->hit());
756 // if (assMonoSimHits.size()==0) continue;
757 // const PSimHit sMonoHit = *(assSimHits.begin());
758 // const Surface * monoSurf = &( tMonoHit->det()->surface() );
759 // if (monoSurf==0) continue;
760 // TSOS monoState = thePropagator->propagate(lastState,*monoSurf);
761 // if (monoState.isValid()==0) continue;
762 
763 // LocalVector monoShitLMom = sMonoHit.momentumAtEntry();
764 // GlobalVector monoShitGMom = monoSurf->toGlobal(monoShitLMom);
765 // LocalPoint monoShitLPos = sMonoHit.localPosition();
766 // GlobalPoint monoShitGPos = monoSurf->toGlobal(monoShitLPos);
767 
768 // GlobalVector monoTsosGMom = monoState.globalMomentum();
769 // GlobalError monoTsosGMEr(monoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3,3));
770 // GlobalPoint monoTsosGPos = monoState.globalPosition();
771 // GlobalError monoTsosGPEr = monoState.cartesianError().position();
772 
773 // GlobalPoint monoRhitGPos = tMonoHit->globalPosition();
774 // GlobalError monoRhitGPEr = tMonoHit->globalPositionError();
775 
776 // double pullGPX_rs_mono = (monoRhitGPos.x()-monoShitGPos.x())/sqrt(monoRhitGPEr.cxx());
777 // double pullGPY_rs_mono = (monoRhitGPos.y()-monoShitGPos.y())/sqrt(monoRhitGPEr.cyy());
778 // double pullGPZ_rs_mono = (monoRhitGPos.z()-monoShitGPos.z())/sqrt(monoRhitGPEr.czz());
779 // // double pullGPX_rs_mono = (monoRhitGPos.x()-monoShitGPos.x())/*/sqrt(monoRhitGPEr.cxx())*/;
780 // // double pullGPY_rs_mono = (monoRhitGPos.y()-monoShitGPos.y())/*/sqrt(monoRhitGPEr.cyy())*/;
781 // // double pullGPZ_rs_mono = (monoRhitGPos.z()-monoShitGPos.z())/*/sqrt(monoRhitGPEr.czz())*/;
782 
783 // title.str("");
784 // title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_mono";
785 // hPullGP_X_rs_mono[title.str()]->Fill( pullGPX_rs_mono );
786 // title.str("");
787 // title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_mono";
788 // hPullGP_Y_rs_mono[title.str()]->Fill( pullGPY_rs_mono );
789 // title.str("");
790 // title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_mono";
791 // hPullGP_Z_rs_mono[title.str()]->Fill( pullGPZ_rs_mono );
792 
793 // double pullGPX_tr_mono = (monoTsosGPos.x()-monoRhitGPos.x())/sqrt(monoTsosGPEr.cxx()+monoRhitGPEr.cxx());
794 // double pullGPY_tr_mono = (monoTsosGPos.y()-monoRhitGPos.y())/sqrt(monoTsosGPEr.cyy()+monoRhitGPEr.cyy());
795 // double pullGPZ_tr_mono = (monoTsosGPos.z()-monoRhitGPos.z())/sqrt(monoTsosGPEr.czz()+monoRhitGPEr.czz());
796 // // double pullGPX_tr_mono = (monoTsosGPos.x()-monoRhitGPos.x())/*/sqrt(monoTsosGPEr.cxx()+monoRhitGPEr.cxx())*/;
797 // // double pullGPY_tr_mono = (monoTsosGPos.y()-monoRhitGPos.y())/*/sqrt(monoTsosGPEr.cyy()+monoRhitGPEr.cyy())*/;
798 // // double pullGPZ_tr_mono = (monoTsosGPos.z()-monoRhitGPos.z())/*/sqrt(monoTsosGPEr.czz()+monoRhitGPEr.czz())*/;
799 
800 // title.str("");
801 // title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_mono";
802 // hPullGP_X_tr_mono[title.str()]->Fill( pullGPX_tr_mono );
803 // title.str("");
804 // title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_mono";
805 // hPullGP_Y_tr_mono[title.str()]->Fill( pullGPY_tr_mono );
806 // title.str("");
807 // title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_mono";
808 // hPullGP_Z_tr_mono[title.str()]->Fill( pullGPZ_tr_mono );
809 
810 // double pullGPX_ts_mono = (monoTsosGPos.x()-monoShitGPos.x())/sqrt(monoTsosGPEr.cxx());
811 // double pullGPY_ts_mono = (monoTsosGPos.y()-monoShitGPos.y())/sqrt(monoTsosGPEr.cyy());
812 // double pullGPZ_ts_mono = (monoTsosGPos.z()-monoShitGPos.z())/sqrt(monoTsosGPEr.czz());
813 // // double pullGPX_ts_mono = (monoTsosGPos.x()-monoShitGPos.x())/*/sqrt(monoTsosGPEr.cxx())*/;
814 // // double pullGPY_ts_mono = (monoTsosGPos.y()-monoShitGPos.y())/*/sqrt(monoTsosGPEr.cyy())*/;
815 // // double pullGPZ_ts_mono = (monoTsosGPos.z()-monoShitGPos.z())/*/sqrt(monoTsosGPEr.czz())*/;
816 
817 // title.str("");
818 // title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_mono";
819 // hPullGP_X_ts_mono[title.str()]->Fill( pullGPX_ts_mono );
820 // title.str("");
821 // title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_mono";
822 // hPullGP_Y_ts_mono[title.str()]->Fill( pullGPY_ts_mono );
823 // title.str("");
824 // title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_mono";
825 // hPullGP_Z_ts_mono[title.str()]->Fill( pullGPZ_ts_mono );
826 
827 // double pullGMX_ts_mono = (monoTsosGMom.x()-monoShitGMom.x())/sqrt(monoTsosGMEr.cxx());
828 // double pullGMY_ts_mono = (monoTsosGMom.y()-monoShitGMom.y())/sqrt(monoTsosGMEr.cyy());
829 // double pullGMZ_ts_mono = (monoTsosGMom.z()-monoShitGMom.z())/sqrt(monoTsosGMEr.czz());
830 // // double pullGMX_ts_mono = (monoTsosGMom.x()-monoShitGMom.x())/*/sqrt(monoTsosGMEr.cxx())*/;
831 // // double pullGMY_ts_mono = (monoTsosGMom.y()-monoShitGMom.y())/*/sqrt(monoTsosGMEr.cyy())*/;
832 // // double pullGMZ_ts_mono = (monoTsosGMom.z()-monoShitGMom.z())/*/sqrt(monoTsosGMEr.czz())*/;
833 
834 // title.str("");
835 // title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_mono";
836 // hPullGM_X_ts_mono[title.str()]->Fill( pullGMX_ts_mono );
837 // title.str("");
838 // title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_mono";
839 // hPullGM_Y_ts_mono[title.str()]->Fill( pullGMY_ts_mono );
840 // title.str("");
841 // title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_mono";
842 // hPullGM_Z_ts_mono[title.str()]->Fill( pullGMZ_ts_mono );
843 
844 // //stereo
845 // LogTrace("TestHits") << "STEREO HIT" << std::endl;
846 // CTTRHp tStereoHit =
847 // theBuilder->build(dynamic_cast<const SiStripMatchedRecHit2D*>((*rhit)->hit())->stereoHit());
848 // if (tStereoHit==0) continue;
849 // vector<PSimHit> assStereoSimHits = hitAssociator->associateHit(*tStereoHit->hit());
850 // if (assStereoSimHits.size()==0) continue;
851 // const PSimHit sStereoHit = *(assSimHits.begin());
852 // const Surface * stereoSurf = &( tStereoHit->det()->surface() );
853 // if (stereoSurf==0) continue;
854 // TSOS stereoState = thePropagator->propagate(lastState,*stereoSurf);
855 // if (stereoState.isValid()==0) continue;
856 
857 // LocalVector stereoShitLMom = sStereoHit.momentumAtEntry();
858 // GlobalVector stereoShitGMom = stereoSurf->toGlobal(stereoShitLMom);
859 // LocalPoint stereoShitLPos = sStereoHit.localPosition();
860 // GlobalPoint stereoShitGPos = stereoSurf->toGlobal(stereoShitLPos);
861 
862 // GlobalVector stereoTsosGMom = stereoState.globalMomentum();
863 // GlobalError stereoTsosGMEr(stereoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3,3));
864 // GlobalPoint stereoTsosGPos = stereoState.globalPosition();
865 // GlobalError stereoTsosGPEr = stereoState.cartesianError().position();
866 
867 // GlobalPoint stereoRhitGPos = tStereoHit->globalPosition();
868 // GlobalError stereoRhitGPEr = tStereoHit->globalPositionError();
869 
870 // double pullGPX_rs_stereo = (stereoRhitGPos.x()-stereoShitGPos.x())/sqrt(stereoRhitGPEr.cxx());
871 // double pullGPY_rs_stereo = (stereoRhitGPos.y()-stereoShitGPos.y())/sqrt(stereoRhitGPEr.cyy());
872 // double pullGPZ_rs_stereo = (stereoRhitGPos.z()-stereoShitGPos.z())/sqrt(stereoRhitGPEr.czz());
873 // // double pullGPX_rs_stereo = (stereoRhitGPos.x()-stereoShitGPos.x())/*/sqrt(stereoRhitGPEr.cxx())*/;
874 // // double pullGPY_rs_stereo = (stereoRhitGPos.y()-stereoShitGPos.y())/*/sqrt(stereoRhitGPEr.cyy())*/;
875 // // double pullGPZ_rs_stereo = (stereoRhitGPos.z()-stereoShitGPos.z())/*/sqrt(stereoRhitGPEr.czz())*/;
876 
877 // title.str("");
878 // title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_stereo";
879 // hPullGP_X_rs_stereo[title.str()]->Fill( pullGPX_rs_stereo );
880 // title.str("");
881 // title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_stereo";
882 // hPullGP_Y_rs_stereo[title.str()]->Fill( pullGPY_rs_stereo );
883 // title.str("");
884 // title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_stereo";
885 // hPullGP_Z_rs_stereo[title.str()]->Fill( pullGPZ_rs_stereo );
886 
887 // double pullGPX_tr_stereo = (stereoTsosGPos.x()-stereoRhitGPos.x())/sqrt(stereoTsosGPEr.cxx()+stereoRhitGPEr.cxx());
888 // double pullGPY_tr_stereo = (stereoTsosGPos.y()-stereoRhitGPos.y())/sqrt(stereoTsosGPEr.cyy()+stereoRhitGPEr.cyy());
889 // double pullGPZ_tr_stereo = (stereoTsosGPos.z()-stereoRhitGPos.z())/sqrt(stereoTsosGPEr.czz()+stereoRhitGPEr.czz());
890 // // double pullGPX_tr_stereo = (stereoTsosGPos.x()-stereoRhitGPos.x())/*/sqrt(stereoTsosGPEr.cxx()+stereoRhitGPEr.cxx())*/;
891 // // double pullGPY_tr_stereo = (stereoTsosGPos.y()-stereoRhitGPos.y())/*/sqrt(stereoTsosGPEr.cyy()+stereoRhitGPEr.cyy())*/;
892 // // double pullGPZ_tr_stereo = (stereoTsosGPos.z()-stereoRhitGPos.z())/*/sqrt(stereoTsosGPEr.czz()+stereoRhitGPEr.czz())*/;
893 
894 // title.str("");
895 // title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_stereo";
896 // hPullGP_X_tr_stereo[title.str()]->Fill( pullGPX_tr_stereo );
897 // title.str("");
898 // title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_stereo";
899 // hPullGP_Y_tr_stereo[title.str()]->Fill( pullGPY_tr_stereo );
900 // title.str("");
901 // title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_stereo";
902 // hPullGP_Z_tr_stereo[title.str()]->Fill( pullGPZ_tr_stereo );
903 
904 // double pullGPX_ts_stereo = (stereoTsosGPos.x()-stereoShitGPos.x())/sqrt(stereoTsosGPEr.cxx());
905 // double pullGPY_ts_stereo = (stereoTsosGPos.y()-stereoShitGPos.y())/sqrt(stereoTsosGPEr.cyy());
906 // double pullGPZ_ts_stereo = (stereoTsosGPos.z()-stereoShitGPos.z())/sqrt(stereoTsosGPEr.czz());
907 // // double pullGPX_ts_stereo = (stereoTsosGPos.x()-stereoShitGPos.x())/*/sqrt(stereoTsosGPEr.cxx())*/;
908 // // double pullGPY_ts_stereo = (stereoTsosGPos.y()-stereoShitGPos.y())/*/sqrt(stereoTsosGPEr.cyy())*/;
909 // // double pullGPZ_ts_stereo = (stereoTsosGPos.z()-stereoShitGPos.z())/*/sqrt(stereoTsosGPEr.czz())*/;
910 
911 // title.str("");
912 // title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_stereo";
913 // hPullGP_X_ts_stereo[title.str()]->Fill( pullGPX_ts_stereo );
914 // title.str("");
915 // title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_stereo";
916 // hPullGP_Y_ts_stereo[title.str()]->Fill( pullGPY_ts_stereo );
917 // title.str("");
918 // title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_stereo";
919 // hPullGP_Z_ts_stereo[title.str()]->Fill( pullGPZ_ts_stereo );
920 
921 // double pullGMX_ts_stereo = (stereoTsosGMom.x()-stereoShitGMom.x())/sqrt(stereoTsosGMEr.cxx());
922 // double pullGMY_ts_stereo = (stereoTsosGMom.y()-stereoShitGMom.y())/sqrt(stereoTsosGMEr.cyy());
923 // double pullGMZ_ts_stereo = (stereoTsosGMom.z()-stereoShitGMom.z())/sqrt(stereoTsosGMEr.czz());
924 // // double pullGMX_ts_stereo = (stereoTsosGMom.x()-stereoShitGMom.x())/*/sqrt(stereoTsosGMEr.cxx())*/;
925 // // double pullGMY_ts_stereo = (stereoTsosGMom.y()-stereoShitGMom.y())/*/sqrt(stereoTsosGMEr.cyy())*/;
926 // // double pullGMZ_ts_stereo = (stereoTsosGMom.z()-stereoShitGMom.z())/*/sqrt(stereoTsosGMEr.czz())*/;
927 
928 // title.str("");
929 // title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_stereo";
930 // hPullGM_X_ts_stereo[title.str()]->Fill( pullGMX_ts_stereo );
931 // title.str("");
932 // title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_stereo";
933 // hPullGM_Y_ts_stereo[title.str()]->Fill( pullGMY_ts_stereo );
934 // title.str("");
935 // title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_stereo";
936 // hPullGM_Z_ts_stereo[title.str()]->Fill( pullGMZ_ts_stereo );
937 // }
938 // }
939 // }
940 // delete hitAssociator;
941 // LogTrace("TestHits") << "end of event" << std::endl;
942 // }
943 
944 void TestHits::endJob() {
945  //file->Write();
946  TDirectory * chi2i = file->mkdir("Chi2_Increment");
947 
948  TDirectory * gp_ts = file->mkdir("GP_TSOS-SimHit");
949  TDirectory * gm_ts = file->mkdir("GM_TSOS-SimHit");
950  TDirectory * gp_tr = file->mkdir("GP_TSOS-RecHit");
951  TDirectory * gp_rs = file->mkdir("GP_RecHit-SimHit");
952 
953  TDirectory * gp_tsx = gp_ts->mkdir("X");
954  TDirectory * gp_tsy = gp_ts->mkdir("Y");
955  TDirectory * gp_tsz = gp_ts->mkdir("Z");
956  TDirectory * gm_tsx = gm_ts->mkdir("X");
957  TDirectory * gm_tsy = gm_ts->mkdir("Y");
958  TDirectory * gm_tsz = gm_ts->mkdir("Z");
959  TDirectory * gp_trx = gp_tr->mkdir("X");
960  TDirectory * gp_try = gp_tr->mkdir("Y");
961  TDirectory * gp_trz = gp_tr->mkdir("Z");
962  TDirectory * gp_rsx = gp_rs->mkdir("X");
963  TDirectory * gp_rsy = gp_rs->mkdir("Y");
964  TDirectory * gp_rsz = gp_rs->mkdir("Z");
965 
966  TDirectory * gp_tsx_mono = gp_ts->mkdir("MONOX");
967  TDirectory * gp_tsy_mono = gp_ts->mkdir("MONOY");
968  TDirectory * gp_tsz_mono = gp_ts->mkdir("MONOZ");
969  TDirectory * gm_tsx_mono = gm_ts->mkdir("MONOX");
970  TDirectory * gm_tsy_mono = gm_ts->mkdir("MONOY");
971  TDirectory * gm_tsz_mono = gm_ts->mkdir("MONOZ");
972  TDirectory * gp_trx_mono = gp_tr->mkdir("MONOX");
973  TDirectory * gp_try_mono = gp_tr->mkdir("MONOY");
974  TDirectory * gp_trz_mono = gp_tr->mkdir("MONOZ");
975  TDirectory * gp_rsx_mono = gp_rs->mkdir("MONOX");
976  TDirectory * gp_rsy_mono = gp_rs->mkdir("MONOY");
977  TDirectory * gp_rsz_mono = gp_rs->mkdir("MONOZ");
978 
979  TDirectory * gp_tsx_stereo = gp_ts->mkdir("STEREOX");
980  TDirectory * gp_tsy_stereo = gp_ts->mkdir("STEREOY");
981  TDirectory * gp_tsz_stereo = gp_ts->mkdir("STEREOZ");
982  TDirectory * gm_tsx_stereo = gm_ts->mkdir("STEREOX");
983  TDirectory * gm_tsy_stereo = gm_ts->mkdir("STEREOY");
984  TDirectory * gm_tsz_stereo = gm_ts->mkdir("STEREOZ");
985  TDirectory * gp_trx_stereo = gp_tr->mkdir("STEREOX");
986  TDirectory * gp_try_stereo = gp_tr->mkdir("STEREOY");
987  TDirectory * gp_trz_stereo = gp_tr->mkdir("STEREOZ");
988  TDirectory * gp_rsx_stereo = gp_rs->mkdir("STEREOX");
989  TDirectory * gp_rsy_stereo = gp_rs->mkdir("STEREOY");
990  TDirectory * gp_rsz_stereo = gp_rs->mkdir("STEREOZ");
991 
992  chi2i->cd();
993  hTotChi2Increment->Write();
994  hProcess_vs_Chi2->Write();
995  hClsize_vs_Chi2->Write();
996  for (int i=0; i!=6; i++)
997  for (int j=0; j!=9; j++){
998  if (i==0 && j>2) break;
999  if (i==1 && j>1) break;
1000  if (i==2 && j>3) break;
1001  if (i==3 && j>2) break;
1002  if (i==4 && j>5) break;
1003  if (i==5 && j>8) break;
1004  chi2i->cd();
1005  title.str("");
1006  title << "Chi2Increment_" << i+1 << "-" << j+1;
1007  hChi2Increment[title.str()]->Write();
1008 
1009  gp_ts->cd();
1010  gp_tsx->cd();
1011  title.str("");
1012  title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts";
1013  hPullGP_X_ts[title.str()]->Write();
1014  gp_tsy->cd();
1015  title.str("");
1016  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts";
1017  hPullGP_Y_ts[title.str()]->Write();
1018  gp_tsz->cd();
1019  title.str("");
1020  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts";
1021  hPullGP_Z_ts[title.str()]->Write();
1022 
1023  gm_ts->cd();
1024  gm_tsx->cd();
1025  title.str("");
1026  title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts";
1027  hPullGM_X_ts[title.str()]->Write();
1028  gm_tsy->cd();
1029  title.str("");
1030  title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts";
1031  hPullGM_Y_ts[title.str()]->Write();
1032  gm_tsz->cd();
1033  title.str("");
1034  title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts";
1035  hPullGM_Z_ts[title.str()]->Write();
1036 
1037  gp_tr->cd();
1038  gp_trx->cd();
1039  title.str("");
1040  title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr";
1041  hPullGP_X_tr[title.str()]->Write();
1042  gp_try->cd();
1043  title.str("");
1044  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr";
1045  hPullGP_Y_tr[title.str()]->Write();
1046  gp_trz->cd();
1047  title.str("");
1048  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr";
1049  hPullGP_Z_tr[title.str()]->Write();
1050 
1051  gp_rs->cd();
1052  gp_rsx->cd();
1053  title.str("");
1054  title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs";
1055  hPullGP_X_rs[title.str()]->Write();
1056  gp_rsy->cd();
1057  title.str("");
1058  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs";
1059  hPullGP_Y_rs[title.str()]->Write();
1060  gp_rsz->cd();
1061  title.str("");
1062  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs";
1063  hPullGP_Z_rs[title.str()]->Write();
1064 
1065  if ( ((i==2||i==4)&&(j==0||j==1)) || (i==3||i==5) ){
1066  //mono
1067  gp_ts->cd();
1068  gp_tsx_mono->cd();
1069  title.str("");
1070  title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts_mono";
1071  hPullGP_X_ts_mono[title.str()]->Write();
1072  gp_tsy_mono->cd();
1073  title.str("");
1074  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts_mono";
1075  hPullGP_Y_ts_mono[title.str()]->Write();
1076  gp_tsz_mono->cd();
1077  title.str("");
1078  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts_mono";
1079  hPullGP_Z_ts_mono[title.str()]->Write();
1080 
1081  gm_ts->cd();
1082  gm_tsx_mono->cd();
1083  title.str("");
1084  title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts_mono";
1085  hPullGM_X_ts_mono[title.str()]->Write();
1086  gm_tsy_mono->cd();
1087  title.str("");
1088  title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts_mono";
1089  hPullGM_Y_ts_mono[title.str()]->Write();
1090  gm_tsz_mono->cd();
1091  title.str("");
1092  title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts_mono";
1093  hPullGM_Z_ts_mono[title.str()]->Write();
1094 
1095  gp_tr->cd();
1096  gp_trx_mono->cd();
1097  title.str("");
1098  title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr_mono";
1099  hPullGP_X_tr_mono[title.str()]->Write();
1100  gp_try_mono->cd();
1101  title.str("");
1102  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr_mono";
1103  hPullGP_Y_tr_mono[title.str()]->Write();
1104  gp_trz_mono->cd();
1105  title.str("");
1106  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr_mono";
1107  hPullGP_Z_tr_mono[title.str()]->Write();
1108 
1109  gp_rs->cd();
1110  gp_rsx_mono->cd();
1111  title.str("");
1112  title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs_mono";
1113  hPullGP_X_rs_mono[title.str()]->Write();
1114  gp_rsy_mono->cd();
1115  title.str("");
1116  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs_mono";
1117  hPullGP_Y_rs_mono[title.str()]->Write();
1118  gp_rsz_mono->cd();
1119  title.str("");
1120  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs_mono";
1121  hPullGP_Z_rs_mono[title.str()]->Write();
1122 
1123  //stereo
1124  gp_ts->cd();
1125  gp_tsx_stereo->cd();
1126  title.str("");
1127  title << "PullGP_X_" << i+1 << "-" << j+1 << "_ts_stereo";
1128  hPullGP_X_ts_stereo[title.str()]->Write();
1129  gp_tsy_stereo->cd();
1130  title.str("");
1131  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_ts_stereo";
1132  hPullGP_Y_ts_stereo[title.str()]->Write();
1133  gp_tsz_stereo->cd();
1134  title.str("");
1135  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_ts_stereo";
1136  hPullGP_Z_ts_stereo[title.str()]->Write();
1137 
1138  gm_ts->cd();
1139  gm_tsx_stereo->cd();
1140  title.str("");
1141  title << "PullGM_X_" << i+1 << "-" << j+1 << "_ts_stereo";
1142  hPullGM_X_ts_stereo[title.str()]->Write();
1143  gm_tsy_stereo->cd();
1144  title.str("");
1145  title << "PullGM_Y_" << i+1 << "-" << j+1 << "_ts_stereo";
1146  hPullGM_Y_ts_stereo[title.str()]->Write();
1147  gm_tsz_stereo->cd();
1148  title.str("");
1149  title << "PullGM_Z_" << i+1 << "-" << j+1 << "_ts_stereo";
1150  hPullGM_Z_ts_stereo[title.str()]->Write();
1151 
1152  gp_tr->cd();
1153  gp_trx_stereo->cd();
1154  title.str("");
1155  title << "PullGP_X_" << i+1 << "-" << j+1 << "_tr_stereo";
1156  hPullGP_X_tr_stereo[title.str()]->Write();
1157  gp_try_stereo->cd();
1158  title.str("");
1159  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_tr_stereo";
1160  hPullGP_Y_tr_stereo[title.str()]->Write();
1161  gp_trz_stereo->cd();
1162  title.str("");
1163  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_tr_stereo";
1164  hPullGP_Z_tr_stereo[title.str()]->Write();
1165 
1166  gp_rs->cd();
1167  gp_rsx_stereo->cd();
1168  title.str("");
1169  title << "PullGP_X_" << i+1 << "-" << j+1 << "_rs_stereo";
1170  hPullGP_X_rs_stereo[title.str()]->Write();
1171  gp_rsy_stereo->cd();
1172  title.str("");
1173  title << "PullGP_Y_" << i+1 << "-" << j+1 << "_rs_stereo";
1174  hPullGP_Y_rs_stereo[title.str()]->Write();
1175  gp_rsz_stereo->cd();
1176  title.str("");
1177  title << "PullGP_Z_" << i+1 << "-" << j+1 << "_rs_stereo";
1178  hPullGP_Z_rs_stereo[title.str()]->Write();
1179  }
1180  }
1181 
1182  file->Close();
1183 }
1184 
1185 //needed by to do the residual for matched hits
1186 //taken from SiStripTrackingRecHitsValid.cc
1187 std::pair<LocalPoint,LocalVector>
1188 TestHits::projectHit( const PSimHit& hit, const StripGeomDetUnit* stripDet, const BoundPlane& plane)
1189 {
1190  const StripTopology& topol = stripDet->specificTopology();
1191  GlobalPoint globalpos= stripDet->surface().toGlobal(hit.localPosition());
1192  LocalPoint localHit = plane.toLocal(globalpos);
1193  //track direction
1194  LocalVector locdir=hit.localDirection();
1195  //rotate track in new frame
1196 
1197  GlobalVector globaldir= stripDet->surface().toGlobal(locdir);
1198  LocalVector dir=plane.toLocal(globaldir);
1199  float scale = -localHit.z() / dir.z();
1200 
1201  LocalPoint projectedPos = localHit + scale*dir;
1202 
1203  float selfAngle = topol.stripAngle( topol.strip( hit.localPosition()));
1204 
1205  LocalVector stripDir( sin(selfAngle), cos(selfAngle), 0); // vector along strip in hit frame
1206 
1207  LocalVector localStripDir( plane.toLocal(stripDet->surface().toGlobal( stripDir)));
1208 
1209  return std::pair<LocalPoint,LocalVector>( projectedPos, localStripDir);
1210 }
1211 
1212 #include "FWCore/Framework/interface/ModuleFactory.h"
1213 #include "FWCore/Framework/interface/MakerMacros.h"
1214 
1215 DEFINE_FWK_MODULE(TestHits);
TestHits::theTCCollection
edm::Handle< TrackCandidateCollection > theTCCollection
Definition: TestHits.h:77
delta
dbl * delta
Definition: mlp_gen.cc:36
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:78
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
i
int i
Definition: DBlmapReader.cc:9
TestHits::projectHit
std::pair< LocalPoint, LocalVector > projectHit(const PSimHit &, const StripGeomDetUnit *, const BoundPlane &)
Definition: TestHits.cc:1188
TestHits::beginRun
virtual void beginRun(edm::Run &run, const edm::EventSetup &)
Definition: TestHits.cc:41
TestHits::hPullGP_Y_tr_mono
std::map< std::string, TH1F * > hPullGP_Y_tr_mono
Definition: TestHits.h:107
TestHits::hPullGP_X_rs_mono
std::map< std::string, TH1F * > hPullGP_X_rs_mono
Definition: TestHits.h:103
TestHits::hPullGP_Z_tr_stereo
std::map< std::string, TH1F * > hPullGP_Z_tr_stereo
Definition: TestHits.h:121
TestHits::hTotChi2Increment
TH1F * hTotChi2Increment
Definition: TestHits.h:94
StripTopology::stripAngle
virtual float stripAngle(float strip) const =0
TrackCandidate::recHits
range recHits() const
Definition: TrackCandidate.h:54
TestHits::theBuilder
edm::ESHandle< TransientTrackingRecHitBuilder > theBuilder
Definition: TestHits.h:75
TestHits::hPullGP_Z_ts_mono
std::map< std::string, TH1F * > hPullGP_Z_ts_mono
Definition: TestHits.h:99
PSimHit::momentumAtEntry
LocalVector momentumAtEntry() const
The momentum of the track that produced the hit, at entry point.
Definition: PSimHit.h:47
TestHits::hPullGM_X_ts_stereo
std::map< std::string, TH1F * > hPullGM_X_ts_stereo
Definition: TestHits.h:113
mag
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
Definition: Basic3DVectorLD.h:142
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
TestHits::hPullGM_X_ts_mono
std::map< std::string, TH1F * > hPullGM_X_ts_mono
Definition: TestHits.h:100
TestHits::builderName
std::string builderName
Definition: TestHits.h:68
TestHits::mineta
double mineta
Definition: TestHits.h:65
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
TrackCandidate::seed
TrajectorySeed const & seed() const
Definition: TrackCandidate.h:56
TrajectoryStateOnSurface::cartesianError
const CartesianTrajectoryError cartesianError() const
Definition: TrajectoryStateOnSurface.h:161
PV3DBase::y
T y() const
Definition: PV3DBase.h:62
TrajectoryStateOnSurface::globalPosition
GlobalPoint globalPosition() const
Definition: TrajectoryStateOnSurface.h:143
TestHits::hPullGP_Y_ts
std::map< std::string, TH1F * > hPullGP_Y_ts
Definition: TestHits.h:82
TestHits::hChi2Increment
std::map< std::string, TH1F * > hChi2Increment
Definition: TestHits.h:93
TrackCandidate::range
std::pair< const_iterator, const_iterator > range
Definition: TrackCandidate.h:27
TestHits::hPullGP_Y_ts_stereo
std::map< std::string, TH1F * > hPullGP_Y_ts_stereo
Definition: TestHits.h:111
TestHits::hPullGP_Z_rs
std::map< std::string, TH1F * > hPullGP_Z_rs
Definition: TestHits.h:89
StripGeomDetUnit::specificTopology
virtual const StripTopology & specificTopology() const
Returns a reference to the strip proxy topology.
Definition: StripGeomDetUnit.cc:19
TestHits::thePropagator
edm::ESHandle< Propagator > thePropagator
Definition: TestHits.h:74
TestHits::hPullGP_Y_rs
std::map< std::string, TH1F * > hPullGP_Y_rs
Definition: TestHits.h:88
TransientTrackingRecHit::RecHitContainer
std::vector< ConstRecHitPointer > RecHitContainer
Definition: TransientTrackingRecHit.h:29
TestHits::hPullGP_X_tr_stereo
std::map< std::string, TH1F * > hPullGP_X_tr_stereo
Definition: TestHits.h:119
StripTopology::strip
virtual float strip(const LocalPoint &) const =0
TestHits::hPullGP_Z_rs_stereo
std::map< std::string, TH1F * > hPullGP_Z_rs_stereo
Definition: TestHits.h:118
TrackCandidate::trajectoryStateOnDet
PTrajectoryStateOnDet const & trajectoryStateOnDet() const
Definition: TrackCandidate.h:52
errorMatrix2Lands_multiChannel.id
tuple id
Definition: errorMatrix2Lands_multiChannel.py:151
TestHits::hitAssociator
TrackerHitAssociator * hitAssociator
Definition: TestHits.h:63
iEvent
int iEvent
Definition: GenABIO.cc:243
AlgebraicSymMatrix33
ROOT::Math::SMatrix< double, 3, 3, ROOT::Math::MatRepSym< double, 3 > > AlgebraicSymMatrix33
Definition: AlgebraicROOTObjects.h:19
TestHits::hPullGP_Y_tr
std::map< std::string, TH1F * > hPullGP_Y_tr
Definition: TestHits.h:91
PSimHit::localPosition
Local3DPoint localPosition() const
Definition: PSimHit.h:44
TestHits::propagatorName
std::string propagatorName
Definition: TestHits.h:67
pileupReCalc_HLTpaths.scale
int scale
Definition: pileupReCalc_HLTpaths.py:111
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:46
GloballyPositioned::toLocal
LocalPoint toLocal(const GlobalPoint &gp) const
Definition: GloballyPositioned.h:118
TestHits::hPullGP_X_rs_stereo
std::map< std::string, TH1F * > hPullGP_X_rs_stereo
Definition: TestHits.h:116
PV3DBase::z
T z() const
Definition: PV3DBase.h:63
query.result
tuple result
Definition: query.py:137
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
heavyFlavorValidationHarvestingSequence_cff.combiner
tuple combiner
Definition: heavyFlavorValidationHarvestingSequence_cff.py:22
j
int j
Definition: DBlmapReader.cc:9
TestHits::endJob
virtual void endJob()
Definition: TestHits.cc:944
TestHits::hPullGP_Y_tr_stereo
std::map< std::string, TH1F * > hPullGP_Y_tr_stereo
Definition: TestHits.h:120
CTTRHp
TransientTrackingRecHit::ConstRecHitPointer CTTRHp
Definition: CkfDebugger.h:44
TestHits::hPullGP_Y_rs_mono
std::map< std::string, TH1F * > hPullGP_Y_rs_mono
Definition: TestHits.h:104
PTrajectoryStateOnDet::detId
unsigned int detId() const
Definition: PTrajectoryStateOnDet.h:69
TestHits::hPullGM_Z_ts
std::map< std::string, TH1F * > hPullGM_Z_ts
Definition: TestHits.h:86
TestHits::hPullGM_Y_ts_mono
std::map< std::string, TH1F * > hPullGM_Y_ts_mono
Definition: TestHits.h:101
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
TestHits::hPullGP_Z_rs_mono
std::map< std::string, TH1F * > hPullGP_Z_rs_mono
Definition: TestHits.h:105
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
TestHits::hPullGP_X_ts_stereo
std::map< std::string, TH1F * > hPullGP_X_ts_stereo
Definition: TestHits.h:110
TestHits::hPullGP_Y_rs_stereo
std::map< std::string, TH1F * > hPullGP_Y_rs_stereo
Definition: TestHits.h:117
CartesianTrajectoryError::matrix
const AlgebraicSymMatrix66 & matrix() const
Definition: CartesianTrajectoryError.h:29
TestHits::theG
edm::ESHandle< TrackerGeometry > theG
Definition: TestHits.h:72
TestHits::fit
edm::ESHandle< TrajectoryFitter > fit
Definition: TestHits.h:76
PSimHit::localDirection
LocalVector localDirection() const
Obsolete. Same as momentumAtEntry().unit(), for backward compatibility.
Definition: PSimHit.h:52
DetId
Definition: DetId.h:20
TestHits::hPullGM_Z_ts_stereo
std::map< std::string, TH1F * > hPullGM_Z_ts_stereo
Definition: TestHits.h:115
TSOS
TrajectoryStateOnSurface TSOS
Definition: TestHits.cc:23
trajectoryStateTransform::transientState
TrajectoryStateOnSurface transientState(const PTrajectoryStateOnDet &ts, const Surface *surface, const MagneticField *field)
Definition: TrajectoryStateTransform.cc:43
TestHits::~TestHits
~TestHits()
Definition: TestHits.cc:39
TestHits::hProcess_vs_Chi2
TH2F * hProcess_vs_Chi2
Definition: TestHits.h:95
CartesianTrajectoryError::position
const GlobalError position() const
Position error submatrix.
Definition: CartesianTrajectoryError.h:46
TestHits::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: TestHits.cc:192
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
TestHits::hPullGP_Z_tr
std::map< std::string, TH1F * > hPullGP_Z_tr
Definition: TestHits.h:92
TestHits::hPullGP_Y_ts_mono
std::map< std::string, TH1F * > hPullGP_Y_ts_mono
Definition: TestHits.h:98
PSimHit::processType
unsigned short processType() const
Definition: PSimHit.h:118
evf::utils::state
char state
Definition: procUtils.cc:75
TestHits::hPullGM_Y_ts_stereo
std::map< std::string, TH1F * > hPullGM_Y_ts_stereo
Definition: TestHits.h:114
TestHits::theMF
edm::ESHandle< MagneticField > theMF
Definition: TestHits.h:73
TestHits::conf_
const edm::ParameterSet conf_
Definition: TestHits.h:62
TestHits::hPullGP_X_tr_mono
std::map< std::string, TH1F * > hPullGP_X_tr_mono
Definition: TestHits.h:106
TestHits::hClsize_vs_Chi2
TH2F * hClsize_vs_Chi2
Definition: TestHits.h:95
TestHits::file
TFile * file
Definition: TestHits.h:79
TrajectoryStateOnSurface::globalMomentum
GlobalVector globalMomentum() const
Definition: TrajectoryStateOnSurface.h:146
GeomDet::surface
const BoundPlane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:35
TestHits::hPullGM_X_ts
std::map< std::string, TH1F * > hPullGM_X_ts
Definition: TestHits.h:84
edm::OwnVector::const_iterator
Definition: OwnVector.h:42
TestHits::hPullGP_X_ts
std::map< std::string, TH1F * > hPullGP_X_ts
Definition: TestHits.h:81
TrackerHitAssociator::associateHit
std::vector< PSimHit > associateHit(const TrackingRecHit &thit)
Definition: TrackerHitAssociator.cc:120
TestHits::fname
std::string fname
Definition: TestHits.h:70
TestHits::srcName
std::string srcName
Definition: TestHits.h:69
TestHits::hPullGP_Z_ts
std::map< std::string, TH1F * > hPullGP_Z_ts
Definition: TestHits.h:83
dir
dbl *** dir
Definition: mlp_gen.cc:35
TestHits::TestHits
TestHits(const edm::ParameterSet &)
Definition: TestHits.cc:28
TestHits::hPullGP_Z_tr_mono
std::map< std::string, TH1F * > hPullGP_Z_tr_mono
Definition: TestHits.h:108
TestHits::hPullGP_X_ts_mono
std::map< std::string, TH1F * > hPullGP_X_ts_mono
Definition: TestHits.h:97
PV3DBase::x
T x() const
Definition: PV3DBase.h:61
GloballyPositioned::position
const PositionType & position() const
Definition: GloballyPositioned.h:41
TestHits::hPullGM_Z_ts_mono
std::map< std::string, TH1F * > hPullGM_Z_ts_mono
Definition: TestHits.h:102
TestHits::hPullGM_Y_ts
std::map< std::string, TH1F * > hPullGM_Y_ts
Definition: TestHits.h:85
TestHits::hPullGP_X_tr
std::map< std::string, TH1F * > hPullGP_X_tr
Definition: TestHits.h:90
edm::Run
Definition: Run.h:33
TestHits::title
std::stringstream title
Definition: TestHits.h:80
SiPixelRecHit
Pixel Reconstructed Hit.
TestHits::maxeta
double maxeta
Definition: TestHits.h:65
TestHits::hPullGP_Z_ts_stereo
std::map< std::string, TH1F * > hPullGP_Z_ts_stereo
Definition: TestHits.h:112
TestHits::hPullGP_X_rs
std::map< std::string, TH1F * > hPullGP_X_rs
Definition: TestHits.h:87
Generated for CMSSW Reference Manual by doxygen 1.8.5