EtMissAnalysis.cc

Go to the documentation of this file.

#include <iostream>
#include <stdio.h>
#include "HepMC/GenEvent.h"
#include "HepMC/IO_GenEvent.h"
#include "HepMC/GenParticle.h"
#include "HepMC/GenVertex.h"
#include "HepMC/IO_AsciiParticles.h"
#include "HepMC/SimpleVector.h"
#include "HepPDT/ParticleData.hh"
#include "CLHEP/Vector/LorentzVector.h"
using namespace std;

// ROOT headers
#include "TH1.h"
#include "TH2.h"
#include "TFile.h"
#include "TMath.h"
#include "TLorentzVector.h"

#include "fastjet/PseudoJet.hh"
#include "fastjet/ClusterSequence.hh"
#include "fastjet/JetDefinition.hh"
#include "fastjet/SISConePlugin.hh"

//EtMiss analysis header
#include "../include/EtMissAnalysis.h"

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

/// empty default constructor
EtMissAnalysis::EtMissAnalysis()
{
}

/// empty default destructor
EtMissAnalysis::~EtMissAnalysis()
{
}

int EtMissAnalysis::Init(double tr_max_eta, double tr_min_pt)
{
  // specify default eta cut
  m_max_eta=tr_max_eta;  
  
  // specify default pt cut
  m_min_pt=tr_min_pt;  
   
  // default Output file name
  m_outputFileName="EtMissAnalysis.root";
  m_outputRootDir="EtMiss";
 
  m_METTruthX =  baseAnalysis::initHist( string("METTruthX"), string("METTruthX"), string("MissingEtTruthX"), 400, -1000, 1000);
  m_METTruthY =  baseAnalysis::initHist( string("METTruthY"), string("METTruthY"), string("MissingEtTruthY"), 400, -1000, 1000);
  m_METTruth =   baseAnalysis::initHist( string("METTruth"), string("METTruth"), string("MissingEtTruth"), 200, 0, 1000);
  m_METSumTruth =  baseAnalysis::initHist( string("METSumTruth"), string("METSumTruth"), string("MissingSumEtTruth"), 200, 0, 1000);
  
  return true; 
}


int EtMissAnalysis::Process(HepMC::GenEvent* hepmcevt){ 

  //calculate the missing energy with function from baseAnalysis
  FindMissingEt(hepmcevt);

  m_METTruthX->Fill(exMissTruth);
  m_METTruthY->Fill(eyMissTruth);
  m_METTruth->Fill(etMissTruth);
  m_METSumTruth->Fill(etsumMissTruth);

  return true;

} //end of Process

---