RecoSelector Class Reference

#include <RecoSelector.h>

Public Member Functions
std::string	GetName ()
void	handleObjects (const edm::Event &)
bool	isSelected (const edm::Event &)
	RecoSelector (const edm::ParameterSet &userCut_params, edm::ConsumesCollector &&iC)
virtual	~RecoSelector ()

Private Attributes
edm::EDGetTokenT < reco::CaloMETCollection >	m_calometSrc
edm::EDGetTokenT < reco::GsfElectronCollection >	m_electronSrc
edm::EDGetTokenT < reco::CaloJetCollection >	m_jetsSrc
edm::EDGetTokenT < reco::MuonCollection >	m_muonSrc
edm::EDGetTokenT < reco::PhotonCollection >	m_photon_token_
std::string	m_photonProducerSrc
std::string	m_photonSrc
std::string	name
double	reco_metMin
double	reco_ptElecMin
double	reco_ptJet1Min
double	reco_ptJet2Min
double	reco_ptMuonMin
double	reco_ptPhotMin
const reco::CaloJetCollection *	theCaloJetCollection
const reco::CaloMETCollection *	theCaloMETCollection
const reco::GsfElectronCollection *	theElectronCollection
const reco::MuonCollection *	theMuonCollection
const reco::PhotonCollection *	thePhotonCollection

Detailed Description

Definition at line 34 of file RecoSelector.h.

Constructor & Destructor Documentation

RecoSelector::RecoSelector	(	const edm::ParameterSet &	userCut_params,
		edm::ConsumesCollector &&	iC
	)

Definition at line 20 of file RecoSelector.cc.

References edm::ParameterSet::getParameter(), HLT_FULL_cff::InputTag, and mergeVDriftHistosByStation::name.

{

  name                   = userCut_params.getParameter<string>("name");
  m_electronSrc      = iC.consumes<GsfElectronCollection>(userCut_params.getParameter<string>("electrons"));
  m_muonSrc          = iC.consumes<MuonCollection>(userCut_params.getParameter<string>("muons"));
  m_jetsSrc          = iC.consumes<CaloJetCollection>(userCut_params.getParameter<string>("jets"));
  m_photonSrc        = userCut_params.getParameter<string>("photons");
  m_photonProducerSrc    = userCut_params.getParameter<string>("photonProducer");
  m_photon_token_    = iC.consumes<PhotonCollection>(edm::InputTag(m_photonProducerSrc, m_photonSrc));
  m_calometSrc           = iC.consumes<CaloMETCollection>(userCut_params.getParameter<string>("calomet"));

  reco_metMin    = userCut_params.getParameter<double>("reco_metMin") ;
  reco_ptJet1Min = userCut_params.getParameter<double>("reco_ptJet1Min");
  reco_ptJet2Min = userCut_params.getParameter<double>("reco_ptJet2Min");
  reco_ptElecMin = userCut_params.getParameter<double>("reco_ptElecMin");
  reco_ptMuonMin = userCut_params.getParameter<double>("reco_ptMuonMin");
  reco_ptPhotMin = userCut_params.getParameter<double>("reco_ptPhotMin");

  edm::LogInfo("RecoSelectorParameters") << endl;
  edm::LogInfo("RecoSelectorParameters") << "UserAnalysis parameters for " << name << " selection:" << endl;
  edm::LogInfo("RecoSelectorParameters") << " reco_metMin    = " << reco_metMin << endl;
  edm::LogInfo("RecoSelectorParameters") << " reco_ptJet1Min = " << reco_ptJet1Min  << endl;
  edm::LogInfo("RecoSelectorParameters") << " reco_ptJet2Min = " << reco_ptJet2Min  << endl;
  edm::LogInfo("RecoSelectorParameters") << " reco_ptElecMin = " << reco_ptElecMin  << endl;
  edm::LogInfo("RecoSelectorParameters") << " reco_ptMuonMin = " << reco_ptMuonMin  << endl;
  edm::LogInfo("RecoSelectorParameters") << " reco_ptPhotMin = " << reco_ptPhotMin  << endl;

}

virtual RecoSelector::~RecoSelector ( )

inlinevirtual

Definition at line 41 of file RecoSelector.h.

{};

Member Function Documentation

string RecoSelector::GetName

(

)

Definition at line 50 of file RecoSelector.cc.

References mergeVDriftHistosByStation::name.

{return name;}

void RecoSelector::handleObjects

(

const edm::Event &

iEvent

)

Definition at line 122 of file RecoSelector.cc.

References edm::Event::getByToken(), and edm::Handle< T >::product().

