---

StandardNtupleMaker.h

Go to the documentation of this file.

//  ====================================================================
//  StandardNtupleMaker.h
//  --------------------------------------------------------------------
//
//  This is the header file for the StandardNtupleMaker class
//  which creates the combined Ntuple based on a list of locations
//  in the TES
//
//  It is based upon AtlfastSTLNtupAlg.h by E.Richter-Was
//
//  Author    : H.T.Phillips, P. Clarke, R. Steward, P. Sherwood, E. Richter-Was
//
//  ====================================================================
#ifndef ATLFAST_STANDARDNTUPLEMAKER_H
#define ATLFAST_STANDARDNTUPLEMAKER_H


// places to look in the TES for input
#define DEFAULT_atlfastEventLocation "/Event/Atlfast"
#define DEFAULT_jetLocation "/Event/AtlfastJets"
#define DEFAULT_electronLocation "/Event/AtlfastIsolatedElectrons"
#define DEFAULT_photonLocation "/Event/AtlfastIsolatedPhotons"
#define DEFAULT_isolatedMuonLocation "/Event/AtlfastIsolatedMuons"
#define DEFAULT_nonIsolatedMuonLocation "/Event/AtlfastNonIsolatedMuons"
#define DEFAULT_trackLocation "/Event/AtlfastTracks"
#define DEFAULT_mcTruthLocation "/Event/McEventCollection"
#define DEFAULT_bPhysicsLocation "/Event/AtlfastBPhysics"
#define DEFAULT_triggerLocation "/Event/AtlfastTrigger"
#define DEFAULT_eventHeaderLocation "/Event/AtlfastEventHeader"

#include <string>

// Framework include files
#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
#include "AtlfastCode/JetRecalibrator.h"
#include "AtlfastCode/TesIO.h"
namespace Atlfast {

class StandardNtupleMaker: public Algorithm {
protected:
  // Items for the column wise n-tuple
  // Miscs
  NTuple::Item<long>              m_isub;
  NTuple::Item<long>              m_njetb, m_njetc, m_njettau;
  NTuple::Item<float>             m_pxmiss, m_pymiss, m_pxnue, m_pynue;
  // Jets
  NTuple::Item<long>              m_njet;
  NTuple::Array<long>             m_codejet;
  NTuple::Array<float>            m_pxjet, m_pyjet, m_pzjet, m_eejet;
  NTuple::Array<float>            m_ptcalo, m_ptbjet, m_ptujet;
  // Electrons
  NTuple::Item<long>              m_nele;
  NTuple::Array<long>             m_codeele;
  NTuple::Array<float>            m_pxele, m_pyele, m_pzele, m_eeele;
  // Muons
  NTuple::Item<long>              m_nmuo;
  NTuple::Array<long>             m_codemuo;
  NTuple::Array<float>            m_pxmuo, m_pymuo, m_pzmuo, m_eemuo;
  NTuple::Item<long>              m_nmux;
  NTuple::Array<long>             m_codemux;
  NTuple::Array<float>            m_pxmux, m_pymux, m_pzmux, m_eemux;
  // Photons
  NTuple::Item<long>              m_npho;
  NTuple::Array<long>             m_codepho;
  NTuple::Array<float>            m_pxpho, m_pypho, m_pzpho, m_eepho;
  // Tracks
  NTuple::Item<long>              m_ntra;
  NTuple::Array<long>             m_kpTruth,   m_kfTruth;
  NTuple::Array<long>             m_kpm1tra, m_kfm1tra;
  NTuple::Array<long>             m_kpm2tra, m_kfm2tra;
  NTuple::Array<long>             m_kpm3tra, m_kfm3tra;
  NTuple::Array<long>             m_kpm4tra, m_kfm4tra;
  NTuple::Array<long>             m_kpm5tra, m_kfm5tra;
  NTuple::Array<long>             m_kpm6tra, m_kfm6tra;
  NTuple::Array<float>            m_d0Track, m_z0Track, m_phiTrack, m_cotTrack, m_ptInvTrack;
  NTuple::Array<float>            m_d0Truth, m_z0Truth, m_phiTruth, m_cotTruth, m_ptInvTruth;
  NTuple::Array<float>            m_corr11tra, m_corr12tra, m_corr13tra, m_corr14tra, m_corr15tra;
  NTuple::Array<float>            m_corr22tra, m_corr23tra, m_corr24tra, m_corr25tra;
  NTuple::Array<float>            m_corr33tra, m_corr34tra, m_corr35tra;
  NTuple::Array<float>            m_corr44tra, m_corr45tra, m_corr55tra;
  // Particle
  NTuple::Item<long>              m_npart;
  NTuple::Array<long>             m_kfpart;
  NTuple::Array<float>            m_pxpart, m_pypart, m_pzpart, m_eepart ;
  // Bphys
  NTuple::Item<long>              m_nbphys;
  NTuple::Array<long>             m_kbphys, m_kfbphys;
  NTuple::Array<long>             m_pxphys, m_pybphys, m_pzbphys, m_eebphys ;
  NTuple::Array<long>             m_vxphys, m_vybphys, m_vzbphys ;


  // Pointer to an ntuple
  NTuple::Tuple* p_nt1;

public:
  StandardNtupleMaker(const std::string& name, ISvcLocator* pSvcLocator); 
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();

private:
  // Counters for event processing
  long   m_nEvents;                    // Number events processed
  TesIO* m_tesIO;

  // member variables pointing to locations in the TES
  std::string m_atlfastEventLocation;
  std::string m_jetLocation;
  std::string m_electronLocation;
  std::string m_isolatedMuonLocation;
  std::string m_nonIsolatedMuonLocation;
  std::string m_photonLocation;
  std::string m_trackLocation;
  std::string m_mcTruthLocation;
  std::string m_bPhysicsLocation;
  std::string m_triggerLocation;  
  std::string m_eventHeaderLocation;  

  //  IJetCalibrator& m_uJetCal;
  //  IJetCalibrator& m_bJetCal;
  JetRecalibrator*  m_jetCal;

  typedef std::vector<HepMC::GenParticle*>            t_MC_particleCollection; 
  typedef std::vector<HepMC::GenParticle*>::iterator  t_MC_particleIterator;
};

} // end of namespace bracket  

#endif 

---