release-05-02-19/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/EventMetaData.h>

#include <framework/dataobjects/EventT0.h>


class TTree;


namespace Belle2 {

  class ECLDigit;

  class ECLCalDigit;

  class Track;

  class ECLCrystalCalib;

  class ECLReferenceCrystalPerCrateCalib;

  class ECLCluster;

  class ECLChannelMapper;


  const static double realNaN = std::numeric_limits<double>::quiet_NaN();

  const static int    intNaN  = std::numeric_limits<int>::quiet_NaN();

  class ECLBhabhaTCollectorModule : public CalibrationCollectorModule {


  public:


    ECLBhabhaTCollectorModule();


    virtual ~ECLBhabhaTCollectorModule();


    void inDefineHisto() override;


    void prepare() override;


    void collect() override;


  private:


    bool m_saveTree;


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


    StoreArray<Track> tracks;

    StoreObjPtr<EventT0> m_eventT0;


    StoreObjPtr<EventMetaData> m_EventMetaData;


    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;

    DBObjPtr<ECLReferenceCrystalPerCrateCalib> m_RefCrystalsCalibDB;

    std::vector<short> m_RefCrystalsCalib;

    TTree* m_dbgTree_electrons = nullptr;

    TTree* m_dbgTree_tracks = nullptr;

    TTree* m_dbgTree_crystals = nullptr;

    TTree* m_dbgTree_event = nullptr;

    TTree* m_dbgTree_allCuts = nullptr;

    TTree* m_dbgTree_evt_allCuts = nullptr;

    TTree* m_dbgTree_crys_allCuts = nullptr;

    /*** 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 = realNaN;

    double m_tree_t0_unc = realNaN;

    double m_tree_t0_ECLclosestCDC = realNaN;

    double m_tree_t0_ECL_minChi2 = realNaN;

    double m_tree_d0 = realNaN;

    double m_tree_z0 = realNaN;

    double m_tree_p = realNaN;

    double m_tree_nCDChits = realNaN;

    double m_tree_clustCrysE_DIV_maxEcrys = realNaN;

    double m_tree_clustCrysE = realNaN;

    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 = realNaN;

    double m_tightTrkZ0 = realNaN;

    double m_looseTrkD0 = realNaN;

    double m_tightTrkD0 = realNaN;

    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_hadronEventT0_TO_bhabhaEventT0_correction;


  };

}


#endif /* ECLBHABHATCOLLECTORMODULE_H */