44 : m_gam_vp( input_vp )
45 , m_pi0_vp( 0 )
46 , Vpi0_list( 0 )
47 , m_low_cut( 0.11 )
48 , m_high_cut( 0.15 ) {
49 int m_Ngam( input_vp.size() );
50 if ( m_Ngam < 2 ) cout << "number of photons less than 2" << endl;
51 for ( int i = 0; i < m_Ngam - 1; i++ )
52 {
53 for ( int j = i + 1; j < m_Ngam; j++ )
54 {
55 double m_temp_mass( ( m_gam_vp[i] + m_gam_vp[j] ).m() );
56 if ( m_temp_mass < m_low_cut || m_temp_mass > m_high_cut ) continue;
57 PionZero m_temp_pi0( m_gam_vp[i], i, m_gam_vp[j], j );
58 Vpi0_list.push_back( m_temp_pi0 );
59 }
60 }
61 m_num_pi0 = Vpi0_list.size();
62}