#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	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

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

ModuleDescription const &	moduleDescription () const

	~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

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::vector< ModuleDescription const * > &modules, 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>
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)

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 86 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 edm::ParameterSetDescription::setUnknown().

◆ 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, runTheMatrix::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, SurveyInfoScenario_cff::seed, SelectedMu_DR, SelectedMu_index, SelectedTrgObj_PtEtaPhiCharge, SkimOnlyElectrons, SkimOnlyMuons, SkipIfNoMuMatch, SLxy_Cut, mathSSE::sqrt(), test_ev, TrackMuDz_Cut, PDWG_EXOHSCP_cff::tracks, Tracks_, TrackSdxy_Cut, TrgExclusionCone, trgresultsToken_, Trk_container, trk_index, TrkObjExclusionCone, UnbiasedWP, KalmanVertexFitter::vertex(), vertex_point, pwdgSkimBPark_cfi::vertices, badGlobalMuonTaggersAOD_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_2018_cff::InputTag, edm::HandleBase::isValid(), dqmiolumiharvest::j, relativeConstraints::keys, trigger::TriggerObject::phi(), trigger::TriggerObject::pt(), and trigobjectsToken_.

Referenced by filter().

Member Data Documentation

◆ beam_x

float LeptonSkimming::beam_x

private

Definition at line 114 of file LeptonSkimming.h.

Referenced by filter().

◆ beam_y

float LeptonSkimming::beam_y

private

Definition at line 114 of file LeptonSkimming.h.

Referenced by filter().

◆ beam_z

float LeptonSkimming::beam_z

private

Definition at line 114 of file LeptonSkimming.h.

Referenced by filter().

◆ beamSpotToken_

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

private

Definition at line 107 of file LeptonSkimming.h.

Referenced by filter().

◆ BiasedWP

double LeptonSkimming::BiasedWP = -100

private

Definition at line 175 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ cleanedObjTracks

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

private

Definition at line 128 of file LeptonSkimming.h.

Referenced by filter().

◆ cleanedPairTracks

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

private

Definition at line 129 of file LeptonSkimming.h.

Referenced by filter().

◆ cleanedTracks

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

private

Definition at line 132 of file LeptonSkimming.h.

Referenced by filter().

◆ conversionsToken_

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

private

Definition at line 108 of file LeptonSkimming.h.

Referenced by filter().

◆ Cosee_Cut

double LeptonSkimming::Cosee_Cut = -1

private

Definition at line 153 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ CoseeK_Cut

double LeptonSkimming::CoseeK_Cut = 0

private

Definition at line 160 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ EarlyStop

bool LeptonSkimming::EarlyStop = false

private

Definition at line 154 of file LeptonSkimming.h.

◆ el_eta

std::vector<float> LeptonSkimming::el_eta

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ el_phi

std::vector<float> LeptonSkimming::el_phi

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ el_pt

std::vector<float> LeptonSkimming::el_pt

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ eleBWPToken_

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

private

Definition at line 102 of file LeptonSkimming.h.

Referenced by filter().

◆ electronsToken_

edm::EDGetToken LeptonSkimming::electronsToken_

private

Definition at line 101 of file LeptonSkimming.h.

Referenced by filter().

◆ eleUnBWPToken_

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

private

Definition at line 103 of file LeptonSkimming.h.

Referenced by filter().

◆ ElTracks

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

private

Definition at line 131 of file LeptonSkimming.h.

Referenced by filter().

◆ ElTrgExclusionCone

double LeptonSkimming::ElTrgExclusionCone = 0

private

Definition at line 169 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ ElTrgMuDz_Cut

double LeptonSkimming::ElTrgMuDz_Cut = 1000

private

Definition at line 172 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ Epair_ObjectId

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

private

Definition at line 118 of file LeptonSkimming.h.

Referenced by filter().

◆ Epair_ObjectIndex

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

private

Definition at line 118 of file LeptonSkimming.h.

Referenced by filter().

◆ Epair_TrkIndex

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

private

Definition at line 118 of file LeptonSkimming.h.

Referenced by filter().

◆ EpairZvtx_Cut

double LeptonSkimming::EpairZvtx_Cut = 10000000

private

Definition at line 157 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ EtaTrack_Cut

double LeptonSkimming::EtaTrack_Cut = 10

private

Definition at line 139 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ HLTFilter_

std::vector<string> LeptonSkimming::HLTFilter_

