CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
TestHits Class Reference

#include <RecoTracker/DebugTools/src/TestHits.cc>

Inheritance diagram for TestHits:
edm::EDAnalyzer edm::EDConsumerBase

Public Member Functions

 TestHits (const edm::ParameterSet &)
 
 ~TestHits () override
 
- Public Member Functions inherited from edm::EDAnalyzer
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzer ()
 
SerialTaskQueue * globalLuminosityBlocksQueue ()
 
SerialTaskQueue * globalRunsQueue ()
 
ModuleDescription const & moduleDescription () const
 
std::string workerType () const
 
 ~EDAnalyzer () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase &&)=default
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
ESProxyIndex const * esGetTokenIndices (edm::Transition iTrans) const
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
void updateLookup (eventsetup::ESRecordsToProxyIndices const &)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void beginRun (edm::Run const &run, const edm::EventSetup &) override
 
void endJob () override
 
std::pair< LocalPoint, LocalVector > projectHit (const PSimHit &, const StripGeomDetUnit *, const BoundPlane &)
 

Private Attributes

std::string builderName
 
TFile * file
 
edm::ESHandle< TrajectoryFitter > fit
 
std::string fname
 
std::map< std::string, TH1F * > hChi2Increment
 
TH2F * hClsize_vs_Chi2
 
TH2F * hProcess_vs_Chi2
 
std::map< std::string, TH1F * > hPullGM_X_ts
 
std::map< std::string, TH1F * > hPullGM_X_ts_mono
 
std::map< std::string, TH1F * > hPullGM_X_ts_stereo
 
std::map< std::string, TH1F * > hPullGM_Y_ts
 
std::map< std::string, TH1F * > hPullGM_Y_ts_mono
 
std::map< std::string, TH1F * > hPullGM_Y_ts_stereo
 
std::map< std::string, TH1F * > hPullGM_Z_ts
 
std::map< std::string, TH1F * > hPullGM_Z_ts_mono
 
std::map< std::string, TH1F * > hPullGM_Z_ts_stereo
 
std::map< std::string, TH1F * > hPullGP_X_rs
 
std::map< std::string, TH1F * > hPullGP_X_rs_mono
 
std::map< std::string, TH1F * > hPullGP_X_rs_stereo
 
std::map< std::string, TH1F * > hPullGP_X_tr
 
std::map< std::string, TH1F * > hPullGP_X_tr_mono
 
std::map< std::string, TH1F * > hPullGP_X_tr_stereo
 
std::map< std::string, TH1F * > hPullGP_X_ts
 
std::map< std::string, TH1F * > hPullGP_X_ts_mono
 
std::map< std::string, TH1F * > hPullGP_X_ts_stereo
 
std::map< std::string, TH1F * > hPullGP_Y_rs
 
std::map< std::string, TH1F * > hPullGP_Y_rs_mono
 
std::map< std::string, TH1F * > hPullGP_Y_rs_stereo
 
std::map< std::string, TH1F * > hPullGP_Y_tr
 
std::map< std::string, TH1F * > hPullGP_Y_tr_mono
 
std::map< std::string, TH1F * > hPullGP_Y_tr_stereo
 
std::map< std::string, TH1F * > hPullGP_Y_ts
 
std::map< std::string, TH1F * > hPullGP_Y_ts_mono
 
std::map< std::string, TH1F * > hPullGP_Y_ts_stereo
 
std::map< std::string, TH1F * > hPullGP_Z_rs
 
std::map< std::string, TH1F * > hPullGP_Z_rs_mono
 
std::map< std::string, TH1F * > hPullGP_Z_rs_stereo
 
std::map< std::string, TH1F * > hPullGP_Z_tr
 
std::map< std::string, TH1F * > hPullGP_Z_tr_mono
 
std::map< std::string, TH1F * > hPullGP_Z_tr_stereo
 
std::map< std::string, TH1F * > hPullGP_Z_ts
 
std::map< std::string, TH1F * > hPullGP_Z_ts_mono
 
std::map< std::string, TH1F * > hPullGP_Z_ts_stereo
 
TH1F * hTotChi2Increment
 
double maxeta
 
double mineta
 
std::string propagatorName
 
std::string srcName
 
edm::ESHandle< TransientTrackingRecHitBuilder > theBuilder
 
edm::ESHandle< TrackerGeometry > theG
 
edm::ESHandle< MagneticField > theMF
 
edm::ESHandle< Propagator > thePropagator
 
edm::Handle< TrackCandidateCollection > theTCCollection
 
std::stringstream title
 
TrackerHitAssociator::Config trackerHitAssociatorConfig_
 

Additional Inherited Members

- Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &)
 
static bool wantsGlobalLuminosityBlocks ()
 
static bool wantsGlobalRuns ()
 
static bool wantsStreamLuminosityBlocks ()
 
static bool wantsStreamRuns ()
 
- Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes ()
 
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto esConsumes (ESInputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 

Detailed Description

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 49 of file TestHits.h.

Constructor & Destructor Documentation

◆ TestHits()

TestHits::TestHits ( const edm::ParameterSet &  iConfig)
explicit

Definition at line 24 of file TestHits.cc.

24  : trackerHitAssociatorConfig_(consumesCollector()) {
25  LogTrace("TestHits") << iConfig << std::endl;
26  propagatorName = iConfig.getParameter<std::string>("Propagator");
27  builderName = iConfig.getParameter<std::string>("TTRHBuilder");
28  srcName = iConfig.getParameter<std::string>("src");
29  fname = iConfig.getParameter<std::string>("Fitter");
30  mineta = iConfig.getParameter<double>("mineta");
31  maxeta = iConfig.getParameter<double>("maxeta");
32 }

References builderName, fname, edm::ParameterSet::getParameter(), LogTrace, maxeta, mineta, propagatorName, srcName, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~TestHits()

TestHits::~TestHits ( )
override

Definition at line 34 of file TestHits.cc.

34 {}

Member Function Documentation

◆ analyze()

void TestHits::analyze ( const edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
)
overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 190 of file TestHits.cc.

190  {
191  //Retrieve tracker topology from geometry
192  edm::ESHandle<TrackerTopology> tTopo;
193  iSetup.get<TrackerTopologyRcd>().get(tTopo);
194 
195  LogTrace("TestHits") << "\nnew event";
196 
197  iEvent.getByLabel(srcName, theTCCollection);
198  TrackerHitAssociator hitAssociator(iEvent, trackerHitAssociatorConfig_);
199 
200  TrajectoryStateCombiner combiner;
201 
202  for (TrackCandidateCollection::const_iterator i = theTCCollection->begin(); i != theTCCollection->end(); i++) {
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  DetId detId(state.detId());
212  TrajectoryStateOnSurface theTSOS =
213  trajectoryStateTransform::transientState(state, &(theG->idToDet(detId)->surface()), theMF.product());
214 
215  if (theTSOS.globalMomentum().eta() > maxeta || theTSOS.globalMomentum().eta() < mineta)
216  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; i != recHitVec.second; i++) {
222  hits.push_back(theBuilder->build(&(*i)));
223  }
224 
225  //call the fitter
226  std::vector<Trajectory> result = fit->fit(theTC->seed(), hits, theTSOS);
227  if (result.empty())
228  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  TransientTrackingRecHit::ConstRecHitPointer rhit = tm->recHit();
235  if ((rhit)->isValid() == 0 && rhit->det() != nullptr)
236  continue;
237  LogTrace("TestHits") << "*****************new hit*****************";
238 
239  int subdetId = rhit->det()->geographicalId().subdetId();
240  DetId id = rhit->det()->geographicalId();
241  int layerId = tTopo->layer(id);
242  LogTrace("TestHits") << "subdetId=" << subdetId << " layerId=" << layerId;
243 
244  double delta = 99999;
245  LocalPoint rhitLPv = rhit->localPosition();
246 
247  std::vector<PSimHit> assSimHits = hitAssociator.associateHit(*(rhit->hit()));
248  if (assSimHits.empty())
249  continue;
250  PSimHit shit;
251  for (std::vector<PSimHit>::const_iterator m = assSimHits.begin(); m < assSimHits.end(); m++) {
252  if ((m->localPosition() - rhitLPv).mag() < delta) {
253  shit = *m;
254  delta = (m->localPosition() - rhitLPv).mag();
255  }
256  }
257 
258  TSOS currentState = tm->forwardPredictedState();
259  if (tm->backwardPredictedState().isValid())
260  currentState = combiner(tm->backwardPredictedState(), tm->forwardPredictedState());
261  TSOS updatedState = tm->updatedState();
262  tchi2 += tm->estimate();
263 
264  //plot chi2 increment
265  double chi2increment = tm->estimate();
266  LogTrace("TestHits") << "tm->estimate()=" << tm->estimate();
267  title.str("");
268  title << "Chi2Increment_" << subdetId << "-" << layerId;
269  hChi2Increment[title.str()]->Fill(chi2increment);
270  hTotChi2Increment->Fill(chi2increment);
271  hProcess_vs_Chi2->Fill(chi2increment, shit.processType());
272  if (dynamic_cast<const SiPixelRecHit*>(rhit->hit()))
273  hClsize_vs_Chi2->Fill(chi2increment, ((const SiPixelRecHit*)(rhit->hit()))->cluster()->size());
274  if (dynamic_cast<const SiStripRecHit2D*>(rhit->hit()))
275  hClsize_vs_Chi2->Fill(chi2increment, ((const SiStripRecHit2D*)(rhit->hit()))->cluster()->amplitudes().size());
276 
277  //test hits
278  const Surface* surf = &((rhit)->det()->surface());
279  LocalVector shitLMom;
280  LocalPoint shitLPos;
281  if (dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())) {
282  double rechitmatchedx = rhit->localPosition().x();
283  double rechitmatchedy = rhit->localPosition().y();
284  double mindist = 999999;
285  double distx, disty;
286  std::pair<LocalPoint, LocalVector> closestPair;
287  const StripGeomDetUnit* stripDet = (StripGeomDetUnit*)((const GluedGeomDet*)(rhit)->det())->stereoDet();
288  const BoundPlane& plane = (rhit)->det()->surface();
289  for (std::vector<PSimHit>::const_iterator m = assSimHits.begin(); m < assSimHits.end(); m++) {
290  //project simhit;
291  std::pair<LocalPoint, LocalVector> hitPair = projectHit((*m), stripDet, plane);
292  distx = fabs(rechitmatchedx - hitPair.first.x());
293  disty = fabs(rechitmatchedy - hitPair.first.y());
294  double dist = distx * distx + disty * disty;
295  if (sqrt(dist) < mindist) {
296  mindist = dist;
297  closestPair = hitPair;
298  }
299  }
300  shitLPos = closestPair.first;
301  shitLMom = closestPair.second;
302  } else {
303  shitLPos = shit.localPosition();
304  shitLMom = shit.momentumAtEntry();
305  }
306  GlobalVector shitGMom = surf->toGlobal(shitLMom);
307  GlobalPoint shitGPos = surf->toGlobal(shitLPos);
308 
309  GlobalVector tsosGMom = currentState.globalMomentum();
310  GlobalError tsosGMEr(currentState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3, 3));
311  GlobalPoint tsosGPos = currentState.globalPosition();
312  GlobalError tsosGPEr = currentState.cartesianError().position();
313 
314  GlobalPoint rhitGPos = (rhit)->globalPosition();
315  GlobalError rhitGPEr = (rhit)->globalPositionError();
316 
317  double pullGPX_rs = (rhitGPos.x() - shitGPos.x()) / sqrt(rhitGPEr.cxx());
318  double pullGPY_rs = (rhitGPos.y() - shitGPos.y()) / sqrt(rhitGPEr.cyy());
319  double pullGPZ_rs = (rhitGPos.z() - shitGPos.z()) / sqrt(rhitGPEr.czz());
320  //double pullGPX_rs = (rhitGPos.x()-shitGPos.x());
321  //double pullGPY_rs = (rhitGPos.y()-shitGPos.y());
322  //double pullGPZ_rs = (rhitGPos.z()-shitGPos.z());
323 
324  LogTrace("TestHits") << "rs" << std::endl;
325  LogVerbatim("TestHits") << "assSimHits.size()=" << assSimHits.size();
326  LogVerbatim("TestHits") << "tsos globalPos =" << tsosGPos;
327  LogVerbatim("TestHits") << "sim hit globalPos=" << shitGPos;
328  LogVerbatim("TestHits") << "rec hit globalPos=" << rhitGPos;
329  LogVerbatim("TestHits") << "geographicalId =" << rhit->det()->geographicalId().rawId();
330  LogVerbatim("TestHits") << "surface position =" << surf->position();
331 
332  title.str("");
333  title << "PullGP_X_" << subdetId << "-" << layerId << "_rs";
334  hPullGP_X_rs[title.str()]->Fill(pullGPX_rs);
335  title.str("");
336  title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs";
337  hPullGP_Y_rs[title.str()]->Fill(pullGPY_rs);
338  title.str("");
339  title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs";
340  hPullGP_Z_rs[title.str()]->Fill(pullGPZ_rs);
341 
342  double pullGPX_tr = (tsosGPos.x() - rhitGPos.x()) / sqrt(tsosGPEr.cxx() + rhitGPEr.cxx());
343  double pullGPY_tr = (tsosGPos.y() - rhitGPos.y()) / sqrt(tsosGPEr.cyy() + rhitGPEr.cyy());
344  double pullGPZ_tr = (tsosGPos.z() - rhitGPos.z()) / sqrt(tsosGPEr.czz() + rhitGPEr.czz());
345  //double pullGPX_tr = (tsosGPos.x()-rhitGPos.x());
346  //double pullGPY_tr = (tsosGPos.y()-rhitGPos.y());
347  //double pullGPZ_tr = (tsosGPos.z()-rhitGPos.z());
348 
349  LogTrace("TestHits") << "tr" << std::endl;
350 
351  title.str("");
352  title << "PullGP_X_" << subdetId << "-" << layerId << "_tr";
353  hPullGP_X_tr[title.str()]->Fill(pullGPX_tr);
354  title.str("");
355  title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr";
356  hPullGP_Y_tr[title.str()]->Fill(pullGPY_tr);
357  title.str("");
358  title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr";
359  hPullGP_Z_tr[title.str()]->Fill(pullGPZ_tr);
360 
361  double pullGPX_ts = (tsosGPos.x() - shitGPos.x()) / sqrt(tsosGPEr.cxx());
362  double pullGPY_ts = (tsosGPos.y() - shitGPos.y()) / sqrt(tsosGPEr.cyy());
363  double pullGPZ_ts = (tsosGPos.z() - shitGPos.z()) / sqrt(tsosGPEr.czz());
364  //double pullGPX_ts = (tsosGPos.x()-shitGPos.x());
365  //double pullGPY_ts = (tsosGPos.y()-shitGPos.y());
366  //double pullGPZ_ts = (tsosGPos.z()-shitGPos.z());
367 
368  LogTrace("TestHits") << "ts1" << std::endl;
369 
370  title.str("");
371  title << "PullGP_X_" << subdetId << "-" << layerId << "_ts";
372  hPullGP_X_ts[title.str()]->Fill(pullGPX_ts);
373  title.str("");
374  title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts";
375  hPullGP_Y_ts[title.str()]->Fill(pullGPY_ts);
376  title.str("");
377  title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts";
378  hPullGP_Z_ts[title.str()]->Fill(pullGPZ_ts);
379 
380  double pullGMX_ts = (tsosGMom.x() - shitGMom.x()) / sqrt(tsosGMEr.cxx());
381  double pullGMY_ts = (tsosGMom.y() - shitGMom.y()) / sqrt(tsosGMEr.cyy());
382  double pullGMZ_ts = (tsosGMom.z() - shitGMom.z()) / sqrt(tsosGMEr.czz());
383  //double pullGMX_ts = (tsosGMom.x()-shitGMom.x());
384  //double pullGMY_ts = (tsosGMom.y()-shitGMom.y());
385  //double pullGMZ_ts = (tsosGMom.z()-shitGMom.z());
386 
387  LogTrace("TestHits") << "ts2" << std::endl;
388 
389  title.str("");
390  title << "PullGM_X_" << subdetId << "-" << layerId << "_ts";
391  hPullGM_X_ts[title.str()]->Fill(pullGMX_ts);
392  title.str("");
393  title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts";
394  hPullGM_Y_ts[title.str()]->Fill(pullGMY_ts);
395  title.str("");
396  title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts";
397  hPullGM_Z_ts[title.str()]->Fill(pullGMZ_ts);
398 
399  if (dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())) {
400  //mono
401  LogTrace("TestHits") << "MONO HIT" << std::endl;
402  auto m = dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())->monoHit();
403  CTTRHp tMonoHit = theBuilder->build(&m);
404  if (tMonoHit == nullptr)
405  continue;
406  vector<PSimHit> assMonoSimHits = hitAssociator.associateHit(*tMonoHit->hit());
407  if (assMonoSimHits.empty())
408  continue;
409  const PSimHit sMonoHit = *(assSimHits.begin());
410  const Surface* monoSurf = &(tMonoHit->det()->surface());
411  if (monoSurf == nullptr)
412  continue;
413  TSOS monoState = thePropagator->propagate(lastState, *monoSurf);
414  if (monoState.isValid() == 0)
415  continue;
416 
417  LocalVector monoShitLMom = sMonoHit.momentumAtEntry();
418  GlobalVector monoShitGMom = monoSurf->toGlobal(monoShitLMom);
419  LocalPoint monoShitLPos = sMonoHit.localPosition();
420  GlobalPoint monoShitGPos = monoSurf->toGlobal(monoShitLPos);
421 
422  GlobalVector monoTsosGMom = monoState.globalMomentum();
423  GlobalError monoTsosGMEr(monoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3, 3));
424  GlobalPoint monoTsosGPos = monoState.globalPosition();
425  GlobalError monoTsosGPEr = monoState.cartesianError().position();
426 
427  GlobalPoint monoRhitGPos = tMonoHit->globalPosition();
428  GlobalError monoRhitGPEr = tMonoHit->globalPositionError();
429 
430  double pullGPX_rs_mono = (monoRhitGPos.x() - monoShitGPos.x()) / sqrt(monoRhitGPEr.cxx());
431  double pullGPY_rs_mono = (monoRhitGPos.y() - monoShitGPos.y()) / sqrt(monoRhitGPEr.cyy());
432  double pullGPZ_rs_mono = (monoRhitGPos.z() - monoShitGPos.z()) / sqrt(monoRhitGPEr.czz());
433  //double pullGPX_rs_mono = (monoRhitGPos.x()-monoShitGPos.x());
434  //double pullGPY_rs_mono = (monoRhitGPos.y()-monoShitGPos.y());
435  //double pullGPZ_rs_mono = (monoRhitGPos.z()-monoShitGPos.z());
436 
437  title.str("");
438  title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_mono";
439  hPullGP_X_rs_mono[title.str()]->Fill(pullGPX_rs_mono);
440  title.str("");
441  title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_mono";
442  hPullGP_Y_rs_mono[title.str()]->Fill(pullGPY_rs_mono);
443  title.str("");
444  title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_mono";
445  hPullGP_Z_rs_mono[title.str()]->Fill(pullGPZ_rs_mono);
446 
447  double pullGPX_tr_mono = (monoTsosGPos.x() - monoRhitGPos.x()) / sqrt(monoTsosGPEr.cxx() + monoRhitGPEr.cxx());
448  double pullGPY_tr_mono = (monoTsosGPos.y() - monoRhitGPos.y()) / sqrt(monoTsosGPEr.cyy() + monoRhitGPEr.cyy());
449  double pullGPZ_tr_mono = (monoTsosGPos.z() - monoRhitGPos.z()) / sqrt(monoTsosGPEr.czz() + monoRhitGPEr.czz());
450  //double pullGPX_tr_mono = (monoTsosGPos.x()-monoRhitGPos.x());
451  //double pullGPY_tr_mono = (monoTsosGPos.y()-monoRhitGPos.y());
452  //double pullGPZ_tr_mono = (monoTsosGPos.z()-monoRhitGPos.z());
453 
454  title.str("");
455  title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_mono";
456  hPullGP_X_tr_mono[title.str()]->Fill(pullGPX_tr_mono);
457  title.str("");
458  title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_mono";
459  hPullGP_Y_tr_mono[title.str()]->Fill(pullGPY_tr_mono);
460  title.str("");
461  title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_mono";
462  hPullGP_Z_tr_mono[title.str()]->Fill(pullGPZ_tr_mono);
463 
464  double pullGPX_ts_mono = (monoTsosGPos.x() - monoShitGPos.x()) / sqrt(monoTsosGPEr.cxx());
465  double pullGPY_ts_mono = (monoTsosGPos.y() - monoShitGPos.y()) / sqrt(monoTsosGPEr.cyy());
466  double pullGPZ_ts_mono = (monoTsosGPos.z() - monoShitGPos.z()) / sqrt(monoTsosGPEr.czz());
467  //double pullGPX_ts_mono = (monoTsosGPos.x()-monoShitGPos.x());
468  //double pullGPY_ts_mono = (monoTsosGPos.y()-monoShitGPos.y());
469  //double pullGPZ_ts_mono = (monoTsosGPos.z()-monoShitGPos.z());
470 
471  title.str("");
472  title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_mono";
473  hPullGP_X_ts_mono[title.str()]->Fill(pullGPX_ts_mono);
474  title.str("");
475  title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_mono";
476  hPullGP_Y_ts_mono[title.str()]->Fill(pullGPY_ts_mono);
477  title.str("");
478  title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_mono";
479  hPullGP_Z_ts_mono[title.str()]->Fill(pullGPZ_ts_mono);
480 
481  double pullGMX_ts_mono = (monoTsosGMom.x() - monoShitGMom.x()) / sqrt(monoTsosGMEr.cxx());
482  double pullGMY_ts_mono = (monoTsosGMom.y() - monoShitGMom.y()) / sqrt(monoTsosGMEr.cyy());
483  double pullGMZ_ts_mono = (monoTsosGMom.z() - monoShitGMom.z()) / sqrt(monoTsosGMEr.czz());
484  //double pullGMX_ts_mono = (monoTsosGMom.x()-monoShitGMom.x());
485  //double pullGMY_ts_mono = (monoTsosGMom.y()-monoShitGMom.y());
486  //double pullGMZ_ts_mono = (monoTsosGMom.z()-monoShitGMom.z());
487 
488  title.str("");
489  title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_mono";
490  hPullGM_X_ts_mono[title.str()]->Fill(pullGMX_ts_mono);
491  title.str("");
492  title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_mono";
493  hPullGM_Y_ts_mono[title.str()]->Fill(pullGMY_ts_mono);
494  title.str("");
495  title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_mono";
496  hPullGM_Z_ts_mono[title.str()]->Fill(pullGMZ_ts_mono);
497 
498  //stereo
499  LogTrace("TestHits") << "STEREO HIT" << std::endl;
500  auto s = dynamic_cast<const SiStripMatchedRecHit2D*>((rhit)->hit())->stereoHit();
501  CTTRHp tStereoHit = theBuilder->build(&s);
502  if (tStereoHit == nullptr)
503  continue;
504  vector<PSimHit> assStereoSimHits = hitAssociator.associateHit(*tStereoHit->hit());
505  if (assStereoSimHits.empty())
506  continue;
507  const PSimHit sStereoHit = *(assSimHits.begin());
508  const Surface* stereoSurf = &(tStereoHit->det()->surface());
509  if (stereoSurf == nullptr)
510  continue;
511  TSOS stereoState = thePropagator->propagate(lastState, *stereoSurf);
512  if (stereoState.isValid() == 0)
513  continue;
514 
515  LocalVector stereoShitLMom = sStereoHit.momentumAtEntry();
516  GlobalVector stereoShitGMom = stereoSurf->toGlobal(stereoShitLMom);
517  LocalPoint stereoShitLPos = sStereoHit.localPosition();
518  GlobalPoint stereoShitGPos = stereoSurf->toGlobal(stereoShitLPos);
519 
520  GlobalVector stereoTsosGMom = stereoState.globalMomentum();
521  GlobalError stereoTsosGMEr(stereoState.cartesianError().matrix().Sub<AlgebraicSymMatrix33>(3, 3));
522  GlobalPoint stereoTsosGPos = stereoState.globalPosition();
523  GlobalError stereoTsosGPEr = stereoState.cartesianError().position();
524 
525  GlobalPoint stereoRhitGPos = tStereoHit->globalPosition();
526  GlobalError stereoRhitGPEr = tStereoHit->globalPositionError();
527 
528  double pullGPX_rs_stereo = (stereoRhitGPos.x() - stereoShitGPos.x()) / sqrt(stereoRhitGPEr.cxx());
529  double pullGPY_rs_stereo = (stereoRhitGPos.y() - stereoShitGPos.y()) / sqrt(stereoRhitGPEr.cyy());
530  double pullGPZ_rs_stereo = (stereoRhitGPos.z() - stereoShitGPos.z()) / sqrt(stereoRhitGPEr.czz());
531  //double pullGPX_rs_stereo = (stereoRhitGPos.x()-stereoShitGPos.x());
532  //double pullGPY_rs_stereo = (stereoRhitGPos.y()-stereoShitGPos.y());
533  //double pullGPZ_rs_stereo = (stereoRhitGPos.z()-stereoShitGPos.z());
534 
535  title.str("");
536  title << "PullGP_X_" << subdetId << "-" << layerId << "_rs_stereo";
537  hPullGP_X_rs_stereo[title.str()]->Fill(pullGPX_rs_stereo);
538  title.str("");
539  title << "PullGP_Y_" << subdetId << "-" << layerId << "_rs_stereo";
540  hPullGP_Y_rs_stereo[title.str()]->Fill(pullGPY_rs_stereo);
541  title.str("");
542  title << "PullGP_Z_" << subdetId << "-" << layerId << "_rs_stereo";
543  hPullGP_Z_rs_stereo[title.str()]->Fill(pullGPZ_rs_stereo);
544 
545  double pullGPX_tr_stereo =
546  (stereoTsosGPos.x() - stereoRhitGPos.x()) / sqrt(stereoTsosGPEr.cxx() + stereoRhitGPEr.cxx());
547  double pullGPY_tr_stereo =
548  (stereoTsosGPos.y() - stereoRhitGPos.y()) / sqrt(stereoTsosGPEr.cyy() + stereoRhitGPEr.cyy());
549  double pullGPZ_tr_stereo =
550  (stereoTsosGPos.z() - stereoRhitGPos.z()) / sqrt(stereoTsosGPEr.czz() + stereoRhitGPEr.czz());
551  //double pullGPX_tr_stereo = (stereoTsosGPos.x()-stereoRhitGPos.x());
552  //double pullGPY_tr_stereo = (stereoTsosGPos.y()-stereoRhitGPos.y());
553  //double pullGPZ_tr_stereo = (stereoTsosGPos.z()-stereoRhitGPos.z());
554 
555  title.str("");
556  title << "PullGP_X_" << subdetId << "-" << layerId << "_tr_stereo";
557  hPullGP_X_tr_stereo[title.str()]->Fill(pullGPX_tr_stereo);
558  title.str("");
559  title << "PullGP_Y_" << subdetId << "-" << layerId << "_tr_stereo";
560  hPullGP_Y_tr_stereo[title.str()]->Fill(pullGPY_tr_stereo);
561  title.str("");
562  title << "PullGP_Z_" << subdetId << "-" << layerId << "_tr_stereo";
563  hPullGP_Z_tr_stereo[title.str()]->Fill(pullGPZ_tr_stereo);
564 
565  double pullGPX_ts_stereo = (stereoTsosGPos.x() - stereoShitGPos.x()) / sqrt(stereoTsosGPEr.cxx());
566  double pullGPY_ts_stereo = (stereoTsosGPos.y() - stereoShitGPos.y()) / sqrt(stereoTsosGPEr.cyy());
567  double pullGPZ_ts_stereo = (stereoTsosGPos.z() - stereoShitGPos.z()) / sqrt(stereoTsosGPEr.czz());
568  //double pullGPX_ts_stereo = (stereoTsosGPos.x()-stereoShitGPos.x());
569  //double pullGPY_ts_stereo = (stereoTsosGPos.y()-stereoShitGPos.y());
570  //double pullGPZ_ts_stereo = (stereoTsosGPos.z()-stereoShitGPos.z());
571 
572  title.str("");
573  title << "PullGP_X_" << subdetId << "-" << layerId << "_ts_stereo";
574  hPullGP_X_ts_stereo[title.str()]->Fill(pullGPX_ts_stereo);
575  title.str("");
576  title << "PullGP_Y_" << subdetId << "-" << layerId << "_ts_stereo";
577  hPullGP_Y_ts_stereo[title.str()]->Fill(pullGPY_ts_stereo);
578  title.str("");
579  title << "PullGP_Z_" << subdetId << "-" << layerId << "_ts_stereo";
580  hPullGP_Z_ts_stereo[title.str()]->Fill(pullGPZ_ts_stereo);
581 
582  double pullGMX_ts_stereo = (stereoTsosGMom.x() - stereoShitGMom.x()) / sqrt(stereoTsosGMEr.cxx());
583  double pullGMY_ts_stereo = (stereoTsosGMom.y() - stereoShitGMom.y()) / sqrt(stereoTsosGMEr.cyy());
584  double pullGMZ_ts_stereo = (stereoTsosGMom.z() - stereoShitGMom.z()) / sqrt(stereoTsosGMEr.czz());
585  //double pullGMX_ts_stereo = (stereoTsosGMom.x()-stereoShitGMom.x());
586  //double pullGMY_ts_stereo = (stereoTsosGMom.y()-stereoShitGMom.y());
587  //double pullGMZ_ts_stereo = (stereoTsosGMom.z()-stereoShitGMom.z());
588 
589  title.str("");
590  title << "PullGM_X_" << subdetId << "-" << layerId << "_ts_stereo";
591  hPullGM_X_ts_stereo[title.str()]->Fill(pullGMX_ts_stereo);
592  title.str("");
593  title << "PullGM_Y_" << subdetId << "-" << layerId << "_ts_stereo";
594  hPullGM_Y_ts_stereo[title.str()]->Fill(pullGMY_ts_stereo);
595  title.str("");
596  title << "PullGM_Z_" << subdetId << "-" << layerId << "_ts_stereo";
597  hPullGM_Z_ts_stereo[title.str()]->Fill(pullGMZ_ts_stereo);
598  }
599  lastState = updatedState;
600  }
601  LogTrace("TestHits") << "traj chi2=" << tchi2;
602  LogTrace("TestHits") << "track chi2=" << result[0].chiSquared();
603  }
604  LogTrace("TestHits") << "end of event" << std::endl;
605 }

