---

TopDecaySubset Class Reference

#include <TopQuarkAnalysis/TopEventProducers/interface/TopDecaySubset.h>

Inheritance diagram for TopDecaySubset:

List of all members.

Public Types
enum	FillMode { kStable, kME }
	supported modes to fill the new vectors of gen particles More...
Public Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &)
	write TopDecaySubset into the event
	TopDecaySubset (const edm::ParameterSet &)
	default constructor
	~TopDecaySubset ()
	default destructor
Private Member Functions
void	addDaughters (int &idx, const reco::GenParticle::const_iterator part, reco::GenParticleCollection &target, bool recursive=true)
	recursively fill vector for all further dacay particles of a given particle
void	addRadiation (int &idx, const reco::GenParticle::const_iterator part, reco::GenParticleCollection &target)
	fill vector including all radiations from quarks originating from W/top
void	clearReferences ()
	clear references
void	fillOutput (edm::Event &evt, const reco::GenParticleCollection &src, FillMode mode)
	fill output collection depending on whether the W boson is contained in the generator listing or not
void	fillReferences (const reco::GenParticleRefProd &refProd, reco::GenParticleCollection &target)
	fill references for output vector
void	fromFullListing (const reco::GenParticleCollection &src, reco::GenParticleCollection &target, const int &partId, FillMode mode)
	fill output vector with full decay chain with intermediate W bosons
void	fromTruncListing (const reco::GenParticleCollection &src, reco::GenParticleCollection &target, const int &partId, FillMode mode)
	fill output vector with full decay chain w/o intermediate W bosons
reco::Particle::LorentzVector	p4 (const reco::GenParticle::const_iterator part, int statusFlag)
	calculate lorentz vector from input
reco::Particle::LorentzVector	p4 (const std::vector< reco::GenParticle >::const_iterator top, int statusFlag)
	calculate lorentz vector from input (dedicated to top reconstruction)
void	printSource (const reco::GenParticleCollection &src)
	print the whole listing if particle with pdgId is contained in the top decay chain
void	printTarget (reco::GenParticleCollection &target)
	print the whole decay chain if particle with pdgId is contained in the top decay chain
bool	wInDecayChain (const reco::GenParticleCollection &src, const int &partId)
	check whether the W boson is contained in the original gen particle listing or not
Private Attributes
bool	addRadiatedGluons_
	add radiated gluons or not?
int	motherPartIdx_
	index in new evt listing of parts with daughters; has to be set to -1 in produce to deliver consistent results!
std::map< int, std::vector< int > >	refs_
	management of daughter indices for fillRefs
edm::InputTag	src_
	input tag for the genParticle source

---

Detailed Description

Definition at line 15 of file TopDecaySubset.h.

---

Member Enumeration Documentation

enum TopDecaySubset::FillMode

supported modes to fill the new vectors of gen particles

Enumerator:

kStable
kME

Definition at line 21 of file TopDecaySubset.h.

{kStable, kME};

---

Constructor & Destructor Documentation

TopDecaySubset::TopDecaySubset

(

const edm::ParameterSet &

cfg

)

[explicit]

default constructor

Definition at line 14 of file TopDecaySubset.cc.

                                                        :
  addRadiatedGluons_( cfg.getParameter<bool>( "addRadiatedGluons" ) ),
  src_( cfg.getParameter<edm::InputTag>( "src" ) )
{
  // produces a set of gen particle collections following
  // the decay branch of top quarks to the first level of 
  // stable decay products
  produces<reco::GenParticleCollection>(); 
  produces<reco::GenParticleCollection>("ME");
}

TopDecaySubset::~TopDecaySubset

(

)

default destructor

Definition at line 25 of file TopDecaySubset.cc.

{
}

---

Member Function Documentation

void TopDecaySubset::addDaughters	(	int &	idx,
		const reco::GenParticle::const_iterator	part,
		reco::GenParticleCollection &	target,
		bool	recursive = true
	)			[private]

recursively fill vector for all further dacay particles of a given particle

Definition at line 425 of file TopDecaySubset.cc.

References ptr, and refs_.

