#include <UEDFilters/UEDMultiLeptonFilter/src/UEDMultiLeptonFilter.cc>

Inheritance diagram for UEDMultiLeptonFilter:

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

	~UEDMultiLeptonFilter ()

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 Member Functions
bool	AcceptEvent ()

void	AllVetoedOff (bool inclusive_message)

virtual void	beginJob () override

virtual void	endJob () override

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

bool	isLepton (HepMC::GenVertex::particles_out_const_iterator part)

bool	isLeptonMinus (HepMC::GenVertex::particles_out_const_iterator part)

bool	isLeptonPlus (HepMC::GenVertex::particles_out_const_iterator part)

void	nLeptons (const std::vector< int > &, int &e, int &mu)

Private Attributes
int	nDidileptons

int	nDiee

int	nDiemu

int	nDileptons

int	nDimumu

int	nFilteredEvents

int	nFour2mu2e

int	nFour4e

int	nFour4mu

int	nFourleptons

int	nProcessed

int	nSSdileptons

int	nSSee

int	nSSemu

int	nSSmumu

int	nTri1mu2e

int	nTri2mu1e

int	nTri3e

int	nTri3mu

int	nTrileptons

int	SSDiEFilter

int	SSDiEMuFilter

int	SSDiEMuVetoedFilter

int	SSDiEVetoedFilter

int	SSDiLepFilter

int	SSDiLepVetoedFilter

int	SSDiMuFilter

int	SSDiMuVetoedFilter

int	UseFilter

int	Vetoed1mu2eFilter

int	Vetoed2mu1eFilter

int	Vetoed2mu2eFilter

int	Vetoed3eFilter

int	Vetoed3muFilter

int	Vetoed4eFilter

int	Vetoed4muFilter

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 49 of file UEDMultiLeptonFilter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 126 of file UEDMultiLeptonFilter.cc.

