CgemSimCheck.h

Go to the documentation of this file.

/*
 * =====================================================================================
 *       Filename:  CgemSimCheck.hh
 *    Description:  Read .rtraw and save it to a .txt
 *        Version:  1.0
 *        Created:  2014年02月10日 14时25分12秒
 *       Revision:  none
 *       Compiler:  gcc
 *         Author:  [email protected] 
 *   Organization:  
 * =====================================================================================
 */
 
#ifndef CGEMSIMCHECK
#define CGEMSIMCHECK
 
/* Header file: BOSS */
#include "GaudiKernel/Algorithm.h"
#include "RawDataProviderSvc/RawDataProviderSvc.h"
 
#include "MdcGeomSvc/IMdcGeomSvc.h"
#include "MdcGeomSvc/MdcGeomSvc.h"
#include "MdcGeomSvc/MdcGeoWire.h"
#include "MdcGeomSvc/MdcGeoLayer.h"
 
#include "CgemGeomSvc/CgemGeomSvc.h"
 
#include "GaudiKernel/NTuple.h"
#include "GaudiKernel/INTupleSvc.h"
 
/* Header file: C++ */
#include <iostream>
#include <fstream>
#include <string>
#include <vector>
 
/* Header file: Geant4 */
#include "G4VHit.hh"
#include "G4THitsCollection.hh"
#include "G4Allocator.hh"
#include "G4ThreeVector.hh"
class CgemSimCheck : public Algorithm
{
    public:
        CgemSimCheck(const std::string& name, ISvcLocator* pSvcLocator);
        ~CgemSimCheck();
 
        StatusCode initialize();
        StatusCode execute();
        StatusCode finalize();
        StatusCode beginRun();
 
    private:
        int     m_F_McTruth;
        fstream m_output_McTruth;
 
        int     m_F_Digi;
        fstream m_output_Digi;
 
        int     m_F_McParticle;
        fstream m_output_McParticle;
 
        string m_postfix;
 
        /* Digi */
        NTuple::Tuple* m_tuple;
        NTuple::Item<int> m_run; 
        NTuple::Item<int> m_evt; 
        NTuple::Item<int> m_nHit; 
        NTuple::Item<int> m_Xnstrip; 
        NTuple::Item<int> m_Vnstrip; 
        NTuple::Array<G4int>    m_ID_track;
        NTuple::Array<G4int>    m_X_ID_strip;
        NTuple::Array<G4int>    m_V_ID_strip;
        NTuple::Array<G4int>    m_X_nhit;
        NTuple::Array<G4int>    m_V_nhit;
        NTuple::Array<G4int>    m_ID_layer;
        NTuple::Array<G4int>    m_ID_sheet;
        NTuple::Array<G4double> m_ID_strip;
        NTuple::Array<G4double> m_globle_time;
        NTuple::Array<G4double> m_E_deposit;
        NTuple::Array<G4double> m_E_X;
        NTuple::Array<G4double> m_E_V;
 
        /* MC */
        NTuple::Tuple* mc_tuple;
        NTuple::Item<int> m_mc_nHit; 
        NTuple::Array<int>         m_mc_ID_track;    /* ID of current track */
        NTuple::Array<int>         m_mc_ID_layer;    /* ID of CgemLayer Hits belong to */
        NTuple::Item<G4int>         m_mc_pdg_code;    /* PDG code of the track particle */
        NTuple::Array<double>      m_mc_global_time; /* Time since the current event began */
        NTuple::Array<double>      m_mc_E_deposit;   /* Total energy deposited during the step */
        NTuple::Array<G4double>      m_mc_L_step;      /* Step length */
        NTuple::Array<double> m_mc_XYZ_pre_X;     /* Position of Pre  Point */
        NTuple::Array<double> m_mc_XYZ_pre_Y;     /* Position of Pre  Point */
        NTuple::Array<double> m_mc_XYZ_pre_Z;     /* Position of Pre  Point */
        NTuple::Array<double> m_mc_XYZ_post_X;    /* Position of Post Point */
        NTuple::Array<double> m_mc_XYZ_post_Y;    /* Position of Post Point */
        NTuple::Array<double> m_mc_XYZ_post_Z;    /* Position of Post Point */
        NTuple::Array<double> m_mc_P_pre_X;       /* Momentum of Pre  Point */
        NTuple::Array<double> m_mc_P_pre_Y;       /* Momentum of Pre  Point */
        NTuple::Array<double> m_mc_P_pre_Z;       /* Momentum of Pre  Point */
        NTuple::Array<double> m_mc_P_post_X;       /* Momentum of Post  Point */
        NTuple::Array<double> m_mc_P_post_Y;       /* Momentum of Post  Point */
        NTuple::Array<double> m_mc_P_post_Z;       /* Momentum of Post  Point */
/*mcparticle*/
        NTuple::Tuple* mcP_tuple;
        NTuple::Item<int>         Nparticle; 
        NTuple::Array<int>        m_trkindex;    
        NTuple::Array<double>     m_mcParticle_x;  
        NTuple::Array<double>     m_mcParticle_y;
        NTuple::Array<double>     m_mcParticle_z;
        NTuple::Array<double>     m_mcParticle_px;
        NTuple::Array<double>     m_mcParticle_py;
        NTuple::Array<double>     m_mcParticle_pz;
        NTuple::Array<double>     m_mcParticle_E;
        NTuple::Array<double>     m_mcParticle_phi;
        NTuple::Array<double>     m_mcParticle_costheta;
        NTuple::Array<double>     m_mcParticle_theta;
        NTuple::Array<double>     m_mcParticle_pt;
};
 
 
#endif