#include <BdecayFilter.h>

Inheritance diagram for BdecayFilter:

Classes
struct	CutStruct

Public Member Functions
	BdecayFilter (const edm::ParameterSet &)

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

	~BdecayFilter ()

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
bool	cuts (const HepMC::GenParticle *jpsi, const CutStruct &cut)

bool	etaInRange (float eta, float etamin, float etamax)

HepMC::GenParticle *	findParticle (HepMC::GenVertex *, const int requested_id)

HepMC::GenEvent::particle_const_iterator	getNextBs (const HepMC::GenEvent::particle_const_iterator start, const HepMC::GenEvent::particle_const_iterator end)

Private Attributes
CutStruct	firstDaughter

int	motherParticle

int	noAccepted

CutStruct	secondDaughter

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 27 of file BdecayFilter.h.

Member Typedef Documentation

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

private

Definition at line 45 of file BdecayFilter.h.

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

private

Definition at line 46 of file BdecayFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 9 of file BdecayFilter.cc.

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

 {
   token_ = consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared"));
   motherParticle = iConfig.getParameter< int >("motherParticle");
 
   firstDaughter.type = iConfig.getParameter< int >("firstDaughter");
   firstDaughter.decayProduct   = iConfig.getParameter< vector<int> >("firstDaughterDecay");
   firstDaughter.etaMin = iConfig.getParameter<double>("firstDaughterDecayEtaMin");
   firstDaughter.etaMax = iConfig.getParameter<double>("firstDaughterDecayEtaMax");
   firstDaughter.ptMin = iConfig.getParameter<double>("firstDaughterDecayPtMin");
 
   secondDaughter.type = iConfig.getParameter< int >("secondDaughter");
   secondDaughter.decayProduct   = iConfig.getParameter< vector<int> >("secondDaughterDecay");
   secondDaughter.etaMin = iConfig.getParameter<double>("secondDaughterDecayEtaMin");
   secondDaughter.etaMax = iConfig.getParameter<double>("secondDaughterDecayEtaMax");
   secondDaughter.ptMin = iConfig.getParameter<double>("secondDaughterDecayPtMin");
 
   noAccepted = 0;
 }

BdecayFilter::~BdecayFilter ( )

Definition at line 30 of file BdecayFilter.cc.

References gather_cfg::cout.

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

Member Function Documentation

bool BdecayFilter::cuts	(	const HepMC::GenParticle *	jpsi,
		const CutStruct &	cut
	)

private

Definition at line 132 of file BdecayFilter.cc.

References funct::abs(), gather_cfg::cout, BdecayFilter::CutStruct::decayProduct, eta, BdecayFilter::CutStruct::etaMax, BdecayFilter::CutStruct::etaMin, i, j, AlCaHLTBitMon_ParallelJobs::p, and BdecayFilter::CutStruct::ptMin.

 {
     cout << "start cuts" << endl;
   HepMC::GenVertex* myVertex = jpsi->end_vertex();
   int numChildren = myVertex->particles_out_size();
   int numDecProd = cut.decayProduct.size();
   std::vector<HepMC::GenParticle*> psiChild;
   // for(std::set<GenParticle*>::const_iterator p = myVertex->particles_out_const_begin();
   for(GenVertex::particles_out_const_iterator p = myVertex->particles_out_const_begin();
       p != myVertex->particles_out_const_end(); p++) 
     psiChild.push_back((*p));
 
   if (numChildren!=numDecProd) return false;
     cout << psiChild[0]->pdg_id()<<" "<<psiChild[1]->pdg_id()<<endl;
 
   for (int i=0; i<numChildren; ++i) {
     bool goodPart = false;
     for (int j=0; j < numChildren; ++j){
       if (abs(psiChild[i]->pdg_id()) == abs(cut.decayProduct[j])) goodPart = true;
     }
     cout << psiChild[i]->momentum().perp() << endl;
     if ( !goodPart || (!etaInRange(psiChild[i]->momentum().eta(), cut.etaMin, cut.etaMax)) || (psiChild[i]->momentum().perp() < cut.ptMin) ) return false;
   }
     cout << "cuts true" << endl;
   return true;
 }

bool BdecayFilter::etaInRange	(	float	eta,
		float	etamin,
		float	etamax
	)

private

Definition at line 159 of file BdecayFilter.cc.

 {
   return ( (etamin < eta) && (eta < etamax) );
 }

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

virtual

