Geant4 Cross Reference |
1 // 1 2 // ******************************************* 3 // * License and Disclaimer 4 // * 5 // * The Geant4 software is copyright of th 6 // * the Geant4 Collaboration. It is provided 7 // * conditions of the Geant4 Software License 8 // * LICENSE and available at http://cern.ch/ 9 // * include a list of copyright holders. 10 // * 11 // * Neither the authors of this software syst 12 // * institutes,nor the agencies providing fin 13 // * work make any representation or warran 14 // * regarding this software system or assum 15 // * use. Please see the license in the file 16 // * for the full disclaimer and the limitatio 17 // * 18 // * This code implementation is the result 19 // * technical work of the GEANT4 collaboratio 20 // * By using, copying, modifying or distri 21 // * any work based on the software) you ag 22 // * use in resulting scientific publicati 23 // * acceptance of all terms of the Geant4 Sof 24 // ******************************************* 25 // 26 /// \file electromagnetic/TestEm12/src/HistoMa 27 /// \brief Implementation of the HistoManager 28 // 29 // 30 // 31 //....oooOO0OOooo........oooOO0OOooo........oo 32 //....oooOO0OOooo........oooOO0OOooo........oo 33 34 #include "HistoManager.hh" 35 36 #include "G4UnitsTable.hh" 37 38 //....oooOO0OOooo........oooOO0OOooo........oo 39 40 HistoManager::HistoManager() 41 { 42 Book(); 43 } 44 45 //....oooOO0OOooo........oooOO0OOooo........oo 46 47 void HistoManager::Book() 48 { 49 // Create or get analysis manager 50 // The choice of analysis technology is done 51 // in HistoManager.hh 52 G4AnalysisManager* analysisManager = G4Analy 53 analysisManager->SetDefaultFileType("root"); 54 analysisManager->SetFileName(fFileName); 55 analysisManager->SetVerboseLevel(1); 56 analysisManager->SetActivation(true); 57 // enable inactivation of histograms 58 59 // Define histograms start values 60 const G4int kMaxHisto = 9; 61 const G4String id[] = {"0", "1", "2", "3", " 62 63 const G4String title[] = { 64 "dummy", // 0 65 "Edep (MeV/mm) along absorber", // 1 66 "total Energy deposited in absorber", // 67 "true track length of the primary particle 68 "true step size of the primary particle", 69 "projected range of the primary particle", 70 "true track length of charged secondaries" 71 "true step size of charged secondaries", 72 "d(E/E0)/d(r/r0) along r/r0" // 8 73 }; 74 75 // Default values (to be reset via /analysis 76 G4int nbins = 100; 77 G4double vmin = 0.; 78 G4double vmax = 100.; 79 80 // Create all histograms as inactivated 81 // as we have not yet set nbins, vmin, vmax 82 for (G4int k = 0; k < kMaxHisto; k++) { 83 G4int ih = analysisManager->CreateH1(id[k] 84 analysisManager->SetH1Activation(ih, false 85 } 86 } 87