Referenced by fromFullListing(), and fromTruncListing().

{
  std::vector<int> daughters;
  int idxBuffer = idx;
  for(reco::GenParticle::const_iterator daughter=part->begin(); daughter!=part->end(); ++daughter){
    std::auto_ptr<reco::GenParticle> ptr( new reco::GenParticle( daughter->threeCharge(), daughter->p4(), daughter->vertex(), daughter->pdgId(), daughter->status(), false) );
    target.push_back( *ptr );
    // increment & push index of daughter
    daughters.push_back( ++idx );
    // continue recursively if desired
    if(recursive){
      addDaughters(idx,daughter,target);  
    }
  }  
  if(daughters.size()) {
    refs_[ idxBuffer ] = daughters;
  }
}

void TopDecaySubset::addRadiation	(	int &	idx,
		const reco::GenParticle::const_iterator	part,
		reco::GenParticleCollection &	target
	)			[private]

fill vector including all radiations from quarks originating from W/top

Definition at line 405 of file TopDecaySubset.cc.

References ptr, refs_, and TopDecayID::stable.

Referenced by fromFullListing(), and fromTruncListing().

{
  std::vector<int> daughters;
  int idxBuffer = idx;
  for(reco::GenParticle::const_iterator daughter=part->begin(); daughter!=part->end(); ++daughter){
    if(daughter->status()<=TopDecayID::stable && daughter->pdgId()!=part->pdgId()){
      // skip comment lines and make sure that
      // the particle is not double counted as 
      // dasughter of itself
      std::auto_ptr<reco::GenParticle> ptr(  new reco::GenParticle( daughter->threeCharge(), daughter->p4(), daughter->vertex(), daughter->pdgId(), daughter->status(), false) );
      target.push_back( *ptr );
      daughters.push_back( ++idx ); //push index of daughter
    }
  }  
  if(daughters.size()) {
    refs_[ idxBuffer ] = daughters;
  }
}

void TopDecaySubset::clearReferences

(

)

[private]

clear references

Definition at line 78 of file TopDecaySubset.cc.

References motherPartIdx_, and refs_.

Referenced by fillOutput().

{
  // clear vector of references 
  refs_.clear();  
  // set idx for mother particles to a start value
  // of -1 (the first entry will raise it to 0)
  motherPartIdx_=-1;
}

void TopDecaySubset::fillOutput	(	edm::Event &	evt,
		const reco::GenParticleCollection &	src,
		FillMode	mode
	)			[private]

fill output collection depending on whether the W boson is contained in the generator listing or not

Definition at line 47 of file TopDecaySubset.cc.

References clearReferences(), fillReferences(), fromFullListing(), fromTruncListing(), edm::Event::getRefBeforePut(), kME, kStable, printTarget(), edm::Event::put(), target, TopDecayID::tID, and wInDecayChain().

Referenced by produce().

{
  // get ref product from the event
  const reco::GenParticleRefProd ref = evt.getRefBeforePut<reco::GenParticleCollection>(); 
  // create target vector
  std::auto_ptr<reco::GenParticleCollection> target( new reco::GenParticleCollection );

  // clear existing refs
  clearReferences();  
  // fill output for top branch
  wInDecayChain(src, TopDecayID::tID) ? fromFullListing (src, *target, TopDecayID::tID, mode): fromTruncListing(src, *target, TopDecayID::tID, mode);
  // fill output for anti-top branch
  wInDecayChain(src,-TopDecayID::tID) ? fromFullListing (src, *target,-TopDecayID::tID, mode): fromTruncListing(src, *target,-TopDecayID::tID, mode);
  // fill references
  fillReferences(ref, *target);
  // print full isting of input collection for debuging
  // with 'TopDecaySubset_printTarget'
  printTarget(*target);

  // write vectors to the event
  switch(mode){
  case kStable:
    evt.put(target);
    break;
  case kME:
    evt.put(target, "ME");
    break;
  }
}

void TopDecaySubset::fillReferences	(	const reco::GenParticleRefProd &	refProd,
		reco::GenParticleCollection &	target
	)			[private]

fill references for output vector

