Inheritance diagram for HSCPFilter:

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

	~HSCPFilter () override

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDFilter () 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)

Private Member Functions
void	beginJob () override

void	endJob () override

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

Private Attributes
double	chi2nMax

double	dedxMaxLeft

double	dedxMin

double	dxyMax

double	dzMax

double	etaMax

double	etaMin

bool	filterFlag

edm::EDGetTokenT< edm::ValueMap< reco::DeDxData > >	input_dedx_collectionToken

edm::EDGetTokenT< reco::MuonCollection >	input_muon_collectionToken

edm::EDGetTokenT< reco::TrackCollection >	input_track_collectionToken

int	ndedxHits

edm::EDGetTokenT< reco::VertexCollection >	recoVertexToken

double	SAMuPtMin

double	trkPtMin

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 51 of file Skim_HSCPFilter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 85 of file Skim_HSCPFilter.cc.

References chi2nMax, dedxMaxLeft, dedxMin, dxyMax, dzMax, etaMax, etaMin, filterFlag, edm::ParameterSet::getParameter(), input_dedx_collectionToken, input_muon_collectionToken, input_track_collectionToken, ndedxHits, recoVertexToken, SAMuPtMin, and trkPtMin.

 {
      filterFlag = iConfig.getParameter< bool >("filter");
 #ifdef THIS_IS_AN_EVENT_EXAMPLE
      pInToken = consumes<ExampleData>(iConfig.getParameter< edm::InputTag >("example"));
 #endif
      recoVertexToken = consumes<reco::VertexCollection>(edm::InputTag("offlinePrimaryVertices"));
      input_muon_collectionToken = consumes<reco::MuonCollection>(iConfig.getParameter< edm::InputTag >("inputMuonCollection"));
      input_track_collectionToken = consumes<reco::TrackCollection>(iConfig.getParameter< edm::InputTag >("inputTrackCollection"));
      input_dedx_collectionToken =  consumes<edm::ValueMap<reco::DeDxData> >(iConfig.getParameter< edm::InputTag >("inputDedxCollection"));
      dedxMin = iConfig.getParameter< double >("dedxMin");
      dedxMaxLeft = iConfig.getParameter< double >("dedxMaxLeft");
      trkPtMin = iConfig.getParameter< double >("trkPtMin");
      etaMin =  iConfig.getParameter< double >("etaMin");
      etaMax =  iConfig.getParameter< double >("etaMax");
      ndedxHits = iConfig.getParameter< int >("ndedxHits");
      chi2nMax = iConfig.getParameter< double >("chi2nMax");
      dxyMax = iConfig.getParameter< double >("dxyMax");
      dzMax = iConfig.getParameter< double >("dzMax");
      SAMuPtMin = iConfig.getParameter< double >("SAMuPtMin");
 }

HSCPFilter::~HSCPFilter ( )

override

Definition at line 108 of file Skim_HSCPFilter.cc.

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

Member Function Documentation

void HSCPFilter::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 212 of file Skim_HSCPFilter.cc.

213 {

214 }

void HSCPFilter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 218 of file Skim_HSCPFilter.cc.

References DEFINE_FWK_MODULE.

218 {

219 }

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

overrideprivate

Definition at line 123 of file Skim_HSCPFilter.cc.

