NtupleReader/src/MyVBFCandidate.cxx

//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// 19.04.2007, AUTHOR: STEFFEN KAISER
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

//:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
//:: IMPLEMENTATION OF METHODS DEFINED IN THE CLASS MyVBFCandidate ::
//:::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::

//::::::::::::::::::
//:: HEADER FILES ::
//::::::::::::::::::

#include "NtupleReader/MyVBFCandidate.h"


//::::::::::::::::::::::::
//:: NAMESPACE SETTINGS ::
//::::::::::::::::::::::::

using namespace std;


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

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

void MyVBFCandidate::init( MyParticlePair tagjet_pair,
                           MyParticlePair lepton_pair,
                           const MyMissingEt & missing_et) {

    has_tagjets=true;
    has_leptons=true;
    cut_map=0;
    m_tagjet_pair = tagjet_pair;
    m_lepton_pair = lepton_pair;
    m_missing_et  = missing_et;

  }

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

//:::::::::::::::::::::
//:: METHOD destruct ::
//:::::::::::::::::::::

void MyVBFCandidate::destruct(void) {

    return;

}

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

//::::::::::::::::::::::::::::::::
//:: METHOD set_has_tagjet_pair ::
//::::::::::::::::::::::::::::::::

void
MyVBFCandidate::set_has_tagjet_pair(bool tagjet_pair_exists) {
    
    has_tagjets = tagjet_pair_exists;
    return ;

}

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

//::::::::::::::::::::::::::::::::
//:: METHOD set_has_lepton_pair ::
//::::::::::::::::::::::::::::::::

void
MyVBFCandidate::set_has_lepton_pair(bool lepton_pair_exists) {

    has_leptons = lepton_pair_exists;
    return ;

}

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

//::::::::::::::::::::::::::::
//:: METHOD has_tagjet_pair ::
//::::::::::::::::::::::::::::

bool
MyVBFCandidate::has_tagjet_pair(void) const {

    return has_tagjets;

}

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

//::::::::::::::::::::::::::::
//:: METHOD has_tagjet_pair ::
//::::::::::::::::::::::::::::

bool
MyVBFCandidate::has_lepton_pair(void) const {

    return has_leptons;

}

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

//::::::::::::::::::::::::
//:: METHOD set_cut_map ::
//::::::::::::::::::::::::

void 
MyVBFCandidate::set_cut_map(long bitmap){

    cut_map = bitmap;
    return ;

}

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

//::::::::::::::::::::::::
//:: METHOD get_cut_map ::
//::::::::::::::::::::::::

long
MyVBFCandidate::get_cut_map(void) const {

    return cut_map;

}

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

//::::::::::::::::::::::::::::
//:: METHOD get_tagjet_pair ::
//::::::::::::::::::::::::::::

MyParticlePair 
MyVBFCandidate::get_tagjet_pair(void) const {

    return m_tagjet_pair;

  }

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

//::::::::::::::::::::::::::::
//:: METHOD get_lepton_pair ::
//::::::::::::::::::::::::::::

MyParticlePair
MyVBFCandidate::get_lepton_pair(void) const {
    
    return m_lepton_pair;

}

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

//::::::::::::::::::::::::::
//:: METHOD get_missinget ::
//::::::::::::::::::::::::::

MyMissingEt
MyVBFCandidate::get_missinget(void) const {
    
    return m_missing_et;

}

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

//::::::::::::::::::::::::::::::::::::::::::::
//:: METHOD get_get_collinear_approximation ::
//::::::::::::::::::::::::::::::::::::::::::::

MyCollinearApproximation
MyVBFCandidate::get_collinear_approximation(void)
    const {MyCollinearApproximation colap;
    if ( has_leptons ) {
        colap = MyCollinearApproximation(m_lepton_pair.get_particle1(),
                                       m_lepton_pair.get_particle2(),
                                       m_missing_et.tlv());
    }

    return colap;

    }

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

//:::::::::::::::::::::::
//:: METHOD get_pt_sum ::
//:::::::::::::::::::::::