References AllVetoedOff(), nDiee, nDiemu, nDileptons, nDimumu, nFilteredEvents, nFour2mu2e, nFour4e, nFour4mu, nFourleptons, nProcessed, nSSdileptons, nSSee, nSSemu, nSSmumu, nTri1mu2e, nTri2mu1e, nTri3e, nTri3mu, nTrileptons, SSDiEFilter, SSDiEMuFilter, SSDiEMuVetoedFilter, SSDiEVetoedFilter, SSDiLepFilter, SSDiLepVetoedFilter, SSDiMuFilter, SSDiMuVetoedFilter, UseFilter, Vetoed1mu2eFilter, Vetoed2mu1eFilter, Vetoed2mu2eFilter, Vetoed3eFilter, Vetoed3muFilter, Vetoed4eFilter, and Vetoed4muFilter.

                                                                          :
 UseFilter(iConfig.getUntrackedParameter("UseFilter",0)),
 SSDiMuFilter(iConfig.getUntrackedParameter("SSDiMuFilter",0)),
 SSDiMuVetoedFilter(iConfig.getUntrackedParameter("SSDiMuVetoedFilter",0)),
 SSDiEFilter(iConfig.getUntrackedParameter("SSDiEFilter",0)),
 SSDiEVetoedFilter(iConfig.getUntrackedParameter("SSDiEVetoedFilter",0)),
 SSDiEMuFilter(iConfig.getUntrackedParameter("SSDiEMuFilter",0)),
 SSDiEMuVetoedFilter(iConfig.getUntrackedParameter("SSDiEMuVetoedFilter",0)),
 SSDiLepFilter(iConfig.getUntrackedParameter("SSDiLepFilter",0)),
 SSDiLepVetoedFilter(iConfig.getUntrackedParameter("SSDiLepVetoedFilter",0)),
 Vetoed3muFilter(iConfig.getUntrackedParameter("Vetoed3muFilter",0)),
 Vetoed2mu1eFilter(iConfig.getUntrackedParameter("Vetoed2mu1eFilter",0)),
 Vetoed1mu2eFilter(iConfig.getUntrackedParameter("Vetoed1mu2eFilter",0)),
 Vetoed3eFilter(iConfig.getUntrackedParameter("Vetoed3eFilter",0)),
 Vetoed4muFilter(iConfig.getUntrackedParameter("Vetoed4muFilter",0)),
 Vetoed2mu2eFilter(iConfig.getUntrackedParameter("Vetoed2mu2eFilter",0)),
 Vetoed4eFilter(iConfig.getUntrackedParameter("Vetoed4eFilter",0))
 {
    //now do what ever initialization is needed
 
    if(UseFilter==0){//std::cout << "************ No Filter ************" << std::endl;
    }else{
 //     std::cout <<"Filters On" << std::endl;
      if(SSDiLepFilter+SSDiMuFilter+SSDiEFilter+SSDiEMuFilter>1){
 //     std::cout << "Bad configuration: more than one non-vetoed filter is On" << std::endl;
 //     std::cout << "You should use only one inclusive filter" << std::endl;
 //     std::cout << "All events will pass: UseFilter=0"<< std::endl;
      UseFilter=0;
      }
     }
 
    if(UseFilter==1){ 
      if(SSDiLepFilter==1){
 //       std::cout << "SSDiLepFilter: On " << std::endl;
        //To avoid multiple counting
        AllVetoedOff(true);
      }
      if(SSDiLepVetoedFilter==1){
        SSDiMuVetoedFilter=1;
        SSDiEMuVetoedFilter=1;   
        SSDiEVetoedFilter=1;
      }
      if(SSDiMuFilter==1){
 //       std::cout <<"SSDiMuFilter: On" << std::endl;
        AllVetoedOff(true);
        SSDiMuVetoedFilter=1;
        Vetoed3muFilter=1;
        Vetoed4muFilter=1;
      }
      if(SSDiEFilter==1){
 //       std::cout <<"SSDiEFilter: On" << std::endl;
        AllVetoedOff(true);
        SSDiEVetoedFilter=1;
        Vetoed3eFilter=1;
        Vetoed4eFilter=1;
      }
      if(SSDiEMuFilter==1){
 //       std::cout <<"SSDiEMuFilter: On" << std::endl;       
        AllVetoedOff(true);
        SSDiEMuVetoedFilter=1;
        Vetoed2mu1eFilter=1;
        Vetoed1mu2eFilter=1;
        Vetoed2mu2eFilter=1; 
      } 
    }
 
 
    nFilteredEvents = 0;
    nProcessed = 0;
 
    nDileptons = 0;
 
    nDimumu = 0;
    nDiee = 0;
    nDiemu = 0;  
 
    nSSdileptons = 0;
 
    nSSmumu = 0;
    nSSee = 0;
    nSSemu = 0;
 
    nTrileptons = 0;
 
    nTri3mu = 0;
    nTri2mu1e = 0;
    nTri1mu2e = 0;
    nTri3e = 0;
 
    nFourleptons = 0;
 
    nFour4mu = 0;
    nFour4e = 0;
    nFour2mu2e = 0;
 
 }

UEDMultiLeptonFilter::~UEDMultiLeptonFilter ( )

Definition at line 244 of file UEDMultiLeptonFilter.cc.

 {
  
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 /*
 
    std::cout << "nProcessed: " << nProcessed << std::endl;
    std::cout << "nFilteredEvents: " << nFilteredEvents << std::endl;
 
    std::cout << "nDileptons: " << nDileptons << std::endl;
 
    std::cout << "nDimumu: " << nDimumu << std::endl;
    std::cout << "nDiemu: " << nDiemu << std::endl;
    std::cout << "nDiee: " << nDiee << std::endl;
 
 
    std::cout << "nSSdileptons: " << nSSdileptons << std::endl;
 
    std::cout << "nSSmumu: " << nSSmumu << std::endl;
    std::cout << "nSSemu: " << nSSemu << std::endl;
    std::cout << "nSSee: " << nSSee << std::endl;   
 
    std::cout << "nTrileptons: " << nTrileptons << std::endl;
 
    std::cout << "nTri3mu: " << nTri3mu << std::endl; 
    std::cout << "nTri2mu1e: " << nTri2mu1e << std::endl;
    std::cout << "nTri1mu2e: " << nTri1mu2e << std::endl;
    std::cout << "nTri3e: " << nTri3e <<std::endl;
 
    std::cout << "nFourleptons: " << nFourleptons << std::endl;
    std::cout << "nFour4mu: " << nFour4mu << std::endl;
    std::cout << "nFour4e: " << nFour4e << std::endl;
    std::cout << "nFour2mu2e: " << nFour2mu2e << std::endl;
 
 */
 
 }

