#include <Configuration/Skimming/src/LeptonSkimming.cc>

Inheritance diagram for LeptonSkimming:

Public Member Functions
	LeptonSkimming (const edm::ParameterSet &)

	~LeptonSkimming () override

Public Member Functions inherited from edm::stream::EDFilter<>
	EDFilter ()=default

	EDFilter (const EDFilter &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDFilter &	operator= (const EDFilter &)=delete

Public Member Functions inherited from edm::stream::EDFilterBase
	EDFilterBase ()

	EDFilterBase (const EDFilterBase &)=delete

ModuleDescription const &	moduleDescription () const

const EDFilterBase &	operator= (const EDFilterBase &)=delete

	~EDFilterBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector < edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector < edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::stream::EDFilterBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	beginStream (edm::StreamID) override

void	endStream () override

bool	filter (edm::Event &, const edm::EventSetup &) override

bool	hltFired (const edm::Event &iEvent, const edm::EventSetup &iSetup, std::vector< string > HLTPath)

std::array< float, 5 >	hltObject (const edm::Event &iEvent, const edm::EventSetup &iSetup, std::vector< string > Seed)

Private Attributes
float	beam_x

float	beam_y

float	beam_z

edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_

edm::ESGetToken< MagneticField, IdealMagneticFieldRecord >	bFieldToken_

double	BiasedWP = -100

std::vector< std::shared_ptr < reco::Track > >	cleanedObjTracks

std::vector< std::shared_ptr < reco::Track > >	cleanedPairTracks

std::vector< std::shared_ptr < reco::Track > >	cleanedTracks

edm::EDGetTokenT < reco::ConversionCollection >	conversionsToken_

double	Cosee_Cut = -1

double	CoseeK_Cut = 0

bool	EarlyStop = false

std::vector< float >	el_eta

std::vector< float >	el_phi

std::vector< float >	el_pt

edm::EDGetTokenT < edm::ValueMap< float > >	eleBWPToken_

edm::EDGetToken	electronsToken_

edm::EDGetTokenT < edm::ValueMap< float > >	eleUnBWPToken_

std::vector< std::shared_ptr < reco::Track > >	ElTracks

double	ElTrgExclusionCone = 0

double	ElTrgMuDz_Cut = 1000

std::vector< unsigned int >	Epair_ObjectId

std::vector< unsigned int >	Epair_ObjectIndex

std::vector< unsigned int >	Epair_TrkIndex

double	EpairZvtx_Cut = 10000000

double	EtaTrack_Cut = 10

std::vector< string >	HLTFilter_

std::vector< string >	HLTPath_

double	Ksdxy_Cut = 1000

double	MaxChi2Track_Cut = 1000

double	MaxMB_Cut = 1000

double	MaxMee_Cut = 1000

double	MinChi2Track_Cut = -1000

double	MinMB_Cut = 0

double	MinMee_Cut = 0

std::vector< float >	muon_eta

std::vector< float >	muon_medium

std::vector< float >	muon_phi

std::vector< float >	muon_pt

std::vector< float >	muon_soft

std::vector< float >	muon_tight

edm::EDGetToken	muonsToken_

std::vector< std::shared_ptr < reco::Track > >	MuTracks

double	MuTrgExclusionCone = 0

double	MuTrgMatchCone = 1000

double	MuTrgMuDz_Cut = 1000

double	MuTrkMinDR_Cut = 0

unsigned int	nel = 0

unsigned int	nmuons = 0

unsigned int	ntracks = 0

std::vector< unsigned int >	object_container

std::vector< unsigned int >	object_id

bool	ObjPtLargerThanTrack = false

double	Probee_Cut = 0

double	ProbeeK_Cut = 0

double	PtB_Cut = 0

double	PtEl_Cut = 0

double	PtKTrack_Cut = 0

double	PtMu_Cut = 0

double	PtTrack_Cut = 0

double	QualMu_Cut = 0

bool	Result = false

bool	SaveOnlyEPairTracks = false

bool	SaveOnlyTracks = false

bool	SaveOutputRoot = true

float	SelectedMu_DR = -1

int	SelectedMu_index

std::vector< float >	SelectedTrgObj_PtEtaPhiCharge

bool	SkimOnlyElectrons = false

bool	SkimOnlyMuons = false

bool	SkipIfNoMuMatch = false

double	SLxy_Cut = 0

std::vector< float >	tempPtEtaPhiM

std::vector< float >	tempXYZ

int	test_ev = 0

std::vector< float >	track_vx

std::vector< float >	track_vy

std::vector< float >	track_vz

double	TrackMuDz_Cut = 100000000

std::vector< unsigned int >	trackObj_container

std::vector< unsigned int >	trackPair_container

edm::EDGetToken	Tracks_

double	TrackSdxy_Cut = 0

double	TrgExclusionCone = -1

edm::EDGetTokenT < edm::TriggerResults >	trgresultsToken_

edm::EDGetTokenT < trigger::TriggerEvent >	trigobjectsToken_

std::vector< unsigned int >	Trk_container

unsigned int	trk_index = 0

double	TrkObjExclusionCone = 0

double	UnbiasedWP = -100

bool	UseOnlyBKeeMCForTriplets = false

reco::TrackBase::Point	vertex_point

edm::EDGetTokenT < reco::VertexCollection >	vtxToken_

float	ZvertexTrg = -1 * std::numeric_limits<float>::max()

Additional Inherited Members
Public Types inherited from edm::stream::EDFilter<>
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Public Types inherited from edm::stream::EDFilterBase
typedef EDFilterAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 84 of file LeptonSkimming.h.

Constructor & Destructor Documentation

LeptonSkimming::LeptonSkimming ( const edm::ParameterSet & iConfig )

explicit

Definition at line 26 of file LeptonSkimming.cc.

References BiasedWP, Cosee_Cut, CoseeK_Cut, ElTrgExclusionCone, ElTrgMuDz_Cut, EpairZvtx_Cut, EtaTrack_Cut, edm::ParameterSet::getParameter(), Ksdxy_Cut, MaxChi2Track_Cut, MaxMB_Cut, MaxMee_Cut, MinChi2Track_Cut, MinMB_Cut, MinMee_Cut, MuTrgExclusionCone, MuTrgMatchCone, MuTrgMuDz_Cut, MuTrkMinDR_Cut, ObjPtLargerThanTrack, Probee_Cut, ProbeeK_Cut, PtB_Cut, PtEl_Cut, PtKTrack_Cut, PtMu_Cut, PtTrack_Cut, QualMu_Cut, environment_file_cff::runParameters, SkimOnlyElectrons, SkimOnlyMuons, SkipIfNoMuMatch, SLxy_Cut, TrackMuDz_Cut, TrackSdxy_Cut, TrgExclusionCone, TrkObjExclusionCone, and UnbiasedWP.

     : electronsToken_(consumes<std::vector<reco::GsfElectron>>(iConfig.getParameter<edm::InputTag>("electrons"))),
       eleBWPToken_(consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("eleBiasedWP"))),
       eleUnBWPToken_(consumes<edm::ValueMap<float>>(iConfig.getParameter<edm::InputTag>("eleUnbiasedWP"))),
       muonsToken_(consumes<std::vector<reco::Muon>>(iConfig.getParameter<edm::InputTag>("muons"))),
       Tracks_(consumes<std::vector<reco::Track>>(iConfig.getParameter<edm::InputTag>("tracks"))),
       vtxToken_(consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("vertices"))),
       beamSpotToken_(consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpot"))),
       conversionsToken_(consumes<reco::ConversionCollection>(iConfig.getParameter<edm::InputTag>("conversions"))),
 
       trgresultsToken_(consumes<edm::TriggerResults>(iConfig.getParameter<edm::InputTag>("triggerresults"))),
       trigobjectsToken_(consumes<trigger::TriggerEvent>(iConfig.getParameter<edm::InputTag>("triggerobjects"))),
       bFieldToken_(esConsumes()),
       HLTFilter_(iConfig.getParameter<vector<string>>("HLTFilter")),
       HLTPath_(iConfig.getParameter<vector<string>>("HLTPath")) {
   edm::ParameterSet runParameters = iConfig.getParameter<edm::ParameterSet>("RunParameters");
   //mu matching with trigger
   MuTrgMatchCone = runParameters.getParameter<double>("MuTrgMatchCone");
   SkipIfNoMuMatch = runParameters.getParameter<bool>("SkipIfNoMuMatch");
   //track cuts
   PtTrack_Cut = runParameters.getParameter<double>("PtTrack_Cut");
   EtaTrack_Cut = runParameters.getParameter<double>("EtaTrack_Cut");
   MinChi2Track_Cut = runParameters.getParameter<double>("MinChi2Track_Cut");
   MaxChi2Track_Cut = runParameters.getParameter<double>("MaxChi2Track_Cut");
   MuTrkMinDR_Cut = runParameters.getParameter<double>("MuTrkMinDR_Cut");
   TrackMuDz_Cut = runParameters.getParameter<double>("TrackMuDz_Cut");
   TrgExclusionCone = runParameters.getParameter<double>("TrgExclusionCone");
   //lepton cuts
   PtMu_Cut = runParameters.getParameter<double>("PtMu_Cut");
   PtEl_Cut = runParameters.getParameter<double>("PtEl_Cut");
   QualMu_Cut = runParameters.getParameter<double>("QualMu_Cut");
   MuTrgExclusionCone = runParameters.getParameter<double>("MuTrgExclusionCone");
   ElTrgExclusionCone = runParameters.getParameter<double>("ElTrgExclusionCone");
   TrkObjExclusionCone = runParameters.getParameter<double>("TrkObjExclusionCone");
   MuTrgMuDz_Cut = runParameters.getParameter<double>("MuTrgMuDz_Cut");
   ElTrgMuDz_Cut = runParameters.getParameter<double>("ElTrgMuDz_Cut");
   BiasedWP = runParameters.getParameter<double>("BiasedWP");
   UnbiasedWP = runParameters.getParameter<double>("UnbiasedWP");
   SkimOnlyMuons = runParameters.getParameter<bool>("SkimOnlyMuons");
   SkimOnlyElectrons = runParameters.getParameter<bool>("SkimOnlyElectrons");
   //pair cuts
   MaxMee_Cut = runParameters.getParameter<double>("MaxMee_Cut");
   MinMee_Cut = runParameters.getParameter<double>("MinMee_Cut");
   Probee_Cut = runParameters.getParameter<double>("Probee_Cut");
   Cosee_Cut = runParameters.getParameter<double>("Cosee_Cut");
   EpairZvtx_Cut = runParameters.getParameter<double>("EpairZvtx_Cut");
   //kaon
   PtKTrack_Cut = runParameters.getParameter<double>("PtKTrack_Cut");
   TrackSdxy_Cut = runParameters.getParameter<double>("TrackSdxy_Cut");
   Ksdxy_Cut = runParameters.getParameter<double>("Ksdxy_Cut");
   //triplet
   MaxMB_Cut = runParameters.getParameter<double>("MaxMB_Cut");
   MinMB_Cut = runParameters.getParameter<double>("MinMB_Cut");
   ProbeeK_Cut = runParameters.getParameter<double>("ProbeeK_Cut");
   CoseeK_Cut = runParameters.getParameter<double>("CoseeK_Cut");
   SLxy_Cut = runParameters.getParameter<double>("SLxy_Cut");
   PtB_Cut = runParameters.getParameter<double>("PtB_Cut");
 
   ObjPtLargerThanTrack = runParameters.getParameter<bool>("ObjPtLargerThanTrack");
 }

LeptonSkimming::~LeptonSkimming ( )

override

Definition at line 87 of file LeptonSkimming.cc.

                                 {
   // do anything here that needs to be done at destruction time
   // (e.g. close files, deallocate resources etc.)
 }

Member Function Documentation

void LeptonSkimming::beginStream ( edm::StreamID )

overrideprivatevirtual

Reimplemented from edm::stream::EDFilterBase.

Definition at line 529 of file LeptonSkimming.cc.

529 {}

void LeptonSkimming::endStream ( )

overrideprivatevirtual

Reimplemented from edm::stream::EDFilterBase.

Definition at line 532 of file LeptonSkimming.cc.

532 {}

void LeptonSkimming::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 567 of file LeptonSkimming.cc.

References edm::ConfigurationDescriptions::addDefault(), submitPVResolutionJobs::desc, and edm::ParameterSetDescription::setUnknown().

                                                                                 {
   //The following says we do not know what parameters are allowed so do no validation
   // Please change this to state exactly what you do use, even if it is no parameters
   edm::ParameterSetDescription desc;
   desc.setUnknown();
   descriptions.addDefault(desc);
 }

bool LeptonSkimming::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivatevirtual

Implements edm::stream::EDFilterBase.

Definition at line 168 of file LeptonSkimming.cc.

References beam_x, beam_y, beam_z, beamSpotToken_, ecalTB2006H4_GenSimDigiReco_cfg::bField, bFieldToken_, BiasedWP, ChiSquaredProbability(), cleanedObjTracks, cleanedPairTracks, cleanedTracks, conversions_cfi::conversions, conversionsToken_, Cosee_Cut, CoseeK_Cut, HLT_FULL_cff::deltaR, el_eta, el_phi, el_pt, eleBWPToken_, HI_PhotonSkim_cff::electrons, electronsToken_, eleUnBWPToken_, ElTracks, ElTrgExclusionCone, ElTrgMuDz_Cut, Epair_ObjectId, Epair_ObjectIndex, Epair_TrkIndex, EpairZvtx_Cut, submitPVValidationJobs::err, reco::TrackBase::eta(), EtaTrack_Cut, edm::Event::getByToken(), edm::EventSetup::getData(), ConversionTools::hasMatchedConversion(), HLTFilter_, hltFired(), hltObject(), HLTPath_, isMediumMuonCustom(), edm::Ref< C, T, F >::isNull(), muon::isSoftMuon(), isTightMuonCustom(), Ksdxy_Cut, SiStripPI::max, MaxChi2Track_Cut, MaxMB_Cut, MaxMee_Cut, MinChi2Track_Cut, MinMB_Cut, MinMee_Cut, RPCpg::mu, muon_eta, muon_medium, muon_phi, muon_pt, muon_soft, muon_tight, patZpeak::muons, muonsToken_, MuTracks, MuTrgExclusionCone, MuTrgMuDz_Cut, nel, nmuons, ntracks, getGTfromDQMFile::obj, object_container, object_id, ObjPtLargerThanTrack, PV3DBase< T, PVType, FrameType >::perp(), Probee_Cut, ProbeeK_Cut, reco::TrackBase::pt(), PtB_Cut, PtEl_Cut, PtKTrack_Cut, PtMu_Cut, PtTrack_Cut, QualMu_Cut, dt_dqm_sourceclient_common_cff::reco, GlobalErrorBase< T, ErrorWeightType >::rerr(), Result, fileCollector::seed, SelectedMu_DR, SelectedMu_index, SelectedTrgObj_PtEtaPhiCharge, SkimOnlyElectrons, SkimOnlyMuons, SkipIfNoMuMatch, SLxy_Cut, mathSSE::sqrt(), test_ev, TrackMuDz_Cut, tracks, Tracks_, TrackSdxy_Cut, TrgExclusionCone, trgresultsToken_, Trk_container, trk_index, TrkObjExclusionCone, UnbiasedWP, KalmanVertexFitter::vertex(), vertex_point, beam_dqm_sourceclient-live_cfg::vertices, vtxToken_, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and ZvertexTrg.

                                                                          {
   using namespace std;
   using namespace edm;
   using namespace reco;
   using namespace trigger;
   // using namespace PhysicsTools;
 
   test_ev++;
   //Get a few collections to apply basic electron ID
   //Get electrons
   edm::Handle<std::vector<reco::GsfElectron>> electrons;
   iEvent.getByToken(electronsToken_, electrons);
   //get bdt values
   edm::Handle<edm::ValueMap<float>> ele_mva_wp_biased;
   iEvent.getByToken(eleBWPToken_, ele_mva_wp_biased);
   edm::Handle<edm::ValueMap<float>> ele_mva_wp_unbiased;
   iEvent.getByToken(eleUnBWPToken_, ele_mva_wp_unbiased);
   //get muons
   edm::Handle<std::vector<reco::Muon>> muons;
   iEvent.getByToken(muonsToken_, muons);
   //Get conversions
   edm::Handle<reco::ConversionCollection> conversions;
   iEvent.getByToken(conversionsToken_, conversions);
   // Get the beam spot
   edm::Handle<reco::BeamSpot> theBeamSpot;
   iEvent.getByToken(beamSpotToken_, theBeamSpot);
   //Get vertices
   edm::Handle<std::vector<reco::Vertex>> vertices;
   iEvent.getByToken(vtxToken_, vertices);
   //continue if there are no vertices
   if (vertices->empty())
     return false;
   edm::Handle<vector<reco::Track>> tracks;
   iEvent.getByToken(Tracks_, tracks);
   edm::Handle<edm::TriggerResults> trigResults;
   iEvent.getByToken(trgresultsToken_, trigResults);
   auto const& bField = iSetup.getData(bFieldToken_);
   KalmanVertexFitter theKalmanFitter(false);
   TransientVertex LLvertex;
 
   // trigger1=0; trigger2=0; trigger3=0; trigger4=0; trigger5=0; trigger6=0;
   nmuons = 0;
   nel = 0;
   ntracks = 0;
 
   SelectedMu_index = -1;
   SelectedMu_DR = std::numeric_limits<float>::max();
   muon_pt.clear();
   muon_eta.clear();
   muon_phi.clear();
   Result = false;
   el_pt.clear();
   el_eta.clear();
   el_phi.clear();
   Trk_container.clear();
   MuTracks.clear();
   ElTracks.clear();
   object_container.clear();
   object_id.clear();
   cleanedTracks.clear();
   Epair_ObjectId.clear();
   muon_soft.clear();
   muon_medium.clear();
   muon_tight.clear();
   Epair_ObjectIndex.clear();
   cleanedObjTracks.clear();
   cleanedPairTracks.clear();
   Epair_ObjectIndex.clear();
   Epair_ObjectId.clear();
   Epair_TrkIndex.clear();
   //internal stuff
   ZvertexTrg = -1 * std::numeric_limits<float>::max();
   std::array<float, 5> SelectedTrgObj_PtEtaPhiCharge{{-999, -999, -999, -999, -999}};
   vertex_point.SetCoordinates(-1 * std::numeric_limits<float>::max(),
                               -1 * std::numeric_limits<float>::max(),
                               -1 * std::numeric_limits<float>::max());
   for (const reco::Vertex& vtx : *vertices) {
     bool isFake = vtx.isFake();
     if (isFake)
       continue;
     vertex_point.SetCoordinates(vtx.x(), vtx.y(), vtx.z());
     if (vertex_point.x() != -1 * std::numeric_limits<float>::max())
       break;
   }
   if (vertex_point.x() == -1 * std::numeric_limits<float>::max())
     return false;
 
   beam_x = theBeamSpot->x0();
   beam_y = theBeamSpot->y0();
   beam_z = theBeamSpot->z0();
   if (!hltFired(iEvent, iSetup, HLTPath_))
     return false;
 
   SelectedTrgObj_PtEtaPhiCharge = hltObject(iEvent, iSetup, HLTFilter_);
   SelectedMu_DR = std::numeric_limits<float>::max();
   MuTracks.clear();
   object_container.clear();
   object_id.clear();
   nmuons = 0;
   for (const reco::Muon& mu : *muons) {
     if (fabs(mu.eta()) > EtaTrack_Cut)
       continue;
     //find triggering muon
     float deltaRmu = deltaR(mu.eta(), mu.phi(), SelectedTrgObj_PtEtaPhiCharge[1], SelectedTrgObj_PtEtaPhiCharge[2]);
     if (deltaRmu < MuTrgMatchCone && SelectedMu_DR > deltaRmu) {
       SelectedMu_DR = deltaR(mu.eta(), mu.phi(), SelectedTrgObj_PtEtaPhiCharge[1], SelectedTrgObj_PtEtaPhiCharge[2]);
       ZvertexTrg = mu.vz();
     }
     //save muons that we want to skim
     if (!SkimOnlyElectrons) {
       bool tight = false, soft = false;
       if (vertices.isValid()) {
         tight = isTightMuonCustom(*(&mu), (*vertices)[0]);
         soft = muon::isSoftMuon(*(&mu), (*vertices)[0]);
       }
       const Track* mutrack = mu.bestTrack();
       muon_medium.push_back(isMediumMuonCustom(*(&mu)));
       muon_tight.push_back(tight);
       muon_soft.push_back(soft);
       muon_pt.push_back(mu.pt());
       muon_eta.push_back(mu.eta());
       muon_phi.push_back(mu.phi());
       auto muTrack = std::make_shared<reco::Track>(*mutrack);
       MuTracks.push_back(muTrack);
       object_container.push_back(nmuons);
       object_id.push_back(13);
       nmuons++;
     }
   }
 
   if (SelectedMu_DR == std::numeric_limits<float>::max() && SkipIfNoMuMatch) {
     return false;
   }
 
   //Save electrons we want to skim
   ElTracks.clear();
   if (!SkimOnlyMuons) {
     unsigned int count_el = -1;
     for (const reco::GsfElectron& el : *electrons) {
       count_el++;
       bool passConvVeto = !ConversionTools::hasMatchedConversion(*(&el), *conversions, theBeamSpot->position());
       if (!passConvVeto)
         continue;
       reco::GsfTrackRef seed = el.gsfTrack();
       if (seed.isNull())
         continue;
       if ((*ele_mva_wp_biased)[seed] < BiasedWP)
         continue;
       if ((*ele_mva_wp_unbiased)[seed] < UnbiasedWP)
         continue;
       const Track* eltrack = el.bestTrack();
       if (fabs(eltrack->eta()) > EtaTrack_Cut)
         continue;
       if (eltrack->pt() < PtEl_Cut)
         continue;
       auto ElTrack = std::make_shared<reco::Track>(*eltrack);
       ElTracks.push_back(ElTrack);
       object_container.push_back(nel);
       el_pt.push_back(eltrack->pt());
       el_eta.push_back(eltrack->eta());
       el_phi.push_back(eltrack->phi());
       nel++;
       object_id.push_back(11);
     }
   }
 
   //Save tracks we want to skim: used both as mu or e candidate
   cleanedTracks.clear();
   trk_index = 0;
   for (const reco::Track& trk : *tracks) {
     if (!trk.quality(Track::highPurity))
       continue;
     if (trk.pt() < PtTrack_Cut)
       continue;
     if (fabs(trk.eta()) > EtaTrack_Cut)
       continue;
     if (trk.charge() == 0)
       continue;
     if (trk.normalizedChi2() > MaxChi2Track_Cut || trk.normalizedChi2() < MinChi2Track_Cut)
       continue;
     if (fabs(trk.dxy()) / trk.dxyError() < TrackSdxy_Cut)
       continue;
     if (SelectedMu_DR < std::numeric_limits<float>::max()) {
       if (fabs(ZvertexTrg - trk.vz()) > TrackMuDz_Cut)
         continue;
       if (deltaR(trk.eta(), trk.phi(), SelectedTrgObj_PtEtaPhiCharge[1], SelectedTrgObj_PtEtaPhiCharge[2]) <
           TrgExclusionCone)
         continue;
     }
     //assignments
     auto cleanTrack = std::make_shared<reco::Track>(trk);
     cleanedTracks.push_back(cleanTrack);
     Trk_container.push_back(trk_index);
     trk_index++;
   }
 
   //create ll combination
 
   // fit track pairs
   cleanedObjTracks.clear();
   cleanedPairTracks.clear();
   TLorentzVector vel1, vel2;
   std::vector<std::shared_ptr<reco::Track>> cleanedObjects;
   //add tracks of objects
   if (!SkimOnlyElectrons) {
     for (auto& vec : MuTracks)
       cleanedObjects.push_back(vec);
   }
 
   if (!SkimOnlyMuons) {
     for (auto& vec : ElTracks)
       cleanedObjects.push_back(vec);
   }
 
   if (cleanedObjects.empty())
     return false;
 
   for (auto& obj : cleanedObjects) {
     auto tranobj = std::make_shared<reco::TransientTrack>(reco::TransientTrack(*obj, &bField));
     unsigned int iobj = &obj - &cleanedObjects[0];
     unsigned int index = object_container.at(iobj);
     float massLep = 0.0005;
     if (object_id.at(iobj) == 13)
       massLep = 0.105;
     //posible cut in mu quality
     if (object_id.at(iobj) == 13 && QualMu_Cut == 1 && !muon_soft.at(index))
       continue;
     if (object_id.at(iobj) == 13 && QualMu_Cut == 2 && !muon_medium.at(index))
       continue;
     if (object_id.at(iobj) == 13 && QualMu_Cut == 3 && !muon_tight.at(index))
       continue;
     // take the corresponding lepton track for lorentz vector
     vel1.SetPtEtaPhiM(obj->pt(), obj->eta(), obj->phi(), massLep);
     if (object_id.at(iobj) == 13 && vel1.Pt() < PtMu_Cut)
       continue;
     for (auto& trk2 : cleanedTracks) {
       unsigned int itrk2 = &trk2 - &cleanedTracks[0];
       //opposite sign
       if (obj->charge() * trk2->charge() == 1)
         continue;
       if (ObjPtLargerThanTrack && vel1.Pt() < trk2->pt())
         continue;
       vel2.SetPtEtaPhiM(trk2->pt(), trk2->eta(), trk2->phi(), massLep);
       //probe side cuts
       if (object_id.at(iobj) == 13 &&
           deltaR(obj->eta(), obj->phi(), SelectedTrgObj_PtEtaPhiCharge[1], SelectedTrgObj_PtEtaPhiCharge[2]) <
               MuTrgExclusionCone)
         continue;
       if (object_id.at(iobj) == 11 &&
           deltaR(obj->eta(), obj->phi(), SelectedTrgObj_PtEtaPhiCharge[1], SelectedTrgObj_PtEtaPhiCharge[2]) <
               ElTrgExclusionCone)
         continue;
       if (SelectedMu_DR < std::numeric_limits<float>::max()) {
         if (object_id.at(iobj) == 13 && fabs(ZvertexTrg - obj->vz()) > MuTrgMuDz_Cut)
           continue;
         if (object_id.at(iobj) == 11 && fabs(ZvertexTrg - obj->vz()) > ElTrgMuDz_Cut)
           continue;
       }
       //additional cuts
       float InvMassLepLep = (vel1 + vel2).M();
       if (InvMassLepLep > MaxMee_Cut || InvMassLepLep < MinMee_Cut)
         continue;
       auto trantrk2 = std::make_shared<reco::TransientTrack>(reco::TransientTrack(*trk2, &bField));
       std::vector<reco::TransientTrack> tempTracks;
       tempTracks.reserve(2);
       tempTracks.push_back(*tranobj);
       tempTracks.push_back(*trantrk2);
       LLvertex = theKalmanFitter.vertex(tempTracks);
       if (!LLvertex.isValid())
         continue;
       if (ChiSquaredProbability(LLvertex.totalChiSquared(), LLvertex.degreesOfFreedom()) < Probee_Cut)
         continue;
       if (ZvertexTrg > -1 * std::numeric_limits<float>::max() &&
           fabs(ZvertexTrg - LLvertex.position().z()) > EpairZvtx_Cut)
         continue;
       GlobalError err = LLvertex.positionError();
       GlobalPoint Dispbeamspot(-1 * ((theBeamSpot->x0() - LLvertex.position().x()) +
                                      (LLvertex.position().z() - theBeamSpot->z0()) * theBeamSpot->dxdz()),
                                -1 * ((theBeamSpot->y0() - LLvertex.position().y()) +
                                      (LLvertex.position().z() - theBeamSpot->z0()) * theBeamSpot->dydz()),
                                0);
       math::XYZVector pperp((vel1 + vel2).Px(), (vel1 + vel2).Py(), 0);
       math::XYZVector vperp(Dispbeamspot.x(), Dispbeamspot.y(), 0.);
       float tempCos = vperp.Dot(pperp) / (vperp.R() * pperp.R());
       if (tempCos < Cosee_Cut)
         continue;
       cleanedObjTracks.push_back(obj);
       cleanedPairTracks.push_back(trk2);
       Epair_ObjectIndex.push_back(object_container.at(iobj));
       Epair_ObjectId.push_back(object_id.at(iobj));
       Epair_TrkIndex.push_back(Trk_container.at(itrk2));
     }
   }
 
   // B recontrsuction
   TLorentzVector vK;
   for (unsigned int iobj = 0; iobj < cleanedObjTracks.size(); iobj++) {
     auto objtrk = cleanedObjTracks.at(iobj);
     auto pairtrk = cleanedPairTracks.at(iobj);
     auto tranobj = std::make_shared<reco::TransientTrack>(reco::TransientTrack(*objtrk, &bField));
     auto tranpair = std::make_shared<reco::TransientTrack>(reco::TransientTrack(*pairtrk, &bField));
     //unsigned int index=Epair_ObjectIndex.at(iobj);
     float massLep = 0.0005;
     if (Epair_ObjectId.at(iobj) == 13)
       massLep = 0.105;
     vel1.SetPtEtaPhiM(objtrk->pt(), objtrk->eta(), objtrk->phi(), massLep);
     vel2.SetPtEtaPhiM(pairtrk->pt(), pairtrk->eta(), pairtrk->phi(), massLep);
     for (auto& trk : cleanedTracks) {
       //reject track corresponding to mu or e
       if (trk->charge() == objtrk->charge() &&
           deltaR(objtrk->eta(), objtrk->phi(), trk->eta(), trk->phi()) < TrkObjExclusionCone)
         continue;
       if (trk->pt() < PtKTrack_Cut)
         continue;
       if (fabs(trk->dxy(vertex_point)) / trk->dxyError() < Ksdxy_Cut)
         continue;
       // skip the track that was used as mu or e
       if (Epair_TrkIndex.at(iobj) == &trk - &cleanedTracks[0])
         continue;
       vK.SetPtEtaPhiM(trk->pt(), trk->eta(), trk->phi(), 0.493);
       //final cuts
       float InvMass = (vel1 + vel2 + vK).M();
       if (InvMass > MaxMB_Cut || InvMass < MinMB_Cut)
         continue;
       if ((vel1 + vel2 + vK).Pt() < PtB_Cut)
         continue;
       auto trantrk = std::make_shared<reco::TransientTrack>(reco::TransientTrack(*trk, &bField));
       std::vector<reco::TransientTrack> tempTracks;
       tempTracks.reserve(3);
       tempTracks.push_back(*tranobj);
       tempTracks.push_back(*tranpair);
       tempTracks.push_back(*trantrk);
       LLvertex = theKalmanFitter.vertex(tempTracks);
       if (!LLvertex.isValid())
         continue;
       if (ChiSquaredProbability(LLvertex.totalChiSquared(), LLvertex.degreesOfFreedom()) < ProbeeK_Cut)
         continue;
       GlobalError err = LLvertex.positionError();
       GlobalPoint Dispbeamspot(-1 * ((theBeamSpot->x0() - LLvertex.position().x()) +
                                      (LLvertex.position().z() - theBeamSpot->z0()) * theBeamSpot->dxdz()),
                                -1 * ((theBeamSpot->y0() - LLvertex.position().y()) +
                                      (LLvertex.position().z() - theBeamSpot->z0()) * theBeamSpot->dydz()),
                                0);
       math::XYZVector pperp((vel1 + vel2 + vK).Px(), (vel1 + vel2 + vK).Py(), 0);
       math::XYZVector vperp(Dispbeamspot.x(), Dispbeamspot.y(), 0.);
       float tempCos = vperp.Dot(pperp) / (vperp.R() * pperp.R());
       if (tempCos < CoseeK_Cut)
         continue;
       if (SLxy_Cut > Dispbeamspot.perp() / sqrt(err.rerr(Dispbeamspot)))
         continue;
       Result = true;
       break;
     }
     if (Result)
       break;
   }
   //decission
   return Result;
 }

bool LeptonSkimming::hltFired	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		std::vector< string >	HLTPath
	)

