heppy::Megajet Class Reference

#include <Megajet.h>

Public Member Functions
std::vector< float >	getAxis1 ()
std::vector< float >	getAxis2 ()
	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)
void	SetMethod (int megajet_association_method)
	~Megajet ()
	Destructor. More...

Private Member Functions
void	Combine ()
	Combine the jets in all the possible pairs of hemispheres. More...
void	CombineGeorgi ()
	Combining the jets in two hemispheres by maximizing (E1-Beta*m1^2/E1 + E2-Beta*m1^2/E2) More...
void	CombineMinEnergyMass ()
	Combining the jets in two hemispheres by minimizing m1^2/E1 + m2^2/E2. More...
void	CombineMinHT ()
void	CombineMinMass ()

Private Attributes
std::vector< float >	Axis1
std::vector< float >	Axis2
std::vector< TLorentzVector >	j1
std::vector< TLorentzVector >	j2
std::vector< TLorentzVector >	jIN
std::vector< TLorentzVector >	jOUT
int	megajet_meth
std::vector< float >	Object_E
std::vector< float >	Object_Px
std::vector< float >	Object_Py
std::vector< float >	Object_Pz
int	status

Detailed Description

Definition at line 21 of file Megajet.h.

Constructor & Destructor Documentation

Megajet() [1/3]

heppy::Megajet::Megajet ( )

inline

Definition at line 35 of file Megajet.h.

{}

Megajet() [2/3]

heppy::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
	)

Definition at line 12 of file Megajet.cc.

References gather_cfg::cout, createfilelist::int, dqmiolumiharvest::j, jIN, Object_E, Object_Px, Object_Py, and Object_Pz.

      : Object_Px(Px_vector),
        Object_Py(Py_vector),
        Object_Pz(Pz_vector),
        Object_E(E_vector),
        megajet_meth(megajet_association_method),
        status(0) {
    if (Object_Px.size() < 2)
      cout << "Error in Megajet: you should provide at least two jets to form Megajets" << endl;
    for (int j = 0; j < (int)jIN.size(); ++j) {
      jIN.push_back(TLorentzVector(Object_Px[j], Object_Py[j], Object_Pz[j], Object_E[j]));
    }
  }

Megajet() [3/3]

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

Definition at line 32 of file Megajet.cc.

References gather_cfg::cout, createfilelist::int, dqmiolumiharvest::j, jIN, Object_E, Object_Px, Object_Py, and Object_Pz.

      : Object_Px(Px_vector), Object_Py(Py_vector), Object_Pz(Pz_vector), Object_E(E_vector), status(0) {
    if (Object_Px.size() < 2)
      cout << "Error in Megajet: you should provide at least two jets to form Megajets" << endl;
    for (int j = 0; j < (int)jIN.size(); ++j) {
      jIN.push_back(TLorentzVector(Object_Px[j], Object_Py[j], Object_Pz[j], Object_E[j]));
    }
  }

~Megajet()

heppy::Megajet::~Megajet ( )

inline

Destructor.

Definition at line 49 of file Megajet.h.

{}

Member Function Documentation

Combine()

void heppy::Megajet::Combine ( )

private

Combine the jets in all the possible pairs of hemispheres.

Definition at line 88 of file Megajet.cc.

References Axis1, Axis2, submitPVResolutionJobs::count, mps_fire::i, createfilelist::int, dqmiolumiharvest::j, j1, j2, and jIN.

Referenced by getAxis1(), and getAxis2().

                        {
    int N_JETS = (int)jIN.size();

    int N_comb = 1;
    for (int i = 0; i < N_JETS; i++) {
      N_comb *= 2;
    }

    // clear some vectors if method Combine() is called again
    if (!j1.empty()) {
      j1.clear();
      j2.clear();
      Axis1.clear();
      Axis2.clear();
    }

    for (int j = 0; j < 5; ++j) {
      Axis1.push_back(0);
      Axis2.push_back(0);
    }

    int j_count;
    for (int i = 1; i < N_comb - 1; i++) {
      TLorentzVector j_temp1, j_temp2;
      int itemp = i;
      j_count = N_comb / 2;
      int count = 0;
      while (j_count > 0) {
        if (itemp / j_count == 1) {
          j_temp1 += jIN[count];
        } else {
          j_temp2 += jIN[count];
        }
        itemp -= j_count * (itemp / j_count);
        j_count /= 2;
        count++;
      }

      j1.push_back(j_temp1);
      j2.push_back(j_temp2);
    }
  }

CombineGeorgi()

void heppy::Megajet::CombineGeorgi ( )

private

Combining the jets in two hemispheres by maximizing (E1-Beta*m1^2/E1 + E2-Beta*m1^2/E2)

Definition at line 186 of file Megajet.cc.

References mps_fire::i, createfilelist::int, j1, j2, jOUT, and status.

Referenced by getAxis1(), and getAxis2().

                              {
    double M_max = -10000;
    // default value (in case none is found)
    TLorentzVector myJ1 = TLorentzVector(0, 0, 0, 0);
    TLorentzVector myJ2 = TLorentzVector(0, 0, 0, 0);
    for (int i = 0; i < (int)j1.size(); i++) {
      int myBeta = 2;
      double M_temp = (j1[i].E() - myBeta * j1[i].M2() / j1[i].E()) + (j2[i].E() - myBeta * j2[i].M2() / j2[i].E());
      if (M_max < -9999 || M_temp > M_max) {
        M_max = M_temp;
        myJ1 = j1[i];
        myJ2 = j2[i];
      }
    }

    //  myJ1.SetPtEtaPhiM(myJ1.Pt(),myJ1.Eta(),myJ1.Phi(),0.0);
    //  myJ2.SetPtEtaPhiM(myJ2.Pt(),myJ2.Eta(),myJ2.Phi(),0.0);

    jOUT.clear();
    if (myJ1.Pt() > myJ2.Pt()) {
      jOUT.push_back(myJ1);
      jOUT.push_back(myJ2);
    } else {
      jOUT.push_back(myJ2);
      jOUT.push_back(myJ1);
    }
    status = 1;
  }

CombineMinEnergyMass()

void heppy::Megajet::CombineMinEnergyMass ( )

private

Combining the jets in two hemispheres by minimizing m1^2/E1 + m2^2/E2.

Definition at line 158 of file Megajet.cc.

References mps_fire::i, createfilelist::int, j1, j2, jOUT, and status.

Referenced by getAxis1(), and getAxis2().

                                     {
    double M_min = -1;
    // default value (in case none is found)
    TLorentzVector myJ1 = TLorentzVector(0, 0, 0, 0);
    TLorentzVector myJ2 = TLorentzVector(0, 0, 0, 0);
    for (int i = 0; i < (int)j1.size(); i++) {
      double M_temp = j1[i].M2() / j1[i].E() + j2[i].M2() / j2[i].E();
      if (M_min < 0 || M_temp < M_min) {
        M_min = M_temp;
        myJ1 = j1[i];
        myJ2 = j2[i];
      }
    }

    //  myJ1.SetPtEtaPhiM(myJ1.Pt(),myJ1.Eta(),myJ1.Phi(),0.0);
    //  myJ2.SetPtEtaPhiM(myJ2.Pt(),myJ2.Eta(),myJ2.Phi(),0.0);

    jOUT.clear();
    if (myJ1.Pt() > myJ2.Pt()) {
      jOUT.push_back(myJ1);
      jOUT.push_back(myJ2);
    } else {
      jOUT.push_back(myJ2);
      jOUT.push_back(myJ1);
    }
    status = 1;
  }

CombineMinHT()

void heppy::Megajet::CombineMinHT ( )

private

Combining the jets in two hemispheres minimizing the difference of HT for the two hemispheres

Definition at line 215 of file Megajet.cc.

References mps_fire::i, createfilelist::int, j1, j2, jOUT, and status.

Referenced by getAxis1(), and getAxis2().

                             {
    double dHT_min = 999999999999999.0;
    // default value (in case none is found)
    TLorentzVector myJ1 = TLorentzVector(0, 0, 0, 0);
    TLorentzVector myJ2 = TLorentzVector(0, 0, 0, 0);
    for (int i = 0; i < (int)j1.size(); i++) {
      double dHT_temp = fabs(j1[i].E() - j2[i].E());
      if (dHT_temp < dHT_min) {
        dHT_min = dHT_temp;
        myJ1 = j1[i];
        myJ2 = j2[i];
      }
    }

    jOUT.clear();
    if (myJ1.Pt() > myJ2.Pt()) {
      jOUT.push_back(myJ1);
      jOUT.push_back(myJ2);
    } else {
      jOUT.push_back(myJ2);
      jOUT.push_back(myJ1);
    }
    status = 1;
  }

CombineMinMass()

void heppy::Megajet::CombineMinMass ( )

private

Combining the jets in two hemispheres minimizing the sum of the invariant masses of the two hemispheres

Definition at line 131 of file Megajet.cc.

References mps_fire::i, createfilelist::int, j1, j2, jOUT, and status.

Referenced by getAxis1(), and getAxis2().

                               {
    double M_min = -1;
    // default value (in case none is found)
    TLorentzVector myJ1 = TLorentzVector(0, 0, 0, 0);
    TLorentzVector myJ2 = TLorentzVector(0, 0, 0, 0);
    for (int i = 0; i < (int)j1.size(); i++) {
      double M_temp = j1[i].M2() + j2[i].M2();
      if (M_min < 0 || M_temp < M_min) {
        M_min = M_temp;
        myJ1 = j1[i];
        myJ2 = j2[i];
      }
    }
    //  myJ1.SetPtEtaPhiM(myJ1.Pt(),myJ1.Eta(),myJ1.Phi(),0.0);
    //  myJ2.SetPtEtaPhiM(myJ2.Pt(),myJ2.Eta(),myJ2.Phi(),0.0);

    jOUT.clear();
    if (myJ1.Pt() > myJ2.Pt()) {
      jOUT.push_back(myJ1);
      jOUT.push_back(myJ2);
    } else {
      jOUT.push_back(myJ2);
      jOUT.push_back(myJ1);
    }
    status = 1;
  }

