Analyze.
346 edm::PSimHitContainer::const_iterator itHit;
353 if (!PxlBrlLowContainer.
isValid()) {
355 <<
"Unable to find TrackerHitsPixelBarrelLowTof in event!";
361 if (!PxlBrlHighContainer.
isValid()) {
363 <<
"Unable to find TrackerHitsPixelBarrelHighTof in event!";
372 if (!PxlFwdLowContainer.
isValid()) {
374 <<
"Unable to find TrackerHitsPixelEndcapLowTof in event!";
380 if (!PxlFwdHighContainer.
isValid()) {
382 <<
"Unable to find TrackerHitsPixelEndcapHighTof in event!";
393 if (!SiTIBLowContainer.
isValid()) {
395 <<
"Unable to find TrackerHitsTIBLowTof in event!";
403 if (!SiTIBHighContainer.
isValid()) {
405 <<
"Unable to find TrackerHitsTIBHighTof in event!";
415 if (!SiTOBLowContainer.
isValid()) {
417 <<
"Unable to find TrackerHitsTOBLowTof in event!";
425 if (!SiTOBHighContainer.
isValid()) {
427 <<
"Unable to find TrackerHitsTOBHighTof in event!";
438 if (!SiTIDLowContainer.
isValid()) {
440 <<
"Unable to find TrackerHitsTIDLowTof in event!";
448 if (!SiTIDHighContainer.
isValid()) {
450 <<
"Unable to find TrackerHitsTIDHighTof in event!";
460 if (!SiTECLowContainer.
isValid()) {
462 <<
"Unable to find TrackerHitsTECLowTof in event!";
470 if (!SiTECHighContainer.
isValid()) {
472 <<
"Unable to find TrackerHitsTECHighTof in event!";
482 if (!G4TrkContainer.
isValid()) {
484 <<
"Unable to find SimTrack in event!";
496 edm::SimTrackContainer::const_iterator itTrk;
497 for (itTrk = G4TrkContainer->begin(); itTrk != G4TrkContainer->end();
502 const CLHEP::HepLorentzVector& G4Trk = CLHEP::HepLorentzVector(itTrk->momentum().x(),
503 itTrk->momentum().y(),
504 itTrk->momentum().z(),
505 itTrk->momentum().e());
506 p =
sqrt(G4Trk[0]*G4Trk[0]+G4Trk[1]*G4Trk[1]+G4Trk[2]*G4Trk[2]);
509 <<
"TrackerTest::INFO: Primary has p = 0 ";
511 double costheta = G4Trk[2]/
p;
515 if (eta>0.0 && eta<=0.5) ir = 0;
516 if (eta>0.5 && eta<=1.0) ir = 1;
517 if (eta>1.0 && eta<=1.5) ir = 2;
518 if (eta>1.5 && eta<=2.0) ir = 3;
519 if (eta>2.0 && eta<=2.5) ir = 4;
522 if (eta>-0.5 && eta<= 0.0) ir = 6;
523 if (eta>-1.0 && eta<=-0.5) ir = 7;
524 if (eta>-1.5 && eta<=-1.0) ir = 8;
525 if (eta>-2.0 && eta<=-1.5) ir = 9;
526 if (eta>-2.5 && eta<=-2.0) ir = 10;
527 if (eta<=-2.5) ir = 11;
537 for (itHit = PxlBrlLowContainer->begin(); itHit != PxlBrlLowContainer->end(); ++itHit) {
540 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
546 h5e[ir]->
Fill(itHit->energyLoss());
547 h5ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
548 h5ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
549 h5ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
550 h5lx[ir]->
Fill(itHit->localPosition().x());
551 h5ly[ir]->
Fill(itHit->localPosition().y());
553 for (itHit = PxlBrlHighContainer->begin(); itHit != PxlBrlHighContainer->end(); ++itHit) {
556 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
563 h5e[ir]->
Fill(itHit->energyLoss());
564 h5ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
565 h5ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
566 h5ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
567 h5lx[ir]->
Fill(itHit->localPosition().x());
568 h5ly[ir]->
Fill(itHit->localPosition().y());
570 for (itHit = PxlFwdLowContainer->begin(); itHit != PxlFwdLowContainer->end(); ++itHit) {
573 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
580 h6e[ir]->
Fill(itHit->energyLoss());
581 h6ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
582 h6ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
583 h6ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
584 h6lx[ir]->
Fill(itHit->localPosition().x());
585 h6ly[ir]->
Fill(itHit->localPosition().y());
587 for (itHit = PxlFwdHighContainer->begin(); itHit != PxlFwdHighContainer->end(); ++itHit) {
590 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
597 h6e[ir]->
Fill(itHit->energyLoss());
598 h6ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
599 h6ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
600 h6ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
601 h6lx[ir]->
Fill(itHit->localPosition().x());
602 h6ly[ir]->
Fill(itHit->localPosition().y());
607 for (itHit = SiTIBLowContainer->begin(); itHit != SiTIBLowContainer->end(); ++itHit) {
610 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
617 h1e[ir]->
Fill(itHit->energyLoss());
618 h1ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
619 h1ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
620 h1ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
621 h1lx[ir]->
Fill(itHit->localPosition().x());
622 h1ly[ir]->
Fill(itHit->localPosition().y());
624 for (itHit = SiTIBHighContainer->begin(); itHit != SiTIBHighContainer->end(); ++itHit) {
627 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
634 h1e[ir]->
Fill(itHit->energyLoss());
635 h1ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
636 h1ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
637 h1ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
638 h1lx[ir]->
Fill(itHit->localPosition().x());
639 h1ly[ir]->
Fill(itHit->localPosition().y());
644 for (itHit = SiTOBLowContainer->begin(); itHit != SiTOBLowContainer->end(); ++itHit) {
647 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
654 h2e[ir]->
Fill(itHit->energyLoss());
655 h2ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
656 h2ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
657 h2ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
658 h2lx[ir]->
Fill(itHit->localPosition().x());
659 h2ly[ir]->
Fill(itHit->localPosition().y());
661 for (itHit = SiTOBHighContainer->begin(); itHit != SiTOBHighContainer->end(); ++itHit) {
664 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
670 h2e[ir]->
Fill(itHit->energyLoss());
671 h2ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
672 h2ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
673 h2ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
674 h2lx[ir]->
Fill(itHit->localPosition().x());
675 h2ly[ir]->
Fill(itHit->localPosition().y());
680 for (itHit = SiTIDLowContainer->begin(); itHit != SiTIDLowContainer->end(); ++itHit) {
683 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
689 h3e[ir]->
Fill(itHit->energyLoss());
690 h3ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
691 h3ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
692 h3ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
693 h3lx[ir]->
Fill(itHit->localPosition().x());
694 h3ly[ir]->
Fill(itHit->localPosition().y());
696 for (itHit = SiTIDHighContainer->begin(); itHit != SiTIDHighContainer->end(); ++itHit) {
699 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
705 h3e[ir]->
Fill(itHit->energyLoss());
706 h3ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
707 h3ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
708 h3ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
709 h3lx[ir]->
Fill(itHit->localPosition().x());
710 h3ly[ir]->
Fill(itHit->localPosition().y());
715 for (itHit = SiTECLowContainer->begin(); itHit != SiTECLowContainer->end(); ++itHit) {
718 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
724 h4e[ir]->
Fill(itHit->energyLoss());
725 h4ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
726 h4ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
727 h4ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
728 h4lx[ir]->
Fill(itHit->localPosition().x());
729 h4ly[ir]->
Fill(itHit->localPosition().y());
731 for (itHit = SiTECHighContainer->begin(); itHit != SiTECHighContainer->end(); ++itHit) {
734 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
740 h4e[ir]->
Fill(itHit->energyLoss());
741 h4ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
742 h4ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
743 h4ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
744 h4lx[ir]->
Fill(itHit->localPosition().x());
745 h4ly[ir]->
Fill(itHit->localPosition().y());
EventNumber_t event() const
MonitorElement * h2lx[12]
edm::InputTag PxlBrlLowSrc_
MonitorElement * h1lx[12]
MonitorElement * h6ex[12]
MonitorElement * h1ey[12]
MonitorElement * h3ly[12]
MonitorElement * h5ez[12]
MonitorElement * h4ex[12]
edm::InputTag SiTIBLowSrc_
MonitorElement * h1ez[12]
MonitorElement * h3ez[12]
edm::InputTag SiTIDLowSrc_
Geom::Phi< T > phi() const
Geom::Theta< T > theta() const
edm::InputTag SiTIBHighSrc_
MonitorElement * h5ey[12]
MonitorElement * h3lx[12]
bool getByType(Handle< PROD > &result) const
edm::InputTag PxlFwdLowSrc_
MonitorElement * h6ly[12]
MonitorElement * h5ex[12]
edm::InputTag SiTECLowSrc_
MonitorElement * h5ly[12]
edm::InputTag SiTECHighSrc_
MonitorElement * h1ex[12]
edm::InputTag SiTOBHighSrc_
MonitorElement * h6ey[12]
Tan< T >::type tan(const T &t)
MonitorElement * h4lx[12]
edm::InputTag SiTOBLowSrc_
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
MonitorElement * h2ex[12]
MonitorElement * h6ez[12]
MonitorElement * h2ez[12]
MonitorElement * h5lx[12]
MonitorElement * h4ey[12]
MonitorElement * h1ly[12]
MonitorElement * h4ly[12]
MonitorElement * h3ex[12]
MonitorElement * h3ey[12]
MonitorElement * h2ey[12]
edm::InputTag PxlBrlHighSrc_
MonitorElement * h2ly[12]
edm::InputTag SiTIDHighSrc_
MonitorElement * h6lx[12]
MonitorElement * h4ez[12]
edm::InputTag PxlFwdHighSrc_