PIDNtupleModule.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2010 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Marko Staric, Jake Bennett                               *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef PIDNTUPLEMODULE_H
#define PIDNTUPLEMODULE_H
 
#include <framework/core/Module.h>
#include <string>
#include <reconstruction/modules/PIDNtuple/PIDTree.h>
#include "TFile.h"
 
namespace Belle2 {
  class PIDNtupleModule : public Module {
 
  public:
 
    PIDNtupleModule();
 
    virtual ~PIDNtupleModule();
 
    virtual void initialize() override;
 
    virtual void beginRun() override;
 
    virtual void event() override;
 
    virtual void endRun() override;
 
    virtual void terminate() override;
 
    void printModuleParams() const;
 
  private:
 
    std::string m_outputFileName; 
    bool m_makeFlat; 
    double m_p1; 
    double m_p2; 
    double m_pMax; 
    double m_norm;  
    double m_value; 
    TFile* m_file;      
    TTree* m_tree;      
    PID::PIDTree m_pid;      
    double momDistribution(double p) const
    {
      double f = 1.0 - exp(-p / m_p1);
      return m_norm * f * f  * f * f * exp(-p / m_p2);
    };
 
  };
 
} // Belle2 namespace
 
#endif