private

Definition at line 92 of file LeptonSkimming.cc.

References edm::Event::getByToken(), dt_dqm_sourceclient_common_cff::reco, AlCaHLTBitMon_QueryRunRegistry::string, trgresultsToken_, edm::TriggerNames::triggerName(), and edm::Event::triggerNames().

Referenced by filter().

                                                                                                             {
   using namespace std;
   using namespace edm;
   using namespace reco;
   using namespace trigger;
 
   edm::Handle<edm::TriggerResults> trigResults;
   iEvent.getByToken(trgresultsToken_, trigResults);
   bool fire = false;
   if (trigResults.failedToGet())
     return false;
   for (unsigned int ip = 0; ip < HLTPath.size(); ip++) {
     int N_Triggers = trigResults->size();
     const edm::TriggerNames& trigName = iEvent.triggerNames(*trigResults);
     for (int i_Trig = 0; i_Trig < N_Triggers; ++i_Trig) {
       if (!trigResults->accept(i_Trig))
         continue;
       const std::string& TrigPath = trigName.triggerName(i_Trig);
       if (TrigPath.find(HLTPath[ip]) != std::string::npos)
         fire = true;
     }
   }
   return fire;
 }

std::array< float, 5 > LeptonSkimming::hltObject	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		std::vector< string >	Seed
	)