private

Definition at line 111 of file LeptonSkimming.h.

Referenced by filter().

◆ HLTPath_

std::vector<string> LeptonSkimming::HLTPath_

private

Definition at line 111 of file LeptonSkimming.h.

Referenced by filter().

◆ Ksdxy_Cut

double LeptonSkimming::Ksdxy_Cut = 1000

private

Definition at line 158 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MaxChi2Track_Cut

double LeptonSkimming::MaxChi2Track_Cut = 1000

private

Definition at line 141 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MaxMB_Cut

double LeptonSkimming::MaxMB_Cut = 1000

private

Definition at line 145 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MaxMee_Cut

double LeptonSkimming::MaxMee_Cut = 1000

private

Definition at line 143 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MinChi2Track_Cut

double LeptonSkimming::MinChi2Track_Cut = -1000

private

Definition at line 140 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MinMB_Cut

double LeptonSkimming::MinMB_Cut = 0

private

Definition at line 146 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MinMee_Cut

double LeptonSkimming::MinMee_Cut = 0

private

Definition at line 151 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ muon_eta

std::vector<float> LeptonSkimming::muon_eta

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_medium

std::vector<float> LeptonSkimming::muon_medium

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_phi

std::vector<float> LeptonSkimming::muon_phi

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_pt

std::vector<float> LeptonSkimming::muon_pt

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_soft

std::vector<float> LeptonSkimming::muon_soft

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ muon_tight

std::vector<float> LeptonSkimming::muon_tight

private

Definition at line 117 of file LeptonSkimming.h.

Referenced by filter().

◆ muonsToken_

edm::EDGetToken LeptonSkimming::muonsToken_

private

Definition at line 104 of file LeptonSkimming.h.

Referenced by filter().

◆ MuTracks

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

private

Definition at line 130 of file LeptonSkimming.h.

Referenced by filter().

◆ MuTrgExclusionCone

double LeptonSkimming::MuTrgExclusionCone = 0

private

Definition at line 168 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MuTrgMatchCone

double LeptonSkimming::MuTrgMatchCone = 1000

private

Definition at line 155 of file LeptonSkimming.h.

Referenced by LeptonSkimming().

◆ MuTrgMuDz_Cut

double LeptonSkimming::MuTrgMuDz_Cut = 1000

private

Definition at line 171 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ MuTrkMinDR_Cut

double LeptonSkimming::MuTrkMinDR_Cut = 0

private

Definition at line 142 of file LeptonSkimming.h.

Referenced by LeptonSkimming().

◆ nel

unsigned int LeptonSkimming::nel = 0

private

Definition at line 136 of file LeptonSkimming.h.

Referenced by filter().

◆ nmuons

unsigned int LeptonSkimming::nmuons = 0

private

Definition at line 136 of file LeptonSkimming.h.

Referenced by filter().

◆ ntracks

unsigned int LeptonSkimming::ntracks = 0

private

Definition at line 136 of file LeptonSkimming.h.

Referenced by filter().

◆ object_container

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

private

Definition at line 134 of file LeptonSkimming.h.

Referenced by filter().

◆ object_id

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

private

Definition at line 134 of file LeptonSkimming.h.

Referenced by filter().

◆ ObjPtLargerThanTrack

bool LeptonSkimming::ObjPtLargerThanTrack = false

private

Definition at line 174 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ Probee_Cut

double LeptonSkimming::Probee_Cut = 0

private

Definition at line 152 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ ProbeeK_Cut

double LeptonSkimming::ProbeeK_Cut = 0

private

Definition at line 159 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtB_Cut

double LeptonSkimming::PtB_Cut = 0

private

Definition at line 164 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtEl_Cut

double LeptonSkimming::PtEl_Cut = 0

private

Definition at line 167 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtKTrack_Cut

double LeptonSkimming::PtKTrack_Cut = 0

private

Definition at line 144 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtMu_Cut

double LeptonSkimming::PtMu_Cut = 0

private

Definition at line 165 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ PtTrack_Cut

double LeptonSkimming::PtTrack_Cut = 0

private

Definition at line 138 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ QualMu_Cut

double LeptonSkimming::QualMu_Cut = 0

private

Definition at line 166 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ Result

bool LeptonSkimming::Result = false

private

Definition at line 113 of file LeptonSkimming.h.

Referenced by filter().

◆ SaveOnlyEPairTracks

