Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/runAndEvent/RE01/src/RE01CalorimeterSD.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/extended/runAndEvent/RE01/src/RE01CalorimeterSD.cc (Version 11.3.0) and /examples/extended/runAndEvent/RE01/src/RE01CalorimeterSD.cc (Version 10.2.p3)


  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 /// \file runAndEvent/RE01/src/RE01Calorimeter     26 /// \file runAndEvent/RE01/src/RE01CalorimeterSD.cc
 27 /// \brief Implementation of the RE01Calorimet     27 /// \brief Implementation of the RE01CalorimeterSD class
 28 //                                                 28 //
                                                   >>  29 // $Id: RE01CalorimeterSD.cc 68761 2013-04-05 12:35:00Z gcosmo $
 29 //                                                 30 //
 30                                                    31 
 31 #include "RE01CalorimeterSD.hh"                    32 #include "RE01CalorimeterSD.hh"
 32                                                << 
 33 #include "RE01CalorimeterHit.hh"                   33 #include "RE01CalorimeterHit.hh"
 34                                                <<  34 #include "G4VPhysicalVolume.hh"
 35 #include "G4LogicalVolume.hh"                      35 #include "G4LogicalVolume.hh"
 36 #include "G4ParticleDefinition.hh"             << 
 37 #include "G4Step.hh"                           << 
 38 #include "G4TouchableHistory.hh"               << 
 39 #include "G4Track.hh"                              36 #include "G4Track.hh"
 40 #include "G4VPhysicalVolume.hh"                <<  37 #include "G4Step.hh"
                                                   >>  38 #include "G4ParticleDefinition.hh"
 41 #include "G4VTouchable.hh"                         39 #include "G4VTouchable.hh"
                                                   >>  40 #include "G4TouchableHistory.hh"
 42 #include "G4ios.hh"                                41 #include "G4ios.hh"
 43                                                    42 
 44 //....oooOO0OOooo........oooOO0OOooo........oo     43 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 45 RE01CalorimeterSD::RE01CalorimeterSD(G4String      44 RE01CalorimeterSD::RE01CalorimeterSD(G4String name)
 46   : G4VSensitiveDetector(name), fCalCollection <<  45   :G4VSensitiveDetector(name), fCalCollection(0),
                                                   >>  46    fNumberOfCellsInZ(20),fNumberOfCellsInPhi(48)
 47 {                                                  47 {
 48   // Initialize data member.                       48   // Initialize data member.
 49   for (G4int j = 0; j < fNumberOfCellsInZ; j++ <<  49   for(G4int j=0;j<fNumberOfCellsInZ;j++)
 50     for (G4int k = 0; k < fNumberOfCellsInPhi; <<  50   for(G4int k=0;k<fNumberOfCellsInPhi;k++)
 51       fCellID[j][k] = -1;                      <<  51   {
 52     }                                          <<  52     fCellID[j][k] = -1;
                                                   >>  53   }
 53   //                                               54   //
 54   G4String HCname;                                 55   G4String HCname;
 55   collectionName.insert(HCname = "calCollectio <<  56   collectionName.insert(HCname="calCollection");
 56 }                                                  57 }
 57                                                    58 
 58 //....oooOO0OOooo........oooOO0OOooo........oo     59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 59 RE01CalorimeterSD::~RE01CalorimeterSD()            60 RE01CalorimeterSD::~RE01CalorimeterSD()
 60 {                                              <<  61 {;}
 61   ;                                            << 
 62 }                                              << 
 63                                                    62 
 64 //....oooOO0OOooo........oooOO0OOooo........oo     63 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 65 void RE01CalorimeterSD::Initialize(G4HCofThisE     64 void RE01CalorimeterSD::Initialize(G4HCofThisEvent* HCE)
 66 {                                                  65 {
 67   fCalCollection = new RE01CalorimeterHitsColl <<  66   fCalCollection = new RE01CalorimeterHitsCollection
 68   for (G4int j = 0; j < fNumberOfCellsInZ; j++ <<  67                       (SensitiveDetectorName,collectionName[0]); 
 69     for (G4int k = 0; k < fNumberOfCellsInPhi; <<  68   for(G4int j=0;j<fNumberOfCellsInZ;j++)
 70       fCellID[j][k] = -1;                      <<  69   for(G4int k=0;k<fNumberOfCellsInPhi;k++)
 71     }                                          <<  70   {
                                                   >>  71     fCellID[j][k] = -1;
                                                   >>  72   }
 72                                                    73 
 73   static G4int HCID = -1;                          74   static G4int HCID = -1;
 74   if (HCID < 0) {                              <<  75   if(HCID<0)
 75     HCID = GetCollectionID(0);                 <<  76   { HCID = GetCollectionID(0); }
 76   }                                            <<  77   HCE->AddHitsCollection( HCID, fCalCollection );
 77   HCE->AddHitsCollection(HCID, fCalCollection) << 
 78 }                                                  78 }
 79                                                    79 
 80 //....oooOO0OOooo........oooOO0OOooo........oo     80 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 81 G4bool RE01CalorimeterSD::ProcessHits(G4Step*  <<  81 G4bool RE01CalorimeterSD::ProcessHits(G4Step*aStep,G4TouchableHistory*ROhist)
 82 {                                                  82 {
 83   //***** RE05CalorimeterSD has been migrated  <<  83   if(!ROhist) return false;
 84   //***** support Readout Geometry in multi-th << 
 85   //***** is assigned to a dedicaed parallel w << 
 86   //***** a G4Step object for the parallel wor << 
 87                                                << 
 88   G4double edep = aStep->GetTotalEnergyDeposit     84   G4double edep = aStep->GetTotalEnergyDeposit();
 89   if (edep == 0.) return false;                <<  85   if(edep==0.) return false;
 90                                                    86 
 91   const G4VTouchable* ROhist = aStep->GetPreSt <<  87   G4VPhysicalVolume* physVol = ROhist->GetVolume();
                                                   >>  88   //ROhist->MoveUpHistory();
                                                   >>  89   //G4VPhysicalVolume* mothVol = ROhist->GetVolume(1);
 92   G4int copyIDinZ = ROhist->GetReplicaNumber()     90   G4int copyIDinZ = ROhist->GetReplicaNumber();
 93   G4int copyIDinPhi = ROhist->GetReplicaNumber     91   G4int copyIDinPhi = ROhist->GetReplicaNumber(1);
 94                                                    92 
 95   if (fCellID[copyIDinZ][copyIDinPhi] == -1) { <<  93   if(fCellID[copyIDinZ][copyIDinPhi]==-1)
 96     RE01CalorimeterHit* calHit =               <<  94   {
 97       new RE01CalorimeterHit(ROhist->GetVolume <<  95     RE01CalorimeterHit* calHit
 98     calHit->SetEdep(edep);                     <<  96       = new RE01CalorimeterHit
                                                   >>  97             (physVol->GetLogicalVolume(),copyIDinZ,copyIDinPhi);
                                                   >>  98     calHit->SetEdep( edep );
 99     G4AffineTransform aTrans = ROhist->GetHist     99     G4AffineTransform aTrans = ROhist->GetHistory()->GetTopTransform();
100     aTrans.Invert();                              100     aTrans.Invert();
101     calHit->SetPos(aTrans.NetTranslation());      101     calHit->SetPos(aTrans.NetTranslation());
102     calHit->SetRot(aTrans.NetRotation());         102     calHit->SetRot(aTrans.NetRotation());
103     calHit->SetTrackInformation(aStep->GetTrac    103     calHit->SetTrackInformation(aStep->GetTrack());
104     G4int icell = fCalCollection->insert(calHi << 104     G4int icell = fCalCollection->insert( calHit );
105     fCellID[copyIDinZ][copyIDinPhi] = icell -     105     fCellID[copyIDinZ][copyIDinPhi] = icell - 1;
106     if (verboseLevel > 0) {                    << 106     if(verboseLevel>0)
107       G4cout << " New Calorimeter Hit on CellI << 107     { G4cout << " New Calorimeter Hit on CellID " 
108     }                                          << 108            << copyIDinZ << " " << copyIDinPhi << G4endl; }
109   }                                               109   }
110   else {                                       << 110   else
                                                   >> 111   { 
111     (*fCalCollection)[fCellID[copyIDinZ][copyI    112     (*fCalCollection)[fCellID[copyIDinZ][copyIDinPhi]]->AddEdep(edep);
112     (*fCalCollection)[fCellID[copyIDinZ][copyI << 113     (*fCalCollection)[fCellID[copyIDinZ][copyIDinPhi]]
113     if (verboseLevel > 0) {                    << 114                              ->SetTrackInformation(aStep->GetTrack());
114       G4cout << " Energy added to CellID " <<  << 115     if(verboseLevel>0)
115     }                                          << 116     { G4cout << " Energy added to CellID " 
                                                   >> 117            << copyIDinZ << " " << copyIDinPhi << G4endl; }
116   }                                               118   }
117                                                   119 
118   return true;                                    120   return true;
119 }                                                 121 }
                                                   >> 122 
120                                                   123