47{
48
49 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
50 analysisManager->SetDefaultFileType("root");
52 analysisManager->SetVerboseLevel(1);
53 analysisManager->SetActivation(true);
54
55
56
57
58 G4int nbins = 100;
59 G4double vmin = 0.;
60 G4double vmax = 100.;
61
62
63 analysisManager->CreateH1("0", "dummy", nbins, vmin, vmax);
64
65 analysisManager->CreateH1("hits per event", "hits per event", nbins, vmin,
66 vmax);
67
68 analysisManager->CreateH1("hits above threshold",
69 "hits per event above threshold", nbins, vmin,
70 vmax);
71
72 analysisManager->CreateH1("scintillation", "scintillation photons per event",
73 nbins, vmin, vmax);
74
75 analysisManager->CreateH1("Cerenkov", "Cerenkov photons per event", nbins,
76 vmin, vmax);
77
78 analysisManager->CreateH1("absorbed", "absorbed photons per event", nbins,
79 vmin, vmax);
80
81 analysisManager->CreateH1("boundary absorbed",
82 "photons absorbed at boundary per event", nbins,
83 vmin, vmax);
84
85 analysisManager->CreateH1(
86 "E dep", "energy deposition in scintillator per event", nbins, vmin, vmax);
87
88
89 for(G4int i = 0; i < analysisManager->GetNofH1s(); ++i)
90 {
91 analysisManager->SetH1Activation(i, false);
92 }
93}