#include <SkimmingForB/LeptonSkimming/plugins/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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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 &	operator= (EDConsumerBase &&)=default

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

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

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_

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<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	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
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

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

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

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...

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

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

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

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 85 of file LeptonSkimming.h.

Constructor & Destructor Documentation

◆ LeptonSkimming()

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

explicit

Definition at line 26 of file LeptonSkimming.cc.

     : 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"))),
       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");
 }

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.

◆ ~LeptonSkimming()

LeptonSkimming::~LeptonSkimming ( )

override

Definition at line 86 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

◆ beginStream()

void LeptonSkimming::beginStream ( edm::StreamID )

overrideprivatevirtual

Reimplemented from edm::stream::EDFilterBase.

Definition at line 529 of file LeptonSkimming.cc.

529 {}

◆ endStream()

void LeptonSkimming::endStream ( )

overrideprivatevirtual

Reimplemented from edm::stream::EDFilterBase.

Definition at line 532 of file LeptonSkimming.cc.

532 {}

◆ fillDescriptions()

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

static

Definition at line 567 of file LeptonSkimming.cc.

                                                                                 {
   //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);
 }

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

◆ filter()

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

overrideprivatevirtual

Implements edm::stream::EDFilterBase.

Definition at line 167 of file LeptonSkimming.cc.

                                                                          {
   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);
   edm::ESHandle<MagneticField> bFieldHandle;
   iSetup.get<IdealMagneticFieldRecord>().get(bFieldHandle);
   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, &(*bFieldHandle)));
     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, &(*bFieldHandle)));
       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, &(*bFieldHandle)));
     auto tranpair = std::make_shared<reco::TransientTrack>(reco::TransientTrack(*pairtrk, &(*bFieldHandle)));
     //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, &(*bFieldHandle)));
       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;
 }

References beam_x, beam_y, beam_z, beamSpotToken_, BiasedWP, ChiSquaredProbability(), cleanedObjTracks, cleanedPairTracks, cleanedTracks, pwdgSkimBPark_cfi::conversions, conversionsToken_, Cosee_Cut, CoseeK_Cut, PbPb_ZMuSkimMuonDPG_cff::deltaR, reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), el_eta, el_phi, el_pt, eleBWPToken_, pwdgSkimBPark_cfi::electrons, electronsToken_, eleUnBWPToken_, ElTracks, ElTrgExclusionCone, ElTrgMuDz_Cut, Epair_ObjectId, Epair_ObjectIndex, Epair_TrkIndex, EpairZvtx_Cut, submitPVResolutionJobs::err, reco::TrackBase::eta(), EtaTrack_Cut, edm::EventSetup::get(), get, ConversionTools::hasMatchedConversion(), muons_cff::highPurity, HLTFilter_, hltFired(), hltObject(), HLTPath_, iEvent, isMediumMuonCustom(), muon::isSoftMuon(), isTightMuonCustom(), TransientVertex::isValid(), Ksdxy_Cut, SiStripPI::max, MaxChi2Track_Cut, MaxMB_Cut, MaxMee_Cut, MinChi2Track_Cut, MinMB_Cut, MinMee_Cut, amptDefaultParameters_cff::mu, muon_eta, muon_medium, muon_phi, muon_pt, muon_soft, muon_tight, PDWG_BPHSkim_cff::muons, muonsToken_, MuTracks, MuTrgExclusionCone, MuTrgMuDz_Cut, nel, nmuons, ntracks, getGTfromDQMFile::obj, object_container, object_id, ObjPtLargerThanTrack, PV3DBase< T, PVType, FrameType >::perp(), reco::TrackBase::phi(), reco::BeamSpot::position(), TransientVertex::position(), TransientVertex::positionError(), Probee_Cut, ProbeeK_Cut, reco::TrackBase::pt(), PtB_Cut, PtEl_Cut, PtKTrack_Cut, PtMu_Cut, PtTrack_Cut, QualMu_Cut, 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, pwdgSkimBPark_cfi::vertices, extraflags_cff::vtx, vtxToken_, PV3DBase< T, PVType, FrameType >::x(), reco::BeamSpot::x0(), PV3DBase< T, PVType, FrameType >::y(), reco::BeamSpot::y0(), PV3DBase< T, PVType, FrameType >::z(), reco::BeamSpot::z0(), and ZvertexTrg.

◆ hltFired()

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

private

Definition at line 91 of file LeptonSkimming.cc.

                                                                                                             {
   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;
 }

References edm::HLTGlobalStatus::accept(), edm::HandleBase::failedToGet(), pwdgSkimBPark_cfi::HLTPath, iEvent, edm::HLTGlobalStatus::size(), AlCaHLTBitMon_QueryRunRegistry::string, trgresultsToken_, and EgHLTOffTrigSelection_cfi::trigName.

