release-05-01-25/doxygen/ECLBhabhaTCollectorModule_8h_source.html

/**************************************************************************

 * BASF2 (Belle Analysis Framework 2)                                     *

 * Copyright(C) 2019 - Belle II Collaboration                             *

 *                                                                        *

 * Author: The Belle II Collaboration                                     *

 * Contributors:                                                          *

 *    Ewan Hill                                                           *

 *    Mikhail Remnev                                                      *

 *                                                                        *

 * This software is provided "as is" without any warranty.                *

 **************************************************************************/


#ifndef ECLBHABHATCOLLECTORMODULE_H

#define ECLBHABHATCOLLECTORMODULE_H


#include <framework/core/Module.h>

#include <ecl/utility/ECLChannelMapper.h>

#include <ecl/utility/ECLTimingUtilities.h>


#include <calibration/CalibrationCollectorModule.h>

#include <framework/database/DBObjPtr.h>

#include <framework/datastore/StoreArray.h>

#include <framework/dataobjects/EventT0.h>


class TTree;


namespace Belle2 {

  class ECLDigit;

  class ECLCalDigit;

  class Track;

  class ECLCrystalCalib;

  class ECLCluster;

  class ECLChannelMapper;


  class ECLBhabhaTCollectorModule : public CalibrationCollectorModule {


  public:


    ECLBhabhaTCollectorModule();


    virtual ~ECLBhabhaTCollectorModule();


    void inDefineHisto();


    void prepare();


    void collect();


  private:


    bool m_saveTree;


    /****** Parameters END ******/


    StoreArray<Track> tracks;

    StoreObjPtr<EventT0> m_eventT0;


    DBObjPtr<ECLCrystalCalib> m_ElectronicsDB;

    std::vector<float> m_Electronics;

    DBObjPtr<ECLCrystalCalib> m_ElectronicsTimeDB;

    std::vector<float> m_ElectronicsTime;

    DBObjPtr<ECLCrystalCalib> m_FlightTimeDB;

    std::vector<float> m_FlightTime;

    DBObjPtr<ECLCrystalCalib> m_PreviousCrystalTimeDB;

    std::vector<float> m_PreviousCrystalTime;

    std::vector<float> m_PreviousCrystalTimeUnc;

    DBObjPtr<ECLCrystalCalib> m_CrateTimeDB;

    std::vector<float> m_CrateTime;

    std::vector<float> m_CrateTimeUnc;

    TTree* m_dbgTree_electrons;

    TTree* m_dbgTree_tracks;

    TTree* m_dbgTree_crystals;

    TTree* m_dbgTree_event;

    TTree* m_dbgTree_allCuts;

    TTree* m_dbgTree_evt_allCuts;

    TTree* m_dbgTree_crys_allCuts;


    /*** tree branches ***/

    /*** See inDefineHisto method for branches description ***/

    int m_tree_evtNum;

    int m_tree_cid;

    double m_tree_phi;

    double m_tree_theta;

    int m_tree_amp;

    double m_tree_en;

    double m_tree_E1Etot;

    double m_tree_E1E2;

    double m_tree_E1p;

    int m_tree_quality;

    double m_tree_timeF;

    double m_tree_time;

    double m_tree_timetsPreviousTimeCalibs;

    double m_tree_t0;

    double m_tree_t0_unc;

    double m_tree_t0_ECLclosestCDC;

    double m_tree_t0_ECL_minChi2;

    double m_tree_d0;

    double m_tree_z0;

    double m_tree_p;

    double m_tree_nCDChits;

    double m_tree_clustCrysE_DIV_maxEcrys;

    double m_tree_clustCrysE;

    double m_tree_enPlus;

    double m_tree_enNeg;

    double m_tree_tClustPos;

    double m_tree_tClustNeg;

    double m_tree_maxEcrystPosClust;

    double m_tree_maxEcrystNegClust;

    double m_tree_tClust;

    double m_tree_ECLCalDigitTime;

    double m_tree_ECLCalDigitE;

    double m_tree_ECLDigitAmplitude;

    int m_charge;

    double m_E_DIV_p;

    double m_massInvTracks;

    /*** tree branches END ***/


    StoreArray<ECLDigit> m_eclDigitArray;

    StoreArray<ECLCalDigit> m_eclCalDigitArray;

    StoreArray<ECLCluster> m_eclClusterArray;

    std::vector<float> m_EperCrys;

    std::vector<int> m_eclCalDigitID;

    std::vector<int> m_eclDigitID;

    /****** Parameters for cuts ******/

    short m_timeAbsMax;

    int m_minCrystal;

    int m_maxCrystal;

    /* d0 and z0 values of the loose and tight tracks*/

    double m_looseTrkZ0;

    double m_tightTrkZ0;

    double m_looseTrkD0;

    double m_tightTrkD0;


    int m_crystalCrate;

    int m_runNum;

    bool m_storeCalib = true;

    // For the energy dependence correction to the time

    // t-t0 = p1 + pow( (p3/(amplitude+p2)), p4 ) + p5*exp(-amplitude/p6)      ("Energy dependence equation")


    std::unique_ptr< Belle2::ECL::ECLTimingUtilities > m_ECLTimeUtil =

      std::make_unique<Belle2::ECL::ECLTimingUtilities>();

    double m_energyDependenceTimeOffsetFitParam_p1 = 0;

    double m_energyDependenceTimeOffsetFitParam_p2 = 88449.;

    double m_energyDependenceTimeOffsetFitParam_p3 = 0.20867E+06;

    double m_energyDependenceTimeOffsetFitParam_p4 = 3.1482;

    double m_energyDependenceTimeOffsetFitParam_p5 = 7.4747;

    double m_energyDependenceTimeOffsetFitParam_p6 = 1279.3;

  };

}


#endif /* ECLBHABHATCOLLECTORMODULE_H */