reco::FlavorHistoryEvent Class Reference

#include <FlavorHistoryEvent.h>

Public Types
typedef std::vector< value_type >	collection_type
typedef collection_type::const_iterator	const_iterator
typedef collection_type::const_pointer	const_pointer
typedef collection_type::const_reference	const_reference
typedef collection_type::const_reverse_iterator	const_reverse_iterator
typedef FlavorHistory::FLAVOR_T	flavor_type
typedef collection_type::iterator	iterator
typedef collection_type::pointer	pointer
typedef collection_type::reference	reference
typedef collection_type::reverse_iterator	reverse_iterator
typedef collection_type::size_type	size_type
typedef FlavorHistory	value_type

Public Member Functions
const_iterator	begin () const
void	cache ()
void	clear ()
double	deltaR () const
const_iterator	end () const
	FlavorHistoryEvent ()
flavor_type	flavorSource () const
unsigned int	highestFlavor () const
bool	isCached () const
unsigned int	nb () const
unsigned int	nc () const
void	push_back (const value_type &v)
const_reverse_iterator	rbegin () const
const_reverse_iterator	rend () const
void	resize (size_t n)
size_type	size () const
	~FlavorHistoryEvent ()

Protected Attributes
bool	cached_
double	dR_
flavor_type	flavorSource_
unsigned int	highestFlavor_
collection_type	histories_
unsigned int	nb_
unsigned int	nc_

Detailed Description

Stores a vector of information about the flavor history of partons as well as a classification of the event.

It will return the following: nb = number of genjets that are matched to b partons nc = number of genjets that are matched to c partons

This can be used to classify the event, for instance, as W+bb (2 b partons), W+c (1 c parton), etc.

Author

: Salvatore Rappoccio (JHU)

Definition at line 84 of file FlavorHistoryEvent.h.

Member Typedef Documentation

typedef std::vector<value_type> reco::FlavorHistoryEvent::collection_type

Definition at line 88 of file FlavorHistoryEvent.h.

typedef collection_type::const_iterator reco::FlavorHistoryEvent::const_iterator

Definition at line 91 of file FlavorHistoryEvent.h.

typedef collection_type::const_pointer reco::FlavorHistoryEvent::const_pointer

Definition at line 95 of file FlavorHistoryEvent.h.

typedef collection_type::const_reference reco::FlavorHistoryEvent::const_reference

Definition at line 97 of file FlavorHistoryEvent.h.

typedef collection_type::const_reverse_iterator reco::FlavorHistoryEvent::const_reverse_iterator

Definition at line 93 of file FlavorHistoryEvent.h.

typedef FlavorHistory::FLAVOR_T reco::FlavorHistoryEvent::flavor_type

Definition at line 98 of file FlavorHistoryEvent.h.

typedef collection_type::iterator reco::FlavorHistoryEvent::iterator

Definition at line 90 of file FlavorHistoryEvent.h.

typedef collection_type::pointer reco::FlavorHistoryEvent::pointer

Definition at line 94 of file FlavorHistoryEvent.h.

typedef collection_type::reference reco::FlavorHistoryEvent::reference

Definition at line 96 of file FlavorHistoryEvent.h.

typedef collection_type::reverse_iterator reco::FlavorHistoryEvent::reverse_iterator

Definition at line 92 of file FlavorHistoryEvent.h.

typedef collection_type::size_type reco::FlavorHistoryEvent::size_type

Definition at line 89 of file FlavorHistoryEvent.h.

typedef FlavorHistory reco::FlavorHistoryEvent::value_type

Definition at line 87 of file FlavorHistoryEvent.h.

Constructor & Destructor Documentation

reco::FlavorHistoryEvent::FlavorHistoryEvent ( )

inline

Definition at line 100 of file FlavorHistoryEvent.h.

References clear().

{ clear(); }

reco::FlavorHistoryEvent::~FlavorHistoryEvent ( )

inline

Definition at line 101 of file FlavorHistoryEvent.h.

{}

Member Function Documentation

const_iterator reco::FlavorHistoryEvent::begin ( void  ) const

inline

Definition at line 144 of file FlavorHistoryEvent.h.

References histories_.

{ return histories_.begin(); }

void FlavorHistoryEvent::cache

(

)

Definition at line 11 of file FlavorHistoryEvent.cc.

