JetFlavourFilter Class Reference

#include <JetFlavourFilter.h>

Inheritance diagram for JetFlavourFilter:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)
	JetFlavourFilter (const edm::ParameterSet &)
	~JetFlavourFilter ()
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 Types
typedef std::vector < HepMC::GenParticle * >	GenPartVect
typedef std::vector < HepMC::GenParticle * > ::const_iterator	GenPartVectIt

Private Member Functions
HepMC::GenParticle *	findParticle (const GenPartVect genPartVect, const int requested_id)
void	printHisto (const HepMC::GenEvent::particle_iterator start, const HepMC::GenEvent::particle_iterator end)

Private Attributes
int	jetType
std::string	label_
int	noAccepted

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 25 of file JetFlavourFilter.h.

Member Typedef Documentation

typedef std::vector< HepMC::GenParticle * > JetFlavourFilter::GenPartVect

private

Definition at line 37 of file JetFlavourFilter.h.

typedef std::vector< HepMC::GenParticle * >::const_iterator JetFlavourFilter::GenPartVectIt

private

Definition at line 38 of file JetFlavourFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 9 of file JetFlavourFilter.cc.

References edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

{
  label_ = iConfig.getUntrackedParameter("moduleLabel",std::string("generator"));
  jetType = iConfig.getParameter< int >("jetType");
  if ((jetType>=1)&&(jetType<=3)) jetType=1;

  noAccepted = 0;
}

JetFlavourFilter::~JetFlavourFilter

(

)

Definition at line 19 of file JetFlavourFilter.cc.

References gather_cfg::cout.

{  
  std::cout << "Total number of accepted events = " << noAccepted << std::endl;
}

Member Function Documentation

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

virtual

Implements edm::EDFilter.

Definition at line 73 of file JetFlavourFilter.cc.

References alignCSCRings::ff, reco::flavour(), edm::Event::getByLabel(), i, AlCaHLTBitMon_ParallelJobs::p, parents, and detailsBasic3DVector::z.

{
  edm::Handle<HepMCProduct> evt;
  iEvent.getByLabel(label_, evt);
  
  const HepMC::GenEvent * generated_event = evt->GetEvent();

  bool event_passed = false;

  vector<int> foundQ;
  HepMC::GenEvent::particle_const_iterator p;
  for (p = generated_event->particles_begin(); p != generated_event->particles_end(); p++) {
    if (((*p)->pdg_id() < 1) && ((*p)->pdg_id() > -10)) { // We have an anti-quark
//       cout << "antiQ "<< (*p)->pdg_id();
      //vector< GenParticle * > parents = (*p)->listParents();
      vector< GenParticle * > parents;
      HepMC::GenVertex* inVertex = (*p)->production_vertex();
      for(GenVertex::particles_in_const_iterator iter = inVertex->particles_in_const_begin();
      iter != inVertex->particles_in_const_end();iter++)
    parents.push_back(*iter);
      
      for (GenPartVectIt z = parents.begin(); z != parents.end(); z++){

    vector< GenParticle * > child;
    HepMC::GenVertex* outVertex = (*z)->end_vertex();
    for(GenVertex::particles_in_const_iterator iter = outVertex->particles_in_const_begin();
        iter != outVertex->particles_in_const_end();iter++)
      child.push_back(*iter);
    
    if (findParticle(child, -(*p)->pdg_id())) foundQ.push_back(-(*p)->pdg_id());
      }
//       cout << " "<< foundQ.size()<<endl;
      
    }

  }
  
  int flavour = 0;
  int ff[6];
  ff[0]=0; ff[1]=0; ff[2]=0; ff[3]=0; ff[4]=0; ff[5]=0; 
  for (vector<int>::iterator i = foundQ.begin(); i != foundQ.end(); i++){
    ++ff[(*i)-1];
    if ((*i)>flavour) flavour = (*i);
  }
  // Is it light quark jet?
  if ((flavour>=0)&&(flavour<=3)) flavour=1;
  // Do we have more than one heavy flavour ?
  if ( (ff[3] && ff[4]) || (ff[3] && ff[5]) || (ff[4] && ff[5]) ) flavour =0;

  if (jetType==flavour) event_passed =true;

//   cout <<"Final flavour: " << ff[0]<<" "<<ff[1]<<" "<<ff[2]<<" "<<ff[3]<<" "<<ff[4]<<" "<<ff[5];
//   if (ff[0]||ff[1]||ff[2])  cout << " light";
//   if (ff[3])  cout << " charm";
//   if (ff[4])  cout << " bottom";
//   if (ff[5])  cout << " top";
//   if ((ff[0]||ff[1]||ff[2])&&((ff[3] + ff[4]+ff[5])))  {cout << " LH";};
//   if ( (ff[3] && ff[4]) || (ff[3] && ff[5]) || (ff[4] && ff[5]) ) {cout <<" notDef";}//printHisto(start,end);//}

//   cout<<" "<<flavour<<event_passed<<endl;

  if (event_passed) noAccepted++;

  delete generated_event; 

  return event_passed;
}

