B5EventAction Class Reference

Event action. More...

#include <B5EventAction.hh>

Inheritance diagram for B5EventAction:

Collaboration diagram for B5EventAction:

Public Member Functions
	B5EventAction ()
virtual	~B5EventAction ()
virtual void	BeginOfEventAction (const G4Event *)
virtual void	EndOfEventAction (const G4Event *)
std::vector< G4double > &	GetEmCalEdep ()
std::vector< G4double > &	GetHadCalEdep ()
Public Member Functions inherited from G4UserEventAction
	G4UserEventAction ()
virtual	~G4UserEventAction ()
virtual void	SetEventManager (G4EventManager *value)

Private Attributes
G4int	fHHC1ID
G4int	fHHC2ID
G4int	fDHC1ID
G4int	fDHC2ID
G4int	fECHCID
G4int	fHCHCID
std::vector< G4double >	fEmCalEdep
std::vector< G4double >	fHadCalEdep

Additional Inherited Members
Protected Attributes inherited from G4UserEventAction
G4EventManager *	fpEventManager

Detailed Description

Event action.

Definition at line 42 of file B5EventAction.hh.

Constructor & Destructor Documentation

B5EventAction()

B5EventAction::B5EventAction

(

)

Definition at line 50 of file B5EventAction.cc.

: G4UserEventAction(), 
  fHHC1ID(-1),
  fHHC2ID(-1),
  fDHC1ID(-1),
  fDHC2ID(-1),
  fECHCID(-1),
  fHCHCID(-1),
  fEmCalEdep(), 
  fHadCalEdep()
{
  // set printing per each event
  G4RunManager::GetRunManager()->SetPrintProgress(1);

  // initialize the vectors
  fEmCalEdep.resize(80, 0.);
  fHadCalEdep.resize(20, 0.);
}

Here is the call graph for this function:

~B5EventAction()

B5EventAction::~B5EventAction ( )

virtual

Definition at line 71 of file B5EventAction.cc.

{}

Member Function Documentation

BeginOfEventAction()

void B5EventAction::BeginOfEventAction ( const G4Event *  )

virtual

Reimplemented from G4UserEventAction.

Definition at line 76 of file B5EventAction.cc.

{
    if (fHHC1ID==-1) {
      G4SDManager* sdManager = G4SDManager::GetSDMpointer();
      fHHC1ID = sdManager->GetCollectionID("hodoscope1/hodoscopeColl");
      fHHC2ID = sdManager->GetCollectionID("hodoscope2/hodoscopeColl");
      fDHC1ID = sdManager->GetCollectionID("chamber1/driftChamberColl");
      fDHC2ID = sdManager->GetCollectionID("chamber2/driftChamberColl");
      fECHCID = sdManager->GetCollectionID("EMcalorimeter/EMcalorimeterColl");
      fHCHCID = sdManager->GetCollectionID("HadCalorimeter/HadCalorimeterColl");
    }
}     

Here is the call graph for this function:

EndOfEventAction()

void B5EventAction::EndOfEventAction ( const G4Event *  event )

virtual

Reimplemented from G4UserEventAction.

Definition at line 91 of file B5EventAction.cc.

