Not complete. More...

#include <NeutrinoWeighting.h>

Inheritance diagram for top::NeutrinoWeighting:

Collaboration diagram for top::NeutrinoWeighting:

Public Member Functions
	NeutrinoWeighting ()

virtual	~NeutrinoWeighting ()

bool	apply (const top::Event &) const override
	This does stuff based on the information in an event. More...

std::string	name () const override
	A human readable name. More...

virtual bool	applyParticleLevel (const top::ParticleLevelEvent &) const
	This does stuff based on the information in a particle level event. More...

Private Member Functions
NWSolution	solveForNeutrinoEta (const TLorentzVector &lepton, const TLorentzVector &bJet, double topMass, int index) const

double	neutrino_weight (const TLorentzVector &, const TLorentzVector &, double, double) const
	Calculate the weight for this combination of particles by comparing with Met. More...

Private Attributes
double	neutrinos [2000][3]
	Sampling points. More...

int	etaSize
	Delta eta sampling size. For the mass we don't need to sample so many points. More...

double	m_bmass
	B jet mass constraint. More...

double	m_wmass
	W jet mass constraint. More...

double	sigmax

double	sigmay

Detailed Description

Not complete.

Work in progress

Definition at line 57 of file NeutrinoWeighting.h.

Constructor & Destructor Documentation

◆ NeutrinoWeighting()

top::NeutrinoWeighting::NeutrinoWeighting ( )

Definition at line 12 of file NeutrinoWeighting.cxx.

                                        : sigmax(10), sigmay(10) {
     double etaStep = 0.2;
  
     //construct the array of eta values to use (cosh and sinh them)
     int index = 0;
  
     for (double eta = -5.0; eta < 5.0001; eta += etaStep) {
       neutrinos[index][0] = eta;
       neutrinos[index][1] = std::sinh(eta);
       neutrinos[index][2] = std::cosh(eta);
       ++index;
     }
  
     etaSize = index;
  
     if (etaSize > 2000) {
       std::cout << "ERROR TOO MANY SAMPLING POINT FOR neutrinos ARRAY" << std::endl;
       exit(1);
     }
   }

◆ ~NeutrinoWeighting()

top::NeutrinoWeighting::~NeutrinoWeighting ( )

virtual

Definition at line 33 of file NeutrinoWeighting.cxx.

33 {

34 }

Member Function Documentation

◆ apply()

bool top::NeutrinoWeighting::apply ( const top::Event & ) const

overridevirtual

This does stuff based on the information in an event.

The idea is that you implement this to return either true or false based on the information held within top::Event. If this returns true then the event is kept. If it returns false then the event is removed.

Parameters

top::Event The current event.

Returns: true if the event should be kept, false otherwise.

Implements top::EventSelectorBase.

Definition at line 36 of file NeutrinoWeighting.cxx.

                                                            {
     TLorentzVector l1;
     TLorentzVector l2;
     TLorentzVector b1;
     TLorentzVector b2;
     double topMass = 172500.;
  
     const double met_ex = event.m_met->mpx();
     const double met_ey = event.m_met->mpy();
  
     for (int i = 0; i < etaSize; ++i) {
       for (int j = 0; j < etaSize; ++j) {
         NWSolution ans1 = solveForNeutrinoEta(l1, b1, topMass, i);
         NWSolution ans2 = solveForNeutrinoEta(l2, b2, topMass, j);
  
         if (ans1.getNumSolutions() > 0 && ans2.getNumSolutions() > 0) {
           const double sol1 = neutrino_weight(ans1.getv1(), ans2.getv1(), met_ex, met_ey);
           const double sol2 = neutrino_weight(ans1.getv1(), ans2.getv2(), met_ex, met_ey);
           const double sol3 = neutrino_weight(ans1.getv2(), ans2.getv1(), met_ex, met_ey);
           const double sol4 = neutrino_weight(ans1.getv2(), ans2.getv2(), met_ex, met_ey);
  
           if (std::isnan(sol1) || std::isnan(sol2) || std::isnan(sol3) ||
               std::isnan(sol4)) std::cout << "One of the quadratic ans is NaN!" << std::endl;
         }
  
         ans1 = solveForNeutrinoEta(l1, b2, topMass, i);
         ans2 = solveForNeutrinoEta(l2, b1, topMass, j);
  
         if (ans1.getNumSolutions() > 0 && ans2.getNumSolutions() > 0) {
           const double sol1 = neutrino_weight(ans1.getv1(), ans2.getv1(), met_ex, met_ey);
           const double sol2 = neutrino_weight(ans1.getv1(), ans2.getv2(), met_ex, met_ey);
           const double sol3 = neutrino_weight(ans1.getv2(), ans2.getv1(), met_ex, met_ey);
           const double sol4 = neutrino_weight(ans1.getv2(), ans2.getv2(), met_ex, met_ey);
  
           if (std::isnan(sol1) || std::isnan(sol2) || std::isnan(sol3) ||
               std::isnan(sol4)) std::cout << "One of the quadratic ans is NaN!" << std::endl;
         }
       }
     }
  
     //normalise?
  
     return true;
   }

◆ applyParticleLevel()

virtual bool top::EventSelectorBase::applyParticleLevel ( const top::ParticleLevelEvent & ) const

inlinevirtualinherited

This does stuff based on the information in a particle level event.

The idea is that you implement this to return either true or false, based on the information held within the top::ParticleLevelEvent. If this function returns true, then the event is kept, otherwise it is removed. The function has a default implementation (which returns true) because it is expected that many EventSelector objects do not operate on ParticleLevelEvent objects.

Parameters

top::ParticleLevelEvent	the current particle level event.
true	if the event should be kept (i.e. it passed the selector criteria), false otherwise.

Reimplemented in top::JetNGhostSelector, top::PrintEventSelector, top::PseudoTopRecoRun, top::NElectronNMuonTightSelector, top::NElectronNMuonSelector, top::NFwdElectronSelector, top::HTSelector, top::OSLeptonTightSelector, top::MLLSelector, top::MWTSelector, top::NElectronTightSelector, top::NFwdElectronTightSelector, top::NMuonTightSelector, top::OSLeptonSelector, top::METMWTSelector, top::METSelector, top::MLLWindow, top::NElectronSelector, top::NJetSelector, top::NMuonSelector, top::NPhotonSelector, top::NSoftMuonSelector, top::NTauSelector, top::SSLeptonTightSelector, top::SSLeptonSelector, top::ParticleLevelSelector, top::RecoLevelSelector, top::NVarRCJetSelector, top::NLargeJetSelector, and top::NRCJetSelector.

Definition at line 73 of file EventSelectorBase.h.

73 {return true;}

◆ name()

std::string top::NeutrinoWeighting::name ( ) const

inlineoverridevirtual

A human readable name.

Mostly used for printing the cut and value to the screen. Must be implemented for each tool.

Implements top::EventSelectorBase.

Definition at line 64 of file NeutrinoWeighting.h.

64 {return "RECO:NEUTRINOWEIGHTING";}

◆ neutrino_weight()

double top::NeutrinoWeighting::neutrino_weight	(	const TLorentzVector &	neutrino1,
		const TLorentzVector &	neutrino2,
		double	met_ex,
		double	met_ey
	)		const

private

Calculate the weight for this combination of particles by comparing with Met.

Weight a solution based on the agreement with missing Et.

Parameters

