include/N3AnaBase/measurements/tof_tools.h

#ifndef __N3AnaBase__tof_tools_h
#define __N3AnaBase__tof_tools_h 1

#include <N3AnaBase/event/ana_event.h>
#include <cmath>
#include <string>
#include <map>
#include <sstream>

#include <N3GeomDef/N3Vec3.h>
#include <N3GeomDef/N3SourceGeom.h>
#include <N3GeomDef/N3ScinGeom.h>


#include <gsl/gsl_cdf.h>

namespace nemo3 {

  using namespace std;


  class tof_tools
  {

  private :
    static const double kC;    
    static const double kMe;   
    static const double kSig;  
    static double kTauSc;      
    static double kTauSc_g;    
    static double kTauEl;      
    static double kVg;         
    static double fwhm_trans;  
  public :
    tof_tools  ();
    tof_tools  ( bool debug_, bool use_energy_loss_ = false, bool timeerr_ = false, bool fluct_ = false );
~tof_tools ();

    void init_from (const datatools::utils::properties &); 

    void set_debug      ( bool debug_ ); 
    void set_run_number ( int run_number_, std::string runtype_ = "" ); 
    void set_ECorr      ( bool energy_loss_ ); 
    void set_TimeErr    ( bool timeerr_ );
    void set_LTFluct    ( bool fluct_ );
    void set_simulation ( bool simu_);
    void set_nodb4simu  ( bool nodb_);
    void set_tau_sc     ( double ts_ ) {kTauSc=ts_;}
    void set_tau_sc_g   ( double ts_ ) {kTauSc_g=ts_;}
    void set_tau_el     ( double ts_ ) {kTauEl=ts_;}
    void set_Vg         ( double vg_ ) {kVg=vg_;}
    void set_run_type   ( std::string rt_ ) {__run_type=rt_;}
    bool get_simulation ();
    string get_run_type () {return __run_type;}

  private:
    float get_tse (int sector_, int iobt_, int fcll_, int block_); 
    float get_ltce (int sector_, int iobt_, int fcll_, int block_); 
    float get_ltfluct (int sector_, int iobt_, int fcll_, int block_); 


  private :
    double Beta ( double ke_, double mass_ ); 

    double T_th ( double ke_, double mass_, double trlen_ ); 

    double ECorrTracking ( double energy_, double trlen_, scin_hit_id & id_ ); 

    double SigmaE ( scin_hit_id & id_, double ke_ ); 

    double SigmaTrlen ( double mass_ ); 
    double SigmaT ( scin_hit_id & id_, double ke_, bool gamma_=false ); 

    double FWHM ( scin_hit_id & id_ ); 
    double g_trlen ( geomtools::vector_3d & vertex_, scin_hit_id & id_ ); 
    double g_trlen ( scin_hit_id& id1_, scin_hit_id& id2_); 

  public :
    bool get_dt (const track & track1_,const track & track2_,
                 double & dt_int_, double & dt_ext_); 

    bool get_dt (const track & track1_,const scin_hit & scin_hit2_,
                 double & dt_int_, double & dt_ext_, bool change_elec_energy_=false,
                 double elec_energy_=0, double elec_length_=0); 

    bool get_dt (const scin_hit & scin_hit1_,const scin_hit & scin_hit2_,
                   double & dt_);

    bool get_chi2 (const track & track1_,const track & track2_,
                   double & internal_chi2_, double & external_chi2_,
                   bool keep_dt_=false); 

    bool get_chi2 (const track & track_,const scin_hit & scin_hit_,
                   double & internal_chi2_, double & external_chi2_,
                   bool keep_dt_=false, bool change_elec_energy_=false,
                   double elec_energy_=0, double elec_length_=0); 

    bool get_chi2 (const scin_hit & scin_hit1_,const scin_hit & scin_hit2_,
                   double &chi2_,
                   bool keep_dt_=false); 

    bool get_chi2_1e (const scin_hit & scin_hit1_,const scin_hit & scin_hit2_,
                      double & internal_chi2_, double & internal_prob_);

    bool get_chi2_2e_1g (const track & track1_, const track & track2_,
                         const scin_hit & scin_hit3_, const bool switch_,
                         double & internal_chi2_, double & internal_prob_);

    bool get_chi2_2e_2g (const track & track1_, const track & track2_,
                         const scin_hit & scin_hit3_, const scin_hit & scin_hit4_,
                         const bool switch_, double & internal_chi2_, double & internal_prob_);

    bool get_chi2_2g (const track & track_, const scin_hit & scin_hit2_,
                      const scin_hit & scin_hit3_, const bool switch_,
                      double & internal_chi2_, double & internal_prob_,
                      bool change_elec_energy_=false, double elec_energy_=0,
                      double elec_length_=0);

    bool get_chi2_3g (const track & track_, const scin_hit & scin_hit2_,
                      const scin_hit & scin_hit3_, const scin_hit & scin_hit4_,
                      const bool switch_, double & internal_chi2_, double & internal_prob_,
                      bool change_elec_energy_=false, double elec_energy_=0,
                      double elec_length_=0);

    static void get_proba (double  internal_chi2_, double  external_chi2_,
                    double & internal_prob_, double & external_prob_,
                           unsigned int freedom_=1); 

    void set_vera (bool vera_){__vera_method=vera_;} 

  private:
    bool   __debug;            
    bool   __use_energy_loss;  
    bool   __use_ltfluct;      
    bool   __use_time_err;     
    bool   __simulation;       
    bool   __vera_method;      
    bool   __nodb4simu;        
    int    __run_number;       
    double __t1;               
    double __t2;               
    double __t_th1_int;        
    double __t_th2_int;        
    double __t_th1_ext;        
    double __t_th2_ext;        
    double __dt_int;           
    double __dt_ext;           
    double __dt;               
    double __e1;               
    double __e2;               
    double __e_last;           
    string __run_type;         


  };
} // end of namespace nemo3

#endif // __N3AnaBase__tof_tools_h

// end of tof_tools.h LvAUKAks