getAxis1()

vector< float > heppy::Megajet::getAxis1

(

)

Definition at line 41 of file Megajet.cc.

References Axis1, Combine(), CombineGeorgi(), CombineMinEnergyMass(), CombineMinHT(), CombineMinMass(), jIN, jOUT, megajet_meth, and status.

                                  {
    if (status != 1) {
      this->Combine();
      if (megajet_meth == 1)
        this->CombineMinMass();
      else if (megajet_meth == 2)
        this->CombineMinHT();
      else if (megajet_meth == 3)
        this->CombineMinEnergyMass();
      else if (megajet_meth == 4)
        this->CombineGeorgi();
      else
        this->CombineMinMass();
    }
    if (jIN.size() > 1) {
      Axis1[0] = jOUT[0].Px() / jOUT[0].P();
      Axis1[1] = jOUT[0].Py() / jOUT[0].P();
      Axis1[2] = jOUT[0].Pz() / jOUT[0].P();
      Axis1[3] = jOUT[0].P();
      Axis1[4] = jOUT[0].E();
    }
    return Axis1;
  }

getAxis2()

vector< float > heppy::Megajet::getAxis2

(

)

Definition at line 64 of file Megajet.cc.

References Axis2, Combine(), CombineGeorgi(), CombineMinEnergyMass(), CombineMinHT(), CombineMinMass(), jIN, jOUT, megajet_meth, and status.

                                  {
    if (status != 1) {
      this->Combine();
      if (megajet_meth == 1)
        this->CombineMinMass();
      else if (megajet_meth == 2)
        this->CombineMinHT();
      else if (megajet_meth == 3)
        this->CombineMinEnergyMass();
      else if (megajet_meth == 4)
        this->CombineGeorgi();
      else
        this->CombineMinMass();
    }
    if (jIN.size() > 1) {
      Axis2[0] = jOUT[1].Px() / jOUT[1].P();
      Axis2[1] = jOUT[1].Py() / jOUT[1].P();
      Axis2[2] = jOUT[1].Pz() / jOUT[1].P();
      Axis2[3] = jOUT[1].P();
      Axis2[4] = jOUT[1].E();
    }
    return Axis2;
  }

SetMethod()

void heppy::Megajet::SetMethod ( int  megajet_association_method )

inline

Definition at line 60 of file Megajet.h.

References megajet_meth, and status.

                                                   {
      megajet_meth = megajet_association_method;
      status = 0;
    }

Member Data Documentation

Axis1

std::vector<float> heppy::Megajet::Axis1

private

Definition at line 88 of file Megajet.h.

Referenced by Combine(), and getAxis1().

Axis2

std::vector<float> heppy::Megajet::Axis2

private

Definition at line 89 of file Megajet.h.

Referenced by Combine(), and getAxis2().

j1

std::vector<TLorentzVector> heppy::Megajet::j1

private

Definition at line 96 of file Megajet.h.

Referenced by Combine(), CombineGeorgi(), CombineMinEnergyMass(), CombineMinHT(), and CombineMinMass().

j2

std::vector<TLorentzVector> heppy::Megajet::j2

private

Definition at line 97 of file Megajet.h.

Referenced by Combine(), CombineGeorgi(), CombineMinEnergyMass(), CombineMinHT(), and CombineMinMass().

jIN

std::vector<TLorentzVector> heppy::Megajet::jIN

private

Definition at line 94 of file Megajet.h.

Referenced by Combine(), getAxis1(), getAxis2(), and Megajet().

jOUT

std::vector<TLorentzVector> heppy::Megajet::jOUT

private

Definition at line 95 of file Megajet.h.

Referenced by CombineGeorgi(), CombineMinEnergyMass(), CombineMinHT(), CombineMinMass(), getAxis1(), and getAxis2().

megajet_meth

int heppy::Megajet::megajet_meth

private

Definition at line 91 of file Megajet.h.

Referenced by getAxis1(), getAxis2(), and SetMethod().

Object_E

std::vector<float> heppy::Megajet::Object_E

private

Definition at line 86 of file Megajet.h.

Referenced by Megajet().

Object_Px

std::vector<float> heppy::Megajet::Object_Px

private

Definition at line 83 of file Megajet.h.

Referenced by Megajet().

Object_Py

std::vector<float> heppy::Megajet::Object_Py

private

Definition at line 84 of file Megajet.h.

Referenced by Megajet().

Object_Pz

std::vector<float> heppy::Megajet::Object_Pz

private

Definition at line 85 of file Megajet.h.

Referenced by Megajet().

status

int heppy::Megajet::status

private

Definition at line 92 of file Megajet.h.

Referenced by dirstructure.Comparison::__make_image(), dirstructure.Comparison::__repr__(), CombineGeorgi(), CombineMinEnergyMass(), CombineMinHT(), CombineMinMass(), getAxis1(), getAxis2(), and SetMethod().