107 const G4Event*
event = G4RunManager::GetRunManager()->GetCurrentEvent();
109 evName += G4UIcommand::ConvertToString( event->GetEventID() );
110 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
113 analysisManager->CreateNtupleIColumn(
"particleID" );
114 analysisManager->CreateNtupleIColumn(
"PID" );
115 analysisManager->CreateNtupleDColumn(
"MC_pX" );
116 analysisManager->CreateNtupleDColumn(
"MC_pY" );
117 analysisManager->CreateNtupleDColumn(
"MC_pZ" );
119 analysisManager->CreateNtupleDColumn(
"tracker_res" );
120 analysisManager->CreateNtupleDColumn(
"tracker_eff" );
121 analysisManager->CreateNtupleDColumn(
"tracker_pX" );
122 analysisManager->CreateNtupleDColumn(
"tracker_pY" );
123 analysisManager->CreateNtupleDColumn(
"tracker_pZ" );
125 analysisManager->CreateNtupleDColumn(
"emcal_res" );
126 analysisManager->CreateNtupleDColumn(
"emcal_eff" );
127 analysisManager->CreateNtupleDColumn(
"emcal_X" );
128 analysisManager->CreateNtupleDColumn(
"emcal_Y" );
129 analysisManager->CreateNtupleDColumn(
"emcal_Z" );
130 analysisManager->CreateNtupleDColumn(
"emcal_E" );
132 analysisManager->CreateNtupleDColumn(
"hcal_res" );
133 analysisManager->CreateNtupleDColumn(
"hcal_eff" );
134 analysisManager->CreateNtupleDColumn(
"hcal_X" );
135 analysisManager->CreateNtupleDColumn(
"hcal_Y" );
136 analysisManager->CreateNtupleDColumn(
"hcal_Z" );
137 analysisManager->CreateNtupleDColumn(
"hcal_E" );
146 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
147 analysisManager->CreateH1(
"Pdiff",
"momentum smeared in tracker", 100, 0.8, 1.2 );
148 analysisManager->SetH1XAxisTitle( 0,
"p_{smeared}/p_{true}" );
149 analysisManager->SetH1YAxisTitle( 0,
"Entries" );
150 analysisManager->CreateH1(
"EMCalEdiff",
"energy smeared in EMCal", 100, 0.8, 1.2 );
151 analysisManager->SetH1XAxisTitle( 1,
"E_{smeared}/E_{true}" );
152 analysisManager->SetH1YAxisTitle( 1,
"Entries" );
153 analysisManager->CreateH1(
"HCalEdiff",
"energy smeared in HCal", 100, 0.0, 2.0 );
154 analysisManager->SetH1XAxisTitle( 2,
"E_{smeared}/E_{true}" );
155 analysisManager->SetH1YAxisTitle( 2,
"Entries" );
161 G4ThreeVector aVector, G4double aResolution,
162 G4double aEfficiency, G4double aEnergy ) {
163 G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
164 switch( aWhatToSave ) {
178 " Wrong particle - trying to save Tracker information of different particle"
189 " Wrong particle - trying to save EMCal information of different particle"
201 " Wrong particle - trying to save HCal information of different particle"
void SaveTrack(SaveType aWhatToSave, G4int aPartID, G4int aPDG, G4ThreeVector aVector, G4double aResolution=0, G4double aEfficiency=1, G4double aEnergy=0)
Saves the information about the particle (track).