CsIStudyModule.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2014 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: beaulieu                                                 *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef CSISTUDYMODULE_H
#define CSISTUDYMODULE_H
 
#include <framework/core/HistoModule.h>
#include <framework/datastore/StoreArray.h>
#include <beast/csi/dataobjects/CsiSimHit.h>
#include <beast/csi/dataobjects/CsiHit.h>
#include <beast/csi/dataobjects/CsiDigiHit.h>
 
//C++/C standard lib elements.
#include <string>
 
 
//ROOT
#include <TH1F.h>
 
namespace Belle2 {
  class CsIStudyModule : public HistoModule {
 
  public:
 
    CsIStudyModule();
 
    virtual ~CsIStudyModule();
 
    virtual void initialize() override;
 
    virtual void beginRun() override;
 
    virtual void endRun() override;
 
 
    virtual void event() override;
 
    virtual void terminate() override;
 
 
    /*
       Histogram definitions such as TH1(), TH2(), TNtuple(), TTree().... are supposed
       to be placed in this function.
    */
    virtual void defineHisto() override;
 
 
    bool* readDPPStatusBits(char packedbits);
 
  private:
 
    StoreArray<CsiSimHit> m_simhits; 
    StoreArray<CsiHit> m_hits; 
    StoreArray<CsiDigiHit> m_digihits; 
    // Parameters
    double m_paramTemplate;  
    int m_nWFSamples
    ; 
    std::string  m_waveformFilename ; 
    //Othe rmembers
    TFile* fWF; 
    // Histograms
    TH1F* h_CrystalEdep;        
    TH1F* h_CrystalSpectrum;    
    TH1F* h_CrystalRadDose;     
    TH1F* h_CrystalRadDoseSH;   
    TH1F* h_NhitCrystal;        
    TH1F* h_LightYieldCrystal;  
    TH1S* h_Waveform;           
    TH1C* h_Gate;               
    TH1F* h_Charge ;            
    TH1F* h_TrueEdep;           
    TH1F* h_Height;           
    TH1F* h_Weight;           
    // Constants
    const double GeVtoJ = 1.6e-10; 
    const double Sampletime = 1000.0; 
    const double usInYr = 1.0e13; 
  };
}
 
#endif /* CSISTUDYMODULE_H */