Definition at line 445 of file TopDecaySubset.cc.

References reco::CompositeRefCandidateT< D >::addDaughter(), edm::errors::InvalidReference, p, and refs_.

Referenced by fillOutput().

{ 
  int idx=0;
  for(reco::GenParticleCollection::iterator p=sel.begin(); p!=sel.end(); ++p, ++idx){
    //find daughter reference vectors in refs_ and add daughters
    std::map<int, std::vector<int> >::const_iterator daughters=refs_.find( idx );
    if( daughters!=refs_.end() ){
      for(std::vector<int>::const_iterator daughter = daughters->second.begin(); 
          daughter!=daughters->second.end(); ++daughter){
        reco::GenParticle* part = dynamic_cast<reco::GenParticle* > (&(*p));
        if(part == 0){
         throw edm::Exception( edm::errors::InvalidReference, "Not a GenParticle" );
        }
        part->addDaughter( reco::GenParticleRef(ref, *daughter) );
        sel[*daughter].addMother( reco::GenParticleRef(ref, idx) );
      }
    }
  }
}

void TopDecaySubset::fromFullListing	(	const reco::GenParticleCollection &	src,
		reco::GenParticleCollection &	target,
		const int &	partId,
		FillMode	mode
	)			[private]

fill output vector with full decay chain with intermediate W bosons

Definition at line 113 of file TopDecaySubset.cc.

References funct::abs(), addDaughters(), addRadiatedGluons_, addRadiation(), TopDecayID::bID, false, TopDecayID::glueID, kME, motherPartIdx_, p4(), refs_, TopDecayID::stable, t, TopDecayID::tauID, getDQMSummary::td, TopDecayID::tID, true, TopDecayID::unfrag, and TopDecayID::WID.

Referenced by fillOutput().

