45{
46
47 auto analysisManager = G4AnalysisManager::Instance();
48 analysisManager->SetDefaultFileType("root");
49
50
51 analysisManager->SetVerboseLevel(1);
52
53
54 analysisManager->SetNtupleMerging(true);
55
56 analysisManager->SetFileName("B5");
57
58
59
60
61 analysisManager
62 ->CreateH1("Chamber1","Drift Chamber 1 # Hits", 50, 0., 50);
63 analysisManager
64 ->CreateH1("Chamber2","Drift Chamber 2 # Hits", 50, 0., 50);
65
66
67 analysisManager
68 ->CreateH2("Chamber1 XY","Drift Chamber 1 X vs Y",
69 50, -1000., 1000, 50, -300., 300.);
70 analysisManager
71 ->CreateH2("Chamber2 XY","Drift Chamber 2 X vs Y",
72 50, -1500., 1500, 50, -300., 300.);
73
74
76 analysisManager->CreateNtuple("B5", "Hits");
77 analysisManager->CreateNtupleIColumn("Dc1Hits");
78 analysisManager->CreateNtupleIColumn("Dc2Hits");
79 analysisManager->CreateNtupleDColumn("ECEnergy");
80 analysisManager->CreateNtupleDColumn("HCEnergy");
81 analysisManager->CreateNtupleDColumn("Time1");
82 analysisManager->CreateNtupleDColumn("Time2");
83 analysisManager
85 analysisManager
87 analysisManager->FinishNtuple();
88 }
89
90
91 analysisManager->SetNtupleFileName(0, "B5ntuple");
92}
std::vector< G4double > & GetHadCalEdep()
std::vector< G4double > & GetEmCalEdep()
EventAction * fEventAction