LeptonSkimming Class Reference

#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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
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 const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 87 of file LeptonSkimming.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 28 of file LeptonSkimming.cc.

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

                                                            :
  electronsToken_(consumes<std::vector<reco::GsfElectron> >(iConfig.getParameter<edm::InputTag>  ("electrons"))),
   eleBWPToken_(consumes<edm::ValueMap<float> >(iConfig.getParameter<edm::InputTag>("eleBiasedWP"))),
  eleUnBWPToken_(consumes<edm::ValueMap<float> >(iConfig.getParameter<edm::InputTag>("eleUnbiasedWP"))),
  muonsToken_(consumes<std::vector<reco::Muon>>(iConfig.getParameter<edm::InputTag>("muons"))),
  Tracks_(consumes<std::vector<reco::Track> >(iConfig.getParameter<edm::InputTag>("tracks"))),
  vtxToken_(consumes<std::vector<reco::Vertex>>(iConfig.getParameter<edm::InputTag>("vertices"))),
  beamSpotToken_(consumes<reco::BeamSpot>(iConfig.getParameter<edm::InputTag>("beamSpot"))),
  conversionsToken_(consumes< reco::ConversionCollection > (iConfig.getParameter<edm::InputTag> ("conversions"))),

   trgresultsToken_(consumes<edm::TriggerResults>(iConfig.getParameter<edm::InputTag> ("triggerresults"))),
  trigobjectsToken_(consumes<trigger::TriggerEvent>(iConfig.getParameter<edm::InputTag> ("triggerobjects"))),
  HLTFilter_(iConfig.getParameter<vector<string> >("HLTFilter")),
  HLTPath_(iConfig.getParameter<vector<string> >("HLTPath"))
  {
  edm::ParameterSet runParameters=iConfig.getParameter<edm::ParameterSet>("RunParameters");     
 //mu matching with trigger
 MuTrgMatchCone=runParameters.getParameter<double>("MuTrgMatchCone");
 SkipIfNoMuMatch=runParameters.getParameter<bool>("SkipIfNoMuMatch");
 //track cuts
 PtTrack_Cut=runParameters.getParameter<double>("PtTrack_Cut");
 EtaTrack_Cut=runParameters.getParameter<double>("EtaTrack_Cut");
 MinChi2Track_Cut=runParameters.getParameter<double>("MinChi2Track_Cut");
 MaxChi2Track_Cut=runParameters.getParameter<double>("MaxChi2Track_Cut");
 MuTrkMinDR_Cut=runParameters.getParameter<double>("MuTrkMinDR_Cut");
 TrackMuDz_Cut=runParameters.getParameter<double>("TrackMuDz_Cut");
 TrgExclusionCone=runParameters.getParameter<double>("TrgExclusionCone");
 //lepton cuts
 PtMu_Cut=runParameters.getParameter<double>("PtMu_Cut");
 PtEl_Cut=runParameters.getParameter<double>("PtEl_Cut");
 QualMu_Cut=runParameters.getParameter<double>("QualMu_Cut");
 MuTrgExclusionCone=runParameters.getParameter<double>("MuTrgExclusionCone");
 ElTrgExclusionCone=runParameters.getParameter<double>("ElTrgExclusionCone");
 TrkObjExclusionCone=runParameters.getParameter<double>("TrkObjExclusionCone");
 MuTrgMuDz_Cut=runParameters.getParameter<double>("MuTrgMuDz_Cut");
 ElTrgMuDz_Cut=runParameters.getParameter<double>("ElTrgMuDz_Cut");
 BiasedWP=runParameters.getParameter<double>("BiasedWP");
 UnbiasedWP=runParameters.getParameter<double>("UnbiasedWP");
 SkimOnlyMuons=runParameters.getParameter<bool>("SkimOnlyMuons");
 SkimOnlyElectrons=runParameters.getParameter<bool>("SkimOnlyElectrons");
 //pair cuts
 MaxMee_Cut=runParameters.getParameter<double>("MaxMee_Cut");
 MinMee_Cut=runParameters.getParameter<double>("MinMee_Cut");
 Probee_Cut=runParameters.getParameter<double>("Probee_Cut");
 Cosee_Cut=runParameters.getParameter<double>("Cosee_Cut");
 EpairZvtx_Cut=runParameters.getParameter<double>("EpairZvtx_Cut");
 //kaon
 PtKTrack_Cut=runParameters.getParameter<double>("PtKTrack_Cut");
 TrackSdxy_Cut=runParameters.getParameter<double>("TrackSdxy_Cut");
 Ksdxy_Cut=runParameters.getParameter<double>("Ksdxy_Cut");
 //triplet
 MaxMB_Cut=runParameters.getParameter<double>("MaxMB_Cut");
 MinMB_Cut=runParameters.getParameter<double>("MinMB_Cut");  
 ProbeeK_Cut=runParameters.getParameter<double>("ProbeeK_Cut");
 CoseeK_Cut=runParameters.getParameter<double>("CoseeK_Cut"); 
 SLxy_Cut=runParameters.getParameter<double>("SLxy_Cut");
 PtB_Cut=runParameters.getParameter<double>("PtB_Cut");
 
 ObjPtLargerThanTrack=runParameters.getParameter<bool>("ObjPtLargerThanTrack");

}

LeptonSkimming::~LeptonSkimming ( )

override

Definition at line 91 of file LeptonSkimming.cc.

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

}

Member Function Documentation

void LeptonSkimming::beginStream ( edm::StreamID  )

overrideprivatevirtual

Reimplemented from edm::stream::EDFilterBase.

Definition at line 436 of file LeptonSkimming.cc.

{
}

void LeptonSkimming::endStream ( )

overrideprivatevirtual

Reimplemented from edm::stream::EDFilterBase.

Definition at line 442 of file LeptonSkimming.cc.

                          {
}

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

static

Definition at line 479 of file LeptonSkimming.cc.

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

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

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

overrideprivate

Definition at line 170 of file LeptonSkimming.cc.

References beam_x, beam_y, beam_z, beamSpotToken_, BiasedWP, ChiSquaredProbability(), cleanedObjTracks, cleanedPairTracks, cleanedTracks, conversions_cfi::conversions, conversionsToken_, Cosee_Cut, CoseeK_Cut, TransientVertex::degreesOfFreedom(), boostedElectronIsolation_cff::deltaR, reco::BeamSpot::dxdz(), reco::BeamSpot::dydz(), el_eta, el_phi, el_pt, eleBWPToken_, nano_cff::electrons, electronsToken_, eleUnBWPToken_, ElTracks, ElTrgExclusionCone, ElTrgMuDz_Cut, Epair_ObjectId, Epair_ObjectIndex, Epair_TrkIndex, EpairZvtx_Cut, reco::TrackBase::eta(), EtaTrack_Cut, edm::EventSetup::get(), edm::Event::getByToken(), ConversionTools::hasMatchedConversion(), muons_cff::highPurity, HLTFilter_, hltFired(), hltObject(), HLTPath_, isMediumMuonCustom(), edm::Ref< C, T, F >::isNull(), muon::isSoftMuon(), isTightMuonCustom(), TransientVertex::isValid(), Ksdxy_Cut, SiStripPI::max, MaxChi2Track_Cut, MaxMB_Cut, MaxMee_Cut, MinChi2Track_Cut, MinMB_Cut, MinMee_Cut, RPCpg::mu, muon_eta, muon_medium, muon_phi, muon_pt, muon_soft, muon_tight, extraflags_cff::muons, muonsToken_, MuTracks, MuTrgExclusionCone, MuTrgMuDz_Cut, nel, nmuons, ntracks, hgcalPlots::obj, object_container, object_id, ObjPtLargerThanTrack, PV3DBase< T, PVType, FrameType >::perp(), 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, GlobalErrorBase< T, ErrorWeightType >::rerr(), Result, SurveyInfoScenario_cff::seed, SelectedMu_DR, SelectedMu_index, SelectedTrgObj_PtEtaPhiCharge, SkimOnlyElectrons, SkimOnlyMuons, SkipIfNoMuMatch, SLxy_Cut, mathSSE::sqrt(), test_ev, TransientVertex::totalChiSquared(), TrackMuDz_Cut, l1t::tracks, Tracks_, TrackSdxy_Cut, TrgExclusionCone, trgresultsToken_, Trk_container, trk_index, TrkObjExclusionCone, UnbiasedWP, KalmanVertexFitter::vertex(), vertex_point, electrons_cff::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.

