PandaRoot Notebook Python styleΒΆ

This notebook provides code snippets how one can run PandaRoot simulations in a notebook.

[2]:
import ROOT
Welcome to JupyROOT 6.22/08

Here we plainly dropped a master macro and made the syntax pythonic. ROOT carries the PandaRoot libraries that provide a dictionary, which have to be accessed through the ROOT. namespace.

[3]:
nEvents = 10
SimEngine = "TGeant3"
BeamMomentum = 6.231552
parAsciiFile = ROOT.TString("all.par")
# TString inputGenerator =
# EvtGen . "xxxxxxxx.dec"
# DPM    . "dpm_xxxxx"
# FTF    . "ftf_xxxxx"
inputGenerator = ROOT.TString("psi2s_Jpsi2pi_Jpsi_mumu.dec")
inputDir = ROOT.TString(ROOT.gSystem.Getenv("VMCWORKDIR")+"/macro/QA/run3/")
ROOT.gRandom.SetSeed(1234)
#-------------------------------------------------------------------------
# -----   Create the Simulation run manager ------------------------------
fRun = ROOT.PndMasterRunSim()
fRun.SetInput(inputGenerator)
fRun.SetInputDir(inputDir)
fRun.SetName(SimEngine)
fRun.SetParamAsciiFile(parAsciiFile)
fRun.SetNumberOfEvents(nEvents)
fRun.SetBeamMom(BeamMomentum)
# -----  Initialization   ------------------------------------------------
fRun.Setup()
# -----   Geometry   -----------------------------------------------------
fRun.CreateGeometry()
# -----   Event generator   ----------------------------------------------
fRun.SetGenerator()
# -----   Add tasks   ----------------------------------------------------
fRun.AddSimTasks()
# -----   Intialise and run   --------------------------------------------
fRun.Init()
fRun.Run(nEvents)
fRun.Finish()
[INFO] PndSensorNameContFact::PndSensorNameContFact()
[INFO] PndSensorNameContFact::setAllContainers(): Added PndSensorNamePar, Match between GeoManager path and SensorId with context: TestDefaultContext
PndFileNameCreator::fFileName psi2s_jpsi2pi_jpsi_mumu 23
[INFO] Media file used: /mnt/work/PandaRoot/geometry/media_pnd.geo
[INFO] Info in (PndGeoHandling::Instance): Making a new instance using the framework.
 -I container name PndEmcGeoPar
ScitTil fGeoH is loading
[INFO]  PndRich:
[INFO] Using PndEvtGenDirect(pbarpSystem, /mnt/work/PandaRoot/macro/QA/run3/psi2s_Jpsi2pi_Jpsi_mumu.dec, 6.23155) generator
<I> PndEvtGenDirect
<I> Particle: pbarpSystem
<I> decfile: /mnt/work/PandaRoot/macro/QA/run3/psi2s_Jpsi2pi_Jpsi_mumu.dec
<I> Beam-Momentum: 6.23155
<I> Rnd Seed: -1
<I> Rnd Seed changed to 1234
EvtGen:Defining EvtPhotosEngine using photonType = gamma
EvtGen:Setting up PHOTOS.
EvtGen:Using EvtGen random number engine also for Photos++

********************************************************************************

                            =========================
                              PHOTOS, Version:  3.52
                              Released at:  4/11/12
                            =========================

                     Photos QED corrections in Particle Decays

           Monte Carlo Program - by E. Barberio, B. van Eijk and Z. Was
           From version 2.09   - by P. Golonka and Z. Was
           From version 3.00   - by N. Davidson, T. Przedzinski and Z. Was
********************************************************************************

                  Internal (default) input parameters:

                    INTERF= 1 ISEC= 0 ITRE= 0 IEXP= 1 IFTOP= 1 IFW= 1
                    ALPHA_QED= 0x1.de3d3945c2bf2p-8 XPHCUT= 0x1.ad7f29abcaf48p-24

                    Option with interference is active
                    Option with exponentiation is active EPSEXP=0x1.a36e2eb1c432dp-14
                    Emision in t tbar production is active
                    Correction wt in decay of W is active

          WARNING:  /HEPEVT/ is not anymore used.

********************************************************************************
 Message from PHCORK(MODCOR):: initialization
 MODOP=5 -- corrects Energy from mass+flow
