EventAction Class Reference

#include <EventAction.hh>

Inheritance diagram for EventAction:

[legend]

Collaboration diagram for EventAction:

[legend]

List of all members.

Public Member Functions
	EventAction (RunAction *)
	~EventAction ()
void	BeginOfEventAction (const G4Event *)
void	EndOfEventAction (const G4Event *)
void	AddEdep (G4double edep)
Private Attributes
RunAction *	m_runAction
G4double	m_edepPerEvent
G4int	m_printModulo

---

Detailed Description

Definition at line 16 of file EventAction.hh.

---

Constructor & Destructor Documentation

EventAction::EventAction

(

RunAction *

run

)

Definition at line 22 of file EventAction.cc.

:m_runAction(run), m_printModulo(5000)
{}

EventAction::~EventAction

(

)

Definition at line 28 of file EventAction.cc.

{}

---

Member Function Documentation

void EventAction::BeginOfEventAction

(

const G4Event *

evt

)

Definition at line 33 of file EventAction.cc.

{
 G4int evtNb = evt->GetEventID();

 //printing survey
 if (evtNb%m_printModulo == 0)
    G4cout << "\n---> Begin of Event: " << evtNb << G4endl;
    
 //energy deposited per event
 m_edepPerEvent = 0.;   
}

void EventAction::EndOfEventAction

(

const G4Event *

event

)

Definition at line 47 of file EventAction.cc.

{
  //energy deposited per event
  //
  if (m_edepPerEvent > 0.) {
   m_runAction->SumEvents(m_edepPerEvent);
  }

  // -- jour4b2:
  // -------------------------------
  // -- Get back collection pointer:
  // -------------------------------
  // -- find back collection, this is made by getting the collection ID:
  G4SDManager*             SDmanager = G4SDManager::GetSDMpointer();
  G4int        astronautCollectionID = SDmanager->GetCollectionID("AstronautCollection");
  AstroHitCollection*           astroCollection(0);
  // -- Collection ID is correct, we get the pointer of the collection:
  if (astronautCollectionID>=0) astroCollection = (AstroHitCollection*)(event->GetHCofThisEvent()->GetHC(astronautCollectionID));
  else G4cout << "Collection `AstronautCollection' not found !!!!" << G4endl;

  // -- Question : why do we need to cast (AstroHitCollection*) two lines above ?

  // --------------------------------------------
  // -- print hits found for the 5 first events:
  // --------------------------------------------
  if (event->GetEventID() < 5)
    {
      G4cout << "Event::EndOfEventAction(...) for event # " << event->GetEventID() << " : " << G4endl;
      // -- get back detector pointer to get astronaut mass:
      DetectorConstruction* detector = (DetectorConstruction*)(G4RunManager::GetRunManager()->GetUserDetectorConstruction());
      G4double mass = detector->GetAstronaut()->GetLogicalVolume()->GetMass();
      if (astroCollection)
        {
          // -- Now print:
          G4cout << " - number of hits : " << astroCollection->entries() << G4endl;
          if  (astroCollection->entries()) G4cout << " - hit list : " << G4endl;
          else                             G4cout << " - hit list : (empty)" << G4endl;
          for (G4int i=0; i<astroCollection->entries(); i++)
            {
              AstroHit* hit = (*astroCollection)[i];
              G4cout << "   k = " << hit->GetSlice() 
                     << " , dose (gray) = "  << hit->GetEnergy()/mass/gray 
                     << " , nb protons : "   << hit->GetNbProton()
                     << " , <eKin>(MeV) = "  << hit->GetProtonMeanEkin()/MeV << G4endl;
            }
        }
    }

  // ------------------------------------------
  // -- fill the RunAction table for SD vector:
  // ------------------------------------------
  if (astroCollection)
    {
      for (G4int i=0; i<astroCollection->entries(); i++)
        {
          AstroHit* hit = (*astroCollection)[i];
          m_runAction->SumDepthDoseSD(hit->GetSlice(),
                                      hit->GetEnergy(),
                                      hit->GetNbProton(),
                                      hit->GetProtonTotalEkin());
        }
    }


  // ----------------------------------------------
  // -- scorer
  // ----------------------------------------------
  // -- read back edep "collection" from scorer:
  G4int floorEdepCollectionID = SDmanager->GetCollectionID("floorScorer/edepFloor");
   
  G4THitsMap<G4double>* eventHitMap = (G4THitsMap<G4double>*)(event->GetHCofThisEvent()->GetHC(floorEdepCollectionID));
  std::map<G4int,G4double*>::iterator itr = eventHitMap->GetMap()->begin();
  G4cout << "Number of edep hits in floor : " << eventHitMap->entries() << G4endl;
  for(; itr != eventHitMap->GetMap()->end(); itr++)
    {
      G4int    key = (itr->first);
      G4double val = *(itr->second);
      G4cout << " key = " << key << " val (MeV) = " << val/MeV << G4endl;
    }



  // -- read back track counter "collection" from scorer:
  G4int floorTrackCounterCollectionID = SDmanager->GetCollectionID("floorScorer/trackCounterFloor");
  
  eventHitMap = (G4THitsMap<G4double>*)(event->GetHCofThisEvent()->GetHC(floorTrackCounterCollectionID));
  itr = eventHitMap->GetMap()->begin();
  G4cout << "Number of track hits in floor : " << eventHitMap->entries() << G4endl;
  for(; itr != eventHitMap->GetMap()->end(); itr++)
    {
      G4int    key = (itr->first);
      G4double val = *(itr->second);
      G4cout << " key = " << key << " , # tracks = " << val << G4endl;
    }

  
}

void EventAction::AddEdep

(

G4double

edep

)

[inline]

Definition at line 26 of file EventAction.hh.

Referenced by SteppingAction::UserSteppingAction().

{m_edepPerEvent += edep;};

---

Member Data Documentation

RunAction* EventAction::m_runAction [private]

Definition at line 26 of file EventAction.hh.

Referenced by EndOfEventAction().

G4double EventAction::m_edepPerEvent [private]

Definition at line 30 of file EventAction.hh.

Referenced by BeginOfEventAction(), and EndOfEventAction().

G4int EventAction::m_printModulo [private]

Definition at line 31 of file EventAction.hh.

Referenced by BeginOfEventAction().

---

The documentation for this class was generated from the following files:

• include/EventAction.hh
• src/EventAction.cc