Referenced by Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::filter(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::find(), Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::setDataAccessor(), and Vispa.Plugins.Browser.BrowserTabController.BrowserTabController::switchCenterView().

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

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

private

Definition at line 100 of file LeptonSkimming.cc.

References edm::HLTGlobalStatus::accept(), edm::HandleBase::failedToGet(), edm::Event::getByToken(), edm::HLTGlobalStatus::size(), AlCaHLTBitMon_QueryRunRegistry::string, trgresultsToken_, edm::TriggerNames::triggerName(), and edm::Event::triggerNames().

Referenced by filter().

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

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

private

Definition at line 120 of file LeptonSkimming.cc.

References a, b, trigger::TriggerObject::eta(), edm::Event::getByToken(), trigger::TriggerEvent::getObjects(), trigger::TriggerObject::id(), edm::HandleBase::isValid(), relativeConstraints::keys, trigger::TriggerObject::phi(), trigger::TriggerObject::pt(), and trigobjectsToken_.

Referenced by filter().

                                                                                                                      {
  using namespace std;  using namespace edm;  using namespace reco;
  using namespace trigger;
 
  edm::Handle<trigger::TriggerEvent> triggerObjectsSummary;
  iEvent.getByToken(trigobjectsToken_ ,triggerObjectsSummary);
  trigger::TriggerObjectCollection selectedObjects;

  std::vector<std::array<float,5>> max_per_trigger; 

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

Member Data Documentation

float LeptonSkimming::beam_x

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

float LeptonSkimming::beam_y

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

float LeptonSkimming::beam_z

private

Definition at line 116 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 109 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::BiasedWP =-100

private

Definition at line 154 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

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

private

Definition at line 128 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 129 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 132 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 110 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::Cosee_Cut =-1

private

Definition at line 144 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::CoseeK_Cut =0

private

Definition at line 147 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::EarlyStop =false

private

Definition at line 145 of file LeptonSkimming.h.

std::vector<float> LeptonSkimming::el_eta

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::el_phi

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::el_pt

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 104 of file LeptonSkimming.h.

Referenced by filter().

edm::EDGetToken LeptonSkimming::electronsToken_

private

Definition at line 103 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 105 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 131 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::ElTrgExclusionCone =0

private

Definition at line 150 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::ElTrgMuDz_Cut =1000

private

Definition at line 152 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

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

private

Definition at line 120 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 120 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 120 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::EpairZvtx_Cut =10000000

private

Definition at line 146 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::EtaTrack_Cut =10

private

Definition at line 138 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

std::vector<string> LeptonSkimming::HLTFilter_

private

Definition at line 113 of file LeptonSkimming.h.

Referenced by filter().

std::vector<string> LeptonSkimming::HLTPath_

private

Definition at line 113 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::Ksdxy_Cut =1000

private

Definition at line 146 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MaxChi2Track_Cut =1000

private

Definition at line 139 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MaxMB_Cut =1000

private

Definition at line 141 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MaxMee_Cut =1000

private

Definition at line 139 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MinChi2Track_Cut =-1000

private

Definition at line 138 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MinMB_Cut =0

private

Definition at line 141 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MinMee_Cut =0

private

Definition at line 144 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

std::vector<float> LeptonSkimming::muon_eta

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_medium

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_phi

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_pt

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_soft

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::muon_tight

private

Definition at line 119 of file LeptonSkimming.h.

Referenced by filter().

edm::EDGetToken LeptonSkimming::muonsToken_

private

Definition at line 106 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 130 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::MuTrgExclusionCone =0

private

Definition at line 150 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MuTrgMatchCone =1000

private

Definition at line 145 of file LeptonSkimming.h.

Referenced by LeptonSkimming().

double LeptonSkimming::MuTrgMuDz_Cut =1000

private

Definition at line 151 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::MuTrkMinDR_Cut =0

private

Definition at line 139 of file LeptonSkimming.h.

Referenced by LeptonSkimming().

unsigned int LeptonSkimming::nel =0

private

Definition at line 136 of file LeptonSkimming.h.

Referenced by filter().

unsigned int LeptonSkimming::nmuons =0

private

Definition at line 136 of file LeptonSkimming.h.

Referenced by filter().

unsigned int LeptonSkimming::ntracks =0

private

Definition at line 136 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 134 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 134 of file LeptonSkimming.h.

Referenced by filter().

bool LeptonSkimming::ObjPtLargerThanTrack =false

private

Definition at line 153 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::Probee_Cut =0

private

Definition at line 144 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::ProbeeK_Cut =0

private

Definition at line 146 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtB_Cut =0

private

Definition at line 149 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtEl_Cut =0

private

Definition at line 149 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtKTrack_Cut =0

private

Definition at line 140 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtMu_Cut =0

private

Definition at line 149 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::PtTrack_Cut =0

private

Definition at line 138 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::QualMu_Cut =0

private

Definition at line 149 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::Result =false

private

Definition at line 115 of file LeptonSkimming.h.

Referenced by filter().

bool LeptonSkimming::SaveOnlyEPairTracks =false

private

Definition at line 142 of file LeptonSkimming.h.

bool LeptonSkimming::SaveOnlyTracks =false

private

Definition at line 142 of file LeptonSkimming.h.

bool LeptonSkimming::SaveOutputRoot =true

private

Definition at line 152 of file LeptonSkimming.h.

float LeptonSkimming::SelectedMu_DR =-1

private

Definition at line 124 of file LeptonSkimming.h.

Referenced by filter().

int LeptonSkimming::SelectedMu_index

private

Definition at line 122 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::SelectedTrgObj_PtEtaPhiCharge

private

Definition at line 122 of file LeptonSkimming.h.

Referenced by filter().

bool LeptonSkimming::SkimOnlyElectrons =false

private

Definition at line 155 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::SkimOnlyMuons =false

private

Definition at line 155 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::SkipIfNoMuMatch =false

private

Definition at line 145 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::SLxy_Cut =0

private

Definition at line 148 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

std::vector<float> LeptonSkimming::tempPtEtaPhiM

private

Definition at line 127 of file LeptonSkimming.h.

std::vector<float> LeptonSkimming::tempXYZ

private

Definition at line 127 of file LeptonSkimming.h.

int LeptonSkimming::test_ev =0

private

Definition at line 156 of file LeptonSkimming.h.

Referenced by filter().

std::vector<float> LeptonSkimming::track_vx

private

Definition at line 118 of file LeptonSkimming.h.

std::vector<float> LeptonSkimming::track_vy

private

Definition at line 118 of file LeptonSkimming.h.

std::vector<float> LeptonSkimming::track_vz

private

Definition at line 118 of file LeptonSkimming.h.

double LeptonSkimming::TrackMuDz_Cut =100000000

private

Definition at line 147 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

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

private

Definition at line 133 of file LeptonSkimming.h.

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

private

Definition at line 133 of file LeptonSkimming.h.

edm::EDGetToken LeptonSkimming::Tracks_

private

Definition at line 107 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::TrackSdxy_Cut =0

private

Definition at line 143 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::TrgExclusionCone =-1

private

Definition at line 148 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

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

private

Definition at line 111 of file LeptonSkimming.h.

Referenced by filter(), and hltFired().

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

private

Definition at line 112 of file LeptonSkimming.h.

Referenced by hltObject().

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

private

Definition at line 134 of file LeptonSkimming.h.

Referenced by filter().

unsigned int LeptonSkimming::trk_index =0

private

Definition at line 125 of file LeptonSkimming.h.

Referenced by filter().

double LeptonSkimming::TrkObjExclusionCone =0

private

Definition at line 151 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

double LeptonSkimming::UnbiasedWP =-100

private

Definition at line 154 of file LeptonSkimming.h.

Referenced by filter(), and LeptonSkimming().

bool LeptonSkimming::UseOnlyBKeeMCForTriplets =false

private

Definition at line 143 of file LeptonSkimming.h.

reco::TrackBase::Point LeptonSkimming::vertex_point

private

Definition at line 126 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 108 of file LeptonSkimming.h.

Referenced by filter().

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

private

Definition at line 125 of file LeptonSkimming.h.

Referenced by filter().