TRGECLModule.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 <framework/core/Module.h>
#include <framework/datastore/StoreArray.h>
#include <framework/database/DBArray.h>
#include <framework/database/DBObjPtr.h>
 
#include "trg/ecl/dataobjects/TRGECLHit.h"
#include "trg/ecl/dataobjects/TRGECLTrg.h"
#include "trg/ecl/dataobjects/TRGECLCluster.h"
#include "trg/ecl/dbobjects/TRGECLETMPara.h"
 
#include <map>
 
namespace Belle2 {

  class TRGECLModule : public Module {
 
 
  public:

    TRGECLModule();

    virtual ~TRGECLModule();

    virtual void initialize() override;

    virtual void beginRun() override;

    virtual void event() override;

    virtual void endRun() override;

    virtual void terminate() override;
 
  public:

    std::string version(void) const;
 
  private: 

    int m_debugLevel = 0;

    std::string m_configFilename;

    double getDBparmap(const std::map<std::string, double>, std::string, double);
 
  protected:

    std::string m_inColName;
    std::string m_eclHitOutColName;
    std::string m_eclTCHitOutColName;

    int m_hitNum = 0;
    int m_hitTCNum = 0;
 
  private:
 
    // double m_timeCPU;        //  CPU time
    int m_nRun = 0;
    int m_nEvent = 0;
    int m_Clustering = 1;
    int m_ClusterLimit = 6;
    int m_EventTiming = 1;
    double m_TimeWindow = 250.0;
    double m_OverlapWindow = 125.0;
    int m_NofTopTC = 3;
    int m_SelectEvent = 1;
    bool m_ConditionDB = true;
    double m_ADCtoEnergy = 0.0;
    std::vector<double> m_TotalEnergy;
    std::vector<double> m_2DBhabhaThresholdFWD;
    std::vector<double> m_2DBhabhaThresholdBWD;
    std::vector<double> m_3DBhabhaVetoThreshold;
    std::vector<int> m_3DBhabhaVetoAngle;
    std::vector<double> m_3DBhabhaSelectionThreshold;
    std::vector<int> m_3DBhabhaSelectionAngle;
    std::vector<int> m_3DBhabhaSelectionPreScale;
    double m_mumuThreshold = 0.0;
    std::vector<int> m_mumuAngle;
    std::vector<double> m_lmlCLELabCut;
    std::vector<double> m_lmlCLECMSCut;
    int m_lml00NCLforMinE = 0;
    int m_lml12NCLforMinE = 0;
    int m_lml13ThetaIdSelection = 0;
    double m_ECLBurstThreshold = 0.0;
    std::vector<double> m_EventTimingQualityThreshold;
    std::vector<int> m_3DBhabhaVetoInTrackThetaRegion;
    std::vector<int> m_taub2bAngleCut;
    double m_taub2bEtotCut = 0.0;
    double m_taub2bCLELabCut = 0.0;
    std::vector<int> m_hie12BhabhaVetoAngle;
    double m_taub2b2EtotCut = 0.0;
    std::vector<int> m_taub2b2AngleCut;
    std::vector<double> m_taub2b2CLELabCut;
    double m_taub2b3EtotCut = 0.0;
    std::vector<int> m_taub2b3AngleCut;
    double m_taub2b3CLEb2bLabCut = 0.0;
    std::vector<double> m_taub2b3CLELabCut;
    double m_hie4LowCLELabCut = 0.0;

    TrgEclMaster* etm;

    StoreArray<TRGECLHit> m_TRGECLHit;
    StoreArray<TRGECLTrg> m_TRGECLTrg;
    StoreArray<TRGECLCluster> m_TRGECLCluster;
 
  };

} // namespace Belle2