References funct::abs(), filterCSVwithJSON::copy, gather_cfg::cout, reco::helpers::FlavorHistoryEventHelper::dR, reco::LeafCandidate::energy(), reco::helpers::FlavorHistoryEventHelper::flavor, reco::FlavorHistory::FLAVOR_NULL, reco::helpers::FlavorHistoryEventHelper::flavorSource, reco::FlavorHistory::flavorSource(), mps_fire::i, edm::Ptr< T >::isNonnull(), reco::ShallowClonePtrCandidate::masterClonePtr(), reco::FlavorHistory::matchedJet(), reco::FlavorHistory::parton(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), and verbose.

                               {
  bool verbose = false;

  if (verbose)
    cout << "----- Caching Flavor History Event  -----" << endl;
  // Set cached to false
  cached_ = false;
  nb_ = nc_ = 0;
  highestFlavor_ = 0;
  dR_ = 0.0;
  flavorSource_ = FlavorHistory::FLAVOR_NULL;

  // get list of flavor --> type --> dR.
  // Will sort this later to determine event classification
  vector<helpers::FlavorHistoryEventHelper> classification;

  // get iterators to the history vector
  const_iterator i = histories_.begin(), ibegin = histories_.begin(), iend = histories_.end();
  // loop over the history vector and count the number of
  // partons of flavors "b" and "c" that have a matched genjet.
  for (; i != iend; ++i) {
    FlavorHistory const& flavHist = *i;
    if (verbose)
      cout << "  Processing flavor history: " << i - ibegin << " = " << endl << flavHist << endl;
    CandidatePtr const& parton = flavHist.parton();
    flavor_type flavorSource = flavHist.flavorSource();

    // Now examine the matched jets to see what the classification should be.
    int pdgId = -1;
    if (parton.isNonnull())
      pdgId = abs(parton->pdgId());
    ShallowClonePtrCandidate const& matchedJet = flavHist.matchedJet();
    // Only count events with a matched genjet
    if (matchedJet.masterClonePtr().isNonnull()) {
      TLorentzVector p41(matchedJet.px(), matchedJet.py(), matchedJet.pz(), matchedJet.energy());
      if (pdgId == 5)
        nb_++;
      if (pdgId == 4)
        nc_++;

      // Get the sister genjet
      ShallowClonePtrCandidate const& sisterJet = i->sisterJet();
      TLorentzVector p42(sisterJet.px(), sisterJet.py(), sisterJet.pz(), sisterJet.energy());

      // Now check the source.
      double dR = -1;
      if (sisterJet.masterClonePtr().isNonnull()) {
        dR = p41.DeltaR(p42);
      }
      // Add to the vector to be sorted later
      if (verbose)
        cout << "Adding classification: pdgId = " << pdgId << ", flavorSource = " << flavorSource << ", dR = " << dR
             << endl;
      classification.push_back(helpers::FlavorHistoryEventHelper(pdgId, flavorSource, dR));
    } else {
      if (verbose)
        cout << "No matched jet found, not adding to classification list" << endl;
    }
  }

  // Sort by priority

  // Priority is:
  //
  //  1. flavor (5 > 4)
  //  2. type:
  //      2a. Flavor decay
  //      2b. Matrix element
  //      2c. Flavor excitation
  //      2d. Gluon splitting
  //  3. delta R (if applicable)
  if (!classification.empty()) {
    std::sort(classification.begin(), classification.end());

    if (verbose) {
      cout << "Writing out list of classifications" << endl;
      copy(classification.begin(), classification.end(), ostream_iterator<helpers::FlavorHistoryEventHelper>(cout, ""));
    }

    helpers::FlavorHistoryEventHelper const& best = *(classification.rbegin());
    dR_ = best.dR;
    highestFlavor_ = best.flavor;
    flavorSource_ = best.flavorSource;
  } else {
    dR_ = -1.0;
    highestFlavor_ = 0;
    flavorSource_ = FlavorHistory::FLAVOR_NULL;
  }

  // now we're cached, can return values quickly
  cached_ = true;
}

void reco::FlavorHistoryEvent::clear ( void  )

inline

Definition at line 158 of file FlavorHistoryEvent.h.

References cached_, dR_, reco::FlavorHistory::FLAVOR_NULL, flavorSource_, highestFlavor_, histories_, nb_, and nc_.

Referenced by FlavorHistoryEvent(), and BeautifulSoup.Tag::setString().

                 {
      cached_ = false;
      histories_.clear();
      nb_ = nc_ = 0;
      dR_ = 0.0;
      highestFlavor_ = 0;
      flavorSource_ = FlavorHistory::FLAVOR_NULL;
    }

