JGJFilter Class Reference

#include <JGJFilter.h>

Inheritance diagram for JGJFilter:

Public Member Functions
	JGJFilter (const edm::ParameterSet &)
	~JGJFilter ()
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 Member Functions
virtual void	beginJob ()
virtual void	endJob ()
virtual bool	filter (edm::Event &, const edm::EventSetup &)

Private Attributes
unsigned int	nAccepted
unsigned int	nEvents

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 32 of file JGJFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 21 of file JGJFilter.cc.

References nEvents.

{
   //Number of events and Number Accepted
   nEvents = 0;
   nAccepted = 0;
}

JGJFilter::~JGJFilter

(

)

Definition at line 29 of file JGJFilter.cc.

References gather_cfg::cout, and nEvents.

{
std::cout << "Total number of tested events = " << nEvents << std::endl;
std::cout << "Total number of accepted events = " << nAccepted << std::endl;
std::cout << "nEfficiency = " << ((double)nAccepted)/((double)nEvents) << std::endl; 
}

Member Function Documentation

void JGJFilter::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 143 of file JGJFilter.cc.

{
}

void JGJFilter::endJob ( void  )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 149 of file JGJFilter.cc.

{
}

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

privatevirtual

Implements edm::EDFilter.

Definition at line 36 of file JGJFilter.cc.

References edm::Event::getByLabel(), metsig::jet, nEvents, and dt_dqm_sourceclient_common_cff::reco.

{
   nEvents++;
   bool accepted = false;

   using namespace edm;
   using namespace reco;
   using namespace std;

   //KT4 jets ***************************************************************************
   Handle<reco::GenJetCollection> o4;
   iEvent.getByLabel("kt4GenJetsjgj",o4);

   if( o4.isValid()) {
 
     float jet_pt[3];
     float jet_eta[3];

     int njets(0);
     for( GenJetCollection::const_iterator jet = o4->begin(); jet != o4->end(); ++jet ) {
       if(njets<4) {
        jet_pt[njets]=jet->pt();
        jet_eta[njets]=jet->eta();
       }
       njets++;
     }

     if(njets>1 && fabs(jet_eta[0]-jet_eta[1])>3.5) accepted = true;

     if(njets>2) {
       if((jet_pt[2]/jet_pt[1])>0.8) {
       if(fabs(jet_eta[0]-jet_eta[2])>3.5) accepted = true;
       if(fabs(jet_eta[1]-jet_eta[2])>3.5) accepted = true;
       }
     }
   } //valid

  
   //KT6 jets ***************************************************************************
   Handle<reco::GenJetCollection> o6;
   iEvent.getByLabel("kt6GenJetsjgj",o6);

   if( o6.isValid()) {
 
     float jet_pt[3];
     float jet_eta[3];

     int njets(0);
     for( GenJetCollection::const_iterator jet = o6->begin(); jet != o6->end(); ++jet ) {
       if(njets<4) {
        jet_pt[njets]=jet->pt();
        jet_eta[njets]=jet->eta();
       }
       njets++;
     }

     if(njets>1 && fabs(jet_eta[0]-jet_eta[1])>3.5) accepted = true;

     if(njets>2) {
       if((jet_pt[2]/jet_pt[1])>0.8) {
       if(fabs(jet_eta[0]-jet_eta[2])>3.5) accepted = true;
       if(fabs(jet_eta[1]-jet_eta[2])>3.5) accepted = true;
       }
     }
   } //valid
  

  //SC5 jets ***************************************************************************
   Handle<reco::GenJetCollection> oo;
   iEvent.getByLabel("sisCone5GenJetsjgj",oo);

   if( oo.isValid()) {
 
     float jet_pt[3];
     float jet_eta[3];

     int njets(0);
     for( GenJetCollection::const_iterator jet = oo->begin(); jet != oo->end(); ++jet ) {
       if(njets<4) {
        jet_pt[njets]=jet->pt();
        jet_eta[njets]=jet->eta();
       }
       njets++;
     }

     if(njets>1 && fabs(jet_eta[0]-jet_eta[1])>3.5) accepted = true;



     if(njets>2) {
       if((jet_pt[2]/jet_pt[1])>0.8) {
       if(fabs(jet_eta[0]-jet_eta[2])>3.5) accepted = true;
       if(fabs(jet_eta[1]-jet_eta[2])>3.5) accepted = true;
       }
     }
   } //valid


   if ( accepted ){
    nAccepted++;
        return true;
   }
   else return false;

}

Member Data Documentation

unsigned int JGJFilter::nAccepted

private

Definition at line 46 of file JGJFilter.h.

unsigned int JGJFilter::nEvents

private

Definition at line 45 of file JGJFilter.h.