References TransientTrackingRecHitBuilder::build(), TrajectoryStateOnSurface::cartesianError(), heavyFlavorValidationHarvestingSequence_cff::combiner, GlobalErrorBase< T, ErrorWeightType >::cxx(), GlobalErrorBase< T, ErrorWeightType >::cyy(), GlobalErrorBase< T, ErrorWeightType >::czz(), dumpMFGeometry_cfg::delta, DetId::det(), PTrajectoryStateOnDet::detId(), edm::EventSetup::get(), get, TrajectoryStateOnSurface::globalMomentum(), TrajectoryStateOnSurface::globalPosition(), hChi2Increment, hClsize_vs_Chi2, TrackQuality_cff::hitAssociator, hfClusterShapes_cfi::hits, hProcess_vs_Chi2, hPullGM_X_ts, hPullGM_X_ts_mono, hPullGM_X_ts_stereo, hPullGM_Y_ts, hPullGM_Y_ts_mono, hPullGM_Y_ts_stereo, hPullGM_Z_ts, hPullGM_Z_ts_mono, hPullGM_Z_ts_stereo, hPullGP_X_rs, hPullGP_X_rs_mono, hPullGP_X_rs_stereo, hPullGP_X_tr, hPullGP_X_tr_mono, hPullGP_X_tr_stereo, hPullGP_X_ts, hPullGP_X_ts_mono, hPullGP_X_ts_stereo, hPullGP_Y_rs, hPullGP_Y_rs_mono, hPullGP_Y_rs_stereo, hPullGP_Y_tr, hPullGP_Y_tr_mono, hPullGP_Y_tr_stereo, hPullGP_Y_ts, hPullGP_Y_ts_mono, hPullGP_Y_ts_stereo, hPullGP_Z_rs, hPullGP_Z_rs_mono, hPullGP_Z_rs_stereo, hPullGP_Z_tr, hPullGP_Z_tr_mono, hPullGP_Z_tr_stereo, hPullGP_Z_ts, hPullGP_Z_ts_mono, hPullGP_Z_ts_stereo, hTotChi2Increment, mps_fire::i, TrackerGeometry::idToDet(), iEvent, TrajectoryStateOnSurface::isValid(), TrackerTopology::layer(), PSimHit::localPosition(), LogTrace, visualization-live-secondInstance_cfg::m, mag(), CartesianTrajectoryError::matrix(), maxeta, mineta, PSimHit::momentumAtEntry(), GloballyPositioned< T >::position(), CartesianTrajectoryError::position(), PSimHit::processType(), edm::ESHandle< T >::product(), projectHit(), Propagator::propagate(), TrackCandidate::recHits(), mps_fire::result, alignCSCRings::s, TrackCandidate::seed(), mathSSE::sqrt(), srcName, GeomDet::surface(), theBuilder, theG, theMF, thePropagator, theTCCollection, title, Surface::toGlobal(), trackerHitAssociatorConfig_, TrackCandidate::trajectoryStateOnDet(), trajectoryStateTransform::transientState(), TrackInfoProducer_cfi::updatedState, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

