#include <PFFilter.h>

Inheritance diagram for PFFilter:

Public Member Functions
virtual void	beginJob ()

bool	checkInput ()

virtual void	endJob ()

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

	PFFilter (const edm::ParameterSet &)

virtual	~PFFilter ()

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

virtual	~EDFilter ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Private Attributes
std::vector< std::string >	collections_

std::vector< int >	doMax_

std::vector< int >	doMin_

std::vector< double >	max_

std::vector< double >	min_

std::vector< std::string >	variables_

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

typedef WorkerT< EDFilter >	WorkerType

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 11 of file PFFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 11 of file PFFilter.cc.

References collections_, doMax_, doMin_, edm::ParameterSet::getParameter(), max_, min_, and variables_.

 {
   collections_= iConfig.getParameter< std::vector<std::string> >("Collections");
   variables_= iConfig.getParameter< std::vector<std::string> >("Variables");
   min_ = iConfig.getParameter< std::vector<double> >("Mins");
   max_ = iConfig.getParameter< std::vector<double> >("Maxs");
   doMin_ = iConfig.getParameter< std::vector<int> >("DoMin");
   doMax_ = iConfig.getParameter< std::vector<int> >("DoMax");
 }

PFFilter::~PFFilter ( )

virtual

Definition at line 21 of file PFFilter.cc.

21 {}

Member Function Documentation

void PFFilter::beginJob ( void )

virtual

Reimplemented from edm::EDFilter.

Definition at line 63 of file PFFilter.cc.

 {
   //std::cout << "FL: beginJob" << std::endl;
 }

bool PFFilter::checkInput ( )

Definition at line 23 of file PFFilter.cc.

References collections_, gather_cfg::cout, doMax_, doMin_, max_, min_, and variables_.

Referenced by filter().

 {
   if (collections_.size()!=min_.size())
   {
     std::cout << "Error: in PFFilter: collections_.size()!=min_.size()" << std::endl;
     std::cout << "collections_.size() = " << collections_.size() << std::endl;
     std::cout << "min_.size() = " << min_.size() << std::endl;
     return false;
   }
   if (collections_.size()!=max_.size())
   {
     std::cout << "Error: in PFFilter: collections_.size()!=max_.size()" << std::endl;
     std::cout << "collections_.size() = " << collections_.size() << std::endl;
     std::cout << "max_.size() = " << max_.size() << std::endl;
     return false;
   }
   if (collections_.size()!=doMin_.size())
   {
     std::cout << "Error: in PFFilter: collections_.size()!=min_.size()" << std::endl;
     std::cout << "collections_.size() = " << collections_.size() << std::endl;
     std::cout << "doMin_.size() = " << doMin_.size() << std::endl;
     return false;
   }
   if (collections_.size()!=doMax_.size())
   {
     std::cout << "Error: in PFFilter: collections_.size()!=min_.size()" << std::endl;
     std::cout << "collections_.size() = " << collections_.size() << std::endl;
     std::cout << "doMax_.size() = " << doMax_.size() << std::endl;
     return false;
   }
   if (collections_.size()!=variables_.size())
   {
     std::cout << "Error: in PFFilter: collections_.size()!=variables_.size()" << std::endl;
     std::cout << "collections_.size() = " << collections_.size() << std::endl;
     std::cout << "variables_.size() = " << variables_.size() << std::endl;
     return false;
   }
   return true;
 }

void PFFilter::endJob ( void )

virtual

Reimplemented from edm::EDFilter.

Definition at line 221 of file PFFilter.cc.

 {
   //std::cout << "FL: endJob" << std::endl;
 }

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

virtual

Implements edm::EDFilter.

Definition at line 68 of file PFFilter.cc.

