Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/advanced/composite_calorimeter/src/CCalRunAction.cc

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /examples/advanced/composite_calorimeter/src/CCalRunAction.cc (Version 11.3.0) and /examples/advanced/composite_calorimeter/src/CCalRunAction.cc (Version 11.0.p4)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //////////////////////////////////////////////     26 ///////////////////////////////////////////////////////////////////////////////
 27 // File: CCalRunAction.cc                          27 // File: CCalRunAction.cc
 28 // Description: A class for providing user act     28 // Description: A class for providing user actions at begin and end of run
 29 //////////////////////////////////////////////     29 ///////////////////////////////////////////////////////////////////////////////
 30 #include "CCalRunAction.hh"                        30 #include "CCalRunAction.hh"
 31                                                    31 
 32 #include "globals.hh"                              32 #include "globals.hh"
 33 #include "G4Run.hh"                                33 #include "G4Run.hh"
 34                                                    34 
 35 #include "G4ios.hh"                                35 #include "G4ios.hh"
 36                                                    36 
 37 #include "G4AnalysisManager.hh"                    37 #include "G4AnalysisManager.hh"
 38 #include "G4Threading.hh"                          38 #include "G4Threading.hh"
 39                                                    39 
 40                                                    40 
 41 CCalRunAction::CCalRunAction()                     41 CCalRunAction::CCalRunAction()
 42 {                                                  42 {
 43   numberOfTimeSlices = 200;                        43   numberOfTimeSlices = 200;
 44   Book();                                          44   Book();
 45 }                                                  45 }
 46                                                    46 
 47                                                    47 
 48 CCalRunAction::~CCalRunAction()                    48 CCalRunAction::~CCalRunAction()
 49 {                                                  49 {
 50 }                                                  50 }
 51                                                    51 
 52                                                    52 
 53 void CCalRunAction::BeginOfRunAction(const G4R     53 void CCalRunAction::BeginOfRunAction(const G4Run* aRun) 
 54 {                                                  54 {
 55   G4cout << "### Run " << aRun->GetRunID() <<      55   G4cout << "### Run " << aRun->GetRunID() << " start." << G4endl;
 56                                                    56 
 57   G4AnalysisManager* analysis = G4AnalysisMana     57   G4AnalysisManager* analysis = G4AnalysisManager::Instance();
 58   analysis->SetDefaultFileType("root");            58   analysis->SetDefaultFileType("root");
 59                                                    59 
 60   //cleanup                                        60   //cleanup
 61   G4int timeHist = analysis->GetH1Id("h300");      61   G4int timeHist = analysis->GetH1Id("h300");
 62   for (G4int i=0; i<numberOfTimeSlices; ++i) {     62   for (G4int i=0; i<numberOfTimeSlices; ++i) {
 63     analysis->GetH1(timeHist+i)->reset();          63     analysis->GetH1(timeHist+i)->reset();
 64   }                                                64   }
 65                                                    65 
 66   // Open an output file                           66   // Open an output file
 67   analysis->OpenFile("ccal");                      67   analysis->OpenFile("ccal");
 68   G4cout << "*********************************     68   G4cout << "********************************************" << G4endl
 69          << "* o/p file ccal"  << G4endl           69          << "* o/p file ccal"  << G4endl
 70          << "*********************************     70          << "********************************************" << G4endl 
 71          << G4endl;                                71          << G4endl;
 72 }                                                  72 }
 73                                                    73 
 74                                                    74 
 75 void CCalRunAction::EndOfRunAction(const G4Run     75 void CCalRunAction::EndOfRunAction(const G4Run* aRun) 
 76 {                                                  76 {
 77   G4cout << "### Run " << aRun->GetRunID() <<      77   G4cout << "### Run " << aRun->GetRunID() << " end." << G4endl;
 78                                                    78 
 79   // Close-out analysis: save histograms and n     79   // Close-out analysis: save histograms and ntuple
 80   G4AnalysisManager* analysisManager = G4Analy     80   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
 81   analysisManager->Write();                        81   analysisManager->Write();
 82   analysisManager->CloseFile();                    82   analysisManager->CloseFile();
 83 }                                                  83 }   
 84                                                    84 
 85                                                    85 
 86 void CCalRunAction::Book()                         86 void CCalRunAction::Book()
 87 {                                                  87 {
 88   G4AnalysisManager* analysisManager = G4Analy     88   G4AnalysisManager* analysisManager = G4AnalysisManager::Instance();
 89   analysisManager->SetVerboseLevel(1);             89   analysisManager->SetVerboseLevel(1);
 90   analysisManager->SetFirstHistoId(1);             90   analysisManager->SetFirstHistoId(1);
 91   analysisManager->SetFirstNtupleId(1);            91   analysisManager->SetFirstNtupleId(1);
 92                                                    92 
 93   // Note: merging ntuples is available only w     93   // Note: merging ntuples is available only with Root output
 94   if ( G4Threading::IsMultithreadedApplication     94   if ( G4Threading::IsMultithreadedApplication() ) analysisManager->SetNtupleMerging(true);
 95                                                    95  
 96   // Create ntuple                                 96   // Create ntuple
 97   analysisManager->CreateNtuple("ntuple1", "Ev     97   analysisManager->CreateNtuple("ntuple1", "Event info");
 98   for (G4int i=0;i<28;++i) {                       98   for (G4int i=0;i<28;++i) {
 99       G4String tupleidString = "hcal" + std::t     99       G4String tupleidString = "hcal" + std::to_string( i );
100       analysisManager->CreateNtupleFColumn( tu    100       analysisManager->CreateNtupleFColumn( tupleidString.c_str() );
101   }                                               101   }
102   for (G4int i=0; i<49; ++i) {                    102   for (G4int i=0; i<49; ++i) {
103       G4String tupleidString = "ecal" + std::t    103       G4String tupleidString = "ecal" + std::to_string( i );
104       analysisManager->CreateNtupleFColumn( tu    104       analysisManager->CreateNtupleFColumn( tupleidString.c_str() );
105   }                                               105   }
106   analysisManager->CreateNtupleFColumn("ELAB")    106   analysisManager->CreateNtupleFColumn("ELAB");
107   analysisManager->CreateNtupleFColumn("XPOS")    107   analysisManager->CreateNtupleFColumn("XPOS");
108   analysisManager->CreateNtupleFColumn("YPOS")    108   analysisManager->CreateNtupleFColumn("YPOS");
109   analysisManager->CreateNtupleFColumn("ZPOS")    109   analysisManager->CreateNtupleFColumn("ZPOS");
110   analysisManager->CreateNtupleFColumn("EDEP")    110   analysisManager->CreateNtupleFColumn("EDEP");
111   analysisManager->CreateNtupleFColumn("EDEC")    111   analysisManager->CreateNtupleFColumn("EDEC");
112   analysisManager->CreateNtupleFColumn("EHDC")    112   analysisManager->CreateNtupleFColumn("EHDC");
113   analysisManager->FinishNtuple();                113   analysisManager->FinishNtuple();
114                                                   114 
115   //                                              115   //
116   //Create histograms. Save the ID of the firs    116   //Create histograms. Save the ID of the first histogram in each block
117   //                                              117   //
118   //Energy deposit in Hcal layers                 118   //Energy deposit in Hcal layers 
119   for (G4int i = 0; i<28; ++i) {                  119   for (G4int i = 0; i<28; ++i) {
120     G4String idString = "h" + std::to_string(     120     G4String idString = "h" + std::to_string( i+100 );
121     G4String ntupletagString = "Energy Deposit    121     G4String ntupletagString = "Energy Deposit in Hcal Layer" + std::to_string( i ) + "  in GeV";
122     analysisManager->CreateH1( idString.c_str(    122     analysisManager->CreateH1( idString.c_str(), ntupletagString.c_str(), 100, 0., 1.0 );    
123   }                                               123   }
124   // Energy deposits in Ecal towers               124   // Energy deposits in Ecal towers
125   for (G4int i = 0; i<49; ++i) {                  125   for (G4int i = 0; i<49; ++i) {
126     G4String idString = "h" + std::to_string(     126     G4String idString = "h" + std::to_string( i+200 );
127     G4String ntupletagString = "Energy Deposit    127     G4String ntupletagString = "Energy Deposit in Ecal Tower" + std::to_string( i ) + "  in GeV";
128     analysisManager->CreateH1( idString.c_str(    128     analysisManager->CreateH1( idString.c_str(), ntupletagString.c_str(), 100, 0., 1.0 );      
129   }                                               129   }
130   // Total energy deposit                         130   // Total energy deposit
131   analysisManager->CreateH1( "h4000", "Total e    131   analysisManager->CreateH1( "h4000", "Total energy deposited in GeV", 100, 0., 100.0 );
132                                                   132 
133   // Time slices                                  133   // Time slices          
134   for (G4int i=0; i<numberOfTimeSlices; ++i){     134   for (G4int i=0; i<numberOfTimeSlices; ++i){
135     G4String idString = "h" + std::to_string(     135     G4String idString = "h" + std::to_string( i+300 );
136     G4String ntupletagString = "Time slice " +    136     G4String ntupletagString = "Time slice " + std::to_string( i ) + " nsec energy profile in GeV";
137     analysisManager->CreateH1( idString.c_str(    137     analysisManager->CreateH1( idString.c_str(), ntupletagString.c_str(), 100, 0., 100.0 );   
138   }                                               138   }
139                                                   139 
140   // Profile of lateral energy deposit in Hcal    140   // Profile of lateral energy deposit in Hcal
141   for (G4int i = 0; i<70; ++i) {                  141   for (G4int i = 0; i<70; ++i) {
142     G4String idString = "h" + std::to_string(     142     G4String idString = "h" + std::to_string( i+500 );
143     G4String ntupletagString = "Lateral energy    143     G4String ntupletagString = "Lateral energy profile at " + std::to_string( i ) + " cm  in GeV";
144     analysisManager->CreateH1( idString.c_str(    144     analysisManager->CreateH1( idString.c_str(), ntupletagString.c_str(), 100, 0., 10.0 );    
145   }                                               145   }
146                                                   146 
147   // Time profile                                 147   // Time profile 
148   analysisManager->CreateH1( "h901", "Time Pro    148   analysisManager->CreateH1( "h901", "Time Profile in Sensitive Detector", 200, 0., 200. );
149   analysisManager->CreateH1( "h902", "Time Pro    149   analysisManager->CreateH1( "h902", "Time Profile in Sensitive+Passive",  200, 0., 200. );
150                                                   150 
151   return;                                         151   return;
152 }                                                 152 }
153                                                   153