PythiaFilterHT Class Reference

#include <IOMC/PythiaFilterHT/src/PythiaFilterHT.cc>

Inheritance diagram for PythiaFilterHT:

Public Member Functions
bool	filter (edm::Event &, const edm::EventSetup &) override
	PythiaFilterHT (const edm::ParameterSet &)
	~PythiaFilterHT () 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 Attributes
edm::EDGetTokenT< edm::HepMCProduct >	label_
double	maxetacut
int	maxnumberofeventsinrun
double	maxpcut
double	maxphicut
double	maxptcut
double	maxrapcut
double	minetacut
double	minhtcut
double	minpcut
double	minphicut
double	minptcut
double	minrapcut
int	motherID
int	particleID
int	processID
double	rapidity
int	status
int	theNumberOfSelected
int	theNumberOfTestedEvt

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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 41 of file PythiaFilterHT.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 9 of file PythiaFilterHT.cc.

References gather_cfg::cout, minhtcut, motherID, theNumberOfSelected, and theNumberOfTestedEvt.

                                                             :
    label_(consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared"))),
        /*minpcut(iConfig.getUntrackedParameter("MinP", 0.)),
    maxpcut(iConfig.getUntrackedParameter("MaxP", 10000.)),
    minptcut(iConfig.getUntrackedParameter("MinPt", 0.)),
    maxptcut(iConfig.getUntrackedParameter("MaxPt", 10000.)),
    minetacut(iConfig.getUntrackedParameter("MinEta", -10.)),
    maxetacut(iConfig.getUntrackedParameter("MaxEta", 10.)),
    minrapcut(iConfig.getUntrackedParameter("MinRapidity", -20.)),
    maxrapcut(iConfig.getUntrackedParameter("MaxRapidity", 20.)),
    minphicut(iConfig.getUntrackedParameter("MinPhi", -3.5)),
    maxphicut(iConfig.getUntrackedParameter("MaxPhi", 3.5)),*/
    minhtcut(iConfig.getUntrackedParameter("MinHT", 0.)),
    motherID(iConfig.getUntrackedParameter("MotherID", 0)) {

        theNumberOfTestedEvt = 0;
        theNumberOfSelected = 0;

        cout << " Cut Definition: " << endl;
        /*cout << " MinP  = " << minpcut << endl;
        cout << " MaxP = " << maxpcut << endl;
        cout << " MinPt = " << minptcut << endl;
        cout << " MaxPt = " << maxptcut << endl;
        cout << " MinEta = " << minetacut << endl;
        cout << " MaxEta = " << maxetacut << endl;
        cout << " MinRapidity = " << minrapcut << endl;
        cout << " MaxRapidity = " << maxrapcut << endl;
        cout << " MinPhi = " << minphicut << endl;
        cout << " MaxPhi = " << maxphicut << endl;*/
        cout << " MinHT = " << minhtcut << endl;
        cout << " MotherID = " << motherID << endl;

}

PythiaFilterHT::~PythiaFilterHT ( )

override

Definition at line 44 of file PythiaFilterHT.cc.

References gather_cfg::cout, theNumberOfSelected, and theNumberOfTestedEvt.

                                {
    std::cout << "Total number of tested events = " << theNumberOfTestedEvt << std::endl;
    std::cout << "Total number of accepted events = " << theNumberOfSelected << std::endl;
}

Member Function Documentation

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

override

Definition at line 55 of file PythiaFilterHT.cc.

References funct::abs(), cms::dd::accepted(), gather_cfg::cout, GenParticle::GenParticle, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), label_, minhtcut, motherID, AlCaHLTBitMon_ParallelJobs::p, theNumberOfSelected, and theNumberOfTestedEvt.

                                                                         {
  
    theNumberOfTestedEvt++;
    if(theNumberOfTestedEvt%1000 == 0) cout << "Number of tested events = " << theNumberOfTestedEvt <<  endl;

    bool accepted = false;
    Handle<HepMCProduct> evt;
    iEvent.getByToken(label_, evt);

    const HepMC::GenEvent * myGenEvent = evt->GetEvent();

    double HT = 0;

    for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end(); ++p ) {

        if ( ( (*p)->status() == 23 ) && ( ( abs( (*p)->pdg_id() ) < 6 ) || ( (*p)->pdg_id() == 21 ) ) ) {
              
            /*rapidity = 0.5*log( ((*p)->momentum().e()+(*p)->momentum().pz()) / ((*p)->momentum().e()-(*p)->momentum().pz()) );

            if ( (*p)->momentum().rho() > minpcut 
                    && (*p)->momentum().rho() < maxpcut
                    && (*p)->momentum().perp() > minptcut 
                    && (*p)->momentum().perp() < maxptcut
                    && (*p)->momentum().eta() > minetacut
                    && (*p)->momentum().eta() < maxetacut 
                    && rapidity > minrapcut
                    && rapidity < maxrapcut 
                    && (*p)->momentum().phi() > minphicut
                    && (*p)->momentum().phi() < maxphicut ) {*/

                if ( motherID == 0 ) {
                   HT += (*p)->momentum().perp();
                } else {
                    HepMC::GenParticle* mother = (*((*p)->production_vertex()->particles_in_const_begin()));
                    if (abs(mother->pdg_id()) == abs(motherID)) {
                        HT += (*p)->momentum().perp();
                    }
                }
            //}         
        }
    }
    if ( HT > minhtcut ) accepted = true;

    if (accepted){
        theNumberOfSelected++;
        cout << "========>  Event preselected " << theNumberOfSelected << " HT = " << HT << endl;
        return true; 
    } else {
        //cout << "========>  Event rejected HT = " << HT << endl;
        return false;
    }

}