private

Definition at line 117 of file LeptonSkimming.cc.

References a, b, trigger::TriggerObject::eta(), edm::Event::getByToken(), trigger::TriggerObject::id(), HLT_FULL_cff::InputTag, dqmiolumiharvest::j, relativeConstraints::keys, trigger::TriggerObject::phi(), trigger::TriggerObject::pt(), dt_dqm_sourceclient_common_cff::reco, and trigobjectsToken_.

Referenced by filter().

                                                                        {
   using namespace std;
   using namespace edm;
   using namespace reco;
   using namespace trigger;
 
   edm::Handle<trigger::TriggerEvent> triggerObjectsSummary;
   iEvent.getByToken(trigobjectsToken_, triggerObjectsSummary);
   trigger::TriggerObjectCollection selectedObjects;
 
   std::vector<std::array<float, 5>> max_per_trigger;
 
   for (unsigned int ipath = 0; ipath < Seed.size(); ipath++) {
     std::vector<std::array<float, 5>> tot_tr_obj_pt_eta_phi;
     if (!triggerObjectsSummary.isValid())
       continue;
     size_t filterIndex = (*triggerObjectsSummary).filterIndex(InputTag(Seed[ipath], "", "HLT"));
     trigger::TriggerObjectCollection allTriggerObjects = triggerObjectsSummary->getObjects();
     if (filterIndex < (*triggerObjectsSummary).sizeFilters()) {
       const trigger::Keys& keys = (*triggerObjectsSummary).filterKeys(filterIndex);
       for (size_t j = 0; j < keys.size(); j++) {
         const trigger::TriggerObject& foundObject = (allTriggerObjects)[keys[j]];
         std::array<float, 5> tr_obj_pt_eta_phi;
         if (fabs(foundObject.id()) != 13)
           continue;
         tr_obj_pt_eta_phi[0] = foundObject.pt();
         tr_obj_pt_eta_phi[1] = foundObject.eta();
         tr_obj_pt_eta_phi[2] = foundObject.phi();
         tr_obj_pt_eta_phi[3] = foundObject.id() / fabs(foundObject.id());
         tot_tr_obj_pt_eta_phi.push_back(tr_obj_pt_eta_phi);
       }
     }
     //take the max per hlt
     if (!tot_tr_obj_pt_eta_phi.empty()) {
       std::sort(tot_tr_obj_pt_eta_phi.begin(),
                 tot_tr_obj_pt_eta_phi.end(),
                 [](const std::array<float, 5>& a, const std::array<float, 5>& b) { return a[0] > b[0]; });
       max_per_trigger.push_back(tot_tr_obj_pt_eta_phi.at(0));
     }
   }
   //we know that at least a trigger fired
   //find the total max
   std::sort(max_per_trigger.begin(),
             max_per_trigger.end(),
             [](const std::array<float, 5>& a, const std::array<float, 5>& b) { return a[0] > b[0]; });
   return max_per_trigger.at(0);
 }

