NeutrinoWeighting.h

Go to the documentation of this file.

/*
   Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
 */
 
#ifndef NEUTRINOWEIGHTING_H_
#define NEUTRINOWEIGHTING_H_
 
#include "TopEventSelectionTools/EventSelectorBase.h"
 
namespace top {
  class Event;
 
  class NWSolution {
  public:
    void setSolutions(int num, const TLorentzVector& a, const TLorentzVector& b) {
      m_v1 = a;
      m_v2 = b;
      m_solutions = num;
    }
 
    void setSolutions(int num) {
      m_solutions = num;
    }
 
    int getNumSolutions() const {
      return m_solutions;
    }
 
    TLorentzVector getv1() const {
      return m_v1;
    }
 
    TLorentzVector getv2() const {
      return m_v2;
    }
 
  private:
    int m_solutions;
 
    TLorentzVector m_v1;
 
    TLorentzVector m_v2;
  };
 
  class NeutrinoWeighting: public EventSelectorBase {
  public:
    NeutrinoWeighting();
    virtual ~NeutrinoWeighting();
 
    bool apply(const top::Event&) const override;
 
    std::string name() const override {return "RECO:NEUTRINOWEIGHTING";}
  private:
    NWSolution solveForNeutrinoEta(const TLorentzVector& lepton, const TLorentzVector& bJet, double topMass,
                                   int index) const;
 
    double neutrino_weight(const TLorentzVector&, const TLorentzVector&, double, double) const;
 
    double neutrinos[2000][3];
 
    int etaSize;
 
    // commented out variables are unused.  experts please check and remove
    // ///Top Quark Mass constraint.
    // double topMass;
 
    double m_bmass;
 
    double m_wmass;
 
    //For weighter.
    double sigmax;
 
    //For weighter.
    double sigmay;
  };
}
 
#endif