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

heppy::Megajet::Megajet ( )

inline

Definition at line 37 of file Megajet.h.

{};

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, 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]));
  }
}

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 25 of file Megajet.cc.

References gather_cfg::cout, createfilelist::int, 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]));
  }
}

heppy::Megajet::~Megajet ( )

inline

Destructor.

Definition at line 44 of file Megajet.h.

References getAxis1(), and getAxis2().

{};

Member Function Documentation

void heppy::Megajet::Combine ( )

private

Combine the jets in all the possible pairs of hemispheres.

Definition at line 73 of file Megajet.cc.

References Axis1, Axis2, KineDebug3::count(), i, createfilelist::int, j, j1, j2, and jIN.

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

                      {
  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);
  }
}

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 171 of file Megajet.cc.

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

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

                           {
  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;
}

void heppy::Megajet::CombineMinEnergyMass ( )

private

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

Definition at line 143 of file Megajet.cc.

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

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

                                   {
  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;
}

void heppy::Megajet::CombineMinHT ( )

private

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

Definition at line 200 of file Megajet.cc.

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

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

                           {
  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;
}

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 116 of file Megajet.cc.

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

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

                             {
  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;
}

vector< float > heppy::Megajet::getAxis1

(

)

Definition at line 36 of file Megajet.cc.

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

Referenced by ~Megajet().

                               {
  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;
}

vector< float > heppy::Megajet::getAxis2

(

)

Definition at line 54 of file Megajet.cc.

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

Referenced by ~Megajet().

                               {
  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;
}

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

inline

Definition at line 55 of file Megajet.h.

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

                                                 {
    megajet_meth = megajet_association_method;
    status = 0;
  }

Member Data Documentation

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

private

Definition at line 85 of file Megajet.h.

Referenced by Combine(), and getAxis1().

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

private

Definition at line 86 of file Megajet.h.

Referenced by Combine(), and getAxis2().

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

private

Definition at line 94 of file Megajet.h.

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

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

private

Definition at line 95 of file Megajet.h.

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

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

private

Definition at line 92 of file Megajet.h.

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

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

private

Definition at line 93 of file Megajet.h.

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

int heppy::Megajet::megajet_meth

private

Definition at line 89 of file Megajet.h.

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

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

private

Definition at line 83 of file Megajet.h.

Referenced by Megajet().

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

private

Definition at line 80 of file Megajet.h.

Referenced by Megajet().

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

private

Definition at line 81 of file Megajet.h.

Referenced by Megajet().

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

private

Definition at line 82 of file Megajet.h.

Referenced by Megajet().

int heppy::Megajet::status

private

Definition at line 90 of file Megajet.h.

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