Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/persistency/gdml/G04/gdml_det.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/persistency/gdml/G04/gdml_det.cc (Version 11.3.0) and /examples/extended/persistency/gdml/G04/gdml_det.cc (Version 9.6.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 /// \file persistency/gdml/G04/gdml_det.cc         26 /// \file persistency/gdml/G04/gdml_det.cc
 27 /// \brief Main program of the persistency/gdm     27 /// \brief Main program of the persistency/gdml/G04 example
 28 //                                                 28 //
 29 //                                                 29 //
                                                   >>  30 // $Id$
 30 //                                                 31 //
 31 //                                                 32 //
 32 // -------------------------------------------     33 // --------------------------------------------------------------
 33 //      GEANT 4 - gdml_det                         34 //      GEANT 4 - gdml_det
 34 //                                                 35 //
 35 // -------------------------------------------     36 // --------------------------------------------------------------
 36                                                    37 
 37 #include "FTFP_BERT.hh"                        <<  38 #include <vector>
 38 #include "G04ActionInitialization.hh"          <<  39 
                                                   >>  40 #include "G4RunManager.hh"
                                                   >>  41 #include "G4UImanager.hh"
                                                   >>  42 
                                                   >>  43 #include "G4LogicalVolumeStore.hh"
                                                   >>  44 #include "G4TransportationManager.hh"
                                                   >>  45 #include "G4SDManager.hh"
                                                   >>  46 
                                                   >>  47 #include "G04PrimaryGeneratorAction.hh"
 39 #include "G04DetectorConstruction.hh"              48 #include "G04DetectorConstruction.hh"
                                                   >>  49 #include "G04PhysicsList.hh"
 40 #include "G04SensitiveDetector.hh"                 50 #include "G04SensitiveDetector.hh"
 41                                                    51 
 42 #include "G4GDMLParser.hh"                     <<  52 #ifdef G4VIS_USE
 43 #include "G4RunManagerFactory.hh"              << 
 44 #include "G4UIExecutive.hh"                    << 
 45 #include "G4UImanager.hh"                      << 
 46 #include "G4VisExecutive.hh"                       53 #include "G4VisExecutive.hh"
                                                   >>  54 #endif
 47                                                    55 
 48 #include <vector>                              <<  56 #ifdef G4UI_USE
                                                   >>  57 #include "G4UIExecutive.hh"
                                                   >>  58 #endif
                                                   >>  59 
                                                   >>  60 #include "G4GDMLParser.hh"
 49                                                    61 
 50 int main(int argc, char** argv)                <<  62 int main(int argc,char **argv)
 51 {                                                  63 {
 52   G4cout << G4endl;                            <<  64    G4cout << G4endl;
 53   G4cout << "Usage: gdml_det <intput_gdml_file <<  65    G4cout << "Usage: gdml_det <intput_gdml_file:mandatory>"
 54   G4cout << G4endl;                            <<  66           << G4endl;
 55                                                <<  67    G4cout << G4endl;
 56   if (argc < 2) {                              <<  68 
 57     G4cout << "Error! Mandatory input file is  <<  69    if (argc<2)
 58     G4cout << G4endl;                          <<  70    {
 59     return -1;                                 <<  71       G4cout << "Error! Mandatory input file is not specified!" << G4endl;
 60   }                                            <<  72       G4cout << G4endl;
 61                                                <<  73       return -1;
 62   // Detect interactive mode (if only one argu <<  74    }
 63   //                                           <<  75 
 64   G4UIExecutive* ui = 0;                       <<  76    G4GDMLParser parser;
 65   if (argc == 2) {                             <<  77    parser.Read(argv[1]);
 66     ui = new G4UIExecutive(argc, argv);        <<  78 
 67   }                                            <<  79    
 68                                                <<  80    G4RunManager* runManager = new G4RunManager;
 69   G4GDMLParser parser;                         <<  81 
 70   parser.Read(argv[1]);                        <<  82    runManager->SetUserInitialization(new G04DetectorConstruction(
 71                                                <<  83                                      parser.GetWorldVolume()));
 72   auto* runManager = G4RunManagerFactory::Crea <<  84    runManager->SetUserInitialization(new G04PhysicsList);
 73                                                <<  85    runManager->SetUserAction(new G04PrimaryGeneratorAction);
 74   runManager->SetUserInitialization(new G04Det <<  86 
 75   runManager->SetUserInitialization(new FTFP_B <<  87    runManager->Initialize();
 76                                                <<  88 
 77   // User action initialization                <<  89    G4UImanager* UImanager = G4UImanager::GetUIpointer();
 78   runManager->SetUserInitialization(new G04Act <<  90 
 79   runManager->Initialize();                    <<  91    //------------------------------------------------ 
 80                                                <<  92    // Sensitive detectors
 81   // Initialize visualization                  <<  93    //------------------------------------------------ 
 82   G4VisManager* visManager = new G4VisExecutiv <<  94    
 83   visManager->Initialize();                    <<  95    G4SDManager* SDman = G4SDManager::GetSDMpointer();
 84                                                <<  96    
 85   // Get the pointer to the User Interface man <<  97    G4String trackerChamberSDname = "Tracker";
 86   G4UImanager* UImanager = G4UImanager::GetUIp <<  98    G04SensitiveDetector* aTrackerSD = new G04SensitiveDetector(trackerChamberSDname);
 87                                                <<  99    SDman->AddNewDetector( aTrackerSD );
 88   // Process macro or start UI session         << 100  
 89   if (!ui)  // batch mode                      << 101    ///////////////////////////////////////////////////////////////////////
 90   {                                            << 102    //
 91     G4String command = "/control/execute ";    << 103    // Example how to retrieve Auxiliary Information for sensitive detector
 92     G4String fileName = argv[2];               << 104    //
 93     UImanager->ApplyCommand(command + fileName << 105    const G4GDMLAuxMapType* auxmap = parser.GetAuxMap();
 94   }                                            << 106    std::cout << "Found " << auxmap->size()
 95   else  // interactive mode                    << 107              << " volume(s) with auxiliary information."
 96   {                                            << 108              << G4endl << G4endl;
 97     UImanager->ApplyCommand("/control/execute  << 109    for(G4GDMLAuxMapType::const_iterator iter=auxmap->begin();
 98     ui->SessionStart();                        << 110        iter!=auxmap->end(); iter++) 
 99     delete ui;                                 << 111    {
100   }                                            << 112      G4cout << "Volume " << ((*iter).first)->GetName()
                                                   >> 113             << " has the following list of auxiliary information: "
                                                   >> 114             << G4endl << G4endl;
                                                   >> 115      for (G4GDMLAuxListType::const_iterator vit=(*iter).second.begin();
                                                   >> 116           vit!=(*iter).second.end(); vit++)
                                                   >> 117      {
                                                   >> 118        std::cout << "--> Type: " << (*vit).type
                                                   >> 119                  << " Value: " << (*vit).value << std::endl;
                                                   >> 120      }
                                                   >> 121    }
                                                   >> 122    G4cout << G4endl;
                                                   >> 123 
                                                   >> 124    // The same as above, but now we are looking for
                                                   >> 125    // sensitive detectors setting them for the volumes
                                                   >> 126 
                                                   >> 127    for(G4GDMLAuxMapType::const_iterator iter=auxmap->begin();
                                                   >> 128        iter!=auxmap->end(); iter++) 
                                                   >> 129    {
                                                   >> 130      G4cout << "Volume " << ((*iter).first)->GetName()
                                                   >> 131             << " has the following list of auxiliary information: "
                                                   >> 132             << G4endl << G4endl;
                                                   >> 133      for (G4GDMLAuxListType::const_iterator vit=(*iter).second.begin();
                                                   >> 134           vit!=(*iter).second.end();vit++)
                                                   >> 135      {
                                                   >> 136        if ((*vit).type=="SensDet")
                                                   >> 137        {
                                                   >> 138          G4cout << "Attaching sensitive detector " << (*vit).value
                                                   >> 139                 << " to volume " << ((*iter).first)->GetName()
                                                   >> 140                 <<  G4endl << G4endl;
                                                   >> 141 
                                                   >> 142          G4VSensitiveDetector* mydet = SDman->FindSensitiveDetector((*vit).value);
                                                   >> 143          if(mydet) 
                                                   >> 144          {
                                                   >> 145            G4LogicalVolume* myvol = (*iter).first;
                                                   >> 146            myvol->SetSensitiveDetector(mydet);
                                                   >> 147          }
                                                   >> 148          else
                                                   >> 149          {
                                                   >> 150            G4cout << (*vit).value << " detector not found" << G4endl;
                                                   >> 151          }
                                                   >> 152        }
                                                   >> 153      }
                                                   >> 154    }
                                                   >> 155    //
                                                   >> 156    // End of Auxiliary Information block
                                                   >> 157    //
                                                   >> 158    ////////////////////////////////////////////////////////////////////////
                                                   >> 159 
                                                   >> 160 #ifdef G4VIS_USE
                                                   >> 161      G4VisManager* visManager = new G4VisExecutive;
                                                   >> 162      visManager->Initialize();
                                                   >> 163 #endif
                                                   >> 164 
                                                   >> 165    if (argc==3)   // batch mode  
                                                   >> 166    {
                                                   >> 167      G4String command = "/control/execute ";
                                                   >> 168      G4String fileName = argv[2];
                                                   >> 169      UImanager->ApplyCommand(command+fileName);
                                                   >> 170    }
                                                   >> 171    else           // interactive mode
                                                   >> 172    {
                                                   >> 173 #ifdef G4UI_USE
                                                   >> 174      G4UIExecutive* ui = new G4UIExecutive(argc, argv);
                                                   >> 175 #ifdef G4VIS_USE
                                                   >> 176      UImanager->ApplyCommand("/control/execute vis.mac");     
                                                   >> 177 #endif
                                                   >> 178      ui->SessionStart();
                                                   >> 179      delete ui;
                                                   >> 180 #endif
                                                   >> 181    }
                                                   >> 182 
                                                   >> 183 #ifdef G4VIS_USE
                                                   >> 184    delete visManager;
                                                   >> 185 #endif
                                                   >> 186    delete runManager;
101                                                   187 
102   delete visManager;                           << 188    return 0;
103   delete runManager;                           << 
104 }                                                 189 }
105                                                   190