MuonTrackValidator Class Reference

#include <MuonTrackValidator.h>

Inheritance diagram for MuonTrackValidator:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	Method called before the event loop. More...
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	Method called at the end of the event loop. More...
	MuonTrackValidator (const edm::ParameterSet &pset)
	Constructor. More...
	~MuonTrackValidator () override
	Destructor. More...
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Private Attributes
edm::EDGetTokenT< SimHitTPAssociationProducer::SimHitTPAssociationList >	_simHitTpMapTag
edm::InputTag	associatormap
bool	BiDirectional_RecoToSim_association
CosmicTrackingParticleSelector	cosmictpSelector
std::string	dirName_
bool	MABH
edm::EDGetTokenT< reco::RecoToSimCollection >	recoToSimCollection_Token
edm::EDGetTokenT< reco::SimToRecoCollection >	simToRecoCollection_Token
TrackingParticleSelector	tpSelector
bool	UseAssociators
bool	useGEMs_
bool	useME0_

Additional Inherited Members
Protected Member Functions inherited from MuonTrackValidatorBase
void	BinLogX (TH1 *h)
void	doProfileX (TH2 *th2, MonitorElement *me)
void	doProfileX (MonitorElement *th2m, MonitorElement *me)
template<typename T >
void	fillPlotNoFlow (MonitorElement *h, T val)
double	getEta (double eta)
double	getPt (double pt)
	MuonTrackValidatorBase (const edm::ParameterSet &pset, edm::ConsumesCollector iC)
	Constructor. More...
	MuonTrackValidatorBase (const edm::ParameterSet &pset)
virtual	~MuonTrackValidatorBase ()(false)
	Destructor. More...
Protected Attributes inherited from MuonTrackValidatorBase
std::vector< std::string >	associators
edm::InputTag	bsSrc
edm::EDGetTokenT< reco::BeamSpot >	bsSrc_Token
std::vector< MonitorElement * >	chi2_vs_eta
std::vector< MonitorElement * >	chi2_vs_nhits
std::vector< MonitorElement * >	chi2_vs_phi
int	cotThetaRes_nbin
double	cotThetaRes_rangeMax
double	cotThetaRes_rangeMin
bool	do_MUOhitsPlots
bool	do_TRKhitsPlots
std::vector< MonitorElement * >	dxypull_vs_eta
int	dxyRes_nbin
double	dxyRes_rangeMax
double	dxyRes_rangeMin
std::vector< MonitorElement * >	dxyres_vs_eta
std::vector< MonitorElement * >	dxyres_vs_pt
std::vector< MonitorElement * >	dzpull_vs_eta
int	dzRes_nbin
double	dzRes_rangeMax
double	dzRes_rangeMin
std::vector< MonitorElement * >	dzres_vs_eta
std::vector< MonitorElement * >	dzres_vs_pt
int	etaRes_nbin
double	etaRes_rangeMax
double	etaRes_rangeMin
std::vector< MonitorElement * >	etares_vs_eta
std::vector< MonitorElement * >	h_assoc2dxy
std::vector< MonitorElement * >	h_assoc2dz
std::vector< MonitorElement * >	h_assoc2eta
std::vector< MonitorElement * >	h_assoc2hit
std::vector< MonitorElement * >	h_assoc2phi
std::vector< MonitorElement * >	h_assoc2pT
std::vector< MonitorElement * >	h_assoc2pu
std::vector< MonitorElement * >	h_assocdxy
std::vector< MonitorElement * >	h_assocdz
std::vector< MonitorElement * >	h_assoceta
std::vector< MonitorElement * >	h_assocFraction
std::vector< MonitorElement * >	h_assochi2
std::vector< MonitorElement * >	h_assochi2_prob
std::vector< MonitorElement * >	h_assochit
std::vector< MonitorElement * >	h_assocphi
std::vector< MonitorElement * >	h_assocpT
std::vector< MonitorElement * >	h_assocpu
std::vector< MonitorElement * >	h_assocRpos
std::vector< MonitorElement * >	h_assocSharedHit
std::vector< MonitorElement * >	h_assocZpos
std::vector< MonitorElement * >	h_charge
std::vector< MonitorElement * >	h_dxypull
std::vector< MonitorElement * >	h_dxypulleta
std::vector< MonitorElement * >	h_dzpull
std::vector< MonitorElement * >	h_dzpulleta
std::vector< MonitorElement * >	h_etaRes
std::vector< MonitorElement * >	h_fakes
std::vector< MonitorElement * >	h_losthits
std::vector< MonitorElement * >	h_misiddxy
std::vector< MonitorElement * >	h_misiddz
std::vector< MonitorElement * >	h_misideta
std::vector< MonitorElement * >	h_misidhit
std::vector< MonitorElement * >	h_misidphi
std::vector< MonitorElement * >	h_misidpT
std::vector< MonitorElement * >	h_misidpu
std::vector< MonitorElement * >	h_nchi2
std::vector< MonitorElement * >	h_nchi2_prob
std::vector< MonitorElement * >	h_nhits
std::vector< MonitorElement * >	h_nlosthits
std::vector< MonitorElement * >	h_nmisslayers_inner
std::vector< MonitorElement * >	h_nmisslayers_outer
std::vector< MonitorElement * >	h_phipull
std::vector< MonitorElement * >	h_phipulleta
std::vector< MonitorElement * >	h_phipullphi
std::vector< MonitorElement * >	h_ptpull
std::vector< MonitorElement * >	h_ptpulleta
std::vector< MonitorElement * >	h_ptpullphi
std::vector< MonitorElement * >	h_PurityVsQuality
std::vector< MonitorElement * >	h_qoverppull
std::vector< MonitorElement * >	h_recodxy
std::vector< MonitorElement * >	h_recodz
std::vector< MonitorElement * >	h_recoeta
std::vector< MonitorElement * >	h_recohit
std::vector< MonitorElement * >	h_recophi
std::vector< MonitorElement * >	h_recopT
std::vector< MonitorElement * >	h_recopu
std::vector< MonitorElement * >	h_simuldxy
std::vector< MonitorElement * >	h_simuldz
std::vector< MonitorElement * >	h_simuleta
std::vector< MonitorElement * >	h_simulhit
std::vector< MonitorElement * >	h_simulphi
std::vector< MonitorElement * >	h_simulpT
std::vector< MonitorElement * >	h_simulpu
std::vector< MonitorElement * >	h_simulRpos
std::vector< MonitorElement * >	h_simulZpos
std::vector< MonitorElement * >	h_thetapull
std::vector< MonitorElement * >	h_thetapulleta
std::vector< MonitorElement * >	h_thetapullphi
std::vector< MonitorElement * >	h_tracks
bool	ignoremissingtkcollection_
std::vector< edm::InputTag >	label
edm::InputTag	label_pileupinfo
edm::InputTag	label_tp_effic
edm::InputTag	label_tp_fake
double	maxCSCHit
double	maxDTHit
double	maxDxy
double	maxDz
double	maxEta
double	maxLayers
double	maxNHit
double	maxPhi
double	maxPixels
double	maxPt
double	maxPU
double	maxRPCHit
double	maxRpos
double	maxZpos
double	minCSCHit
double	minDTHit
double	minDxy
double	minDz
double	minEta
double	minLayers
double	minNHit
double	minPhi
double	minPixels
double	minPt
double	minPU
double	minRPCHit
double	minRpos
double	minZpos
edm::ParameterSet	muonHistoParameters
std::vector< MonitorElement * >	nCSChits_vs_eta
std::vector< MonitorElement * >	nDThits_vs_eta
std::vector< MonitorElement * >	nGEMhits_vs_eta
std::vector< MonitorElement * >	nhits_vs_eta
std::vector< MonitorElement * >	nhits_vs_phi
int	nintCSCHit
int	nintDTHit
int	nintDxy
int	nintDz
int	nintEta
int	nintLayers
int	nintNHit
int	nintPhi
int	nintPixels
int	nintPt
int	nintPU
int	nintRPCHit
int	nintRpos
int	nintZpos
std::vector< MonitorElement * >	nlosthits_vs_eta
std::vector< MonitorElement * >	nME0hits_vs_eta
std::vector< MonitorElement * >	nPixel_LayersWithMeas_vs_eta
std::vector< MonitorElement * >	nRecHits_vs_nSimHits
std::vector< MonitorElement * >	nRPChits_vs_eta
std::vector< MonitorElement * >	nTRK_LayersWithMeas_vs_eta
std::string	out
std::string	parametersDefiner
std::vector< MonitorElement * >	phipull_vs_eta
std::vector< MonitorElement * >	phipull_vs_phi
int	phiRes_nbin
double	phiRes_rangeMax
double	phiRes_rangeMin
std::vector< MonitorElement * >	phires_vs_eta
std::vector< MonitorElement * >	phires_vs_phi
std::vector< MonitorElement * >	phires_vs_pt
edm::EDGetTokenT< std::vector< PileupSummaryInfo > >	pileupinfo_Token
std::vector< MonitorElement * >	ptpull_vs_eta
std::vector< MonitorElement * >	ptpull_vs_phi
int	ptRes_nbin
double	ptRes_rangeMax
double	ptRes_rangeMin
std::vector< MonitorElement * >	ptres_vs_eta
std::vector< MonitorElement * >	ptres_vs_phi
std::vector< MonitorElement * >	ptres_vs_pt
edm::ESHandle< MagneticField >	theMF
std::vector< MonitorElement * >	thetaCotres_vs_eta
std::vector< MonitorElement * >	thetaCotres_vs_pt
std::vector< MonitorElement * >	thetapull_vs_eta
std::vector< MonitorElement * >	thetapull_vs_phi
edm::EDGetTokenT< TrackingParticleCollection >	tp_effic_Token
edm::EDGetTokenT< TrackingParticleCollection >	tp_fake_Token
std::vector< edm::EDGetTokenT< edm::View< reco::Track > > >	track_Collection_Token
bool	useFabsEta
bool	useInvPt
bool	useLogPt
bool	usemuon
bool	usetracker

