CMS 3D CMS Logo

List of all members | Public Member Functions | Static Public Member Functions | Private Types | Private Attributes
PackedCandidateTrackValidator Class Reference
Inheritance diagram for PackedCandidateTrackValidator:
one::DQMEDAnalyzer< T > one::dqmimplementation::DQMBaseClass< T... >

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
 PackedCandidateTrackValidator (const edm::ParameterSet &pset)
 
 ~PackedCandidateTrackValidator () override
 
- Public Member Functions inherited from one::DQMEDAnalyzer< T >
 DQMEDAnalyzer ()=default
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
 
 DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
 
 ~DQMEDAnalyzer () override=default
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 

Private Types

enum  {
  sf_AllTracks = 0, sf_AssociatedToPC = 1, sf_PCIsCharged = 2, sf_PCHasTrack = 3,
  sf_PCIsNotElectron = 4, sf_PCHasHits = 5, sf_PCNdofNot0 = 6, sf_NoMissingInnerHits = 7
}
 

Private Attributes

bool debug_
 
MonitorElementh_diffCharge
 
PackedValueCheck< Float16Helper > h_diffCovDszDsz
 
PackedValueCheck< Float16Helper > h_diffCovDxyDsz
 
PackedValueCheck< Float16Helper > h_diffCovDxyDxy
 
PackedValueCheck< LogIntHelper > h_diffCovLambdaDsz
 
PackedValueCheck< LogIntHelper > h_diffCovLambdaLambda
 
PackedValueCheck< LogIntHelper > h_diffCovPhiDxy
 
PackedValueCheck< LogIntHelper > h_diffCovPhiPhi
 
PackedValueCheck< LogIntHelper > h_diffCovQoverpQoverp
 
MonitorElementh_diffDszError
 
PackedValueCheck< Float16Helper > h_diffDxyAssocPV
 
MonitorElementh_diffDxyError
 
MonitorElementh_diffDxyPV
 
PackedValueCheck< Float16Helper > h_diffDzAssocPV
 
MonitorElementh_diffDzError
 
MonitorElementh_diffDzPV
 
MonitorElementh_diffEta
 
MonitorElementh_diffEtaError
 
MonitorElementh_diffHitPatternHasValidHitInFirstPixelBarrel
 
MonitorElementh_diffHitPatternNumberOfLostInnerHits
 
MonitorElementh_diffHitPatternNumberOfValidHits
 
MonitorElementh_diffHitPatternNumberOfValidPixelHits
 
MonitorElementh_diffHitPatternPixelLayersWithMeasurement
 
MonitorElementh_diffHitPatternStripLayersWithMeasurement
 
MonitorElementh_diffHitPatternTrackerLayersWithMeasurement
 
MonitorElementh_diffIsHighPurity
 
MonitorElementh_diffLostInnerHits
 
MonitorElementh_diffNdof
 
MonitorElementh_diffNormalizedChi2
 
MonitorElementh_diffNumberOfHits
 
MonitorElementh_diffNumberOfPixelHits
 
MonitorElementh_diffNumberOfPixelLayers
 
MonitorElementh_diffNumberOfStripLayers
 
MonitorElementh_diffPhi
 
MonitorElementh_diffPt
 
MonitorElementh_diffPtError
 
MonitorElementh_diffTrackDxyAssocPV
 
MonitorElementh_diffTrackDxyError
 
MonitorElementh_diffTrackDzAssocPV
 
MonitorElementh_diffTrackDzError
 
MonitorElementh_diffVx
 
MonitorElementh_diffVy
 
MonitorElementh_diffVz
 
MonitorElementh_numberHitsOverMax
 
MonitorElementh_numberLayersOverMax
 
MonitorElementh_numberPixelHitsOverMax
 
MonitorElementh_numberPixelLayersOverMax
 
MonitorElementh_numberStripHitsOverMax
 
MonitorElementh_numberStripLayersOverMax
 
MonitorElementh_selectionFlow
 
std::string rootFolder_
 
edm::EDGetTokenT< edm::View< reco::Track > > tracksToken_
 
edm::EDGetTokenT< edm::Association< pat::PackedCandidateCollection > > trackToPackedCandidateToken_
 
edm::EDGetTokenT< reco::VertexCollectionverticesToken_
 

Detailed Description

Definition at line 434 of file PackedCandidateTrackValidator.cc.

Member Enumeration Documentation

anonymous enum
private

Constructor & Destructor Documentation

PackedCandidateTrackValidator::PackedCandidateTrackValidator ( const edm::ParameterSet pset)

Definition at line 528 of file PackedCandidateTrackValidator.cc.

528  :
529  tracksToken_(consumes<edm::View<reco::Track>>(iConfig.getUntrackedParameter<edm::InputTag>("tracks"))),
530  verticesToken_(consumes<reco::VertexCollection>(iConfig.getUntrackedParameter<edm::InputTag>("vertices"))),
531  trackToPackedCandidateToken_(consumes<edm::Association<pat::PackedCandidateCollection>>(iConfig.getUntrackedParameter<edm::InputTag>("trackToPackedCandidateAssociation"))),
532  rootFolder_(iConfig.getUntrackedParameter<std::string>("rootFolder")),
533  debug_(iConfig.getUntrackedParameter<bool>("debug")),
534  h_diffDxyAssocPV(RangeAbs(0.001)),
535  h_diffDzAssocPV(RangeAbs(0.001)),
536  h_diffCovQoverpQoverp(Range(-1e-6, 0.13), -15, 0), // despite of ceil in pack, there is rounding in double->float
537  h_diffCovLambdaLambda(Range(-1e-6, 0.13), -20, -5), // despite of ceil in pack, there is rounding in double->float
538  h_diffCovLambdaDsz(RangeAbs(0.13), -17, -4),
539  h_diffCovPhiPhi(RangeAbs(0.13), -15, 0),
540  h_diffCovPhiDxy(RangeAbs(0.13), -17, -4),
541  h_diffCovDxyDxy(RangeAbs(0.001)),
542  h_diffCovDxyDsz(RangeAbs(0.001)),
543  h_diffCovDszDsz(RangeAbs(0.001))
544 {}
PackedValueCheck< LogIntHelper > h_diffCovLambdaDsz
PackedValueCheck< Float16Helper > h_diffDzAssocPV
PackedValueCheck< Float16Helper > h_diffDxyAssocPV
edm::EDGetTokenT< edm::View< reco::Track > > tracksToken_
PixelRecoRange< float > Range
edm::EDGetTokenT< edm::Association< pat::PackedCandidateCollection > > trackToPackedCandidateToken_
PackedValueCheck< LogIntHelper > h_diffCovLambdaLambda
PackedValueCheck< LogIntHelper > h_diffCovPhiDxy
PackedValueCheck< Float16Helper > h_diffCovDxyDxy
PackedValueCheck< LogIntHelper > h_diffCovPhiPhi
edm::EDGetTokenT< reco::VertexCollection > verticesToken_
PackedValueCheck< LogIntHelper > h_diffCovQoverpQoverp
PackedValueCheck< Float16Helper > h_diffCovDxyDsz
PackedValueCheck< Float16Helper > h_diffCovDszDsz
PackedCandidateTrackValidator::~PackedCandidateTrackValidator ( )
override

Definition at line 546 of file PackedCandidateTrackValidator.cc.

546 {}

Member Function Documentation

