MuonMiniAOD Class Reference

#include <MuonMiniAOD.h>

Inheritance diagram for MuonMiniAOD:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	Inizialize parameters for histo binning. More...
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	MuonMiniAOD (const edm::ParameterSet &)
	Constructor. More...
	~MuonMiniAOD () 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 Member Functions
bool	PassesCut_A (edm::View< pat::Muon >::const_iterator, reco::Vertex, TString)
bool	PassesCut_B (edm::View< pat::Muon >::const_iterator, reco::Vertex, TString)

Private Attributes
bool	doPVCheck_
std::string	metname
edm::ParameterSet	parameters
edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotLabel_
edm::EDGetTokenT< edm::View< pat::Muon > >	theMuonCollectionLabel_
edm::EDGetTokenT< reco::VertexCollection >	theVertexLabel_
std::vector< MonitorElement * >	workingPoints

Detailed Description

Definition at line 29 of file MuonMiniAOD.h.

Constructor & Destructor Documentation

MuonMiniAOD::MuonMiniAOD

(

const edm::ParameterSet &

pSet

)

Constructor.

Definition at line 11 of file MuonMiniAOD.cc.

                                                    {
  parameters = pSet;

  // the services:
  
  theMuonCollectionLabel_ = consumes<edm::View<pat::Muon> >  (parameters.getParameter<edm::InputTag>("MuonCollection"));
  theVertexLabel_          = consumes<reco::VertexCollection>(parameters.getParameter<edm::InputTag>("VertexLabel"));
  theBeamSpotLabel_        = mayConsume<reco::BeamSpot>      (parameters.getParameter<edm::InputTag>("BeamSpotLabel"));

}

MuonMiniAOD::~MuonMiniAOD ( )

override

Destructor.

Definition at line 23 of file MuonMiniAOD.cc.

{ }

Member Function Documentation

void MuonMiniAOD::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Inizialize parameters for histo binning.

Definition at line 87 of file MuonMiniAOD.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), relativeConstraints::error, edm::Event::getByToken(), edm::HandleBase::isValid(), LogTrace, metname, extraflags_cff::muons, reco::BeamSpot::position(), position, reco::BeamSpot::sigmaZ(), and mvaElectronID_Fall17_iso_V2_cff::workingPoints.

                                                                            {
  LogTrace(metname)<<"[MuonMiniAOD] Analyze the mu";
  
 
  // Take the muon container
  edm::Handle<edm::View<pat::Muon> > muons; 
  iEvent.getByToken(theMuonCollectionLabel_,muons);
 
  //Vertex information
  edm::Handle<reco::VertexCollection> vertex;
  iEvent.getByToken(theVertexLabel_, vertex);  


  if(!muons.isValid()) return;

  reco::Vertex::Point posVtx;
  reco::Vertex::Error errVtx;
  unsigned int theIndexOfThePrimaryVertex = 999.;
  if (!vertex.isValid()) {
    LogTrace(metname) << "[EfficiencyAnalyzer] Could not find vertex collection" << std::endl;
    for (unsigned int ind=0; ind<vertex->size(); ++ind) {
      if ( (*vertex)[ind].isValid() && !((*vertex)[ind].isFake()) ) {
    theIndexOfThePrimaryVertex = ind;
    break;
      }
    }
  }

  if (theIndexOfThePrimaryVertex<100) {
    posVtx = ((*vertex)[theIndexOfThePrimaryVertex]).position();
    errVtx = ((*vertex)[theIndexOfThePrimaryVertex]).error();
  }   
  else {
    LogInfo("RecoMuonValidator") << "reco::PrimaryVertex not found, use BeamSpot position instead\n";
    
    edm::Handle<reco::BeamSpot> recoBeamSpotHandle;
    iEvent.getByToken(theBeamSpotLabel_,recoBeamSpotHandle);
    reco::BeamSpot bs = *recoBeamSpotHandle;
    
    posVtx = bs.position();
    errVtx(0,0) = bs.BeamWidthX();
    errVtx(1,1) = bs.BeamWidthY();
    errVtx(2,2) = bs.sigmaZ();
  }
    


  const reco::Vertex thePrimaryVertex(posVtx,errVtx);

  for (edm::View<pat::Muon>::const_iterator muon1 = muons->begin(); muon1 != muons->end(); ++muon1){

    for (std::vector<MonitorElement*>::iterator monitor = workingPoints.begin();
     monitor != workingPoints.end(); ++monitor){
      int Pass_A = 0;
      int Pass_B = 0;
      if (PassesCut_A(muon1, thePrimaryVertex,(*monitor)->getName()))
    Pass_A = 1;
      else
    Pass_A = 2;
      if (PassesCut_B(muon1, thePrimaryVertex,(*monitor)->getName()))
    Pass_B = 1;
      else
    Pass_B = 2;

      (*monitor)->Fill(Pass_A,Pass_B);
    }

  }
}

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