Detailed Description

Class that produces histograms to validate Muon Track Reconstruction performances

Definition at line 17 of file MuonTrackValidator.h.

Constructor & Destructor Documentation

MuonTrackValidator::MuonTrackValidator ( const edm::ParameterSet &  pset )

inline

Constructor.

Definition at line 20 of file MuonTrackValidator.h.

References _simHitTpMapTag, ctfWithMaterialTrackMCMatch_cfi::associator, associatormap, MuonTrackValidatorBase::associators, BiDirectional_RecoToSim_association, MuonTrackValidatorBase::bsSrc, MuonTrackValidatorBase::bsSrc_Token, cosmictpSelector, dirName_, edm::ParameterSet::dump(), edm::ParameterSet::getParameter(), MuonTrackValidatorBase::label, MuonTrackValidatorBase::label_pileupinfo, MuonTrackValidatorBase::label_tp_effic, MuonTrackValidatorBase::label_tp_fake, MABH, MuonTrackValidatorBase::pileupinfo_Token, recoToSimCollection_Token, simToRecoCollection_Token, AlCaHLTBitMon_QueryRunRegistry::string, MuonTrackValidatorBase::tp_effic_Token, MuonTrackValidatorBase::tp_fake_Token, tpSelector, MuonTrackValidatorBase::track_Collection_Token, UseAssociators, useGEMs_, useME0_, MuonTrackValidatorBase::usemuon, and MuonTrackValidatorBase::usetracker.

                                                 :MuonTrackValidatorBase(pset){
    dirName_ = pset.getParameter<std::string>("dirName");
    associatormap = pset.getParameter< edm::InputTag >("associatormap");
    UseAssociators = pset.getParameter< bool >("UseAssociators");
    useGEMs_ = pset.getParameter< bool >("useGEMs");
    useME0_ = pset.getParameter< bool >("useME0");
    edm::ParameterSet tpset = pset.getParameter<edm::ParameterSet>("muonTPSelector");
    tpSelector = TrackingParticleSelector(tpset.getParameter<double>("ptMin"),
                                          tpset.getParameter<double>("ptMax"),
                      tpset.getParameter<double>("minRapidity"),
                      tpset.getParameter<double>("maxRapidity"),
                      tpset.getParameter<double>("tip"),
                      tpset.getParameter<double>("lip"),
                      tpset.getParameter<int>("minHit"),
                      tpset.getParameter<bool>("signalOnly"),
                      tpset.getParameter<bool>("intimeOnly"),
                      tpset.getParameter<bool>("chargedOnly"),
                      tpset.getParameter<bool>("stableOnly"),
                      tpset.getParameter<std::vector<int> >("pdgId"));

    cosmictpSelector = CosmicTrackingParticleSelector(tpset.getParameter<double>("ptMin"),
                              tpset.getParameter<double>("minRapidity"),
                              tpset.getParameter<double>("maxRapidity"),
                              tpset.getParameter<double>("tip"),
                              tpset.getParameter<double>("lip"),
                              tpset.getParameter<int>("minHit"),
                              tpset.getParameter<bool>("chargedOnly"),
                              tpset.getParameter<std::vector<int> >("pdgId"));
    
    BiDirectional_RecoToSim_association = pset.getParameter<bool>("BiDirectional_RecoToSim_association");

    // dump cfg parameters
    edm::LogVerbatim("MuonTrackValidator") << "constructing MuonTrackValidator: " << pset.dump();
    
    // Declare consumes (also for the base class)
    bsSrc_Token = consumes<reco::BeamSpot>(bsSrc);
    tp_effic_Token = consumes<TrackingParticleCollection>(label_tp_effic);
    tp_fake_Token = consumes<TrackingParticleCollection>(label_tp_fake);
    pileupinfo_Token = consumes<std::vector<PileupSummaryInfo> >(label_pileupinfo);
    for (unsigned int www=0;www<label.size();www++){
      track_Collection_Token.push_back(consumes<edm::View<reco::Track> >(label[www]));
    }
    simToRecoCollection_Token = consumes<reco::SimToRecoCollection>(associatormap);
    recoToSimCollection_Token = consumes<reco::RecoToSimCollection>(associatormap);

    _simHitTpMapTag = mayConsume<SimHitTPAssociationProducer::SimHitTPAssociationList>(pset.getParameter<edm::InputTag>("simHitTpMapTag"));

    MABH = false;
    if (!UseAssociators) {
      // flag MuonAssociatorByHits
      if (associators[0] == "MuonAssociationByHits") MABH = true;
      // reset string associators to the map label
      associators.clear();
      associators.push_back(associatormap.label());
      edm::LogVerbatim("MuonTrackValidator") << "--> associators reset to: " <<associators[0];
    } else {
      for (auto const& associator :associators) {
        consumes<reco::TrackToTrackingParticleAssociator>(edm::InputTag(associator));
      }
    }
    
    // inform on which SimHits will be counted
    if (usetracker) edm::LogVerbatim("MuonTrackValidator")
      <<"\n usetracker = TRUE : Tracker SimHits WILL be counted";
    else edm::LogVerbatim("MuonTrackValidator")
      <<"\n usetracker = FALSE : Tracker SimHits WILL NOT be counted";
    if (usemuon) edm::LogVerbatim("MuonTrackValidator")
      <<" usemuon = TRUE : Muon SimHits WILL be counted";
    else edm::LogVerbatim("MuonTrackValidator")
      <<" usemuon = FALSE : Muon SimHits WILL NOT be counted"<<std::endl;
    
    // loop over the reco::Track collections to validate: check for inconsistent input settings
    for (unsigned int www=0;www<label.size();www++) {
      std::string recoTracksLabel = label[www].label();
      std::string recoTracksInstance = label[www].instance();
      
      // tracks with hits only on tracker
      if (recoTracksLabel=="generalTracks" ||
      (recoTracksLabel.find("cutsRecoTracks") != std::string::npos) ||
      recoTracksLabel=="ctfWithMaterialTracksP5LHCNavigation" ||
      recoTracksLabel=="hltL3TkTracksFromL2" ||
      (recoTracksLabel=="hltL3Muons" && recoTracksInstance=="L2Seeded"))
    {
      if (usemuon) {
        edm::LogWarning("MuonTrackValidator")
          <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
          <<"\n with usemuon == true"<<"\n ---> please change to usemuon == false ";
      }
      if (!usetracker) {
        edm::LogWarning("MuonTrackValidator")
          <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
          <<"\n with usetracker == false"<<"\n ---> please change to usetracker == true ";
      } 
    }
      
      // tracks with hits only on muon detectors
      else if (recoTracksLabel=="standAloneMuons" ||
           recoTracksLabel=="standAloneSETMuons" ||
           recoTracksLabel=="cosmicMuons" ||
           recoTracksLabel=="hltL2Muons")
    {
      if (usetracker) {
        edm::LogWarning("MuonTrackValidator")
          <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
          <<"\n with usetracker == true"<<"\n ---> please change to usetracker == false ";
      }
      if (!usemuon) {
        edm::LogWarning("MuonTrackValidator")
          <<"\n*** WARNING : inconsistent input tracksTag = "<<label[www]
          <<"\n with usemuon == false"<<"\n ---> please change to usemuon == true ";
      }
    }
      
    } // for (unsigned int www=0;www<label.size();www++)
  }

MuonTrackValidator::~MuonTrackValidator ( )

inlineoverride

Destructor.

Definition at line 137 of file MuonTrackValidator.h.

References analyze(), and bookHistograms().

{ }

Member Function Documentation

void MuonTrackValidator::analyze ( const edm::Event &  event, const edm::EventSetup &  setup  )

override

Method called before the event loop.

Method called once per event

Definition at line 217 of file MuonTrackValidator.cc.