EvtGen:Defining EvtPythiaEngine: data tables defined in /cvmfs/fairsoft.gsi.de/debian10/fairsoft/apr21/share/Pythia8/xmldoc
EvtGen:Pythia 6 codes in decay files will be converted to Pythia 8 codes
EvtGen:Using EvtGen random engine for Pythia 8 as well
EvtGen:Creating generic Pythia generator

 *------------------------------------------------------------------------------------*
 |                                                                                    |
 |  *------------------------------------------------------------------------------*  |
 |  |                                                                              |  |
 |  |                                                                              |  |
 |  |   PPP   Y   Y  TTTTT  H   H  III    A      Welcome to the Lund Monte Carlo!  |  |
 |  |   P  P   Y Y     T    H   H   I    A A     This is PYTHIA version 8.303      |  |
 |  |   PPP     Y      T    HHHHH   I   AAAAA    Last date of change:  1 Sep 2020  |  |
 |  |   P       Y      T    H   H   I   A   A                                      |  |
 |  |   P       Y      T    H   H  III  A   A    Now is 10 Aug 2022 at 14:56:07    |  |
 |  |                                                                              |  |
 |  |   Christian Bierlich;  Department of Astronomy and Theoretical Physics,      |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: christian.bierlich@thep.lu.se                                   |  |
 |  |   Nishita Desai;  Department of Theoretical Physics, Tata Institute,         |  |
 |  |      Homi Bhabha Road, Mumbai 400005, India;                                 |  |
 |  |      e-mail: desai@theory.tifr.res.in                                        |  |
 |  |   Leif Gellersen;  Department of Astronomy and Theoretical Physics,          |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: leif.gellersen@thep.lu.se                                       |  |
 |  |   Ilkka Helenius;  Department of Physics, University of Jyvaskyla,           |  |
 |  |      P.O. Box 35, FI-40014 University of Jyvaskyla, Finland;                 |  |
 |  |      e-mail: ilkka.m.helenius@jyu.fi                                         |  |
 |  |   Philip Ilten;  Department of Physics,                                      |  |
 |  |      University of Cincinnati, Cincinnati, OH 45221, USA;                    |  |
 |  |      School of Physics and Astronomy,                                        |  |
 |  |      University of Birmingham, Birmingham, B152 2TT, UK;                     |  |
 |  |      e-mail: philten@cern.ch                                                 |  |
 |  |   Leif Lonnblad;  Department of Astronomy and Theoretical Physics,           |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: leif.lonnblad@thep.lu.se                                        |  |
 |  |   Stephen Mrenna;  Computing Division, Simulations Group,                    |  |
 |  |      Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA;  |  |
 |  |      e-mail: mrenna@fnal.gov                                                 |  |
 |  |   Stefan Prestel;  Department of Astronomy and Theoretical Physics,          |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: stefan.prestel@thep.lu.se                                       |  |
 |  |   Christine O. Rasmussen;  Department of Astronomy and Theoretical Physics,  |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: christine.rasmussen@thep.lu.se                                  |  |
 |  |   Torbjorn Sjostrand;  Department of Astronomy and Theoretical Physics,      |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: torbjorn@thep.lu.se                                             |  |
 |  |   Peter Skands;  School of Physics and Astronomy,                            |  |
 |  |      Monash University, PO Box 27, 3800 Melbourne, Australia;                |  |
 |  |      e-mail: peter.skands@monash.edu                                         |  |
 |  |   Marius Utheim;  Department of Astronomy and Theoretical Physics,           |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: marius.utheim@thep.lu.se                                        |  |
 |  |                                                                              |  |
 |  |   The main program reference is 'An Introduction to PYTHIA 8.2',             |  |
 |  |   T. Sjostrand et al, Comput. Phys. Commun. 191 (2015) 159                   |  |
 |  |   [arXiv:1410.3012 [hep-ph]]                                                 |  |
 |  |                                                                              |  |
 |  |   The main physics reference is the 'PYTHIA 6.4 Physics and Manual',         |  |
 |  |   T. Sjostrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 [hep-ph/0603175]  |  |
 |  |                                                                              |  |
 |  |   An archive of program versions and documentation is found on the web:      |  |
 |  |   http://www.thep.lu.se/Pythia                                               |  |
 |  |                                                                              |  |
 |  |   This program is released under the GNU General Public Licence version 2.   |  |
 |  |   Please respect the MCnet Guidelines for Event Generator Authors and Users. |  |
 |  |                                                                              |  |
 |  |   Disclaimer: this program comes without any guarantees.                     |  |
 |  |   Beware of errors and use common sense when interpreting results.           |  |
 |  |                                                                              |  |
 |  |   Copyright (C) 2020 Torbjorn Sjostrand                                      |  |
 |  |                                                                              |  |
 |  |                                                                              |  |
 |  *------------------------------------------------------------------------------*  |
 |                                                                                    |
 *------------------------------------------------------------------------------------*

EvtGen:Creating alias Pythia generator

 *------------------------------------------------------------------------------------*
 |                                                                                    |
 |  *------------------------------------------------------------------------------*  |
 |  |                                                                              |  |
 |  |                                                                              |  |
 |  |   PPP   Y   Y  TTTTT  H   H  III    A      Welcome to the Lund Monte Carlo!  |  |
 |  |   P  P   Y Y     T    H   H   I    A A     This is PYTHIA version 8.303      |  |
 |  |   PPP     Y      T    HHHHH   I   AAAAA    Last date of change:  1 Sep 2020  |  |
 |  |   P       Y      T    H   H   I   A   A                                      |  |
 |  |   P       Y      T    H   H  III  A   A    Now is 10 Aug 2022 at 14:56:07    |  |
 |  |                                                                              |  |
 |  |   Christian Bierlich;  Department of Astronomy and Theoretical Physics,      |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: christian.bierlich@thep.lu.se                                   |  |
 |  |   Nishita Desai;  Department of Theoretical Physics, Tata Institute,         |  |
 |  |      Homi Bhabha Road, Mumbai 400005, India;                                 |  |
 |  |      e-mail: desai@theory.tifr.res.in                                        |  |
 |  |   Leif Gellersen;  Department of Astronomy and Theoretical Physics,          |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: leif.gellersen@thep.lu.se                                       |  |
 |  |   Ilkka Helenius;  Department of Physics, University of Jyvaskyla,           |  |
 |  |      P.O. Box 35, FI-40014 University of Jyvaskyla, Finland;                 |  |
 |  |      e-mail: ilkka.m.helenius@jyu.fi                                         |  |
 |  |   Philip Ilten;  Department of Physics,                                      |  |
 |  |      University of Cincinnati, Cincinnati, OH 45221, USA;                    |  |
 |  |      School of Physics and Astronomy,                                        |  |
 |  |      University of Birmingham, Birmingham, B152 2TT, UK;                     |  |
 |  |      e-mail: philten@cern.ch                                                 |  |
 |  |   Leif Lonnblad;  Department of Astronomy and Theoretical Physics,           |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: leif.lonnblad@thep.lu.se                                        |  |
 |  |   Stephen Mrenna;  Computing Division, Simulations Group,                    |  |
 |  |      Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA;  |  |
 |  |      e-mail: mrenna@fnal.gov                                                 |  |
 |  |   Stefan Prestel;  Department of Astronomy and Theoretical Physics,          |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: stefan.prestel@thep.lu.se                                       |  |
 |  |   Christine O. Rasmussen;  Department of Astronomy and Theoretical Physics,  |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: christine.rasmussen@thep.lu.se                                  |  |
 |  |   Torbjorn Sjostrand;  Department of Astronomy and Theoretical Physics,      |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: torbjorn@thep.lu.se                                             |  |
 |  |   Peter Skands;  School of Physics and Astronomy,                            |  |
 |  |      Monash University, PO Box 27, 3800 Melbourne, Australia;                |  |
 |  |      e-mail: peter.skands@monash.edu                                         |  |
 |  |   Marius Utheim;  Department of Astronomy and Theoretical Physics,           |  |
 |  |      Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden;                |  |
 |  |      e-mail: marius.utheim@thep.lu.se                                        |  |
 |  |                                                                              |  |
 |  |   The main program reference is 'An Introduction to PYTHIA 8.2',             |  |
 |  |   T. Sjostrand et al, Comput. Phys. Commun. 191 (2015) 159                   |  |
 |  |   [arXiv:1410.3012 [hep-ph]]                                                 |  |
 |  |                                                                              |  |
 |  |   The main physics reference is the 'PYTHIA 6.4 Physics and Manual',         |  |
 |  |   T. Sjostrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 [hep-ph/0603175]  |  |
 |  |                                                                              |  |
 |  |   An archive of program versions and documentation is found on the web:      |  |
 |  |   http://www.thep.lu.se/Pythia                                               |  |
 |  |                                                                              |  |
 |  |   This program is released under the GNU General Public Licence version 2.   |  |
 |  |   Please respect the MCnet Guidelines for Event Generator Authors and Users. |  |
 |  |                                                                              |  |
 |  |   Disclaimer: this program comes without any guarantees.                     |  |
 |  |   Beware of errors and use common sense when interpreting results.           |  |
 |  |                                                                              |  |
 |  |   Copyright (C) 2020 Torbjorn Sjostrand                                      |  |
 |  |                                                                              |  |
 |  |                                                                              |  |
 |  *------------------------------------------------------------------------------*  |
 |                                                                                    |
 *------------------------------------------------------------------------------------*

EvtGen:Initializing EvtGen
EvtGen:Storing known decay models
EvtGen:Defining EvtBcToNPi model: Bc -> V + npi and Bc -> P + npi decays
from A.V. Berezhnoy, A.K. Likhoded, A.V. Luchinsky: Phys.Rev.D 82, 014012 (2010) and arXiV:1104.0808.
EvtGen:Main decay file name  :/mnt/work/PandaRoot/pgenerators/EvtGen/EvtGen/Private/DECAY.DEC
EvtGen:PDT table file name   :/mnt/work/PandaRoot/pgenerators/EvtGen/EvtGen/Private/evt.pdl
EvtGen:In readDecayFile, reading:/mnt/work/PandaRoot/pgenerators/EvtGen/EvtGen/Private/DECAY.DEC
EvtGen:Warning, sum of branching fractions for omega is 0.994
EvtGen:rescaled to one!
EvtGen:Warning, sum of branching fractions for phi is 0.999
EvtGen:rescaled to one!
EvtGen:Warning, sum of branching fractions for f_2 is 0.999
EvtGen:rescaled to one!
EvtGen:Warning, sum of branching fractions for chi_c1 is 0.996
EvtGen:rescaled to one!
EvtGen:Warning, sum of branching fractions for Upsilon(2S) is 0.992
EvtGen:rescaled to one!
EvtGen:Warning, sum of branching fractions for Upsilon(3S) is 0.999
EvtGen:rescaled to one!
EvtGen:Mixing type integer set to 1
EvtGen:Done initializing EvtGen
EvtGen:In readDecayFile, reading:/mnt/work/PandaRoot/macro/QA/run3/psi2s_Jpsi2pi_Jpsi_mumu.dec
EvtGen:As requested, PHOTOS will be turned off.
EvtGen:Redefined decay of J/psi
EvtGen:Given allowed decays, resetting minMass pbarpSystem 3.100 to 3.375

############# Generating with following conditions:

incident 4-mom : (7.240, 0, 0, 6.232), m = 3.686

######################


[INFO] Using no Vertex smearing
[INFO] FairRootFileSink initialized.
[INFO]  - pnd_0
[INFO]     - ./psi2s_jpsi2pi_jpsi_mumu_sim.root
[INFO] ==============  FairRunSim: Initialising simulation run ==============
[INFO] FairGeoMedia: Read media
[INFO] PndSensorNameContFact::createContainer()
[INFO] Created PndSensorNamePar - PndSensorNamePar -

*************************************************************
     initialisation for run id 1660136168
*************************************************************
-I- FairRunTimeDB::InitContainer() PndEmcGeoPar
[ERROR] init() PndEmcGeoPar not initialized
-I- FairRunTimeDB::InitContainer() PndSensorNamePar
[ERROR] init() PndSensorNamePar not initialized
 -I container name PndEmcDigiPar
 -I container name PndEmcDigiNonuniformityPar

*************************************************************
     initialisation for run id 1660136168
*************************************************************
-I- FairRunTimeDB::InitContainer() PndEmcGeoPar
[ERROR] init() PndEmcGeoPar not initialized
-I- FairRunTimeDB::InitContainer() PndSensorNamePar
[ERROR] init() PndSensorNamePar not initialized
-I- FairRunTimeDB::InitContainer() PndEmcDigiPar
PndEmcDigiPar initialized from Ascii file
-I- FairRunTimeDB::InitContainer() PndEmcDigiNonuniformityPar
[ERROR] init() PndEmcDigiNonuniformityPar not initialized
[INFO] PndFieldMap: Reading field map from ROOT file  /mnt/work/PandaRoot/input/fieldmaps/TransMap.0890.root
[INFO] PndFieldMap: Reading field map from ROOT file  /mnt/work/PandaRoot/input/fieldmaps/DipoleMap1.0890.root
[INFO] PndFieldMap: Reading field map from ROOT file  /mnt/work/PandaRoot/input/fieldmaps/DipoleMap2.0890.root
[INFO] PndFieldMap: Reading field map from ROOT file  /mnt/work/PandaRoot/input/fieldmaps/SolenoidMap1.root
[INFO] PndFieldMap: Reading field map from ROOT file  /mnt/work/PandaRoot/input/fieldmaps/SolenoidMap2.root
[INFO] PndFieldMap: Reading field map from ROOT file  /mnt/work/PandaRoot/input/fieldmaps/SolenoidMap3.root
[INFO] PndFieldMap: Reading field map from ROOT file  /mnt/work/PandaRoot/input/fieldmaps/SolenoidMap4.root
[WARN] FairGenericVMCConfig::Setup() Using gConfig.C macro DEPRACATED.
[WARN] Check FairRoot/examples/common/gconfig/ for current YAML implementation.
[INFO] -------------- Standard Config is called ------------------------------------
[INFO] ---VMCWORKDIR path for Configuration (g3libs.C) is used: /mnt/work/PandaRoot/gconfig/

=============================================================
 Virtual Monte Carlo Library
 Version 1.0.p3 ( 21 July 2020 )
=============================================================

=============================================================
 Geant3 Virtual Monte Carlo
 Version 3.8 ( 17 February 2021 )
 WWW : https://vmc-project.github.io/
=============================================================


 MZSTOR.  ZEBRA table base TAB(0) in /MZCC/ at adr   281008583    10BFD9C7 HEX

 MZSTOR.  Initialize Store  0  in /GCBANK/
          with Store/Table at absolute adrs   275562925   281008583
                                        HEX    106CC1AD    10BFD9C7
                                        HEX    FFACE5D2           0
                              relative adrs    -5446190           0
          with     1 Str. in     2 Links in   5300 Low words in 4999970 words.
          This store has a fence of   16 words.

 MZLOGL.  Set Log Level 0 for store  0
1*****  GEANT Version  3.21/11 Released on 100298
0*****  Correction Cradle Version  0.1100

 MZDIV.   Initialize Division  Constant  in Store  0
          NW/NWMAX=   20004000000,  MODE/KIND=  1  2
          Division 20 initialized.

 MZLINK.  Initialize Link Area  /GCLINK/  for Store  0 NL/NS=    20    20

 MZLINK.  Initialize Link Area  /GCSLNK/  for Store  0 NL/NS=   100   100
PndTargetGenerator::ReadDensityFile():
PndTargetGenerator::Init():
<I> PndPipe - Using geometry /mnt/work/PandaRoot/geometry/beampipe_201309.root
 -I container name PndGeoSttPar
[INFO]  STT total number of tubes: 4542
[INFO]  PndStt::ConstructGeometry : STT total number of tubes: 0
fwendcap & bwendcap flags == 1 / 1

 ====== EMC 2)::  ConstructROOTGeometry() ===
 =============================================
File name = /mnt/work/PandaRoot/geometry/emc_module12_2018v1.root

 ====== EMC::  ConstructRootGeometry() m3a ===
 =============================================
File name = /mnt/work/PandaRoot/geometry/emc_module3_2012_new.root

 ====== EMC::  ConstructRootGeometry() m4a ===
 =============================================
fgeoName3:: /mnt/work/PandaRoot/geometry/emc_module4_2017.root
File name Bw1= /mnt/work/PandaRoot/geometry/emc_module4_2017.root
PndEmc::ConstructRootGeomMod4: halfLength = 2.190000e+01, fullLength = 4.380000e+01, z-pos = -7.790000e+01, distTargetBWEC = 5.600000e+01

 ====== EMC::  ConstructRootGeometry() m5a ===
 =============================================
fgeoName4:: /mnt/work/PandaRoot/geometry/emc_module5_fsc.root
File name Fsc= /mnt/work/PandaRoot/geometry/emc_module5_fsc.root
[INFO]  =======  DRC::  ConstructGeometry()  ========
 -I container name PndGeoFtsPar
