AnaTools/src/StopCTauWeight.cc

#include "OSUT3Analysis/AnaTools/interface/StopCTauWeight.h"

StopCTauWeight::StopCTauWeight (double currentCTau, double targetCTau) :
  currentCTau_ (currentCTau),
  targetCTau_ (targetCTau)
{
}

StopCTauWeight::StopCTauWeight (double currentCTau, double targetCTau, edm::InputTag stops) :
  currentCTau_ (currentCTau),
  targetCTau_ (targetCTau),
  stops_ (stops)
{
}

StopCTauWeight::~StopCTauWeight ()
{
}

double
StopCTauWeight::operator[] (const BNstopCollection &stops)
{
  double weight, ctau0, ctau1;

  if (stops.size () != 2)
    cout << "Wrong number of stops!: " << stops.size () << endl;
  ctau0 = stops.at (0).ctau;
  ctau1 = stops.at (1).ctau;

  weight = (currentCTau_ / targetCTau_) * (currentCTau_ / targetCTau_);
  weight *= exp (-(ctau0 + ctau1) * ((1 / targetCTau_) - (1 / currentCTau_)));

  return weight;
}

double
StopCTauWeight::operator[] (const edm::Event &event)
{
  edm::Handle<BNstopCollection> stops;
  event.getByLabel (stops_, stops);

  return (*this)[*stops.product ()];
}
Revision:	1.1
Committed:	Fri May 3 20:09:35 2013 UTC (12 years ago) by ahart
Content type:	text/plain
Branch:	MAIN
CVS Tags:	V02-03-02, V02-03-01, V02-03-00, V02-02-00, V02-01-01, V02-01-00, V01-01-00, V01-00-01, V01-00-00, V00-01-00, HEAD
Log Message:	Moved here from DisplacedSUSY/Configuration/src/CTauWeight.cc so we can get rid of the DISPLACED_SUSY macro.
#	User	Rev	Content
1	ahart	1.1	#include "OSUT3Analysis/AnaTools/interface/StopCTauWeight.h"
2
3			StopCTauWeight::StopCTauWeight (double currentCTau, double targetCTau) :
4			currentCTau_ (currentCTau),
5			targetCTau_ (targetCTau)
6			{
7			}
8
9			StopCTauWeight::StopCTauWeight (double currentCTau, double targetCTau, edm::InputTag stops) :
10			currentCTau_ (currentCTau),
11			targetCTau_ (targetCTau),
12			stops_ (stops)
13			{
14			}
15
16			StopCTauWeight::~StopCTauWeight ()
17			{
18			}
19
20			double
21			StopCTauWeight::operator[] (const BNstopCollection &stops)
22			{
23			double weight, ctau0, ctau1;
24
25			if (stops.size () != 2)
26			cout << "Wrong number of stops!: " << stops.size () << endl;
27			ctau0 = stops.at (0).ctau;
28			ctau1 = stops.at (1).ctau;
29
30			weight = (currentCTau_ / targetCTau_) * (currentCTau_ / targetCTau_);
31			weight = exp (-(ctau0 + ctau1) ((1 / targetCTau_) - (1 / currentCTau_)));
32
33			return weight;
34			}
35
36			double
37			StopCTauWeight::operator[] (const edm::Event &event)
38			{
39			edm::Handle<BNstopCollection> stops;
40			event.getByLabel (stops_, stops);
41
42			return (this)[stops.product ()];
43			}