NtupleReader/src/MyPileUpExtractedEvent.cxx

//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// 02.08.2007, AUTHOR: OLIVER KORTNER
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

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

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

#include "NtupleReader/MyPileUpExtractedEvent.h"

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

using namespace std;

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

  //:::::::::::::::::
  //:: CONSTRUCTOR ::
  //:::::::::::::::::

MyPileUpExtractedEvent::MyPileUpExtractedEvent(const MyEvent & event) {

    //---------------
    //-- VARIABLES --
    //---------------

    int vertex(0);  // index of the vertex with highest total
                    // transverse track energy
    double Et_vert; // transverse vertex energy

    //--------------------------------------------------------------------
    //-- LOOK FOR THE VERTEX WITH HIGHEST TOTAL TRANSVERSE TRACK ENERGY --
    //--------------------------------------------------------------------
    
    // no primary vertex has been reconstructed //
    if (event.primary_vertices()->size()==0) {
        vertex = -1;
        m_vertex = NULL;
    }

    // if there is a primary vertex, find the right one //
    if (vertex!=-1) {
        
        Et_vert = (*(event.primary_vertices()))[0]->Et();
        
        for (unsigned int k=1; k<event.primary_vertices()->size(); k++) {
            if (Et_vert<(*(event.primary_vertices()))[k]->Et()) {
                Et_vert=(*(event.primary_vertices()))[k]->Et();
                vertex = k;
            }
        }
        
        m_vertex = (*(event.primary_vertices()))[vertex]; 
    }

    int truth_vertex(0);
    Et_vert=0;

    if (event.truth_vertices()->size()==0) {
        truth_vertex = -1;
        m_truth_vertex = NULL;
    }

    // if there is a truth vertex, find the right one //
    if (truth_vertex!=-1) {
        
        Et_vert = (*(event.truth_vertices()))[0]->Et();
        
        for (unsigned int k=1; k<event.truth_vertices()->size(); k++) {
            if (Et_vert<(*(event.truth_vertices()))[k]->Et()) {
                Et_vert=(*(event.truth_vertices()))[k]->Et();
                truth_vertex = k;
            }
        }
        
        m_truth_vertex = (*(event.truth_vertices()))[truth_vertex]; 
    }
    

    //---------------------------------------------------
    //-- COLLECT ELECTRONS AND PHOTON FROM THAT VERTEX --
    //---------------------------------------------------

    // electrons //
    m_rec_electron.clear();
    for (unsigned int k=0; k<(event.reco_electrons()).size(); k++) {
        
        if ((event.reco_electrons())[k]->vertexIndex()==vertex) {
        
            m_rec_electron.push_back((event.reco_electrons())[k]);
        }
    }
   
    // muons //
    m_rec_muon.clear();
    std::vector<MyParticle*> muons = event.reco_muons("all");
    
    for (unsigned int k=0; k<muons.size(); k++) {
        
        if ( muons[k]->vertexIndex() == vertex ){
            m_rec_muon.push_back(muons[k]);
        }
    }
    
    
    //-----------------------------------
    //-- COPY THE OTHER EVENT CONTENTS --
    //-----------------------------------

    m_rec_jet        = event.reco_jets();
    m_truth_jet      = event.truth_jets();
    m_track_jet      = event.track_jets();
    m_track_jet      = event.truth_track_jets();
    m_rec_taujet     = event.reco_taujets();
    m_met_manager    = event.met_manager();
    m_trigger        = event.trigger();
    m_truth_particle = event.truth_particles();
        
    m_dataset_info = &MyDatasetInfo::Instance();
    
}

---