Implements edm::EDFilter.

Definition at line 76 of file BdecayFilter.cc.

References gather_cfg::cout, ParticleFlowFastSim_cff::cuts, GenParticle::GenParticle, edm::Event::getByToken(), jpsi, and phi().

 {
   edm::Handle<HepMCProduct> evt;
   iEvent.getByToken(token_, evt);
   
   const HepMC::GenEvent * generated_event = evt->GetEvent();
   //cout << "Start\n";
 
   bool event_passed = false;
   HepMC::GenEvent::particle_const_iterator bs = getNextBs(generated_event->particles_begin(), generated_event->particles_end());
   while (bs!=  generated_event->particles_end() ) {
 
     // vector< GenParticle * > bsChild = (*bs)->listChildren();
 
     //***
     HepMC::GenVertex* outVertex = (*bs)->end_vertex();
     //***
     
     GenParticle * jpsi = 0;
     GenParticle * phi = 0;
     // cout << "bs size "<<bsChild.size()<<endl;
     //***
     int numChildren = outVertex->particles_out_size();
     cout<< "B children "<<numChildren<<endl;
     //***
     
     /*    if ((bsChild.size()==2) && ((jpsi = findParticle(bsChild, 443))!=0) && 
       ((phi = findParticle(bsChild, 333))!=0)) {
       cout << bsChild[0]->momentum()<<" "<<bsChild[0]->momentum().eta()
       <<" "<<bsChild[1]->momentum()<<" "<<bsChild[1]->momentum().eta()<<endl;
     */
     
     //***
     if( (numChildren==2) && ((jpsi = findParticle(outVertex, firstDaughter.type))!=0) && 
     ((phi = findParticle(outVertex, secondDaughter.type))!=0)) {
       
       cout << jpsi->momentum().rho() <<" "<<jpsi->momentum().eta() <<" "<<phi->momentum().rho()<<" "<<phi->momentum().eta()<<endl;
       //cout <<"bs dec trouve"<<endl;
       if (cuts(phi, secondDaughter) && cuts(jpsi, firstDaughter)) {
         cout <<"decay found"<<endl;
         event_passed = true;
         break;
       }
     }
     bs = getNextBs(++bs, generated_event->particles_end());
   }
   
   if (event_passed) noAccepted++;
   cout << "End filter\n";
   
   delete generated_event; 
   
   return event_passed;
 }

HepMC::GenParticle * BdecayFilter::findParticle	(	HepMC::GenVertex *	vertex,
		const int	requested_id
	)

private

Definition at line 49 of file BdecayFilter.cc.

References funct::abs(), gather_cfg::cout, and AlCaHLTBitMon_ParallelJobs::p.

 {
   // for(std::set<GenParticle*>::const_iterator p = vertex->particles_out_const_begin(); 
   for(GenVertex::particles_out_const_iterator  p = vertex->particles_out_const_begin(); 
       p != vertex->particles_out_const_end(); p++)
     {
       int event_particle_id = abs( (*p)->pdg_id() );
       cout << "particle Id: "<<event_particle_id<<"\n";
       if (requested_id == event_particle_id) return *p;
     }
   return 0;
 }

HepMC::GenEvent::particle_const_iterator BdecayFilter::getNextBs	(	const HepMC::GenEvent::particle_const_iterator	start,
		const HepMC::GenEvent::particle_const_iterator	end
	)

private

Definition at line 63 of file BdecayFilter.cc.

References funct::abs(), end, and AlCaHLTBitMon_ParallelJobs::p.

 {
   HepMC::GenEvent::particle_const_iterator p;
   for (p = start; p != end; p++) 
     {
       int event_particle_id = abs( (*p)->pdg_id() );
 //   cout << "search "<<event_particle_id<<"\n";
       if (event_particle_id == motherParticle) return p;
     }
   return p;  
 }

Member Data Documentation

CutStruct BdecayFilter::firstDaughter

private

Definition at line 61 of file BdecayFilter.h.

int BdecayFilter::motherParticle

private

Definition at line 65 of file BdecayFilter.h.

int BdecayFilter::noAccepted

private

Definition at line 64 of file BdecayFilter.h.

CutStruct BdecayFilter::secondDaughter

private

Definition at line 61 of file BdecayFilter.h.

edm::EDGetTokenT<edm::HepMCProduct> BdecayFilter::token_

private

Definition at line 63 of file BdecayFilter.h.

Classes

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