HadronSaturation.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <string>
#include <iostream>
#include <vector>
#include <fstream>
 
#include <TFile.h>
#include <TTree.h>
#include <TRandom.h>
#include <TMath.h>
#include <TFitter.h>
 
#include <framework/database/DBObjPtr.h>
#include <cdc/utilities/CDCDedxHadSat.h>
 
namespace Belle2 {
  class HadronSaturation {
 
  public:
 
    HadronSaturation();
 
    HadronSaturation(double alpha, double gamma, double delta, double power, double ratio, int cosbins);
 
    virtual ~HadronSaturation() { clear(); };
 
    void fillSample(TString infilename);
 
    void printEvents(int firstevent, int nevents);
 
    void fitSaturation();
 
    void clear();
 
    void setCosBins(int nbins) { m_cosbins = nbins; }
 
    double myFunction(double alpha, double gamma, double delta, double power, double ratio);
 
    static void minuitFunction(int&, double*, double& result, double* para, int);
 
  private:
 
    int m_cosbins;  
    double m_alpha; 
    double m_gamma; 
    double m_delta; 
    double m_power; 
    double m_ratio; 
    std::vector< double > m_dedx;      
    std::vector< double > m_dedxerror; 
    std::vector< double > m_betagamma; 
    std::vector< double > m_costheta;  
  };
}