ChargeExchangeMC/src/CexmcSimpleEnergyDeposit.cc

0001 //
0002 // ********************************************************************
0003 // * License and Disclaimer                                           *
0004 // *                                                                  *
0005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
0006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
0007 // * conditions of the Geant4 Software License,  included in the file *
0008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
0009 // * include a list of copyright holders.                             *
0010 // *                                                                  *
0011 // * Neither the authors of this software system, nor their employing *
0012 // * institutes,nor the agencies providing financial support for this *
0013 // * work  make  any representation or  warranty, express or implied, *
0014 // * regarding  this  software system or assume any liability for its *
0015 // * use.  Please see the license in the file  LICENSE  and URL above *
0016 // * for the full disclaimer and the limitation of liability.         *
0017 // *                                                                  *
0018 // * This  code  implementation is the result of  the  scientific and *
0019 // * technical work of the GEANT4 collaboration.                      *
0020 // * By using,  copying,  modifying or  distributing the software (or *
0021 // * any work based  on the software)  you  agree  to acknowledge its *
0022 // * use  in  resulting  scientific  publications,  and indicate your *
0023 // * acceptance of all terms of the Geant4 Software license.          *
0024 // ********************************************************************
0025 //
0026 /*
0027  * ============================================================================
0028  *
0029  *       Filename:  CexmcSimpleEnergyDeposit.cc
0030  *
0031  *    Description:  simple energy deposit scorer
0032  *
0033  *        Version:  1.0
0034  *        Created:  14.11.2009 12:48:22
0035  *       Revision:  none
0036  *       Compiler:  gcc
0037  *
0038  *         Author:  Alexey Radkov (),
0039  *        Company:  PNPI
0040  *
0041  * ============================================================================
0042  */
0043
0044 #include <G4UnitsTable.hh>
0045 #include <G4HCofThisEvent.hh>
0046 #include "CexmcSimpleEnergyDeposit.hh"
0047
0048
0049 CexmcSimpleEnergyDeposit::CexmcSimpleEnergyDeposit( const G4String &  name ) :
0050     CexmcPrimitiveScorer( name ), hcId( -1 )
0051 {
0052 }
0053
0054
0055 G4int  CexmcSimpleEnergyDeposit::GetIndex( G4Step * )
0056 {
0057     return 0;
0058 }
0059
0060
0061 G4bool  CexmcSimpleEnergyDeposit::ProcessHits( G4Step *  step,
0062                                                G4TouchableHistory * )
0063 {
0064     G4double  energyDeposit( step->GetTotalEnergyDeposit() );
0065
0066     if ( energyDeposit == 0. )
0067         return false;
0068
0069     eventMap->add( GetIndex( step ), energyDeposit );
0070
0071     return true;
0072 }
0073
0074
0075 void  CexmcSimpleEnergyDeposit::Initialize( G4HCofThisEvent *  hcOfEvent )
0076 {
0077     eventMap = new CexmcEnergyDepositCollection( detector->GetName(),
0078                                                  primitiveName );
0079
0080     if ( hcId < 0 )
0081         hcId = GetCollectionID( 0 );
0082
0083     hcOfEvent->AddHitsCollection( hcId, eventMap );
0084 }
0085
0086
0087 void  CexmcSimpleEnergyDeposit::EndOfEvent( G4HCofThisEvent * )
0088 {
0089     if ( GetVerboseLevel() > 0 )
0090         PrintAll();
0091 }
0092
0093
0094 void  CexmcSimpleEnergyDeposit::clear( void )
0095 {
0096     eventMap->clear();
0097 }
0098
0099
0100 void  CexmcSimpleEnergyDeposit::DrawAll( void )
0101 {
0102 }
0103
0104
0105 void  CexmcSimpleEnergyDeposit::PrintAll( void )
0106 {
0107     G4int   nmbOfEntries( eventMap->entries() );
0108
0109     if ( nmbOfEntries == 0 )
0110         return;
0111
0112     PrintHeader( nmbOfEntries );
0113
0114     /* index is always 0 */
0115     for ( CexmcEnergyDepositCollectionData::iterator
0116                          itr( eventMap->GetMap()->begin() );
0117                                      itr != eventMap->GetMap()->end(); ++itr )
0118     {
0119         G4cout << "       energy deposit: " <<
0120                 G4BestUnit( *( itr->second ), "Energy" ) << G4endl;
0121     }
0122 }
0123