{

  //Get the electrons
  Handle<GsfElectronCollection> theElectronCollectionHandle; 
  iEvent.getByToken(m_electronSrc, theElectronCollectionHandle);
  theElectronCollection = theElectronCollectionHandle.product();

  //Get the Muons
  Handle<MuonCollection> theMuonCollectionHandle; 
  iEvent.getByToken(m_muonSrc, theMuonCollectionHandle);
  theMuonCollection = theMuonCollectionHandle.product();

  //Get the Photons
  Handle<PhotonCollection> thePhotonCollectionHandle; 
  iEvent.getByToken(m_photon_token_, thePhotonCollectionHandle);
  thePhotonCollection = thePhotonCollectionHandle.product();

  //Get the CaloJets
  Handle<CaloJetCollection> theCaloJetCollectionHandle;
  iEvent.getByToken(m_jetsSrc, theCaloJetCollectionHandle);
  theCaloJetCollection = theCaloJetCollectionHandle.product();

  //Get the CaloMET
  Handle<CaloMETCollection> theCaloMETCollectionHandle;
  iEvent.getByToken(m_calometSrc, theCaloMETCollectionHandle);
  theCaloMETCollection = theCaloMETCollectionHandle.product();

}

bool RecoSelector::isSelected

(

const edm::Event &

iEvent

)

Definition at line 52 of file RecoSelector.cc.

References reco::GsfElectron::ecalDrivenSeed(), metsig::electron, i, LogDebug, metsig::muon, reco::LeafCandidate::p(), reco::LeafCandidate::pt(), and reco::GsfElectron::superCluster().

{


  this->handleObjects(iEvent);


  bool TotalCutPassed = false;

  bool ElectronCutPassed = false;
  for(unsigned int i=0; i<theElectronCollection->size(); i++) {
    GsfElectron electron = (*theElectronCollection)[i];
    if (electron.ecalDrivenSeed()) {
      float elenergy = electron.superCluster()->energy();
      float elpt = electron.pt() * elenergy / electron.p();
      if(elpt>reco_ptElecMin)  ElectronCutPassed = true;
      LogDebug("RecoSelectorCuts") << "elpt = " << elpt << endl;
    }
    else {
      //      float elenergy = 0;
      //      float elpt = 0;
      ElectronCutPassed = false;
    }
  }

  bool MuonCutPassed = false;
  for(unsigned int i=0; i<theMuonCollection->size(); i++) {
    Muon muon = (*theMuonCollection)[i];
    float muonpt = muon.pt();
    if(muonpt>reco_ptMuonMin) MuonCutPassed = true;
    LogDebug("RecoSelectorCuts") << "muonpt = " << muonpt << endl;
  }
  
  bool PhotonCutPassed = false;
  for(unsigned int i=0; i<thePhotonCollection->size(); i++) {
    Photon photon = (*thePhotonCollection)[i];
    float photonpt = photon.pt();
    if(photonpt>reco_ptPhotMin) PhotonCutPassed = true;
    LogDebug("RecoSelectorCuts") << "photonpt = " << photonpt << endl;
  }


  bool JetCutPassed = false;  
  if(theCaloJetCollection->size()>1) {
    if((*theCaloJetCollection)[0].pt() > reco_ptJet1Min &&
       (*theCaloJetCollection)[1].pt() > reco_ptJet2Min ) JetCutPassed = true;
  }

  bool MetCutPassed = false;  
  const CaloMET theCaloMET = theCaloMETCollection->front();
  float calomet = theCaloMET.pt();
  LogDebug("RecoSelectorCuts") << "calomet = " << calomet << endl;
  if(calomet > reco_metMin) MetCutPassed = true;

  edm::LogInfo("RecoSelectorPassed") << "ElectronCutPassed = " << (int)ElectronCutPassed << endl;
  edm::LogInfo("RecoSelectorPassed") << "MuonCutPassed     = " << (int)MuonCutPassed << endl;
  edm::LogInfo("RecoSelectorPassed") << "PhotonCutPassed   = " << (int)PhotonCutPassed << endl;
  edm::LogInfo("RecoSelectorPassed") << "JetCutPassed      = " << (int)JetCutPassed << endl;
  edm::LogInfo("RecoSelectorPassed") << "MetCutPassed      = " << (int)MetCutPassed << endl;


  if(
     (ElectronCutPassed || MuonCutPassed) &&
     PhotonCutPassed                      &&
     JetCutPassed                         &&
     MetCutPassed      )   TotalCutPassed = true;

  return TotalCutPassed;
}

Member Data Documentation

edm::EDGetTokenT<reco::CaloMETCollection> RecoSelector::m_calometSrc

private

Definition at line 58 of file RecoSelector.h.

edm::EDGetTokenT<reco::GsfElectronCollection> RecoSelector::m_electronSrc

