#include <ComphepSingletopFilterPy8.h>

Inheritance diagram for ComphepSingletopFilterPy8:

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

	~ComphepSingletopFilterPy8 () 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 Member Functions
void	beginJob () override

void	endJob () override

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

Private Attributes
int	hardLep

int	pass22

int	pass23

double	ptsep

int	read22

int	read23

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 19 of file ComphepSingletopFilterPy8.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 37 of file ComphepSingletopFilterPy8.cc.

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

                                                                                    :
  token_(consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared")))
 {   
     ptsep = iConfig.getParameter<double>("pTSep");
 }

ComphepSingletopFilterPy8::~ComphepSingletopFilterPy8 ( )

override

Definition at line 43 of file ComphepSingletopFilterPy8.cc.

43 {}

Member Function Documentation

void ComphepSingletopFilterPy8::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 45 of file ComphepSingletopFilterPy8.cc.

References hardLep, pass22, pass23, read22, and read23.

 { 
     read22 = read23 = 0;
     pass22 = pass23 = 0;
     hardLep = 23; //identifies the "hard part" in Pythia8
  }

void ComphepSingletopFilterPy8::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 52 of file ComphepSingletopFilterPy8.cc.

References gather_cfg::cout, cmsPerfStripChart::format, pass22, pass23, read22, and read23.

 {
     cout << "Proc:     2-->2     2-->3     Total" << endl;
     cout << boost::format("Read: %9d %9d %9d") % read22 % read23 % (read22+read23)
          << endl;
     cout << boost::format("Pass: %9d %9d %9d") % pass22 % pass23 % (pass22+pass23)
          << endl;
 }

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

overrideprivate

Definition at line 62 of file ComphepSingletopFilterPy8.cc.

References funct::abs(), MessageLogger_cfi::cerr, class-composition::children, gather_cfg::cout, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), runTauDisplay::gp, hardLep, parents, pass22, pass23, HiggsValidation_cfi::pdg_id, EnergyCorrector::pt, ptsep, read22, read23, and token_.

 {
   
   
 edm::Handle<edm::HepMCProduct> evt;
 //iEvent.getByLabel("generator","unsmeared", evt);
 iEvent.getByToken(token_, evt);
 const HepMC::GenEvent * myEvt = evt->GetEvent();
 
 int id_bdec=0, id_lJet=0, id_b_from_top=0, id_lep = 0;
 //vars for lepton top
 const   GenParticle * gp_clep = nullptr;
 GenVertex * gv_hard = nullptr;
 
 const   GenParticle * gp_lep_FSR = nullptr;
 GenVertex * gv_lep = nullptr;
 vector<const GenParticle *> vgp_lep;
 
 //vars for add b from top
 GenVertex * gv = nullptr;
 const GenParticle * gp = nullptr;
 vector<const GenParticle *> vgp_bsec;
 
 //vars for light Jet (light q)
 GenVertex * gv_lJet= nullptr;
 const GenParticle * gplJet= nullptr;
 vector<const GenParticle *> vgp_lJet;
 
 //vars for b from top
 GenVertex * gv_b_t= nullptr;
 const GenParticle * gp_b_t= nullptr;
 vector<const GenParticle *> vgp_b_t;
 
 
 //Run through all particles in myEvt, if lepton found saves it gp_clep
     for (GenEvent::particle_const_iterator it = myEvt->particles_begin(); it != myEvt->particles_end(); ++it)
    {
     if (abs((*it)->status()) == hardLep) 
       {
             int abs_id = abs((*it)->pdg_id());//11=e, -11=E, 13=mu, -13=Mu, 15=l(tau), -15=L
             if (abs_id==11 || abs_id==13 || abs_id==15) 
             {
             gp_clep = *it;
             id_lep = (*it)->pdg_id();
             gv_lep=(*it)->production_vertex();
         // Lepton FSR
          while (gv_lep) 
           {  
               gp_lep_FSR = nullptr;
           for (GenVertex::particle_iterator it = gv_lep->particles_begin(children); it != gv_lep->particles_end(children); ++it) 
             {
             if ((*it)->pdg_id() == id_lep) 
             {
                   if (!gp_lep_FSR || (*it)->momentum().perp2() > gp_lep_FSR->momentum().perp2()) 
                   {
                   gp_lep_FSR = *it;
                   }
               }
               }      
               if (gp_lep_FSR) 
               {
               gv_lep = gp_lep_FSR->end_vertex();
               vgp_lep.push_back(gp_lep_FSR); //vertex of 
               } 
               else 
               {
               gv_lep = nullptr; //exits the "while" loop
               }
           }
                break;
             }
         }
    }
 // Goes to lepton production vertex
     gv_hard = gp_clep->production_vertex(); 
     
     if (! gp_clep) 
     {
         cout << "ERROR: ComphepSingletopFilterPy8: no charged lepton" << endl;
         return false;
     }
 
 
 //Run through lepton production_vertex 
     for (GenVertex::particle_iterator it = gv_hard->particles_begin(children); it != gv_hard->particles_end(children); ++it) 
    {
         int pdg_id = (*it)->pdg_id(); // KF in Pythia
         int abs_id = abs(pdg_id);
 
 //selection of light quark among particles in primary vertex
     if(abs_id<5) 
     {
         id_lJet = (*it)->pdg_id();
         gv_lJet=(*it)->production_vertex();
         
          while (gv_lJet) 
           {  
               gplJet = nullptr;
           for (GenVertex::particle_iterator it = gv_lJet->particles_begin(children); it != gv_lJet->particles_end(children); ++it) 
             {
             if ((*it)->pdg_id() == id_lJet) 
             {
                   if (!gplJet || (*it)->momentum().perp2() > gplJet->momentum().perp2()) 
                   {
                   gplJet = *it;
                   }
               }
               }
               
               if (gplJet) 
               {
               gv_lJet = gplJet->end_vertex();
               vgp_lJet.push_back(gplJet); //vertex of 
               } 
               else 
               {
               gv_lJet = nullptr; //exits the "while" loop
               }
           }       
     }
 
     if (abs(pdg_id) == 5) // 5 = b
        {
             if (pdg_id * (gp_clep->pdg_id()) < 0) 
            { //b is from top 
                 id_bdec = pdg_id;
 
         id_b_from_top = (*it)->pdg_id();
         gv_b_t = (*it)->production_vertex();
         
          while (gv_b_t) 
           {  
               gp_b_t = nullptr;
           for (GenVertex::particle_iterator it = gv_b_t->particles_begin(children); it != gv_b_t->particles_end(children); ++it) 
             {
             if ((*it)->pdg_id() == id_b_from_top) 
                 {
                   if (!gp_b_t || (*it)->momentum().perp2() > gp_b_t->momentum().perp2()) 
                   {
                   gp_b_t = *it;
                   }
                 }   
             }
               
               if (gp_b_t) 
               {
               gv_b_t = gp_b_t->end_vertex();
               vgp_b_t.push_back(gp_b_t); //vertex of 
               } 
               else 
               {
               gv_b_t = nullptr; //exits the "while" loop
               }
           }  
            } 
            else  
            {//If process 2-3, then aditional b in the initial state fills
                 vgp_bsec.push_back(*it); 
            }
         }
     }
 
 
 bool process22 = (vgp_bsec.empty()); //if there is no aditional b-quark in primary vexrtes, then it is tq-process (2->2) 
 
     if (process22) 
     {
         for (GenVertex::particle_iterator it = gv_hard->particles_begin(parents); it != gv_hard->particles_end(parents); ++it)  //Among parents of lepton production vertex(primary vertex) we find b-quark.
        {
 
             if ((*it)->pdg_id() == id_bdec) 
             {
                 gv = (*it)->production_vertex();
                break;
             }
         }
         if (! gv) 
         {
             cerr << "ERROR: ComphepSingletopFilterPy8: HepMC inconsistency (! gv)" << endl;
             myEvt->print();
             return false;
         }
     } 
    else 
     {
         gv = vgp_bsec.back()->end_vertex();
     }
 
 //##Correction for GV vertex, because of additional gluons in ISR
 bool WeFoundAdditional_b_quark = false;
 int loopCount = 0, gv_loopCount = 0;
     while (WeFoundAdditional_b_quark != true) 
     {     
         for (GenVertex::particle_iterator it = gv->particles_begin(children); it != gv->particles_end(children); ++it) 
        {
            if ((*it)->pdg_id() == -id_bdec) 
            { //we found right vertex of ISR gluon splitting!
          WeFoundAdditional_b_quark = true;//gv = (*it)->production_vertex();
            }
        }
        if (WeFoundAdditional_b_quark == false)
        {
         for (GenVertex::particle_iterator it = gv->particles_begin(parents); it != gv->particles_end(parents); ++it) 
          {
            if ((*it)->pdg_id() == id_bdec) 
         {
           gv = (*it)->production_vertex();
 
         }
          }
        }
       loopCount++;
 
       if (loopCount > 100)//loop protection, nothing more 
       {
 cerr << "ERROR: ComphepSingletopFilterPy8: HepMC inconsistency (No add b vertex found)" << endl;
     break;
       }
     }
 
 
 
     while (gv) 
     {
         gp = nullptr;
 
         for (GenVertex::particle_iterator it = gv->particles_begin(children); it != gv->particles_end(children); ++it) 
        {
            if ((*it)->pdg_id() == -id_bdec) 
            {
                 if (!gp || (*it)->momentum().perp2() > gp->momentum().perp2()) 
                 {
                     gp = *it;
             
                 }
             }
         }
         if (gp) 
         {
             gv = gp->end_vertex();
             vgp_bsec.push_back(gp); //vertex of 
         } 
         else 
         {
             gv = nullptr; //exits the "while" loop
         }
         gv_loopCount++;
     }
 
 
     if (vgp_bsec.empty()) 
     {
         cerr << "ERROR: ComphepSingletopFilterPy8: HepMC inconsistency (vgp_bsec.size() == 0)" << endl;
         return false;
     }
     
     double pt = vgp_bsec.back()->momentum().perp();
    // double eta = vgp_bsec.back()->momentum().eta();  
    bool pass;
     if (process22) 
    {
         read22 += 1;
         pass = pt < ptsep;
         if (pass) pass22 += 1;
     }
     else
     {
         read23 += 1;
         pass = ptsep <= pt;
         if (pass) pass23 += 1;
     }
  return pass;
 }