StandardNtupleMaker.h

Go to the documentation of this file.

#ifndef ATLFAST_STANDARDNTUPLEMAKER_H
#define ATLFAST_STANDARDNTUPLEMAKER_H

#include <string>

#include "GaudiKernel/Algorithm.h"
#include "GaudiKernel/NTuple.h"
#include "AtlfastUtils/TesIO.h"
namespace Atlfast {
  using std::string;
  class StandardNtupleMaker: public Algorithm {
  protected:
    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;
    NTuple::Item<float>             m_mcWeight0, m_mcWeight1, m_mcWeight2;
    NTuple::Item<long>              m_njet;
    NTuple::Array<long>             m_jetlabel;
    NTuple::Array<float>            m_bjetdR, m_cjetdR, m_taujetdR;  
    NTuple::Array<float>            m_pxjet, m_pyjet, m_pzjet, m_eejet;
    NTuple::Array<float>            m_ptcalo, m_ucorr;
    NTuple::Array<float>            m_btag, m_bcorr;
    NTuple::Array<float>            m_tautag, m_taucorr;
    NTuple::Array<float>            m_tau1ptag, m_tau3ptag, m_tau1p3pcorr;
    NTuple::Item<long>              m_nele;
    NTuple::Array<long>             m_codeele;
    NTuple::Array<float>            m_pxele, m_pyele, m_pzele, m_eeele;
    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;
    NTuple::Item<long>              m_npho;
    NTuple::Array<long>             m_codepho;
    NTuple::Array<float>            m_pxpho, m_pypho, m_pzpho, m_eepho;
    /*
    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;
    */
    NTuple::Item<long>              m_npart;
    NTuple::Array<long>             m_kppart, m_kspart, m_kfpart;
    NTuple::Array<long>             m_kpmoth, m_kfmoth;
    NTuple::Array<float>            m_pxpart, m_pypart, m_pzpart, m_eepart ;
    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); 
    ~StandardNtupleMaker();
    StatusCode initialize();
    StatusCode execute();
    StatusCode finalize();
    
  private:
    std::pair<int, int> getParentCodes(MCparticleCollection::const_iterator&);
    long   m_nEvents;                    // Number events processed
    TesIO* m_tesIO;
    //holds the Stotegate MC location
    std::string m_mcLocation;    
    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;  
    
    typedef std::vector<HepMC::GenParticle*>            t_MC_particleCollection; 
    typedef std::vector<HepMC::GenParticle*>::iterator  t_MC_particleIterator;
  };
  
} // end of namespace bracket  

#endif 

---