CDCCosmicAnalysisModule.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 <string>
 
#include <framework/datastore/StoreArray.h>
#include <framework/datastore/StoreObjPtr.h>
#include <framework/dataobjects/EventMetaData.h>
#include <framework/dataobjects/EventT0.h>
#include <mdst/dataobjects/Track.h>
#include <mdst/dataobjects/TrackFitResult.h>
#include <tracking/dataobjects/RecoTrack.h>
#include <TTree.h>
#include <TFile.h>
#include <TVector3.h>
 
namespace Belle2 {
 
  namespace CDC {
 
 
 

 
    class CDCCosmicAnalysisModule : public Module {
 
    public:
      CDCCosmicAnalysisModule();

      virtual ~CDCCosmicAnalysisModule() override;
 

      void initialize() override;

      void beginRun() override;
      void event() override;
      void terminate() override;
 
    private:
      StoreObjPtr<EventMetaData> m_EventMetaData;   
      StoreObjPtr<EventT0> m_eventTimeStoreObject;  
      StoreArray<Track> m_Tracks;                   
      StoreArray<RecoTrack> m_RecoTracks;           
      StoreArray<TrackFitResult> m_TrackFitResults; 
      std::string m_trackArrayName;           
      std::string m_recoTrackArrayName ;       
      std::string m_trackFitResultArrayName;  
      std::string m_relRecoTrackTrackName;     
      std::string m_outputFileName;           
      std::string m_treeName;    
 
 
      TFile* tfile = nullptr;  
      TTree* tree = nullptr;             
      int    run = 0;              
      double ndf1 = 0.0;             
      double ndf2 = 0.0;             
      double Pval1 = 0.0;            
      double Pval2 = 0.0;            
      double D01 = 0.0;              
      double D02 = 0.0;              
      double Phi01 = 0.0;            
      double Phi02 = 0.0;            
      //      double Om1;              /**< Omega of 1st track. */
      //      double Om2;              /**< Omega of 2nd track. */
      double Z01 = 0.0;              
      double Z02 = 0.0;              
      double tanLambda1 = 0.0;       
      double tanLambda2 = 0.0;       
      double eD01 = 0.0;             
      double eD02 = 0.0;             
      double eOm1 = 0.0;             
      double eOm2 = 0.0;             
      double ePhi01 = 0.0;           
      double ePhi02 = 0.0;           
      double eZ01 = 0.0;             
      double eZ02 = 0.0;             
      double etanL1 = 0.0;           
      double etanL2 = 0.0;           
      double Pt1 = 0.0;              
      double Pt2 = 0.0;              
 
      TVector3 posSeed1;                
      TVector3 posSeed2;                
 
      TVector3 Mom1;                    
      TVector3 Mom2;                    
      double Omega1 = 0.0;                    
      double Omega2 = 0.0;                    
      double evtT0 = 0.0;                       
      short charge = 0;                     
 
      bool m_bField = false; 
      bool m_eventT0Extraction = false;
      bool m_phi0InRad = false; 
      bool m_storeTrackParErrors = false;
    };
  }
}