Ćwiczenie 1
a)
Proszę zdefiniować własności optyczne PbWO4 jako przykład
implementacji proszę wykorzystać kod ze strony GEANT'a
GEANT4
User's
Guide
b) skompilować i uruchomić program
sprawdzając, czy własności optyczne PBWO4 zostały poprawnie
zaimplementowane
Wskazówki:
- Własności optyczne materiału zapisane są jako elementy
tablicy G4MaterialPropertiesTable
przypisanej
do obiektu klasy G4Material
za pomocą funkcji SetMaterialPropertiesTable().
Wiersze tablicy G4MaterialPropertiesTable
mają postać pary [słowo-kluczowe
;
wartości] i dodajemy je do tablicy za pomocą funkcji AddConstProperty() gdy wartość jest
jedna lub AddProperty gdy wartość zmienia się wraz z energią
fotonu. Własności scyntylacyjne określamy za
pomocą słów kluczowych:
"SCINTILLATIONYIELD"
średnia
ilość emitowanych fotonów na jednostkę deponowanej
energii w scyntylatorze
"RESOLUTIONSCALE"
współczynnik zwiększający dyspersje
rozkładu liczby fotonów sigma=RES*sqrt(N)
"FASTCOMPONENT"
rozkład widmowy
"szybkich" fotonów
"FASTTIMECONSTANT"
stała czasowa eksponencjalnego rozkładu czasowego
emisji fotonów szybkich
"SLOWCOMPONENT"
i
"SLOWTIMECONSTANT" dla
składowej powolnej.
"YIELDRATIO"
względny udział fotonów szybkich
- Aby sprawdzić jakie wartości są zapisane w tablicy MaterialPropertiesTable można
użyć
funkcji DumpTable()
Ćwiczenie 2
Proszę
aktywować
proces produkcji fotonów optycznych w PbWO4 dla wszystkich
cząstek których ten proces dotyczy:
Wskazówki:
-
Proces ten reprezentowany jest przez klasę G4Scintillation
- Aktywacja tego procesu może wyglądać identycznie jak
aktywacja procesu rozpadu G4Decay w funkcji MyPhysicsList::ConstructGeneral omawianej na
wykładzie. Aby zminimalizować ilość zmian w kodzie można ten
proces aktywować w funkcji ConstructGeneral.
- Proces G4Scintillation ma podobnie jak
proces G4Decay komponentę AtRest oraz PostStep
- Uwaga: 1 GeV
elektron może wyprodukować niewiele mniej niż 200
tysięcy fotonów optycznych w scyntylatorze. Ich
wizualizacja jest prawie niemożliwa. Dlatego przed
uruchomieniem programu proszę zmniejszyć ilość
produkowanych fotonów w scyntylatorze do np.
1/MeV.
- Proszę sprawdzić czy program się kompiluje i czy podczas
wizualizacji widać wszystkie obiekty jak na rysunku
poniżej: