GenFilterInfo Class Reference

#include <GenFilterInfo.h>

Public Member Functions
double	filterEfficiency (int idwtup=+3) const
double	filterEfficiencyError (int idwtup=+3) const
	GenFilterInfo ()
	GenFilterInfo (unsigned int, unsigned int)
	GenFilterInfo (unsigned int, unsigned int, unsigned int, unsigned int, double, double, double, double)
	GenFilterInfo (const GenFilterInfo &)
bool	mergeProduct (GenFilterInfo const &other)
unsigned int	numEventsPassed () const
unsigned int	numEventsTotal () const
unsigned int	numEventsTried () const
unsigned int	numPassNegativeEvents () const
unsigned int	numPassPositiveEvents () const
unsigned int	numTotalNegativeEvents () const
unsigned int	numTotalPositiveEvents () const
double	sumFailWeights () const
double	sumFailWeights2 () const
double	sumPassWeights () const
double	sumPassWeights2 () const
double	sumWeights () const
double	sumWeights2 () const
virtual	~GenFilterInfo ()

Private Attributes
unsigned int	numPassNegativeEvents_
unsigned int	numPassPositiveEvents_
unsigned int	numTotalNegativeEvents_
unsigned int	numTotalPositiveEvents_
double	sumPassWeights2_
double	sumPassWeights_
double	sumTotalWeights2_
double	sumTotalWeights_

Detailed Description

Definition at line 11 of file GenFilterInfo.h.

Constructor & Destructor Documentation

GenFilterInfo::GenFilterInfo

(

)

Definition at line 11 of file GenFilterInfo.cc.

                             :
  numPassPositiveEvents_(0),
  numPassNegativeEvents_(0),
  numTotalPositiveEvents_(0),
  numTotalNegativeEvents_(0),
  sumPassWeights_(0.),
  sumPassWeights2_(0.),
  sumTotalWeights_(0.),
  sumTotalWeights2_(0.)
{
}

GenFilterInfo::GenFilterInfo	(	unsigned int	tried,
		unsigned int	pass
	)

Definition at line 23 of file GenFilterInfo.cc.

                                                                  :
  numPassPositiveEvents_(pass),
  numPassNegativeEvents_(0),
  numTotalPositiveEvents_(tried),
  numTotalNegativeEvents_(0),
  sumPassWeights_(pass),
  sumPassWeights2_(pass),
  sumTotalWeights_(tried),
  sumTotalWeights2_(tried)
{
}

GenFilterInfo::GenFilterInfo	(	unsigned int	passp,
		unsigned int	passn,
		unsigned int	totalp,
		unsigned int	totaln,
		double	passw,
		double	passw2,
		double	totalw,
		double	totalw2
	)

Definition at line 35 of file GenFilterInfo.cc.

                                                                             :
  numPassPositiveEvents_(passp),
  numPassNegativeEvents_(passn),
  numTotalPositiveEvents_(totalp),
  numTotalNegativeEvents_(totaln),
  sumPassWeights_(passw),
  sumPassWeights2_(passw2),
  sumTotalWeights_(totalw),
  sumTotalWeights2_(totalw2)
{
}

GenFilterInfo::GenFilterInfo

(

const GenFilterInfo &

other

)

Definition at line 48 of file GenFilterInfo.cc.

                                                      :
  numPassPositiveEvents_(other.numPassPositiveEvents_),
  numPassNegativeEvents_(other.numPassNegativeEvents_),
  numTotalPositiveEvents_(other.numTotalPositiveEvents_),
  numTotalNegativeEvents_(other.numTotalNegativeEvents_),
  sumPassWeights_(other.sumPassWeights_),
  sumPassWeights2_(other.sumPassWeights2_),
  sumTotalWeights_(other.sumTotalWeights_),
  sumTotalWeights2_(other.sumTotalWeights2_)
{
}

GenFilterInfo::~GenFilterInfo ( )

virtual

Definition at line 60 of file GenFilterInfo.cc.

{
}

Member Function Documentation

double GenFilterInfo::filterEfficiency

(

int

idwtup = +3

)

const

Definition at line 82 of file GenFilterInfo.cc.

References alignCSCRings::e, numEventsPassed(), numEventsTotal(), sumPassWeights(), and sumWeights().

Referenced by GenFilterEfficiencyAnalyzer::endJob(), GenXSecAnalyzer::endJob(), and GenXSecAnalyzer::endRun().

                                                       {
  double eff = -1;
  switch(idwtup) {
  case 3: case -3:
    eff = numEventsTotal() > 0 ? (double) numEventsPassed() / (double) numEventsTotal(): -1;
    break;
  default:
    eff = sumWeights() > 1e-6? sumPassWeights()/sumWeights(): -1;
    break;
  }
  return eff;

}

double GenFilterInfo::filterEfficiencyError

(

int

idwtup = +3

)

const

Definition at line 96 of file GenFilterInfo.cc.

References cuy::denominator, alignCSCRings::e, cuy::numerator, numEventsTotal(), numPassNegativeEvents(), numPassPositiveEvents(), numTotalNegativeEvents(), numTotalPositiveEvents(), mathSSE::sqrt(), sumFailWeights(), sumFailWeights2(), sumPassWeights(), sumPassWeights2(), and sumWeights().

Referenced by GenFilterEfficiencyAnalyzer::endJob(), GenXSecAnalyzer::endJob(), and GenXSecAnalyzer::endRun().

                                                            {

  double efferr = -1;
  switch(idwtup) {
  case 3: case -3:
    {
      double effp  = numTotalPositiveEvents() > 1e-6? 
    (double)numPassPositiveEvents()/(double)numTotalPositiveEvents():0;
      double effp_err2 = numTotalPositiveEvents() > 1e-6?
    (1-effp)*effp/(double)numTotalPositiveEvents(): 0;

      double effn  = numTotalNegativeEvents() > 1e-6? 
    (double)numPassNegativeEvents()/(double)numTotalNegativeEvents():0;
      double effn_err2 = numTotalNegativeEvents() > 1e-6? 
    (1-effn)*effn/(double)numTotalNegativeEvents(): 0;

      efferr = numEventsTotal() > 0 ? sqrt ( 
                        ((double)numTotalPositiveEvents()*(double)numTotalPositiveEvents()*effp_err2 + 
                         (double)numTotalNegativeEvents()*(double)numTotalNegativeEvents()*effn_err2)
                        /(double)numEventsTotal()/(double)numEventsTotal()
                         ) : -1;
      break;
    }
  default:
    {
      double denominator = sumWeights()*sumWeights()*sumWeights()*sumWeights();
      double numerator =
    sumPassWeights2() * sumFailWeights()* sumFailWeights() +
    sumFailWeights2() * sumPassWeights()* sumPassWeights();
      efferr = denominator>1e-6? 
    sqrt(numerator/denominator):-1;
      break;
    }
  }

  return efferr;
}

bool GenFilterInfo::mergeProduct

(

GenFilterInfo const &

other

)

Definition at line 64 of file GenFilterInfo.cc.

References numPassNegativeEvents_, numPassPositiveEvents_, numTotalNegativeEvents_, numTotalPositiveEvents_, sumPassWeights2_, sumPassWeights_, sumTotalWeights2_, and sumTotalWeights_.

Referenced by GenFilterEfficiencyAnalyzer::endLuminosityBlock(), and GenXSecAnalyzer::endLuminosityBlock().

{
  // merging two GenFilterInfos means that the numerator and
  // denominator from the original product need to besummed with
  // those in the product we are going to merge

  numPassPositiveEvents_ += other.numPassPositiveEvents_;
  numPassNegativeEvents_ += other.numPassNegativeEvents_;
  numTotalPositiveEvents_ += other.numTotalPositiveEvents_;
  numTotalNegativeEvents_ += other.numTotalNegativeEvents_;
  sumPassWeights_        += other.sumPassWeights_;
  sumPassWeights2_       += other.sumPassWeights2_;
  sumTotalWeights_        += other.sumTotalWeights_;
  sumTotalWeights2_       += other.sumTotalWeights2_;

  return true;
}

unsigned int GenFilterInfo::numEventsPassed ( ) const

inline

Definition at line 24 of file GenFilterInfo.h.

References numPassNegativeEvents_, and numPassPositiveEvents_.

Referenced by GenXSecAnalyzer::endJob(), and filterEfficiency().

{ return fmax(0, (numPassPositiveEvents_ - numPassNegativeEvents_));}

unsigned int GenFilterInfo::numEventsTotal ( ) const

inline

Definition at line 25 of file GenFilterInfo.h.

References numTotalNegativeEvents_, and numTotalPositiveEvents_.

Referenced by GenXSecAnalyzer::endJob(), filterEfficiency(), and filterEfficiencyError().

{ return fmax(0, (numTotalPositiveEvents_ - numTotalNegativeEvents_));}

unsigned int GenFilterInfo::numEventsTried ( ) const

inline

Definition at line 23 of file GenFilterInfo.h.

References numTotalNegativeEvents_, and numTotalPositiveEvents_.

{ return (numTotalPositiveEvents_ + numTotalNegativeEvents_);}

unsigned int GenFilterInfo::numPassNegativeEvents ( ) const