Member Function Documentation

bool UEDMultiLeptonFilter::AcceptEvent ( )

private

Definition at line 513 of file UEDMultiLeptonFilter.cc.

References nFilteredEvents.

Referenced by filter().

                                  {
      nFilteredEvents++;
 //     std::cout << "Event accepted" << std::endl;
      return true;
 }

void UEDMultiLeptonFilter::AllVetoedOff ( bool inclusive_message )

private

Definition at line 224 of file UEDMultiLeptonFilter.cc.

References SSDiEMuVetoedFilter, SSDiEVetoedFilter, SSDiLepVetoedFilter, SSDiMuVetoedFilter, Vetoed1mu2eFilter, Vetoed2mu1eFilter, Vetoed2mu2eFilter, Vetoed3eFilter, Vetoed3muFilter, Vetoed4eFilter, and Vetoed4muFilter.

Referenced by UEDMultiLeptonFilter().

                                                         {
 
 //      if(inclusive_message)std::cout << "All Vetoed Filters Off" << std::endl;
 
       SSDiMuVetoedFilter=0;
       SSDiEVetoedFilter=0;
       SSDiEMuVetoedFilter=0;
       SSDiLepVetoedFilter=0;
 
       Vetoed3muFilter=0;
       Vetoed2mu1eFilter=0;
       Vetoed1mu2eFilter=0;
       Vetoed3eFilter=0;
 
       Vetoed4muFilter=0;
       Vetoed2mu2eFilter=0;
       Vetoed4eFilter=0;
 
 }

void UEDMultiLeptonFilter::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 474 of file UEDMultiLeptonFilter.cc.

475 {

476 }

void UEDMultiLeptonFilter::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDFilter.

Definition at line 480 of file UEDMultiLeptonFilter.cc.

480 {

481 }

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

overrideprivatevirtual

Implements edm::EDFilter.

Definition at line 291 of file UEDMultiLeptonFilter.cc.