TLorentzVector
MyVBFCandidate::get_pt_sum(void) const {

    TLorentzVector tmp_tlv;
    tmp_tlv = m_lepton_pair.get_particle1()->tlv();
    tmp_tlv += m_lepton_pair.get_particle2()->tlv();
    tmp_tlv += m_tagjet_pair.get_particle1()->tlv();
    tmp_tlv += m_tagjet_pair.get_particle2()->tlv();
    tmp_tlv += m_missing_et.tlv();

    return tmp_tlv;

}

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

//:::::::::::::::::::::::::::::
//:: METHOD get_lj_DeltaEta1 ::
//:::::::::::::::::::::::::::::

double
MyVBFCandidate::get_lj_DeltaEta1(void) const {
    if ( !has_tagjets || !has_leptons ){
        cerr<<"Error in MyVBFCandidate::get_lj_AbsDeltaEta1: "<<endl
                <<"You can not call this function "
                <<"if the candidate has no tagjet or lepton pair!"<<endl;
        return -9999.;
    }

    double jet1_eta = m_tagjet_pair.get_particle1()->tlv().Eta();
    double jet2_eta = m_tagjet_pair.get_particle2()->tlv().Eta();
    double lepton1_eta = m_lepton_pair.get_particle1()->tlv().Eta();
    double lepton2_eta = m_lepton_pair.get_particle2()->tlv().Eta();

    return ( TMath::Min(lepton1_eta, lepton2_eta) - TMath::Min(jet1_eta,jet2_eta) );

}

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

//:::::::::::::::::::::::::::::
//:: METHOD get_lj_DeltaEta2 ::
//:::::::::::::::::::::::::::::

double
MyVBFCandidate::get_lj_DeltaEta2(void) const {
    if ( !has_tagjets || !has_leptons ){
        cerr<<"Error in MyVBFCandidate::get_lj_AbsDeltaEta2: "<<endl
            <<"You can not call this function "
            <<"if the candidate has no tagjet or lepton pair!"<<endl;
        return -9999.;
    }

    double jet1_eta = m_tagjet_pair.get_particle1()->tlv().Eta();
    double jet2_eta = m_tagjet_pair.get_particle2()->tlv().Eta();
    double lepton1_eta = m_lepton_pair.get_particle1()->tlv().Eta();
    double lepton2_eta = m_lepton_pair.get_particle2()->tlv().Eta();

    return ( TMath::Max(jet1_eta,jet2_eta) - TMath::Max(lepton1_eta, lepton2_eta) );

}

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

//:::::::::::::::::::::::::::::::
//:: METHOD get_lj_DeltaEtaMin ::
//:::::::::::::::::::::::::::::::

double
MyVBFCandidate::get_lj_DeltaEtaMin(void) const {

    return TMath::Min(get_lj_DeltaEta1(),get_lj_DeltaEta2());

}

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

//:::::::::::::::::::::::::::::::
//:: METHOD get_lj_DeltaEtaMax ::
//:::::::::::::::::::::::::::::::

double
MyVBFCandidate::get_lj_DeltaEtaMax(void) const {

    return TMath::Max(get_lj_DeltaEta1(),get_lj_DeltaEta2());

}

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

//::::::::::::::::::
//:: METHOD Print ::
//::::::::::::::::::

void 
MyVBFCandidate::Print(string option) {

    cout << "VBFCandidate: " << endl;
    m_tagjet_pair.Print(option);
    m_lepton_pair.Print(option);
    if(option=="v"){
        m_missing_et.Print(option); 
    }

  }


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

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

// TString
// MyVBFCandidate::get_lepton_flavors(void) {
    
//     int pdgId1 = m_lepton_pair.get_particle1()->pdgId();
//     int pdgId2 = m_lepton_pair.get_particle2()->pdgId();
    
//     TString type1 = pdgIdToString(pdgId1);
//     TString type2 = pdgIdToString(pdgId2);
    
//     TString lp_type = type1+type2;
    
//     if (lp_type == "mue"){
//         lp_type = "emu";
//     }
    
//     return lp_type;
    
//   }

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

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

// TString
// MyVBFCandidate::pdgIdToString(int pdgId){

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

    
//   }

---