void PackedCandidateTrackValidator::analyze ( const edm::Event iEvent,
const edm::EventSetup iSetup 
)
override

Definition at line 658 of file PackedCandidateTrackValidator.cc.

References funct::abs(), reco::TrackBase::charge(), reco::TrackBase::chi2(), reco::TrackBase::covariance(), debug_, DEFINE_FWK_MODULE, reco::deltaPhi(), reco::TrackBase::dszError(), reco::TrackBase::dxy(), reco::TrackBase::dxyError(), reco::TrackBase::dz(), reco::TrackBase::dzError(), MillePedeFileConverter_cfg::e, reco::TrackBase::eta(), reco::TrackBase::etaError(), MonitorElement::Fill(), edm::Event::getByToken(), h_diffCharge, h_diffCovDszDsz, h_diffCovDxyDsz, h_diffCovDxyDxy, h_diffCovLambdaDsz, h_diffCovLambdaLambda, h_diffCovPhiDxy, h_diffCovPhiPhi, h_diffCovQoverpQoverp, h_diffDszError, h_diffDxyAssocPV, h_diffDxyError, h_diffDxyPV, h_diffDzAssocPV, h_diffDzError, h_diffDzPV, h_diffEta, h_diffEtaError, h_diffHitPatternHasValidHitInFirstPixelBarrel, h_diffHitPatternNumberOfLostInnerHits, h_diffHitPatternNumberOfValidHits, h_diffHitPatternNumberOfValidPixelHits, h_diffHitPatternPixelLayersWithMeasurement, h_diffHitPatternStripLayersWithMeasurement, h_diffHitPatternTrackerLayersWithMeasurement, h_diffIsHighPurity, h_diffLostInnerHits, h_diffNdof, h_diffNormalizedChi2, h_diffNumberOfHits, h_diffNumberOfPixelHits, h_diffNumberOfPixelLayers, h_diffNumberOfStripLayers, h_diffPhi, h_diffPt, h_diffPtError, h_diffTrackDxyAssocPV, h_diffTrackDxyError, h_diffTrackDzAssocPV, h_diffTrackDzError, h_diffVx, h_diffVy, h_diffVz, h_numberHitsOverMax, h_numberLayersOverMax, h_numberPixelHitsOverMax, h_numberPixelLayersOverMax, h_numberStripHitsOverMax, h_numberStripLayersOverMax, h_selectionFlow, reco::HitPattern::hasValidHitInPixelLayer(), reco::TrackBase::highPurity, reco::TrackBase::hitPattern(), mps_fire::i, reco::TrackBase::i_dsz, reco::TrackBase::i_dxy, reco::TrackBase::i_lambda, reco::TrackBase::i_phi, reco::TrackBase::i_qoverp, edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::isNull(), edm::Ref< C, T, F >::key(), reco::HitPattern::MISSING_INNER_HITS, reco::HitPattern::MISSING_OUTER_HITS, reco::TrackBase::momentum(), pat::PackedCandidate::moreLostInnerHits, reco::TrackBase::ndof(), pat::PackedCandidate::noLostInnerHits, reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfLostHits(), reco::HitPattern::numberOfValidHits(), reco::HitPattern::numberOfValidPixelHits(), reco::HitPattern::numberOfValidStripHits(), pat::PackedCandidate::oneLostInnerHit, reco::TrackBase::phi(), reco::TrackBase::phiError(), GeomDetEnumerators::PixelBarrel, reco::HitPattern::pixelLayersWithMeasurement(), reco::Vertex::position(), reco::TrackBase::pt(), reco::TrackBase::ptError(), MetAnalyzer::pv(), reco::TrackBase::quality(), reco::TrackBase::referencePoint(), sf_AllTracks, sf_AssociatedToPC, sf_NoMissingInnerHits, sf_PCHasHits, sf_PCHasTrack, sf_PCIsCharged, sf_PCIsNotElectron, sf_PCNdofNot0, reco::HitPattern::stripLayersWithMeasurement(), HiIsolationCommonParameters_cff::track, reco::HitPattern::trackerLayersWithMeasurement(), pat::PackedCandidate::trackPixelHitsMask, l1t::tracks, tracksToken_, pat::PackedCandidate::trackStripHitsMask, trackToPackedCandidateToken_, heppy_batch::val, pat::PackedCandidate::validHitInFirstPixelBarrelLayer, electrons_cff::vertices, verticesToken_, reco::TrackBase::vx(), reco::TrackBase::vy(), and reco::TrackBase::vz().

658  {
660  iEvent.getByToken(tracksToken_, htracks);
661  const auto& tracks = *htracks;
662 
664  iEvent.getByToken(verticesToken_, hvertices);
665  const auto& vertices = *hvertices;
666 
667  if(vertices.empty())
668  return;
669  const reco::Vertex& pv = vertices[0];
670 
673  const auto& trackToPackedCandidate = *hassoc;
674 
675  for(size_t i=0; i<tracks.size(); ++i) {
676  auto trackPtr = tracks.ptrAt(i);
677  const reco::Track& track = *trackPtr;
679 
680  pat::PackedCandidateRef pcRef = trackToPackedCandidate[trackPtr];
681  if(pcRef.isNull()) {
682  continue;
683  }
685 
686  // Filter out neutral PackedCandidates, some of them may have track associated, and for those the charge comparison fails
687  if(pcRef->charge() == 0) {
688  continue;
689  }
691 
692  const reco::Track *trackPcPtr = pcRef->bestTrack();
693  if(!trackPcPtr) {
694  continue;
695  }
697 
698  // Filter out electrons to avoid comparisons to PackedCandidates with GsfTrack
699  if(std::abs(pcRef->pdgId()) == 11) {
700  continue;
701  }
703 
704  // Filter out PackedCandidate-tracks with no hits, as they won't have their details filled
705  const reco::Track& trackPc = *trackPcPtr;
706  if(trackPc.hitPattern().numberOfValidHits() == 0) {
707  continue;
708  }
710 
711  auto slimmedVertexRef = pcRef->vertexRef();
712  const reco::Vertex& pcVertex = vertices[slimmedVertexRef.key()];
713 
714  fillNoFlow(h_diffVx, trackPc.vx() - track.vx());
715  fillNoFlow(h_diffVy, trackPc.vy() - track.vy());
716  fillNoFlow(h_diffVz, trackPc.vz() - track.vz());
717 
718  // PackedCandidate recalculates the ndof in unpacking as
719  // (nhits+npixelhits-5), but some strip hits may have dimension 2.
720  // If PackedCandidate has ndof=0, the resulting normalizedChi2
721  // will be 0 too. Hence, the comparison makes sense only for those
722  // PackedCandidates that have ndof != 0.
723  double diffNormalizedChi2 = 0;
724  if(trackPc.ndof() != 0) {
726  diffNormalizedChi2 = trackPc.normalizedChi2() - track.normalizedChi2();
727  fillNoFlow(h_diffNormalizedChi2, diffNormalizedChi2);
728  }
729  fillNoFlow(h_diffNdof, trackPc.ndof() - track.ndof());
730 
731  auto diffCharge = trackPc.charge() - track.charge();
732  fillNoFlow(h_diffCharge, diffCharge);
733  int diffHP = static_cast<int>(trackPc.quality(reco::TrackBase::highPurity)) - static_cast<int>(track.quality(reco::TrackBase::highPurity));
734  fillNoFlow(h_diffIsHighPurity, diffHP);
735 
736  const auto diffPt = diffRelative(trackPc.pt(), track.pt());
737  const auto diffPhi = reco::deltaPhi(trackPc.phi(), track.phi());
738  fillNoFlow(h_diffPt , diffPt);
739  fillNoFlow(h_diffEta, trackPc.eta() - track.eta());
740  fillNoFlow(h_diffPhi, diffPhi);
741 
742  const auto diffDxyAssocPV = h_diffDxyAssocPV.fill(pcRef->dxy( ) , track.dxy(pcVertex.position()), [](double value){return value*100.;});
743  const auto diffDzAssocPV = h_diffDzAssocPV.fill( pcRef->dzAssociatedPV(), track.dz(pcVertex.position()) , [](double value){return value*100.;});
744  const auto diffDxyPV = diffRelative(pcRef->dxy(pv.position()), track.dxy(pv.position()));
745  const auto diffDzPV = diffRelative(pcRef->dz(pv.position()) , track.dz(pv.position()));
746  fillNoFlow(h_diffDxyPV , diffDxyPV);
747  fillNoFlow(h_diffDzPV , diffDzPV);
748  fillNoFlow(h_diffTrackDxyAssocPV, diffRelative(trackPc.dxy(pcVertex.position()), track.dxy(pcVertex.position())));
749  fillNoFlow(h_diffTrackDzAssocPV , diffRelative(trackPc.dz(pcVertex.position()) , track.dz(pcVertex.position())));
750 
751  auto fillCov1 = [&](auto& hlp, const int i, const int j) {
752  return hlp.fill(trackPc.covariance(i, j), track.covariance(i, j));
753  };
754  auto fillCov2 = [&](auto& hlp, const int i, const int j, std::function<double(double)> modifyPack) {
755  return hlp.fill(trackPc.covariance(i, j), track.covariance(i, j), modifyPack);
756  };
757  auto fillCov3 = [&](auto& hlp, const int i, const int j, std::function<double(double)> modifyPack, std::function<double(double)> modifyUnpack) {
758  return hlp.fill(trackPc.covariance(i, j), track.covariance(i, j), modifyPack, modifyUnpack);
759  };
760 
761  const auto pcPt = pcRef->pt();
762  const auto diffCovQoverpQoverp = fillCov3(h_diffCovQoverpQoverp, reco::TrackBase::i_qoverp, reco::TrackBase::i_qoverp, [=](double val){return val*pcPt*pcPt;}, [=](double val){return val/pcPt/pcPt;});
763  const auto diffCovLambdaLambda = fillCov1(h_diffCovLambdaLambda, reco::TrackBase::i_lambda, reco::TrackBase::i_lambda);
764  const auto diffCovLambdaDsz = fillCov1(h_diffCovLambdaDsz, reco::TrackBase::i_lambda, reco::TrackBase::i_dsz);
765  const auto diffCovPhiPhi = fillCov3(h_diffCovPhiPhi, reco::TrackBase::i_phi, reco::TrackBase::i_phi, [=](double val){return val*pcPt*pcPt;}, [=](double val){return val/pcPt/pcPt;});
766  const auto diffCovPhiDxy = fillCov1(h_diffCovPhiDxy, reco::TrackBase::i_phi, reco::TrackBase::i_dxy);
767  const auto diffCovDxyDxy = fillCov2(h_diffCovDxyDxy, reco::TrackBase::i_dxy, reco::TrackBase::i_dxy, [](double value){return value*10000.;});
768  const auto diffCovDxyDsz = fillCov2(h_diffCovDxyDsz, reco::TrackBase::i_dxy, reco::TrackBase::i_dsz, [](double value){return value*10000.;});
769  const auto diffCovDszDsz = fillCov2(h_diffCovDszDsz, reco::TrackBase::i_dsz, reco::TrackBase::i_dsz, [](double value){return value*10000.;});
770 
771  if(isInRange(diffCovDszDsz.status())) {
772  fillNoFlow(h_diffDszError, diffRelative(pcRef->dzError(), track.dszError()));
773  fillNoFlow(h_diffDzError, diffRelative(pcRef->dzError(), track.dzError()));
774  fillNoFlow(h_diffTrackDzError, diffRelative(trackPc.dzError(), track.dzError()));
775  }
776  if(isInRange(diffCovDxyDxy.status())) {
777  fillNoFlow(h_diffDxyError, diffRelative(pcRef->dxyError(), track.dxyError()));
778  fillNoFlow(h_diffTrackDxyError, diffRelative(trackPc.dxyError(), track.dxyError()));
779  }
780  fillNoFlow(h_diffPtError , diffRelative(trackPc.ptError() , track.ptError() ));
781  fillNoFlow(h_diffEtaError , diffRelative(trackPc.etaError() , track.etaError() ));
782 
783 
784 
785  // For the non-HitPattern ones, take into account the PackedCandidate packing precision
786  const auto trackNumberOfHits = track.hitPattern().numberOfValidHits();
787  const auto trackNumberOfPixelHits = track.hitPattern().numberOfValidPixelHits();
788  const auto trackNumberOfStripHits = track.hitPattern().numberOfValidStripHits();
789  const auto pcNumberOfHits = pcRef->numberOfHits();
790  const auto pcNumberOfPixelHits = pcRef->numberOfPixelHits();
791  const auto pcNumberOfStripHits = pcNumberOfHits - pcNumberOfPixelHits;
792  const auto trackNumberOfLayers = track.hitPattern().trackerLayersWithMeasurement();
793  const auto trackNumberOfPixelLayers = track.hitPattern().pixelLayersWithMeasurement();
794  const auto trackNumberOfStripLayers = track.hitPattern().stripLayersWithMeasurement();
795  const auto pcNumberOfLayers = pcRef->trackerLayersWithMeasurement();
796  const auto pcNumberOfPixelLayers = pcRef->pixelLayersWithMeasurement();
797  const auto pcNumberOfStripLayers = pcRef->stripLayersWithMeasurement();
798 
799  // layer number overflow (should be zero)
800  const int pixelLayerOverflow = trackNumberOfPixelLayers > pat::PackedCandidate::trackPixelHitsMask ? trackNumberOfPixelLayers - pat::PackedCandidate::trackPixelHitsMask : 0;
801  const int stripLayerOverflow = trackNumberOfStripLayers > pat::PackedCandidate::trackStripHitsMask ? trackNumberOfStripLayers - pat::PackedCandidate::trackStripHitsMask : 0;
803 
804  // hit overflow (should also be zero)
805  const int pixelOverflow = trackNumberOfPixelHits - pcNumberOfPixelLayers > pat::PackedCandidate::trackPixelHitsMask ? trackNumberOfPixelHits - pcNumberOfPixelLayers - pat::PackedCandidate::trackPixelHitsMask : 0;
806  const int stripOverflow = trackNumberOfStripHits - pcNumberOfStripLayers > pat::PackedCandidate::trackStripHitsMask ? trackNumberOfStripHits - pcNumberOfStripLayers - pat::PackedCandidate::trackStripHitsMask : 0;
807  const int hitsOverflow = trackNumberOfHits - pcNumberOfLayers > (pat::PackedCandidate::trackPixelHitsMask+pat::PackedCandidate::trackStripHitsMask) ? trackNumberOfHits - pcNumberOfLayers - (pat::PackedCandidate::trackPixelHitsMask+pat::PackedCandidate::trackStripHitsMask) : 0;
808  // PackedCandidate counts overflow pixel hits as strip
809  const int pixelInducedStripOverflow = (trackNumberOfStripHits+pixelOverflow-pcNumberOfStripLayers) > pat::PackedCandidate::trackStripHitsMask ? (trackNumberOfStripHits+pixelOverflow-stripOverflow-pcNumberOfStripLayers) - pat::PackedCandidate::trackStripHitsMask : 0;
810  h_numberPixelLayersOverMax->Fill(pixelLayerOverflow);
811  h_numberStripLayersOverMax->Fill(stripLayerOverflow);
812  h_numberLayersOverMax->Fill(layerOverflow);
813  h_numberPixelHitsOverMax->Fill(pixelOverflow);
814  h_numberStripHitsOverMax->Fill(stripOverflow);
815  h_numberHitsOverMax->Fill(hitsOverflow);
816 
817  int diffNumberOfPixelHits = 0;
818  int diffNumberOfHits = 0;
819  int diffNumberOfPixelLayers = 0;
820  int diffNumberOfStripLayers = 0;
821  if(pixelLayerOverflow) {
822  diffNumberOfPixelLayers = pcNumberOfPixelLayers - pat::PackedCandidate::trackPixelHitsMask;
823  }
824  else {
825  diffNumberOfPixelLayers = pcNumberOfPixelLayers - trackNumberOfPixelLayers;
826  }
827  if(pixelOverflow) {
828  diffNumberOfPixelHits = pcNumberOfPixelHits - pcNumberOfPixelLayers - pat::PackedCandidate::trackPixelHitsMask;
829  }
830  else {
831  diffNumberOfPixelHits = pcNumberOfPixelHits - trackNumberOfPixelHits;
832  }
833  if(stripLayerOverflow) {
834  diffNumberOfStripLayers = pcNumberOfStripLayers - pat::PackedCandidate::trackStripHitsMask;
835  }
836  else {
837  diffNumberOfStripLayers = pcNumberOfStripLayers - trackNumberOfStripLayers;
838  }
839  if(stripOverflow || pixelInducedStripOverflow || pixelOverflow) {
840  int diffNumberOfStripHits = 0;
841  if(stripOverflow || pixelInducedStripOverflow) {
842  diffNumberOfStripHits = pcNumberOfStripHits - pat::PackedCandidate::trackStripHitsMask;
843  }
844  else if(pixelOverflow) {
845  diffNumberOfStripHits = (pcNumberOfStripHits - pixelOverflow) - trackNumberOfStripHits;
846  }
847 
848  diffNumberOfHits = diffNumberOfPixelHits + diffNumberOfStripHits;
849  }
850  else {
851  diffNumberOfHits = pcNumberOfHits - trackNumberOfHits;
852  }
853 
854 
855 
856  fillNoFlow(h_diffNumberOfPixelHits, diffNumberOfPixelHits);
857  fillNoFlow(h_diffNumberOfHits, diffNumberOfHits);
858  fillNoFlow(h_diffNumberOfPixelLayers, diffNumberOfPixelLayers);
859  fillNoFlow(h_diffNumberOfStripLayers, diffNumberOfStripLayers);
860 
861  int diffLostInnerHits = 0;
862  const auto trackLostInnerHits = track.hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS);
863  switch(pcRef->lostInnerHits()) {
866  diffLostInnerHits = -trackLostInnerHits;
867  break;
869  diffLostInnerHits = 1-trackLostInnerHits;
870  break;
872  diffLostInnerHits = trackLostInnerHits>=2 ? 0 : 2-trackLostInnerHits;
873  break;
874  }
875  fillNoFlow(h_diffLostInnerHits, diffLostInnerHits);
876 
877  // For HitPattern ones, calculate the full diff (i.e. some differences are expected)
878  auto diffHitPatternPixelLayersWithMeasurement = trackPc.hitPattern().pixelLayersWithMeasurement() - trackNumberOfPixelLayers;
879  fillNoFlow(h_diffHitPatternPixelLayersWithMeasurement, diffHitPatternPixelLayersWithMeasurement);
880  auto diffHitPatternStripLayersWithMeasurement = trackPc.hitPattern().stripLayersWithMeasurement() - trackNumberOfStripLayers;
881  fillNoFlow(h_diffHitPatternStripLayersWithMeasurement, diffHitPatternStripLayersWithMeasurement);
882  auto diffHitPatternTrackerLayersWithMeasurement = trackPc.hitPattern().trackerLayersWithMeasurement() - trackNumberOfLayers;
883  fillNoFlow(h_diffHitPatternTrackerLayersWithMeasurement, diffHitPatternTrackerLayersWithMeasurement);
884  auto diffHitPatternNumberOfValidPixelHits = trackPc.hitPattern().numberOfValidPixelHits() - trackNumberOfPixelHits;
885  fillNoFlow(h_diffHitPatternNumberOfValidPixelHits, diffHitPatternNumberOfValidPixelHits);
886  auto diffHitPatternNumberOfValidHits = trackPc.hitPattern().numberOfValidHits() - trackNumberOfHits;
887  fillNoFlow(h_diffHitPatternNumberOfValidHits, diffHitPatternNumberOfValidHits);
889 
890  // hasValidHitInFirstPixelBarrel is set only if numberOfLostHits(MISSING_INNER_HITS) == 0
891  int diffHitPatternHasValidHitInFirstPixelBarrel = 0;
894  diffHitPatternHasValidHitInFirstPixelBarrel = static_cast<int>(trackPc.hitPattern().hasValidHitInPixelLayer(PixelSubdetector::SubDetector::PixelBarrel, 1)) - static_cast<int>(track.hitPattern().hasValidHitInPixelLayer(PixelSubdetector::SubDetector::PixelBarrel, 1));
895  fillNoFlow(h_diffHitPatternHasValidHitInFirstPixelBarrel, diffHitPatternHasValidHitInFirstPixelBarrel);
896  }
897 
898  // Print warning if there are differences outside the expected range
899  if(debug_ &&
900  (diffNormalizedChi2 < -1 || diffNormalizedChi2 > 0 || diffCharge != 0 || diffHP != 0
901  || std::abs(diffPhi) > 5e-4
902  || diffDxyAssocPV.outsideExpectedRange()
903  || diffDzAssocPV.outsideExpectedRange()
904  || std::abs(diffDxyPV) > 0.05 || std::abs(diffDzPV) > 0.05
905  || diffCovQoverpQoverp.outsideExpectedRange() || diffCovLambdaLambda.outsideExpectedRange()
906  || diffCovLambdaDsz.outsideExpectedRange() || diffCovPhiPhi.outsideExpectedRange()
907  || diffCovPhiDxy.outsideExpectedRange() || diffCovDxyDxy.outsideExpectedRange()
908  || diffCovDxyDsz.outsideExpectedRange() || diffCovDszDsz.outsideExpectedRange()
909  || diffNumberOfPixelHits != 0 || diffNumberOfHits != 0 || diffLostInnerHits != 0
910  || diffHitPatternHasValidHitInFirstPixelBarrel != 0
911  )
912  ) {
913 
914  edm::LogInfo("PackedCandidateTrackValidator") << "Track " << i << " pt " << track.pt() << " eta " << track.eta() << " phi " << track.phi() << " chi2 " << track.chi2() << " ndof " << track.ndof()
915  << "\n"
916  << " ptError " << track.ptError() << " etaError " << track.etaError() << " phiError " << track.phiError() << " dxyError " << track.dxyError() << " dzError " << track.dzError()
917  << "\n"
918  << " refpoint " << track.referencePoint() << " momentum " << track.momentum()
919  << "\n"
920  << " dxy " << track.dxy() << " dz " << track.dz()
921  << "\n"
922  << " " << TrackAlgoPrinter(track)
923  << " lost inner hits " << trackLostInnerHits
924  << " lost outer hits " << track.hitPattern().numberOfLostHits(reco::HitPattern::MISSING_OUTER_HITS)
925  << " hitpattern " << HitPatternPrinter(track)
926  << " \n"
927  << " PC " << pcRef.id() << ":" << pcRef.key() << " track pt " << trackPc.pt() << " eta " << trackPc.eta() << " phi " << trackPc.phi() << " (PC " << pcRef->phi() << ") chi2 " << trackPc.chi2() << " ndof " << trackPc.ndof() << " pdgId " << pcRef->pdgId() << " mass " << pcRef->mass()
928  << "\n"
929  << " ptError " << trackPc.ptError() << " etaError " << trackPc.etaError() << " phiError " << trackPc.phiError()
930  << "\n"
931  << " pc.vertex " << pcRef->vertex() << " momentum " << pcRef->momentum() << " track " << trackPc.momentum()
932  << "\n"
933  << " dxy " << trackPc.dxy() << " dz " << trackPc.dz() << " pc.dz " << pcRef->dz()
934  << " dxyError " << trackPc.dxyError() << " dzError " << trackPc.dzError()
935  << "\n"
936  << " dxy(PV) " << trackPc.dxy(pv.position()) << " dz(PV) " << trackPc.dz(pv.position())
937  << " dxy(assocPV) " << trackPc.dxy(pcVertex.position()) << " dz(assocPV) " << trackPc.dz(pcVertex.position())
938  << "\n"
939  << " (diff PackedCandidate track)"
940  << " highPurity " << diffHP << " " << trackPc.quality(reco::TrackBase::highPurity) << " " << track.quality(reco::TrackBase::highPurity)
941  << " charge " << diffCharge << " " << trackPc.charge() << " " << track.charge()
942  << " normalizedChi2 " << diffNormalizedChi2 << " " << trackPc.normalizedChi2() << " " << track.normalizedChi2()
943  << "\n "
944  << " numberOfAllHits " << diffNumberOfHits << " " << pcNumberOfHits << " " << trackNumberOfHits
945  << " numberOfPixelHits " << diffNumberOfPixelHits << " " << pcNumberOfPixelHits << " " << trackNumberOfPixelHits
946  << " numberOfStripHits # " << pcNumberOfStripHits << " " << trackNumberOfStripHits
947  << "\n "
948  << " hitPattern.numberOfValidPixelHits " << diffHitPatternNumberOfValidPixelHits << " " << trackPc.hitPattern().numberOfValidPixelHits() << " " << track.hitPattern().numberOfValidPixelHits()
949  << " hitPattern.numberOfValidHits " << diffHitPatternNumberOfValidHits << " " << trackPc.hitPattern().numberOfValidHits() << " " << track.hitPattern().numberOfValidHits()
950  << " hitPattern.hasValidHitInFirstPixelBarrel " << diffHitPatternHasValidHitInFirstPixelBarrel << " " << trackPc.hitPattern().hasValidHitInPixelLayer(PixelSubdetector::SubDetector::PixelBarrel, 1) << " " << track.hitPattern().hasValidHitInPixelLayer(PixelSubdetector::SubDetector::PixelBarrel, 1)
951  << "\n "
952  << " lostInnerHits " << diffLostInnerHits << " " << pcRef->lostInnerHits() << " #"
953  << " phi (5e-4) " << diffPhi << " " << trackPc.phi() << " " << track.phi()
954  << "\n "
955  << " dxy(assocPV) " << diffDxyAssocPV
956  << "\n "
957  << " dz(assocPV) " << diffDzAssocPV
958  << "\n "
959  << " dxy(PV) (0.05) " << diffDxyPV << " " << pcRef->dxy(pv.position()) << " " << track.dxy(pv.position())
960  << "\n "
961  << " dz(PV) (0.05) " << diffDzPV << " " << pcRef->dz(pv.position()) << " " << track.dz(pv.position())
962  << "\n "
963  << " cov(qoverp, qoverp) " << diffCovQoverpQoverp
964  << "\n "
965  << " cov(lambda, lambda) " << diffCovLambdaLambda
966  << "\n "
967  << " cov(lambda, dsz) " << diffCovLambdaDsz
968  << "\n "
969  << " cov(phi, phi) " << diffCovPhiPhi
970  << "\n "
971  << " cov(phi, dxy) " << diffCovPhiDxy
972  << "\n "
973  << " cov(dxy, dxy) " << diffCovDxyDxy
974  << "\n "
975  << " cov(dxy, dsz) " << diffCovDxyDsz
976  << "\n "
977  << " cov(dsz, dsz) " << diffCovDszDsz;
978  }
979  }
980 }
constexpr double deltaPhi(double phi1, double phi2)
Definition: deltaPhi.h:22
const Point & referencePoint() const
Reference point on the track.
Definition: TrackBase.h:720
PackedValueCheck< LogIntHelper > h_diffCovLambdaDsz
PackedValueCheck< Float16Helper > h_diffDzAssocPV
PackedValueCheck< Float16Helper > h_diffDxyAssocPV
bool hasValidHitInPixelLayer(enum PixelSubdetector::SubDetector, uint16_t layer) const
Definition: HitPattern.cc:365
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
Definition: TrackBase.h:600
edm::EDGetTokenT< edm::View< reco::Track > > tracksToken_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
double dxyError() const
error on dxy
Definition: TrackBase.h:847
int numberOfValidHits() const
Definition: HitPattern.h:896
edm::EDGetTokenT< edm::Association< pat::PackedCandidateCollection > > trackToPackedCandidateToken_
double etaError() const
error on eta
Definition: TrackBase.h:835
int numberOfValidStripHits() const
Definition: HitPattern.h:931
double phi() const
azimuthal angle of momentum vector
Definition: TrackBase.h:684
key_type key() const
Accessor for product key.
Definition: Ref.h:263
int pixelLayersWithMeasurement() const
Definition: HitPattern.cc:538
const Vector & momentum() const
track momentum vector
Definition: TrackBase.h:714
int trackerLayersWithMeasurement() const
Definition: HitPattern.cc:557
const Point & position() const
position
Definition: Vertex.h:109
ProductID id() const
Accessor for product ID.
Definition: Ref.h:257
void Fill(long long x)
PackedValueCheck< LogIntHelper > h_diffCovLambdaLambda
double dszError() const
error on dsz
Definition: TrackBase.h:859
double eta() const
pseudorapidity of momentum vector
Definition: TrackBase.h:690
PackedValueCheck< LogIntHelper > h_diffCovPhiDxy
double chi2() const
chi-squared of the fit
Definition: TrackBase.h:588
PackedValueCheck< Float16Helper > h_diffCovDxyDxy
PackedValueCheck< LogIntHelper > h_diffCovPhiPhi
double ndof() const
number of degrees of freedom of the fit
Definition: TrackBase.h:594
CovarianceMatrix covariance() const
return track covariance matrix
Definition: TrackBase.h:782
double pt() const
track transverse momentum
Definition: TrackBase.h:660
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
Definition: TrackBase.h:814
double phiError() const
error on phi
Definition: TrackBase.h:841
def pv(vc)
Definition: MetAnalyzer.py:7
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
Definition: value.py:1
bool isNull() const
Checks for null.
Definition: Ref.h:248
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
Definition: TrackBase.h:648
double dzError() const
error on dz
Definition: TrackBase.h:865
double vz() const
z coordinate of the reference point on track
Definition: TrackBase.h:708
edm::EDGetTokenT< reco::VertexCollection > verticesToken_
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
Definition: TrackBase.h:485
int stripLayersWithMeasurement() const
Definition: HitPattern.h:1138
bool quality(const TrackQuality) const
Track quality.
Definition: TrackBase.h:549
int numberOfLostHits(HitCategory category) const
Definition: HitPattern.h:990
double vy() const
y coordinate of the reference point on track
Definition: TrackBase.h:702
PackedValueCheck< LogIntHelper > h_diffCovQoverpQoverp
PackedValueCheck< Float16Helper > h_diffCovDxyDsz
int numberOfValidPixelHits() const
Definition: HitPattern.h:916
PackedValueCheck< Float16Helper > h_diffCovDszDsz
int charge() const
track electric charge
Definition: TrackBase.h:606
double dxy() const
dxy parameter. (This is the transverse impact parameter w.r.t. to (0,0,0) ONLY if refPoint is close t...
Definition: TrackBase.h:630
double vx() const
x coordinate of the reference point on track
Definition: TrackBase.h:696
void PackedCandidateTrackValidator::bookHistograms ( DQMStore::IBooker iBooker,
edm::Run const &  ,
edm::EventSetup const &   
)
override

Definition at line 560 of file PackedCandidateTrackValidator.cc.

References DQMStore::IBooker::book1D(), constexpr, h_diffCharge, h_diffCovDszDsz, h_diffCovDxyDsz, h_diffCovDxyDxy, h_diffCovLambdaDsz, h_diffCovLambdaLambda, h_diffCovPhiDxy, h_diffCovPhiPhi, h_diffCovQoverpQoverp, h_diffDszError, h_diffDxyAssocPV, h_diffDxyError, h_diffDxyPV, h_diffDzAssocPV, h_diffDzError, h_diffDzPV, h_diffEta, h_diffEtaError, h_diffHitPatternHasValidHitInFirstPixelBarrel, h_diffHitPatternNumberOfLostInnerHits, h_diffHitPatternNumberOfValidHits, h_diffHitPatternNumberOfValidPixelHits, h_diffHitPatternPixelLayersWithMeasurement, h_diffHitPatternStripLayersWithMeasurement, h_diffHitPatternTrackerLayersWithMeasurement, h_diffIsHighPurity, h_diffLostInnerHits, h_diffNdof, h_diffNormalizedChi2, h_diffNumberOfHits, h_diffNumberOfPixelHits, h_diffNumberOfPixelLayers, h_diffNumberOfStripLayers, h_diffPhi, h_diffPt, h_diffPtError, h_diffTrackDxyAssocPV, h_diffTrackDxyError, h_diffTrackDzAssocPV, h_diffTrackDzError, h_diffVx, h_diffVy, h_diffVz, h_numberHitsOverMax, h_numberLayersOverMax, h_numberPixelHitsOverMax, h_numberPixelLayersOverMax, h_numberStripHitsOverMax, h_numberStripLayersOverMax, h_selectionFlow, rootFolder_, MonitorElement::setBinLabel(), and DQMStore::IBooker::setCurrentFolder().

560  {
561  iBooker.setCurrentFolder(rootFolder_);
562 
563  h_selectionFlow = iBooker.book1D("selectionFlow", "Track selection flow", 8, -0.5, 7.5);
564  h_selectionFlow->setBinLabel(1, "All tracks");
565  h_selectionFlow->setBinLabel(2, "Associated to PackedCandidate");
566  h_selectionFlow->setBinLabel(3, "PC is charged"),
567  h_selectionFlow->setBinLabel(4, "PC has track");
568  h_selectionFlow->setBinLabel(5, "PC is not electron");
569  h_selectionFlow->setBinLabel(6, "PC has hits");
570  h_selectionFlow->setBinLabel(7, "PC ndof != 0");
571  h_selectionFlow->setBinLabel(8, "Track: no missing inner hits");
572 
573  constexpr int diffBins = 50;
574 
575  h_diffVx = iBooker.book1D("diffVx", "PackedCandidate::bestTrack() - reco::Track in vx()", diffBins, -0.2, 0.2); // not equal
576  h_diffVy = iBooker.book1D("diffVy", "PackedCandidate::bestTrack() - reco::Track in vy()", diffBins, -0.2, 0.2); // not equal
577  h_diffVz = iBooker.book1D("diffVz", "PackedCandidate::bestTrack() - reco::Track in vz()", diffBins, -0.4, 0.4); // not equal
578 
579  h_diffNormalizedChi2 = iBooker.book1D("diffNormalizedChi2", "PackedCandidate::bestTrack() - reco::Track in normalizedChi2()", 30, -1.5, 1.5); // expected difference in -1...0
580  h_diffNdof = iBooker.book1D("diffNdof", "PackedCandidate::bestTrack() - reco::Track in ndof()", 33, -30.5, 2.5); // to monitor the difference
581 
582  h_diffCharge = iBooker.book1D("diffCharge", "PackedCandidate::bestTrack() - reco::Track in charge()", 5, -2.5, 2.5); // expect equality
583  h_diffIsHighPurity = iBooker.book1D("diffIsHighPurity", "PackedCandidate::bestTrack() - reco::Track in quality(highPurity)", 3, -1.5, 1.5); // expect equality
584 
585  h_diffPt = iBooker.book1D("diffPt", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in pt()", diffBins, -1.1, 1.1); // not equal, keep
586  h_diffEta = iBooker.book1D("diffEta", "PackedCandidate::bestTrack() - reco::Track in eta()", diffBins, -0.001, 0.001); // not equal, keep
587  h_diffPhi = iBooker.book1D("diffPhi", "PackedCandidate::bestTrack() - reco::Track in phi()", diffBins, -0.0005, 0.0005); // expect equality within precision
588 
589  h_diffDxyAssocPV.book(iBooker, "diffDxyAssocPV", "(PackedCandidate::dxy() - reco::Track::dxy(assocPV))/reco::Track",
590  40, -0.001, 0.001, // expect equality within precision
591  50, -0.5, 0.5);
592  h_diffDzAssocPV.book(iBooker, "diffDzAssocPV","(PackedCandidate::dzAssociatedPV() - reco::Track::dz(assocPV))/reco::Track",
593  40, -0.001, 0.001, // expect equality within precision
594  50, -0.5, 0.5);
595  h_diffDxyPV = iBooker.book1D("diffDxyPV", "(PackedCandidate::dxy(PV) - reco::Track::dxy(PV))/reco::Track", diffBins, -0.01, 0.01); // expect equality within precision (worse than assocPV)
596  h_diffDzPV = iBooker.book1D("diffDzPV", "(PackedCandidate::dz(PV) - reco::Track::dz(PV))/reco::Track", diffBins, -0.01, 0.01); // expect equality wihtin precision (worse than assocPV)
597  h_diffTrackDxyAssocPV = iBooker.book1D("diffTrackDxyAssocPV", "(PackedCandidate::bestTrack()::dxy(assocPV)) - reco::Track::dxy(assocPV))/reco::Track", diffBins, -0.01, 0.01); // not equal
598  h_diffTrackDzAssocPV = iBooker.book1D("diffTrackDzAssocPV", "(PackedCandidate::bestTrack()::dz(assocPV)) - reco::Track::dz(assocPV))/reco::Track", diffBins, -0.01, 0.01); // not equal
599 
600  h_diffCovQoverpQoverp.book(iBooker, "diffCovQoverpQoverp", "(PackedCandidate::bestTrack() - reco::Track)/reco::track in cov(qoverp, qoverp)",
601  40, -0.05, 0.15, // expect equality within precision (worst precision is exp(1/128*15) =~ 12 %
602  50, -0.5, 0.5);
603  h_diffCovLambdaLambda.book(iBooker, "diffCovLambdaLambda", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in cov(lambda, lambda)",
604  40, -0.05, 0.15, // expect equality within precision worst precision is exp(1/128*(20-5)) =~ 12 % (multiplied by pt^2 in packing & unpacking)
605  50, -0.5, 0.5);
606  h_diffCovLambdaDsz.book(iBooker, "diffCovLambdaDsz", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in cov(lambda, dsz)",
607  60, -0.15, 0.15, // expect equality within precision, worst precision is exp(1/128*(17-4) =~ 11 %
608  50, -1, 1);
609  h_diffCovPhiPhi.book(iBooker, "diffCovPhiPhi", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in cov(phi, phi)",
610  40, -0.05, 0.15, // expect equality within precision worst precision is exp(1/128*(20-5)) =~ 12 % (multiplied by pt^2 in packing & unpacking)
611  50, -0.5, 0.5);
612  h_diffCovPhiDxy.book(iBooker, "diffCovPhiDxy", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in cov(phi, dxy)",
613  60, -0.15, 0.15, // expect equality within precision, wors precision is exp(1/128)*(17-4) =~ 11 %
614  50, -1, 1);
615  h_diffCovDxyDxy.book(iBooker, "diffCovDxyDxy", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in cov(dxy, dxy)",
616  40, -0.001, 0.001,
617  50, -0.1, 0.1);
618  h_diffCovDxyDsz.book(iBooker, "diffCovDxyDsz", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in cov(dxy, dsz)",
619  40, -0.001, 0.001, // expect equality within precision
620  50, -0.5, 0.5);
621  h_diffCovDszDsz.book(iBooker, "diffCovDszDsz", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in cov(dsz, dsz)",
622  40, -0.001, 0.001, // expect equality within precision
623  50, -0.1, 0.1);
624 
625  h_diffDxyError = iBooker.book1D("diffDxyError", "(PackedCandidate::dxyError() - reco::Track::dxyError())/reco::Track", 40, -0.001, 0.001); // expect equality within precision
626  h_diffDszError = iBooker.book1D("diffDszError", "(PackedCandidate::dzError() - reco::Track::dszError())/reco::Track", 40, -0.001, 0.001); // ideally, not equal, but for now they are
627  h_diffDzError = iBooker.book1D("diffDzError", "(PackedCandidate::dzError() - reco::Track::dzError())/reco::Track", 40, -0.001, 0.001); // expect equality within precision (not currently the case)
628 
629  h_diffTrackDxyError = iBooker.book1D("diffTrackDxyError", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in dxyError()", 40, -0.001, 0.001); // expect equality within precision
630  h_diffTrackDzError = iBooker.book1D("diffTrackDzError", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in dzError()", 40, -0.05, 0.05); // not equal
631 
632  h_diffPtError = iBooker.book1D("diffPtError", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in ptError()", diffBins, -1.1, 1.1); // not equal
633  h_diffEtaError = iBooker.book1D("diffEtaError", "(PackedCandidate::bestTrack() - reco::Track)/reco::Track in etaError()", 60, -0.15, 0.15); // not equal
634 
635 
636  h_diffNumberOfPixelLayers = iBooker.book1D("diffNumberOfPixelLayers", "PackedCandidate::pixelLayersWithMeasurement() - reco::Track::hitPattern::pixelLayersWithMeasurement()", 5, -2.5, 2.5); // expect equality
637  h_diffNumberOfStripLayers = iBooker.book1D("diffNumberOfStripLayers", "PackedCandidate::stripLayersWithMeasurement() - reco::Track::hitPattern::stripLayersWithMeasurement()", 5, -2.5, 2.5); // expect equality
638  h_diffNumberOfPixelHits = iBooker.book1D("diffNumberOfPixelHits", "PackedCandidate::numberOfPixelHits() - reco::Track::hitPattern::numberOfValidPixelHits()", 5, -2.5, 2.5); // expect equality
639  h_diffNumberOfHits = iBooker.book1D("diffNumberOfHits", "PackedCandidate::numberHits() - reco::Track::hitPattern::numberOfValidHits()", 5, -2.5, 2.5); // expect equality
640  h_diffLostInnerHits = iBooker.book1D("diffLostInnerHits", "PackedCandidate::lostInnerHits() - reco::Track::hitPattern::numberOfLostHits(MISSING_INNER_HITS)", 5, -2.5, 2.5); // expect equality
641 
642  h_diffHitPatternPixelLayersWithMeasurement = iBooker.book1D("diffHitPatternPixelLayersWithMeasurement", "PackedCandidate::bestTrack() - reco::Track in hitPattern::pixelLayersWithMeasurement()", 13, -10.5, 2.5); // not equal
643  h_diffHitPatternStripLayersWithMeasurement = iBooker.book1D("diffHitPatternStripLayersWithMeasurement", "PackedCandidate::bestTrack() - reco::Track in hitPattern::stripLayersWithMeasurement()", 13, -10.5, 2.5); // not equal
644  h_diffHitPatternTrackerLayersWithMeasurement = iBooker.book1D("diffHitPatternTrackerLayersWithMeasurement", "PackedCandidate::bestTrack() - reco::Track in hitPattern::trackerLayersWithMeasurement()", 13, -10.5, 2.5); // not equal
645  h_diffHitPatternNumberOfValidPixelHits = iBooker.book1D("diffHitPatternNumberOfValidPixelHits", "PackedCandidate::bestTrack() - reco::Track in hitPattern::numberOfValidPixelHits()", 13, -10.5, 2.5); // not equal
646  h_diffHitPatternNumberOfValidHits = iBooker.book1D("diffHitPatternNumberOfValidHits", "PackedCandidate::bestTrack() - reco::Track in hitPattern::numberOfValidHits()", 13, -10.5, 2.5); // not equal
647  h_diffHitPatternNumberOfLostInnerHits = iBooker.book1D("diffHitPatternNumberOfLostPixelHits", "PackedCandidate::bestTrack() - reco::Track in hitPattern::numberOfLostHits(MISSING_INNER_HITS)", 13, -10.5, 2.5); // not equal
648  h_diffHitPatternHasValidHitInFirstPixelBarrel = iBooker.book1D("diffHitPatternHasValidHitInFirstPixelBarrel", "PackedCandidate::bestTrack() - reco::Track in hitPattern::hasValidHitInFirstPixelBarrel", 3, -1.5, 1.5); // expect equality
649 
650  h_numberPixelLayersOverMax = iBooker.book1D("numberPixelLayersOverMax", "Number of pixel layers over the maximum of PackedCandidate", 10, 0, 10);
651  h_numberStripLayersOverMax = iBooker.book1D("numberStripLayersOverMax", "Number of strip layers over the maximum of PackedCandidate", 10, 0, 10);
652  h_numberLayersOverMax = iBooker.book1D("numberLayersOverMax", "Number of layers over the maximum of PackedCandidate", 20, 0, 20);
653  h_numberPixelHitsOverMax = iBooker.book1D("numberPixelHitsOverMax", "Number of pixel hits over the maximum of PackedCandidate", 10, 0, 10);
654  h_numberStripHitsOverMax = iBooker.book1D("numberStripHitsOverMax", "Number of strip hits over the maximum of PackedCandidate", 10, 0, 10);
655  h_numberHitsOverMax = iBooker.book1D("numberHitsOverMax", "Number of hits over the maximum of PackedCandidate", 20, 0, 20);
656 }
PackedValueCheck< LogIntHelper > h_diffCovLambdaDsz
PackedValueCheck< Float16Helper > h_diffDzAssocPV
PackedValueCheck< Float16Helper > h_diffDxyAssocPV
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
PackedValueCheck< LogIntHelper > h_diffCovLambdaLambda
void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:268
PackedValueCheck< LogIntHelper > h_diffCovPhiDxy
PackedValueCheck< Float16Helper > h_diffCovDxyDxy
PackedValueCheck< LogIntHelper > h_diffCovPhiPhi
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:106
PackedValueCheck< LogIntHelper > h_diffCovQoverpQoverp
PackedValueCheck< Float16Helper > h_diffCovDxyDsz
PackedValueCheck< Float16Helper > h_diffCovDszDsz
#define constexpr
void PackedCandidateTrackValidator::fillDescriptions ( edm::ConfigurationDescriptions descriptions)
static

Definition at line 548 of file PackedCandidateTrackValidator.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::addUntracked(), and AlCaHLTBitMon_QueryRunRegistry::string.

548  {
550 
551  desc.addUntracked<edm::InputTag>("tracks", edm::InputTag("generalTracks"));
552  desc.addUntracked<edm::InputTag>("vertices", edm::InputTag("offlinePrimaryVertices"));
553  desc.addUntracked<edm::InputTag>("trackToPackedCandidateAssociation", edm::InputTag("packedPFCandidates"));
554  desc.addUntracked<std::string>("rootFolder", "Tracking/PackedCandidate");
555  desc.addUntracked<bool>("debug", false);
556 
557  descriptions.add("packedCandidateTrackValidator", desc);
558 }
ParameterDescriptionBase * addUntracked(U const &iLabel, T const &value)
void add(std::string const &label, ParameterSetDescription const &psetDescription)

Member Data Documentation

bool PackedCandidateTrackValidator::debug_
private

Definition at line 451 of file PackedCandidateTrackValidator.cc.

Referenced by analyze().

MonitorElement* PackedCandidateTrackValidator::h_diffCharge
private

Definition at line 472 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<Float16Helper> PackedCandidateTrackValidator::h_diffCovDszDsz
private

Definition at line 493 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<Float16Helper> PackedCandidateTrackValidator::h_diffCovDxyDsz
private

Definition at line 492 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<Float16Helper> PackedCandidateTrackValidator::h_diffCovDxyDxy
private

Definition at line 491 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<LogIntHelper> PackedCandidateTrackValidator::h_diffCovLambdaDsz
private

Definition at line 488 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<LogIntHelper> PackedCandidateTrackValidator::h_diffCovLambdaLambda
private

Definition at line 487 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<LogIntHelper> PackedCandidateTrackValidator::h_diffCovPhiDxy
private

Definition at line 490 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<LogIntHelper> PackedCandidateTrackValidator::h_diffCovPhiPhi
private

Definition at line 489 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<LogIntHelper> PackedCandidateTrackValidator::h_diffCovQoverpQoverp
private

Definition at line 486 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffDszError
private

Definition at line 496 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<Float16Helper> PackedCandidateTrackValidator::h_diffDxyAssocPV
private

Definition at line 478 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffDxyError
private

Definition at line 495 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffDxyPV
private

Definition at line 480 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

PackedValueCheck<Float16Helper> PackedCandidateTrackValidator::h_diffDzAssocPV
private

Definition at line 479 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffDzError
private

Definition at line 497 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffDzPV
private

Definition at line 481 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffEta
private

Definition at line 476 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffEtaError
private

Definition at line 503 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffHitPatternHasValidHitInFirstPixelBarrel
private

Definition at line 518 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffHitPatternNumberOfLostInnerHits
private

Definition at line 517 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffHitPatternNumberOfValidHits
private

Definition at line 516 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffHitPatternNumberOfValidPixelHits
private

Definition at line 515 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffHitPatternPixelLayersWithMeasurement
private

Definition at line 512 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffHitPatternStripLayersWithMeasurement
private

Definition at line 514 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffHitPatternTrackerLayersWithMeasurement
private

Definition at line 513 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffIsHighPurity
private

Definition at line 473 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffLostInnerHits
private

Definition at line 510 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffNdof
private

Definition at line 470 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffNormalizedChi2
private

Definition at line 469 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffNumberOfHits
private

Definition at line 509 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffNumberOfPixelHits
private

Definition at line 508 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffNumberOfPixelLayers
private

Definition at line 506 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffNumberOfStripLayers
private

Definition at line 507 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffPhi
private

Definition at line 477 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffPt
private

Definition at line 475 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffPtError
private

Definition at line 502 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffTrackDxyAssocPV
private

Definition at line 483 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffTrackDxyError
private

Definition at line 499 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffTrackDzAssocPV
private

Definition at line 484 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffTrackDzError
private

Definition at line 500 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffVx
private

Definition at line 465 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffVy
private

Definition at line 466 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_diffVz
private

Definition at line 467 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_numberHitsOverMax
private

Definition at line 525 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_numberLayersOverMax
private

Definition at line 522 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_numberPixelHitsOverMax
private

Definition at line 523 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_numberPixelLayersOverMax
private

Definition at line 520 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_numberStripHitsOverMax
private

Definition at line 524 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_numberStripLayersOverMax
private

Definition at line 521 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* PackedCandidateTrackValidator::h_selectionFlow
private

Definition at line 463 of file PackedCandidateTrackValidator.cc.

Referenced by analyze(), and bookHistograms().

std::string PackedCandidateTrackValidator::rootFolder_
private

Definition at line 450 of file PackedCandidateTrackValidator.cc.

Referenced by bookHistograms().

edm::EDGetTokenT<edm::View<reco::Track> > PackedCandidateTrackValidator::tracksToken_
private

Definition at line 446 of file PackedCandidateTrackValidator.cc.

Referenced by analyze().

edm::EDGetTokenT<edm::Association<pat::PackedCandidateCollection> > PackedCandidateTrackValidator::trackToPackedCandidateToken_
private

Definition at line 448 of file PackedCandidateTrackValidator.cc.

Referenced by analyze().

edm::EDGetTokenT<reco::VertexCollection> PackedCandidateTrackValidator::verticesToken_
private

Definition at line 447 of file PackedCandidateTrackValidator.cc.

Referenced by analyze().