TTbar_P4Violation Class Reference

#include <TTbar_P4Violation.h>

Inheritance diagram for TTbar_P4Violation:

Public Member Functions
	TTbar_P4Violation (const edm::ParameterSet &)
	~TTbar_P4Violation ()
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 ()

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)

Private Member Functions
virtual void	beginJob ()
virtual void	endJob ()
virtual bool	filter (edm::Event &, const edm::EventSetup &)

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 43 of file TTbar_P4Violation.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 of file TTbar_P4Violation.cc.

{
   //now do what ever initialization is needed

}

TTbar_P4Violation::~TTbar_P4Violation

(

)

Definition at line 21 of file TTbar_P4Violation.cc.

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

}

Member Function Documentation

void TTbar_P4Violation::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 109 of file TTbar_P4Violation.cc.

{
}

void TTbar_P4Violation::endJob ( void  )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 115 of file TTbar_P4Violation.cc.

                          {
}

void TTbar_P4Violation::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 120 of file TTbar_P4Violation.cc.

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

                                                                              {
  //The following says we do not know what parameters are allowed so do no validation
  // Please change this to state exactly what you do use, even if it is no parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

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

privatevirtual

Implements edm::EDFilter.

Definition at line 36 of file TTbar_P4Violation.cc.

References funct::false, TtGenEvtProducer_cfi::genEvt, edm::Event::getByLabel(), and funct::true.

{
   using namespace edm;

   // --- get TopQuarkAnalysis TtGenEvent
   Handle<TtGenEvent> genEvt;
   iEvent.getByLabel("genEvt", genEvt);

   if(!genEvt.isValid()) return false;

   const reco::GenParticle*        top = 0;
   const reco::GenParticle*    antitop = 0;
   const reco::GenParticle*     bottom = 0;
   const reco::GenParticle* antibottom = 0;
   const reco::GenParticle*      Wplus = 0;
   const reco::GenParticle*      Wmin  = 0;

   top        = genEvt->top();
   antitop    = genEvt->topBar();
   bottom     = genEvt->b();
   antibottom = genEvt->bBar();
   Wplus      = genEvt->wPlus();
   Wmin       = genEvt->wMinus();
   

   if(top && antitop && bottom && antibottom && Wplus && Wmin){
     const reco::Particle::LorentzVector     topP4 =     bottom->p4() + Wplus->p4() ;
     const reco::Particle::LorentzVector antitopP4 = antibottom->p4() + Wmin ->p4() ;
     
     double tolerance = 0.1 ;
     
     bool     topViolated = false ;
     bool antitopViolated = false ;
     
     if ( (top->p4().px() - topP4.px() > tolerance) || 
      (top->p4().py() - topP4.py() > tolerance) ||
      (top->p4().pz() - topP4.pz() > tolerance) ||
      (top->p4().e () - topP4.e () > tolerance) )  {
       
       topViolated = true ;
       
       //printf( "momentum not conserved for top:\n" ) ;
       //printf( " %5.5f\t %5.5f \t %5.5f \t %5.5f \n",      top->p4().px(),     top->p4().py(),     top->p4().pz(),     top->p4().e() ) ;
       //printf( " %5.5f\t %5.5f \t %5.5f \t %5.5f \n",          topP4.px(),         topP4.py(),         topP4.pz(),         topP4.e() ) ;
     }
     
     if ( (antitop->p4().px() - antitopP4.px() > tolerance) || 
      (antitop->p4().py() - antitopP4.py() > tolerance) ||
      (antitop->p4().pz() - antitopP4.pz() > tolerance) ||
      (antitop->p4().e () - antitopP4.e () > tolerance) )  {
       
       antitopViolated = true ;
       
       //printf( "momentum not conserved for anti-top:\n" ) ;
       //printf( " %5.5f\t %5.5f \t %5.5f \t %5.5f \n ", antitop->p4().px(), antitop->p4().py(), antitop->p4().pz(), antitop->p4().e() ) ;
       //printf( " %5.5f\t %5.5f \t %5.5f \t %5.5f \n ",     antitopP4.px(),     antitopP4.py(),     antitopP4.pz(),     antitopP4.e() ) ;
     }
     
     return (topViolated || antitopViolated);
     
     
     // GOSSIE temp
     bool     bottomOK = true ;
     bool antibottomOK = true ;
     if ( fabs(    bottom->p4().pz()) < 1. )     bottomOK = false ;
     if ( fabs(antibottom->p4().pz()) < 1. ) antibottomOK = false ;
     return (bottomOK && antibottomOK);
   }
   return false;
}