inline

Definition at line 30 of file GenFilterInfo.h.

References numPassNegativeEvents_.

Referenced by GenXSecAnalyzer::endJob(), and filterEfficiencyError().

{ return numPassNegativeEvents_;}

unsigned int GenFilterInfo::numPassPositiveEvents ( ) const

inline

Definition at line 27 of file GenFilterInfo.h.

References numPassPositiveEvents_.

Referenced by GenXSecAnalyzer::endJob(), and filterEfficiencyError().

{ return numPassPositiveEvents_;}

unsigned int GenFilterInfo::numTotalNegativeEvents ( ) const

inline

Definition at line 31 of file GenFilterInfo.h.

References numTotalNegativeEvents_.

Referenced by GenXSecAnalyzer::endJob(), and filterEfficiencyError().

{ return numTotalNegativeEvents_;}

unsigned int GenFilterInfo::numTotalPositiveEvents ( ) const

inline

Definition at line 28 of file GenFilterInfo.h.

References numTotalPositiveEvents_.

Referenced by GenXSecAnalyzer::endJob(), and filterEfficiencyError().

{ return numTotalPositiveEvents_;}

double GenFilterInfo::sumFailWeights ( ) const

inline

Definition at line 37 of file GenFilterInfo.h.

References sumPassWeights_, and sumTotalWeights_.

Referenced by filterEfficiencyError().

{ return sumTotalWeights_ - sumPassWeights_;}

double GenFilterInfo::sumFailWeights2 ( ) const

inline

Definition at line 38 of file GenFilterInfo.h.

References sumPassWeights2_, and sumTotalWeights2_.

Referenced by filterEfficiencyError().

{ return sumTotalWeights2_ - sumPassWeights2_;}

double GenFilterInfo::sumPassWeights ( ) const

inline

Definition at line 34 of file GenFilterInfo.h.

References sumPassWeights_.

Referenced by GenFilterEfficiencyAnalyzer::endJob(), GenXSecAnalyzer::endJob(), filterEfficiency(), and filterEfficiencyError().

{ return sumPassWeights_;}

double GenFilterInfo::sumPassWeights2 ( ) const

inline

Definition at line 35 of file GenFilterInfo.h.

References sumPassWeights2_.

Referenced by filterEfficiencyError().

{ return sumPassWeights2_;}

double GenFilterInfo::sumWeights ( ) const

inline

Definition at line 40 of file GenFilterInfo.h.

References sumTotalWeights_.

Referenced by GenFilterEfficiencyAnalyzer::endJob(), GenXSecAnalyzer::endJob(), filterEfficiency(), and filterEfficiencyError().

{ return sumTotalWeights_;}

double GenFilterInfo::sumWeights2 ( ) const

inline

Definition at line 41 of file GenFilterInfo.h.

References sumTotalWeights2_.

Referenced by GenXSecAnalyzer::endJob(), and GenXSecAnalyzer::endRun().

{ return sumTotalWeights2_;}

Member Data Documentation

unsigned int GenFilterInfo::numPassNegativeEvents_

private

Definition at line 51 of file GenFilterInfo.h.

Referenced by mergeProduct(), numEventsPassed(), and numPassNegativeEvents().

unsigned int GenFilterInfo::numPassPositiveEvents_

private

Definition at line 50 of file GenFilterInfo.h.

Referenced by mergeProduct(), numEventsPassed(), and numPassPositiveEvents().

unsigned int GenFilterInfo::numTotalNegativeEvents_

private

Definition at line 53 of file GenFilterInfo.h.

Referenced by mergeProduct(), numEventsTotal(), numEventsTried(), and numTotalNegativeEvents().

unsigned int GenFilterInfo::numTotalPositiveEvents_

private

Definition at line 52 of file GenFilterInfo.h.

Referenced by mergeProduct(), numEventsTotal(), numEventsTried(), and numTotalPositiveEvents().

double GenFilterInfo::sumPassWeights2_

private

Definition at line 56 of file GenFilterInfo.h.

Referenced by mergeProduct(), sumFailWeights2(), and sumPassWeights2().

double GenFilterInfo::sumPassWeights_

private

Definition at line 55 of file GenFilterInfo.h.

Referenced by mergeProduct(), sumFailWeights(), and sumPassWeights().

double GenFilterInfo::sumTotalWeights2_

private

Definition at line 58 of file GenFilterInfo.h.

Referenced by mergeProduct(), sumFailWeights2(), and sumWeights2().

double GenFilterInfo::sumTotalWeights_

private

Definition at line 57 of file GenFilterInfo.h.

Referenced by mergeProduct(), sumFailWeights(), and sumWeights().