Member Data Documentation

edm::EDGetTokenT<edm::HepMCProduct> PythiaFilterHT::label_

private

Definition at line 51 of file PythiaFilterHT.h.

Referenced by SequenceTypes.DummyModule::__repr__(), Mixins._Labelable::_findDependencies(), filter(), and Mixins._Labelable::setLabel().

double PythiaFilterHT::maxetacut

private

Definition at line 59 of file PythiaFilterHT.h.

int PythiaFilterHT::maxnumberofeventsinrun

private

Definition at line 72 of file PythiaFilterHT.h.

double PythiaFilterHT::maxpcut

private

Definition at line 54 of file PythiaFilterHT.h.

double PythiaFilterHT::maxphicut

private

Definition at line 62 of file PythiaFilterHT.h.

double PythiaFilterHT::maxptcut

private

Definition at line 57 of file PythiaFilterHT.h.

double PythiaFilterHT::maxrapcut

private

Definition at line 61 of file PythiaFilterHT.h.

double PythiaFilterHT::minetacut

private

Definition at line 58 of file PythiaFilterHT.h.

double PythiaFilterHT::minhtcut

private

Definition at line 56 of file PythiaFilterHT.h.

Referenced by filter(), and PythiaFilterHT().

double PythiaFilterHT::minpcut

private

Definition at line 53 of file PythiaFilterHT.h.

double PythiaFilterHT::minphicut

private

Definition at line 63 of file PythiaFilterHT.h.

double PythiaFilterHT::minptcut

private

Definition at line 55 of file PythiaFilterHT.h.

double PythiaFilterHT::minrapcut

private

Definition at line 60 of file PythiaFilterHT.h.

int PythiaFilterHT::motherID

private

Definition at line 68 of file PythiaFilterHT.h.

Referenced by filter(), and PythiaFilterHT().

int PythiaFilterHT::particleID

private

Definition at line 52 of file PythiaFilterHT.h.

int PythiaFilterHT::processID

private

Definition at line 69 of file PythiaFilterHT.h.

double PythiaFilterHT::rapidity

private

Definition at line 65 of file PythiaFilterHT.h.

int PythiaFilterHT::status

private

Definition at line 67 of file PythiaFilterHT.h.

Referenced by dirstructure.Comparison::__make_image(), and dirstructure.Comparison::__repr__().

int PythiaFilterHT::theNumberOfSelected

private

Definition at line 71 of file PythiaFilterHT.h.

Referenced by filter(), PythiaFilterHT(), and ~PythiaFilterHT().

int PythiaFilterHT::theNumberOfTestedEvt

private

Definition at line 70 of file PythiaFilterHT.h.

Referenced by filter(), PythiaFilterHT(), and ~PythiaFilterHT().