References accept(), AcceptEvent(), alignCSCRings::e, edm::EventSetup::get(), edm::Event::getByLabel(), i, isLepton(), isLeptonMinus(), isLeptonPlus(), RPCpg::mu, nDiee, nDiemu, nDileptons, nDimumu, nFilteredEvents, nFour2mu2e, nFour4e, nFour4mu, nFourleptons, nLeptons(), nProcessed, nSSdileptons, nSSee, nSSemu, nSSmumu, nTri1mu2e, nTri2mu1e, nTri3e, nTri3mu, nTrileptons, SSDiEMuVetoedFilter, SSDiEVetoedFilter, SSDiLepFilter, SSDiMuVetoedFilter, UseFilter, Vetoed1mu2eFilter, Vetoed2mu1eFilter, Vetoed2mu2eFilter, Vetoed3eFilter, Vetoed3muFilter, Vetoed4eFilter, and Vetoed4muFilter.

 {
 
    using namespace edm;
 
    std::map<std::string, int> accept; 
 
 //   std::cout << "===================== New event ========================: " << nFilteredEvents  << std::endl;
 
    nProcessed++;
    if(UseFilter==0){
      nFilteredEvents++;
      //std::cout << "Event accepted" << std::endl;
      return true;
    }
    std::vector<int> leptonIDminus;
    std::vector<int> leptonIDplus;
 
    Handle<HepMCProduct> mcEventHandle;
 
    iEvent.getByLabel("generator","unsmeared",mcEventHandle);
 
    const HepMC::GenEvent* mcEvent = mcEventHandle->GetEvent() ;
 
    HepMC::GenEvent::particle_const_iterator i;
    for(i = mcEvent->particles_begin(); i!= mcEvent->particles_end(); i++){
 
    bool D_part = (fabs((*i)->pdg_id())-5100000>0 && fabs((*i)->pdg_id())-5100000<40);
    bool S_part = (fabs((*i)->pdg_id())-6100000>0 && fabs((*i)->pdg_id())-6100000<40 );
      if(D_part || S_part){
 
        //debug std::cout << "UED particle ID: " << (*i)->pdg_id() << " status: "<< (*i)->status()<< std::endl;
        HepMC::GenVertex* vertex = (*i)->end_vertex();
        if(vertex!=0){
        for(HepMC::GenVertex::particles_out_const_iterator part = vertex->particles_out_const_begin();
      part != vertex->particles_out_const_end(); part++ ){
           //debug std::cout << "   Outgoing particle id :"<< (*part)->pdg_id() << " status:  " << (*part)->status() << std::endl;
             if(isLepton(part)){
                 if((*part)->status()==1){//debug std::cout << "Final State Lepton " << std::endl;
                 }
                 HepMC::GenVertex* lepton_vertex = (*part)->end_vertex();
                 if(lepton_vertex!=0){
                 for(HepMC::GenVertex::particles_out_const_iterator lepton_part = lepton_vertex->particles_out_const_begin();
      lepton_part != lepton_vertex->particles_out_const_end(); lepton_part++ ){
                         //debug std::cout << "      Part Id: "<< (*lepton_part)->pdg_id() << " status: "<< (*lepton_part)->status()<< std::endl;
                         if((*part)->pdg_id() == (*lepton_part)->pdg_id()){//debug std::cout << "         Part pt: " << (*part)->momentum().perp() << " mu 1: " << (*lepton_part)->momentum().perp() << std::endl;
                         }
                         //std::cout << "         Part px: " << (*part)->momentum().px() << " end_vtx px: " << (*mu_part)->momentum().px() << std::endl;
                         if((*part)->pdg_id() == (*lepton_part)->pdg_id() && (*lepton_part)->status()==1){
                             if(isLeptonPlus(lepton_part))leptonIDplus.push_back((*lepton_part)->pdg_id());
                             if(isLeptonMinus(lepton_part))leptonIDminus.push_back((*lepton_part)->pdg_id());
                           }
                    }
                  }
               }
           }
           }
       }
      }
 
         if((leptonIDplus.size()==1 && leptonIDminus.size()==1) || (leptonIDplus.size()==1 && leptonIDminus.size()==1)){
         //Opposite sign only
         nDileptons++; 
           int e=0;
           int mu=0;
 
           nLeptons(leptonIDplus,e,mu);
           nLeptons(leptonIDminus,e,mu);
           if(mu==2 && e==0)nDimumu++;
           if(mu==1 && e==1)nDiemu++;
           if(mu==0 && e==2)nDiee++;
 
         }
 
 
      if(leptonIDplus.size()>1 || leptonIDminus.size()>1){
 //  std::cout << "SS dilepton" << std::endl;
         //Check dimuon inclusive filter
         nSSdileptons++;
         if(SSDiLepFilter==1)return AcceptEvent();
         //Vetoed decays
         if((leptonIDplus.size()==2 && leptonIDminus.size()==0) || (leptonIDplus.size()==0 && leptonIDminus.size()==2)){
 
           int e=0;
           int mu=0;      
 
           nLeptons(leptonIDplus,e,mu);
           if(mu==2 && e==0){
            nSSmumu++;
            if(SSDiMuVetoedFilter==1)return AcceptEvent();
           }
           if(mu==1 && e==1){
            nSSemu++;
            if(SSDiEMuVetoedFilter==1)return AcceptEvent();
           }
           if(mu==0 && e==2){
        nSSee++;
        if(SSDiEVetoedFilter==1)return AcceptEvent(); 
       }
 
           e=0;
           mu=0;
 
           nLeptons(leptonIDminus,e,mu);
           if(mu==2 && e==0){
            nSSmumu++;
            if(SSDiMuVetoedFilter==1)return AcceptEvent();           
           }
           if(mu==1 && e==1){
        nSSemu++;
            if(SSDiEMuVetoedFilter==1)return AcceptEvent();
           }           
           if(mu==0 && e==2){
            nSSee++;
            if(SSDiEVetoedFilter==1)return AcceptEvent();
           }
         }        
 
      }
 
      if(leptonIDplus.size()+leptonIDminus.size()==4){
 //        std::cout << "Four lepton" << std::endl;
         nFourleptons++;
         int e=0;
         int mu=0;
         nLeptons(leptonIDplus,e,mu);
         nLeptons(leptonIDminus,e,mu);
         if(mu==4 && e==0){
      nFour4mu++;
          if(Vetoed4muFilter==1)return AcceptEvent();         
         }
         if(mu==2 && e==2){
          nFour2mu2e++;
          if(Vetoed2mu2eFilter==1)return AcceptEvent();
     }
         if(mu==0 && e==4){
      nFour4e++;
          if(Vetoed4eFilter==1)return AcceptEvent();  
     }
 
      }
      if(leptonIDplus.size()+leptonIDminus.size()==3){
 //        std::cout << "Three lepton " << std::endl;
         nTrileptons++;
         int e=0;
         int mu=0;
         nLeptons(leptonIDplus,e,mu);
         nLeptons(leptonIDminus,e,mu);
         if(mu==3 && e==0){
      nTri3mu++;
          if(Vetoed3muFilter==1)return AcceptEvent();     
     }
         if(mu==2 && e==1){
          nTri2mu1e++;
          if(Vetoed2mu1eFilter==1)return AcceptEvent();
         }
         if(mu==1 && e==2){
      nTri1mu2e++;
          if(Vetoed1mu2eFilter==1)return AcceptEvent();
     }
         if(mu==0 && e==3){
      nTri3e++;
          if(Vetoed3eFilter==1)return AcceptEvent();
     }
      }
 
 
 
 
 #ifdef THIS_IS_AN_EVENT_EXAMPLE
    Handle<ExampleData> pIn;
    iEvent.getByLabel("example",pIn);
 #endif
 
 #ifdef THIS_IS_AN_EVENTSETUP_EXAMPLE
    ESHandle<SetupData> pSetup;
    iSetup.get<SetupRecord>().get(pSetup);
 #endif
    return false;
 }