References chi2nMax, dedxMaxLeft, dedxMin, SoftLeptonByDistance_cfi::distance, PVValHelper::dxy, dxyMax, PVValHelper::dz, dzMax, etaMax, etaMin, filterFlag, edm::EventSetup::get(), edm::Event::getByToken(), mps_fire::i, input_dedx_collectionToken, input_muon_collectionToken, input_track_collectionToken, edm::Ref< C, T, F >::isNull(), ndedxHits, position, edm::Handle< T >::product(), ecalDetailedTimeRecHit_cfi::recoVertex, recoVertexToken, SAMuPtMin, mathSSE::sqrt(), and trkPtMin.

 {
    using namespace edm;
 #ifdef THIS_IS_AN_EVENT_EXAMPLE
    Handle<ExampleData> pIn;
    iEvent.getByToken(pInToken,pIn);
 #endif
 
 #ifdef THIS_IS_AN_EVENTSETUP_EXAMPLE
    ESHandle<SetupData> pSetup;
    iSetup.get<SetupRecord>().get(pSetup);
 #endif
 
    using namespace reco;
 
    edm::Handle<reco::VertexCollection> recoVertexHandle;
    iEvent.getByToken(recoVertexToken, recoVertexHandle);
    reco::VertexCollection recoVertex = *recoVertexHandle;
 
    if(!filterFlag) return true;
 
    if(recoVertex.empty()) return false;
 
    using reco::MuonCollection;
 
    Handle<MuonCollection> muTracks;
    iEvent.getByToken(input_muon_collectionToken,muTracks);
    const reco::MuonCollection muonC = *(muTracks.product());
    for(unsigned int i=0; i<muonC.size(); i++){
       reco::MuonRef muon  = reco::MuonRef( muTracks, i );
       if(!muon->standAloneMuon().isNull()) {
          TrackRef SATrack = muon->standAloneMuon();
          if(SATrack->pt()>SAMuPtMin) return true;
       }
 
    }
 
 
 
 
 
 
    using reco::TrackCollection;
    Handle<TrackCollection> tkTracks;
    iEvent.getByToken(input_track_collectionToken,tkTracks);
    const reco::TrackCollection tkTC = *(tkTracks.product());
 
    Handle<ValueMap<DeDxData> >          dEdxTrackHandle;
    iEvent.getByToken(input_dedx_collectionToken, dEdxTrackHandle);
    const ValueMap<DeDxData> dEdxTrack = *dEdxTrackHandle.product();
 
    for(size_t i=0; i<tkTracks->size(); i++){
 
       reco::TrackRef trkRef = reco::TrackRef(tkTracks, i);
 
 
       if(trkRef->pt()>trkPtMin && trkRef->eta()<etaMax && trkRef->eta()>etaMin && trkRef->normalizedChi2()<chi2nMax){
 
            double dz  = trkRef->dz (recoVertex[0].position());
            double dxy = trkRef->dxy(recoVertex[0].position());
            double distancemin =sqrt(dxy*dxy+dz*dz);
            int closestvertex=0;
            for(unsigned int i=1;i<recoVertex.size();i++){
               dz  = trkRef->dz (recoVertex[i].position());
               dxy = trkRef->dxy(recoVertex[i].position());
               double distance = sqrt(dxy*dxy+dz*dz);
               if(distance < distancemin ){
                  distancemin = distance;
                  closestvertex=i;
               }
            }
 
            dz  = trkRef->dz (recoVertex[closestvertex].position());
            dxy = trkRef->dxy(recoVertex[closestvertex].position());
 
            if(fabs(dz)<dzMax && fabs(dxy)<dxyMax ){
 
              double dedx = dEdxTrack[trkRef].dEdx();
               int dedxnhits  = dEdxTrack[trkRef].numberOfMeasurements();
               if((dedx >dedxMin || dedx<dedxMaxLeft) && dedxnhits > ndedxHits) return true;
            }
         }
    }
    return false;
 
 }

Member Data Documentation

double HSCPFilter::chi2nMax

private

Definition at line 69 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

double HSCPFilter::dedxMaxLeft

private

Definition at line 69 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

double HSCPFilter::dedxMin

private

Definition at line 69 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

double HSCPFilter::dxyMax

private

Definition at line 69 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

double HSCPFilter::dzMax

private

Definition at line 69 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

double HSCPFilter::etaMax

private

Definition at line 69 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

double HSCPFilter::etaMin

private

Definition at line 69 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

bool HSCPFilter::filterFlag

private

Definition at line 60 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

edm::EDGetTokenT<edm::ValueMap<reco::DeDxData> > HSCPFilter::input_dedx_collectionToken

private

Definition at line 67 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

edm::EDGetTokenT<reco::MuonCollection> HSCPFilter::input_muon_collectionToken

private

Definition at line 65 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

edm::EDGetTokenT<reco::TrackCollection> HSCPFilter::input_track_collectionToken

private

Definition at line 66 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

int HSCPFilter::ndedxHits

private

Definition at line 68 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

edm::EDGetTokenT<reco::VertexCollection> HSCPFilter::recoVertexToken

private

Definition at line 64 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

double HSCPFilter::SAMuPtMin

private

Definition at line 69 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

double HSCPFilter::trkPtMin

private

Definition at line 69 of file Skim_HSCPFilter.cc.

Referenced by filter(), and HSCPFilter().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation