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;

};
}

#endif