SVDPositionErrorScaleFactorImporterModule.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 <framework/core/Module.h>
#include <framework/datastore/StoreArray.h>
#include <framework/database/Database.h>
#include <framework/database/IntervalOfValidity.h>
 
#include <svd/dataobjects/SVDHistograms.h>
#include <svd/dataobjects/SVDCluster.h>
#include <svd/dataobjects/SVDTrueHit.h>
 
#include <string>
#include <TFile.h>
#include <TTree.h>
#include <TBranch.h>
#include <TH1F.h>
 
namespace Belle2 {
  class SVDPositionErrorScaleFactorImporterModule : public Module {
 
  public:
 
    SVDPositionErrorScaleFactorImporterModule();
 
    virtual void initialize() override;
 
    virtual void beginRun() override;
 
    virtual void event() override;
 
    virtual void endRun() override;
 
  private:
 
    /* user-defined parameters */
    std::string m_rootFileName;   
    std::string m_posAlgorithm = "CoGOnly"; 
    bool m_noOutliers = false; 
    float m_min = -10; 
    float m_max = 10; 
    float m_nBins = 100; 
    std::string m_uniqueID = "PositionErrorScaleFactors_TEST"; 
    StoreArray<SVDCluster> m_clusters; 
    StoreArray<SVDTrueHit> m_truehits; 
    //* ROOT file related parameters
    TFile* m_rootFilePtr = nullptr; 
    TTree* m_tree = nullptr; 
    //tree variables
    int m_exp = -1; 
    int m_run = -1; 
    int m_layer = -1; 
    int m_ladder = -1; 
    int m_sensor = -1; 
    int m_side = -1; 
    int m_size = -1; 
    float m_clsCharge = -1; 
    float m_clsTime = -1; 
    float m_clsPos = -1; 
    float m_clsErr = -1; 
    float m_clsResid = -1; 
    float m_clsPull = -1; 
    double oneSigma(TH1F*);
 
    //CLUSTER PULLS
    SVDHistograms<TH1F>* m_hClsPullSize1 = nullptr; 
    SVDHistograms<TH1F>* m_hClsPullSize2 = nullptr; 
    SVDHistograms<TH1F>* m_hClsPullSize3 = nullptr; 
    SVDHistograms<TH1F>* m_hClsPullSize4 = nullptr; 
    SVDHistograms<TH1F>* m_hClsPullSize5 = nullptr; 
    const int maxSize = 5; 
    TH1F* m_hL3Pulls[5][2] = {}; 
    TH1F* m_hFWPulls[5][2] = {}; 
    TH1F* m_hBWPulls[5][2] = {}; 
    TH1F* m_hORPulls[5][2] = {}; 
  };
}