bool LeptonSkimming::SaveOnlyEPairTracks = false

private

Definition at line 148 of file LeptonSkimming.h.

◆ SaveOnlyTracks

bool LeptonSkimming::SaveOnlyTracks = false

private

Definition at line 147 of file LeptonSkimming.h.

◆ SaveOutputRoot

bool LeptonSkimming::SaveOutputRoot = true

private

Definition at line 173 of file LeptonSkimming.h.

◆ SelectedMu_DR

float LeptonSkimming::SelectedMu_DR = -1

private

Definition at line 123 of file LeptonSkimming.h.

Referenced by filter().

◆ SelectedMu_index

int LeptonSkimming::SelectedMu_index

private

Definition at line 121 of file LeptonSkimming.h.

Referenced by filter().

◆ SelectedTrgObj_PtEtaPhiCharge

std::vector<float> LeptonSkimming::SelectedTrgObj_PtEtaPhiCharge

private

Definition at line 120 of file LeptonSkimming.h.

Referenced by filter().

◆ SkimOnlyElectrons

bool LeptonSkimming::SkimOnlyElectrons = false

private

Definition at line 178 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ SkimOnlyMuons

bool LeptonSkimming::SkimOnlyMuons = false

private

Definition at line 177 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ SkipIfNoMuMatch

bool LeptonSkimming::SkipIfNoMuMatch = false

private

Definition at line 156 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ SLxy_Cut

double LeptonSkimming::SLxy_Cut = 0

private

Definition at line 163 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ tempPtEtaPhiM

std::vector<float> LeptonSkimming::tempPtEtaPhiM

private

Definition at line 127 of file LeptonSkimming.h.

◆ tempXYZ

std::vector<float> LeptonSkimming::tempXYZ

private

Definition at line 127 of file LeptonSkimming.h.

◆ test_ev

int LeptonSkimming::test_ev = 0

private

Definition at line 179 of file LeptonSkimming.h.

Referenced by filter().

◆ track_vx

std::vector<float> LeptonSkimming::track_vx

private

Definition at line 116 of file LeptonSkimming.h.

◆ track_vy

std::vector<float> LeptonSkimming::track_vy

private

Definition at line 116 of file LeptonSkimming.h.

◆ track_vz

std::vector<float> LeptonSkimming::track_vz

private

Definition at line 116 of file LeptonSkimming.h.

◆ TrackMuDz_Cut

double LeptonSkimming::TrackMuDz_Cut = 100000000

private

Definition at line 161 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ trackObj_container

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

private

Definition at line 133 of file LeptonSkimming.h.

◆ trackPair_container

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

private

Definition at line 133 of file LeptonSkimming.h.

◆ Tracks_

edm::EDGetToken LeptonSkimming::Tracks_

private

Definition at line 105 of file LeptonSkimming.h.

Referenced by filter().

◆ TrackSdxy_Cut

double LeptonSkimming::TrackSdxy_Cut = 0

private

Definition at line 149 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ TrgExclusionCone

double LeptonSkimming::TrgExclusionCone = -1

private

Definition at line 162 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ trgresultsToken_

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

private

Definition at line 109 of file LeptonSkimming.h.

Referenced by filter(), and hltFired().

◆ trigobjectsToken_

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

private

Definition at line 110 of file LeptonSkimming.h.

Referenced by hltObject().

◆ Trk_container

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

private

Definition at line 134 of file LeptonSkimming.h.

Referenced by filter().

◆ trk_index

unsigned int LeptonSkimming::trk_index = 0

private

Definition at line 125 of file LeptonSkimming.h.

Referenced by filter().

◆ TrkObjExclusionCone

double LeptonSkimming::TrkObjExclusionCone = 0

private

Definition at line 170 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ UnbiasedWP

double LeptonSkimming::UnbiasedWP = -100

private

Definition at line 176 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

◆ UseOnlyBKeeMCForTriplets

bool LeptonSkimming::UseOnlyBKeeMCForTriplets = false

private

Definition at line 150 of file LeptonSkimming.h.

◆ vertex_point

reco::TrackBase::Point LeptonSkimming::vertex_point

private

Definition at line 126 of file LeptonSkimming.h.

Referenced by filter().

◆ vtxToken_

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

private

Definition at line 106 of file LeptonSkimming.h.

Referenced by filter().

◆ ZvertexTrg

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

private

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