neutrino1	The first neutrino (actually the order doesn't matter).
neutrino2	The second neutrino.
met_ex	The x component of missing Et.
met_ey	The y component of missing Et.

Returns: The weight of this solution based on the input neutrinos.

Definition at line 124 of file NeutrinoWeighting.cxx.

                                                                                 {
     const double dx = met_ex - neutrino1.Px() - neutrino2.Px();
     const double dy = met_ey - neutrino1.Py() - neutrino2.Py();
  
     return exp(-dx * dx / (2. * sigmax * sigmax) - dy * dy / (2. * sigmay * sigmay));
   }

◆ solveForNeutrinoEta()

NWSolution top::NeutrinoWeighting::solveForNeutrinoEta	(	const TLorentzVector &	lepton,
		const TLorentzVector &	bJet,
		double	topMass,
		int	index
	)		const

private

Parameters

lepton	The lepton 4 vector.
bJet	The b jet 4 vector.
topMass	The top mass is a constraint in the spin correlation calculation.
index	Gives the position in the neutrino eta array of this neutrino.

Definition at line 81 of file NeutrinoWeighting.cxx.

                                                                                      {
     double Wmass2 = m_wmass * m_wmass;
     double bmass = m_bmass;
  
     double Elprime = lepton.E() * neutrinos[index][2] - lepton.Pz() * neutrinos[index][1];
     double Ebprime = bJet.E() * neutrinos[index][2] - bJet.Pz() * neutrinos[index][1];
  
     double A = (lepton.Py() * Ebprime - bJet.Py() * Elprime) / (bJet.Px() * Elprime - lepton.Px() * Ebprime);
     double B = (Elprime * (topMass * topMass - Wmass2 - bmass * bmass - 2. * lepton * bJet) - Ebprime * Wmass2) /
                (2. * (lepton.Px() * Ebprime - bJet.Px() * Elprime));
  
     double par1 = (lepton.Px() * A + lepton.Py()) / Elprime;
     double C = A * A + 1. - par1 * par1;
     double par2 = (Wmass2 / 2. + lepton.Px() * B) / Elprime;
     double D = 2. * (A * B - par2 * par1);
     double F = B * B - par2 * par2;
     double det = D * D - 4. * C * F;
  
     NWSolution sol;
  
     sol.setSolutions(0);
  
     if (det > 0.) {
       double tmp = std::sqrt(det) / (2. * C);
       double py1 = -D / (2. * C) + tmp;
       double py2 = -D / (2. * C) - tmp;
       double px1 = A * py1 + B;
       double px2 = A * py2 + B;
       double pT2_1 = px1 * px1 + py1 * py1;
       double pT2_2 = px2 * px2 + py2 * py2;
       double pz1 = std::sqrt(pT2_1) * neutrinos[index][1];
       double pz2 = std::sqrt(pT2_2) * neutrinos[index][1];
  
       TLorentzVector a1(px1, py1, pz1, sqrt(pT2_1 + pz1* pz1));
       TLorentzVector a2(px2, py2, pz2, sqrt(pT2_2 + pz2* pz2));
  
       sol.setSolutions(2, a1, a2);
     }
  
     return sol;
   }

Member Data Documentation

◆ etaSize

int top::NeutrinoWeighting::etaSize

private

Delta eta sampling size. For the mass we don't need to sample so many points.

Definition at line 93 of file NeutrinoWeighting.h.

◆ m_bmass

double top::NeutrinoWeighting::m_bmass

private

B jet mass constraint.

Definition at line 100 of file NeutrinoWeighting.h.

◆ m_wmass

double top::NeutrinoWeighting::m_wmass

private

W jet mass constraint.

Definition at line 103 of file NeutrinoWeighting.h.

◆ neutrinos

double top::NeutrinoWeighting::neutrinos[2000][3]

private

Sampling points.

Definition at line 90 of file NeutrinoWeighting.h.

◆ sigmax

double top::NeutrinoWeighting::sigmax

private

Definition at line 106 of file NeutrinoWeighting.h.

◆ sigmay

double top::NeutrinoWeighting::sigmay

private

Definition at line 109 of file NeutrinoWeighting.h.

The documentation for this class was generated from the following files:

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ NeutrinoWeighting()

◆ ~NeutrinoWeighting()

Member Function Documentation

◆ apply()

◆ applyParticleLevel()

◆ name()

◆ neutrino_weight()

◆ solveForNeutrinoEta()

Member Data Documentation

◆ etaSize

◆ m_bmass

◆ m_wmass

◆ neutrinos

◆ sigmax

◆ sigmay