{
    G4HCofThisEvent* hce = event->GetHCofThisEvent();
    if (!hce) 
    {
        G4ExceptionDescription msg;
        msg << "No hits collection of this event found." << G4endl; 
        G4Exception("B5EventAction::EndOfEventAction()",
                    "B5Code001", JustWarning, msg);
        return;
    }   


    // Get hits collections 
    B5HodoscopeHitsCollection* hHC1 
      = static_cast<B5HodoscopeHitsCollection*>(hce->GetHC(fHHC1ID));
      
    B5HodoscopeHitsCollection* hHC2 
      = static_cast<B5HodoscopeHitsCollection*>(hce->GetHC(fHHC2ID));
      
    B5DriftChamberHitsCollection* dHC1 
      = static_cast<B5DriftChamberHitsCollection*>(hce->GetHC(fDHC1ID));
      
    B5DriftChamberHitsCollection* dHC2 
      = static_cast<B5DriftChamberHitsCollection*>(hce->GetHC(fDHC2ID));
      
    B5EmCalorimeterHitsCollection* ecHC 
      = static_cast<B5EmCalorimeterHitsCollection*>(hce->GetHC(fECHCID));
      
    B5HadCalorimeterHitsCollection* hcHC 
      = static_cast<B5HadCalorimeterHitsCollection*>(hce->GetHC(fHCHCID));
      
    if ( (!hHC1) || (!hHC2) || (!dHC1) || (!dHC2) || (!ecHC) || (!hcHC) ) 
    {
        G4ExceptionDescription msg;
        msg << "Some of hits collections of this event not found." << G4endl; 
        G4Exception("B5EventAction::EndOfEventAction()",
                    "B5Code001", JustWarning, msg);
        return;
    }   
    
    //
    // Fill histograms & ntuple
    // 
    
    // Get analysis manager
    G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
 
    // Fill histograms
 
    G4int n_hit = dHC1->entries();
    analysisManager->FillH1(0, n_hit);

    for (G4int i=0;i<n_hit;i++)
    {
       B5DriftChamberHit* hit = (*dHC1)[i];
       G4ThreeVector localPos = hit->GetLocalPos();
       analysisManager->FillH2(0, localPos.x(), localPos.y());
    }
 
    n_hit = dHC2->entries();
    analysisManager->FillH1(1, n_hit);

    for (G4int i=0;i<n_hit;i++)
    {
       B5DriftChamberHit* hit = (*dHC2)[i];
       G4ThreeVector localPos = hit->GetLocalPos();
       analysisManager->FillH2(1, localPos.x(), localPos.y());
    }
    
        
    // Fill ntuple
    
    // Dc1Hits
    analysisManager->FillNtupleIColumn(0, dHC1->entries());
    // Dc2Hits
    analysisManager->FillNtupleIColumn(1, dHC1->entries());
    
    // ECEnergy
    G4int totalEmHit = 0;
    G4double totalEmE = 0.;
    for (G4int i=0;i<80;i++)
    {
        B5EmCalorimeterHit* hit = (*ecHC)[i];
        G4double eDep = hit->GetEdep();
        if (eDep>0.)
        {
            totalEmHit++;
            totalEmE += eDep;
        }
        fEmCalEdep[i] = eDep;
    }
    analysisManager->FillNtupleDColumn(2, totalEmE);

    // HCEnergy
    G4int totalHadHit = 0;
    G4double totalHadE = 0.;
    for (G4int i=0;i<20;i++)
    {
        B5HadCalorimeterHit* hit = (*hcHC)[i];
        G4double eDep = hit->GetEdep();
        if (eDep>0.)
        {
            totalHadHit++;
            totalHadE += eDep;
        }
        fHadCalEdep[i] = eDep;
    }
    analysisManager->FillNtupleDColumn(3, totalHadE);

    // Time 1
    for (G4int i=0;i<hHC1->entries();i++)
    {
      analysisManager->FillNtupleDColumn(4,(*hHC1)[i]->GetTime());
    }
      
    // Time 2
    for (G4int i=0;i<hHC2->entries();i++)
    {
      analysisManager->FillNtupleDColumn(5,(*hHC2)[i]->GetTime());
    }
    
    analysisManager->AddNtupleRow();  
    
    //
    // Print diagnostics
    // 
    
    G4int printModulo = G4RunManager::GetRunManager()->GetPrintProgress();
    if ( printModulo==0 || event->GetEventID() % printModulo != 0) return;
    
    G4PrimaryParticle* primary = event->GetPrimaryVertex(0)->GetPrimary(0);
    G4cout << G4endl
           << ">>> Event " << event->GetEventID() << " >>> Simulation truth : "
           << primary->GetG4code()->GetParticleName()
           << " " << primary->GetMomentum() << G4endl;
    
    // Hodoscope 1
    n_hit = hHC1->entries();
    G4cout << "Hodoscope 1 has " << n_hit << " hits." << G4endl;
    for (G4int i=0;i<n_hit;i++)
    {
        B5HodoscopeHit* hit = (*hHC1)[i];
        hit->Print();
    }

    // Hodoscope 2
    n_hit = hHC2->entries();
    G4cout << "Hodoscope 2 has " << n_hit << " hits." << G4endl;
    for (G4int i=0;i<n_hit;i++)
    {
        B5HodoscopeHit* hit = (*hHC2)[i];
        hit->Print();
    }

    // Drift chamber 1
    n_hit = dHC1->entries();
    G4cout << "Drift Chamber 1 has " << n_hit << " hits." << G4endl;
    for (G4int i2=0;i2<5;i2++)
    {
        for (G4int i=0;i<n_hit;i++)
        {
            B5DriftChamberHit* hit = (*dHC1)[i];
            if (hit->GetLayerID()==i2) hit->Print();
        }
    }

    // Drift chamber 2
    n_hit = dHC2->entries();
    G4cout << "Drift Chamber 2 has " << n_hit << " hits." << G4endl;
    for (G4int i2=0;i2<5;i2++)
    {
        for (G4int i=0;i<n_hit;i++)
        {
            B5DriftChamberHit* hit = (*dHC2)[i];
            if (hit->GetLayerID()==i2) hit->Print();
        }
    }

    // EM calorimeter
    G4cout << "EM Calorimeter has " << totalEmHit << " hits. Total Edep is "
    << totalEmE/MeV << " (MeV)" << G4endl;

    // Had calorimeter
    G4cout << "Hadron Calorimeter has " << totalHadHit << " hits. Total Edep is "
    << totalHadE/MeV << " (MeV)" << G4endl;
}

Here is the call graph for this function:

GetEmCalEdep()

std::vector<G4double>& B5EventAction::GetEmCalEdep ( )

inline

Definition at line 51 of file B5EventAction.hh.

{ return fEmCalEdep; }

Here is the caller graph for this function:

GetHadCalEdep()

std::vector<G4double>& B5EventAction::GetHadCalEdep ( )

inline

Definition at line 52 of file B5EventAction.hh.

{ return fEmCalEdep; }

Here is the caller graph for this function:

Member Data Documentation

fDHC1ID

G4int B5EventAction::fDHC1ID

private

Definition at line 57 of file B5EventAction.hh.

fDHC2ID

G4int B5EventAction::fDHC2ID

private

Definition at line 58 of file B5EventAction.hh.

fECHCID

G4int B5EventAction::fECHCID

private

Definition at line 59 of file B5EventAction.hh.

fEmCalEdep

std::vector<G4double> B5EventAction::fEmCalEdep

private

Definition at line 61 of file B5EventAction.hh.

fHadCalEdep

std::vector<G4double> B5EventAction::fHadCalEdep

private

Definition at line 62 of file B5EventAction.hh.

fHCHCID

G4int B5EventAction::fHCHCID

private

Definition at line 60 of file B5EventAction.hh.

fHHC1ID

G4int B5EventAction::fHHC1ID

private

Definition at line 55 of file B5EventAction.hh.

fHHC2ID

G4int B5EventAction::fHHC2ID

private

Definition at line 56 of file B5EventAction.hh.

---

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

• B5EventAction.hh
• B5EventAction.cc