Referenced by filter().

◆ hltObject()

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

private

Definition at line 116 of file LeptonSkimming.cc.

                                                                        {
   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);
 }

References a, b, trigger::TriggerObject::eta(), trigger::TriggerEvent::getObjects(), trigger::TriggerObject::id(), iEvent, HLT_FULL_cff::InputTag, edm::HandleBase::isValid(), dqmiolumiharvest::j, relativeConstraints::keys, trigger::TriggerObject::phi(), trigger::TriggerObject::pt(), jetUpdater_cfi::sort, and trigobjectsToken_.

Referenced by filter().

Member Data Documentation

◆ beam_x

float LeptonSkimming::beam_x

private

Definition at line 113 of file LeptonSkimming.h.

Referenced by filter().

◆ beam_y

float LeptonSkimming::beam_y

private

Definition at line 113 of file LeptonSkimming.h.

Referenced by filter().

◆ beam_z

float LeptonSkimming::beam_z

private

Definition at line 113 of file LeptonSkimming.h.

Referenced by filter().

◆ beamSpotToken_

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

private

Definition at line 106 of file LeptonSkimming.h.

Referenced by filter().

◆ BiasedWP

double LeptonSkimming::BiasedWP = -100

private

Definition at line 174 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ cleanedObjTracks

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

private

Definition at line 127 of file LeptonSkimming.h.

Referenced by filter().

◆ cleanedPairTracks

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

private

Definition at line 128 of file LeptonSkimming.h.

Referenced by filter().

◆ cleanedTracks

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

private

Definition at line 131 of file LeptonSkimming.h.

Referenced by filter().

◆ conversionsToken_

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

private

Definition at line 107 of file LeptonSkimming.h.

Referenced by filter().

◆ Cosee_Cut

double LeptonSkimming::Cosee_Cut = -1

private

Definition at line 152 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ CoseeK_Cut

double LeptonSkimming::CoseeK_Cut = 0

private

Definition at line 159 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ EarlyStop

bool LeptonSkimming::EarlyStop = false

private

Definition at line 153 of file LeptonSkimming.h.

◆ el_eta

std::vector<float> LeptonSkimming::el_eta

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

◆ el_phi

std::vector<float> LeptonSkimming::el_phi

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

◆ el_pt

std::vector<float> LeptonSkimming::el_pt

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

◆ eleBWPToken_

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

private

Definition at line 101 of file LeptonSkimming.h.

Referenced by filter().

◆ electronsToken_

edm::EDGetToken LeptonSkimming::electronsToken_

private

Definition at line 100 of file LeptonSkimming.h.

Referenced by filter().

◆ eleUnBWPToken_

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

private

Definition at line 102 of file LeptonSkimming.h.

Referenced by filter().

◆ ElTracks

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

private

Definition at line 130 of file LeptonSkimming.h.

Referenced by filter().

◆ ElTrgExclusionCone

double LeptonSkimming::ElTrgExclusionCone = 0

private

Definition at line 168 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ ElTrgMuDz_Cut

double LeptonSkimming::ElTrgMuDz_Cut = 1000

private

Definition at line 171 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ Epair_ObjectId

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

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ Epair_ObjectIndex

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

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ Epair_TrkIndex

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

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ EpairZvtx_Cut

double LeptonSkimming::EpairZvtx_Cut = 10000000

private

Definition at line 156 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ EtaTrack_Cut

double LeptonSkimming::EtaTrack_Cut = 10

private

Definition at line 138 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ HLTFilter_

std::vector<string> LeptonSkimming::HLTFilter_

private

Definition at line 110 of file LeptonSkimming.h.

Referenced by filter().

◆ HLTPath_

std::vector<string> LeptonSkimming::HLTPath_

private

Definition at line 110 of file LeptonSkimming.h.

Referenced by filter().

◆ Ksdxy_Cut

double LeptonSkimming::Ksdxy_Cut = 1000

private

Definition at line 157 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MaxChi2Track_Cut

double LeptonSkimming::MaxChi2Track_Cut = 1000

private

Definition at line 140 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MaxMB_Cut

double LeptonSkimming::MaxMB_Cut = 1000

private

Definition at line 144 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MaxMee_Cut

double LeptonSkimming::MaxMee_Cut = 1000

private

Definition at line 142 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MinChi2Track_Cut

double LeptonSkimming::MinChi2Track_Cut = -1000

private

Definition at line 139 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MinMB_Cut

double LeptonSkimming::MinMB_Cut = 0

private

Definition at line 145 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MinMee_Cut

double LeptonSkimming::MinMee_Cut = 0

private

Definition at line 150 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ muon_eta

std::vector<float> LeptonSkimming::muon_eta

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_medium