Member Data Documentation

float LeptonSkimming::beam_x

private

Definition at line 113 of file LeptonSkimming.h.

Referenced by filter().

float LeptonSkimming::beam_y

private

Definition at line 113 of file LeptonSkimming.h.

Referenced by filter().

float LeptonSkimming::beam_z

private

Definition at line 113 of file LeptonSkimming.h.

Referenced by filter().

edm::EDGetTokenT<reco::BeamSpot> LeptonSkimming::beamSpotToken_

private

Definition at line 105 of file LeptonSkimming.h.

Referenced by filter().

edm::ESGetToken<MagneticField, IdealMagneticFieldRecord> LeptonSkimming::bFieldToken_

private

Definition at line 109 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::BiasedWP = -100

private

Definition at line 174 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

std::vector<std::shared_ptr<reco::Track> > LeptonSkimming::cleanedObjTracks

private

Definition at line 127 of file LeptonSkimming.h.

Referenced by filter().

std::vector<std::shared_ptr<reco::Track> > LeptonSkimming::cleanedPairTracks

private

Definition at line 128 of file LeptonSkimming.h.

Referenced by filter().

std::vector<std::shared_ptr<reco::Track> > LeptonSkimming::cleanedTracks

private

Definition at line 131 of file LeptonSkimming.h.

