Analyze.
335 edm::LogInfo(
"EventInfo") <<
" Run = " <<
e.id().run() <<
" Event = " <<
e.id().event();
338 edm::PSimHitContainer::const_iterator itHit;
345 if (!PxlBrlLowContainer.
isValid()) {
346 edm::LogError(
"TrackerHitAnalyzer::analyze") <<
"Unable to find TrackerHitsPixelBarrelLowTof in event!";
352 if (!PxlBrlHighContainer.
isValid()) {
353 edm::LogError(
"TrackerHitAnalyzer::analyze") <<
"Unable to find TrackerHitsPixelBarrelHighTof in event!";
362 if (!PxlFwdLowContainer.
isValid()) {
363 edm::LogError(
"TrackerHitAnalyzer::analyze") <<
"Unable to find TrackerHitsPixelEndcapLowTof in event!";
369 if (!PxlFwdHighContainer.
isValid()) {
370 edm::LogError(
"TrackerHitAnalyzer::analyze") <<
"Unable to find TrackerHitsPixelEndcapHighTof in event!";
379 if (!SiTIBLowContainer.
isValid()) {
380 edm::LogError(
"TrackerHitProducer::analyze") <<
"Unable to find TrackerHitsTIBLowTof in event!";
387 if (!SiTIBHighContainer.
isValid()) {
388 edm::LogError(
"TrackerHitProducer::analyze") <<
"Unable to find TrackerHitsTIBHighTof in event!";
397 if (!SiTOBLowContainer.
isValid()) {
398 edm::LogError(
"TrackerHitProducer::analyze") <<
"Unable to find TrackerHitsTOBLowTof in event!";
405 if (!SiTOBHighContainer.
isValid()) {
406 edm::LogError(
"TrackerHitProducer::analyze") <<
"Unable to find TrackerHitsTOBHighTof in event!";
416 if (!SiTIDLowContainer.
isValid()) {
417 edm::LogError(
"TrackerHitProducer::analyze") <<
"Unable to find TrackerHitsTIDLowTof in event!";
424 if (!SiTIDHighContainer.
isValid()) {
425 edm::LogError(
"TrackerHitProducer::analyze") <<
"Unable to find TrackerHitsTIDHighTof in event!";
434 if (!SiTECLowContainer.
isValid()) {
435 edm::LogError(
"TrackerHitProducer::analyze") <<
"Unable to find TrackerHitsTECLowTof in event!";
442 if (!SiTECHighContainer.
isValid()) {
443 edm::LogError(
"TrackerHitProducer::analyze") <<
"Unable to find TrackerHitsTECHighTof in event!";
453 if (!G4TrkContainer.
isValid()) {
454 edm::LogError(
"TrackerHitAnalyzer::analyze") <<
"Unable to find SimTrack in event!";
462 edm::SimTrackContainer::const_iterator itTrk;
463 for (itTrk = G4TrkContainer->begin(); itTrk != G4TrkContainer->end(); ++itTrk) {
465 double eta = 0,
p = 0;
466 const CLHEP::HepLorentzVector &G4Trk = CLHEP::HepLorentzVector(
467 itTrk->momentum().x(), itTrk->momentum().y(), itTrk->momentum().z(), itTrk->momentum().e());
468 p =
sqrt(G4Trk[0] * G4Trk[0] + G4Trk[1] * G4Trk[1] + G4Trk[2] * G4Trk[2]);
470 edm::LogError(
"TrackerHitAnalyzer::analyze") <<
"TrackerTest::INFO: Primary has p = 0 ";
472 double costheta = G4Trk[2] /
p;
476 if (
eta > 0.0 &&
eta <= 0.5)
478 if (
eta > 0.5 &&
eta <= 1.0)
480 if (
eta > 1.0 &&
eta <= 1.5)
482 if (
eta > 1.5 &&
eta <= 2.0)
484 if (
eta > 2.0 &&
eta <= 2.5)
489 if (
eta > -0.5 &&
eta <= 0.0)
491 if (
eta > -1.0 &&
eta <= -0.5)
493 if (
eta > -1.5 &&
eta <= -1.0)
495 if (
eta > -2.0 &&
eta <= -1.5)
497 if (
eta > -2.5 &&
eta <= -2.0)
514 for (itHit = PxlBrlLowContainer->begin(); itHit != PxlBrlLowContainer->end(); ++itHit) {
517 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
523 h5e[ir]->
Fill(itHit->energyLoss());
524 h5ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
525 h5ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
526 h5ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
527 h5lx[ir]->
Fill(itHit->localPosition().x());
528 h5ly[ir]->
Fill(itHit->localPosition().y());
530 for (itHit = PxlBrlHighContainer->begin(); itHit != PxlBrlHighContainer->end(); ++itHit) {
533 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
539 h5e[ir]->
Fill(itHit->energyLoss());
540 h5ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
541 h5ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
542 h5ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
543 h5lx[ir]->
Fill(itHit->localPosition().x());
544 h5ly[ir]->
Fill(itHit->localPosition().y());
546 for (itHit = PxlFwdLowContainer->begin(); itHit != PxlFwdLowContainer->end(); ++itHit) {
549 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
555 h6e[ir]->
Fill(itHit->energyLoss());
556 h6ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
557 h6ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
558 h6ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
559 h6lx[ir]->
Fill(itHit->localPosition().x());
560 h6ly[ir]->
Fill(itHit->localPosition().y());
562 for (itHit = PxlFwdHighContainer->begin(); itHit != PxlFwdHighContainer->end(); ++itHit) {
565 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
571 h6e[ir]->
Fill(itHit->energyLoss());
572 h6ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
573 h6ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
574 h6ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
575 h6lx[ir]->
Fill(itHit->localPosition().x());
576 h6ly[ir]->
Fill(itHit->localPosition().y());
582 for (itHit = SiTIBLowContainer->begin(); itHit != SiTIBLowContainer->end(); ++itHit) {
585 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
591 h1e[ir]->
Fill(itHit->energyLoss());
592 h1ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
593 h1ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
594 h1ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
595 h1lx[ir]->
Fill(itHit->localPosition().x());
596 h1ly[ir]->
Fill(itHit->localPosition().y());
598 for (itHit = SiTIBHighContainer->begin(); itHit != SiTIBHighContainer->end(); ++itHit) {
601 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
607 h1e[ir]->
Fill(itHit->energyLoss());
608 h1ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
609 h1ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
610 h1ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
611 h1lx[ir]->
Fill(itHit->localPosition().x());
612 h1ly[ir]->
Fill(itHit->localPosition().y());
617 for (itHit = SiTOBLowContainer->begin(); itHit != SiTOBLowContainer->end(); ++itHit) {
620 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
626 h2e[ir]->
Fill(itHit->energyLoss());
627 h2ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
628 h2ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
629 h2ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
630 h2lx[ir]->
Fill(itHit->localPosition().x());
631 h2ly[ir]->
Fill(itHit->localPosition().y());
633 for (itHit = SiTOBHighContainer->begin(); itHit != SiTOBHighContainer->end(); ++itHit) {
636 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
642 h2e[ir]->
Fill(itHit->energyLoss());
643 h2ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
644 h2ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
645 h2ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
646 h2lx[ir]->
Fill(itHit->localPosition().x());
647 h2ly[ir]->
Fill(itHit->localPosition().y());
652 for (itHit = SiTIDLowContainer->begin(); itHit != SiTIDLowContainer->end(); ++itHit) {
655 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
661 h3e[ir]->
Fill(itHit->energyLoss());
662 h3ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
663 h3ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
664 h3ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
665 h3lx[ir]->
Fill(itHit->localPosition().x());
666 h3ly[ir]->
Fill(itHit->localPosition().y());
668 for (itHit = SiTIDHighContainer->begin(); itHit != SiTIDHighContainer->end(); ++itHit) {
671 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
677 h3e[ir]->
Fill(itHit->energyLoss());
678 h3ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
679 h3ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
680 h3ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
681 h3lx[ir]->
Fill(itHit->localPosition().x());
682 h3ly[ir]->
Fill(itHit->localPosition().y());
687 for (itHit = SiTECLowContainer->begin(); itHit != SiTECLowContainer->end(); ++itHit) {
690 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
696 h4e[ir]->
Fill(itHit->energyLoss());
697 h4ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
698 h4ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
699 h4ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
700 h4lx[ir]->
Fill(itHit->localPosition().x());
701 h4ly[ir]->
Fill(itHit->localPosition().y());
703 for (itHit = SiTECHighContainer->begin(); itHit != SiTECHighContainer->end(); ++itHit) {
706 GlobalPoint gpos = det->toGlobal(itHit->localPosition());
712 h4e[ir]->
Fill(itHit->energyLoss());
713 h4ex[ir]->
Fill(itHit->entryPoint().x() - itHit->exitPoint().x());
714 h4ey[ir]->
Fill(itHit->entryPoint().y() - itHit->exitPoint().y());
715 h4ez[ir]->
Fill(std::fabs(itHit->entryPoint().z() - itHit->exitPoint().z()));
716 h4lx[ir]->
Fill(itHit->localPosition().x());
717 h4ly[ir]->
Fill(itHit->localPosition().y());