std::vector<float> LeptonSkimming::muon_medium

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_phi

std::vector<float> LeptonSkimming::muon_phi

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_pt

std::vector<float> LeptonSkimming::muon_pt

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_soft

std::vector<float> LeptonSkimming::muon_soft

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_tight

std::vector<float> LeptonSkimming::muon_tight

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

◆ muonsToken_

edm::EDGetToken LeptonSkimming::muonsToken_

private

Definition at line 103 of file LeptonSkimming.h.

Referenced by filter().

◆ MuTracks

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

private

Definition at line 129 of file LeptonSkimming.h.

Referenced by filter().

◆ MuTrgExclusionCone

double LeptonSkimming::MuTrgExclusionCone = 0

private

Definition at line 167 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MuTrgMatchCone

double LeptonSkimming::MuTrgMatchCone = 1000

private

Definition at line 154 of file LeptonSkimming.h.

Referenced by LeptonSkimming().

◆ MuTrgMuDz_Cut

double LeptonSkimming::MuTrgMuDz_Cut = 1000

private

Definition at line 170 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MuTrkMinDR_Cut

double LeptonSkimming::MuTrkMinDR_Cut = 0

private

Definition at line 141 of file LeptonSkimming.h.

Referenced by LeptonSkimming().

◆ nel

unsigned int LeptonSkimming::nel = 0

private

Definition at line 135 of file LeptonSkimming.h.

Referenced by filter().

◆ nmuons

unsigned int LeptonSkimming::nmuons = 0

private

Definition at line 135 of file LeptonSkimming.h.

Referenced by filter().

◆ ntracks

unsigned int LeptonSkimming::ntracks = 0

private

Definition at line 135 of file LeptonSkimming.h.

Referenced by filter().

◆ object_container

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

private

Definition at line 133 of file LeptonSkimming.h.

Referenced by filter().

◆ object_id

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

private

Definition at line 133 of file LeptonSkimming.h.

Referenced by filter().

◆ ObjPtLargerThanTrack

bool LeptonSkimming::ObjPtLargerThanTrack = false

private

Definition at line 173 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ Probee_Cut

double LeptonSkimming::Probee_Cut = 0

private

Definition at line 151 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ ProbeeK_Cut

double LeptonSkimming::ProbeeK_Cut = 0

private

Definition at line 158 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtB_Cut

double LeptonSkimming::PtB_Cut = 0

private

Definition at line 163 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtEl_Cut

double LeptonSkimming::PtEl_Cut = 0

private

Definition at line 166 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtKTrack_Cut

double LeptonSkimming::PtKTrack_Cut = 0

private

Definition at line 143 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtMu_Cut

double LeptonSkimming::PtMu_Cut = 0

private

Definition at line 164 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtTrack_Cut

double LeptonSkimming::PtTrack_Cut = 0

private

Definition at line 137 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ QualMu_Cut

double LeptonSkimming::QualMu_Cut = 0

private

Definition at line 165 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ Result

bool LeptonSkimming::Result = false

private

Definition at line 112 of file LeptonSkimming.h.

Referenced by filter().

◆ SaveOnlyEPairTracks

bool LeptonSkimming::SaveOnlyEPairTracks = false

private

Definition at line 147 of file LeptonSkimming.h.

◆ SaveOnlyTracks

bool LeptonSkimming::SaveOnlyTracks = false

private

Definition at line 146 of file LeptonSkimming.h.

◆ SaveOutputRoot

bool LeptonSkimming::SaveOutputRoot = true

private

Definition at line 172 of file LeptonSkimming.h.

◆ SelectedMu_DR

float LeptonSkimming::SelectedMu_DR = -1

private

Definition at line 122 of file LeptonSkimming.h.

Referenced by filter().

◆ SelectedMu_index

int LeptonSkimming::SelectedMu_index

private

Definition at line 120 of file LeptonSkimming.h.

Referenced by filter().

◆ SelectedTrgObj_PtEtaPhiCharge

std::vector<float> LeptonSkimming::SelectedTrgObj_PtEtaPhiCharge

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

◆ SkimOnlyElectrons

bool LeptonSkimming::SkimOnlyElectrons = false

private

Definition at line 177 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ SkimOnlyMuons

bool LeptonSkimming::SkimOnlyMuons = false

private

Definition at line 176 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ SkipIfNoMuMatch

bool LeptonSkimming::SkipIfNoMuMatch = false

private

Definition at line 155 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ SLxy_Cut

double LeptonSkimming::SLxy_Cut = 0

private

Definition at line 162 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ tempPtEtaPhiM

std::vector<float> LeptonSkimming::tempPtEtaPhiM

private

Definition at line 126 of file LeptonSkimming.h.

◆ tempXYZ