Referenced by filter().

edm::EDGetTokenT<reco::ConversionCollection> LeptonSkimming::conversionsToken_

private

Definition at line 106 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::Cosee_Cut = -1

private

Definition at line 152 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::CoseeK_Cut = 0

private

Definition at line 159 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::EarlyStop = false

private

Definition at line 153 of file LeptonSkimming.h.

std::vector<float> LeptonSkimming::el_eta

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::el_phi

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::el_pt

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

edm::EDGetTokenT<edm::ValueMap<float> > LeptonSkimming::eleBWPToken_

private

Definition at line 100 of file LeptonSkimming.h.

Referenced by filter().

edm::EDGetToken LeptonSkimming::electronsToken_

private

Definition at line 99 of file LeptonSkimming.h.

Referenced by filter().

edm::EDGetTokenT<edm::ValueMap<float> > LeptonSkimming::eleUnBWPToken_

private

Definition at line 101 of file LeptonSkimming.h.

Referenced by filter().

std::vector<std::shared_ptr<reco::Track> > LeptonSkimming::ElTracks

private

Definition at line 130 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::ElTrgExclusionCone = 0

private

Definition at line 168 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::ElTrgMuDz_Cut = 1000

private

Definition at line 171 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

std::vector<unsigned int> LeptonSkimming::Epair_ObjectId

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