HepMC::GenParticle * JetFlavourFilter::findParticle ( const GenPartVect  genPartVect, const int  requested_id  )

private

Definition at line 25 of file JetFlavourFilter.cc.

References AlCaHLTBitMon_ParallelJobs::p.

{
  for (GenPartVectIt p = genPartVect.begin(); p != genPartVect.end(); p++)
    {
      if (requested_id == (*p)->pdg_id()) return *p;
    }
  return 0;
}

void JetFlavourFilter::printHisto ( const HepMC::GenEvent::particle_iterator  start, const HepMC::GenEvent::particle_iterator  end  )

private

Definition at line 38 of file JetFlavourFilter.cc.

References gather_cfg::cout, end, AlCaHLTBitMon_ParallelJobs::p, parents, and detailsBasic3DVector::z.

{
  HepMC::GenEvent::particle_iterator p;
  for (p = start; p != end; p++) 
    {
      //vector< GenParticle * > parents = (*p)->listParents();
      vector< GenParticle * > parents;
      HepMC::GenVertex* inVertex = (*p)->production_vertex();
      for(std::vector<GenParticle*>::const_iterator iter = inVertex->particles_in_const_begin();
      // for(GenVertex::particles_in_const_iterator iter = inVertex->particles_in_const_begin();
      iter != inVertex->particles_in_const_end();iter++)
    parents.push_back(*iter);
      
      cout << "isC "<<(*p)->pdg_id()<<" status "<<(*p)->status()<<" Parents: "<<parents.size()<<" - ";
      for (GenPartVectIt z = parents.begin(); z != parents.end(); z++){
    cout << (*z)->pdg_id()<<" ";
      }

      //vector< GenParticle * > child = (*p)->listChildren();
      vector< GenParticle * > child;
      HepMC::GenVertex* outVertex = (*p)->end_vertex();
      for(GenVertex::particles_in_const_iterator iter = outVertex->particles_in_const_begin();
      iter != outVertex->particles_in_const_end();iter++)
    child.push_back(*iter);
      
      cout << " - Child: "<<child.size()<<" - ";
      for (GenPartVectIt z = child.begin(); z != child.end(); z++){
    cout << (*z)->pdg_id()<<" ";
      }
      cout<<"\n";
    }
}

Member Data Documentation

int JetFlavourFilter::jetType

private

Definition at line 46 of file JetFlavourFilter.h.

std::string JetFlavourFilter::label_

private

Definition at line 48 of file JetFlavourFilter.h.

Referenced by Mixins._Labelable::_findDependencies(), and Mixins._TypedParameterizable::dumpPython().

int JetFlavourFilter::noAccepted

private

Definition at line 49 of file JetFlavourFilter.h.