{
  int statusFlag;
  // determine status flag of the new 
  // particle depending on the FillMode
  mode == kME ? statusFlag=3 : statusFlag=2;

  for(reco::GenParticleCollection::const_iterator t=src.begin(); t!=src.end(); ++t){
    if( t->status() == TopDecayID::unfrag && t->pdgId()==partId ){ 
      // if particle is top or anti-top depending 
      // on the function's argument partId 
      std::auto_ptr<reco::GenParticle> topPtr( new reco::GenParticle( t->threeCharge(), p4( t, statusFlag), t->vertex(), t->pdgId(), statusFlag, false ) );
      target.push_back( *topPtr );
      ++motherPartIdx_;
      // keep the top index for the map to manage the daughter refs
      int iTop=motherPartIdx_; 
      std::vector<int> topDaughters;
      // define the W boson index (to be set later) for the map to 
      // manage the daughter refs
      int iW = 0;
      std::vector<int> wDaughters;
      //iterate over top daughters
      for(reco::GenParticle::const_iterator td=t->begin(); td!=t->end(); ++td){
        if( td->status()==TopDecayID::unfrag && abs( td->pdgId() )<=TopDecayID::bID ){ 
          // if particle is beauty or other quark 
          std::auto_ptr<reco::GenParticle> bPtr( new reco::GenParticle( td->threeCharge(), p4( td, statusFlag ), td->vertex(), td->pdgId(), statusFlag, false ) );
          target.push_back( *bPtr );      
          // increment & push index of the top daughter
          topDaughters.push_back( ++motherPartIdx_ ); 
          if(addRadiatedGluons_){
            addRadiation(motherPartIdx_,td,target); 
          }
        }
        if( td->status()==TopDecayID::unfrag && abs( td->pdgId() )==TopDecayID::WID ){ 
          // if particle is is W boson
          std::auto_ptr<reco::GenParticle> wPtr(  new reco::GenParticle( td->threeCharge(), p4( td, statusFlag), td->vertex(), td->pdgId(), statusFlag, true ) );
          target.push_back( *wPtr );
          // increment & push index of the top daughter
          topDaughters.push_back( ++motherPartIdx_ );
          // keep the W idx for the map
          iW=motherPartIdx_; 
          if(addRadiatedGluons_){
            addRadiation(motherPartIdx_,td,target); 
          }

          // iterate over W daughters
          for(reco::GenParticle::const_iterator wd=td->begin(); wd!=td->end(); ++wd){
            // make sure the W daughter is of status unfrag and not the W itself
            if( wd->status()==TopDecayID::unfrag && !(abs(wd->pdgId())==TopDecayID::WID) ) {
              std::auto_ptr<reco::GenParticle> qPtr( new reco::GenParticle( wd->threeCharge(), p4( wd, statusFlag ), wd->vertex(), wd->pdgId(), statusFlag, false) );
              target.push_back( *qPtr );
              // increment & push index of the top daughter
              wDaughters.push_back( ++motherPartIdx_ );
              if(addRadiatedGluons_){
                addRadiation(motherPartIdx_,wd,target); 
              }
              if( wd->status()==TopDecayID::unfrag && abs( wd->pdgId() )==TopDecayID::tauID ){ 
                // add tau daughters if the particle is a tau pass
                // the daughter of the tau which is of status 2
                addDaughters(motherPartIdx_,wd->begin(),target); 
              }
            } 
          }
        }
        if(td->status()==TopDecayID::stable && ( td->pdgId()==TopDecayID::glueID || abs(td->pdgId())<TopDecayID::bID)){
          // collect additional radiation from the top 
          std::auto_ptr<reco::GenParticle> radPtr( new reco::GenParticle( td->threeCharge(), td->p4(), td->vertex(), td->pdgId(), statusFlag, false ) );
          target.push_back( *radPtr );  
          // increment & push index of the top daughter  
          topDaughters.push_back( ++motherPartIdx_ );
        }
      }
      // add potential sisters of the top quark;
      // only for top to prevent double counting
      if(t->numberOfMothers()>0 && t->pdgId()==TopDecayID::tID){
        for(reco::GenParticle::const_iterator ts = t->mother()->begin(); ts!=t->mother()->end(); ++ts){
          // loop over all daughters of the top mother i.e.
          // the two top quarks and their potential sisters
          if(abs(ts->pdgId())!=t->pdgId()){
            // add all further particles
            // but the two top quarks 
            reco::GenParticle* cand = new reco::GenParticle( ts->threeCharge(), ts->p4(), ts->vertex(), ts->pdgId(), ts->status(), false );
            std::auto_ptr<reco::GenParticle> sPtr( cand );
            target.push_back( *sPtr );
            // increment & add 1. generation of daughters
            addDaughters(++motherPartIdx_,ts->begin(),target,false);
          }
        }
      }
      // fill the map for the administration 
      // of daughter indices
      refs_[ iTop ] = topDaughters;
      refs_[ iW   ] = wDaughters; 
    }
  }
}

void TopDecaySubset::fromTruncListing	(	const reco::GenParticleCollection &	src,
		reco::GenParticleCollection &	target,
		const int &	partId,
		FillMode	mode
	)			[private]

fill output vector with full decay chain w/o intermediate W bosons

Definition at line 211 of file TopDecaySubset.cc.

References funct::abs(), addDaughters(), addRadiatedGluons_, addRadiation(), TopDecayID::bID, reco::Particle::charge(), false, TopDecayID::glueID, kME, motherPartIdx_, reco::Particle::p4(), p4(), refs_, TopDecayID::stable, t, TopDecayID::tauID, getDQMSummary::td, reco::Particle::threeCharge(), TopDecayID::tID, true, TopDecayID::unfrag, and reco::Particle::vertex().

Referenced by fillOutput().

{
  int statusFlag;
  // determine status flag of the new 
  // particle depending on the FillMode
  mode == kME ? statusFlag=3 : statusFlag=2;

  // needed for W reconstruction from 
  // daughters
  reco::Particle::Point wVtx;
  reco::Particle::Charge wQ=0;
  reco::Particle::LorentzVector wP4;
  for(reco::GenParticleCollection::const_iterator t=src.begin(); t!=src.end(); ++t){
    if( t->status() == TopDecayID::unfrag && t->pdgId()==partId ){ 
      // if particle is top or anti-top 
      std::auto_ptr<reco::GenParticle> topPtr( new reco::GenParticle( t->threeCharge(), p4( t, statusFlag), t->vertex(), t->pdgId(), statusFlag, false ) );
      target.push_back( *topPtr );
      ++motherPartIdx_;
      // keep the top index for the map to manage the daughter refs
      int iTop=motherPartIdx_; 
      std::vector<int> topDaughters;
      // define the W boson index (to be set later) for the map to 
      // manage the daughter refs
      int iW = 0;
      std::vector<int> wDaughters;
      
      //iterate over top daughters
      for(reco::GenParticle::const_iterator td=t->begin(); td!=t->end(); ++td){
        if( td->status()==TopDecayID::unfrag && abs( td->pdgId() )==TopDecayID::bID ){ 
          // if particle is beauty or other quark 
          std::auto_ptr<reco::GenParticle> bPtr( new reco::GenParticle( td->threeCharge(), p4( td, statusFlag ), td->vertex(), td->pdgId(), statusFlag, false ) );
          target.push_back( *bPtr );      
          // increment & push index of the top daughter
          topDaughters.push_back( ++motherPartIdx_ ); 
          if(addRadiatedGluons_){
            addRadiation(motherPartIdx_,td,target); 
          }
        }

        // count all 4-vectors but the b quark 
        // as daughters of the W boson
        if( td->status()==TopDecayID::unfrag && abs( td->pdgId() )!=TopDecayID::bID ){
          reco::GenParticle* cand = new reco::GenParticle( td->threeCharge(), p4( td, statusFlag ), td->vertex(), td->pdgId(), statusFlag, false);
          std::auto_ptr<reco::GenParticle> qPtr( cand );
          target.push_back( *qPtr );
          // increment & push index of the top daughter
          wDaughters.push_back( ++motherPartIdx_ );
          if(addRadiatedGluons_){
            if(abs(td->pdgId())<TopDecayID::tID){
              addRadiation(motherPartIdx_,td,target); 
            }
          }

          // reconstruct the quantities of the W boson from its daughters; take 
          // care of the non-trivial charge definition of quarks/leptons and of 
          // the non-integer charge of the quarks
          if( fabs(td->pdgId() )<TopDecayID::tID ) //quark
            wQ += ((td->pdgId()>0)-(td->pdgId()<0))*abs(cand->threeCharge());
          else //lepton
            wQ +=-((td->pdgId()>0)-(td->pdgId()<0))*abs(cand->charge());
          wVtx=cand->vertex();        
          wP4+=cand->p4();

          if( td->status()==TopDecayID::unfrag && abs( td->pdgId() )==TopDecayID::tauID ){ 
            // add tau daughters if the particle is a tau pass
            // the daughter of the tau which is of status 2
            addDaughters(motherPartIdx_,td->begin(),target); 
          }
        }
        if( (td+1)==t->end() ){
          // the reconstruction of the W boson 
          // candidate; now add it to the list
          std::auto_ptr<reco::GenParticle> wPtr(  new reco::GenParticle( td->threeCharge(), p4( td, statusFlag), td->vertex(), td->pdgId(), statusFlag, true ) );
          target.push_back( *wPtr );
          // increment & push index of the top daughter
          topDaughters.push_back( ++motherPartIdx_ );
          // keep the W idx for the map
          iW=motherPartIdx_; 
          if(addRadiatedGluons_){
            addRadiation(motherPartIdx_,td,target); 
          }
          
          // reset the quantities of the W boson 
          // for the next reconstruction cycle
          wQ  = 0;
          wVtx= reco::Particle::Point(0, 0, 0);
          wP4 = reco::Particle::LorentzVector(0, 0, 0, 0);
        }
        if(td->status()==TopDecayID::stable && ( td->pdgId()==TopDecayID::glueID || abs(td->pdgId())<TopDecayID::bID)){
          // collect additional radiation from the top 
          std::auto_ptr<reco::GenParticle> radPtr( new reco::GenParticle( td->threeCharge(), td->p4(), td->vertex(), td->pdgId(), statusFlag, false ) );
          target.push_back( *radPtr );  
          // increment & push index of the top daughter  
          topDaughters.push_back( ++motherPartIdx_ );
        }
      }
      // add potential sisters of the top quark;
      // only for top to prevent double counting
      if(t->numberOfMothers()>0 && t->pdgId()==TopDecayID::tID){
        for(reco::GenParticle::const_iterator ts = t->mother()->begin(); ts!=t->mother()->end(); ++ts){
          // loop over all daughters of the top mother i.e.
          // the two top quarks and their potential sisters
          if(abs(ts->pdgId())!=t->pdgId()){
            // add all further particles
            // but the two top quarks 
            reco::GenParticle* cand = new reco::GenParticle( ts->threeCharge(), ts->p4(), ts->vertex(), ts->pdgId(), ts->status(), false );
            std::auto_ptr<reco::GenParticle> sPtr( cand );
            target.push_back( *sPtr );
            // increment & add 1. generation of daughters
            addDaughters(++motherPartIdx_,ts->begin(),target,false);
          }
        }
      }
      // fill the map for the administration 
      // of daughter indices
      refs_[ iTop ] = topDaughters;
      refs_[ iW   ] = wDaughters; 
    }
  }
}