std::vector<float> LeptonSkimming::tempXYZ

private

Definition at line 126 of file LeptonSkimming.h.

◆ test_ev

int LeptonSkimming::test_ev = 0

private

Definition at line 178 of file LeptonSkimming.h.

Referenced by filter().

◆ track_vx

std::vector<float> LeptonSkimming::track_vx

private

Definition at line 115 of file LeptonSkimming.h.

◆ track_vy

std::vector<float> LeptonSkimming::track_vy

private

Definition at line 115 of file LeptonSkimming.h.

◆ track_vz

std::vector<float> LeptonSkimming::track_vz

private

Definition at line 115 of file LeptonSkimming.h.

◆ TrackMuDz_Cut

double LeptonSkimming::TrackMuDz_Cut = 100000000

private

Definition at line 160 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ trackObj_container

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

private

Definition at line 132 of file LeptonSkimming.h.

◆ trackPair_container

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

private

Definition at line 132 of file LeptonSkimming.h.

◆ Tracks_

edm::EDGetToken LeptonSkimming::Tracks_

private

Definition at line 104 of file LeptonSkimming.h.

Referenced by filter().

◆ TrackSdxy_Cut

double LeptonSkimming::TrackSdxy_Cut = 0

private

Definition at line 148 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ TrgExclusionCone

double LeptonSkimming::TrgExclusionCone = -1

private

Definition at line 161 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ trgresultsToken_

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

private

Definition at line 108 of file LeptonSkimming.h.

Referenced by filter(), and hltFired().

◆ trigobjectsToken_

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

private

Definition at line 109 of file LeptonSkimming.h.

Referenced by hltObject().

◆ Trk_container

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

private

Definition at line 133 of file LeptonSkimming.h.

Referenced by filter().

◆ trk_index

unsigned int LeptonSkimming::trk_index = 0

private

Definition at line 124 of file LeptonSkimming.h.

Referenced by filter().

◆ TrkObjExclusionCone

double LeptonSkimming::TrkObjExclusionCone = 0

private

Definition at line 169 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ UnbiasedWP

double LeptonSkimming::UnbiasedWP = -100

private

Definition at line 175 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ UseOnlyBKeeMCForTriplets

bool LeptonSkimming::UseOnlyBKeeMCForTriplets = false

private

Definition at line 149 of file LeptonSkimming.h.

◆ vertex_point

reco::TrackBase::Point LeptonSkimming::vertex_point

private

Definition at line 125 of file LeptonSkimming.h.

Referenced by filter().

◆ vtxToken_

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

private

Definition at line 105 of file LeptonSkimming.h.

Referenced by filter().

◆ ZvertexTrg

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

◆ LeptonSkimming()

◆ ~LeptonSkimming()

Member Function Documentation

◆ beginStream()

◆ endStream()

◆ fillDescriptions()

◆ filter()

◆ hltFired()

◆ hltObject()

Member Data Documentation

◆ beam_x

◆ beam_y

◆ beam_z

◆ beamSpotToken_

◆ BiasedWP

◆ cleanedObjTracks

◆ cleanedPairTracks

◆ cleanedTracks

◆ conversionsToken_

◆ Cosee_Cut

◆ CoseeK_Cut

◆ EarlyStop

◆ el_eta

◆ el_phi

◆ el_pt

◆ eleBWPToken_

◆ electronsToken_

◆ eleUnBWPToken_

◆ ElTracks

◆ ElTrgExclusionCone

◆ ElTrgMuDz_Cut

◆ Epair_ObjectId

◆ Epair_ObjectIndex

◆ Epair_TrkIndex

◆ EpairZvtx_Cut

◆ EtaTrack_Cut

◆ HLTFilter_

◆ HLTPath_

◆ Ksdxy_Cut

◆ MaxChi2Track_Cut

◆ MaxMB_Cut

◆ MaxMee_Cut

◆ MinChi2Track_Cut

◆ MinMB_Cut

◆ MinMee_Cut

◆ muon_eta

◆ muon_medium

◆ muon_phi

◆ muon_pt

◆ muon_soft

◆ muon_tight

◆ muonsToken_

◆ MuTracks

◆ MuTrgExclusionCone

◆ MuTrgMatchCone

◆ MuTrgMuDz_Cut

◆ MuTrkMinDR_Cut

◆ nel

◆ nmuons

◆ ntracks

◆ object_container

◆ object_id

◆ ObjPtLargerThanTrack

◆ Probee_Cut

◆ ProbeeK_Cut

◆ PtB_Cut

◆ PtEl_Cut

◆ PtKTrack_Cut

◆ PtMu_Cut

◆ PtTrack_Cut

◆ QualMu_Cut

◆ Result

◆ SaveOnlyEPairTracks