ETMiss.h

Go to the documentation of this file.

#ifndef _ETMISS
#define _ETMISS

#include <stdio.h>
#include <math.h>

#include <TObject.h>
#include <TROOT.h>

#include "classes/DelphesClasses.h"


class ETMiss {
 public:

    ETMiss(MissingET* met) {
        double missingET_ET = met->MET;
        double missingET_Phi = met->Phi;
        double missingET_Ex = missingET_ET*cos(missingET_Phi);
        double missingET_Ey = missingET_ET*sin(missingET_Phi);
        
        // Smearing of ET due to pile up set at 20GeV
        double deltaPT = 20.0;
        
        double x = rand() /(RAND_MAX+1.);
        double y = rand() /(RAND_MAX+1.);
        double uni_gauss = sqrt(-2.*log(x))*cos(2.*3.1415*y);
        double uni_normal = rand() / (RAND_MAX+1.);

        double smear_x = deltaPT*uni_gauss*cos(2.*3.1415*uni_normal);
        double smear_y = deltaPT*uni_gauss*sin(2.*3.1415*uni_normal);

        double new_ET = sqrt(pow(missingET_Ex+smear_x, 2)+pow(missingET_Ey+smear_y, 2));
        content.SetPxPyPzE(missingET_Ex+smear_x, missingET_Ey+smear_y, 0, new_ET);
        PT = content.Pt();
        Eta = content.Eta();
        Phi = content.Phi();    
    }
    

    void addMuons(std::vector<Muon*> muons) {
        TLorentzVector sum = TLorentzVector(0,0,0,0);
        for(int i = 0; i < muons.size(); i++)
            sum += muons[i]->P4();
      
        double new_ET = sqrt(pow(content.Px()-sum.Px(), 2)+pow(content.Py()-sum.Py(), 2));
        content.SetPxPyPzE(content.Px()-sum.Px(), content.Py()-sum.Py(), 0, new_ET);
        PT = content.Pt();
        Eta = content.Eta();
        Phi = content.Phi();
    }
    
    Float_t PT; 
    Float_t Eta; 
    Float_t Phi; 
    TLorentzVector P4() { return content; } 

 private:
    TLorentzVector content;
};

#endif