HeavyQuarkFromMPIFilterAlgo Class Reference

#include <HeavyQuarkFromMPIFilterAlgo.h>

Public Member Functions
bool	filter (const edm::Event &iEvent)
bool	hasMPIAncestor (const reco::GenParticle *)
	HeavyQuarkFromMPIFilterAlgo (const edm::ParameterSet &)
	~HeavyQuarkFromMPIFilterAlgo ()

Private Attributes
edm::InputTag	genParSource_
int	HeavyQuarkFlavour

Detailed Description

Definition at line 23 of file HeavyQuarkFromMPIFilterAlgo.h.

Constructor & Destructor Documentation

HeavyQuarkFromMPIFilterAlgo::HeavyQuarkFromMPIFilterAlgo

(

const edm::ParameterSet &

iConfig

)

Definition at line 10 of file HeavyQuarkFromMPIFilterAlgo.cc.

References edm::ParameterSet::getParameter().

                                                                                       { 
  HeavyQuarkFlavour=(int)iConfig.getParameter<int>("HQFlavour");
  genParSource_=iConfig.getParameter<edm::InputTag>("genParSource");
}

HeavyQuarkFromMPIFilterAlgo::~HeavyQuarkFromMPIFilterAlgo

(

)

Definition at line 15 of file HeavyQuarkFromMPIFilterAlgo.cc.

                                                          {
}

Member Function Documentation

bool HeavyQuarkFromMPIFilterAlgo::filter

(

const edm::Event &

iEvent

)

Definition at line 21 of file HeavyQuarkFromMPIFilterAlgo.cc.

References funct::abs(), edm::RefVector< C, T, F >::begin(), edm::RefVector< C, T, F >::end(), edm::Event::getByLabel(), reco::CompositeRefCandidateT< D >::motherRefVector(), reco::LeafCandidate::pdgId(), and reco::LeafCandidate::status().

                                                                {
  bool Veto=false;
  Handle<reco::GenParticleCollection> genParsHandle;
  iEvent.getByLabel(genParSource_,genParsHandle);
  reco::GenParticleCollection genPars=*genParsHandle;
//  std::cout<<" in the filter HQF="<<HeavyQuarkFlavour<<std::endl;
  bool fromMPI=false;
  for (uint32_t ig=0;ig<genPars.size();ig++) {
    reco::GenParticle gp=genPars[ig];
    if (abs(gp.pdgId())==HeavyQuarkFlavour) {
        if(gp.status()>=30 && gp.status()<40){
//      cout<<"Found a B with status of 3X (eta="<<gp.eta()<<")==>fromMPI=true"<<endl;
        fromMPI=true;
    }
    if(gp.status()>=40 && fromMPI==false){
//      cout<<"Found a B with status >= 4X (eta="<<gp.eta()<<")==>Need to check ancestors!"<<endl;
        const reco::GenParticleRefVector& mothers = gp.motherRefVector();
//      cout<<"Note: it has "<<mothers.size()<<" mothers"<<endl;
        for( reco::GenParticleRefVector::const_iterator im = mothers.begin(); im!=mothers.end(); ++im) {
            const reco::GenParticle& part = **im;
//          cout<<"--->Going to a mother, having eta="<<part.eta()<<endl;
            if( hasMPIAncestor( &part) ){
//              cout<<"------>Found one ancestor with status of 3X (eta="<<gp.eta()<<")==>fromMPI=true"<<endl;
                fromMPI=true;
            }
        } 
    }
    if(fromMPI)Veto=true;
        else Veto=false;
    }
  }
//  cout<<"RETURN "<<Veto<<endl;
  return Veto;
}

bool HeavyQuarkFromMPIFilterAlgo::hasMPIAncestor

(

const reco::GenParticle *

particle

)

Definition at line 56 of file HeavyQuarkFromMPIFilterAlgo.cc.

References edm::RefVector< C, T, F >::begin(), edm::RefVector< C, T, F >::end(), reco::CompositeRefCandidateT< D >::motherRefVector(), and reco::LeafCandidate::status().

                                                                                  {
    if( particle->status() >=30 && particle->status()<40 ){
//      cout<<"------->Mother found with eta="<<particle->eta()<<" and status "<<particle->status()<<"==> returning true!"<<endl;
        return true;
    }
//  cout<<"------->in hasMPIAncestor, current particle has eta="<<particle->eta()<<" and status!=3X ==> checking next mothers"<<endl;
    const reco::GenParticleRefVector& mothers = particle->motherRefVector();
//   cout<<"------>Note: it has "<<mothers.size()<<" mothers"<<endl;
    for( reco::GenParticleRefVector::const_iterator im = mothers.begin(); im!=mothers.end(); ++im) {
        const reco::GenParticle& part = **im;
        if( hasMPIAncestor( &part )){
//          cout<<"--------->Found one ancestor with status of 3X ==>returning true"<<endl;
            return true;
        }
    }
//  cout<<"------>No 30's ancestor found ==>returning false"<<endl;
    return false;
}

Member Data Documentation

edm::InputTag HeavyQuarkFromMPIFilterAlgo::genParSource_

private

Definition at line 33 of file HeavyQuarkFromMPIFilterAlgo.h.

int HeavyQuarkFromMPIFilterAlgo::HeavyQuarkFlavour

private

Definition at line 32 of file HeavyQuarkFromMPIFilterAlgo.h.