Histo_BW.h

Go to the documentation of this file.

#ifndef HISTO_BW_LINESHAPE_HH
#define HISTO_BW_LINESHAPE_HH
// author: Philippe d'Argent (p.dargent@cern.ch)
// status: 19 March 2015 GMT
#include <complex>

#include "TH1D.h"
#include "TFile.h"

#include "Mint/ILineshape.h"
#include "Mint/BW_BW.h"
#include "Mint/NamedParameter.h"
#include "Mint/NamedParameterBase.h"
#include "Mint/ResonancePropertiesList.h"
#include "Mint/FitParDependent.h"
#include "ResonancePropertiesFitRef.h"
#include "ResonanceProperties.h"

using namespace MINT;
using namespace std;

class Histo_BW  : public BW_BW, virtual public ILineshape
{
 public:
  Histo_BW( const AssociatedDecayTree& tree, const std::string& namePrefix )
    : BW_BW(tree, namePrefix)
    , _runningMassHist(0)
    , _runningWidthHist(0)
  {
    if( !std::getenv("MINT2") ){
      std::cout
    << "ERROR: Mint2 environment variable not set."
    << "       Energy-dependent mass/width histograms cannot be found."
    << std::endl;
      exit(1);
    }

    const std::string hist_path =
      static_cast<std::string>(std::getenv("MINT2"))+"/share/";

    NamedParameter<string> fileNameMass(("RunningMass_"+ namePrefix+ (BW_BW::resonanceProperties()->nameFromPid(abs(mumsPID())))),(hist_path+"running_mass/"+"RunningMass_"+ (BW_BW::resonanceProperties()->nameFromPid(abs(mumsPID()))) + ".root"));
    NamedParameter<string> fileNameWidth(("RunningWidth_"+ namePrefix + (BW_BW::resonanceProperties()->nameFromPid(abs(mumsPID())))),(hist_path+"running_width/"+"RunningWidth_"+ (BW_BW::resonanceProperties()->nameFromPid(abs(mumsPID()))) + ".root"));

    TFile* const f_m = TFile::Open(((string)fileNameMass).c_str());
    _runningMassHist = get_mass_histo(f_m, "RunningMass");
    TFile* const f_w = TFile::Open(((string)fileNameWidth).c_str());
    _runningWidthHist = get_width_histo(f_w, "RunningWidth");
  }

  virtual ~Histo_BW(){}

  virtual std::string name() const
  { return "Histo_BW("+prefix()+_theDecay.oneLiner() +")"; }

 protected:
  TH1D* _runningMassHist;
  TH1D* _runningWidthHist;
  double runningMass2();
  virtual double GofM(); 
  virtual std::complex<double> BreitWigner();
  TH1D* get_width_histo(TFile* f, const std::string& hname); 
  TH1D* get_mass_histo(TFile* f, const std::string& hname); 

  TH1D* producePhaseSpaceHist();
};

#endif //HISTO_BW_LINESHAPE_HH