std::vector<unsigned int> LeptonSkimming::Epair_ObjectIndex

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

std::vector<unsigned int> LeptonSkimming::Epair_TrkIndex

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::EpairZvtx_Cut = 10000000

private

Definition at line 156 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::EtaTrack_Cut = 10

private

Definition at line 138 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

std::vector<string> LeptonSkimming::HLTFilter_

private

Definition at line 110 of file LeptonSkimming.h.

Referenced by filter().

std::vector<string> LeptonSkimming::HLTPath_

private

Definition at line 110 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::Ksdxy_Cut = 1000

private

Definition at line 157 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MaxChi2Track_Cut = 1000

private

Definition at line 140 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MaxMB_Cut = 1000

private

Definition at line 144 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MaxMee_Cut = 1000

private

Definition at line 142 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MinChi2Track_Cut = -1000

private

Definition at line 139 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MinMB_Cut = 0

private

Definition at line 145 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MinMee_Cut = 0

private

Definition at line 150 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

std::vector<float> LeptonSkimming::muon_eta

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_medium

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_phi

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_pt

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_soft

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_tight

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

edm::EDGetToken LeptonSkimming::muonsToken_

private

Definition at line 102 of file LeptonSkimming.h.

Referenced by filter().

std::vector<std::shared_ptr<reco::Track> > LeptonSkimming::MuTracks

