Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/examples/extended/eventgenerator/HepMC/HepMCEx01/src/HepMCG4Interface.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/eventgenerator/HepMC/HepMCEx01/src/HepMCG4Interface.cc (Version 11.3.0) and /examples/extended/eventgenerator/HepMC/HepMCEx01/src/HepMCG4Interface.cc (Version 8.0.p1)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                    <<   3 // * DISCLAIMER                                                       *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th <<   5 // * The following disclaimer summarizes all the specific disclaimers *
  6 // * the Geant4 Collaboration.  It is provided <<   6 // * of contributors to this software. The specific disclaimers,which *
  7 // * conditions of the Geant4 Software License <<   7 // * govern, are listed with their locations in:                      *
  8 // * LICENSE and available at  http://cern.ch/ <<   8 // *   http://cern.ch/geant4/license                                  *
  9 // * include a list of copyright holders.      << 
 10 // *                                                9 // *                                                                  *
 11 // * Neither the authors of this software syst     10 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     11 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     12 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     13 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file  <<  14 // * use.                                                             *
 16 // * for the full disclaimer and the limitatio << 
 17 // *                                               15 // *                                                                  *
 18 // * This  code  implementation is the result  <<  16 // * This  code  implementation is the  intellectual property  of the *
 19 // * technical work of the GEANT4 collaboratio <<  17 // * GEANT4 collaboration.                                            *
 20 // * By using,  copying,  modifying or  distri <<  18 // * By copying,  distributing  or modifying the Program (or any work *
 21 // * any work based  on the software)  you  ag <<  19 // * based  on  the Program)  you indicate  your  acceptance of  this *
 22 // * use  in  resulting  scientific  publicati <<  20 // * statement, and all its terms.                                    *
 23 // * acceptance of all terms of the Geant4 Sof << 
 24 // *******************************************     21 // ********************************************************************
 25 //                                                 22 //
 26 /// \file eventgenerator/HepMC/HepMCEx01/src/H <<  23 // ====================================================================
 27 /// \brief Implementation of the HepMCG4Interf << 
 28 //                                                 24 //
                                                   >>  25 //   HepMCG4Interface.cc
                                                   >>  26 //   $Id: HepMCG4Interface.cc,v 1.2 2003/12/09 15:30:52 gunter Exp $
 29 //                                                 27 //
 30                                                <<  28 // ====================================================================
 31 #include "HepMCG4Interface.hh"                     29 #include "HepMCG4Interface.hh"
 32                                                    30 
 33 #include "G4Event.hh"                          <<  31 #include "G4RunManager.hh"
 34 #include "G4LorentzVector.hh"                      32 #include "G4LorentzVector.hh"
 35 #include "G4PhysicalConstants.hh"              <<  33 #include "G4Event.hh"
 36 #include "G4PrimaryParticle.hh"                    34 #include "G4PrimaryParticle.hh"
 37 #include "G4PrimaryVertex.hh"                      35 #include "G4PrimaryVertex.hh"
 38 #include "G4RunManager.hh"                     << 
 39 #include "G4SystemOfUnits.hh"                  << 
 40 #include "G4TransportationManager.hh"              36 #include "G4TransportationManager.hh"
 41                                                    37 
 42 //....oooOO0OOooo........oooOO0OOooo........oo <<  38 ////////////////////////////////////
 43 HepMCG4Interface::HepMCG4Interface() : hepmcEv <<  39 HepMCG4Interface::HepMCG4Interface()
                                                   >>  40   : hepmcEvent(0)
                                                   >>  41 ////////////////////////////////////
                                                   >>  42 {
                                                   >>  43 }
 44                                                    44 
 45 //....oooOO0OOooo........oooOO0OOooo........oo <<  45 /////////////////////////////////////
 46 HepMCG4Interface::~HepMCG4Interface()              46 HepMCG4Interface::~HepMCG4Interface()
                                                   >>  47 /////////////////////////////////////
 47 {                                                  48 {
 48   delete hepmcEvent;                               49   delete hepmcEvent;
 49 }                                                  50 }
 50                                                    51 
 51 //....oooOO0OOooo........oooOO0OOooo........oo <<  52 /////////////////////////////////////////////////////////
 52 G4bool HepMCG4Interface::CheckVertexInsideWorl <<  53 G4bool HepMCG4Interface::CheckVertexInsideWorld
                                                   >>  54                          (const G4ThreeVector& pos) const
                                                   >>  55 /////////////////////////////////////////////////////////
 53 {                                                  56 {
 54   G4Navigator* navigator =                     <<  57   G4Navigator* navigator= G4TransportationManager::GetTransportationManager()
 55     G4TransportationManager::GetTransportation <<  58                                                  -> GetNavigatorForTracking();
 56                                                    59 
 57   G4VPhysicalVolume* world = navigator->GetWor <<  60   G4VPhysicalVolume* world= navigator-> GetWorldVolume();
 58   G4VSolid* solid = world->GetLogicalVolume()- <<  61   G4VSolid* solid= world-> GetLogicalVolume()-> GetSolid();
 59   EInside qinside = solid->Inside(pos);        <<  62   EInside qinside= solid-> Inside(pos);
 60                                                <<  63 
 61   if (qinside != kInside)                      <<  64   if( qinside != kInside) return false;
 62     return false;                              <<  65   else return true;
 63   else                                         << 
 64     return true;                               << 
 65 }                                                  66 }
 66                                                    67 
 67 //....oooOO0OOooo........oooOO0OOooo........oo <<  68 ////////////////////////////////////////////////////////////////
 68 void HepMCG4Interface::HepMC2G4(const HepMC::G <<  69 void HepMCG4Interface::HepMC2G4(const HepMC::GenEvent* hepmcevt, 
                                                   >>  70         G4Event* g4event)
                                                   >>  71 ////////////////////////////////////////////////////////////////
 69 {                                                  72 {
 70   for (HepMC::GenEvent::vertex_const_iterator  <<  73   for(HepMC::GenEvent::vertex_const_iterator vitr= hepmcevt->vertices_begin();
 71        vitr != hepmcevt->vertices_end(); ++vit <<  74       vitr != hepmcevt->vertices_end(); ++vitr ) { // loop for vertex ...
 72   {  // loop for vertex ...                    << 
 73                                                    75 
 74     // real vertex?                                76     // real vertex?
 75     G4bool qvtx = false;                       <<  77     G4bool qvtx=false;
 76     for (HepMC::GenVertex::particle_iterator p <<  78     for (HepMC::GenVertex::particle_iterator 
 77          pitr != (*vitr)->particles_end(HepMC: <<  79      pitr= (*vitr)->particles_begin(HepMC::children);
 78     {                                          <<  80    pitr != (*vitr)->particles_end(HepMC::children); ++pitr) {
 79       if (!(*pitr)->end_vertex() && (*pitr)->s <<  81 
 80         qvtx = true;                           <<  82       if (!(*pitr)->end_vertex() && (*pitr)->status()==1) {
 81         break;                                 <<  83   qvtx=true;
                                                   >>  84   break;
 82       }                                            85       }
 83     }                                              86     }
 84     if (!qvtx) continue;                           87     if (!qvtx) continue;
 85                                                    88 
 86     // check world boundary                        89     // check world boundary
 87     HepMC::FourVector pos = (*vitr)->position( <<  90     G4LorentzVector xvtx= (*vitr)-> position();
 88     G4LorentzVector xvtx(pos.x(), pos.y(), pos <<  91     if (! CheckVertexInsideWorld(xvtx.vect()*mm)) continue;
 89     if (!CheckVertexInsideWorld(xvtx.vect() *  << 
 90                                                    92 
 91     // create G4PrimaryVertex and associated G     93     // create G4PrimaryVertex and associated G4PrimaryParticles
 92     G4PrimaryVertex* g4vtx =                   <<  94     G4PrimaryVertex* g4vtx= 
 93       new G4PrimaryVertex(xvtx.x() * mm, xvtx. <<  95       new G4PrimaryVertex(xvtx.x()*mm, xvtx.y()*mm, xvtx.z()*mm, 
 94                                                <<  96         xvtx.t()*mm/c_light);
 95     for (HepMC::GenVertex::particle_iterator v <<  97 
 96          vpitr != (*vitr)->particles_end(HepMC <<  98     for (HepMC::GenVertex::particle_iterator 
 97     {                                          <<  99      vpitr= (*vitr)->particles_begin(HepMC::children);
 98       if ((*vpitr)->status() != 1) continue;   << 100    vpitr != (*vitr)->particles_end(HepMC::children); ++vpitr) {
 99                                                << 101 
100       G4int pdgcode = (*vpitr)->pdg_id();      << 102       if( (*vpitr)->status() != 1 ) continue;
101       pos = (*vpitr)->momentum();              << 103 
102       G4LorentzVector p(pos.px(), pos.py(), po << 104       G4int pdgcode= (*vpitr)-> pdg_id();
103       G4PrimaryParticle* g4prim =              << 105       G4LorentzVector p= (*vpitr)-> momentum();
104         new G4PrimaryParticle(pdgcode, p.x() * << 106       G4PrimaryParticle* g4prim= 
                                                   >> 107   new G4PrimaryParticle(pdgcode, p.x()*GeV, p.y()*GeV, p.z()*GeV);
105                                                   108 
106       g4vtx->SetPrimary(g4prim);               << 109       g4vtx-> SetPrimary(g4prim);
107     }                                             110     }
108     g4event->AddPrimaryVertex(g4vtx);          << 111     g4event-> AddPrimaryVertex(g4vtx);
109   }                                            << 112   } 
110 }                                              << 113 } 
                                                   >> 114 
111                                                   115 
112 //....oooOO0OOooo........oooOO0OOooo........oo << 116 ///////////////////////////////////////////////////////
113 HepMC::GenEvent* HepMCG4Interface::GenerateHep    117 HepMC::GenEvent* HepMCG4Interface::GenerateHepMCEvent()
                                                   >> 118 ///////////////////////////////////////////////////////
114 {                                                 119 {
115   HepMC::GenEvent* aevent = new HepMC::GenEven << 120   HepMC::GenEvent* aevent= new HepMC::GenEvent();
116   return aevent;                                  121   return aevent;
117 }                                                 122 }
118                                                   123 
119 //....oooOO0OOooo........oooOO0OOooo........oo << 124 //////////////////////////////////////////////////////////////
120 void HepMCG4Interface::GeneratePrimaryVertex(G    125 void HepMCG4Interface::GeneratePrimaryVertex(G4Event* anEvent)
                                                   >> 126 //////////////////////////////////////////////////////////////
121 {                                                 127 {
122   // delete previous event object                 128   // delete previous event object
123   delete hepmcEvent;                              129   delete hepmcEvent;
124                                                   130 
125   // generate next event                          131   // generate next event
126   hepmcEvent = GenerateHepMCEvent();           << 132   hepmcEvent= GenerateHepMCEvent();
127   if (!hepmcEvent) {                           << 133   if(! hepmcEvent) {
128     G4cout << "HepMCInterface: no generated pa << 134     G4cout << "HepMCInterface: no generated particles. run terminated..." 
129     G4RunManager::GetRunManager()->AbortRun(); << 135            << G4endl;
                                                   >> 136     G4RunManager::GetRunManager()-> AbortRun();
130     return;                                       137     return;
131   }                                               138   }
132   HepMC2G4(hepmcEvent, anEvent);                  139   HepMC2G4(hepmcEvent, anEvent);
133 }                                                 140 }
                                                   >> 141 
134                                                   142