References funct::abs(), ctfWithMaterialTrackMCMatch_cfi::associator, muonValidationHeavyIons_cff::associatormap, MuonErrorMatrixAnalyzer_cfi::associatorName, simPFProducer_cfi::associators, edm::AssociationMap< Tag >::begin(), MuonTrackValidator_cfi::BiDirectional_RecoToSim_association, reco::TrackBase::charge(), reco::TrackBase::chi2(), CosmicParametersDefinerForTP::clone(), funct::cos(), histoParameters_cff::do_MUOhitsPlots, histoParameters_cff::do_TRKhitsPlots, TrackingParticleIP::dxy(), reco::TrackBase::dxy(), reco::TrackBase::dxyError(), TrackingParticleIP::dz(), reco::TrackBase::dz(), reco::TrackBase::dzError(), edm::AssociationMap< Tag >::end(), reco::TrackBase::eta(), edm::AssociationMap< Tag >::find(), dedxEstimators_cff::fraction, edm::EventSetup::get(), edm::Ref< C, T, F >::get(), edm::Event::getByToken(), reco::TrackBase::hitPattern(), mps_fire::i, createfilelist::int, edm::RefToBase< T >::key(), edm::Ref< C, T, F >::key(), label, LogTrace, M_PI, associators_cff::MABH, reco::HitPattern::MISSING_INNER_HITS, reco::HitPattern::MISSING_OUTER_HITS, ParametersDefinerForTP::momentum(), TrackingParticle::momentum(), reco::TrackBase::ndof(), reco::TrackBase::normalizedChi2(), reco::HitPattern::numberOfLostHits(), reco::TrackBase::numberOfLostHits(), reco::TrackBase::numberOfValidHits(), reco::HitPattern::numberOfValidMuonCSCHits(), reco::HitPattern::numberOfValidMuonDTHits(), reco::HitPattern::numberOfValidMuonGEMHits(), reco::HitPattern::numberOfValidMuonME0Hits(), reco::HitPattern::numberOfValidMuonRPCHits(), MuonTrackValidator_cfi::parametersDefiner, reco::TrackBase::phi(), muonME0PseudoDigis_cfi::phiError, reco::TrackBase::phiError(), reco::HitPattern::pixelLayersWithMeasurement(), reco::BeamSpot::position(), edm::AssociationMap< Tag >::post_insert(), edm::Handle< T >::product(), reco::TrackBase::pt(), reco::TrackBase::ptError(), reco::TrackBase::qoverp(), reco::TrackBase::qoverpError(), jets_cff::quality, funct::sin(), mathSSE::sqrt(), reco::TrackBase::theta(), reco::TrackBase::thetaError(), RecoMuonValidator_cfi::tpSelector, HiIsolationCommonParameters_cff::track, findElectronsInSiStrips_cfi::trackCollection, reco::HitPattern::trackerLayersWithMeasurement(), TrackAlgoCompareUtil_cfi::UseAssociators, histoParameters_cff::usemuon, histoParameters_cff::usetracker, ParametersDefinerForTP::vertex(), TrackingParticle::vertex(), and w.