double reco::FlavorHistoryEvent::deltaR ( ) const

inline

Definition at line 122 of file FlavorHistoryEvent.h.

References dR_, and isCached().

                          {
      if (isCached())
        return dR_;
      else
        return -1.0;
    }

const_iterator reco::FlavorHistoryEvent::end ( void  ) const

inline

Definition at line 145 of file FlavorHistoryEvent.h.

References histories_.

Referenced by Types.LuminosityBlockRange::cppID(), and Types.EventRange::cppID().

{ return histories_.end(); }

flavor_type reco::FlavorHistoryEvent::flavorSource ( ) const

inline

Definition at line 134 of file FlavorHistoryEvent.h.

References reco::FlavorHistory::FLAVOR_NULL, flavorSource_, and isCached().

                                     {
      if (isCached())
        return flavorSource_;
      else
        return FlavorHistory::FLAVOR_NULL;
    }

unsigned int reco::FlavorHistoryEvent::highestFlavor ( ) const

inline

Definition at line 128 of file FlavorHistoryEvent.h.

References highestFlavor_, and isCached().

                                       {
      if (isCached())
        return highestFlavor_;
      else
        return 0;
    }

bool reco::FlavorHistoryEvent::isCached ( ) const

inline

Definition at line 105 of file FlavorHistoryEvent.h.

References cached_.

Referenced by deltaR(), flavorSource(), highestFlavor(), nb(), and nc().

{ return cached_; }

unsigned int reco::FlavorHistoryEvent::nb ( ) const

inline

Definition at line 108 of file FlavorHistoryEvent.h.

References isCached(), and nb_.

                            {
      if (isCached())
        return nb_;
      else
        return 0;
    }

unsigned int reco::FlavorHistoryEvent::nc ( ) const

inline

Definition at line 114 of file FlavorHistoryEvent.h.

References isCached(), and nc_.

                            {
      if (isCached())
        return nc_;
      else
        return 0;
    }

void reco::FlavorHistoryEvent::push_back ( const value_type &  v )

inline

Definition at line 150 of file FlavorHistoryEvent.h.

References cached_, and histories_.

                                        {
      cached_ = false;
      histories_.push_back(v);
    }

const_reverse_iterator reco::FlavorHistoryEvent::rbegin ( ) const

inline

Definition at line 146 of file FlavorHistoryEvent.h.

References histories_.

{ return histories_.rbegin(); }

const_reverse_iterator reco::FlavorHistoryEvent::rend ( ) const

inline

Definition at line 147 of file FlavorHistoryEvent.h.

References histories_.

{ return histories_.rend(); }

void reco::FlavorHistoryEvent::resize ( size_t  n )

inline

Definition at line 154 of file FlavorHistoryEvent.h.

References cached_, and histories_.

                          {
      cached_ = false;
      histories_.resize(n);
    }

size_type reco::FlavorHistoryEvent::size ( void  ) const

inline

Definition at line 143 of file FlavorHistoryEvent.h.

References histories_.

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

{ return histories_.size(); }

Member Data Documentation

bool reco::FlavorHistoryEvent::cached_

protected

Definition at line 169 of file FlavorHistoryEvent.h.

Referenced by clear(), isCached(), push_back(), and resize().

double reco::FlavorHistoryEvent::dR_

protected

Definition at line 172 of file FlavorHistoryEvent.h.

Referenced by clear(), and deltaR().

flavor_type reco::FlavorHistoryEvent::flavorSource_

protected

Definition at line 174 of file FlavorHistoryEvent.h.

Referenced by clear(), and flavorSource().

unsigned int reco::FlavorHistoryEvent::highestFlavor_

protected

Definition at line 173 of file FlavorHistoryEvent.h.

Referenced by clear(), and highestFlavor().

collection_type reco::FlavorHistoryEvent::histories_

protected

Definition at line 168 of file FlavorHistoryEvent.h.

Referenced by begin(), clear(), end(), push_back(), rbegin(), rend(), resize(), and size().

unsigned int reco::FlavorHistoryEvent::nb_

protected

Definition at line 170 of file FlavorHistoryEvent.h.

Referenced by clear(), and nb().

unsigned int reco::FlavorHistoryEvent::nc_

protected

Definition at line 171 of file FlavorHistoryEvent.h.

Referenced by clear(), and nc().