NtupleReader/NtupleReader/MyPhoton.ixx

//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// 23.10.2006, AUTHOR: OLIVER KORTNER
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

//::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
//:: IMPLEMENTATION OF INLINE METHODS DEFINED IN THE CLASS MyPhoton ::
//::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::

//*****************************************************************************

  //:::::::::::::::::
  //:: METHOD init ::
  //:::::::::::::::::

inline void MyPhoton::init(const TLorentzVector & p,
                             const double & echarge,
                             const bool & egamma_flag,
                             const double & isem,
                             const double & neuralnet,
                             const int & index,
                             const double et_cone[4]) {

    m_tlv = p;
    m_charge = echarge;
    m_egamma_flag = egamma_flag;
    m_isEM = (int)isem;
    m_NeuralNet = neuralnet;
    m_index = index;
    m_type = "photon";
    m_pdgId = (int)m_charge*22;
    for (unsigned int k=0; k<4; k++) {
        if (et_cone!=0) {
            m_Et_in_cone[k] = et_cone[k];
        } else {
            m_Et_in_cone[k] = 0.0;
        }
    }
    m_matchingflag = 0;

    return;

  }

//*****************************************************************************

  //:::::::::::::::::::::::::
  //:: METHOD e_gamma_flag ::
  //:::::::::::::::::::::::::

inline bool MyPhoton::e_gamma_flag(void) const {

    return m_egamma_flag;

  }

//*****************************************************************************

  //::::::::::::::::::
  //:: METHOD is_em ::
  //::::::::::::::::::

inline int MyPhoton::is_em(void) const {

    return m_isEM;

  }
//*****************************************************************************

  //::::::::::::::::::::::
  //:: METHOD NeuralNet ::
  //::::::::::::::::::::::

inline double MyPhoton::NeuralNet(void) const {

    return m_NeuralNet;

  }

//*****************************************************************************

  //:::::::::::::::::::::::
  //:: METHOD Et_in_cone ::
  //:::::::::::::::::::::::

inline double MyPhoton::Et_in_cone(const int & delta_R_index) const {

    // CHECK INDEX //

    if (delta_R_index<0 || delta_R_index>4) {
        std::cerr << "\n"
                  << "Class MyPhoton, method Et_in_cone: ERROR!\n"
                  << "Wrong delta_R_index " << delta_R_index
                  << ", must be >=0 and <=3.\n";
        exit(1);
    }

    // RETURN THE REQUESTED ENERGY //

    return m_Et_in_cone[delta_R_index];

  }

//*****************************************************************************

  //:::::::::::::::::::::::::::::
  //:: METHOD set_e_gamma_flag ::
  //:::::::::::::::::::::::::::::

inline void MyPhoton::set_e_gamma_flag(const bool & flag) {

    m_egamma_flag = flag;
    return;

  }

//*****************************************************************************

  //:::::::::::::::::::::::::::
  //:: METHOD set_Et_in_cone ::
  //:::::::::::::::::::::::::::

inline void MyPhoton::set_Et_in_cone(const double et_in_cone[4]) {

    for (int k=0; k<4; k++) {
        m_Et_in_cone[k] = et_in_cone[k];
    }

    return;

  }

//*****************************************************************************

  //::::::::::::::::::::::::::
  //:: METHOD PrintParticle ::
  //::::::::::::::::::::::::::

inline void MyPhoton::PrintParticle(std::string option) {
    
    if(option=="v"){
        std::cout << "   is_em:                       " 
                  <<  m_isEM << std::endl;
    }
    
    return;
  }

---