private

Definition at line 129 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::MuTrgExclusionCone = 0

private

Definition at line 167 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MuTrgMatchCone = 1000

private

Definition at line 154 of file LeptonSkimming.h.

Referenced by LeptonSkimming().

double LeptonSkimming::MuTrgMuDz_Cut = 1000

private

Definition at line 170 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MuTrkMinDR_Cut = 0

private

Definition at line 141 of file LeptonSkimming.h.

Referenced by LeptonSkimming().

unsigned int LeptonSkimming::nel = 0

private

Definition at line 135 of file LeptonSkimming.h.

Referenced by filter().

unsigned int LeptonSkimming::nmuons = 0

private

Definition at line 135 of file LeptonSkimming.h.

Referenced by filter().

unsigned int LeptonSkimming::ntracks = 0

private

Definition at line 135 of file LeptonSkimming.h.

Referenced by filter().

std::vector<unsigned int> LeptonSkimming::object_container

private

Definition at line 133 of file LeptonSkimming.h.

Referenced by filter().

std::vector<unsigned int> LeptonSkimming::object_id

private

Definition at line 133 of file LeptonSkimming.h.

Referenced by filter().

bool LeptonSkimming::ObjPtLargerThanTrack = false

private

Definition at line 173 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::Probee_Cut = 0

