Megajet.h

Go to the documentation of this file.

//-------------------------------------------------------
// Description:  Razor Megajet calculator
// Authors:      Maurizio Pierini, CERN
// Authors:      Christopher Rogan, Javier Duarte, Caltech
// Authors:      Emanuele Di Marco, CERN
//-------------------------------------------------------


#ifndef Megajet_h
#define Megajet_h

#include <vector>
#include <iostream>
#include <TLorentzVector.h>

namespace heppy {

  class Megajet {
  public:
    // There are 2 constructors:
    // 1. Constructor taking as argument vectors of Px, Py, Pz and E of the objects in the event and the hemisphere
    // association method,
    // 2. Constructor taking as argument vectors of Px, Py, Pz and E of the objects in the event.
    // The hemisphere association method should then be defined by SetMethod(association_method).
    //
    // Hemisphere association method:
    //   1: minimum sum of the invariant masses of the two hemispheres
    //   2: minimum difference of HT for the two hemispheres
    //   3: minimum m1^2/E1 + m2^2/E2
    //   4: Georgi distance: maximum (E1-Beta*m1^2/E1 + E2-Beta*m1^2/E2)

    Megajet(){};

    Megajet(std::vector<float> Px_vector,
            std::vector<float> Py_vector,
            std::vector<float> Pz_vector,
            std::vector<float> E_vector,
            int megajet_association_method);

    Megajet(std::vector<float> Px_vector,
            std::vector<float> Py_vector,
            std::vector<float> Pz_vector,
            std::vector<float> E_vector);

    ~Megajet(){};

    // return Nx, Ny, Nz, P, E of the axis of group 1
    std::vector<float> getAxis1();
    // return Nx, Ny, Nz, P, E of the axis of group 2
    std::vector<float> getAxis2();

    // where Nx, Ny, Nz are the direction cosines e.g. Nx = Px/P,
    // P is the momentum, E is the energy

    // set or overwrite the seed and association methods
    void SetMethod(int megajet_association_method) {
      megajet_meth = megajet_association_method;
      status = 0;
    }

  private:
    void CombineMinMass();

    void CombineMinHT();

    void CombineMinEnergyMass();

    void CombineGeorgi();

    void Combine();

    std::vector<float> Object_Px;
    std::vector<float> Object_Py;
    std::vector<float> Object_Pz;
    std::vector<float> Object_E;

    std::vector<float> Axis1;
    std::vector<float> Axis2;

    int megajet_meth;
    int status;

    std::vector<TLorentzVector> jIN;
    std::vector<TLorentzVector> jOUT;
    std::vector<TLorentzVector> j1;
    std::vector<TLorentzVector> j2;
  };
}  // namespace heppy

#endif