bool UEDMultiLeptonFilter::isLepton ( HepMC::GenVertex::particles_out_const_iterator part )

private

Definition at line 484 of file UEDMultiLeptonFilter.cc.

Referenced by filter().

                                                                              {
 
   return (fabs((*part)->pdg_id())==13 || fabs((*part)->pdg_id())==11);
 
 } 

bool UEDMultiLeptonFilter::isLeptonMinus ( HepMC::GenVertex::particles_out_const_iterator part )

private

Definition at line 496 of file UEDMultiLeptonFilter.cc.

Referenced by filter().

                                                                                   {
   return ((*part)->pdg_id()==13 || (*part)->pdg_id()==11);
 }

bool UEDMultiLeptonFilter::isLeptonPlus ( HepMC::GenVertex::particles_out_const_iterator part )

private

Definition at line 491 of file UEDMultiLeptonFilter.cc.

Referenced by filter().

                                                                                  {
   return ((*part)->pdg_id()==-13 || (*part)->pdg_id()==-11);  
 }

void UEDMultiLeptonFilter::nLeptons	(	const std::vector< int > &	leptons_id,
		int &	e,
		int &	mu
	)

private

Definition at line 501 of file UEDMultiLeptonFilter.cc.

References funct::abs(), and i.

Referenced by filter().

                                                                               {
   int nentries = (int)leptons_id.size();
   for(int i=0; i<nentries; i++){
    if(abs(leptons_id.at(i))==11)e++;
    if(abs(leptons_id.at(i))==13)mu++;
   }
 
 //  std::cout <<"mu "<< mu << std::endl;
 //  std::cout <<"e "<< e <<std::endl; 
 }