release-05-02-19/doxygen/SVDROIFinderAnalysisModule_8h_source.html

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

 * BASF2 (Belle Analysis Framework 2)                                     *

 * Copyright(C) 2011 - Belle II Collaboration                             *

 *                                                                        *

 * Author: The Belle II Collaboration                                     *

 * Contributors: Giulia Casarosa                                          *

 *                                                                        *

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

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


#pragma once


#include <framework/core/Module.h>

#include <mdst/dataobjects/MCParticle.h>

#include <tracking/dataobjects/RecoTrack.h>

#include <tracking/dataobjects/ROIid.h>

#include <tracking/dataobjects/SVDIntercept.h>

#include <svd/dataobjects/SVDShaperDigit.h>

#include <string>

#include <TFile.h>

#include <TH1F.h>

#include <TH2F.h>

#include <TGraphErrors.h>


namespace Belle2 {

  class SVDROIFinderAnalysisModule : public Module {


  public:


    SVDROIFinderAnalysisModule();


    ~SVDROIFinderAnalysisModule();


    void initialize() override;


    void beginRun() override;


    void event() override;


    void endRun() override;


    void terminate() override;


  private:


    StoreArray<SVDShaperDigit> m_shapers;

    StoreArray<ROIid> m_ROIs;

    StoreArray<RecoTrack> m_trackList;

    StoreArray<SVDIntercept> m_SVDIntercepts;

    StoreArray<MCParticle> m_mcParticles;

    bool m_isSimulation;

    std::string m_shapersName;

    std::string m_recoTrackListName;

    std::string m_SVDInterceptListName;

    std::string m_ROIListName;

    TFile* m_rootFilePtr;

    std::string m_rootFileName;

    bool m_writeToRoot;

    int m_rootEvent;

    //graphs & histos - results

    Double_t pt[6]  = {0.05, 0.15, 0.25, 0.4, 0.75, 1.5};

    Double_t ptErr[6] = { 0.05, 0.05, 0.05, 0.1, 0.25, 0.5};

    TGraphErrors* m_gEff2;

    TGraphErrors* m_gEff;

    TH1F* m_h1DigitsPerParticle;

    TH1F* m_h1RecoTracksPerParticle;

    TH1F* m_h1digiIn;

    TH1F* m_h1digiOut2;

    TH1F* m_h1digiOut3;

    TH1F* m_h1digiOut4;

    TH1F* m_h1digiOut5;

    TH1F* m_h1TrackOneDigiIn;

    TH1F* m_h1nnotINtrack2;

    TH1F* m_h1nnotINtrack3;

    TH1F* m_h1nnotINtrack4;

    TH1F* m_h1nnotINtrack5;

    TH1F* m_h1Track;

    TH1F* m_h1Track_pt;

    TH1F* m_h1Track_phi;

    TH1F* m_h1Track_lambda;

    TH1F* m_h1Track_cosTheta;

    TH1F* m_h1Track_pVal;

    TH1F* m_h1Track_nSVDhits;

    TH1F* m_h1Track_nCDChits;

    TH1F* m_h1INtrack1;

    TH1F* m_h1INtrack1_pt;

    TH1F* m_h1INtrack1_phi;

    TH1F* m_h1INtrack1_lambda;

    TH1F* m_h1INtrack1_cosTheta;

    TH1F* m_h1INtrack1_pVal;

    TH1F* m_h1INtrack1_nSVDhits;

    TH1F* m_h1INtrack1_nCDChits;

    TH1F* m_h1notINtrack5;

    TH1F* m_h1notINtrack5_pt;

    TH1F* m_h1notINtrack5_phi;

    TH1F* m_h1notINtrack5_lambda;

    TH1F* m_h1notINtrack5_cosTheta;

    TH1F* m_h1notINtrack5_pVal;

    TH1F* m_h1notINtrack5_nSVDhits;

    TH1F* m_h1notINtrack5_nCDChits;

    //fill digits inside ROIs

    TH1F* m_h1PullU;

    TH1F* m_h1PullV;

    TH2F* m_h2sigmaUphi;

    TH2F* m_h2sigmaVphi;

    TH1F* m_h1ResidU;

    TH1F* m_h1ResidV;

    TH1F* m_h1SigmaU;

    TH1F* m_h1SigmaV;

    TH1F* m_h1GlobalTime;

    //fill digits outside2 ROIs

    TH2F* m_h2sigmaUphi_out2;

    TH2F* m_h2sigmaVphi_out2;

    TH1F* m_h1ResidU_out2;

    TH1F* m_h1ResidV_out2;

    TH1F* m_h1SigmaU_out2;

    TH1F* m_h1SigmaV_out2;

    TH1F* m_h1GlobalTime_out2;

    //fill digits outside3 ROIs

    TH2F* m_h2sigmaUphi_out3;

    TH2F* m_h2sigmaVphi_out3;

    TH1F* m_h1ResidU_out3;

    TH1F* m_h1ResidV_out3;

    TH1F* m_h1SigmaU_out3;

    TH1F* m_h1SigmaV_out3;

    TH1F* m_h1GlobalTime_out3;

    //fill digits outside4 ROIs

    TH2F* m_h2sigmaUphi_out4;

    TH2F* m_h2sigmaVphi_out4;

    TH1F* m_h1SigmaU_out4;

    TH1F* m_h1SigmaV_out4;

    TH1F* m_h1GlobalTime_out4;

    //fill digits outside5 ROIs

    TH1F* m_h1GlobalTime_out5;

    //ROI stuff

    TH2F* m_h2ROIbottomLeft;

    TH2F* m_h2ROItopRight;

    TH2F* m_h2ROIuMinMax;

    TH2F* m_h2ROIvMinMax;

    TH1F* m_h1totROIs;

    TH1F* m_h1okROIs;

    TH1F* m_h1totUstrips;

    TH1F* m_h1totVstrips;

    TH1F* m_h1effPerTrack;

    //variables

    double m_globalTime;

    double m_coorU;

    double m_coorV;

    double m_sigmaU;

    double m_sigmaV;

    int m_vxdID;

    double m_coormc;

    int m_idmc;

    int m_vxdIDmc;

    double m_pTmc;

    double m_momXmc;

    double m_momYmc;

    double m_momZmc;

    double m_thetamc;

    double m_costhetamc;

    double m_phimc;

    double m_lambdamc;

    unsigned int Ntrack;

    unsigned int NtrackHit;

    unsigned int n_notINtrack2;

    unsigned int n_notINtrack3;

    unsigned int n_notINtrack4;

    unsigned int n_notINtrack5;

    unsigned int n_rois;

    unsigned int n_OKrois;

    unsigned int m_nGoodROIs;

    unsigned int n_intercepts;

    unsigned int n_tracks;

    unsigned int n_tracksWithDigits;

    unsigned int n_tracksWithDigitsInROI;

    unsigned int n_svdDigit;

    unsigned int n_svdDigitInROI;

    unsigned int n_notINdigit2;

    unsigned int n_notINdigit3;

    unsigned int n_notINdigit4;

    unsigned int n_notINdigit5;

    unsigned int nsvdDigit[6];

    unsigned int nsvdDigitInROI[6];

    unsigned int nnotINdigit2[6];

    unsigned int nnotINdigit3[6];

    unsigned int nnotINdigit4[6];

    unsigned int nnotINdigit5[6];

    unsigned int TrackOneDigiIn[6];

    unsigned int nnotINtrack2[6];

    unsigned int nnotINtrack3[6];

    unsigned int nnotINtrack4[6];

    unsigned int nnotINtrack5[6];

  };


}