#include <PionZeroList.h>

Classes
class	f_less

Public Member Functions
	PionZeroList (Vp4 &input_vp)

void	set_cut (int index, double cut)

void	sort ()

void	reduce ()

void	print ()

void	refresh ()

int	get_num_pi0 ()

Vp4	get_pi0_list ()

Detailed Description

Definition at line 32 of file PionZeroList.h.

Constructor & Destructor Documentation

◆ PionZeroList()

PionZeroList::PionZeroList ( Vp4 & input_vp )

Definition at line 46 of file PionZeroList.cxx.

                                       :  
  m_gam_vp(input_vp), 
  m_pi0_vp(0),
  Vpi0_list(0),
  m_low_cut(0.11),
  m_high_cut(0.15){
  int m_Ngam(input_vp.size());
  if(m_Ngam<2) cout << "number of photons less than 2" << endl;
  for(int i=0; i<m_Ngam-1; i++){
    for(int j=i+1; j<m_Ngam; j++){
      double m_temp_mass((m_gam_vp[i]+m_gam_vp[j]).m());
      if(m_temp_mass<m_low_cut || m_temp_mass>m_high_cut) continue;
      PionZero m_temp_pi0(m_gam_vp[i], i, m_gam_vp[j], j);
      Vpi0_list.push_back(m_temp_pi0);
    }
  }
  m_num_pi0 = Vpi0_list.size();
}

Member Function Documentation

◆ get_num_pi0()

int PionZeroList::get_num_pi0 ( )

inline

Definition at line 40 of file PionZeroList.h.

40{ return m_num_pi0; }

◆ get_pi0_list()

Vp4 PionZeroList::get_pi0_list ( )

inline

Definition at line 41 of file PionZeroList.h.

41{ return m_pi0_vp; }

◆ print()

void PionZeroList::print ( )

Definition at line 108 of file PionZeroList.cxx.

                        {
  std::cout << "number of pion0 in the list " << m_num_pi0 << std::endl;
  for(int i=0; i<Vpi0_list.size(); i++){
    std::cout << "i= " << i << endl;
    std::cout << "pi0 mass " << Vpi0_list[i].get_mass() << endl; 
    std::cout << "pi0 goodness " << Vpi0_list[i].get_goodness() << endl;
    std::cout << "pi0 four momentum " << Vpi0_list[i].get_pi0_vp() << endl;
    std::cout << "fir gam index " << Vpi0_list[i].get_index(0) << endl;
    std::cout << "sec gam index " << Vpi0_list[i].get_index(1) << endl;
    std::cout << "fir gam mom " << Vpi0_list[i].get_gam_vp(0) << endl;
    std::cout << "sec gam mom " << Vpi0_list[i].get_gam_vp(1) << endl;
  }
}

◆ reduce()

void PionZeroList::reduce ( )

Definition at line 84 of file PionZeroList.cxx.

                         {
  std::vector<int> gam_index_v(0);
  std::vector<PionZero> temp_pi0_list(0);
  for(int i=0; i<Vpi0_list.size(); i++){
    int temp_one(Vpi0_list[i].get_index(0));
    int temp_two(Vpi0_list[i].get_index(1));
    if(find(gam_index_v.begin(), gam_index_v.end(), temp_one) != gam_index_v.end() ) continue;
    if(find(gam_index_v.begin(), gam_index_v.end(), temp_two) != gam_index_v.end() ) continue;
    gam_index_v.push_back(temp_one);
    gam_index_v.push_back(temp_two);
    temp_pi0_list.push_back(Vpi0_list[i]);
  }
  Vpi0_list = temp_pi0_list;
  refresh();
}

◆ refresh()

void PionZeroList::refresh ( )

Definition at line 101 of file PionZeroList.cxx.

                          {
  m_num_pi0 = Vpi0_list.size();
  m_pi0_vp.clear();
  for(int i=0; i<m_num_pi0; i++) m_pi0_vp.push_back(Vpi0_list[i].get_pi0_vp());
}

Referenced by reduce(), set_cut(), and sort().

◆ set_cut()

void PionZeroList::set_cut	(	int	index,
		double	cut )

Definition at line 66 of file PionZeroList.cxx.

                                               {
  if(index < 1) m_low_cut = cut;
  else m_high_cut = cut;
  std::vector<PionZero>::iterator m_ind;
  for(m_ind=Vpi0_list.begin(); m_ind!=Vpi0_list.end(); m_ind++){
    double m_temp_mass((*m_ind).get_mass());
    if(m_temp_mass<m_low_cut||m_temp_mass>m_high_cut) Vpi0_list.erase(m_ind); 
  }
  refresh();
}

◆ sort()

void PionZeroList::sort ( )

Definition at line 78 of file PionZeroList.cxx.

                       { // sort the pion0 list respect to goodness
  std::sort(Vpi0_list.begin(), Vpi0_list.end(), f_less());
  refresh();
}

The documentation for this class was generated from the following files:

Classes

Public Member Functions