TotalKinematicsFilter Class Reference

#include <TotalKinematicsFilter.h>

Inheritance diagram for TotalKinematicsFilter:

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

Private Attributes
edm::InputTag	src_
double	tolerance_
bool	verbose_

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 19 of file TotalKinematicsFilter.h.

Constructor & Destructor Documentation

TotalKinematicsFilter::TotalKinematicsFilter ( const edm::ParameterSet &  iPSet )

explicit

Definition at line 5 of file TotalKinematicsFilter.cc.

                                                                        :  
  src_(iPSet.getParameter<edm::InputTag>("src")),
  tolerance_(iPSet.getParameter<double>("tolerance")),
  verbose_(iPSet.getUntrackedParameter<bool>("verbose",false))
{    
}

TotalKinematicsFilter::~TotalKinematicsFilter

(

)

Definition at line 12 of file TotalKinematicsFilter.cc.

{}

Member Function Documentation

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

privatevirtual

Implements edm::EDFilter.

Definition at line 14 of file TotalKinematicsFilter.cc.

References abs, gather_cfg::cout, genParticleCandidates2GenParticles_cfi::genParticles, edm::Event::getByLabel(), i, nPart(), mathSSE::sqrt(), src_, tolerance_, and verbose_.

{ 

  float nEcms = 0.;
  unsigned int nInit = 0;

  std::vector<float> p4tot(4,0.);
  unsigned int nPart = 0;

  // Gather information on the reco::GenParticle collection
  edm::Handle<reco::GenParticleCollection> genParticles;
  iEvent.getByLabel(src_, genParticles );
  
  for (reco::GenParticleCollection::const_iterator iter=genParticles->begin();iter!=genParticles->end();++iter){
    if ( nInit < 3 && (*iter).status() == 3 && (*iter).pdgId() == 2212 ) {
      nInit++;
      nEcms += (*iter).energy();
    }
    if ( (*iter).status() == 1) { 
      nPart++;
      if (verbose_) {
        std::cout << "Status 1 part # " << std::setw(4) << std::fixed << nPart 
                  << std::setw(14) << std::fixed << (*iter).pdgId() 
                  << std::setw(14) << std::fixed << (*iter).px() 
                  << std::setw(14) << std::fixed << (*iter).py() 
                  << std::setw(14) << std::fixed << (*iter).pz() << std::endl;
      }
      p4tot[0] += (*iter).px();
      p4tot[1] += (*iter).py();
      p4tot[2] += (*iter).pz();
      p4tot[3] += std::sqrt( (*iter).px()*(*iter).px() + 
                           (*iter).py()*(*iter).py() + 
                           (*iter).pz()*(*iter).pz() + 
                           (*iter).mass()*(*iter).mass()) ; 
    }
  }

  if ( verbose_ ) {
    std::cout << "Initial sqrt(s) = " << nEcms << std::endl;
    for (unsigned int i=0; i<4; i++) {
      std::cout << "p4tot["<<i<<"] = " << p4tot[i] << std::endl;
    }
  }

  bool pass = true;
  if ( std::abs(p4tot[0]) > tolerance_ || std::abs(p4tot[1]) > tolerance_ || std::abs(p4tot[2]) > tolerance_ || std::abs(p4tot[3]-nEcms) > tolerance_ ) 
    { pass = false; }

  return pass;
  
}

Member Data Documentation

edm::InputTag TotalKinematicsFilter::src_

private

Definition at line 30 of file TotalKinematicsFilter.h.

Referenced by filter().

double TotalKinematicsFilter::tolerance_

private

Definition at line 31 of file TotalKinematicsFilter.h.

Referenced by filter().

bool TotalKinematicsFilter::verbose_

private

Definition at line 32 of file TotalKinematicsFilter.h.

Referenced by filter().