[INFO] TGeometry will be imported to VMC

0*** GSTMED *** Warning, medium =     4, value of EPSIL= 0.000E+00 reset to 1 micron
[INFO] PndDrc: constructing optical geometry
[INFO] PndDrc: geometry construction finished
 ====================================================
 =======  Rich::  ConstructOpticalGeometry()  ========
GetGeometryFileName() = /mnt/work/PandaRoot/geometry/rich_v313.root 313
fGeoVersion = 313
nOpt.size() = 3
 =======  RICH::ConstructOpGeometry -> Finished! ======
[INFO] FairMCApplication::InitGeometry: 0
[INFO]  Initializing PndStt()
 -I- Initializing PndSdsDetector()
[INFO] Branch: MVDPoint not registered!
[WARN]  PndSdsDetector: New branch MVDPoint created!
 -I- Initializing PndGemDetector()

************************************************************
PndEmc::SetSpecialPhysicsCuts():
   using special physics cuts ...
************************************************************
 -I- Initializing PndSciT()
 Fairdetector geht
 -I- Initialized PndSciT()
[INFO] PndDrc: intialization started
[INFO] PndDrc: Cherenkov Propagation is OFF
[INFO] PndDrc: intialization successfull
 -I- Initializing PndFts()
 -I- Initializing PndFtof()
there is gGeoManager
list of sensitives has 1 entries
 -I- PndRich: Switching OFF Cherenkov Propagation
GetGeometryFileName() = /mnt/work/PandaRoot/geometry/rich_v313.root 313
Mid ref. index = 1.050
[INFO] Initialize Tasks--------------------------
[INFO] PndSensorNamePar::FillMap() - There are 2809 entries.
[INFO] PndSensorNamePar::FillMap() done
 -I- PndEmcHitProducer INITIALIZATION *********************
[INFO]  PndEmcHitProducer: Using nonuniform light output
HitProducer has EnergyHitThreshold of 0.000001 GeV and Use_nonuniformity 1
[INFO]  PndEmcHitProducer: Initialization successful
[INFO] Simulation RunID: 1660136168

 Calculating cross section tables, see gphysi.dat for more information

 Cross section calculation concluded successfully
[INFO] Monte Carlo Engine Initialisation with: TGeant3TGeo
[INFO] FairGenericVMCConfig::SetupPostInit() OFF.
 **** GTRIGI: IEVENT=      1 IDEVT=      1 Random Seeds =       1234           0
 **** GTRIGI: IEVENT=      2 IDEVT=      2 Random Seeds =       1234           0
 **** GTRIGI: IEVENT=      3 IDEVT=      3 Random Seeds =       1234           0
 **** GTRIGI: IEVENT=      4 IDEVT=      4 Random Seeds =       1234           0
 **** GTRIGI: IEVENT=      5 IDEVT=      5 Random Seeds =       1234           0
 **** GTRIGI: IEVENT=      6 IDEVT=      6 Random Seeds =       1234           0
 **** GTRIGI: IEVENT=      7 IDEVT=      7 Random Seeds =       1234           0
 **** GTRIGI: IEVENT=      8 IDEVT=      8 Random Seeds =       1234           0
 **** GTRIGI: IEVENT=      9 IDEVT=      9 Random Seeds =       1234           0
 **** GTRIGI: IEVENT=     10 IDEVT=     10 Random Seeds =       1234           0
=========================================================
PndEmcHitProducer::FinishTask
*********************************************************
Read points # 114651
Produc hits# 893, threshold# 0.000
Hits above threshhod#263
*********************************************************
[INFO] ***  PndEmcGeoPar written to ROOT file   version: 1
[INFO] ***  PndSensorNamePar written to ROOT file   version: 1
[INFO] ***  PndEmcDigiPar written to ROOT file   version: 1
[INFO] ***  PndEmcDigiNonuniformityPar written to ROOT file   version: 1
[INFO] ***  FairBaseParSet written to ROOT file   version: 1
[INFO] ***  FairGeoParSet written to ROOT file   version: 1
[INFO] ***  PndMultiFieldPar written to ROOT file   version: 1
[INFO] ***  PndGeoPassivePar written to ROOT file   version: 1
[INFO] ***  PndGeoSttPar written to ROOT file   version: 1
[INFO] ***  PndGeoFtsPar written to ROOT file   version: 1
[INFO] ***  PndGeoSciTPar written to ROOT file   version: 1
[INFO] ***  PndGeoFtofPar written to ROOT file   version: 1
[INFO] ***  PndRichGeoPar written to ROOT file   version: 1
PndMasterRunSim::Finish(): Tasks that ran just now:
file is not open
Opening file: ./psi2s_jpsi2pi_jpsi_mumu_sim.root


Generated Events = 10
 - PndGeoHandling (PndGeoHandling)
 - PndEventCounterTask (Event Counter)
 - PndMasterSimTask (Master Sim Task)
   - PndEmcHitProducer (Ideal EMC hit Producer)

<DartMeasurement name="MaxMemory" type="numeric/double">831.688</DartMeasurement>
<DartMeasurement name="CpuLoad" type="numeric/double">0.987</DartMeasurement>

Output file is          ./psi2s_jpsi2pi_jpsi_mumu_sim.root
Parameter ROOT file is  ./psi2s_jpsi2pi_jpsi_mumu_par.root
Parameter ASCII file is all.par
Real time 16.156 s, CPU time 15.950s
CPU usage 98.724%
Max Memory 831.688 MB
Macro finished successfully.
 Message from PHCORK(MODCOR):: initialization
 MODOP=5 -- corrects Energy from mass+flow
 Message from PHOTOS: IPHQRK(MODCOR):: (re)initialization
 MODOP=1 -- blocks emission from light quarks: DEFAULT
 Message from PHOTOS: IPHEKL(MODCOR):: (re)initialization
 MODOP=2 -- blocks emission in pi0 to gamma e+e-: DEFAULT
 MODOP=2 -- blocks emission in Kl  to gamma e+e-: DEFAULT
Info in <TGeoManager::TGeoManager>: Geometry FAIRGeom, FAIR geometry created
Error in <FairRuntimeDb::initContainers()>: Error occured during initialization
Error in <FairRuntimeDb::initContainers()>: Error occured during initialization
Info in <TGeoManager::SetTopVolume>: Top volume is cave. Master volume is cave
Info in <TGeoNavigator::BuildCache>: --- Maximum geometry depth set to 100
Info in <TGeoManager::CheckGeometry>: Fixing runtime shapes...
Info in <TGeoManager::CheckGeometry>: ...Nothing to fix
Info in <TGeoManager::CloseGeometry>: Counting nodes...
Info in <TGeoManager::Voxelize>: Voxelizing...
Info in <TGeoManager::CloseGeometry>: Building cache...
Info in <TGeoManager::CountLevels>: max level = 12, max placements = 4550
Info in <TGeoManager::CloseGeometry>: 2029436 nodes/ 5737 volume UID's in FAIR geometry
Info in <TGeoManager::CloseGeometry>: ----------------modeler ready----------------
Warning in <TGeant3TGeo::DefineOpSurface>: Called for surface LensSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::DefineOpSurface>: Called for surface MirrSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::DefineOpSurface>: Called for surface EVSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::DefineOpSurface>: Called for surface BlackSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::DefineOpSurface>: Called for surface DrcChamferSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetMaterialProperty>: Called for material surface  BlackSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetMaterialProperty>: Called for material surface  MirrSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface BarMirrorSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface DrcBar_ChamferSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface DrcBar_ChamferSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface DrcBar_ChamferSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface Lens1AirSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface Lens2AirSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface Lens1AirSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface Lens2AirSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface Lens1AirSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface Lens2AirSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetSkinSurface>: Called for skin surface AirMirrorSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::DefineOpSurface>: Called for surface RichMirrSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetMaterialProperty>: Called for material surface  RichMirrSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetBorderSurface>: Called for border surface BarRichMirrorSurface. Not applicable in Geant3 - setting is ignored.
Warning in <TGeant3TGeo::SetSkinSurface>: Called for skin surface RichAirMirrorSurface. Not applicable in Geant3 - setting is ignored.

Plotting with ROOT is still working. To show a canvas use its Draw() function.

[9]:
myfile = ROOT.TFile("psi2s_jpsi2pi_jpsi_mumu_sim.root")
mytree = myfile.Get("pndsim")
can = ROOT.TCanvas()
mytree.Draw("MCTrack.fStartZ","")
can.Draw()
../../_images/Jupyter_notebooks_PandaRootPySim_5_0.png
[ ]: