Integrate geant4star with mainline

StRoot/StBFChain/BigFullChain.h

@@ -1211,6 +1211,7 @@ Bfc_st BFC[] = { // standard chains
   {"McEvOut"     ,""  ,"","StMcEvent,Tree"                       ,"","","Write StMcEvent to StTree",kFALSE},
   {"EvOut"       ,""  ,"","Tree"                                   ,"","","Write StEvent to StTree",kFALSE},
   {"GeantOut"    ,""  ,"","Tree"                                ,"","","Write g2t tables to StTree",kFALSE},
+  {"Geant4Out"   ,""  ,"","Tree"                        ,"","","Write g2t tables to StTree from G4",kFALSE},
   {"Simu"        ,""  ,"","" ,"","","Simulated Data, force to use Db time stamp from used geometry",kFALSE},
   {"HitsBranch"  ,""  ,"",""  ,"","","take out points from dst branch and put them into HitsBranch",kFALSE},
   {"paw"         ,""  ,"",""                                      ,"","","Allocate memory for pawc",kFALSE},
@@ -1376,12 +1377,64 @@ Bfc_st BFC[] = { // standard chains
   {"------------","-----------","-----------","------------------------------------------","","","",kFALSE},
   {"Generators  ","-----------","-----------","------------------------------------------","","","",kFALSE},
   {"------------","-----------","-----------","------------------------------------------","","","",kFALSE},
-  {"stargen",     "","", "gen_T,sim_T"/*+++*/,                     "", "libVMC.so,libStarGeneratorUtil.so,libStarGeneratorEvent.so,libStarGeneratorBase.so,libStarGeneratorFilt.so,libMathMore.so","STAR Generator BASE",false},
+  {"stargen.filt", "","", "gen_T,sim_T"/*+++*/,                     "", "libVMC.so,libStarGeneratorFilt.so,libMathMore.so","STAR Generator BASE",false},
+  {"stargen.util", "","", "gen_T,sim_T"/*+++*/,                     "", "libVMC.so,libStarGeneratorUtil.so,libMathMore.so","STAR Generator BASE",false},
+  {"stargen.event","","", "gen_T,sim_T"/*+++*/,                     "", "libVMC.so,libStarGeneratorEvent.so,libMathMore.so","STAR Generator BASE",false},
+  {"stargen.base", "","", "gen_T,sim_T"/*+++*/,                     "", "libVMC.so,libStarGeneratorBase.so,libMathMore.so","STAR Generator BASE",false},
+  {"stargen",     "","",  "stargen.util,stargen.event,stargen.base"/*+++*/,                     "", "","STAR Generator BASE",false},
+  {"pythia6.4.28","","","stargen", "", "Pythia6_4_28.so",       "Load Pythia 6.4.28 generator", false },
   {"pythia8.1.86","","","stargen", "", "Pythia8_1_86.so",       "Load Pythia 8.1.86 generator", false },
   {"pythia8.2.35","","","stargen", "", "Pythia8_2_35.so",       "Load Pythia 8.1.86 generator", false },
   {"hijing1.383" ,"","","stargen", "", "Hijing1_383.so",        "Load Hijing  1.383 generator", false },
   {"kinematics"  ,"","","stargen", "", "Kinematics.so",         "Load STAR Particle Gun", false },
   {"genreader"   ,"","","stargen", "", "StarGenEventReader.so", "Load STAR Gen Event Reader", false },
+  {"------------","-----------","-----------","------------------------------------------","","","",kFALSE},
+  {"Geant4star  ","-----------","-----------","------------------------------------------","","","",kFALSE},
+  {"------------","-----------","-----------","------------------------------------------","","","",kFALSE},
+  {"g4clhep",     "","", "",                      "", "libG4clhep.so","Load G4 libs",false},
+  {"g4global",    "","", "g4clhep",               "", "libG4global.so","Load G4 libs",false},
+  {"g4intercoms", "","", "g4global",              "", "libG4intercoms.so","Load G4 libs",false},
+  {"g4materials", "","", "g4intercoms",           "", "libG4materials.so","Load G4 libs",false},
+  {"g4interfaces","","", "g4intercoms",           "", "libG4interfaces.so","Load G4 libs",false},
+  {"g4graphics",  "","", "g4intercoms",           "", "libG4graphics_reps.so","Load G4 libs",false},
+  {"g4geometry",  "","", "g4graphics,g4materials","", "libG4geometry.so","Load G4 libs",false},
+  {"g4particles", "","", "g4geometry",            "", "libG4particles.so", "Load G4",false},
+  {"g4track"    , "","", "g4particles",           "", "libG4track.so",     "Load G4",false},
+  {"g4digits",    "","", "g4track",               "", "libG4zlib.so,libG4digits_hits.so", "Load G4",false},
+  {"g4processes", "","", "g4digits",              "", "libG4processes.so", "Load G4",false},
+
+  {"g4tracking" , "","", "g4processes",           "", "libG4tracking.so",  "Load G4",false},
+  {"g4event",     "","", "g4tracking",            "", "libG4event.so",       "Load G4",false},
+  {"g4run",       "", "","g4event",               "", "libG4run.so",         "Load G4",false},
+  {"g4physics",   "", "", "g4run",                "", "libG4physicslists.so", "Load G4",false},
+
+  {"g4modeling",  "", "", "g4physics",            "", "libG4modeling.so", "Load G4", false},
+  {"g4persist",   "", "", "",                     "", "libG4persistency.so", "Geant4 persistency",false},
+  {"g4errprop", "", "", "",                     "", "libG4error_propagation.so", "Geant4 error propagation",false},
+  {"g4readout", "", "", "",                     "", "libG4readout.so","G4 readout",false},
+  {"g4parmodels", "", "", "",                   "", "libG4parmodels.so", "G4 par models", false},
+  {"g4analysis",  "","", "g4errprop,g4readout,g4parmodels",                 "", "libG4analysis.so",false}, 
+
+  {"g4visual",    "", "", "g4modeling",           "", "libG4vis_management.so", "Load G4", false},  
+  {"g4ascii",     "", "", "g4visual",             "", "libG4Tree.so", "G4 ASCII Tree", false},
+  {"g4vrml",      "", "", "g4visual",             "", "libG4VRML.so", "G4 ASCII Tree", false},
+  {"g4dawn",      "", "", "g4visual",             "", "libG4FR.so",   "G4 DAWN", false},
+  {"g4heprep",    "", "", "g4visual",             "", "libG4zlib.so,libG4visHepRep.so",   "G4 DAWN", false},
+  {"g4raytracer", "", "", "g4visual",             "", "libG4RayTracer.so",   "G4 DAWN", false},
+  {"g4gmocren",   "", "", "g4visual",             "", "libG4GMocren.so", "G4 ASCII Tree", false},
+  {"g4opengl",    "", "", "g4visual",             "", "libG4gl2ps.so,libG4OpenGL.so", "G4 ASCII Tree", false},
+  {"g4geant3",    "", "", "",                     "", "libG3toG4.so", "Load g3 to g4 support", false },
+  {"geant4",      "","", "g4physics,g4interfaces,g4visual","","","Load G4 libs", false},
+
+  {"g4root",      "","", "g4analysis",       "", "libg4root.so", "Load g4root support", false },
+  {"geant4vmc",   "","", "geant4,g4root,g4ascii,g4dawn,g4heprep,g4raytracer,g4vrml,g4gmocren,g4opengl,g4persist", "", 
+                         "libVMC.so,libgeant4vmc.so", "Load G4 VMC libs", false},
+
+  {"ag4ml",       "","", "g4geometry, agml"            ,"","", "AgML + g4",false},
+  {"stargen",     "","", "gen_T,sim_T"/*+++*/,                     "", "libStarGeneratorUtil.so,libStarGeneratorEvent.so,libStarGeneratorBase.so,libStarGeneratorFilt.so,libMathMore.so","STAR Generator BASE",false},
+  {"geant4mk",    "","", "stargen,geant4vmc", "", "StGeant4Maker.so,StarMagField.so", "Load G4 VMC libs", false},
+  {"geant321", "", "", "",                 "", "/star/simu/simu/jwebb/GEANT3/lib64/libgeant321.so", "Load G3 VMC library", false},
+
   {"------------","-----------","-----------","------------------------------------------","","","",kFALSE},
   {"I/O Makers  ","-----------","-----------","------------------------------------------","","","",kFALSE},
   {"------------","-----------","-----------","------------------------------------------","","","",kFALSE},
@@ -1587,18 +1640,6 @@ Bfc_st BFC[] = { // standard chains
                                                                                                     kFALSE},
   {"fstUtil"  , "", "", "",          "",       "StFstUtil",    "Fst Utilities",                     
                                                                                                     kFALSE},
-
-  {"fstFastSim","fstChain","","StMcEvent,StEvent","StFstFastSimMaker","StFstSimMaker","FST fast simulator", 
-                                                                                                    kFALSE},
-  {"fstRawHit", "", "fstChain", "fstUtil,fstDb","StFstRawHitMaker", "StFstRawHitMaker","FST raw hit maker", 
-                                                                                                    kFALSE},
-  {"fstCluster", "", "fstChain", "fstRawHit","StFstClusterMaker",  "StFstClusterMaker","FST Cluster maker", 
-                                                                                                    kFALSE},
-  {"fstHit",     "", "fstChain", "event,fstCluster", "StFstHitMaker",       "StFstHitMaker","FST Hit Maker",
-                                                                                                    kFALSE},
-  {"fstUtil"  , "", "fstChain", "",          "",       "StFstUtil",    "Fst Utilities",                     
-                                                                                                    kFALSE},
-
   {"ssddat"      ,"","","ssd_daq"                             ,"","","SSD full chain for Real Data",kFALSE},
   {"sstdat"      ,"","","sst_daq"                             ,"","","SST full chain for Real Data",kFALSE},
   {"ssd_daq","","","ssdCalDb,svt_T,-sls,-spa,ssdUtil","StSsdDaqMaker","StSsdDaqMaker","... SSD Daq",kFALSE},

StRoot/StBFChain/StBFChain.cxx

@@ -1942,6 +1942,7 @@ void StBFChain::SetTreeOptions()
     treeMk->IntoBranch("McEventBranch","StMcEvent");
   }
   if (GetOption("GeantOut")) treeMk->IntoBranch("geantBranch","geant");
+  if (GetOption("Geant4Out")) treeMk->IntoBranch("geantBranch","geant4star");
   if (GetOption("AllEvent")) {
     if (GetOption("fzin")   ||
 	GetOption("ntin")   ||

StRoot/StBTofUtil/StVpdSimConfig.h

@@ -11,6 +11,7 @@
 
 #include <iostream>
 #include <fstream>
+#include "StMessMgr.h"
 #include <vector>
 #include "St_db_Maker/St_db_Maker.h"
 #include "tables/St_vpdSimParams_Table.h"

StRoot/StChain/StMaker.cxx

@@ -17,7 +17,7 @@
  * The method will be described.
  *                                                                     
  */
-#define STAR_LOGGER 1
+#define STAR_LOGGER 1 
 #include "StMaker.h"
 #include "Stiostream.h"
 #include <stdio.h>

StRoot/StGeant4Maker/AgMLBTofVolumeId.h

@@ -0,0 +1,38 @@
+#ifndef __AgMLBtofVolumeId_h__
+#define __AgMLBtofVolumeId_h__
+
+#include <StarVMC/StarAgmlLib/AgMLExtension.h>
+#include <StMessMgr.h>
+
+class AgMLBtofVolumeId : public AgMLVolumeId {
+public:
+
+  AgMLBtofVolumeId(){}
+  // Applies to btog.version = 8 with btog.choice =13 (run 13 onwards)
+
+  virtual int id( int* numbv ) const { 
+
+    int rileft = numbv[0];
+    int sector = numbv[1];
+    int module = numbv[2];
+    int layer  = numbv[3];
+
+    // Adjust for GMT modules
+    if ( rileft==1 ) {
+      if ( sector== 8 )   module += 4;
+      if ( sector==23 )   module += 4;
+    }
+    else if ( rileft==2 ) {
+      if ( sector == 33 ) module += 4;
+      if ( sector == 48 ) module += 4;
+    }
+
+    int _id = layer + 10 * (module + 100 * (sector+100*rileft) );
+
+    return _id;
+
+  };
+};
+
+
+#endif

StRoot/StGeant4Maker/AgMLEEmcVolumeId.h

@@ -0,0 +1,70 @@
+#ifndef __AgMLEEmcVolumeId_h__
+#define __AgMLEEmcVolumeId_h__
+
+#include <StarVMC/StarAgmlLib/AgMLExtension.h>
+#include <TGeoNavigator.h>
+#include <StMessMgr.h>
+#include <TLorentzVector.h>
+#include <TVirtualMC.h>
+#include <TMath.h>
+
+#include <TString.h>
+#include <StMessMgr.h>
+
+class AgMLEEmcVolumeId : public AgMLVolumeId {
+
+  const int onoff    = 1;
+  const int fillmode = 3;
+  const int sectormap[2][6] = { 
+    { 4, 5, 6, 7, 8, 9}, 
+    {10,11,12, 1, 2, 3} 
+  };
+
+
+public:
+
+  AgMLEEmcVolumeId() : AgMLVolumeId() { /* nada */ };
+
+  virtual int id( int* numbv ) const { 
+
+    // cd == ESCI
+    int rileft = onoff;
+    int shift  = 0;
+
+    int wheel   = numbv[0];  
+    int section = numbv[1];
+    int idx     = numbv[2];
+    int phi30d  = sectormap[wheel-1][idx-1]; // sector number
+    int subsec  = numbv[3]; // subsection (i.e. layer) within each section
+    int sublay  = numbv[4]; // ???
+    int phi     = numbv[5]; // phibin in sector
+    int eta     = numbv[6];
+
+    int depth   = subsec + 3 * ( section - 1 );
+
+    //    int  volumeid = 100000*rileft + 1000*(5*(phi30d-1)+phi) + 10*eta + depth;
+
+    int volumeid = 100000 * rileft 
+                 +   1000 * ( 5 * ( phi30d-1 ) + phi )
+                 +     10 *    eta
+                 +           depth;
+
+    // LOG_INFO << Form("AgML EEmc VolumdId = wheel=%i section=%i idx=%i phi30d=%i subsection=%i %i %i %i %i",
+    // 		     numbv[0],
+    // 		     numbv[1],
+    // 		     numbv[2],
+    // 		     phi30d,
+    // 		     numbv[3],
+    // 		     numbv[4],
+    // 		     numbv[5], 
+    // 		     numbv[6],
+    // 		     volumeid ) << endm;    
+
+
+    return volumeid;
+
+  };
+};
+
+
+#endif

StRoot/StGeant4Maker/AgMLEmcVolumeId.h

@@ -0,0 +1,66 @@
+#ifndef __AgMLEmcVolumeId_h__
+#define __AgMLEmcVolumeId_h__
+
+#include <StarVMC/StarAgmlLib/AgMLExtension.h>
+#include <TGeoNavigator.h>
+#include <StMessMgr.h>
+#include <TLorentzVector.h>
+#include <TVirtualMC.h>
+#include <TMath.h>
+
+#include <TString.h>
+#include <StMessMgr.h>
+
+class AgMLEmcVolumeId : public AgMLVolumeId {
+
+  const int    numberOfEtaRings = 20;
+
+public:
+
+  AgMLEmcVolumeId() : AgMLVolumeId() { /* nada */ };
+
+  virtual int id( int* numbv ) const { 
+
+
+    int rileft = numbv[0]; // east barrel vs west barrel
+    int phi    = numbv[1]; // module
+    int superl = numbv[2]; // before / after SMD layer
+
+    TLorentzVector _direction(0,0,0,0);
+    TVirtualMC::GetMC()->TrackPosition( _direction );
+
+    double xg[4], xl[4];
+    _direction.GetXYZT( xg );
+    double pseudoRapidity = TMath::Abs(_direction.Eta());        
+
+    // Obtain local coordinates from global coordinates
+    TVirtualMC::GetMC()->Gmtod( xg, xl, 1 );
+
+    // Calculate the eta ring and submodule
+    int eta_tow = ( pseudoRapidity * double(numberOfEtaRings) ) + 1.0;
+    int phi_sub = ( xl[1]>= 0 )? 1 : 0; // (-13,0)=0, (0,13)=1
+
+    //    LOG_INFO << Form("eta=%f xlocal=%f %f %f | phi_sub=%i eta_tow=%i", pseudoRapidity, xl[0], xl[1], xl[2],phi_sub,eta_tow) << endm;
+
+    int volumeid = -999;
+
+    if ( rileft==1 ) {
+      phi = 60 - phi + 1;
+      if ( phi_sub==0 ) {
+	phi_sub=2;
+      }
+    }
+    else {
+      phi = 60+phi;
+      phi_sub = phi_sub + 1;      
+    }
+
+    volumeid=10000000*rileft+100000*eta_tow+100*phi+10*phi_sub+superl;
+
+    return volumeid;
+
+  };
+};
+
+
+#endif

StRoot/StGeant4Maker/AgMLEpdVolumeId.h

@@ -0,0 +1,40 @@
+#ifndef __AgMLEpdVolumeId_h__
+#define __AgMLEpdVolumeId_h__
+
+#include <StarVMC/StarAgmlLib/AgMLExtension.h>
+#include <StMessMgr.h>
+
+class AgMLEpdVolumeId : public AgMLVolumeId {
+public:
+
+  AgMLEpdVolumeId(){}
+
+  virtual int id( int* numbv ) const { 
+
+    int epdm = numbv[0]; // 1=east, 2=west
+    int epss = numbv[1]; // 1 for PP1, 2 for PP2, PP-postion 1'o,2'o clock etc
+    int epdt = numbv[2]; // 1:T1 trap, 2:T1 Triangular, 3:T2 Thin, 4:T3 Thick
+
+    /*
+
+      " EPD volume_id " 
+      " 100,000 : east or west "
+      "   1,000 : Position clock wise, 1 to 12 "
+      "      10 : Tile number 1 to 31, refer EPD Drupal page"
+      "       1 : 1 T1 trap or T2 thin; 0 T1 triangular or T2 thick
+
+     */
+
+    int _id = 0;
+    _id += 100000 * epdm;
+    _id +=   1000 * epss;
+    _id +=     10 * ( epdt%2 + epdt/2 );
+    _id +=      1 * ( epdt%2 );
+
+    return _id;
+
+  };
+};
+
+
+#endif

StRoot/StGeant4Maker/AgMLFstVolumeId.h

@@ -0,0 +1,28 @@
+#ifndef __AgMLFstVolumeId_h__
+#define __AgMLFstVolumeId_h__
+
+#include <StarVMC/StarAgmlLib/AgMLExtension.h>
+#include <StMessMgr.h>
+
+class AgMLFstVolumeId : public AgMLVolumeId {
+public:
+
+  AgMLFstVolumeId(){}
+
+  virtual int id( int* numbv ) const { 
+
+    static const int mapping[] = { 2, 3, 1 };
+
+    int disk   = numbv[0];
+    int wedge  = numbv[1];
+    assert( numbv[2]>0 && numbv[2]<=3 );
+    int sensor = mapping[ numbv[2]-1 ];
+
+
+    return 1000*disk + 10*wedge + sensor;
+
+  };
+};
+
+
+#endif

StRoot/StGeant4Maker/AgMLHcaVolumeId.h

@@ -0,0 +1,23 @@
+#ifndef __AgMLHcaVolumeId_h__
+#define __AgMLHcaVolumeId_h__
+
+#include <StarVMC/StarAgmlLib/AgMLExtension.h>
+#include <StMessMgr.h>
+
+class AgMLHcaVolumeId : public AgMLVolumeId {
+public:
+
+  AgMLHcaVolumeId(){}
+
+  virtual int id( int* numbv ) const { 
+
+    int mod = numbv[0]; // module
+    int tow = numbv[1]; // tower
+
+    return 1000*mod + tow;
+
+  };
+};
+
+
+#endif