References checkInput(), collections_, gather_cfg::cout, delta, doMax_, doMin_, reco::Candidate::et(), reco::Candidate::eta(), edm::Event::getByLabel(), M_PI, max_, min_, reco::Candidate::phi(), edm::Handle< T >::product(), findQualityFiles::size, createPayload::skip, and variables_.

 {
   //std::cout << "FL: filter" << std::endl;
   //std::cout << "FL: Mins = " << min_ << std::endl;
 
   if (!checkInput()) return true; // no filtering !
 
   bool skip=false;
 
   for (unsigned int varc=0;varc<collections_.size();++varc)
   {
     //std::cout << "FL: var[" << varc << "] = " << collections_[varc] << std::endl;
     //std::cout << "FL: var[0] = " << collections_[0] << std::endl;
 
     // if the collection is collection1-collection2:
     const unsigned int minuspos=collections_[varc].find("-");
     if (minuspos<collections_[varc].size())
     {
       std::string collection1;
       collection1.assign(collections_[varc],0,minuspos);
       //std::cout << "collection1 = " << collection1 << std::endl;
       std::string collection2;
       collection2.assign(collections_[varc],minuspos+1,collections_[varc].size());
       //std::cout << "collection2 = " << collection2 << std::endl;
 
       const edm::View<reco::Candidate> *var1;
       edm::Handle< edm::View<reco::Candidate> > var1_hnd;
       bool isVar1 = iEvent.getByLabel(collection1, var1_hnd);
       if ( !isVar1 ) { 
     std::cout << "Warning : no " << collection1 << " in input !" << std::endl;
     return false;
       }
       var1 = var1_hnd.product();
       const reco::Candidate *var10 = &(*var1)[0];
       //std::cout << "FL: var10.pt = " << var10->et() << std::endl;
       double coll_var1;
       if (variables_[varc]=="et") coll_var1=var10->et();
       else if (variables_[varc]=="phi") coll_var1=var10->phi();
       else if (variables_[varc]=="eta") coll_var1=var10->eta();
       else
       {
     std::cout << "Error: PFFilter: variable unknown: " << variables_[varc] << std::endl;
     return true;
       }
       //std::cout << "FL: coll_var1[" << variables_[varc] << "] = " << coll_var1 << std::endl;
 
       const edm::View<reco::Candidate> *var2;
       edm::Handle< edm::View<reco::Candidate> > var2_hnd;
       bool isVar2 = iEvent.getByLabel(collection2, var2_hnd);
       if ( !isVar2 ) { 
     std::cout << "Warning : no " << collection2 << " in input !" << std::endl;
     return false;
       }
       var2 = var2_hnd.product();
       const reco::Candidate *var20 = &(*var2)[0];
       //std::cout << "FL: var20.pt = " << var20->et() << std::endl;
       double coll_var2;
       if (variables_[varc]=="et") coll_var2=var20->et();
       else if (variables_[varc]=="phi") coll_var2=var20->phi();
       else if (variables_[varc]=="eta") coll_var2=var20->eta();
       else
       {
     std::cout << "Error: PFFilter: variable unknown: " << variables_[varc] << std::endl;
     return true;
       }
       //std::cout << "FL: coll_var2[" << variables_[varc] << "] = " << coll_var2 << std::endl;
       //std::cout << "FL: max_[varc] = " << max_[varc] << std::endl;
       //std::cout << "FL: min_[varc] = " << min_[varc] << std::endl;
 
       // Delta computation
       double delta=coll_var1-coll_var2;
       if (variables_[varc]=="phi")
       {
     if (coll_var1 > M_PI) coll_var1 -= ceil((coll_var1 - M_PI) / (2 * M_PI)) * 2 * M_PI;
     if (coll_var1 <= - M_PI) coll_var1 += ceil((coll_var1 + M_PI) / (-2. * M_PI)) * 2. * M_PI;
     if (coll_var2 > M_PI) coll_var2 -= ceil((coll_var2 - M_PI) / (2 * M_PI)) * 2 * M_PI;
     if (coll_var2 <= - M_PI) coll_var2 += ceil((coll_var2 + M_PI) / (-2. * M_PI)) * 2 * M_PI;
 
     double deltaphi=-999.0;
     if (fabs(coll_var1 - coll_var2)<M_PI)
     {
       deltaphi=(coll_var1-coll_var2);
     }
     else
     {
       if ((coll_var1-coll_var2)>0.0)
       {
         deltaphi=(2*M_PI - fabs(coll_var1 - coll_var2));
       }
       else
       {
         deltaphi=-(2*M_PI - fabs(coll_var1 - coll_var2));
       }
     }
     delta=deltaphi;
       }
 
       // cuts
       if (doMin_[varc] && doMax_[varc] && max_[varc]<min_[varc])
       {
     if (delta>max_[varc] && delta<min_[varc]) skip=true;
       }
       else
       {
     if (doMin_[varc] && delta<min_[varc]) skip=true;
     if (doMax_[varc] && delta>max_[varc]) skip=true;
       }
       //std::cout << "skip = " << skip << std::endl;
     }
     else
     {
       // get the variable:
       const edm::View<reco::Candidate> *var0;
       edm::Handle< edm::View<reco::Candidate> > var0_hnd;
       bool isVar0 = iEvent.getByLabel(collections_[varc], var0_hnd);
       if ( !isVar0 ) { 
     std::cout << "Warning : no " << collections_[varc] << " in input !" << std::endl;
     return false;
       }
       var0 = var0_hnd.product();
       const reco::Candidate *var00 = &(*var0)[0];
       //std::cout << "FL: var00.pt = " << var00->et() << std::endl;
       double coll_var;
       if (variables_[varc]=="et") coll_var=var00->et();
       else if (variables_[varc]=="phi") coll_var=var00->phi();
       else if (variables_[varc]=="eta") coll_var=var00->eta();
       else
       {
     std::cout << "Error: PFFilter: variable unknown: " << variables_[varc] << std::endl;
     return true;
       }
       //std::cout << "FL: coll_var[" << variables_[varc] << "] = " << coll_var << std::endl;
       //std::cout << "FL: max_[varc] = " << max_[varc] << std::endl;
       //std::cout << "FL: min_[varc] = " << min_[varc] << std::endl;
 
       // cuts
       if (doMin_[varc] && doMax_[varc] && max_[varc]<min_[varc])
       {
     if (coll_var>max_[varc] && coll_var<min_[varc]) skip=true;
       }
       else
       {
     if (doMin_[varc] && coll_var<min_[varc]) skip=true;
     if (doMax_[varc] && coll_var>max_[varc]) skip=true;
       }
       //std::cout << "skip = " << skip << std::endl;
     }
   }
   //std::cout << "final skip = " << skip << std::endl;
   return !skip;
 }