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 ]

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