Referenced by ~MuonTrackValidator().

                                                                                  {
  using namespace reco;

  edm::LogInfo("MuonTrackValidator") << "\n====================================================" << "\n"
                     << "Analyzing new event" << "\n"
                     << "====================================================\n" << "\n";

  edm::Handle< std::vector<PileupSummaryInfo> > puinfoH;
  int PU_NumInteractions(-1);

  edm::ESHandle<ParametersDefinerForTP> Lhc_parametersDefinerTP;
  std::unique_ptr<ParametersDefinerForTP> Cosmic_parametersDefinerTP;

  if(parametersDefiner=="LhcParametersDefinerForTP") {
    setup.get<TrackAssociatorRecord>().get(parametersDefiner, Lhc_parametersDefinerTP);

    // PileupSummaryInfo is contained only in collision events
    event.getByToken(pileupinfo_Token,puinfoH);
    for (std::vector<PileupSummaryInfo>::const_iterator puInfoIt = puinfoH->begin(); puInfoIt != puinfoH->end(); ++puInfoIt) {
      if (puInfoIt->getBunchCrossing()==0) {
    PU_NumInteractions = puInfoIt->getPU_NumInteractions();
    break;
      }
    }

  }
  else if(parametersDefiner=="CosmicParametersDefinerForTP") {
    edm::ESHandle<CosmicParametersDefinerForTP>  _Cosmic_parametersDefinerTP;
    setup.get<TrackAssociatorRecord>().get(parametersDefiner, _Cosmic_parametersDefinerTP);

     //Since we modify the object, we must clone it
    Cosmic_parametersDefinerTP = _Cosmic_parametersDefinerTP->clone();

    edm::Handle<SimHitTPAssociationProducer::SimHitTPAssociationList> simHitsTPAssoc;
    //warning: make sure the TP collection used in the map is the same used here
    event.getByToken(_simHitTpMapTag,simHitsTPAssoc);
    Cosmic_parametersDefinerTP->initEvent(simHitsTPAssoc);
    cosmictpSelector.initEvent(simHitsTPAssoc);
  }
  else {
    edm::LogError("MuonTrackValidator")
      << "Unexpected label: parametersDefiner = "<< parametersDefiner.c_str() << "\n";
  }

  edm::Handle<TrackingParticleCollection> TPCollectionHeff ;
  event.getByToken(tp_effic_Token,TPCollectionHeff);
  TrackingParticleCollection const & tPCeff = *(TPCollectionHeff.product());

  edm::Handle<TrackingParticleCollection> TPCollectionHfake ;
  event.getByToken(tp_fake_Token,TPCollectionHfake);
  
  edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
  event.getByToken(bsSrc_Token,recoBeamSpotHandle);
  reco::BeamSpot bs = *recoBeamSpotHandle;

  std::vector<const reco::TrackToTrackingParticleAssociator*> associator;
  if (UseAssociators) {
    edm::Handle<reco::TrackToTrackingParticleAssociator> theAssociator;
    for (auto const& associatorName : associators) {
      event.getByLabel(associatorName,theAssociator);
      associator.push_back( theAssociator.product() );
    }
  }
  
  int w=0;
  for (unsigned int ww=0;ww<associators.size();ww++) {
    for (unsigned int www=0;www<label.size();www++) {
      //
      //get collections from the event
      //
      edm::Handle<edm::View<Track> > trackCollection;
      unsigned int trackCollectionSize = 0;

      reco::RecoToSimCollection recSimColl;
      reco::SimToRecoCollection simRecColl;

      // account for missing track collections (HLT case)
      if (!event.getByToken(track_Collection_Token[www], trackCollection) && ignoremissingtkcollection_) {
    recSimColl.post_insert();
    simRecColl.post_insert();
      }

      //associate tracks to TrackingParticles
      else {
    trackCollectionSize = trackCollection->size();

    if(UseAssociators) {
      edm::LogVerbatim("MuonTrackValidator") << "Analyzing "
                            << label[www].process()<<":"
                            << label[www].label()<<":"
                            << label[www].instance()<<" with "
                            << associators[ww].c_str() <<"\n";
      
      LogTrace("MuonTrackValidator") << "Calling associateRecoToSim method" << "\n";
      recSimColl = associator[ww]->associateRecoToSim(trackCollection,
                                    TPCollectionHfake);
      LogTrace("MuonTrackValidator") << "Calling associateSimToReco method" << "\n";
      simRecColl = associator[ww]->associateSimToReco(trackCollection,
                                    TPCollectionHeff);
    } else {
      edm::LogVerbatim("MuonTrackValidator") << "Analyzing "
                            << label[www].process()<<":"
                            << label[www].label()<<":"
                            << label[www].instance()<<" with "
                            << associatormap.process()<<":"
                            << associatormap.label()<<":"
                            << associatormap.instance()<<"\n";
      
      Handle<reco::SimToRecoCollection > simtorecoCollectionH;
      event.getByToken(simToRecoCollection_Token,simtorecoCollectionH);
      simRecColl = *simtorecoCollectionH.product();
      
      Handle<reco::RecoToSimCollection > recotosimCollectionH;
      event.getByToken(recoToSimCollection_Token,recotosimCollectionH);
      recSimColl = *recotosimCollectionH.product();
    }   
      }
       
      //
      //fill simulation histograms
      //
      edm::LogVerbatim("MuonTrackValidator") << "\n# of TrackingParticles: " << tPCeff.size() << "\n";
      int ats = 0;
      int st = 0;
      for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
    bool TP_is_matched = false;
        bool isChargeOK = true;
    double quality = 0.;
    
    TrackingParticleRef tpr(TPCollectionHeff, i);
    TrackingParticle* tp = const_cast<TrackingParticle*>(tpr.get());
    
    TrackingParticle::Vector momentumTP;
    TrackingParticle::Point vertexTP;
    double dxySim = 0;
    double dzSim = 0;
    
    //If the TrackingParticle is collision-like, get the momentum and vertex at production state
    //and the impact parameters w.r.t. PCA
    if(parametersDefiner=="LhcParametersDefinerForTP")
      {
        LogTrace("MuonTrackValidator") <<"TrackingParticle "<< i;
        if(! tpSelector(*tp)) continue;
        momentumTP = tp->momentum();
        vertexTP = tp->vertex();
        TrackingParticle::Vector momentum = Lhc_parametersDefinerTP->momentum(event,setup,tpr);
        TrackingParticle::Point vertex = Lhc_parametersDefinerTP->vertex(event,setup,tpr);
        dxySim = TrackingParticleIP::dxy(vertex, momentum, bs.position());
        dzSim = TrackingParticleIP::dz(vertex, momentum, bs.position());
      }
    //for cosmics get the momentum and vertex at PCA
    else if(parametersDefiner=="CosmicParametersDefinerForTP")
      {
        edm::LogVerbatim("MuonTrackValidator") <<"TrackingParticle "<< i;
        if(! cosmictpSelector(tpr,&bs,event,setup)) continue;   
        momentumTP = Cosmic_parametersDefinerTP->momentum(event,setup,tpr);
        vertexTP = Cosmic_parametersDefinerTP->vertex(event,setup,tpr);
        dxySim = TrackingParticleIP::dxy(vertexTP, momentumTP, bs.position());
        dzSim = TrackingParticleIP::dz(vertexTP, momentumTP, bs.position());
      }
    edm::LogVerbatim("MuonTrackValidator") <<"--------------------Selected TrackingParticle #"<<tpr.key();
    edm::LogVerbatim("MuonTrackValidator") <<"momentumTP: pt = "<<sqrt(momentumTP.perp2())<<", pz = "<<momentumTP.z() 
                           <<", \t vertexTP: radius = "<<sqrt(vertexTP.perp2())<<  ",  z = "<<vertexTP.z() <<"\n";
    st++;

    double TPeta = momentumTP.eta();
    double xTPeta = getEta(TPeta);            // may be |eta| in histos according to useFabsEta
        double TPpt = sqrt(momentumTP.perp2());
    double xTPpt = getPt(TPpt);               // may be 1/pt in histos according to useInvPt
        double TPphi = momentumTP.phi();
    double TPrpos = sqrt(vertexTP.perp2());
    double TPzpos = vertexTP.z();

    // Number of counted SimHits depend on the selection of tracker and muon detectors (via cfg parameters) 
    int nSimHits = 0;
    if (usetracker && usemuon) {
      nSimHits= tpr.get()->numberOfHits();
    }
    else if (!usetracker && usemuon) {
      nSimHits= tpr.get()->numberOfHits() - tpr.get()->numberOfTrackerHits();
    }
    else if (usetracker && !usemuon) {
      nSimHits=tpr.get()->numberOfTrackerHits();
    }
    edm::LogVerbatim("MuonTrackValidator") << "\t N simhits = "<< nSimHits<<"\n";

    int assoc_recoTrack_NValidHits = 0;
    if(simRecColl.find(tpr) != simRecColl.end()) {
      auto const & rt = simRecColl[tpr];
      if (!rt.empty()) {
        RefToBase<Track> assoc_recoTrack = rt.begin()->first;
        edm::LogVerbatim("MuonTrackValidator")<<"-----------------------------associated Track #"<<assoc_recoTrack.key();
        TP_is_matched = true;
        ats++;
        if (assoc_recoTrack->charge() != tpr->charge()) isChargeOK = false;
        quality = rt.begin()->second;
        assoc_recoTrack_NValidHits = assoc_recoTrack->numberOfValidHits();
        edm::LogVerbatim("MuonTrackValidator") << "TrackingParticle #" <<tpr.key()
                           << " with pt=" << sqrt(momentumTP.perp2())
                           << " associated with quality:" << quality <<"\n";
      }
    } else {
      edm::LogVerbatim("MuonTrackValidator")
        << "TrackingParticle #" << tpr.key()
        << " with pt,eta,phi: "
        << sqrt(momentumTP.perp2()) << " , "<< momentumTP.eta() << " , "<< momentumTP.phi() << " , "
        << " NOT associated to any reco::Track" << "\n";
    }

    // histos for efficiency vs eta

    fillPlotNoFlow(h_simuleta[w], xTPeta);
    if (TP_is_matched) {
      fillPlotNoFlow(h_assoceta[w], xTPeta);
      if (!isChargeOK) fillPlotNoFlow(h_misideta[w], xTPeta);
    }
      
    // histos for efficiency vs phi 
    fillPlotNoFlow(h_simulphi[w], TPphi);
    if (TP_is_matched) {
      fillPlotNoFlow(h_assocphi[w], TPphi);
      if (!isChargeOK) fillPlotNoFlow(h_misidphi[w], TPphi);
    }
    
    // histos for efficiency vs pT
    fillPlotNoFlow(h_simulpT[w], xTPpt);
    if (TP_is_matched) {
      fillPlotNoFlow(h_assocpT[w], xTPpt);
      if (!isChargeOK) fillPlotNoFlow(h_misidpT[w], xTPpt);
    }
    
    // histos for efficiency vs dxy
    fillPlotNoFlow(h_simuldxy[w], dxySim);
    if (TP_is_matched) {
      fillPlotNoFlow(h_assocdxy[w], dxySim);
      if (!isChargeOK) fillPlotNoFlow(h_misiddxy[w], dxySim);
    }

    // histos for efficiency vs dz
    fillPlotNoFlow(h_simuldz[w], dzSim);
    if (TP_is_matched) {
      fillPlotNoFlow(h_assocdz[w], dzSim);
      if (!isChargeOK) fillPlotNoFlow(h_misiddz[w], dzSim);
    }

    // histos for efficiency vs Radius
    fillPlotNoFlow(h_simulRpos[w], TPrpos);
    if (TP_is_matched) fillPlotNoFlow(h_assocRpos[w], TPrpos);
    
    // histos for efficiency vs z position
    fillPlotNoFlow(h_simulZpos[w], TPzpos);
    if (TP_is_matched) fillPlotNoFlow(h_assocZpos[w], TPzpos);

    // histos for efficiency vs Number of Hits
    fillPlotNoFlow(h_simulhit[w], nSimHits );
    if (TP_is_matched) {
      fillPlotNoFlow(h_assochit[w], nSimHits );
      nRecHits_vs_nSimHits[w]->Fill(nSimHits, assoc_recoTrack_NValidHits);

      // charge misid is more useful w.r.t. nRecHits (filled after)
      //if (!isChargeOK) fillPlotNoFlow(h_misidhit[w], nSimHits);
    }
    
    // histos for efficiency vs PileUp
    fillPlotNoFlow(h_simulpu[w], PU_NumInteractions);
    if (TP_is_matched) {
      fillPlotNoFlow(h_assocpu[w], PU_NumInteractions);
      if (!isChargeOK) fillPlotNoFlow(h_misidpu[w], PU_NumInteractions);
    }
    
      } // End for (TrackingParticleCollection::size_type i=0; i<tPCeff.size(); i++){
      
      //
      //fill reconstructed track histograms
      //
      edm::LogVerbatim("MuonTrackValidator") << "\n# of reco::Tracks with "
                         << label[www].process()<<":"
                         << label[www].label()<<":"
                         << label[www].instance()
                         << ": " << trackCollectionSize << "\n";

      int at = 0;
      int rT = 0;
      for(edm::View<Track>::size_type i=0; i<trackCollectionSize; ++i) {
        bool Track_is_matched = false;
        bool isChargeOK = true;
    RefToBase<Track> track(trackCollection, i);
    rT++;
    
    std::vector<std::pair<TrackingParticleRef, double> > tp;
    TrackingParticleRef tpr;

    // new logic (bidirectional)
    if (BiDirectional_RecoToSim_association) {
      edm::LogVerbatim("MuonTrackValidator")<<"----------------------------------------Track #"<< track.key();
      
      if(recSimColl.find(track) != recSimColl.end()) {
        tp = recSimColl[track]; 
        if (!tp.empty()) {
          tpr = tp.begin()->first;  
          // RtS and StR must associate the same pair !
          if(simRecColl.find(tpr) != simRecColl.end()) {
        auto const & assoc_track_checkback = simRecColl[tpr].begin()->first;
        
        if ( assoc_track_checkback.key() == track.key() ) {
          edm::LogVerbatim("MuonTrackValidator")<<"------------------associated TrackingParticle #"<<tpr.key();
          Track_is_matched = true;
          at++;
          if (track->charge() != tpr->charge()) isChargeOK = false;
          double Purity = tp.begin()->second;
          double Quality = simRecColl[tpr].begin()->second;
          edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
                             << " associated with purity:" << Purity <<"\n";
          if (MABH) h_PurityVsQuality[w]->Fill(Quality,Purity);
        }
          }
        }
      }
      
      if (!Track_is_matched)
        edm::LogVerbatim("MuonTrackValidator")
          << "reco::Track #" << track.key() << " with pt=" << track->pt() << " NOT associated to any TrackingParticle" << "\n";
    }
    // old logic, valid for cosmics 2-legs reco (bugged for collision scenario)
    else {
      if(recSimColl.find(track) != recSimColl.end()){
        tp = recSimColl[track];
        if (!tp.empty()) {
          tpr = tp.begin()->first;
          Track_is_matched = true;
          at++;
          if (track->charge() != tpr->charge()) isChargeOK = false;
          edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
                             << " associated with quality:" << tp.begin()->second <<"\n";
        }
      } else {
        edm::LogVerbatim("MuonTrackValidator") << "reco::Track #" << track.key() << " with pt=" << track->pt()
                           << " NOT associated to any TrackingParticle" << "\n";    
      }
    }

    // 
    int nRecHits = track->numberOfValidHits();
    edm::LogVerbatim("MuonTrackValidator") << "\t N valid rechits = "<< nRecHits <<"\n";
    
    double etaRec = track->eta();
    double xetaRec = getEta(etaRec);

    double ptRec = track->pt();
    double xptRec = getPt(ptRec);

    double qoverpRec = track->qoverp();
    double phiRec = track->phi();
    double thetaRec = track->theta();
    double dxyRec = track->dxy(bs.position());
    double dzRec = track->dz(bs.position()); 

    double qoverpError = track->qoverpError();
    double ptError = track->ptError();
        double thetaError = track->thetaError();
    double phiError = track->phiError();
    double dxyError = track->dxyError();
    double dzError = track->dzError();

    // histos for fake rate vs eta
    fillPlotNoFlow(h_recoeta[w], xetaRec);
    if (Track_is_matched) {
      fillPlotNoFlow(h_assoc2eta[w], xetaRec);
    }

    // histos for fake rate vs phi
    fillPlotNoFlow(h_recophi[w], phiRec);
    if (Track_is_matched) {
      fillPlotNoFlow(h_assoc2phi[w], phiRec);
    }
    
    // histos for fake rate vs pT
    fillPlotNoFlow(h_recopT[w], xptRec);
    if (Track_is_matched) {
      fillPlotNoFlow(h_assoc2pT[w], xptRec);
    }
    
    // histos for fake rate vs dxy
    fillPlotNoFlow(h_recodxy[w], dxyRec);
    if (Track_is_matched) {
      fillPlotNoFlow(h_assoc2dxy[w], dxyRec);
    }
    
    // histos for fake rate vs dz
    fillPlotNoFlow(h_recodz[w], dzRec);
    if (Track_is_matched) {
      fillPlotNoFlow(h_assoc2dz[w], dzRec);
    }

    // histos for fake rate vs Number of RecHits in track 
    fillPlotNoFlow(h_recohit[w], nRecHits);
    if (Track_is_matched) {
      fillPlotNoFlow(h_assoc2hit[w], nRecHits);
      // charge misid w.r.t. nRecHits
      if (!isChargeOK) fillPlotNoFlow(h_misidhit[w], nRecHits);
    }

    // histos for fake rate vs Number of PU interactions
    fillPlotNoFlow(h_recopu[w], PU_NumInteractions);
    if (Track_is_matched) {
      fillPlotNoFlow(h_assoc2pu[w], PU_NumInteractions);
    }

    // Fill other histos
    TrackingParticle* tpp = const_cast<TrackingParticle*>(tpr.get());
    // TrackingParticle parameters at point of closest approach to the beamline
    TrackingParticle::Vector momentumTP;
    TrackingParticle::Point vertexTP;
    
    if (parametersDefiner=="LhcParametersDefinerForTP") {
      // following reco plots are made only from tracks associated to selected signal TPs 
      if (! (Track_is_matched && tpSelector(*tpp)) ) continue;
      else {
        momentumTP = Lhc_parametersDefinerTP->momentum(event,setup,tpr) ;
        vertexTP = Lhc_parametersDefinerTP->vertex(event,setup,tpr);
      }
    }
    else if (parametersDefiner=="CosmicParametersDefinerForTP") {
      // following reco plots are made only from tracks associated to selected signal TPs 
      if (! (Track_is_matched && cosmictpSelector(tpr,&bs,event,setup)) ) continue;
      else {
        momentumTP = Cosmic_parametersDefinerTP->momentum(event,setup,tpr) ;
        vertexTP = Cosmic_parametersDefinerTP->vertex(event,setup,tpr);     
      }
    }
    
    if (associators[ww]=="trackAssociatorByChi2"){
      //association chi2
      double assocChi2 = -tp.begin()->second;//in association map is stored -chi2
      h_assochi2[www]->Fill(assocChi2);
      h_assochi2_prob[www]->Fill(TMath::Prob((assocChi2)*5,5));
    }
    else if (associators[ww]=="trackAssociatorByHits"){
      double fraction = tp.begin()->second;
      h_assocFraction[www]->Fill(fraction);
      h_assocSharedHit[www]->Fill(fraction*nRecHits);
    }
    
    h_charge[w]->Fill(track->charge());

    // Hits
    h_nhits[w]->Fill(nRecHits);
    nhits_vs_eta[w]->Fill(xetaRec, nRecHits);
    nhits_vs_phi[w]->Fill(phiRec, nRecHits);

    if (do_MUOhitsPlots) {
      nDThits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonDTHits());
      nCSChits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonCSCHits());
      nRPChits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonRPCHits()); 
      if(useGEMs_) nGEMhits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonGEMHits());
      if(useME0_) nME0hits_vs_eta[w]->Fill(xetaRec,track->hitPattern().numberOfValidMuonME0Hits());
    }

    if (do_TRKhitsPlots) {
      nTRK_LayersWithMeas_vs_eta[w]->Fill(xetaRec,track->hitPattern().trackerLayersWithMeasurement());
      nPixel_LayersWithMeas_vs_eta[w]->Fill(xetaRec,track->hitPattern().pixelLayersWithMeasurement());
      h_nlosthits[w]->Fill(track->numberOfLostHits());
      h_nmisslayers_inner[w]->Fill(track->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_INNER_HITS));
      h_nmisslayers_outer[w]->Fill(track->hitPattern().numberOfLostHits(reco::HitPattern::MISSING_OUTER_HITS));
      nlosthits_vs_eta[w]->Fill(xetaRec,track->numberOfLostHits());
    }

    // normalized chi2
    h_nchi2[w]->Fill(track->normalizedChi2());
    h_nchi2_prob[w]->Fill(TMath::Prob(track->chi2(),(int)track->ndof()));
    chi2_vs_nhits[w]->Fill(nRecHits, track->normalizedChi2());
    chi2_vs_eta[w]->Fill(xetaRec, track->normalizedChi2());
    chi2_vs_phi[w]->Fill(phiRec,track->normalizedChi2());
    
    double ptSim = sqrt(momentumTP.perp2());
    double xptSim = getPt(ptSim);
    double qoverpSim = tpr->charge() /
      sqrt(momentumTP.x()*momentumTP.x()+momentumTP.y()*momentumTP.y()+momentumTP.z()*momentumTP.z());
    double etaSim = momentumTP.eta();
    double thetaSim = momentumTP.theta();
    double phiSim = momentumTP.phi();
    double dxySim = TrackingParticleIP::dxy(vertexTP, momentumTP, bs.position());
    double dzSim = TrackingParticleIP::dz(vertexTP, momentumTP, bs.position());
    
    double etares = etaRec - etaSim;
    double ptRelRes = (ptRec - ptSim) / ptSim;    // relative residual -> resolution
    double ptPull = (ptRec - ptSim) / ptError; 
    double qoverpPull = (qoverpRec-qoverpSim) / qoverpError;
    double thetaPull = (thetaRec - thetaSim) / thetaError;
    double phiDiff = phiRec - phiSim;
    if (abs(phiDiff) > M_PI) {
      if (phiDiff >0.) phiDiff = phiDiff - 2.*M_PI;
      else phiDiff = phiDiff + 2.*M_PI;
    }
    double phiPull = phiDiff / phiError;
    double dxyPull = (dxyRec-dxySim) / dxyError;
    double dzPull = (dzRec-dzSim) / dzError;
    
    h_etaRes[w]->Fill(etares);
    etares_vs_eta[w]->Fill(xetaRec, etares);
    
    ptres_vs_eta[w]->Fill(xetaRec,ptRelRes);
    ptres_vs_pt[w]->Fill(xptSim,ptRelRes);
    ptres_vs_phi[w]->Fill(phiRec,ptRelRes);
    h_ptpull[w]->Fill(ptPull);
    ptpull_vs_eta[w]->Fill(xetaRec,ptPull);
    ptpull_vs_phi[w]->Fill(phiRec,ptPull);
    h_qoverppull[w]->Fill(qoverpPull);

    thetaCotres_vs_eta[w]->Fill(xetaRec, cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
    thetaCotres_vs_pt[w]->Fill(xptSim, cos(thetaRec)/sin(thetaRec) - cos(thetaSim)/sin(thetaSim));
    h_thetapull[w]->Fill(thetaPull);
    thetapull_vs_eta[w]->Fill(xetaRec,thetaPull);
    thetapull_vs_phi[w]->Fill(phiRec,thetaPull);

    phires_vs_eta[w]->Fill(xetaRec,phiDiff);
    phires_vs_pt[w]->Fill(xptSim,phiDiff);
    phires_vs_phi[w]->Fill(phiRec,phiDiff);
    h_phipull[w]->Fill(phiPull);
    phipull_vs_eta[w]->Fill(xetaRec,phiPull);
    phipull_vs_phi[w]->Fill(phiRec,phiPull);

    dxyres_vs_eta[w]->Fill(xetaRec,dxyRec-dxySim);
    dxyres_vs_pt[w]->Fill(xptSim,dxyRec-dxySim);
    h_dxypull[w]->Fill(dxyPull);
    dxypull_vs_eta[w]->Fill(xetaRec,dxyPull);

    dzres_vs_eta[w]->Fill(xetaRec,dzRec-dzSim);
    dzres_vs_pt[w]->Fill(xptSim,dzRec-dzSim);
    h_dzpull[w]->Fill(dzPull);  
    dzpull_vs_eta[w]->Fill(xetaRec,dzPull);

    double contrib_Qoverp = qoverpPull*qoverpPull/5;
    double contrib_dxy = dxyPull*dxyPull/5;
    double contrib_dz = dzPull*dzPull/5;
    double contrib_theta = thetaPull*thetaPull/5;
    double contrib_phi = phiPull*phiPull/5;
    double assoChi2 = contrib_Qoverp+contrib_dxy+contrib_dz+contrib_theta+contrib_phi;

    edm::LogVerbatim("MuonTrackValidator") << "normalized Chi2 (track 5-dofs matching) = "<<assoChi2 <<"\n"
                           << "\t contrib_Qoverp = " << contrib_Qoverp << "\n"
                           << "\t contrib_theta = " << contrib_theta << "\n"
                           << "\t contrib_phi = " << contrib_phi << "\n"
                           << "\t contrib_dxy = " << contrib_dxy << "\n"
                           << "\t contrib_dz = " << contrib_dz << "\n";
    
    edm::LogVerbatim("MuonTrackValidator") << "ptRec = " << ptRec << "\n"
                           << "etaRec = " << etaRec << "\n"
                           << "qoverpRec = " << qoverpRec << "\n"
                           << "thetaRec = " << thetaRec << "\n"
                           << "phiRec = " << phiRec << "\n"
                           << "dxyRec = " << dxyRec << "\n"
                           << "dzRec = " << dzRec << "\n"
                           << "" << "\n"
                           << "qoverpError = " << qoverpError << "\n"
                           << "thetaError = " << thetaError << "\n"
                           << "phiError = " << phiError << "\n"
                           << "dxyError = " << dxyError << "\n"
                           << "dzError = " << dzError << "\n"
                           << "" << "\n"
                           << "ptSim = " << ptSim << "\n"
                           << "etaSim = " << etaSim << "\n"
                           << "qoverpSim = " << qoverpSim << "\n"
                           << "thetaSim = " << thetaSim << "\n"
                           << "phiSim = " << phiSim << "\n"
                           << "dxySim = " << dxySim << "\n"
                           << "dzSim = " << dzSim << "\n";
      } // End of for(edm::View<Track>::size_type i=0; i<trackCollectionSize; ++i) {
      
      h_tracks[w]->Fill(at);
      h_fakes[w]->Fill(rT-at);
      edm::LogVerbatim("MuonTrackValidator") << "Total Simulated: " << st << "\n"
                         << "Total Associated (simToReco): " << ats << "\n"
                         << "Total Reconstructed: " << rT << "\n"
                         << "Total Associated (recoToSim): " << at << "\n"
                         << "Total Fakes: " << rT-at << "\n";
       w++;
    } // End of for (unsigned int www=0;www<label.size();www++){
  } //END of for (unsigned int ww=0;ww<associators.size();ww++){
}