◆ beginRun()

void TestHits::beginRun ( edm::Run const &  run,
const edm::EventSetup &  iSetup 
)
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 36 of file TestHits.cc.

36  {
37  iSetup.get<TrackerDigiGeometryRecord>().get(theG);
38  iSetup.get<IdealMagneticFieldRecord>().get(theMF);
39  iSetup.get<TrackingComponentsRecord>().get(propagatorName, thePropagator);
40  iSetup.get<TransientRecHitRecord>().get(builderName, theBuilder);
41  iSetup.get<TrajectoryFitter::Record>().get(fname, fit);
42 
43  file = new TFile("testhits.root", "recreate");
44  for (int i = 0; i != 6; i++)
45  for (int j = 0; j != 9; j++) {
46  if (i == 0 && j > 2)
47  break;
48  if (i == 1 && j > 1)
49  break;
50  if (i == 2 && j > 3)
51  break;
52  if (i == 3 && j > 2)
53  break;
54  if (i == 4 && j > 5)
55  break;
56  if (i == 5 && j > 8)
57  break;
58  title.str("");
59  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts";
60  hPullGP_X_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
61  title.str("");
62  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts";
63  hPullGP_Y_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
64  title.str("");
65  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts";
66  hPullGP_Z_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
67  title.str("");
68  title << "Chi2Increment_" << i + 1 << "-" << j + 1;
69  hChi2Increment[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, 0, 100);
70 
71  title.str("");
72  title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts";
73  hPullGM_X_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
74  title.str("");
75  title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts";
76  hPullGM_Y_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
77  title.str("");
78  title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts";
79  hPullGM_Z_ts[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
80 
81  title.str("");
82  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr";
83  hPullGP_X_tr[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
84  title.str("");
85  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr";
86  hPullGP_Y_tr[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
87  title.str("");
88  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr";
89  hPullGP_Z_tr[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
90 
91  title.str("");
92  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs";
93  hPullGP_X_rs[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
94  title.str("");
95  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs";
96  hPullGP_Y_rs[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
97  title.str("");
98  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs";
99  hPullGP_Z_rs[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
100 
101  if (((i == 2 || i == 4) && (j == 0 || j == 1)) || (i == 3 || i == 5)) {
102  //mono
103  title.str("");
104  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts_mono";
105  hPullGP_X_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
106  title.str("");
107  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts_mono";
108  hPullGP_Y_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
109  title.str("");
110  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts_mono";
111  hPullGP_Z_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
112 
113  title.str("");
114  title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts_mono";
115  hPullGM_X_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
116  title.str("");
117  title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts_mono";
118  hPullGM_Y_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
119  title.str("");
120  title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts_mono";
121  hPullGM_Z_ts_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
122 
123  title.str("");
124  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr_mono";
125  hPullGP_X_tr_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
126  title.str("");
127  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr_mono";
128  hPullGP_Y_tr_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
129  title.str("");
130  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr_mono";
131  hPullGP_Z_tr_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
132 
133  title.str("");
134  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs_mono";
135  hPullGP_X_rs_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
136  title.str("");
137  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs_mono";
138  hPullGP_Y_rs_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
139  title.str("");
140  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs_mono";
141  hPullGP_Z_rs_mono[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
142 
143  //stereo
144  title.str("");
145  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts_stereo";
146  hPullGP_X_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
147  title.str("");
148  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts_stereo";
149  hPullGP_Y_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
150  title.str("");
151  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts_stereo";
152  hPullGP_Z_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
153 
154  title.str("");
155  title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts_stereo";
156  hPullGM_X_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
157  title.str("");
158  title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts_stereo";
159  hPullGM_Y_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
160  title.str("");
161  title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts_stereo";
162  hPullGM_Z_ts_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
163 
164  title.str("");
165  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr_stereo";
166  hPullGP_X_tr_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
167  title.str("");
168  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr_stereo";
169  hPullGP_Y_tr_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
170  title.str("");
171  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr_stereo";
172  hPullGP_Z_tr_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
173 
174  title.str("");
175  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs_stereo";
176  hPullGP_X_rs_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
177  title.str("");
178  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs_stereo";
179  hPullGP_Y_rs_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
180  title.str("");
181  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs_stereo";
182  hPullGP_Z_rs_stereo[title.str()] = new TH1F(title.str().c_str(), title.str().c_str(), 1000, -50, 50);
183  }
184  }
185  hTotChi2Increment = new TH1F("TotChi2Increment", "TotChi2Increment", 1000, 0, 100);
186  hProcess_vs_Chi2 = new TH2F("Process_vs_Chi2", "Process_vs_Chi2", 1000, 0, 100, 17, -0.5, 16.5);
187  hClsize_vs_Chi2 = new TH2F("Clsize_vs_Chi2", "Clsize_vs_Chi2", 1000, 0, 100, 17, -0.5, 16.5);
188 }

References builderName, file, fname, edm::EventSetup::get(), get, hChi2Increment, hClsize_vs_Chi2, hProcess_vs_Chi2, hPullGM_X_ts, hPullGM_X_ts_mono, hPullGM_X_ts_stereo, hPullGM_Y_ts, hPullGM_Y_ts_mono, hPullGM_Y_ts_stereo, hPullGM_Z_ts, hPullGM_Z_ts_mono, hPullGM_Z_ts_stereo, hPullGP_X_rs, hPullGP_X_rs_mono, hPullGP_X_rs_stereo, hPullGP_X_tr, hPullGP_X_tr_mono, hPullGP_X_tr_stereo, hPullGP_X_ts, hPullGP_X_ts_mono, hPullGP_X_ts_stereo, hPullGP_Y_rs, hPullGP_Y_rs_mono, hPullGP_Y_rs_stereo, hPullGP_Y_tr, hPullGP_Y_tr_mono, hPullGP_Y_tr_stereo, hPullGP_Y_ts, hPullGP_Y_ts_mono, hPullGP_Y_ts_stereo, hPullGP_Z_rs, hPullGP_Z_rs_mono, hPullGP_Z_rs_stereo, hPullGP_Z_tr, hPullGP_Z_tr_mono, hPullGP_Z_tr_stereo, hPullGP_Z_ts, hPullGP_Z_ts_mono, hPullGP_Z_ts_stereo, hTotChi2Increment, mps_fire::i, dqmiolumiharvest::j, propagatorName, theBuilder, theG, theMF, thePropagator, and title.

◆ endJob()

void TestHits::endJob ( void  )
overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 947 of file TestHits.cc.

947  {
948  //file->Write();
949  TDirectory* chi2i = file->mkdir("Chi2_Increment");
950 
951  TDirectory* gp_ts = file->mkdir("GP_TSOS-SimHit");
952  TDirectory* gm_ts = file->mkdir("GM_TSOS-SimHit");
953  TDirectory* gp_tr = file->mkdir("GP_TSOS-RecHit");
954  TDirectory* gp_rs = file->mkdir("GP_RecHit-SimHit");
955 
956  TDirectory* gp_tsx = gp_ts->mkdir("X");
957  TDirectory* gp_tsy = gp_ts->mkdir("Y");
958  TDirectory* gp_tsz = gp_ts->mkdir("Z");
959  TDirectory* gm_tsx = gm_ts->mkdir("X");
960  TDirectory* gm_tsy = gm_ts->mkdir("Y");
961  TDirectory* gm_tsz = gm_ts->mkdir("Z");
962  TDirectory* gp_trx = gp_tr->mkdir("X");
963  TDirectory* gp_try = gp_tr->mkdir("Y");
964  TDirectory* gp_trz = gp_tr->mkdir("Z");
965  TDirectory* gp_rsx = gp_rs->mkdir("X");
966  TDirectory* gp_rsy = gp_rs->mkdir("Y");
967  TDirectory* gp_rsz = gp_rs->mkdir("Z");
968 
969  TDirectory* gp_tsx_mono = gp_ts->mkdir("MONOX");
970  TDirectory* gp_tsy_mono = gp_ts->mkdir("MONOY");
971  TDirectory* gp_tsz_mono = gp_ts->mkdir("MONOZ");
972  TDirectory* gm_tsx_mono = gm_ts->mkdir("MONOX");
973  TDirectory* gm_tsy_mono = gm_ts->mkdir("MONOY");
974  TDirectory* gm_tsz_mono = gm_ts->mkdir("MONOZ");
975  TDirectory* gp_trx_mono = gp_tr->mkdir("MONOX");
976  TDirectory* gp_try_mono = gp_tr->mkdir("MONOY");
977  TDirectory* gp_trz_mono = gp_tr->mkdir("MONOZ");
978  TDirectory* gp_rsx_mono = gp_rs->mkdir("MONOX");
979  TDirectory* gp_rsy_mono = gp_rs->mkdir("MONOY");
980  TDirectory* gp_rsz_mono = gp_rs->mkdir("MONOZ");
981 
982  TDirectory* gp_tsx_stereo = gp_ts->mkdir("STEREOX");
983  TDirectory* gp_tsy_stereo = gp_ts->mkdir("STEREOY");
984  TDirectory* gp_tsz_stereo = gp_ts->mkdir("STEREOZ");
985  TDirectory* gm_tsx_stereo = gm_ts->mkdir("STEREOX");
986  TDirectory* gm_tsy_stereo = gm_ts->mkdir("STEREOY");
987  TDirectory* gm_tsz_stereo = gm_ts->mkdir("STEREOZ");
988  TDirectory* gp_trx_stereo = gp_tr->mkdir("STEREOX");
989  TDirectory* gp_try_stereo = gp_tr->mkdir("STEREOY");
990  TDirectory* gp_trz_stereo = gp_tr->mkdir("STEREOZ");
991  TDirectory* gp_rsx_stereo = gp_rs->mkdir("STEREOX");
992  TDirectory* gp_rsy_stereo = gp_rs->mkdir("STEREOY");
993  TDirectory* gp_rsz_stereo = gp_rs->mkdir("STEREOZ");
994 
995  chi2i->cd();
996  hTotChi2Increment->Write();
997  hProcess_vs_Chi2->Write();
998  hClsize_vs_Chi2->Write();
999  for (int i = 0; i != 6; i++)
1000  for (int j = 0; j != 9; j++) {
1001  if (i == 0 && j > 2)
1002  break;
1003  if (i == 1 && j > 1)
1004  break;
1005  if (i == 2 && j > 3)
1006  break;
1007  if (i == 3 && j > 2)
1008  break;
1009  if (i == 4 && j > 5)
1010  break;
1011  if (i == 5 && j > 8)
1012  break;
1013  chi2i->cd();
1014  title.str("");
1015  title << "Chi2Increment_" << i + 1 << "-" << j + 1;
1016  hChi2Increment[title.str()]->Write();
1017 
1018  gp_ts->cd();
1019  gp_tsx->cd();
1020  title.str("");
1021  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts";
1022  hPullGP_X_ts[title.str()]->Write();
1023  gp_tsy->cd();
1024  title.str("");
1025  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts";
1026  hPullGP_Y_ts[title.str()]->Write();
1027  gp_tsz->cd();
1028  title.str("");
1029  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts";
1030  hPullGP_Z_ts[title.str()]->Write();
1031 
1032  gm_ts->cd();
1033  gm_tsx->cd();
1034  title.str("");
1035  title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts";
1036  hPullGM_X_ts[title.str()]->Write();
1037  gm_tsy->cd();
1038  title.str("");
1039  title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts";
1040  hPullGM_Y_ts[title.str()]->Write();
1041  gm_tsz->cd();
1042  title.str("");
1043  title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts";
1044  hPullGM_Z_ts[title.str()]->Write();
1045 
1046  gp_tr->cd();
1047  gp_trx->cd();
1048  title.str("");
1049  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr";
1050  hPullGP_X_tr[title.str()]->Write();
1051  gp_try->cd();
1052  title.str("");
1053  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr";
1054  hPullGP_Y_tr[title.str()]->Write();
1055  gp_trz->cd();
1056  title.str("");
1057  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr";
1058  hPullGP_Z_tr[title.str()]->Write();
1059 
1060  gp_rs->cd();
1061  gp_rsx->cd();
1062  title.str("");
1063  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs";
1064  hPullGP_X_rs[title.str()]->Write();
1065  gp_rsy->cd();
1066  title.str("");
1067  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs";
1068  hPullGP_Y_rs[title.str()]->Write();
1069  gp_rsz->cd();
1070  title.str("");
1071  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs";
1072  hPullGP_Z_rs[title.str()]->Write();
1073 
1074  if (((i == 2 || i == 4) && (j == 0 || j == 1)) || (i == 3 || i == 5)) {
1075  //mono
1076  gp_ts->cd();
1077  gp_tsx_mono->cd();
1078  title.str("");
1079  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts_mono";
1080  hPullGP_X_ts_mono[title.str()]->Write();
1081  gp_tsy_mono->cd();
1082  title.str("");
1083  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts_mono";
1084  hPullGP_Y_ts_mono[title.str()]->Write();
1085  gp_tsz_mono->cd();
1086  title.str("");
1087  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts_mono";
1088  hPullGP_Z_ts_mono[title.str()]->Write();
1089 
1090  gm_ts->cd();
1091  gm_tsx_mono->cd();
1092  title.str("");
1093  title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts_mono";
1094  hPullGM_X_ts_mono[title.str()]->Write();
1095  gm_tsy_mono->cd();
1096  title.str("");
1097  title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts_mono";
1098  hPullGM_Y_ts_mono[title.str()]->Write();
1099  gm_tsz_mono->cd();
1100  title.str("");
1101  title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts_mono";
1102  hPullGM_Z_ts_mono[title.str()]->Write();
1103 
1104  gp_tr->cd();
1105  gp_trx_mono->cd();
1106  title.str("");
1107  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr_mono";
1108  hPullGP_X_tr_mono[title.str()]->Write();
1109  gp_try_mono->cd();
1110  title.str("");
1111  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr_mono";
1112  hPullGP_Y_tr_mono[title.str()]->Write();
1113  gp_trz_mono->cd();
1114  title.str("");
1115  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr_mono";
1116  hPullGP_Z_tr_mono[title.str()]->Write();
1117 
1118  gp_rs->cd();
1119  gp_rsx_mono->cd();
1120  title.str("");
1121  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs_mono";
1122  hPullGP_X_rs_mono[title.str()]->Write();
1123  gp_rsy_mono->cd();
1124  title.str("");
1125  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs_mono";
1126  hPullGP_Y_rs_mono[title.str()]->Write();
1127  gp_rsz_mono->cd();
1128  title.str("");
1129  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs_mono";
1130  hPullGP_Z_rs_mono[title.str()]->Write();
1131 
1132  //stereo
1133  gp_ts->cd();
1134  gp_tsx_stereo->cd();
1135  title.str("");
1136  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_ts_stereo";
1137  hPullGP_X_ts_stereo[title.str()]->Write();
1138  gp_tsy_stereo->cd();
1139  title.str("");
1140  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_ts_stereo";
1141  hPullGP_Y_ts_stereo[title.str()]->Write();
1142  gp_tsz_stereo->cd();
1143  title.str("");
1144  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_ts_stereo";
1145  hPullGP_Z_ts_stereo[title.str()]->Write();
1146 
1147  gm_ts->cd();
1148  gm_tsx_stereo->cd();
1149  title.str("");
1150  title << "PullGM_X_" << i + 1 << "-" << j + 1 << "_ts_stereo";
1151  hPullGM_X_ts_stereo[title.str()]->Write();
1152  gm_tsy_stereo->cd();
1153  title.str("");
1154  title << "PullGM_Y_" << i + 1 << "-" << j + 1 << "_ts_stereo";
1155  hPullGM_Y_ts_stereo[title.str()]->Write();
1156  gm_tsz_stereo->cd();
1157  title.str("");
1158  title << "PullGM_Z_" << i + 1 << "-" << j + 1 << "_ts_stereo";
1159  hPullGM_Z_ts_stereo[title.str()]->Write();
1160 
1161  gp_tr->cd();
1162  gp_trx_stereo->cd();
1163  title.str("");
1164  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_tr_stereo";
1165  hPullGP_X_tr_stereo[title.str()]->Write();
1166  gp_try_stereo->cd();
1167  title.str("");
1168  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_tr_stereo";
1169  hPullGP_Y_tr_stereo[title.str()]->Write();
1170  gp_trz_stereo->cd();
1171  title.str("");
1172  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_tr_stereo";
1173  hPullGP_Z_tr_stereo[title.str()]->Write();
1174 
1175  gp_rs->cd();
1176  gp_rsx_stereo->cd();
1177  title.str("");
1178  title << "PullGP_X_" << i + 1 << "-" << j + 1 << "_rs_stereo";
1179  hPullGP_X_rs_stereo[title.str()]->Write();
1180  gp_rsy_stereo->cd();
1181  title.str("");
1182  title << "PullGP_Y_" << i + 1 << "-" << j + 1 << "_rs_stereo";
1183  hPullGP_Y_rs_stereo[title.str()]->Write();
1184  gp_rsz_stereo->cd();
1185  title.str("");
1186  title << "PullGP_Z_" << i + 1 << "-" << j + 1 << "_rs_stereo";
1187  hPullGP_Z_rs_stereo[title.str()]->Write();
1188  }
1189  }
1190 
1191  file->Close();
1192 }

References file, hChi2Increment, hClsize_vs_Chi2, hProcess_vs_Chi2, hPullGM_X_ts, hPullGM_X_ts_mono, hPullGM_X_ts_stereo, hPullGM_Y_ts, hPullGM_Y_ts_mono, hPullGM_Y_ts_stereo, hPullGM_Z_ts, hPullGM_Z_ts_mono, hPullGM_Z_ts_stereo, hPullGP_X_rs, hPullGP_X_rs_mono, hPullGP_X_rs_stereo, hPullGP_X_tr, hPullGP_X_tr_mono, hPullGP_X_tr_stereo, hPullGP_X_ts, hPullGP_X_ts_mono, hPullGP_X_ts_stereo, hPullGP_Y_rs, hPullGP_Y_rs_mono, hPullGP_Y_rs_stereo, hPullGP_Y_tr, hPullGP_Y_tr_mono, hPullGP_Y_tr_stereo, hPullGP_Y_ts, hPullGP_Y_ts_mono, hPullGP_Y_ts_stereo, hPullGP_Z_rs, hPullGP_Z_rs_mono, hPullGP_Z_rs_stereo, hPullGP_Z_tr, hPullGP_Z_tr_mono, hPullGP_Z_tr_stereo, hPullGP_Z_ts, hPullGP_Z_ts_mono, hPullGP_Z_ts_stereo, hTotChi2Increment, mps_fire::i, dqmiolumiharvest::j, and title.

◆ projectHit()

std::pair< LocalPoint, LocalVector > TestHits::projectHit ( const PSimHit &  hit,
const StripGeomDetUnit *  stripDet,
const BoundPlane &  plane 
)
private

Definition at line 1196 of file TestHits.cc.

1198  {
1199  const StripTopology& topol = stripDet->specificTopology();
1200  GlobalPoint globalpos = stripDet->surface().toGlobal(hit.localPosition());
1201  LocalPoint localHit = plane.toLocal(globalpos);
1202  //track direction
1203  LocalVector locdir = hit.localDirection();
1204  //rotate track in new frame
1205 
1206  GlobalVector globaldir = stripDet->surface().toGlobal(locdir);
1207  LocalVector dir = plane.toLocal(globaldir);
1208  float scale = -localHit.z() / dir.z();
1209 
1210  LocalPoint projectedPos = localHit + scale * dir;
1211 
1212  float selfAngle = topol.stripAngle(topol.strip(hit.localPosition()));
1213 
1214  LocalVector stripDir(sin(selfAngle), cos(selfAngle), 0); // vector along strip in hit frame
1215 
1216  LocalVector localStripDir(plane.toLocal(stripDet->surface().toGlobal(stripDir)));
1217 
1218  return std::pair<LocalPoint, LocalVector>(projectedPos, localStripDir);
1219 }

References funct::cos(), DeadROC_duringRun::dir, Scenarios_cff::scale, funct::sin(), StripGeomDetUnit::specificTopology(), StripTopology::strip(), StripTopology::stripAngle(), GeomDet::surface(), Surface::toGlobal(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

Member Data Documentation

◆ builderName

std::string TestHits::builderName
private

Definition at line 66 of file TestHits.h.

Referenced by beginRun(), and TestHits().

◆ file

TFile* TestHits::file
private

Definition at line 77 of file TestHits.h.

Referenced by beginRun(), and endJob().

◆ fit

edm::ESHandle<TrajectoryFitter> TestHits::fit
private

Definition at line 74 of file TestHits.h.

Referenced by trackingPlots.Iteration::modules().

◆ fname

std::string TestHits::fname
private

Definition at line 68 of file TestHits.h.

Referenced by beginRun(), and TestHits().

◆ hChi2Increment

std::map<std::string, TH1F*> TestHits::hChi2Increment
private

Definition at line 91 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hClsize_vs_Chi2

TH2F * TestHits::hClsize_vs_Chi2
private

Definition at line 93 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hProcess_vs_Chi2

TH2F* TestHits::hProcess_vs_Chi2
private

Definition at line 93 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGM_X_ts

std::map<std::string, TH1F*> TestHits::hPullGM_X_ts
private

Definition at line 82 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGM_X_ts_mono

std::map<std::string, TH1F*> TestHits::hPullGM_X_ts_mono
private

Definition at line 98 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGM_X_ts_stereo

std::map<std::string, TH1F*> TestHits::hPullGM_X_ts_stereo
private

Definition at line 111 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGM_Y_ts

std::map<std::string, TH1F*> TestHits::hPullGM_Y_ts
private

Definition at line 83 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGM_Y_ts_mono

std::map<std::string, TH1F*> TestHits::hPullGM_Y_ts_mono
private

Definition at line 99 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGM_Y_ts_stereo

std::map<std::string, TH1F*> TestHits::hPullGM_Y_ts_stereo
private

Definition at line 112 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGM_Z_ts

std::map<std::string, TH1F*> TestHits::hPullGM_Z_ts
private

Definition at line 84 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGM_Z_ts_mono

std::map<std::string, TH1F*> TestHits::hPullGM_Z_ts_mono
private

Definition at line 100 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGM_Z_ts_stereo

std::map<std::string, TH1F*> TestHits::hPullGM_Z_ts_stereo
private

Definition at line 113 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_X_rs

std::map<std::string, TH1F*> TestHits::hPullGP_X_rs
private

Definition at line 85 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_X_rs_mono

std::map<std::string, TH1F*> TestHits::hPullGP_X_rs_mono
private

Definition at line 101 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_X_rs_stereo

std::map<std::string, TH1F*> TestHits::hPullGP_X_rs_stereo
private

Definition at line 114 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_X_tr

std::map<std::string, TH1F*> TestHits::hPullGP_X_tr
private

Definition at line 88 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_X_tr_mono

std::map<std::string, TH1F*> TestHits::hPullGP_X_tr_mono
private

Definition at line 104 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_X_tr_stereo

std::map<std::string, TH1F*> TestHits::hPullGP_X_tr_stereo
private

Definition at line 117 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_X_ts

std::map<std::string, TH1F*> TestHits::hPullGP_X_ts
private

Definition at line 79 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_X_ts_mono

std::map<std::string, TH1F*> TestHits::hPullGP_X_ts_mono
private

Definition at line 95 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_X_ts_stereo

std::map<std::string, TH1F*> TestHits::hPullGP_X_ts_stereo
private

Definition at line 108 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Y_rs

std::map<std::string, TH1F*> TestHits::hPullGP_Y_rs
private

Definition at line 86 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Y_rs_mono

std::map<std::string, TH1F*> TestHits::hPullGP_Y_rs_mono
private

Definition at line 102 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Y_rs_stereo

std::map<std::string, TH1F*> TestHits::hPullGP_Y_rs_stereo
private

Definition at line 115 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Y_tr

std::map<std::string, TH1F*> TestHits::hPullGP_Y_tr
private

Definition at line 89 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Y_tr_mono

std::map<std::string, TH1F*> TestHits::hPullGP_Y_tr_mono
private

Definition at line 105 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Y_tr_stereo

std::map<std::string, TH1F*> TestHits::hPullGP_Y_tr_stereo
private

Definition at line 118 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Y_ts

std::map<std::string, TH1F*> TestHits::hPullGP_Y_ts
private

Definition at line 80 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Y_ts_mono

std::map<std::string, TH1F*> TestHits::hPullGP_Y_ts_mono
private

Definition at line 96 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Y_ts_stereo

std::map<std::string, TH1F*> TestHits::hPullGP_Y_ts_stereo
private

Definition at line 109 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Z_rs

std::map<std::string, TH1F*> TestHits::hPullGP_Z_rs
private

Definition at line 87 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Z_rs_mono

std::map<std::string, TH1F*> TestHits::hPullGP_Z_rs_mono
private

Definition at line 103 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Z_rs_stereo

std::map<std::string, TH1F*> TestHits::hPullGP_Z_rs_stereo
private

Definition at line 116 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Z_tr

std::map<std::string, TH1F*> TestHits::hPullGP_Z_tr
private

Definition at line 90 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Z_tr_mono

std::map<std::string, TH1F*> TestHits::hPullGP_Z_tr_mono
private

Definition at line 106 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Z_tr_stereo

std::map<std::string, TH1F*> TestHits::hPullGP_Z_tr_stereo
private

Definition at line 119 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Z_ts

std::map<std::string, TH1F*> TestHits::hPullGP_Z_ts
private

Definition at line 81 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Z_ts_mono

std::map<std::string, TH1F*> TestHits::hPullGP_Z_ts_mono
private

Definition at line 97 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hPullGP_Z_ts_stereo

std::map<std::string, TH1F*> TestHits::hPullGP_Z_ts_stereo
private

Definition at line 110 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ hTotChi2Increment

TH1F* TestHits::hTotChi2Increment
private

Definition at line 92 of file TestHits.h.

Referenced by analyze(), beginRun(), and endJob().

◆ maxeta

double TestHits::maxeta
private

Definition at line 63 of file TestHits.h.

Referenced by analyze(), and TestHits().

◆ mineta

double TestHits::mineta
private

Definition at line 63 of file TestHits.h.

Referenced by analyze(), and TestHits().

◆ propagatorName

std::string TestHits::propagatorName
private

Definition at line 65 of file TestHits.h.

Referenced by beginRun(), and TestHits().

◆ srcName

std::string TestHits::srcName
private

Definition at line 67 of file TestHits.h.

Referenced by analyze(), and TestHits().

◆ theBuilder

edm::ESHandle<TransientTrackingRecHitBuilder> TestHits::theBuilder
private

Definition at line 73 of file TestHits.h.

Referenced by analyze(), and beginRun().

◆ theG

edm::ESHandle<TrackerGeometry> TestHits::theG
private

Definition at line 70 of file TestHits.h.

Referenced by analyze(), and beginRun().

◆ theMF

edm::ESHandle<MagneticField> TestHits::theMF
private

Definition at line 71 of file TestHits.h.

Referenced by analyze(), and beginRun().

◆ thePropagator

edm::ESHandle<Propagator> TestHits::thePropagator
private

Definition at line 72 of file TestHits.h.

Referenced by analyze(), and beginRun().

◆ theTCCollection

edm::Handle<TrackCandidateCollection> TestHits::theTCCollection
private

Definition at line 75 of file TestHits.h.

Referenced by analyze().

◆ title

std::stringstream TestHits::title
private

Definition at line 78 of file TestHits.h.

Referenced by analyze(), beginRun(), endJob(), and Formatter.SimpleHTMLFormatter::headers().

◆ trackerHitAssociatorConfig_

TrackerHitAssociator::Config TestHits::trackerHitAssociatorConfig_
private

Definition at line 61 of file TestHits.h.

Referenced by analyze().

Vector3DBase< float, LocalTag >
TrackCandidate::trajectoryStateOnDet
PTrajectoryStateOnDet const & trajectoryStateOnDet() const
Definition: TrackCandidate.h:56
TrackerGeometry::idToDet
const TrackerGeomDet * idToDet(DetId) const override
Definition: TrackerGeometry.cc:193
TestHits::hPullGP_Y_rs
std::map< std::string, TH1F * > hPullGP_Y_rs
Definition: TestHits.h:86
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
AlgebraicSymMatrix33
ROOT::Math::SMatrix< double, 3, 3, ROOT::Math::MatRepSym< double, 3 > > AlgebraicSymMatrix33
Definition: AlgebraicROOTObjects.h:21
mps_fire.i
i
Definition: mps_fire.py:355
TestHits::hPullGP_Y_tr_stereo
std::map< std::string, TH1F * > hPullGP_Y_tr_stereo
Definition: TestHits.h:118
TestHits::theMF
edm::ESHandle< MagneticField > theMF
Definition: TestHits.h:71
TestHits::file
TFile * file
Definition: TestHits.h:77
TestHits::hPullGP_X_tr_mono
std::map< std::string, TH1F * > hPullGP_X_tr_mono
Definition: TestHits.h:104
TestHits::fname
std::string fname
Definition: TestHits.h:68
hfClusterShapes_cfi.hits
hits
Definition: hfClusterShapes_cfi.py:5
TrajectoryStateOnSurface::cartesianError
const CartesianTrajectoryError cartesianError() const
Definition: TrajectoryStateOnSurface.h:71
PV3DBase::x
T x() const
Definition: PV3DBase.h:59
DetId::det
constexpr Detector det() const
get the detector field from this detid
Definition: DetId.h:46
TestHits::trackerHitAssociatorConfig_
TrackerHitAssociator::Config trackerHitAssociatorConfig_
Definition: TestHits.h:61
TrajectoryStateOnSurface::globalPosition
GlobalPoint globalPosition() const
Definition: TrajectoryStateOnSurface.h:65
TransientRecHitRecord
Definition: TransientRecHitRecord.h:14
TestHits::hPullGP_Y_ts
std::map< std::string, TH1F * > hPullGP_Y_ts
Definition: TestHits.h:80
TestHits::hPullGP_Z_ts
std::map< std::string, TH1F * > hPullGP_Z_ts
Definition: TestHits.h:81
TestHits::hTotChi2Increment
TH1F * hTotChi2Increment
Definition: TestHits.h:92
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:71
Surface
Definition: Surface.h:36
TrackerTopology::layer
unsigned int layer(const DetId &id) const
Definition: TrackerTopology.cc:47
SiStripRecHit2D
Definition: SiStripRecHit2D.h:7
TestHits::title
std::stringstream title
Definition: TestHits.h:78
TestHits::builderName
std::string builderName
Definition: TestHits.h:66
TransientTrackingRecHitBuilder::build
virtual RecHitPointer build(const TrackingRecHit *p) const =0
build a tracking rechit from an existing rechit
TestHits::hPullGM_Y_ts
std::map< std::string, TH1F * > hPullGM_Y_ts
Definition: TestHits.h:83
TrackQuality_cff.hitAssociator
hitAssociator
Definition: TrackQuality_cff.py:10
TestHits::hPullGP_Z_tr_mono
std::map< std::string, TH1F * > hPullGP_Z_tr_mono
Definition: TestHits.h:106
SiPixelRecHit
Our base class.
Definition: SiPixelRecHit.h:23
TestHits::hClsize_vs_Chi2
TH2F * hClsize_vs_Chi2
Definition: TestHits.h:93
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
TestHits::hPullGP_Z_tr_stereo
std::map< std::string, TH1F * > hPullGP_Z_tr_stereo
Definition: TestHits.h:119
PV3DBase::z
T z() const
Definition: PV3DBase.h:61
TestHits::theTCCollection
edm::Handle< TrackCandidateCollection > theTCCollection
Definition: TestHits.h:75
TestHits::maxeta
double maxeta
Definition: TestHits.h:63
TestHits::hPullGM_X_ts_mono
std::map< std::string, TH1F * > hPullGM_X_ts_mono
Definition: TestHits.h:98
IdealMagneticFieldRecord
Definition: IdealMagneticFieldRecord.h:11
DetId
Definition: DetId.h:17
GeomDet::surface
const Plane & surface() const
The nominal surface of the GeomDet.
Definition: GeomDet.h:37
TrajectoryStateOnSurface
Definition: TrajectoryStateOnSurface.h:16
alignCSCRings.s
s
Definition: alignCSCRings.py:92
GlobalErrorBase::cyy
T cyy() const
Definition: GlobalErrorBase.h:101
TestHits::hPullGM_Y_ts_stereo
std::map< std::string, TH1F * > hPullGM_Y_ts_stereo
Definition: TestHits.h:112
PTrajectoryStateOnDet::detId
unsigned int detId() const
Definition: PTrajectoryStateOnDet.h:65
funct::cos
Cos< T >::type cos(const T &t)
Definition: Cos.h:22
CTTRHp
TransientTrackingRecHit::ConstRecHitPointer CTTRHp
Definition: CkfDebugger.h:39
PSimHit::processType
unsigned short processType() const
Definition: PSimHit.h:120
edm::EventSetup::get
T get() const
Definition: EventSetup.h:73
TestHits::hPullGM_Y_ts_mono
std::map< std::string, TH1F * > hPullGM_Y_ts_mono
Definition: TestHits.h:99
TrackingRecHit::RecHitContainer
std::vector< ConstRecHitPointer > RecHitContainer
Definition: TrackingRecHit.h:31
TestHits::mineta
double mineta
Definition: TestHits.h:63
visualization-live-secondInstance_cfg.m
m
Definition: visualization-live-secondInstance_cfg.py:72
TrackerDigiGeometryRecord
Definition: TrackerDigiGeometryRecord.h:15
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:19
TestHits::hPullGP_X_tr
std::map< std::string, TH1F * > hPullGP_X_tr
Definition: TestHits.h:88
TestHits::hPullGM_X_ts
std::map< std::string, TH1F * > hPullGM_X_ts
Definition: TestHits.h:82
Surface::toGlobal
GlobalPoint toGlobal(const Point2DBase< Scalar, LocalTag > lp) const
Definition: Surface.h:79
TestHits::hPullGP_X_ts_stereo
std::map< std::string, TH1F * > hPullGP_X_ts_stereo
Definition: TestHits.h:108
edm::ESHandle< TrackerTopology >
PSimHit::localPosition
Local3DPoint localPosition() const
Definition: PSimHit.h:52
TrajectoryFitterRecord
Definition: TrajectoryFitterRecord.h:12
TestHits::hPullGP_Z_ts_mono
std::map< std::string, TH1F * > hPullGP_Z_ts_mono
Definition: TestHits.h:97
Point3DBase< float, LocalTag >
TestHits::hPullGP_Z_ts_stereo
std::map< std::string, TH1F * > hPullGP_Z_ts_stereo
Definition: TestHits.h:110
GlobalErrorBase::cxx
T cxx() const
Definition: GlobalErrorBase.h:97
StripTopology::strip
virtual float strip(const LocalPoint &) const =0
GluedGeomDet
Definition: GluedGeomDet.h:7
edm::OwnVector::const_iterator
Definition: OwnVector.h:41
TrackCandidate::range
std::pair< const_iterator, const_iterator > range
Definition: TrackCandidate.h:27
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
Scenarios_cff.scale
scale
Definition: Scenarios_cff.py:2186
TrajectoryStateCombiner
Definition: TrajectoryStateCombiner.h:13
TestHits::hPullGP_Y_tr_mono
std::map< std::string, TH1F * > hPullGP_Y_tr_mono
Definition: TestHits.h:105
TestHits::theG
edm::ESHandle< TrackerGeometry > theG
Definition: TestHits.h:70
TrackCandidate
Definition: TrackCandidate.h:23
TrackCandidate::recHits
range recHits() const
Definition: TrackCandidate.h:58
TestHits::hPullGP_Y_tr
std::map< std::string, TH1F * > hPullGP_Y_tr
Definition: TestHits.h:89
TestHits::hPullGP_X_rs_mono
std::map< std::string, TH1F * > hPullGP_X_rs_mono
Definition: TestHits.h:101
TestHits::hChi2Increment
std::map< std::string, TH1F * > hChi2Increment
Definition: TestHits.h:91
TestHits::hPullGP_Z_rs_stereo
std::map< std::string, TH1F * > hPullGP_Z_rs_stereo
Definition: TestHits.h:116
dumpMFGeometry_cfg.delta
delta
Definition: dumpMFGeometry_cfg.py:25
TrackCandidate::seed
TrajectorySeed const & seed() const
Definition: TrackCandidate.h:60
PV3DBase::y
T y() const
Definition: PV3DBase.h:60
TrackingRecHit::ConstRecHitPointer
std::shared_ptr< TrackingRecHit const > ConstRecHitPointer
Definition: TrackingRecHit.h:25
StripTopology::stripAngle
virtual float stripAngle(float strip) const =0
iEvent
int iEvent
Definition: GenABIO.cc:224
trajectoryStateTransform::transientState
TrajectoryStateOnSurface transientState(const PTrajectoryStateOnDet &ts, const Surface *surface, const MagneticField *field)
Definition: TrajectoryStateTransform.cc:35
edm::LogVerbatim
Definition: MessageLogger.h:297
TestHits::projectHit
std::pair< LocalPoint, LocalVector > projectHit(const PSimHit &, const StripGeomDetUnit *, const BoundPlane &)
Definition: TestHits.cc:1196
TestHits::hPullGM_Z_ts_stereo
std::map< std::string, TH1F * > hPullGM_Z_ts_stereo
Definition: TestHits.h:113
GlobalErrorBase< double, ErrorMatrixTag >
TestHits::srcName
std::string srcName
Definition: TestHits.h:67
Propagator::propagate
TrajectoryStateOnSurface propagate(STA const &state, SUR const &surface) const
Definition: Propagator.h:50
get
#define get
TestHits::hPullGP_Y_rs_stereo
std::map< std::string, TH1F * > hPullGP_Y_rs_stereo
Definition: TestHits.h:115
CartesianTrajectoryError::matrix
const AlgebraicSymMatrix66 & matrix() const
Definition: CartesianTrajectoryError.h:28
TestHits::hPullGP_Z_rs
std::map< std::string, TH1F * > hPullGP_Z_rs
Definition: TestHits.h:87
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
TrackerHitAssociator
Definition: TrackerHitAssociator.h:55
GloballyPositioned::position
const PositionType & position() const
Definition: GloballyPositioned.h:36
TrajectoryStateOnSurface::globalMomentum
GlobalVector globalMomentum() const
Definition: TrajectoryStateOnSurface.h:66
mag
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
Definition: Basic3DVectorLD.h:127
TestHits::hPullGP_X_rs_stereo
std::map< std::string, TH1F * > hPullGP_X_rs_stereo
Definition: TestHits.h:114
TestHits::propagatorName
std::string propagatorName
Definition: TestHits.h:65
TrackInfoProducer_cfi.updatedState
updatedState
Definition: TrackInfoProducer_cfi.py:6
PSimHit::momentumAtEntry
LocalVector momentumAtEntry() const
The momentum of the track that produced the hit, at entry point.
Definition: PSimHit.h:55
TestHits::thePropagator
edm::ESHandle< Propagator > thePropagator
Definition: TestHits.h:72
BoundPlane
TestHits::hPullGM_Z_ts_mono
std::map< std::string, TH1F * > hPullGM_Z_ts_mono
Definition: TestHits.h:100
TestHits::hPullGP_Y_ts_mono
std::map< std::string, TH1F * > hPullGP_Y_ts_mono
Definition: TestHits.h:96
TestHits::theBuilder
edm::ESHandle< TransientTrackingRecHitBuilder > theBuilder
Definition: TestHits.h:73
GlobalErrorBase::czz
T czz() const
Definition: GlobalErrorBase.h:107
TestHits::hPullGP_X_ts_mono
std::map< std::string, TH1F * > hPullGP_X_ts_mono
Definition: TestHits.h:95
StripGeomDetUnit::specificTopology
virtual const StripTopology & specificTopology() const
Returns a reference to the strip proxy topology.
Definition: StripGeomDetUnit.cc:17
TestHits::hPullGP_X_rs
std::map< std::string, TH1F * > hPullGP_X_rs
Definition: TestHits.h:85
TestHits::hPullGP_Y_rs_mono
std::map< std::string, TH1F * > hPullGP_Y_rs_mono
Definition: TestHits.h:102
TestHits::hPullGP_X_ts
std::map< std::string, TH1F * > hPullGP_X_ts
Definition: TestHits.h:79
mps_fire.result
result
Definition: mps_fire.py:303
TestHits::hPullGM_Z_ts
std::map< std::string, TH1F * > hPullGM_Z_ts
Definition: TestHits.h:84
TestHits::hPullGP_X_tr_stereo
std::map< std::string, TH1F * > hPullGP_X_tr_stereo
Definition: TestHits.h:117
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:671
TrackerTopologyRcd
Definition: TrackerTopologyRcd.h:10
PTrajectoryStateOnDet
Definition: PTrajectoryStateOnDet.h:10
heavyFlavorValidationHarvestingSequence_cff.combiner
combiner
Definition: heavyFlavorValidationHarvestingSequence_cff.py:107
TestHits::hProcess_vs_Chi2
TH2F * hProcess_vs_Chi2
Definition: TestHits.h:93
PSimHit
Definition: PSimHit.h:15
TestHits::hPullGP_Z_tr
std::map< std::string, TH1F * > hPullGP_Z_tr
Definition: TestHits.h:90
dqmiolumiharvest.j
j
Definition: dqmiolumiharvest.py:66
TestHits::hPullGM_X_ts_stereo
std::map< std::string, TH1F * > hPullGM_X_ts_stereo
Definition: TestHits.h:111
StripTopology
Definition: StripTopology.h:11
TestHits::hPullGP_Y_ts_stereo
std::map< std::string, TH1F * > hPullGP_Y_ts_stereo
Definition: TestHits.h:109
TestHits::hPullGP_Z_rs_mono
std::map< std::string, TH1F * > hPullGP_Z_rs_mono
Definition: TestHits.h:103
CartesianTrajectoryError::position
const GlobalError position() const
Position error submatrix.
Definition: CartesianTrajectoryError.h:41
fit
Definition: CombinedChiSquaredLikelihood.h:6
TrajectoryStateOnSurface::isValid
bool isValid() const
Definition: TrajectoryStateOnSurface.h:54
hit
Definition: SiStripHitEffFromCalibTree.cc:88
DeadROC_duringRun.dir
dir
Definition: DeadROC_duringRun.py:23
StripGeomDetUnit
Definition: StripGeomDetUnit.h:15
TrackingComponentsRecord
Definition: TrackingComponentsRecord.h:12
Generated for CMSSW Reference Manual by doxygen 1.8.16