GEANT4学习第三节 B2中的SD 敏感探测器 2 hit信息输出

本文链接：https://blog.csdn.net/weixin_42449640/article/details/118195883

GEANT4学习第三节 B2中的SD 敏感探测器 2 Hit信息输出

之前的博文看了B2a这个例子的设置的类，接下来看输出信息的类。

B2a的输出信息类只有两个：

B2RunAction
B2EventAction

B2RunAction里面没什么可看的，主要还是看一下B2EventAction

1.B2EventAction

在B2EventAction的EndOfEventAction里面，例子实质上已经写了一个输出hit总数的例子，如下：

G4VHitsCollection* hc = event->GetHCofThisEvent()->GetHC(0);
G4cout << "    "  
       << hc->GetSize() << " hits stored in this event" << G4endl;

这里的

hc

就是我们之前写的hitCollection.

2.输出每个事件的沉积能量

单个事件沉积能量的输出，我们需要返回B2TrackerSD::EndOfEvent函数中去进行。

我们将这个函数改成如下：

void B2TrackerSD::EndOfEvent(G4HCofThisEvent*)
{
  //if ( verboseLevel>1 ) { 
     G4int nofHits = fHitsCollection->entries();

     G4cout << G4endl
            << "-------->Hits Collection: in this event they are " << nofHits 
            << " hits in the tracker chambers: " << G4endl;
     //for ( G4int i=0; i<nofHits; i++ ) (*fHitsCollection)[i]->Print();
     G4double edepofThisEvent=0.;

     for ( G4int i=0; i<nofHits; i++ ) {
     edepofThisEvent+=(*fHitsCollection)[i]->GetEdep();
     //G4cout <<"SAY: "<<(*fHitsCollection)[i]->GetTrackID()<< G4endl;  
     };
     G4cout <<edepofThisEvent<<" Mev energy deposit in this event"<< G4endl ;
  //}
}

这样就可以输出了。

### Run 0 starts.
--> Event 0 starts.

-------->Hits Collection: in this event they are 124 hits in the tracker chambers: 
0.751416 Mev energy deposit in this event
>>> Event: 0
    125 trajectories stored in this event.
    124 hits stored in this event

-------->Hits Collection: in this event they are 100 hits in the tracker chambers: 
0.640863 Mev energy deposit in this event
>>> Event: 1
    102 trajectories stored in this event.
    100 hits stored in this event

-------->Hits Collection: in this event they are 67 hits in the tracker chambers: 
0.444898 Mev energy deposit in this event
>>> Event: 2
    88 trajectories stored in this event.
    67 hits stored in this event

-------->Hits Collection: in this event they are 101 hits in the tracker chambers: 
0.629 Mev energy deposit in this event
>>> Event: 3
    104 trajectories stored in this event.
    101 hits stored in this event

-------->Hits Collection: in this event they are 0 hits in the tracker chambers: 
0 Mev energy deposit in this event
>>> Event: 4
    1066 trajectories stored in this event.
    0 hits stored in this event

-------->Hits Collection: in this event they are 0 hits in the tracker chambers: 
0 Mev energy deposit in this event
>>> Event: 5
    637 trajectories stored in this event.
    0 hits stored in this event

-------->Hits Collection: in this event they are 158 hits in the tracker chambers: 
0.967425 Mev energy deposit in this event
>>> Event: 6
    143 trajectories stored in this event.
    158 hits stored in this event

-------->Hits Collection: in this event they are 104 hits in the tracker chambers: 
0.699531 Mev energy deposit in this event
>>> Event: 7
    107 trajectories stored in this event.
    104 hits stored in this event

-------->Hits Collection: in this event they are 104 hits in the tracker chambers: 
0.595938 Mev energy deposit in this event
>>> Event: 8
    128 trajectories stored in this event.
    104 hits stored in this event

-------->Hits Collection: in this event they are 128 hits in the tracker chambers: 
0.868534 Mev energy deposit in this event
>>> Event: 9
    109 trajectories stored in this event.
    128 hits stored in this event