private

Definition at line 151 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::ProbeeK_Cut = 0

private

Definition at line 158 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtB_Cut = 0

private

Definition at line 163 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtEl_Cut = 0

private

Definition at line 166 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtKTrack_Cut = 0

private

Definition at line 143 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtMu_Cut = 0

private

Definition at line 164 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtTrack_Cut = 0

private

Definition at line 137 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::QualMu_Cut = 0

private

Definition at line 165 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::Result = false

private

Definition at line 112 of file LeptonSkimming.h.

Referenced by filter().

bool LeptonSkimming::SaveOnlyEPairTracks = false

private

Definition at line 147 of file LeptonSkimming.h.

bool LeptonSkimming::SaveOnlyTracks = false

private

Definition at line 146 of file LeptonSkimming.h.

bool LeptonSkimming::SaveOutputRoot = true

private

Definition at line 172 of file LeptonSkimming.h.

float LeptonSkimming::SelectedMu_DR = -1

private

Definition at line 122 of file LeptonSkimming.h.

Referenced by filter().

int LeptonSkimming::SelectedMu_index

private

Definition at line 120 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::SelectedTrgObj_PtEtaPhiCharge

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

bool LeptonSkimming::SkimOnlyElectrons = false

private

Definition at line 177 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::SkimOnlyMuons = false

private

Definition at line 176 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::SkipIfNoMuMatch = false

private

Definition at line 155 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::SLxy_Cut = 0

private

Definition at line 162 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

std::vector<float> LeptonSkimming::tempPtEtaPhiM

private

Definition at line 126 of file LeptonSkimming.h.

std::vector<float> LeptonSkimming::tempXYZ

private

Definition at line 126 of file LeptonSkimming.h.

int LeptonSkimming::test_ev = 0

private

Definition at line 178 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::track_vx

private

Definition at line 115 of file LeptonSkimming.h.

std::vector<float> LeptonSkimming::track_vy

private

Definition at line 115 of file LeptonSkimming.h.

std::vector<float> LeptonSkimming::track_vz

private

Definition at line 115 of file LeptonSkimming.h.

double LeptonSkimming::TrackMuDz_Cut = 100000000

private

Definition at line 160 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

std::vector<unsigned int> LeptonSkimming::trackObj_container

private

Definition at line 132 of file LeptonSkimming.h.

std::vector<unsigned int> LeptonSkimming::trackPair_container

private

Definition at line 132 of file LeptonSkimming.h.

edm::EDGetToken LeptonSkimming::Tracks_

private

Definition at line 103 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::TrackSdxy_Cut = 0

private

Definition at line 148 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::TrgExclusionCone = -1

private

Definition at line 161 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

edm::EDGetTokenT<edm::TriggerResults> LeptonSkimming::trgresultsToken_

private

Definition at line 107 of file LeptonSkimming.h.

Referenced by filter(), and hltFired().

edm::EDGetTokenT<trigger::TriggerEvent> LeptonSkimming::trigobjectsToken_

private

Definition at line 108 of file LeptonSkimming.h.

Referenced by hltObject().

std::vector<unsigned int> LeptonSkimming::Trk_container

private

Definition at line 133 of file LeptonSkimming.h.

Referenced by filter().

unsigned int LeptonSkimming::trk_index = 0

private

Definition at line 124 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::TrkObjExclusionCone = 0

private

Definition at line 169 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::UnbiasedWP = -100

private

Definition at line 175 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::UseOnlyBKeeMCForTriplets = false

private

Definition at line 149 of file LeptonSkimming.h.

reco::TrackBase::Point LeptonSkimming::vertex_point

private

Definition at line 125 of file LeptonSkimming.h.

Referenced by filter().

edm::EDGetTokenT<reco::VertexCollection> LeptonSkimming::vtxToken_

private

Definition at line 104 of file LeptonSkimming.h.

Referenced by filter().

float LeptonSkimming::ZvertexTrg = -1 * std::numeric_limits<float>::max()

private

Definition at line 123 of file LeptonSkimming.h.

Referenced by filter().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation