NtupleReader/NtupleReader/MyTools.h

//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
// 13.12.2006, AUTHOR: all
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

#ifndef MyToolsH
#define MyToolsH

//:::::::::::::::::::
//:: CLASS MyTools ::
//:::::::::::::::::::


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

// STL //
#include <iostream> 
//#include <fstream> 

//#include <string>
#include <vector>

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

// HiggsAnalysis //
#include "MyParticle.h"
#include "MyTauJet.h"
#include "MyParticleJet.h"

class MyTools {
    
public:
    // Constructors //
    MyTools(void) {
        init();
    }
    
    ~MyTools(void) {
        destruct();
    }
    
  // Methods //

    // change weight of jet
    double change_btag_weight(double current_btag_weight, double delta_weight=0);
    //recalculated all btag-weight of jets with pdgID 'flavour' a particle vector 'jets'
    //    You have to do your weight cut on the changed vector!!!
    void scale_b_tags( std::vector<MyParticle*> * jets, int flavour, double delta_bweight=0);
    
    //recalib tau et Joern GK 15.5.2007//
    double tauRec_ETcorr(double et_in, double eta_in, int ntrk);
    //input: et_in =ethadcalib+etemcalib => Et of tau candidate
    //       eta_in                      => eta of tau candidate    
    //       ntrk                        => # of assosiated tracks
    //return = recalibrated et, llh,eta,phi not changing.
    //recalib ALL taus//
    std::vector<MyParticle*> correct_tau_et( std::vector<MyParticle*> tau_vector);
    
    // matching // 

    MyParticle* find_matching_particle( MyParticle* particle, 
                                        std::vector<MyParticle*> & particle_vector,
                                        double deltaR_cut = 0.1,
                                        bool bool_matchingflag = true,
                                        bool bool_usematching = true); 
    MyParticle* find_matching_particle( TLorentzVector particle_tlv, 
                                        std::vector<MyParticle*> & particle_vector,
                                        double deltaR_cut = 0.1,
                                        bool bool_matchingflag = true,
                                        bool bool_usematching = true); 
        
    Int_t nb_multiple_matches(void);

    void clear_nb_multiple_matches(void);

    void clear_vector(std::vector<MyParticle*> & particle_vector);


    // overlap removal //

    std::vector<MyParticle*> 
    remove_overlap( std::vector<MyParticle*> vec_part_tb_rm,
                    std::vector<MyParticle*> vec_part_rm_with,
                    double deltaR_cut);


    std::vector<MyParticle*>
    delete_doubles(std::vector<MyParticle*> p_vector, double delta_r = 0.1);


    //check if as particle is stored twice   .
    //std::vector<MyParticle*> remove_own_overlap(std::vector<MyParticle*> particle_vector, 
        //                                          double radius );


    // vector sorting //

    std::vector<MyParticle*>
    sort_particles(std::vector<MyParticle*> p_vector);

    void check_order(std::vector<MyParticle*> p_vector);
    //append one vector to another
    std::vector<MyParticle*> append_vector(std::vector<MyParticle*> first_vector, 
                                           std::vector<MyParticle*> second_vector);


    // pt and eta cut 
    
    std::vector<MyParticle*> 
    pt_eta_cut( const std::vector<MyParticle*> p_vector,
                const double pt_min,  const double pt_max,
                const double eta_min, const double eta_max);

    // et and eta cut 
    
    std::vector<MyParticle*> 
    et_eta_cut( const std::vector<MyParticle*> p_vector,
                const double et_min,  const double et_max,
                const double eta_min, const double eta_max);

    std::vector<MyParticle*> 
    particles_from_vertex( const std::vector<MyParticle*> & p_vector, 
                           const int vertex);
    


    static void PrintVector( std::vector<MyParticle*> p_vector, 
                             std::string option = "");

private:
  // private variables and methods //
    
    // matching // 
    Int_t m_nb_multiple_matches; // nb of multiple matches
    
    void init();
         // initialization method (variables as defined
         // in the constructor of this class)
    void destruct(void);
         // destruction method
    
};

#endif

---