private

Definition at line 52 of file RecoSelector.h.

edm::EDGetTokenT<reco::CaloJetCollection> RecoSelector::m_jetsSrc

private

Definition at line 54 of file RecoSelector.h.

edm::EDGetTokenT<reco::MuonCollection> RecoSelector::m_muonSrc

private

Definition at line 53 of file RecoSelector.h.

edm::EDGetTokenT<reco::PhotonCollection> RecoSelector::m_photon_token_

private

Definition at line 57 of file RecoSelector.h.

std::string RecoSelector::m_photonProducerSrc

private

Definition at line 56 of file RecoSelector.h.

std::string RecoSelector::m_photonSrc

private

Definition at line 55 of file RecoSelector.h.

std::string RecoSelector::name

private

Definition at line 51 of file RecoSelector.h.

Referenced by ElectronMVAID.ElectronMVAID::__call__(), dirstructure.Directory::__create_pie_image(), dqm_interfaces.DirID::__eq__(), dirstructure.Directory::__get_full_path(), dirstructure.Comparison::__get_img_name(), dataset.Dataset::__getDataType(), dataset.Dataset::__getFileInfoList(), cuy.divideElement::__init__(), cuy.plotElement::__init__(), cuy.additionElement::__init__(), cuy.superimposeElement::__init__(), cuy.graphElement::__init__(), dirstructure.Comparison::__make_image(), core.autovars.NTupleVariable::__repr__(), core.autovars.NTupleObjectType::__repr__(), core.autovars.NTupleObject::__repr__(), core.autovars.NTupleCollection::__repr__(), dirstructure.Directory::__repr__(), dqm_interfaces.DirID::__repr__(), dirstructure.Comparison::__repr__(), config.CFG::__str__(), counter.Counter::__str__(), average.Average::__str__(), core.autovars.NTupleObjectType::addSubObjects(), core.autovars.NTupleObjectType::addVariables(), core.autovars.NTupleObjectType::allVars(), dirstructure.Directory::calcStats(), validation.Sample::digest(), python.rootplot.utilities.Hist::divide(), python.rootplot.utilities.Hist::divide_wilson(), core.autovars.NTupleVariable::fillBranch(), core.autovars.NTupleObject::fillBranches(), core.autovars.NTupleCollection::fillBranchesScalar(), core.autovars.NTupleCollection::fillBranchesVector(), core.autovars.NTupleCollection::get_cpp_declaration(), core.autovars.NTupleCollection::get_cpp_wrapper_class(), core.autovars.NTupleCollection::get_py_wrapper_class(), utils.StatisticalTest::get_status(), production_tasks.Task::getname(), dataset.CMSDataset::getPrimaryDatasetEntries(), dataset.PrivateDataset::getPrimaryDatasetEntries(), VIDSelectorBase.VIDSelectorBase::initialize(), core.autovars.NTupleVariable::makeBranch(), core.autovars.NTupleObject::makeBranches(), core.autovars.NTupleCollection::makeBranchesScalar(), core.autovars.NTupleCollection::makeBranchesVector(), dirstructure.Directory::print_report(), dataset.BaseDataset::printInfo(), dataset.Dataset::printInfo(), production_tasks.MonitorJobs::run(), python.rootplot.utilities.Hist::TGraph(), python.rootplot.utilities.Hist::TH1F(), Vispa.Views.PropertyView.Property::valueChanged(), counter.Counter::write(), and average.Average::write().

double RecoSelector::reco_metMin

private

Definition at line 59 of file RecoSelector.h.

double RecoSelector::reco_ptElecMin

private

Definition at line 62 of file RecoSelector.h.

double RecoSelector::reco_ptJet1Min

private

Definition at line 60 of file RecoSelector.h.

double RecoSelector::reco_ptJet2Min

private

Definition at line 61 of file RecoSelector.h.

double RecoSelector::reco_ptMuonMin

private

Definition at line 63 of file RecoSelector.h.

double RecoSelector::reco_ptPhotMin

private

Definition at line 64 of file RecoSelector.h.

const reco::CaloJetCollection* RecoSelector::theCaloJetCollection

private

Definition at line 69 of file RecoSelector.h.

const reco::CaloMETCollection* RecoSelector::theCaloMETCollection

private

Definition at line 70 of file RecoSelector.h.

const reco::GsfElectronCollection* RecoSelector::theElectronCollection

private

Definition at line 66 of file RecoSelector.h.

const reco::MuonCollection* RecoSelector::theMuonCollection

private

Definition at line 67 of file RecoSelector.h.

const reco::PhotonCollection* RecoSelector::thePhotonCollection

private

Definition at line 68 of file RecoSelector.h.