| Message: Re: G4VPrimitive Scorer-Surface Current | Not Logged In (login) |
|
Hi.
This is based on my previous question My aim is to get the number of photons entering the scorer. I used G4PSPassageCurrent and G4PSFlatSurfaceCurrent in my code:
G4SDManager* manager = G4SDManager::GetSDMpointer(); G4String MFDname = "MyDetector"; G4MultiFunctionalDetector* MFD = new G4MultiFunctionalDetector(MFDname);
manager->AddNewDetector(MFD); logicHoles->SetSensitiveDetector(MFD);
//filter G4String fltName, particleName;
G4SDParticleFilter* gammafilter= new G4SDParticleFilter(fltName="gammaFilter", particleName="gamma");
G4PSPassageCurrent* scorer1 = new G4PSPassageCurrent("Current", fCurrent_In);
G4PSFlatSurfaceCurrent* scorer2 = new G4PSFlatSurfaceCurrent("FlatCurrent", fCurrent_In);
scorer1->SetFilter(gammafilter); scorer2->SetFilter(gammafilter);
MFD->RegisterPrimitive(scorer1); MFD->RegisterPrimitive(scorer2);
I run for 10000000, and here is the result:
---------------------------------------------
Number of Events Processed: 10000000
----------------------------------------------
# Cell Current FlatCurrent
0 21 7
1 22 9
2 34 7
3 23 10
4 23 11
5 15 8
6 16 8
7 22 10
8 26 16
9 27 13
10 20 15
Does this mean that there are 21 photons/mm3 for scorer1 and 7/mm2 photons in scorer 2 for #Cell 0, and so on? Inputs and suggestion are highly appreciated. Many thanks, Suhaili.
|
| Inline Depth: | Outline Depth: | Add message: |
|
to: |