void MuonTrackValidator::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &  setup  )

override

Method called at the end of the event loop.

Definition at line 25 of file MuonTrackValidator.cc.

References patPFMETCorrections_cff::algo, trackingPlots::assoc, simPFProducer_cfi::associators, DQMStore::IBooker::book1D(), DQMStore::IBooker::book2D(), DQMStore::IBooker::cd(), histoParameters_cff::cotThetaRes_nbin, histoParameters_cff::cotThetaRes_rangeMax, histoParameters_cff::cotThetaRes_rangeMin, TrackerOfflineValidation_Dqm_cff::dirName, histoParameters_cff::do_MUOhitsPlots, histoParameters_cff::do_TRKhitsPlots, histoParameters_cff::dxyRes_nbin, histoParameters_cff::dxyRes_rangeMax, histoParameters_cff::dxyRes_rangeMin, histoParameters_cff::dzRes_nbin, histoParameters_cff::dzRes_rangeMax, histoParameters_cff::dzRes_rangeMin, histoParameters_cff::etaRes_nbin, histoParameters_cff::etaRes_rangeMax, histoParameters_cff::etaRes_rangeMin, edm::InputTag::instance(), label, edm::InputTag::label(), associators_cff::MABH, histoParameters_cff::maxCSCHit, histoParameters_cff::maxDTHit, gather_cfg::maxDxy, cutBasedMuonId_MuonPOG_V0_cff::maxDz, maxEta, histoParameters_cff::maxLayers, histoParameters_cff::maxNHit, trackingParticleSelector_cfi::maxPhi, histoParameters_cff::maxPixels, MuonErrorMatrixAnalyzer_cfi::maxPt, pileupFilter_cff::maxPU, histoParameters_cff::maxRPCHit, histoParameters_cff::maxRpos, HGVHistoProducerAlgoBlock_cfi::maxZpos, histoParameters_cff::minCSCHit, histoParameters_cff::minDTHit, histoParameters_cff::minDxy, histoParameters_cff::minDz, trackingPOGFilters_cfi::minEta, dt1dClusters_cfi::minLayers, histoParameters_cff::minNHit, trackingParticleSelector_cfi::minPhi, histoParameters_cff::minPixels, lostTracks_cfi::minPt, pileupFilter_cff::minPU, histoParameters_cff::minRPCHit, histoParameters_cff::minRpos, HGVHistoProducerAlgoBlock_cfi::minZpos, histoParameters_cff::nintCSCHit, histoParameters_cff::nintDTHit, histoParameters_cff::nintDxy, histoParameters_cff::nintDz, HGVHistoProducerAlgoBlock_cfi::nintEta, histoParameters_cff::nintLayers, histoParameters_cff::nintNHit, HGVHistoProducerAlgoBlock_cfi::nintPhi, histoParameters_cff::nintPixels, HGVHistoProducerAlgoBlock_cfi::nintPt, histoParameters_cff::nintPU, histoParameters_cff::nintRPCHit, histoParameters_cff::nintRpos, HGVHistoProducerAlgoBlock_cfi::nintZpos, histoParameters_cff::phiRes_nbin, histoParameters_cff::phiRes_rangeMax, histoParameters_cff::phiRes_rangeMin, edm::InputTag::process(), histoParameters_cff::ptRes_nbin, histoParameters_cff::ptRes_rangeMax, histoParameters_cff::ptRes_rangeMin, python.rootplot.root2matplotlib::replace(), DQMStore::IBooker::setCurrentFolder(), and TrackValidationHeavyIons_cff::useLogPt.

Referenced by ~MuonTrackValidator().

                                                                                                           {

  for (unsigned int ww=0;ww<associators.size();ww++){
    for (unsigned int www=0;www<label.size();www++){

      ibooker.cd();
      InputTag algo = label[www];
      string dirName=dirName_;
      if (algo.process()!="")
    dirName+=algo.process()+"_";
      if(algo.label()!="")
    dirName+=algo.label();
      if(algo.instance()!="")
        dirName+=("_"+algo.instance());
      if (dirName.find("Tracks")<dirName.length()){
    dirName.replace(dirName.find("Tracks"),6,"Trks");
      }
      if (dirName.find("UpdatedAtVtx")<dirName.length()){
        dirName.replace(dirName.find("UpdatedAtVtx"),12,"UpdAtVtx");
      }
      string assoc= associators[ww];
      if (assoc.find("tpToTkmuTrackAssociation")<assoc.length()){
        dirName+="_TkAsso";
      }
      std::replace(dirName.begin(), dirName.end(), ':', '_');
      ibooker.setCurrentFolder(dirName);

      h_tracks.push_back( ibooker.book1D("Ntracks","Number of reconstructed tracks",100,-0.5,99.5) );
      h_fakes.push_back( ibooker.book1D("Nfakes","Number of fake reco tracks",20,-0.5,19.5) );
      h_charge.push_back( ibooker.book1D("Ncharge","track charge",3,-1.5,1.5) );

      h_recoeta.push_back( ibooker.book1D("num_reco_eta","N of reco track vs eta",nintEta,minEta,maxEta) );
      h_assoceta.push_back( ibooker.book1D("num_assoSimToReco_eta","N of associated tracks (simToReco) vs eta",nintEta,minEta,maxEta) );
      h_assoc2eta.push_back( ibooker.book1D("num_assoRecoToSim_eta","N of associated (recoToSim) tracks vs eta",nintEta,minEta,maxEta) );
      h_simuleta.push_back( ibooker.book1D("num_simul_eta","N of simulated tracks vs eta",nintEta,minEta,maxEta) );
      h_misideta.push_back( ibooker.book1D("num_chargemisid_eta","N of associated (simToReco) tracks with charge misID vs eta",nintEta,minEta,maxEta) );

      h_recopT.push_back( ibooker.book1D("num_reco_pT","N of reco track vs pT",nintPt,minPt,maxPt) );
      h_assocpT.push_back( ibooker.book1D("num_assoSimToReco_pT","N of associated tracks (simToReco) vs pT",nintPt,minPt,maxPt) );
      h_assoc2pT.push_back( ibooker.book1D("num_assoRecoToSim_pT","N of associated (recoToSim) tracks vs pT",nintPt,minPt,maxPt) );
      h_simulpT.push_back( ibooker.book1D("num_simul_pT","N of simulated tracks vs pT",nintPt,minPt,maxPt) );
      h_misidpT.push_back( ibooker.book1D("num_chargemisid_pT","N of associated (simToReco) tracks with charge misID vs pT",nintPt,minPt,maxPt) );      

      h_recophi.push_back( ibooker.book1D("num_reco_phi","N of reco track vs phi",nintPhi,minPhi,maxPhi) );
      h_assocphi.push_back( ibooker.book1D("num_assoSimToReco_phi","N of associated tracks (simToReco) vs phi",nintPhi,minPhi,maxPhi) );
      h_assoc2phi.push_back( ibooker.book1D("num_assoRecoToSim_phi","N of associated (recoToSim) tracks vs phi",nintPhi,minPhi,maxPhi) );
      h_simulphi.push_back( ibooker.book1D("num_simul_phi","N of simulated tracks vs phi",nintPhi,minPhi,maxPhi) );
      h_misidphi.push_back( ibooker.book1D("num_chargemisid_phi","N of associated (simToReco) tracks with charge misID vs phi",nintPhi,minPhi,maxPhi) );

      h_recohit.push_back( ibooker.book1D("num_reco_hit","N of reco tracks vs N SimHits",nintNHit,minNHit,maxNHit) );
      h_assochit.push_back( ibooker.book1D("num_assoSimToReco_hit","N of associated tracks (simToReco) vs N SimHits",nintNHit,minNHit,maxNHit) );
      h_assoc2hit.push_back( ibooker.book1D("num_assoRecoToSim_hit","N of associated (recoToSim) tracks vs N Rechits",nintNHit,minNHit,maxNHit) );
      h_simulhit.push_back( ibooker.book1D("num_simul_hit","N of simulated tracks vs N SimHits",nintNHit,minNHit,maxNHit) );
      h_misidhit.push_back( ibooker.book1D("num_chargemisid_hit","N of associated (recoToSim) tracks with charge misID vs N RecHits",nintNHit,minNHit,maxNHit) );

      h_recodxy.push_back( ibooker.book1D("num_reco_dxy","N of reco track vs dxy",nintDxy,minDxy,maxDxy) );
      h_assocdxy.push_back( ibooker.book1D("num_assoSimToReco_dxy","N of associated tracks (simToReco) vs dxy",nintDxy,minDxy,maxDxy) );
      h_assoc2dxy.push_back( ibooker.book1D("num_assoRecoToSim_dxy","N of associated (recoToSim) tracks vs dxy",nintDxy,minDxy,maxDxy) );
      h_simuldxy.push_back( ibooker.book1D("num_simul_dxy","N of simulated tracks vs dxy",nintDxy,minDxy,maxDxy) );
      h_misiddxy.push_back( ibooker.book1D("num_chargemisid_dxy","N of associated (simToReco) tracks with charge misID vs dxy",nintDxy,minDxy,maxDxy) );           
      h_recodz.push_back( ibooker.book1D("num_reco_dz","N of reco track vs dz",nintDz,minDz,maxDz) );
      h_assocdz.push_back( ibooker.book1D("num_assoSimToReco_dz","N of associated tracks (simToReco) vs dz",nintDz,minDz,maxDz) );
      h_assoc2dz.push_back( ibooker.book1D("num_assoRecoToSim_dz","N of associated (recoToSim) tracks vs dz",nintDz,minDz,maxDz) );
      h_simuldz.push_back( ibooker.book1D("num_simul_dz","N of simulated tracks vs dz",nintDz,minDz,maxDz) );
      h_misiddz.push_back( ibooker.book1D("num_chargemisid_dz","N of associated (simToReco) tracks with charge misID vs dz",nintDz,minDz,maxDz) );

      h_assocRpos.push_back( ibooker.book1D("num_assoSimToReco_Rpos","N of associated tracks (simToReco) vs Radius",nintRpos,minRpos,maxRpos) );
      h_simulRpos.push_back( ibooker.book1D("num_simul_Rpos","N of simulated tracks vs Radius",nintRpos,minRpos,maxRpos) );

      h_assocZpos.push_back( ibooker.book1D("num_assoSimToReco_Zpos","N of associated tracks (simToReco) vs Z",nintZpos,minZpos,maxZpos) );
      h_simulZpos.push_back( ibooker.book1D("num_simul_Zpos","N of simulated tracks vs Z",nintZpos,minZpos,maxZpos) );

      h_recopu.push_back( ibooker.book1D("num_reco_pu","N of reco track vs pu",nintPU,minPU,maxPU) );
      h_assocpu.push_back( ibooker.book1D("num_assoSimToReco_pu","N of associated tracks (simToReco) vs pu",nintPU,minPU,maxPU) );
      h_assoc2pu.push_back( ibooker.book1D("num_assoRecoToSim_pu","N of associated (recoToSim) tracks vs pu",nintPU,minPU,maxPU) );
      h_simulpu.push_back( ibooker.book1D("num_simul_pu","N of simulated tracks vs pu",nintPU,minPU,maxPU) );
      h_misidpu.push_back( ibooker.book1D("num_chargemisid_pu","N of associated (simToReco) charge misIDed tracks vs pu",nintPU,minPU,maxPU) );

      h_nchi2.push_back( ibooker.book1D("chi2","Track normalized #chi^{2}", 80, 0., 20.) );
      h_nchi2_prob.push_back( ibooker.book1D("chi2prob", "Probability of track normalized #chi^{2}",100,0.,1.) );

      chi2_vs_nhits.push_back( ibooker.book2D("chi2_vs_nhits","#chi^{2} vs nhits",nintNHit,minNHit,maxNHit,20,0.,10.) );
      chi2_vs_eta.push_back( ibooker.book2D("chi2_vs_eta","chi2_vs_eta",nintEta,minEta,maxEta, 40, 0., 20. ));
      chi2_vs_phi.push_back( ibooker.book2D("chi2_vs_phi","#chi^{2} vs #phi",nintPhi,minPhi,maxPhi, 40, 0., 20. ) );

      h_nhits.push_back( ibooker.book1D("nhits", "Number of hits per track", nintNHit,minNHit,maxNHit ) );
      nhits_vs_eta.push_back( ibooker.book2D("nhits_vs_eta","Number of Hits vs eta",nintEta,minEta,maxEta,nintNHit,minNHit,maxNHit) );
      nhits_vs_phi.push_back( ibooker.book2D("nhits_vs_phi","#hits vs #phi",nintPhi,minPhi,maxPhi,nintNHit,minNHit,maxNHit) );
      
      if (do_MUOhitsPlots) {
    nDThits_vs_eta.push_back( ibooker.book2D("nDThits_vs_eta","Number of DT hits vs eta",nintEta,minEta,maxEta,nintDTHit,minDTHit,maxDTHit) );
    nCSChits_vs_eta.push_back( ibooker.book2D("nCSChits_vs_eta","Number of CSC hits vs eta",nintEta,minEta,maxEta,nintCSCHit,minCSCHit,maxCSCHit) );
    nRPChits_vs_eta.push_back( ibooker.book2D("nRPChits_vs_eta","Number of RPC hits vs eta",nintEta,minEta,maxEta,nintRPCHit,minRPCHit,maxRPCHit) );
    if (useGEMs_)
      nGEMhits_vs_eta.push_back( ibooker.book2D("nGEMhits_vs_eta","Number of GEM hits vs eta",nintEta,minEta,maxEta,nintNHit,minNHit,maxNHit) );
    if(useME0_)
      nME0hits_vs_eta.push_back( ibooker.book2D("nME0hits_vs_eta","Number of ME0 hits vs eta",nintEta,minEta,maxEta,nintNHit,minNHit,maxNHit) );

      }

      if (do_TRKhitsPlots) {
    nTRK_LayersWithMeas_vs_eta.push_back(ibooker.book2D("nTRK_LayersWithMeas_vs_eta","# TRK Layers with measurement vs eta",
                                nintEta,minEta,maxEta,nintLayers,minLayers,maxLayers) );
    nPixel_LayersWithMeas_vs_eta.push_back(ibooker.book2D("nPixel_LayersWithMeas_vs_eta","Number of Pixel Layers with measurement vs eta",
                                  nintEta,minEta,maxEta,nintPixels,minPixels,maxPixels) );
    h_nmisslayers_inner.push_back( ibooker.book1D("nTRK_misslayers_inner", "Number of missing inner TRK layers", nintLayers,minLayers,maxLayers ) );
    h_nmisslayers_outer.push_back( ibooker.book1D("nTRK_misslayers_outer", "Number of missing outer TRK layers", nintLayers,minLayers,maxLayers ) );
    h_nlosthits.push_back( ibooker.book1D("nlosthits", "Number of lost hits per track", 6,-0.5,5.5 ) );
    nlosthits_vs_eta.push_back( ibooker.book2D("nlosthits_vs_eta","Number of lost hits per track vs eta",nintEta,minEta,maxEta,6,-0.5,5.5) );
      }

      ptres_vs_eta.push_back(ibooker.book2D("ptres_vs_eta","p_{T} Relative Residual vs #eta",
                        nintEta,minEta,maxEta, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
      ptres_vs_phi.push_back( ibooker.book2D("ptres_vs_phi","p_{T} Relative Residual vs #phi",
                         nintPhi,minPhi,maxPhi, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
      ptres_vs_pt.push_back(ibooker.book2D("ptres_vs_pt","p_{T} Relative Residual vs p_{T}",
                       nintPt,minPt,maxPt, ptRes_nbin, ptRes_rangeMin, ptRes_rangeMax));
      h_ptpull.push_back( ibooker.book1D("ptpull", "p_{T} Pull", 100, -10., 10.) );
      ptpull_vs_eta.push_back(ibooker.book2D("ptpull_vs_eta","p_{T} Pull vs #eta",nintEta,minEta,maxEta,100,-10.,10.));
      ptpull_vs_phi.push_back(ibooker.book2D("ptpull_vs_phi","p_{T} Pull vs #phi",nintPhi,minPhi,maxPhi,100,-10.,10.));
      h_qoverppull.push_back( ibooker.book1D("qoverppull","q/p Pull", 100, -10., 10.) );

      h_etaRes.push_back( ibooker.book1D("etaRes", "#eta residual", etaRes_nbin, etaRes_rangeMin, etaRes_rangeMax) );
      etares_vs_eta.push_back( ibooker.book2D("etares_vs_eta","#eta Residual vs #eta",
                          nintEta,minEta,maxEta, etaRes_nbin, etaRes_rangeMin, etaRes_rangeMax) );

      thetaCotres_vs_eta.push_back(ibooker.book2D("thetaCotres_vs_eta","cot(#theta) Residual vs #eta",
                          nintEta,minEta,maxEta,cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
      thetaCotres_vs_pt.push_back(ibooker.book2D("thetaCotres_vs_pt","cot(#theta) Residual vs p_{T}",
                         nintPt,minPt,maxPt, cotThetaRes_nbin, cotThetaRes_rangeMin, cotThetaRes_rangeMax));
      h_thetapull.push_back( ibooker.book1D("thetapull","#theta Pull",100,-10.,10.) );
      thetapull_vs_eta.push_back(ibooker.book2D("thetapull_vs_eta","#theta Pull vs #eta",nintEta,minEta,maxEta,100,-10,10));
      thetapull_vs_phi.push_back(ibooker.book2D("thetapull_vs_phi","#theta Pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));

      phires_vs_eta.push_back(ibooker.book2D("phires_vs_eta","#phi Residual vs #eta",
                         nintEta,minEta,maxEta, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
      phires_vs_pt.push_back(ibooker.book2D("phires_vs_pt","#phi Residual vs p_{T}",
                        nintPt,minPt,maxPt, phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
      phires_vs_phi.push_back(ibooker.book2D("phires_vs_phi","#phi Residual vs #phi",
                         nintPhi,minPhi,maxPhi,phiRes_nbin, phiRes_rangeMin, phiRes_rangeMax));
      h_phipull.push_back( ibooker.book1D("phipull","#phi Pull",100,-10.,10.) );
      phipull_vs_eta.push_back(ibooker.book2D("phipull_vs_eta","#phi Pull vs #eta",nintEta,minEta,maxEta,100,-10,10));
      phipull_vs_phi.push_back(ibooker.book2D("phipull_vs_phi","#phi Pull vs #phi",nintPhi,minPhi,maxPhi,100,-10,10));

      dxyres_vs_eta.push_back(ibooker.book2D("dxyres_vs_eta","dxy Residual vs #eta",
                         nintEta,minEta,maxEta,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
      dxyres_vs_pt.push_back( ibooker.book2D("dxyres_vs_pt","dxy Residual vs p_{T}",
                         nintPt,minPt,maxPt,dxyRes_nbin, dxyRes_rangeMin, dxyRes_rangeMax));
      h_dxypull.push_back( ibooker.book1D("dxypull","dxy Pull",100,-10.,10.) );
      dxypull_vs_eta.push_back(ibooker.book2D("dxypull_vs_eta","dxy Pull vs #eta",nintEta,minEta,maxEta,100,-10,10));

      dzres_vs_eta.push_back(ibooker.book2D("dzres_vs_eta","dz Residual vs #eta",
                        nintEta,minEta,maxEta,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
      dzres_vs_pt.push_back(ibooker.book2D("dzres_vs_pt","dz Residual vs p_{T}",
                           nintPt,minPt,maxPt,dzRes_nbin, dzRes_rangeMin, dzRes_rangeMax));
      h_dzpull.push_back( ibooker.book1D("dzpull","dz Pull",100,-10.,10.) );
      dzpull_vs_eta.push_back(ibooker.book2D("dzpull_vs_eta","dz Pull vs #eta",nintEta,minEta,maxEta,100,-10,10));

      nRecHits_vs_nSimHits.push_back( ibooker.book2D("nRecHits_vs_nSimHits","nRecHits vs nSimHits",nintNHit,minNHit,maxNHit, nintNHit,minNHit,maxNHit )); 
     
      if (MABH) {
    h_PurityVsQuality.push_back
      (ibooker.book2D("PurityVsQuality","Purity vs Quality (MABH)",20,0.01,1.01,20,0.01,1.01) );
      }

      if (associators[ww]=="trackAssociatorByChi2"){
    h_assochi2.push_back( ibooker.book1D("assocChi2","track association #chi^{2}",1000,0.,100.) );
    h_assochi2_prob.push_back(ibooker.book1D("assocChi2_prob","probability of association #chi^{2}",100,0.,1.));
      } else if (associators[ww]=="trackAssociatorByHits"){
    h_assocFraction.push_back( ibooker.book1D("assocFraction","fraction of shared hits",22,0.,1.1) );
    h_assocSharedHit.push_back(ibooker.book1D("assocSharedHit","number of shared hits",41,-0.5,40.5));
      }
      
      if (useLogPt) {
    BinLogX(h_simulpT.back()->getTH1F());
    BinLogX(h_assocpT.back()->getTH1F());
    BinLogX(h_recopT.back()->getTH1F());
    BinLogX(h_assoc2pT.back()->getTH1F());
    BinLogX(h_misidpT.back()->getTH1F());

    BinLogX(phires_vs_pt.back()->getTH2F());
    BinLogX(thetaCotres_vs_pt.back()->getTH2F());
    BinLogX(dxyres_vs_pt.back()->getTH2F());
    BinLogX(dzres_vs_pt.back()->getTH2F());
    BinLogX(ptres_vs_pt.back()->getTH2F());
      }

    }   //for (unsigned int www=0;www<label.size();www++)
  }     //for (unsigned int ww=0;ww<associators.size();ww++)

}

Member Data Documentation

edm::EDGetTokenT<SimHitTPAssociationProducer::SimHitTPAssociationList> MuonTrackValidator::_simHitTpMapTag

private

Definition at line 152 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

edm::InputTag MuonTrackValidator::associatormap

private

Definition at line 149 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::BiDirectional_RecoToSim_association

private

Definition at line 164 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

CosmicTrackingParticleSelector MuonTrackValidator::cosmictpSelector

private

Definition at line 161 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

std::string MuonTrackValidator::dirName_

private

Definition at line 148 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::MABH

private

Definition at line 166 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

edm::EDGetTokenT<reco::RecoToSimCollection> MuonTrackValidator::recoToSimCollection_Token

private

Definition at line 151 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

edm::EDGetTokenT<reco::SimToRecoCollection> MuonTrackValidator::simToRecoCollection_Token

private

Definition at line 150 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

TrackingParticleSelector MuonTrackValidator::tpSelector

private

Definition at line 160 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::UseAssociators

private

Definition at line 154 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::useGEMs_

private

Definition at line 155 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().

bool MuonTrackValidator::useME0_

private

Definition at line 156 of file MuonTrackValidator.h.

Referenced by MuonTrackValidator().