#include <BsJpsiPhiFilter.h>

Inheritance diagram for BsJpsiPhiFilter:

Classes
struct	CutStruct

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

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

	~BsJpsiPhiFilter () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDFilter () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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	hadronCuts

CutStruct	leptonCuts

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
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

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

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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 BsJpsiPhiFilter.h.

Member Typedef Documentation

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

private

Definition at line 44 of file BsJpsiPhiFilter.h.

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

private

Definition at line 45 of file BsJpsiPhiFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 9 of file BsJpsiPhiFilter.cc.

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

 {
   token_ = consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared"));
   hadronCuts.type = iConfig.getParameter< int >("hadronType");
   hadronCuts.etaMin = iConfig.getParameter<double>("hadronEtaMin");
   hadronCuts.etaMax = iConfig.getParameter<double>("hadronEtaMax");
   hadronCuts.ptMin = iConfig.getParameter<double>("hadronPtMin");
   leptonCuts.type = iConfig.getParameter< int >("leptonType");
   leptonCuts.etaMin = iConfig.getParameter<double>("leptonEtaMin");
   leptonCuts.etaMax = iConfig.getParameter<double>("leptonEtaMax");
   leptonCuts.ptMin = iConfig.getParameter<double>("leptonPtMin");
 
   noAccepted = 0;
 }

BsJpsiPhiFilter::~BsJpsiPhiFilter ( )

override

Definition at line 25 of file BsJpsiPhiFilter.cc.

References gather_cfg::cout.

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

Member Function Documentation

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

private

Definition at line 148 of file BsJpsiPhiFilter.cc.

References funct::abs(), gather_cfg::cout, PVValHelper::eta, BsJpsiPhiFilter::CutStruct::etaMax, BsJpsiPhiFilter::CutStruct::etaMin, AlCaHLTBitMon_ParallelJobs::p, HiggsValidation_cfi::pdg_id, BsJpsiPhiFilter::CutStruct::ptMin, and BsJpsiPhiFilter::CutStruct::type.

 {
   HepMC::GenVertex* myVertex = jpsi->end_vertex();
   int numChildren = myVertex->particles_out_size();
   std::vector<HepMC::GenParticle*> psiChild;
   for(std::vector<GenParticle*>::const_iterator p = myVertex->particles_out_const_begin();
       p != myVertex->particles_out_const_end(); p++) 
     psiChild.push_back((*p));
   
   if(numChildren>1) {
     cout << psiChild[0]->pdg_id()<<" "<<psiChild[1]->pdg_id()<<endl;
   
     if (psiChild.size()==2 && (abs(psiChild[0]->pdg_id()) == cut.type) &&
     (abs(psiChild[1]->pdg_id()) == cut.type))
       {
     cout << psiChild[0]->momentum().rho()<<" "<<psiChild[0]->momentum().eta()
          <<" "<<psiChild[1]->momentum().rho()<<" "<<psiChild[1]->momentum().eta()<<endl;
     return ( (etaInRange(psiChild[0]->momentum().eta(), cut.etaMin, cut.etaMax)) &&
          (etaInRange(psiChild[1]->momentum().eta(), cut.etaMin, cut.etaMax)) &&
          (psiChild[0]->momentum().perp()> cut.ptMin) &&
          (psiChild[1]->momentum().perp()> cut.ptMin));
       }
     return false;
   }
   return false;
 }

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

private

Definition at line 175 of file BsJpsiPhiFilter.cc.

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

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

override

Definition at line 72 of file BsJpsiPhiFilter.cc.

References gather_cfg::cout, egammaForCoreTracking_cff::cuts, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), and jpsi.

 {
   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 = nullptr;
     GenParticle * phi = nullptr;
     // cout << "bs size "<<bsChild.size()<<endl;
     //***
     int numChildren = outVertex->particles_out_size();
     cout<< "bs size "<<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, 443))!=nullptr) && 
     ((phi = findParticle(outVertex, 333))!=nullptr)) {
       
       cout << jpsi->momentum().rho()<<" "<<jpsi->momentum().eta()
        <<" "<<phi->momentum().rho() <<" "<<phi->momentum().eta()<<endl;
       cout <<"bs dec trouve"<<endl;
       if (cuts(phi, hadronCuts) && cuts(jpsi, leptonCuts)) {
         cout <<"OK trouve"<<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 * BsJpsiPhiFilter::findParticle	(	HepMC::GenVertex *	vertex,
		const int	requested_id
	)

private

Definition at line 44 of file BsJpsiPhiFilter.cc.

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

 {
   for(std::vector<GenParticle*>::const_iterator p = vertex->particles_out_const_begin();
       p != vertex->particles_out_const_end(); p++)
     {
       int event_particle_id = abs( (*p)->pdg_id() );
       cout << "isC "<<event_particle_id<<"\n";
       if (requested_id == event_particle_id) return *p;
     }
   return nullptr;
 }

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

private

Definition at line 58 of file BsJpsiPhiFilter.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 == 531) return p;
     }
   return p;  
 }

Member Data Documentation

CutStruct BsJpsiPhiFilter::hadronCuts

private

Definition at line 60 of file BsJpsiPhiFilter.h.

CutStruct BsJpsiPhiFilter::leptonCuts

private

Definition at line 60 of file BsJpsiPhiFilter.h.

int BsJpsiPhiFilter::noAccepted

private

Definition at line 63 of file BsJpsiPhiFilter.h.

edm::EDGetTokenT<edm::HepMCProduct> BsJpsiPhiFilter::token_

private

Definition at line 62 of file BsJpsiPhiFilter.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