override

Definition at line 24 of file MuonMiniAOD.cc.

References DQMStore::IBooker::book2D(), DQMStore::IBooker::cd(), DQMStore::IBooker::setCurrentFolder(), and mvaElectronID_Fall17_iso_V2_cff::workingPoints.

                                              {
    
  ibooker.cd();
  ibooker.setCurrentFolder("Muons_miniAOD/MuonMiniAOD");
  

  workingPoints.push_back(ibooker.book2D("tightMuons" ,"Tight Muons"  ,2,1,3,2,1,3));
  workingPoints.push_back(ibooker.book2D("mediumMuons","Medium Muons" ,2,1,3,2,1,3));
  workingPoints.push_back(ibooker.book2D("looseMuons" ,"Loose Muons"  ,2,1,3,2,1,3));
  workingPoints.push_back(ibooker.book2D("highPtMuons","High Pt Muons",2,1,3,2,1,3));
  workingPoints.push_back(ibooker.book2D("softMuons"  ,"Soft Muons"   ,2,1,3,2,1,3));

  for (std::vector<MonitorElement*>::iterator monitor = workingPoints.begin();
       monitor != workingPoints.end(); ++monitor){
    (*monitor)-> setBinLabel(1,"Pass",1);
    (*monitor) -> setBinLabel(2,"No Pass",1);
    (*monitor) -> setBinLabel(1,"Pass",2);
    (*monitor) -> setBinLabel(2,"No Pass",2);
  }

}

bool MuonMiniAOD::PassesCut_A ( edm::View< pat::Muon >::const_iterator  muon1, reco::Vertex  thePrimaryVertex, TString  WorkingPoint  )

private

Definition at line 48 of file MuonMiniAOD.cc.

References muon::isHighPtMuon(), muon::isLooseMuon(), muon::isMediumMuon(), muon::isSoftMuon(), and muon::isTightMuon().

                                                                                                                    {

  if (WorkingPoint == "tightMuons")
    return muon::isTightMuon(*muon1,thePrimaryVertex);
  else if (WorkingPoint == "mediumMuons")
    return muon::isMediumMuon(*muon1);
  else if (WorkingPoint == "looseMuons")
    return muon::isLooseMuon(*muon1);
  else if (WorkingPoint == "highPtMuons")
    return muon::isHighPtMuon(*muon1,thePrimaryVertex);
  else if (WorkingPoint == "softMuons")
    return muon::isSoftMuon(*muon1,thePrimaryVertex);
  else{
    LogInfo("RecoMuonValidator") << "[MuonMiniAOD]: MuonMiniAOD. Unknown WP, returning false.\n";
    return false;
  }
  
}

bool MuonMiniAOD::PassesCut_B ( edm::View< pat::Muon >::const_iterator  muon1, reco::Vertex  thePrimaryVertex, TString  WorkingPoint  )

private

Definition at line 67 of file MuonMiniAOD.cc.

References muon::isHighPtMuon(), muon::isLooseMuon(), muon::isMediumMuon(), muon::isSoftMuon(), and muon::isTightMuon().

                                                                                                                    {
  
  if (WorkingPoint == "tightMuons")
    return muon1 -> isTightMuon(thePrimaryVertex);
  else if (WorkingPoint == "mediumMuons")
    return muon1 -> isMediumMuon();
  else if (WorkingPoint == "looseMuons")
    return muon1 -> isLooseMuon();
  else if (WorkingPoint == "highPtMuons")
    return muon1 -> isHighPtMuon(thePrimaryVertex);
  else if (WorkingPoint == "softMuons")
    return muon1 -> isSoftMuon(thePrimaryVertex);
  else{
    LogInfo("RecoMuonValidator") << "[MuonMiniAOD]: MuonMiniAOD. Unknown WP, returning false.\n";
    return false;
  }
   
}

Member Data Documentation

bool MuonMiniAOD::doPVCheck_

private

Definition at line 60 of file MuonMiniAOD.h.

std::string MuonMiniAOD::metname

private

Definition at line 50 of file MuonMiniAOD.h.

edm::ParameterSet MuonMiniAOD::parameters

private

Definition at line 46 of file MuonMiniAOD.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

edm::EDGetTokenT<reco::BeamSpot> MuonMiniAOD::theBeamSpotLabel_

private

Definition at line 62 of file MuonMiniAOD.h.

edm::EDGetTokenT< edm::View<pat::Muon> > MuonMiniAOD::theMuonCollectionLabel_

private

Definition at line 48 of file MuonMiniAOD.h.

edm::EDGetTokenT<reco::VertexCollection> MuonMiniAOD::theVertexLabel_

private

Definition at line 61 of file MuonMiniAOD.h.

std::vector<MonitorElement*> MuonMiniAOD::workingPoints

private

Definition at line 53 of file MuonMiniAOD.h.