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;
514 if ( G4Trk[0] != 0 || G4Trk[1] != 0)
phi = atan2(G4Trk[1],G4Trk[0]);
516 if (eta>0.0 && eta<=0.5) ir = 0;
517 if (eta>0.5 && eta<=1.0) ir = 1;
518 if (eta>1.0 && eta<=1.5) ir = 2;
519 if (eta>1.5 && eta<=2.0) ir = 3;
520 if (eta>2.0 && eta<=2.5) ir = 4;
523 if (eta>-0.5 && eta<= 0.0) ir = 6;
524 if (eta>-1.0 && eta<=-0.5) ir = 7;
525 if (eta>-1.5 && eta<=-1.0) ir = 8;
526 if (eta>-2.0 && eta<=-1.5) ir = 9;
527 if (eta>-2.5 && eta<=-2.0) ir = 10;
528 if (eta<=-2.5) ir = 11;
538 for (itHit = PxlBrlLowContainer->begin(); itHit != PxlBrlLowContainer->end(); ++itHit) {
541 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
547 h5e[ir]->
Fill(itHit->energyLoss());
548 h5ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
549 h5ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
550 h5ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
551 h5lx[ir]->
Fill(itHit->localPosition().x());
552 h5ly[ir]->
Fill(itHit->localPosition().y());
554 for (itHit = PxlBrlHighContainer->begin(); itHit != PxlBrlHighContainer->end(); ++itHit) {
557 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
564 h5e[ir]->
Fill(itHit->energyLoss());
565 h5ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
566 h5ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
567 h5ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
568 h5lx[ir]->
Fill(itHit->localPosition().x());
569 h5ly[ir]->
Fill(itHit->localPosition().y());
571 for (itHit = PxlFwdLowContainer->begin(); itHit != PxlFwdLowContainer->end(); ++itHit) {
574 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
581 h6e[ir]->
Fill(itHit->energyLoss());
582 h6ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
583 h6ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
584 h6ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
585 h6lx[ir]->
Fill(itHit->localPosition().x());
586 h6ly[ir]->
Fill(itHit->localPosition().y());
588 for (itHit = PxlFwdHighContainer->begin(); itHit != PxlFwdHighContainer->end(); ++itHit) {
591 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
598 h6e[ir]->
Fill(itHit->energyLoss());
599 h6ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
600 h6ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
601 h6ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
602 h6lx[ir]->
Fill(itHit->localPosition().x());
603 h6ly[ir]->
Fill(itHit->localPosition().y());
608 for (itHit = SiTIBLowContainer->begin(); itHit != SiTIBLowContainer->end(); ++itHit) {
611 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
618 h1e[ir]->
Fill(itHit->energyLoss());
619 h1ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
620 h1ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
621 h1ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
622 h1lx[ir]->
Fill(itHit->localPosition().x());
623 h1ly[ir]->
Fill(itHit->localPosition().y());
625 for (itHit = SiTIBHighContainer->begin(); itHit != SiTIBHighContainer->end(); ++itHit) {
628 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
635 h1e[ir]->
Fill(itHit->energyLoss());
636 h1ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
637 h1ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
638 h1ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
639 h1lx[ir]->
Fill(itHit->localPosition().x());
640 h1ly[ir]->
Fill(itHit->localPosition().y());
645 for (itHit = SiTOBLowContainer->begin(); itHit != SiTOBLowContainer->end(); ++itHit) {
648 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
655 h2e[ir]->
Fill(itHit->energyLoss());
656 h2ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
657 h2ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
658 h2ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
659 h2lx[ir]->
Fill(itHit->localPosition().x());
660 h2ly[ir]->
Fill(itHit->localPosition().y());
662 for (itHit = SiTOBHighContainer->begin(); itHit != SiTOBHighContainer->end(); ++itHit) {
665 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
671 h2e[ir]->
Fill(itHit->energyLoss());
672 h2ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
673 h2ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
674 h2ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
675 h2lx[ir]->
Fill(itHit->localPosition().x());
676 h2ly[ir]->
Fill(itHit->localPosition().y());
681 for (itHit = SiTIDLowContainer->begin(); itHit != SiTIDLowContainer->end(); ++itHit) {
684 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
690 h3e[ir]->
Fill(itHit->energyLoss());
691 h3ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
692 h3ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
693 h3ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
694 h3lx[ir]->
Fill(itHit->localPosition().x());
695 h3ly[ir]->
Fill(itHit->localPosition().y());
697 for (itHit = SiTIDHighContainer->begin(); itHit != SiTIDHighContainer->end(); ++itHit) {
700 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
706 h3e[ir]->
Fill(itHit->energyLoss());
707 h3ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
708 h3ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
709 h3ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
710 h3lx[ir]->
Fill(itHit->localPosition().x());
711 h3ly[ir]->
Fill(itHit->localPosition().y());
716 for (itHit = SiTECLowContainer->begin(); itHit != SiTECLowContainer->end(); ++itHit) {
719 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
725 h4e[ir]->
Fill(itHit->energyLoss());
726 h4ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
727 h4ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
728 h4ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
729 h4lx[ir]->
Fill(itHit->localPosition().x());
730 h4ly[ir]->
Fill(itHit->localPosition().y());
732 for (itHit = SiTECHighContainer->begin(); itHit != SiTECHighContainer->end(); ++itHit) {
735 GlobalPoint gpos=det->toGlobal(itHit->localPosition());
741 h4e[ir]->
Fill(itHit->energyLoss());
742 h4ex[ir]->
Fill(itHit->entryPoint().x()-itHit->exitPoint().x());
743 h4ey[ir]->
Fill(itHit->entryPoint().y()-itHit->exitPoint().y());
744 h4ez[ir]->
Fill(std::fabs(itHit->entryPoint().z()-itHit->exitPoint().z()));
745 h4lx[ir]->
Fill(itHit->localPosition().x());
746 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]
Log< T >::type log(const T &t)
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_