NtupleReader/NtupleReader/MyTauJet.h

//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// 20.11.2006, AUTHOR: THIES EHRICH
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

#ifndef MyTauJetH
#define MyTauJetH

//:::::::::::::::::::::::::
//:: CLASS MyTauJet ::
//:::::::::::::::::::::::::


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

// standard C++ //
#include <iostream>
#include <string>
#include <vector>
#include "math.h"

// ROOT //
#include "TLorentzVector.h"

// MyParticle base class //
#include "MyParticle.h"
#include "NtupleReader/MyDatasetInfo.h"

class MyTauJet : public MyParticle {

public:
    // Constructors //
    MyTauJet(void) {
        init( TLorentzVector(0.0, 0.0, 0.0, 0.0), 
              0.0, 0, 0, 0.0, 0.0, 0.0, 0.0, 0.0, 
              0.0, 0.0 ,0.0, 0.0 ,0.0, 0.0, 0,
              TLorentzVector(0.0,0.0,0.0,0.0),
              TLorentzVector(0.0,0.0,0.0,0.0),
              TLorentzVector(0.0,0.0,0.0,0.0),
              0,0,0,0,0,0);
    }
    MyTauJet(const TLorentzVector & p,
             const double & charge,
             const int & rec_flag,
             const int & index,
             const double & emRadius,
             const double & isolationFraction,
             const double & centralityFraction,
             const double & stripWidth2,
             const double & nStripCells,
             const double & llh,
             const double & lowPtTauEleDiscriminant,
             const double & tauENeuralNetwork,
             const double & tauJetNeuralnetwork,
             const double & etHadCalib,
             const double & etEMCalib,
             const int & nb_tracks,
             const TLorentzVector & track_1,
             const TLorentzVector & track_2,
             const TLorentzVector & track_3,
             const int & track_TRTHits_1,
             const int & track_TRTHighThresholdHits_1,
             const int & track_TRTHits_2,
             const int & track_TRTHighThresholdHits_2,
             const int & track_TRTHits_3,
             const int & track_TRTHighThresholdHits_3
             ){

        init(p, charge, rec_flag, index, emRadius,
             isolationFraction,centralityFraction,stripWidth2,nStripCells,
             llh, lowPtTauEleDiscriminant, tauENeuralNetwork, tauJetNeuralnetwork,
             etHadCalib, etEMCalib, nb_tracks, track_1, track_2, track_3,
             track_TRTHits_1, track_TRTHighThresholdHits_1,
             track_TRTHits_2, track_TRTHighThresholdHits_2,
             track_TRTHits_3, track_TRTHighThresholdHits_3);
    }

    // Methods //
    // get methods //
    inline double llh(void) const;
    inline double emRadius(void) const;
    inline double isolationFraction(void) const;
    inline double centralityFraction(void) const;
    inline double stripWidth2(void) const;
    inline double nStripCells(void) const;
    inline double lowPtTauEleDiscriminant(void) const;
    inline double tauENeuralNetwork(void) const;
    inline double tauJetNeuralnetwork(void) const;
    inline int reconstruction_flag(void) const;

    inline double ethadcalib(void) const;
    inline double etemcalib(void) const;
    inline int nb_tracks(void) const;
    inline TLorentzVector track_1(void) const;
    inline TLorentzVector track_2(void) const;
    inline TLorentzVector track_3(void) const;
    inline TLorentzVector tlv_corr(void) const;
    int track_TRTHits_1(void) const;
    int track_TRTHighThresholdHits_1(void) const;
    int track_TRTHits_2(void) const;
    int track_TRTHighThresholdHits_2(void) const;
    int track_TRTHits_3(void) const;
    int track_TRTHighThresholdHits_3(void) const;

    // set methods //
    inline void set_reconstruction_flag(const int & flag);

    inline MyTauJet* Clone(void) const;

    inline void PrintParticle(std::string option = "");

    private:
    
    // The MyDatasetInfo object 
    MyDatasetInfo *m_dataset_info;
    // taujet data //
    int m_rec_flag; // reconstruction flag (see the definition above)
    double m_emRadius;
    double m_isolationFraction;
    double m_centralityFraction;
    double m_stripWidth2;
    double m_nStripCells;
    double m_llh;       // logLikelihoodRatio from tauRec
    double m_lowPtTauEleDiscriminant;   // 
    double m_tauENeuralNetwork; // 
    double m_tauJetNeuralnetwork;       // 
    double m_etHadCalib;        //  tauRec
    double m_etEMCalib; //  tauRec
    int m_nb_tracks;    //  nb of charged tracks from taujets (pions)
    TLorentzVector m_track_1;   //  pt of first track
    TLorentzVector m_track_2;   //  pt of second track
    TLorentzVector m_track_3;   //  pt of third track
    int m_track_TRTHits_1;
    int m_track_TRTHighThresholdHits_1;
    int m_track_TRTHits_2;
    int m_track_TRTHighThresholdHits_2;
    int m_track_TRTHits_3;
    int m_track_TRTHighThresholdHits_3;

    // private methods //
    inline void init(const TLorentzVector & p,
                     const double & charge,
                     const int & rec_flag,
                     const int & index,
                     const double & emRadius ,
                     const double & isolationFraction,
                     const double & centralityFraction,
                     const double & stripWidth2,
                     const double & nStripCells,
                     const double & llh,
                     const double & lowPtTauEleDiscriminant,
                     const double & tauENeuralNetwork,
                     const double & tauJetNeuralnetwork,
                     const double & etHadCalib,
                     const double & etEMCalib,
                     const int & nb_tracks,
                     const TLorentzVector & track_1,
                     const TLorentzVector & track_2,
                     const TLorentzVector & track_3,
                     const int & track_TRTHits_1,
                     const int & track_TRTHighThresholdHits_1,
                     const int & track_TRTHits_2,
                     const int & track_TRTHighThresholdHits_2,
                     const int & track_TRTHits_3,
                     const int & track_TRTHighThresholdHits_3
                     ) ;

    // initialization method;
    // tcharge: charge of the taujet;
    // rec_flag: reconstruction flag:
    //    = 0: unknown algorithm,
    //    = 1: cone-4 algorithm,
    //    = 2: cone-7 algorithm,
    //    = 3: kt algorithm;
    // index = index in original vector
    // llh = likelihood value
    // 

};

//::::::::::::::::::::::::::::::::::::::
//:: IMPLEMENTATION OF INLINE METHODS ::
//::::::::::::::::::::::::::::::::::::::

#include "MyTauJet.ixx"

#endif

---