Inheritance diagram for ChargedHadronTrackResolutionFilter:

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

	~ChargedHadronTrackResolutionFilter ()

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

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

ModuleDescription const &	moduleDescription () const

virtual	~EDFilterBase ()

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual bool	filter (edm::StreamID iID, edm::Event &, const edm::EventSetup &) const override

Private Attributes
const bool	debug_

const double	metSignifMin_

const double	p1_

const double	p2_

const double	p3_

const double	ptMin_

const bool	taggingMode_

edm::EDGetTokenT < reco::PFCandidateCollection >	tokenPFCandidates_

edm::EDGetTokenT < reco::PFMETCollection >	tokenPFMET_

Additional Inherited Members
Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Definition at line 27 of file ChargedHadronTrackResolutionFilter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 60 of file ChargedHadronTrackResolutionFilter.cc.

   : tokenPFCandidates_ ( consumes<reco::PFCandidateCollection>(iConfig.getParameter<edm::InputTag> ("PFCandidates")  ))
   , tokenPFMET_ (        consumes<reco::PFMETCollection>(iConfig.getParameter<edm::InputTag> ("PFMET")  )) 
   , taggingMode_          ( iConfig.getParameter<bool>    ("taggingMode")         )
   , ptMin_                ( iConfig.getParameter<double>        ("ptMin")         )
   , metSignifMin_         ( iConfig.getParameter<double>        ("MetSignifMin")  )
   , p1_                   ( iConfig.getParameter<double>        ("p1")            )
   , p2_                   ( iConfig.getParameter<double>        ("p2")            )
   , p3_                   ( iConfig.getParameter<double>        ("p3")            )
   , debug_                ( iConfig.getParameter<bool>          ("debug")         )
 {
   produces<bool>();
 }

ChargedHadronTrackResolutionFilter::~ChargedHadronTrackResolutionFilter ( )

Definition at line 74 of file ChargedHadronTrackResolutionFilter.cc.

74 { }

Member Function Documentation

bool ChargedHadronTrackResolutionFilter::filter	(	edm::StreamID	iID,
		edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)		const

overrideprivatevirtual

Implements edm::global::EDFilterBase.

Definition at line 83 of file ChargedHadronTrackResolutionFilter.cc.

References funct::abs(), gather_cfg::cout, debug_, edm::Event::getByToken(), i, edm::Ref< C, T, F >::isNull(), metSignifMin_, P, p1_, p2_, p3_, reco::LeafCandidate::pdgId(), reco::tau::pfCandidates(), pfMET_cfi::pfMET, funct::pow(), reco::tau::disc::Pt(), ptMin_, edm::Event::put(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), mathSSE::sqrt(), taggingMode_, tokenPFCandidates_, tokenPFMET_, and reco::PFCandidate::trackRef().

 {
   using namespace std;
   using namespace edm;
 
   Handle<reco::PFCandidateCollection>     pfCandidates;
   iEvent.getByToken(tokenPFCandidates_,pfCandidates);
   Handle<reco::PFMETCollection>     pfMET;
   iEvent.getByToken(tokenPFMET_,pfMET);
 
   bool foundBadTrack = false;
   if ( debug_ ) cout << "starting loop over pfCandidates" << endl;
 
   for ( unsigned i=0; i<pfCandidates->size(); ++i ) {
 
     const reco::PFCandidate & cand = (*pfCandidates)[i];
   
     if ( std::abs(cand.pdgId()) == 211 ) {
 
       if (cand.trackRef().isNull()) continue;
 
       const reco::TrackRef trackref = cand.trackRef();
       const double Pt = trackref->pt();
       const double DPt = trackref->ptError();
       if (Pt < ptMin_) continue;
       if ( debug_ ) cout << "charged hadron track pT > " << Pt << " GeV - " << " dPt: " << DPt << " GeV - algorithm: "  << trackref->algo() << std::endl;
 
       const double P = trackref->p();
       
       const unsigned int LostHits = trackref->numberOfLostHits();
 
       if ( (DPt/Pt) > (p1_ * sqrt(p2_*p2_/P+p3_*p3_) / (1.+LostHits)) ) {
         
         const double MET_px = pfMET->begin()->px();
         const double MET_py = pfMET->begin()->py();
         const double MET_et = pfMET->begin()->et();
         const double MET_sumEt = pfMET->begin()->sumEt();
         const double hadron_px = cand.px();
         const double hadron_py = cand.py();
         if (MET_sumEt == 0) continue;
         const double MET_signif = MET_et/MET_sumEt;
         const double MET_et_corr = sqrt( (MET_px + hadron_px)*(MET_px + hadron_px) + (MET_py + hadron_py)*(MET_py + hadron_py) );
         const double MET_signif_corr = MET_et_corr/MET_sumEt;
         if ( debug_ ) std::cout << "MET signif before: " << MET_signif << " - after: " << MET_signif_corr << " - reduction factor: " << MET_signif/MET_signif_corr << endl;
 
         if ( (MET_signif/MET_signif_corr) > metSignifMin_ ) {
 
           foundBadTrack = true;
 
           if ( debug_ ) {
             cout << cand << endl;
             cout << "charged hadron \t" << "track pT = " << Pt << " +/- " << DPt;
             cout << endl;
             cout << "MET: " << MET_et << " MET phi: " << pfMET->begin()->phi()<<
               " MET sumet: " << MET_sumEt << " MET significance: " << MET_et/MET_sumEt << endl;
             cout << "MET_px: " << MET_px << " MET_py: " << MET_py << " hadron_px: " << hadron_px << " hadron_py: " << hadron_py << endl;
             cout << "corrected: " << sqrt( pow((MET_px + hadron_px),2) + pow((MET_py + hadron_py),2)) << endl;
           }
         }
       }
     }
   } // end loop over PF candidates
 
   bool pass = !foundBadTrack;
 
   iEvent.put( std::auto_ptr<bool>(new bool(pass)) );
 
   return taggingMode_ || pass;
 }