reco::Particle::LorentzVector TopDecaySubset::p4	(	const reco::GenParticle::const_iterator	part,
		int	statusFlag
	)			[private]

calculate lorentz vector from input

Definition at line 374 of file TopDecaySubset.cc.

References funct::abs(), p, p4(), TopDecayID::stable, TopDecayID::unfrag, and TopDecayID::WID.

{
  // calculate the four vector for all top daughters from 
  // their daughters including additional radiation 
  if(statusFlag==TopDecayID::unfrag){
    // return 4 momentum as it is
    return part->p4();
  }
  reco::Particle::LorentzVector vec;
  for(reco::GenParticle::const_iterator p=part->begin(); p!=part->end(); ++p){
    if( p->status()<=TopDecayID::stable && 
        abs(p->pdgId ())==TopDecayID::WID){
      vec=p->p4();
    }
    else{
      if(p->status()<=TopDecayID::stable){
        // sum up the p4 of all stable particles 
        // (of status 1 or 2)
        vec+=p->p4();
      }
      else 
        if( p->status()==TopDecayID::unfrag)
          // if the particle is unfragmented (i.e.
          // status 3) decend by one level
          vec+=p4(p, statusFlag);   
    }
  }
  return vec;
}

reco::Particle::LorentzVector TopDecaySubset::p4	(	const std::vector< reco::GenParticle >::const_iterator	top,
		int	statusFlag
	)			[private]

calculate lorentz vector from input (dedicated to top reconstruction)

Definition at line 333 of file TopDecaySubset.cc.

References funct::abs(), p, TopDecayID::unfrag, and TopDecayID::WID.

Referenced by fromFullListing(), fromTruncListing(), and p4().

{
  // calculate the four vector for top/anti-top quarks from 
  // the W boson and the b quark plain or including all 
  // additional radiation depending on switch 'plain'
  if(statusFlag==TopDecayID::unfrag){
    // return 4 momentum as it is
    return top->p4();
  }
  reco::Particle::LorentzVector vec;
  for(reco::GenParticle::const_iterator p=top->begin(); p!=top->end(); ++p){
    if( p->status() == TopDecayID::unfrag ){
      // decend by one level for each
      // status 3 particle on the way
      vec+=p4( p, statusFlag );
    }
    else{ 
      if( abs(top->pdgId())==TopDecayID::WID ){
        // in case of a W boson skip the status 
        // 2 particle to prevent double counting
        if( abs(p->pdgId())!=TopDecayID::WID ) 
          vec+=p->p4();
      } 
      else{
        // add all four vectors for each stable
        // particle (status 1 or 2) on the way 
        // else
        vec+=p->p4();
        if( vec.mass()-top->mass()>0 ){
          // continue adding up gluons and qqbar pairs on the top 
          // line until the nominal top mass is reached; then break 
          // in order to prevent picking up virtualities
          break;
        }
      }
    }
  }
  return vec;
}

void TopDecaySubset::printSource

(

const reco::GenParticleCollection &

src

)

[private]

print the whole listing if particle with pdgId is contained in the top decay chain

Definition at line 466 of file TopDecaySubset.cc.

References d, kMAX, funct::log(), p, t, and TopDecayID::tID.

Referenced by produce().

{
  edm::LogVerbatim log("TopDecaySubset_printSource");
  log << "\n   idx   pdg   stat      px          py         pz             mass          daughter pdg's  "
      << "\n===========================================================================================\n";

  for(unsigned int t=0; t<src.size(); ++t){
    if( src[t].pdgId()==TopDecayID::tID ){
      // restrict to the top in order 
      // to have it shown only once 
      int idx=0;
      for(reco::GenParticleCollection::const_iterator p=src.begin(); p!=src.end(); ++p, ++idx){
        // loop the top daughters
        log << std::right << std::setw( 5) << idx
            << std::right << std::setw( 7) << src[idx].pdgId()
            << std::right << std::setw( 5) << src[idx].status() << "  "
            << std::right << std::setw(10) << std::setprecision( 6 ) << src[idx].p4().x() << "  "       
            << std::right << std::setw(10) << std::setprecision( 6 ) << src[idx].p4().y() << "  "       
            << std::right << std::setw(10) << std::setprecision( 6 ) << src[idx].p4().z() << "  "       
            << std::right << std::setw(15) << std::setprecision( 6 ) << src[idx].p4().mass() 
            << "   ";
        // search for potential daughters; if they exits 
        // print the daughter to the screen in the last 
        // column of the table separated by ','
        TString pdgIds;
        unsigned int jdx=0;
        for(reco::GenParticle::const_iterator d=p->begin(); d!=p->end(); ++d, ++jdx){
          if(jdx<kMAX){
            pdgIds+=d->pdgId();
            if(d+1 != p->end()){
              pdgIds+= ",";
            }
          }
          else{
            pdgIds+="...(";
            pdgIds+= p->numberOfDaughters();
            pdgIds+=")";
            break;
          }
        }
        if(idx>0){
          log << std::setfill( ' ' ) << std::right << std::setw(15) << pdgIds; 
          log << "\n";
        }
        else{
          log << std::setfill( ' ' ) << std::right << std::setw(15) << "-\n";
        }
      }
    }
  }
}

void TopDecaySubset::printTarget

(

reco::GenParticleCollection &

target

)

[private]

print the whole decay chain if particle with pdgId is contained in the top decay chain

Definition at line 519 of file TopDecaySubset.cc.

References d, kMAX, funct::log(), p, refs_, t, and TopDecayID::tID.

Referenced by fillOutput().

{
  edm::LogVerbatim log("TopDecaySubset_printTarget");
  log << "\n   idx   pdg   stat      px          py         pz             mass          daughter pdg's  "
      << "\n===========================================================================================\n";

  for(unsigned int t=0; t<sel.size(); ++t){
    if( sel[t].pdgId()==TopDecayID::tID ){
      // restrict to the top in order 
      // to have it shown only once      
      int idx=0;
      for(reco::GenParticleCollection::iterator p=sel.begin(); p!=sel.end(); ++p, ++idx){
        // loop the top daughters
        log << std::right << std::setw( 5) << idx
            << std::right << std::setw( 7) << sel[idx].pdgId()
            << std::right << std::setw( 5) << sel[idx].status() << "  "
            << std::right << std::setw(10) << std::setprecision( 6 ) << sel[idx].p4().x() << "  "       
            << std::right << std::setw(10) << std::setprecision( 6 ) << sel[idx].p4().y() << "  "       
            << std::right << std::setw(10) << std::setprecision( 6 ) << sel[idx].p4().z() << "  "       
            << std::right << std::setw(15) << std::setprecision( 6 ) << sel[idx].p4().mass() 
            << "   ";
        // search for potential daughters; if they exits 
        // print the daughter to the screen in the last 
        // column of the table separated by ','
        TString pdgIds;
        std::map<int, std::vector<int> >::const_iterator daughters=refs_.find( idx );
        if( daughters!=refs_.end() ){
          unsigned int jdx=0;
          for(std::vector<int>::const_iterator d = daughters->second.begin(); d!=daughters->second.end(); ++d, ++jdx){
            if(jdx<kMAX){
              pdgIds+=sel[*d].pdgId();
              if(d+1 != daughters->second.end()){
                pdgIds+= ",";
              }
            }
            else{
              pdgIds+="...(";
              pdgIds+= daughters->second.size();
              pdgIds+=")";
              break;
            }
          }
          log << std::setfill( ' ' ) << std::right << std::setw(15) << pdgIds; 
          log << "\n";
        }
        else{
          log << std::setfill( ' ' ) << std::right << std::setw(15) << "-\n";
        }
      }
    }
  }
}

void TopDecaySubset::produce	(	edm::Event &	evt,
		const edm::EventSetup &	setup
	)			[virtual]

write TopDecaySubset into the event

Implements edm::EDProducer.

Definition at line 30 of file TopDecaySubset.cc.

References fillOutput(), edm::Event::getByLabel(), kME, kStable, printSource(), HLT_VtxMuL3::src, and src_.

{     
  // get source collection
  edm::Handle<reco::GenParticleCollection> src;
  evt.getByLabel(src_, src);

  // print full listing of input collection for debuging
  // with 'TopDecaySubset_printSource'
  printSource(*src);

  // fill output vectors with references
  fillOutput (evt, *src, kME      );
  fillOutput (evt, *src, kStable  );
}

bool TopDecaySubset::wInDecayChain	(	const reco::GenParticleCollection &	src,
		const int &	partId
	)			[private]

check whether the W boson is contained in the original gen particle listing or not

Definition at line 88 of file TopDecaySubset.cc.

References funct::abs(), lhef::isContained(), t, getDQMSummary::td, TopDecayID::unfrag, and TopDecayID::WID.

Referenced by fillOutput().

{
  bool isContained=false;
  for(reco::GenParticleCollection::const_iterator t=src.begin(); t!=src.end(); ++t){
    if( t->status() == TopDecayID::unfrag && t->pdgId()==partId ){ 
      for(reco::GenParticle::const_iterator td=t->begin(); td!=t->end(); ++td){
        if( td->status()==TopDecayID::unfrag && abs( td->pdgId() )==TopDecayID::WID ){ 
          isContained=true;
          break;
        }
      }
    }
  }
  if( !isContained ){
    edm::LogWarning( "decayChain" )
      << " W boson is not contained in decay chain in the gen particle listing \n"
      << " of the input collection. A status 2 equivalent W candidate will be  \n"
      << " re-reconstructed from the W daughters but the hadronization of the  \n"
      << " W boson might be screwed up. Contact an expert for the generator in \n"
      << " use to assure that what you are doing makes sense.";     
  }
  return isContained;
}

---

Member Data Documentation

bool TopDecaySubset::addRadiatedGluons_ [private]

add radiated gluons or not?

Definition at line 76 of file TopDecaySubset.h.

Referenced by fromFullListing(), and fromTruncListing().

int TopDecaySubset::motherPartIdx_ [private]

index in new evt listing of parts with daughters; has to be set to -1 in produce to deliver consistent results!

Definition at line 72 of file TopDecaySubset.h.

Referenced by clearReferences(), fromFullListing(), and fromTruncListing().

std::map<int,std::vector<int> > TopDecaySubset::refs_ [private]

management of daughter indices for fillRefs

Definition at line 74 of file TopDecaySubset.h.

Referenced by addDaughters(), addRadiation(), clearReferences(), fillReferences(), fromFullListing(), fromTruncListing(), and printTarget().

edm::InputTag TopDecaySubset::src_ [private]

input tag for the genParticle source

Definition at line 78 of file TopDecaySubset.h.

Referenced by produce().

---

The documentation for this class was generated from the following files:

• TopQuarkAnalysis/TopEventProducers/interface/TopDecaySubset.h
• TopQuarkAnalysis/TopEventProducers/src/TopDecaySubset.cc

---