event.hh

Go to the documentation of this file.

 
#ifndef EVENT_HH
#define EVENT_HH
 
#include <vector>
 
#include "G4RunManager.hh"
#include "G4UserEventAction.hh"
#include "G4Event.hh"
#include "G4AnalysisManager.hh"
 
#include "globalsettings.hh"
#include "hit.hh"
#include "generator.hh"
 
class MyEventAction : public G4UserEventAction
{
public:
    MyEventAction() = default; 
    ~MyEventAction() override = default; 
    void BeginOfEventAction(const G4Event *event) override;
    void EndOfEventAction(const G4Event *event) override;
 
    inline void SetArrivalandFirstInteraction(G4double newtimein, G4ThreeVector newposin, G4double newtimeinter, G4ThreeVector newposinter)
    {
        if(newtimein < fTimeIn)
        {
            fTimeIn = newtimein;
            fPosXIn = newposin.x();
            fPosYIn = newposin.y();
            fPosZIn = newposin.z();
 
            fTimeFirstInter = newtimeinter;
            fPosXFirstInter = newposinter.x();
            fPosYFirstInter = newposinter.y();
            fPosZFirstInter = newposinter.z();
        }
    }
    inline void AddEdep(G4double edep) { fEdep += edep; }
    inline void SetMaxEdep(G4double edepondx, G4ThreeVector maxedeppos)
    {
        if(edepondx > fMaxEdep)
        {
            fMaxEdep = edepondx;
            fMaxEdepPos = maxedeppos;
        }
    };
 
 
    inline void SetDecayTriggerSi(G4bool trg) { fDecayTriggerSi = trg; }
    inline void SetCosmicTriggerUp(G4bool trg) { fCosmicTriggerUp = trg; }
    inline void SetCosmicTriggerBottom(G4bool trg) { fCosmicTriggerBottom = trg; }
 
 
    // Primary's data
    G4double fTimeIn, 
             fPosXIn, 
             fPosYIn, 
             fPosZIn, 
             fTimeFirstInter, 
             fPosXFirstInter, 
             fPosYFirstInter, 
             fPosZFirstInter; 
    // Crystal's data
    G4double      fEdep; 
    G4double      fMaxEdep; 
    G4ThreeVector fMaxEdepPos; 
    // Detectors' data
    G4int                 fHitsNum_F, 
                          fHitsNum_B; 
    std::vector<G4double> fT_F, 
                          fX_F, 
                          fY_F, 
                          fT_B, 
                          fX_B, 
                          fY_B; 
    std::vector<G4int>    fHitsNum_F_Ch,
                          fChannel_F, 
                          fHitsNum_B_Ch,
                          fChannel_B; 
    // Intercalibration triggers
    G4double fTimeOfDecay;
    G4bool   fDecayTriggerSi,
             fCosmicTriggerUp,
             fCosmicTriggerBottom;
};
 
#endif  // EVENT_HH