#include <FlavorHistorySelectorUtil.h>

Public Types
typedef reco::FlavorHistory::FLAVOR_T	flavor_type

typedef std::vector< int >	flavor_vector

Public Member Functions
	FlavorHistorySelectorUtil (unsigned int flavor, unsigned int noutput, flavor_vector const &flavorSource, double minDR, double maxDR, bool verbose)

bool	select (unsigned int nb, unsigned int nc, unsigned int highestFlavor, FlavorHistory::FLAVOR_T flavorSource, double dr) const

	~FlavorHistorySelectorUtil ()

Private Attributes
int	flavor_

flavor_vector	flavorSource_

double	maxDR_

double	minDR_

int	noutput_

bool	verbose_

Detailed Description

Definition at line 54 of file FlavorHistorySelectorUtil.h.

Member Typedef Documentation

typedef reco::FlavorHistory::FLAVOR_T reco::FlavorHistorySelectorUtil::flavor_type

Definition at line 56 of file FlavorHistorySelectorUtil.h.

typedef std::vector<int> reco::FlavorHistorySelectorUtil::flavor_vector

Definition at line 57 of file FlavorHistorySelectorUtil.h.

Constructor & Destructor Documentation

FlavorHistorySelectorUtil::FlavorHistorySelectorUtil	(	unsigned int	flavor,
		unsigned int	noutput,
		flavor_vector const &	flavorSource,
		double	minDR,
		double	maxDR,
		bool	verbose
	)

Definition at line 8 of file FlavorHistorySelectorUtil.cc.

References maxDR_, and minDR_.

                                              :
   flavor_(flavor),
   noutput_(noutput),
   flavorSource_(flavorSource),
   minDR_(minDR),
   maxDR_(maxDR),
   verbose_(verbose)
 {
   
   // Deal with the case if minDR == maxDR, just increment maxDR by epsilon
   if ( maxDR_ - minDR_ <= 0.001 ) maxDR_ += 0.001;
 }

reco::FlavorHistorySelectorUtil::~FlavorHistorySelectorUtil ( )

inline

Definition at line 65 of file FlavorHistorySelectorUtil.h.

65 {}

Member Function Documentation

bool FlavorHistorySelectorUtil::select	(	unsigned int	nb,
		unsigned int	nc,
		unsigned int	highestFlavor,
		FlavorHistory::FLAVOR_T	flavorSource,
		double	dr
	)		const

Definition at line 26 of file FlavorHistorySelectorUtil.cc.

References reco::FlavorHistory::bQuarkId, gather_cfg::cout, reco::FlavorHistory::cQuarkId, spr::find(), flavor_, flavorSource_, maxDR_, minDR_, noutput_, and verbose_.

Referenced by FlavorHistoryFilter::filter().

 {
   
   // Print out some information about this event
   if ( verbose_ ) {
     cout << "Looking at flavor history event: " << endl;
     cout << "source   = " << flavorSource << endl;
     cout << "nbjet    = " << nb << endl;
     cout << "ncjet    = " << nc << endl;
     cout << "flavor   = " << highestFlavor << endl;
     cout << "dr       = " << dr << endl;
   }
 
   // First check that the highest flavor in the event is what this
   // filter is checking. Otherwise we need to fail the event,
   // since it should be handled by another filter
   if ( highestFlavor > static_cast<unsigned int>(flavor_) ) {
     if ( verbose_ ) cout << "Rejecting event, highest flavor is " << highestFlavor << endl;
     return false;
   }
 
   // Next check that the flavor source is one of the desired ones
   vector<int>::const_iterator iflavorSource = find( flavorSource_.begin(), flavorSource_.end(), static_cast<int>(flavorSource) );
   if ( iflavorSource == flavorSource_.end() ) {
     if ( verbose_ ) cout << "Rejecting event, didn't find flavor source " << static_cast<int>(flavorSource) << endl;
     return false;
   }
   
   // If we are examining b quarks
   if ( flavor_ == reco::FlavorHistory::bQuarkId ) {
     // if we have no b quarks, return false
     if ( nb <= 0 ) {
       if ( verbose_ ) cout << "Rejecting event, nb = " << nb << endl;
       return false;
     }
     // here, nb > 0
     else {
       // if we want 1 b, require nb == 1
       if ( noutput_ == 1 && nb == 1 ) {
     if ( verbose_ ) cout << "Accepting event" << endl;
     return true;
       }
       // if we want 2 b, then look at delta R
       else if ( noutput_ > 1 && nb > 1 ) {
     // If dr is within the range we want, pass.
     // Otherwise, fail.
     if ( verbose_ ) cout << "Want multiples, dr = " << dr << endl;
     return ( dr >= minDR_ && dr < maxDR_ );
       }
       // otherwise return false
       else {
     if ( verbose_ ) cout << "Rejecting event, isn't output = 1 + nb = 1, or output > 0 and delta R in proper range" << endl;
     return false;
       }
     }// end if nb > 0
     
   } // end if flavor is b quark
 
   // If we are examining c quarks
   else if ( flavor_ == reco::FlavorHistory::cQuarkId ) {
     // make sure there are no b quarks in the event.
     // If there are, another filter should handle it.
     if ( nb > 0 ) return false;
     
     // if we have no c quarks, return false
     if ( nc <= 0 ) return false;
     // here, nc > 0
     else {
       // if we want 1 c, require nc == 1
       if ( noutput_ == 1 && nc == 1 ) {
     return true;
       }
       // if we want 2 c, then look at delta R
       else if ( noutput_ > 1 && nc > 1 ) {
     // If dr is within the range we want, pass.
     // Otherwise, fail.
     return ( dr >= minDR_ && dr < maxDR_ );
       }
       // otherwise return false
       else {
     return false;
       }
     }// end if nc > 0
     
   }
   // Otherwise return false
   else {
     if ( verbose_ ) cout << "Something is weird, flavor is " << flavor_ << endl;
     return false;
   }
 }