#include <JetFlavourCutFilter.h>

Inheritance diagram for JetFlavourCutFilter:

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

	JetFlavourCutFilter (const edm::ParameterSet &)

	~JetFlavourCutFilter ()

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

ModuleDescription const &	moduleDescription () const

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

int	noAccepted

edm::EDGetTokenT < edm::HepMCProduct >	token_

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

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

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 25 of file JetFlavourCutFilter.h.

Member Typedef Documentation

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

private

Definition at line 37 of file JetFlavourCutFilter.h.

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

private

Definition at line 38 of file JetFlavourCutFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 9 of file JetFlavourCutFilter.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_25ns14e33_v1_cff::InputTag, HLT_25ns14e33_v1_cff::jetType, and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   token_ = consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared"));
   jetType = iConfig.getParameter< int >("jetType");
   if ((jetType>=1)&&(jetType<=3)) jetType=1;
 
   noAccepted = 0;
 }

JetFlavourCutFilter::~JetFlavourCutFilter ( )

Definition at line 19 of file JetFlavourCutFilter.cc.

References gather_cfg::cout.

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

Member Function Documentation

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

virtual

Implements edm::EDFilter.

Definition at line 72 of file JetFlavourCutFilter.cc.

References class-composition::child, alignCSCRings::ff, reco::flavour(), edm::Event::getByToken(), i, HLT_25ns14e33_v1_cff::jetType, AlCaHLTBitMon_ParallelJobs::p, and parents.

 {
   edm::Handle<HepMCProduct> evt;
   iEvent.getByToken(token_, 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(std::vector<GenParticle*>::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(std::vector<GenParticle*>::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 * JetFlavourCutFilter::findParticle	(	const GenPartVect &	genPartVect,
		const int	requested_id
	)

private

Definition at line 25 of file JetFlavourCutFilter.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 JetFlavourCutFilter::printHisto	(	const HepMC::GenEvent::particle_iterator	start,
		const HepMC::GenEvent::particle_iterator	end
	)

private

Definition at line 38 of file JetFlavourCutFilter.cc.

References class-composition::child, gather_cfg::cout, end, AlCaHLTBitMon_ParallelJobs::p, and parents.

 {
   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();
       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(std::vector<GenParticle*>::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 JetFlavourCutFilter::jetType

private

Definition at line 46 of file JetFlavourCutFilter.h.

int JetFlavourCutFilter::noAccepted

private

Definition at line 49 of file JetFlavourCutFilter.h.

edm::EDGetTokenT<edm::HepMCProduct> JetFlavourCutFilter::token_

private

Definition at line 48 of file JetFlavourCutFilter.h.

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation