ParticleDecayDrawer Class Reference

Inheritance diagram for ParticleDecayDrawer:

Public Member Functions
	ParticleDecayDrawer (const edm::ParameterSet &)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
	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
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Private Member Functions
bool	accept (const reco::Candidate &, const std::list< const reco::Candidate * > &) const
	accept candidate More...
void	analyze (const edm::Event &, const edm::EventSetup &) override
std::string	decay (const reco::Candidate &, std::list< const reco::Candidate * > &) const
bool	hasValidDaughters (const reco::Candidate &) const
	has valid daughters in the chain More...
std::string	printP4 (const reco::Candidate &) const
	print 4 momenta More...
bool	select (const reco::Candidate &) const
	select candidate More...

Private Attributes
edm::ESHandle< ParticleDataTable >	pdt_
bool	printP4_
	print parameters More...
bool	printPtEtaPhi_
bool	printVertex_
edm::EDGetTokenT< edm::View < reco::Candidate > >	srcToken_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Static Public Member Functions inherited from edm::EDAnalyzer
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

Definition at line 12 of file ParticleDecayDrawer.cc.

Constructor & Destructor Documentation

ParticleDecayDrawer::ParticleDecayDrawer

(

const edm::ParameterSet &

cfg

)

Definition at line 44 of file ParticleDecayDrawer.cc.

                                                                   :
  srcToken_( consumes<edm::View<reco::Candidate> >(cfg.getParameter<InputTag>( "src" ) ) ),
  printP4_( cfg.getUntrackedParameter<bool>( "printP4", false ) ),
  printPtEtaPhi_( cfg.getUntrackedParameter<bool>( "printPtEtaPhi", false ) ),
  printVertex_( cfg.getUntrackedParameter<bool>( "printVertex", false ) ) {
}

Member Function Documentation

bool ParticleDecayDrawer::accept ( const reco::Candidate &  , const std::list< const reco::Candidate * > &    ) const

private

accept candidate

Definition at line 51 of file ParticleDecayDrawer.cc.

References trackerHits::c, spr::find(), and select().

Referenced by Vispa.Gui.BoxContentDialog.BoxContentDialog::apply(), Vispa.Plugins.ConfigEditor.ToolDialog.ToolDialog::apply(), and decay().

                                                                                                      {
  if( find( skip.begin(), skip.end(), & c ) != skip.end() ) return false;
  return select( c );
}

void ParticleDecayDrawer::analyze ( const edm::Event &  event, const edm::EventSetup &  es  )

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 68 of file ParticleDecayDrawer.cc.

References gather_cfg::cout, decay(), edm::EventSetup::getData(), j, m, reco::Candidate::mother(), n, AlCaHLTBitMon_ParallelJobs::p, pdt_, select(), runGlobalFakeInputProducer::skip, and srcToken_.

                                                                              {
  es.getData( pdt_ );
  Handle<View<Candidate> > particles;
  event.getByToken( srcToken_, particles );
  list<const Candidate *> skip;
  vector<const Candidate *> nodes, moms;
  for( View<Candidate>::const_iterator p = particles->begin();
       p != particles->end(); ++ p ) {
    if( p->numberOfMothers() > 1 ) {
      if ( select( * p ) ) {
    skip.push_back( & * p );
    nodes.push_back( & * p );
    for( size_t j = 0; j < p->numberOfMothers(); ++ j ) {
      const Candidate * mom = p->mother( j );
      const Candidate * grandMom;
      while ( ( grandMom = mom->mother() ) != 0 )
        mom = grandMom;
      if ( select( * mom ) ) {
        moms.push_back( mom );
      }
    }
      }
    }
  }
  cout << "-- decay: --" << endl;
  if( moms.size() > 0 ) {
    if ( moms.size() > 1 )
      for( size_t m = 0; m < moms.size(); ++ m ) {
    string dec = decay( * moms[ m ], skip );
    if ( ! dec.empty() )
      cout << "{ " << dec << " } ";
      }
    else
      cout << decay( * moms[ 0 ], skip );
  }
  if ( nodes.size() > 0 ) {
    cout << "-> ";
    if ( nodes.size() > 1 ) {
      for( size_t n = 0; n < nodes.size(); ++ n ) {
    skip.remove( nodes[ n ] );
    string dec = decay( * nodes[ n ], skip );
    if ( ! dec.empty() ) {
      if ( dec.find( "->", 0 ) != string::npos )
        cout << " ( " << dec << " )";
      else
        cout << " " << dec;
        }
      }
    } else {
      skip.remove( nodes[ 0 ] );
      cout << decay( * nodes[ 0 ], skip );
    }
  }
  cout << endl;
}

string ParticleDecayDrawer::decay ( const reco::Candidate &  , std::list< const reco::Candidate * > &    ) const

private

Definition at line 132 of file ParticleDecayDrawer.cc.

References accept(), trackerHits::c, reco::Candidate::daughter(), spr::find(), i, reco::Candidate::numberOfDaughters(), dbtoconf::out, reco::Candidate::pdgId(), pdt_, printP4(), and runGlobalFakeInputProducer::skip.

Referenced by analyze().

                                                          {
  string out;
  if ( find( skip.begin(), skip.end(), & c ) != skip.end() )
    return out;
  skip.push_back( & c );


  int id = c.pdgId();
  const ParticleData * pd = pdt_->particle( id );
  assert( pd != 0 );
  out += ( pd->name() + printP4( c ) );

  size_t validDau = 0, ndau = c.numberOfDaughters();
  for( size_t i = 0; i < ndau; ++ i )
    if ( accept( * c.daughter( i ), skip ) )
      ++ validDau;
  if ( validDau == 0 ) return out;

  out += " ->";

  for( size_t i = 0; i < ndau; ++ i ) {
    const Candidate * d = c.daughter( i );
    if ( accept( * d, skip ) ) {
      string dec = decay( * d, skip );
      if ( dec.find( "->", 0 ) != string::npos )
    out += ( " ( " + dec + " )" );
      else
    out += ( " " + dec );
    }
  }
  return out;
}

bool ParticleDecayDrawer::hasValidDaughters ( const reco::Candidate &  c ) const

private

has valid daughters in the chain

Definition at line 60 of file ParticleDecayDrawer.cc.

References reco::Candidate::daughter(), i, reco::Candidate::numberOfDaughters(), and select().

                                                                           {
  size_t ndau = c.numberOfDaughters();
  for( size_t i = 0; i < ndau; ++ i )
    if ( select( * c.daughter( i ) ) )
      return true;
  return false;
}

string ParticleDecayDrawer::printP4 ( const reco::Candidate &  c ) const

private

print 4 momenta

Definition at line 124 of file ParticleDecayDrawer.cc.

References gather_cfg::cout, reco::Candidate::energy(), reco::Candidate::eta(), reco::Candidate::phi(), printP4_, printPtEtaPhi_, printVertex_, reco::Candidate::pt(), reco::Candidate::px(), reco::Candidate::py(), reco::Candidate::pz(), reco::Candidate::vx(), reco::Candidate::vy(), and reco::Candidate::vz().

Referenced by decay().

                                                               {
  ostringstream cout;
  if ( printP4_ ) cout << " (" << c.px() << ", " << c.py() << ", " << c.pz() << "; " << c.energy() << ")";
  if ( printPtEtaPhi_ ) cout << " [" << c.pt() << ": " << c.eta() << ", " << c.phi() << "]";
  if ( printVertex_ ) cout << " {" << c.vx() << ", " << c.vy() << ", " << c.vz() << "}";
  return cout.str();
}

bool ParticleDecayDrawer::select ( const reco::Candidate &  c ) const

private

select candidate

Definition at line 56 of file ParticleDecayDrawer.cc.

References reco::Candidate::status().

Referenced by accept(), analyze(), and hasValidDaughters().

                                                                {
  return c.status() == 3;
}

Member Data Documentation

edm::ESHandle<ParticleDataTable> ParticleDecayDrawer::pdt_

private

Definition at line 19 of file ParticleDecayDrawer.cc.

Referenced by analyze(), and decay().

bool ParticleDecayDrawer::printP4_

private

print parameters

Definition at line 21 of file ParticleDecayDrawer.cc.

Referenced by printP4().

bool ParticleDecayDrawer::printPtEtaPhi_

private

Definition at line 21 of file ParticleDecayDrawer.cc.

Referenced by printP4().

bool ParticleDecayDrawer::printVertex_

private

Definition at line 21 of file ParticleDecayDrawer.cc.

Referenced by printP4().

edm::EDGetTokenT<edm::View<reco::Candidate> > ParticleDecayDrawer::srcToken_

private

Definition at line 17 of file ParticleDecayDrawer.cc.

Referenced by analyze().