NtupleReader/NtupleReader/MyParticlePair.ixx

//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// 10.12.2006, AUTHOR: STEFFEN KAISER
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

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

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

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

inline void MyParticlePair::init(MyParticle* particle1,
                                 MyParticle* particle2) {

    m_particle1 = particle1;
    m_particle2 = particle2;

    return;

}

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

//:::::::::::::::::::::::::
//:: METHOD get_particle1 ::
//:::::::::::::::::::::::::

inline MyParticle* MyParticlePair::get_particle1(void) const {

    return m_particle1;

}

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

//::::::::::::::::::::::::::
//:: METHOD get_particle2 ::
//::::::::::::::::::::::::::

inline MyParticle* MyParticlePair::get_particle2(void) const {

    return m_particle2;

}

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

//::::::::::::::::
//:: METHOD tlv ::
//::::::::::::::::

inline TLorentzVector MyParticlePair::tlv(void) const {

    return (TLorentzVector)(m_particle1->tlv() + m_particle2->tlv());

}

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

//::::::::::::::::::::
//:: METHOD EtaProd ::
//::::::::::::::::::::

inline double MyParticlePair::EtaProd(void) const {

        return (double)(m_particle1->tlv().Eta() * m_particle2->tlv().Eta());

}

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

//:::::::::::::::::::::
//:: METHOD DeltaPhi ::
//:::::::::::::::::::::

inline double MyParticlePair::DeltaPhi(void) const {

    return (double) (m_particle1->tlv()).DeltaPhi(m_particle2->tlv());

}

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

//::::::::::::::::::::::::
//:: METHOD AbsDeltaPhi ::
//::::::::::::::::::::::::

inline double MyParticlePair::AbsDeltaPhi(void) const {

    return TMath::Abs( DeltaPhi() ) ;

}

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

//:::::::::::::::::::
//:: METHOD DeltaR ::
//:::::::::::::::::::

inline double MyParticlePair::DeltaR(void) const {

    return (double) (m_particle1->tlv()).DeltaR(m_particle2->tlv());

}

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

//::::::::::::::
//:: METHOD M ::
//::::::::::::::

inline double MyParticlePair::M(void) const {

    return (double) (m_particle1->tlv()+m_particle2->tlv()).M();

}

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

//:::::::::::::::::::::
//:: METHOD DeltaEta ::
//:::::::::::::::::::::

inline double MyParticlePair::DeltaEta(void) const {

    double eta1 =  m_particle1->tlv().Eta();
    double eta2 =  m_particle2->tlv().Eta();

    return (eta1 - eta2) ;

}

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

//::::::::::::::::::::::::
//:: METHOD AbsDeltaEta ::
//::::::::::::::::::::::::

inline double MyParticlePair::AbsDeltaEta(void) const {

    return TMath::Abs( DeltaEta() ) ;

}

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

//::::::::::::::::::::
//:: METHOD DeltaPt ::
//::::::::::::::::::::

inline double MyParticlePair::DeltaPt(void) const {

    double pt1 =  m_particle1->tlv().Pt();
    double pt2 =  m_particle2->tlv().Pt();

    return (pt1 - pt2) ;

}

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

//:::::::::::::::::::::::
//:: METHOD AbsDeltaPt ::
//:::::::::::::::::::::::

inline double MyParticlePair::AbsDeltaPt(void) const {

    return TMath::Abs( DeltaPt() ) ;

}


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

//:::::::::::::::::::::::::::::::
//:: METHOD get_lepton_flavors ::
//:::::::::::::::::::::::::::::::

TString
MyParticlePair::get_lepton_flavors(void) {

    int pdgId1 = m_particle1->pdgId();
    int pdgId2 = m_particle2->pdgId();

    TString type1 = pdgIdToString(pdgId1);
    TString type2 = pdgIdToString(pdgId2);

    TString lp_type = type1+type2;

    if (lp_type == "mue"){
        lp_type = "emu";
    }

    if(lp_type.Contains("unknown")){
        std::cout << "WARNING: get_lepton_flavors applied on particle pair "
             << "without two leptons!" << std::endl;
    }

    return lp_type;

  }

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

//::::::::::::::::::::::::::
//:: METHOD pdgIdToString ::
//::::::::::::::::::::::::::

TString
MyParticlePair::pdgIdToString(int pdgId){

    if (TMath::Abs(pdgId)==11){
        return TString("e");
    }
    else if(TMath::Abs(pdgId)==13){
        return TString("mu");
    }
    else{
        return TString("unknown");
            //std::cerr << "ERROR: Wrong lepton type in lepon pair pdgId = " << pdgId 
            //<< "." << std::endl;
            //exit(0);
    }


  }

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

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

inline void MyParticlePair::Print(std::string option) {

    std::cout << "ParticlePair: " << std::endl;
    m_particle1->Print(option);
    m_particle2->Print(option);

}

---