{ "cells": [ { "cell_type": "markdown", "id": "2126cfd4", "metadata": {}, "source": [ "# PandaRoot Notebook Python style\n", "\n", "This notebook provides code snippets how one can run PandaRoot simulations in a notebook." ] }, { "cell_type": "code", "execution_count": 2, "id": "3bdfe179", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Welcome to JupyROOT 6.22/08\n" ] } ], "source": [ "import ROOT" ] }, { "cell_type": "markdown", "id": "283ea6ab", "metadata": {}, "source": [ "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." ] }, { "cell_type": "code", "execution_count": 3, "id": "4ec2fbcb", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[INFO] PndSensorNameContFact::PndSensorNameContFact()\n", "[INFO] PndSensorNameContFact::setAllContainers(): Added PndSensorNamePar, Match between GeoManager path and SensorId with context: TestDefaultContext\n", "PndFileNameCreator::fFileName psi2s_jpsi2pi_jpsi_mumu 23\n", "[INFO] Media file used: /mnt/work/PandaRoot/geometry/media_pnd.geo\n", "[INFO] Info in (PndGeoHandling::Instance): Making a new instance using the framework.\n", " -I container name PndEmcGeoPar\n", "ScitTil fGeoH is loading\n", "[INFO] PndRich:\n", "[INFO] Using PndEvtGenDirect(pbarpSystem, /mnt/work/PandaRoot/macro/QA/run3/psi2s_Jpsi2pi_Jpsi_mumu.dec, 6.23155) generator\n", " PndEvtGenDirect\n", " Particle: pbarpSystem\n", " decfile: /mnt/work/PandaRoot/macro/QA/run3/psi2s_Jpsi2pi_Jpsi_mumu.dec\n", " Beam-Momentum: 6.23155\n", " Rnd Seed: -1\n", " Rnd Seed changed to 1234\n", "EvtGen:Defining EvtPhotosEngine using photonType = gamma\n", "EvtGen:Setting up PHOTOS.\n", "EvtGen:Using EvtGen random number engine also for Photos++\n", "\n", "********************************************************************************\n", "\n", " =========================\n", " PHOTOS, Version: 3.52\n", " Released at: 4/11/12\n", " =========================\n", "\n", " Photos QED corrections in Particle Decays\n", "\n", " Monte Carlo Program - by E. Barberio, B. van Eijk and Z. Was\n", " From version 2.09 - by P. Golonka and Z. Was\n", " From version 3.00 - by N. Davidson, T. Przedzinski and Z. Was\n", "********************************************************************************\n", "\n", " Internal (default) input parameters: \n", "\n", " INTERF= 1 ISEC= 0 ITRE= 0 IEXP= 1 IFTOP= 1 IFW= 1\n", " ALPHA_QED= 0x1.de3d3945c2bf2p-8 XPHCUT= 0x1.ad7f29abcaf48p-24\n", "\n", " Option with interference is active\n", " Option with exponentiation is active EPSEXP=0x1.a36e2eb1c432dp-14\n", " Emision in t tbar production is active\n", " Correction wt in decay of W is active\n", "\n", " WARNING: /HEPEVT/ is not anymore used.\n", "\n", "********************************************************************************\n", " Message from PHCORK(MODCOR):: initialization\n", " MODOP=5 -- corrects Energy from mass+flow\n", "EvtGen:Defining EvtPythiaEngine: data tables defined in /cvmfs/fairsoft.gsi.de/debian10/fairsoft/apr21/share/Pythia8/xmldoc\n", "EvtGen:Pythia 6 codes in decay files will be converted to Pythia 8 codes\n", "EvtGen:Using EvtGen random engine for Pythia 8 as well\n", "EvtGen:Creating generic Pythia generator\n", "\n", " *------------------------------------------------------------------------------------* \n", " | | \n", " | *------------------------------------------------------------------------------* | \n", " | | | | \n", " | | | | \n", " | | PPP Y Y TTTTT H H III A Welcome to the Lund Monte Carlo! | | \n", " | | P P Y Y T H H I A A This is PYTHIA version 8.303 | | \n", " | | PPP Y T HHHHH I AAAAA Last date of change: 1 Sep 2020 | | \n", " | | P Y T H H I A A | | \n", " | | P Y T H H III A A Now is 10 Aug 2022 at 14:56:07 | | \n", " | | | | \n", " | | Christian Bierlich; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: christian.bierlich@thep.lu.se | | \n", " | | Nishita Desai; Department of Theoretical Physics, Tata Institute, | | \n", " | | Homi Bhabha Road, Mumbai 400005, India; | | \n", " | | e-mail: desai@theory.tifr.res.in | | \n", " | | Leif Gellersen; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: leif.gellersen@thep.lu.se | | \n", " | | Ilkka Helenius; Department of Physics, University of Jyvaskyla, | | \n", " | | P.O. Box 35, FI-40014 University of Jyvaskyla, Finland; | | \n", " | | e-mail: ilkka.m.helenius@jyu.fi | | \n", " | | Philip Ilten; Department of Physics, | | \n", " | | University of Cincinnati, Cincinnati, OH 45221, USA; | | \n", " | | School of Physics and Astronomy, | | \n", " | | University of Birmingham, Birmingham, B152 2TT, UK; | | \n", " | | e-mail: philten@cern.ch | | \n", " | | Leif Lonnblad; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: leif.lonnblad@thep.lu.se | | \n", " | | Stephen Mrenna; Computing Division, Simulations Group, | | \n", " | | Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA; | | \n", " | | e-mail: mrenna@fnal.gov | | \n", " | | Stefan Prestel; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: stefan.prestel@thep.lu.se | | \n", " | | Christine O. Rasmussen; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: christine.rasmussen@thep.lu.se | | \n", " | | Torbjorn Sjostrand; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: torbjorn@thep.lu.se | | \n", " | | Peter Skands; School of Physics and Astronomy, | | \n", " | | Monash University, PO Box 27, 3800 Melbourne, Australia; | | \n", " | | e-mail: peter.skands@monash.edu | | \n", " | | Marius Utheim; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: marius.utheim@thep.lu.se | | \n", " | | | | \n", " | | The main program reference is 'An Introduction to PYTHIA 8.2', | | \n", " | | T. Sjostrand et al, Comput. Phys. Commun. 191 (2015) 159 | | \n", " | | [arXiv:1410.3012 [hep-ph]] | | \n", " | | | | \n", " | | The main physics reference is the 'PYTHIA 6.4 Physics and Manual', | | \n", " | | T. Sjostrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 [hep-ph/0603175] | | \n", " | | | | \n", " | | An archive of program versions and documentation is found on the web: | | \n", " | | http://www.thep.lu.se/Pythia | | \n", " | | | | \n", " | | This program is released under the GNU General Public Licence version 2. | | \n", " | | Please respect the MCnet Guidelines for Event Generator Authors and Users. | | \n", " | | | | \n", " | | Disclaimer: this program comes without any guarantees. | | \n", " | | Beware of errors and use common sense when interpreting results. | | \n", " | | | | \n", " | | Copyright (C) 2020 Torbjorn Sjostrand | | \n", " | | | | \n", " | | | | \n", " | *------------------------------------------------------------------------------* | \n", " | | \n", " *------------------------------------------------------------------------------------* \n", "\n", "EvtGen:Creating alias Pythia generator\n", "\n", " *------------------------------------------------------------------------------------* \n", " | | \n", " | *------------------------------------------------------------------------------* | \n", " | | | | \n", " | | | | \n", " | | PPP Y Y TTTTT H H III A Welcome to the Lund Monte Carlo! | | \n", " | | P P Y Y T H H I A A This is PYTHIA version 8.303 | | \n", " | | PPP Y T HHHHH I AAAAA Last date of change: 1 Sep 2020 | | \n", " | | P Y T H H I A A | | \n", " | | P Y T H H III A A Now is 10 Aug 2022 at 14:56:07 | | \n", " | | | | \n", " | | Christian Bierlich; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: christian.bierlich@thep.lu.se | | \n", " | | Nishita Desai; Department of Theoretical Physics, Tata Institute, | | \n", " | | Homi Bhabha Road, Mumbai 400005, India; | | \n", " | | e-mail: desai@theory.tifr.res.in | | \n", " | | Leif Gellersen; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: leif.gellersen@thep.lu.se | | \n", " | | Ilkka Helenius; Department of Physics, University of Jyvaskyla, | | \n", " | | P.O. Box 35, FI-40014 University of Jyvaskyla, Finland; | | \n", " | | e-mail: ilkka.m.helenius@jyu.fi | | \n", " | | Philip Ilten; Department of Physics, | | \n", " | | University of Cincinnati, Cincinnati, OH 45221, USA; | | \n", " | | School of Physics and Astronomy, | | \n", " | | University of Birmingham, Birmingham, B152 2TT, UK; | | \n", " | | e-mail: philten@cern.ch | | \n", " | | Leif Lonnblad; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: leif.lonnblad@thep.lu.se | | \n", " | | Stephen Mrenna; Computing Division, Simulations Group, | | \n", " | | Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA; | | \n", " | | e-mail: mrenna@fnal.gov | | \n", " | | Stefan Prestel; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: stefan.prestel@thep.lu.se | | \n", " | | Christine O. Rasmussen; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: christine.rasmussen@thep.lu.se | | \n", " | | Torbjorn Sjostrand; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: torbjorn@thep.lu.se | | \n", " | | Peter Skands; School of Physics and Astronomy, | | \n", " | | Monash University, PO Box 27, 3800 Melbourne, Australia; | | \n", " | | e-mail: peter.skands@monash.edu | | \n", " | | Marius Utheim; Department of Astronomy and Theoretical Physics, | | \n", " | | Lund University, Solvegatan 14A, SE-223 62 Lund, Sweden; | | \n", " | | e-mail: marius.utheim@thep.lu.se | | \n", " | | | | \n", " | | The main program reference is 'An Introduction to PYTHIA 8.2', | | \n", " | | T. Sjostrand et al, Comput. Phys. Commun. 191 (2015) 159 | | \n", " | | [arXiv:1410.3012 [hep-ph]] | | \n", " | | | | \n", " | | The main physics reference is the 'PYTHIA 6.4 Physics and Manual', | | \n", " | | T. Sjostrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 [hep-ph/0603175] | | \n", " | | | | \n", " | | An archive of program versions and documentation is found on the web: | | \n", " | | http://www.thep.lu.se/Pythia | | \n", " | | | | \n", " | | This program is released under the GNU General Public Licence version 2. | | \n", " | | Please respect the MCnet Guidelines for Event Generator Authors and Users. | | \n", " | | | | \n", " | | Disclaimer: this program comes without any guarantees. | | \n", " | | Beware of errors and use common sense when interpreting results. | | \n", " | | | | \n", " | | Copyright (C) 2020 Torbjorn Sjostrand | | \n", " | | | | \n", " | | | | \n", " | *------------------------------------------------------------------------------* | \n", " | | \n", " *------------------------------------------------------------------------------------* \n", "\n", "EvtGen:Initializing EvtGen\n", "EvtGen:Storing known decay models\n", "EvtGen:Defining EvtBcToNPi model: Bc -> V + npi and Bc -> P + npi decays\n", "from A.V. Berezhnoy, A.K. Likhoded, A.V. Luchinsky: Phys.Rev.D 82, 014012 (2010) and arXiV:1104.0808.\n", "EvtGen:Main decay file name :/mnt/work/PandaRoot/pgenerators/EvtGen/EvtGen/Private/DECAY.DEC\n", "EvtGen:PDT table file name :/mnt/work/PandaRoot/pgenerators/EvtGen/EvtGen/Private/evt.pdl\n", "EvtGen:In readDecayFile, reading:/mnt/work/PandaRoot/pgenerators/EvtGen/EvtGen/Private/DECAY.DEC\n", "EvtGen:Warning, sum of branching fractions for omega is 0.994\n", "EvtGen:rescaled to one! \n", "EvtGen:Warning, sum of branching fractions for phi is 0.999\n", "EvtGen:rescaled to one! \n", "EvtGen:Warning, sum of branching fractions for f_2 is 0.999\n", "EvtGen:rescaled to one! \n", "EvtGen:Warning, sum of branching fractions for chi_c1 is 0.996\n", "EvtGen:rescaled to one! \n", "EvtGen:Warning, sum of branching fractions for Upsilon(2S) is 0.992\n", "EvtGen:rescaled to one! \n", "EvtGen:Warning, sum of branching fractions for Upsilon(3S) is 0.999\n", "EvtGen:rescaled to one! \n", "EvtGen:Mixing type integer set to 1\n", "EvtGen:Done initializing EvtGen\n", "EvtGen:In readDecayFile, reading:/mnt/work/PandaRoot/macro/QA/run3/psi2s_Jpsi2pi_Jpsi_mumu.dec\n", "EvtGen:As requested, PHOTOS will be turned off.\n", "EvtGen:Redefined decay of J/psi\n", "EvtGen:Given allowed decays, resetting minMass pbarpSystem 3.100 to 3.375\n", "\n", "############# Generating with following conditions:\n", "\n", "incident 4-mom : (7.240, 0, 0, 6.232), m = 3.686\n", "\n", "######################\n", "\n", "\n", "[INFO] Using no Vertex smearing\n", "[INFO] FairRootFileSink initialized.\n", "[INFO] - pnd_0\n", "[INFO] - ./psi2s_jpsi2pi_jpsi_mumu_sim.root\n", "[INFO] ============== FairRunSim: Initialising simulation run ==============\n", "[INFO] FairGeoMedia: Read media \n", "[INFO] PndSensorNameContFact::createContainer()\n", "[INFO] Created PndSensorNamePar - PndSensorNamePar - \n", "\n", "************************************************************* \n", " initialisation for run id 1660136168\n", "************************************************************* \n", "-I- FairRunTimeDB::InitContainer() PndEmcGeoPar\n", "[ERROR] init() PndEmcGeoPar not initialized\n", "-I- FairRunTimeDB::InitContainer() PndSensorNamePar\n", "[ERROR] init() PndSensorNamePar not initialized\n", " -I container name PndEmcDigiPar\n", " -I container name PndEmcDigiNonuniformityPar\n", "\n", "************************************************************* \n", " initialisation for run id 1660136168\n", "************************************************************* \n", "-I- FairRunTimeDB::InitContainer() PndEmcGeoPar\n", "[ERROR] init() PndEmcGeoPar not initialized\n", "-I- FairRunTimeDB::InitContainer() PndSensorNamePar\n", "[ERROR] init() PndSensorNamePar not initialized\n", "-I- FairRunTimeDB::InitContainer() PndEmcDigiPar\n", "PndEmcDigiPar initialized from Ascii file\n", "-I- FairRunTimeDB::InitContainer() PndEmcDigiNonuniformityPar\n", "[ERROR] init() PndEmcDigiNonuniformityPar not initialized\n", "[INFO] PndFieldMap: Reading field map from ROOT file /mnt/work/PandaRoot/input/fieldmaps/TransMap.0890.root\n", "[INFO] PndFieldMap: Reading field map from ROOT file /mnt/work/PandaRoot/input/fieldmaps/DipoleMap1.0890.root\n", "[INFO] PndFieldMap: Reading field map from ROOT file /mnt/work/PandaRoot/input/fieldmaps/DipoleMap2.0890.root\n", "[INFO] PndFieldMap: Reading field map from ROOT file /mnt/work/PandaRoot/input/fieldmaps/SolenoidMap1.root\n", "[INFO] PndFieldMap: Reading field map from ROOT file /mnt/work/PandaRoot/input/fieldmaps/SolenoidMap2.root\n", "[INFO] PndFieldMap: Reading field map from ROOT file /mnt/work/PandaRoot/input/fieldmaps/SolenoidMap3.root\n", "[INFO] PndFieldMap: Reading field map from ROOT file /mnt/work/PandaRoot/input/fieldmaps/SolenoidMap4.root\n", "[WARN] FairGenericVMCConfig::Setup() Using gConfig.C macro DEPRACATED.\n", "[WARN] Check FairRoot/examples/common/gconfig/ for current YAML implementation.\n", "[INFO] -------------- Standard Config is called ------------------------------------\n", "[INFO] ---VMCWORKDIR path for Configuration (g3libs.C) is used: /mnt/work/PandaRoot/gconfig/\n", "\n", "=============================================================\n", " Virtual Monte Carlo Library\n", " Version 1.0.p3 ( 21 July 2020 )\n", "=============================================================\n", "\n", "=============================================================\n", " Geant3 Virtual Monte Carlo \n", " Version 3.8 ( 17 February 2021 )\n", " WWW : https://vmc-project.github.io/\n", "=============================================================\n", "\n", "\n", " MZSTOR. ZEBRA table base TAB(0) in /MZCC/ at adr 281008583 10BFD9C7 HEX\n", "\n", " MZSTOR. Initialize Store 0 in /GCBANK/\n", " with Store/Table at absolute adrs 275562925 281008583\n", " HEX 106CC1AD 10BFD9C7\n", " HEX FFACE5D2 0\n", " relative adrs -5446190 0\n", " with 1 Str. in 2 Links in 5300 Low words in 4999970 words.\n", " This store has a fence of 16 words.\n", "\n", " MZLOGL. Set Log Level 0 for store 0\n", "1***** GEANT Version 3.21/11 Released on 100298\n", "0***** Correction Cradle Version 0.1100\n", "\n", " MZDIV. Initialize Division Constant in Store 0\n", " NW/NWMAX= 20004000000, MODE/KIND= 1 2\n", " Division 20 initialized.\n", "\n", " MZLINK. Initialize Link Area /GCLINK/ for Store 0 NL/NS= 20 20\n", "\n", " MZLINK. Initialize Link Area /GCSLNK/ for Store 0 NL/NS= 100 100\n", "PndTargetGenerator::ReadDensityFile():\n", "PndTargetGenerator::Init():\n", " PndPipe - Using geometry /mnt/work/PandaRoot/geometry/beampipe_201309.root\n", " -I container name PndGeoSttPar\n", "[INFO] STT total number of tubes: 4542\n", "[INFO] PndStt::ConstructGeometry : STT total number of tubes: 0\n", "fwendcap & bwendcap flags == 1 / 1\n", " \n", " ====== EMC 2):: ConstructROOTGeometry() === \n", " ============================================= \n", "File name = /mnt/work/PandaRoot/geometry/emc_module12_2018v1.root\n", " \n", " ====== EMC:: ConstructRootGeometry() m3a === \n", " ============================================= \n", "File name = /mnt/work/PandaRoot/geometry/emc_module3_2012_new.root\n", " \n", " ====== EMC:: ConstructRootGeometry() m4a === \n", " ============================================= \n", "fgeoName3:: /mnt/work/PandaRoot/geometry/emc_module4_2017.root\n", "File name Bw1= /mnt/work/PandaRoot/geometry/emc_module4_2017.root\n", "PndEmc::ConstructRootGeomMod4: halfLength = 2.190000e+01, fullLength = 4.380000e+01, z-pos = -7.790000e+01, distTargetBWEC = 5.600000e+01\n", " \n", " ====== EMC:: ConstructRootGeometry() m5a === \n", " ============================================= \n", "fgeoName4:: /mnt/work/PandaRoot/geometry/emc_module5_fsc.root\n", "File name Fsc= /mnt/work/PandaRoot/geometry/emc_module5_fsc.root\n", "[INFO] ======= DRC:: ConstructGeometry() ======== \n", " -I container name PndGeoFtsPar\n", "[INFO] TGeometry will be imported to VMC\n", "\n", "0*** GSTMED *** Warning, medium = 4, value of EPSIL= 0.000E+00 reset to 1 micron\n", "[INFO] PndDrc: constructing optical geometry\n", "[INFO] PndDrc: geometry construction finished\n", " ==================================================== \n", " ======= Rich:: ConstructOpticalGeometry() ======== \n", "GetGeometryFileName() = /mnt/work/PandaRoot/geometry/rich_v313.root 313\n", "fGeoVersion = 313\n", "nOpt.size() = 3\n", " ======= RICH::ConstructOpGeometry -> Finished! ====== \n", "[INFO] FairMCApplication::InitGeometry: 0\n", "[INFO] Initializing PndStt()\n", " -I- Initializing PndSdsDetector()\n", "[INFO] Branch: MVDPoint not registered!\n", "[WARN] PndSdsDetector: New branch MVDPoint created!\n", " -I- Initializing PndGemDetector()\n", "\n", "************************************************************\n", "PndEmc::SetSpecialPhysicsCuts():\n", " using special physics cuts ...\n", "************************************************************\n", " -I- Initializing PndSciT()\n", " Fairdetector geht\n", " -I- Initialized PndSciT()\n", "[INFO] PndDrc: intialization started\n", "[INFO] PndDrc: Cherenkov Propagation is OFF\n", "[INFO] PndDrc: intialization successfull\n", " -I- Initializing PndFts()\n", " -I- Initializing PndFtof()\n", "there is gGeoManager\n", "list of sensitives has 1 entries\n", " -I- PndRich: Switching OFF Cherenkov Propagation\n", "GetGeometryFileName() = /mnt/work/PandaRoot/geometry/rich_v313.root 313\n", "Mid ref. index = 1.050\n", "[INFO] Initialize Tasks--------------------------\n", "[INFO] PndSensorNamePar::FillMap() - There are 2809 entries.\n", "[INFO] PndSensorNamePar::FillMap() done\n", " -I- PndEmcHitProducer INITIALIZATION *********************\n", "[INFO] PndEmcHitProducer: Using nonuniform light output\n", "HitProducer has EnergyHitThreshold of 0.000001 GeV and Use_nonuniformity 1\n", "[INFO] PndEmcHitProducer: Initialization successful\n", "[INFO] Simulation RunID: 1660136168\n", "\n", " Calculating cross section tables, see gphysi.dat for more information\n", "\n", " Cross section calculation concluded successfully\n", "[INFO] Monte Carlo Engine Initialisation with: TGeant3TGeo\n", "[INFO] FairGenericVMCConfig::SetupPostInit() OFF.\n", " **** GTRIGI: IEVENT= 1 IDEVT= 1 Random Seeds = 1234 0\n", " **** GTRIGI: IEVENT= 2 IDEVT= 2 Random Seeds = 1234 0\n", " **** GTRIGI: IEVENT= 3 IDEVT= 3 Random Seeds = 1234 0\n", " **** GTRIGI: IEVENT= 4 IDEVT= 4 Random Seeds = 1234 0\n", " **** GTRIGI: IEVENT= 5 IDEVT= 5 Random Seeds = 1234 0\n", " **** GTRIGI: IEVENT= 6 IDEVT= 6 Random Seeds = 1234 0\n", " **** GTRIGI: IEVENT= 7 IDEVT= 7 Random Seeds = 1234 0\n", " **** GTRIGI: IEVENT= 8 IDEVT= 8 Random Seeds = 1234 0\n", " **** GTRIGI: IEVENT= 9 IDEVT= 9 Random Seeds = 1234 0\n", " **** GTRIGI: IEVENT= 10 IDEVT= 10 Random Seeds = 1234 0\n", "=========================================================\n", "PndEmcHitProducer::FinishTask\n", "*********************************************************\n", "Read points # 114651\n", "Produc hits# 893, threshold# 0.000\n", "Hits above threshhod#263\n", "*********************************************************\n", "[INFO] *** PndEmcGeoPar written to ROOT file version: 1\n", "[INFO] *** PndSensorNamePar written to ROOT file version: 1\n", "[INFO] *** PndEmcDigiPar written to ROOT file version: 1\n", "[INFO] *** PndEmcDigiNonuniformityPar written to ROOT file version: 1\n", "[INFO] *** FairBaseParSet written to ROOT file version: 1\n", "[INFO] *** FairGeoParSet written to ROOT file version: 1\n", "[INFO] *** PndMultiFieldPar written to ROOT file version: 1\n", "[INFO] *** PndGeoPassivePar written to ROOT file version: 1\n", "[INFO] *** PndGeoSttPar written to ROOT file version: 1\n", "[INFO] *** PndGeoFtsPar written to ROOT file version: 1\n", "[INFO] *** PndGeoSciTPar written to ROOT file version: 1\n", "[INFO] *** PndGeoFtofPar written to ROOT file version: 1\n", "[INFO] *** PndRichGeoPar written to ROOT file version: 1\n", "PndMasterRunSim::Finish(): Tasks that ran just now:\n", "file is not open\n", "Opening file: ./psi2s_jpsi2pi_jpsi_mumu_sim.root\n", "\n", "\n", "Generated Events = 10\n", " - PndGeoHandling (PndGeoHandling)\n", " - PndEventCounterTask (Event Counter)\n", " - PndMasterSimTask (Master Sim Task)\n", " - PndEmcHitProducer (Ideal EMC hit Producer)\n", "\n", "831.688\n", "0.987\n", "\n", "Output file is\t\t./psi2s_jpsi2pi_jpsi_mumu_sim.root\n", "Parameter ROOT file is\t./psi2s_jpsi2pi_jpsi_mumu_par.root\n", "Parameter ASCII file is\tall.par\n", "Real time 16.156 s, CPU time 15.950s\n", "CPU usage 98.724%\n", "Max Memory 831.688 MB\n", "Macro finished successfully.\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ " Message from PHCORK(MODCOR):: initialization\n", " MODOP=5 -- corrects Energy from mass+flow\n", " Message from PHOTOS: IPHQRK(MODCOR):: (re)initialization\n", " MODOP=1 -- blocks emission from light quarks: DEFAULT\n", " Message from PHOTOS: IPHEKL(MODCOR):: (re)initialization\n", " MODOP=2 -- blocks emission in pi0 to gamma e+e-: DEFAULT\n", " MODOP=2 -- blocks emission in Kl to gamma e+e-: DEFAULT\n", "Info in : Geometry FAIRGeom, FAIR geometry created\n", "Error in : Error occured during initialization\n", "Error in : Error occured during initialization\n", "Info in : Top volume is cave. Master volume is cave\n", "Info in : --- Maximum geometry depth set to 100\n", "Info in : Fixing runtime shapes...\n", "Info in : ...Nothing to fix\n", "Info in : Counting nodes...\n", "Info in : Voxelizing...\n", "Info in : Building cache...\n", "Info in : max level = 12, max placements = 4550\n", "Info in : 2029436 nodes/ 5737 volume UID's in FAIR geometry\n", "Info in : ----------------modeler ready----------------\n", "Warning in : Called for surface LensSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for surface MirrSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for surface EVSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for surface BlackSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for surface DrcChamferSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for material surface BlackSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for material surface MirrSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface BarMirrorSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface DrcBar_ChamferSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface DrcBar_ChamferSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface DrcBar_ChamferSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface Lens1AirSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface Lens2AirSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface Lens1AirSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface Lens2AirSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface Lens1AirSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface Lens2AirSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for skin surface AirMirrorSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for surface RichMirrSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for material surface RichMirrSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for border surface BarRichMirrorSurface. Not applicable in Geant3 - setting is ignored.\n", "Warning in : Called for skin surface RichAirMirrorSurface. Not applicable in Geant3 - setting is ignored.\n" ] } ], "source": [ "nEvents = 10\n", "SimEngine = \"TGeant3\"\n", "BeamMomentum = 6.231552\n", "parAsciiFile = ROOT.TString(\"all.par\")\n", "# TString inputGenerator =\n", "# EvtGen . \"xxxxxxxx.dec\"\n", "# DPM . \"dpm_xxxxx\"\n", "# FTF . \"ftf_xxxxx\"\n", "inputGenerator = ROOT.TString(\"psi2s_Jpsi2pi_Jpsi_mumu.dec\")\n", "inputDir = ROOT.TString(ROOT.gSystem.Getenv(\"VMCWORKDIR\")+\"/macro/QA/run3/\")\n", "ROOT.gRandom.SetSeed(1234)\n", "#-------------------------------------------------------------------------\n", "# ----- Create the Simulation run manager ------------------------------\n", "fRun = ROOT.PndMasterRunSim()\n", "fRun.SetInput(inputGenerator)\n", "fRun.SetInputDir(inputDir)\n", "fRun.SetName(SimEngine)\n", "fRun.SetParamAsciiFile(parAsciiFile)\n", "fRun.SetNumberOfEvents(nEvents)\n", "fRun.SetBeamMom(BeamMomentum)\n", "# ----- Initialization ------------------------------------------------\n", "fRun.Setup()\n", "# ----- Geometry -----------------------------------------------------\n", "fRun.CreateGeometry()\n", "# ----- Event generator ----------------------------------------------\n", "fRun.SetGenerator()\n", "# ----- Add tasks ----------------------------------------------------\n", "fRun.AddSimTasks()\n", "# ----- Intialise and run --------------------------------------------\n", "fRun.Init()\n", "fRun.Run(nEvents)\n", "fRun.Finish()" ] }, { "cell_type": "markdown", "id": "f9bf2b01", "metadata": {}, "source": [ "Plotting with ROOT is still working. To show a canvas use its `Draw()` function." ] }, { "cell_type": "code", "execution_count": 9, "id": "70e8977e", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAArgAAAHYCAIAAAApvgy/AAAABmJLR0QAAAAAAAD5Q7t/AAAgAElEQVR4nO3dbZKjSJouUBibfQHTe6nZRQPLuLWXbsHKuD/eTh8PkCsUESiE0DlWViYhJAGZKX/wz3pZlgoA4Jr/evYBAADHJSgAAEWCAgBQJCgAAEWCAgBQJCgAAEWCAgBQ9N/PPgB4JdM0DcMQD766Q3ppnuemadq2jaerV0tW++8oPvZLHx5HuzqkYRimaYqzi6dt2+56pMUjub3Do48BTm4B7tb3ffzD6fv+6g5RRlZVdblcrm5fSbulT/79f61f/fDL5RJvaZpmtWUl7ZB2W12W78k/5NPr9vOvgzen6QG+o3SfOs/zdmPbtnGfnf/bixKu67rYZxiG/NUIFqti9VEn80Vx7pfLJV2EOItVeGqaZp7nVZVDOt+ffHvXdemrV9dtdQGBnxMU4MuiCNxuj0JxVUSl2vhVthiGIVXRl77oF6ruvy0dW5xX0zSrE4nt4zj+7nH956vjD+g46Qpel6AAXxZl5LaAv1rNcKMHwO3uDlflfSC2XQTaP65+Xbwl3P7SYRhK8WUYhij7V/tczTR936fYNAxDFN7RleGew96ebNpy4wirPzUcpQYR4Gt2aL6AtxHtBakEWr1aVVXf92mftHHV6HCPKF8/PYBU27/9p7360m1V/Oq9qz1LnTC2X5EO5nb/g6vHdvtHaXuyq7O48UWl4we+So0CfMe2D13c4K5urB/U5b7ruigI8y/Ny+loHMkb8ud5zsvOqtAoEN0p+r4v3a8vf8rvZVni89u2jfK767q6rlcVBvkbU/qJHeIr8qOKHVZvz092mqbIDXGyV4+/utYOAnyboADfsW19SKVmvtvVjT+3KgijqM6/ZfWNqbEgbYnCflUkf5oSSqZpSslpHMdIDJ82cGxP5GpbzP2lfmrdMCQSdiQowHdESZzflEcRe3W33a0KztVN/DRN+YHFS6tji8ECq2yR+hB875CiViA1EMzz3HXdjU/bHvaNnhyfSmetawLsS1CAb8qL3qvtDsmNG9xpmna5/Y0ugXVd13X9vSGIaaKkH4abKP5TC8U4jrdPPz/s0uDSe743dWA88lAReEWCAnxTPpyh1MQQW66Wf6Hrup9PLdC27TiOUaURjfffuKtO8yLknRvuUZqBcRiGqw0c+Q4RDn5y2CF1TZASYHeCAvxIVHdfbXcIVzvohTQDwQ+PIc3mlMrs/Ou2W9L27W5RVH8pu6RJC7ZuF9tx6fLD/h5dE+ChBAX4vggHUciVirqocsgnE0yiPP5h//wbnShzq7I8Fa4raQjD/Ud14/Rvt8hU1yanuvNL87ekwPHV9wL3EBTg+4Y/izxV5eKwbds0W3O030evvbquqz1qy1O3yihlo9V/dYcd9QR1XadxiVG4Xv3q7YyK0YegdADpIqSxkdOfJaPS8lfbzw+rMZzxpffEhbRPqvxoCz79KOAT+0/NAOe1mkxpuTa9z3afpbB20Y1JgT6dcOnGhzdNs52RaXsA+ZRQq++Kt6+mRVp93dUDvn2C6RhWMzXlh5TvcPWLlqzmYLmjFqF0hYE71Yv6OvgtaYzDI25248PTJ6+eXt1nd5+e4PYAtlvSpM6ffotZleAXCAoAQJE+CgBAkaAAABQJCgBAkaAAABQJCgBA0X8/+wAAOJsbM3RxHHcOexQUANifsfcHd3+Y0/QAABQJCgCc3GqtVL5EUADg5PLlx1aGYTAX+G36KADwvtQ0fEpQAOAtTH/EwmOxJdZkz+sVrq5tlrbEbvHSds98t1jt7AzVFc9dvBKA8zla4VJVVayEntZDjzXQV0ufL38WTG+aJh6kpdLTG1cP0tNYtz1fzz299Iwz/tz9B6aPAgDnN8/zsizTNC3L0jRNWqk8MkF6mnabpqnv+3Ec8w+5XC7xCasPrD42YUTsmKbpcrlUr78euqAAwPml+/uqqtq2jRaHlXEc892igM+L+dQSkaocrkqhoW3bpmlWaePlCAoAnF8q42+b57nO5C/dSAY//94j05kRAP6jaZpVS8EJSvofEhQAoKr+1BnkySCNcfjS58R4h/R4l2N7Ik0PAFBVfzoz5j0MxnH8Ro1C13X5B+ZDIV6RGgUA3lekgbqul2Vp27bv+1TMV1UVwxa+qmma1L9h25bxcurFAl8A7CrK3Wcfxb3SLEk3ttxpGIZxHGNs5Pc+4dfc/2ekRgGAt7Ytzn9ewB85InzV2YLC/QtsAwCfOltQqKrqheq7AE7pbe/ZzlSRkLxSM9I9XqthDOCU3jYovBZ9FAB4GvdsB3d/mDthUCidvL+1APBVJwwKAgEA7MXMjABA0QlrFAA4ppiPaLv92zXBwzBM0/Tq6ykMw9C27WrERDqv1Uul7au3x2fucniCAgC/6qvzIk/T1HXd1TBxjuGIEZ5Wi1GN4xiLVI3jeLlc4tW2bed53m6fpilt352gAMCv2rF0396Iv5Yo+LfbV+EgclKkgbR9GIaUn6KO4UE1K/ooAHAIdV3HAs11Xdd1HWspRXVCenX4Iz1dVb/Xf6SN0zSljUdLFbEM1Wrjam3raZryNJC2p+sT/39QdUKlRgGAX7a98U2FX9d1ccccde9RYdD3fay0VP1ZuLlpmtgtbrLjvXVdx/bYLc2/13VdWsKx67p416+c6OfiqFb9NqLUL/VFWIloVVXVPM8pHqVah30s53K+MwJ4OaWf4u3dc14SVVXV933+IZfLZVmWKPtjY9w3rz5wtc/q7fnHXi6X2HgoqxNPdQNN08TjeDXOMe2Zv5TyweVyie2fnub9xaUaBQB+1VIe43DPffDVOva4+Y4b9Hxj27ZN04zjGI9XOxxT1JGkqxSVK9HIEpUrUQMR1yGuWH5J05nuVXFywj4KdcGzjwuA39P3fRoREPfi4zimrg9HlioSwmo0RFQtLMtyIweU+kh+zwmDwo1qFgBOKQ0ESNLGuMOOLoFxR/7MA73DjWqV1Xmlx9s5GHbs23jCoADAkU0bP//MFBTiaT6zU9d1x69FyEV9QLoswzCkUj/aIOJxdNKMx/M856Mk8qc7uLMvw6s43xkBvJzST3GpM2PqdZh3wcufpt2iWn71gfF4NY9Teu/qS/NugwexParVMV/dnl+H1bnfc473F5f/GT1yGmk8DCR1/fd247L89ftHAm/iET/FaRzgp7tV12rvS9uP7Kvn8qVzvP/P6CjFampSWm282k81TbKxfYugwFZd/72KBdstwI78FB/f/X9Gh+ijME1TjF3JN9Z1HS1M0U81bY9JK6s/k20dZ94MADifJweFqDCIgj+XBoamqStT5UHMgB39X9JkWwDAIzy/RuHqTNfVxyk1YjLLajMDdszl+fhjBIA39eSgEP0MrvZOmOc5tuePt2+vHrZeFgBw0Cmc06Sb0U0hn2Dr00kkvjEJo043AHDV85seropahLQcSD7FxKfuHBi6Gk4KwKPFhPpXV498icmVd1EatZe259cnevLlbn94jArc94APGhTGcUy1CGkZjGozf/UrjosFeHPboPAmvc3SYL22bVcD+tJAv2mauq5LlyhfR/seXdftPkf1QYNCVSj+t9NZ/8qxALCbVUn2Pr/kMRlz3PTnA/ri/zHQL1rY02DAeHrndNcPum0+aFDIL1NVVeM45utp5kMld1z3AoBHi2FueYE3DMNq7FvceYe88Mu3500YseZTtF+s3nI0pXaB0nKR26c3Pnme54eUid9o0X+E6uMU38vHfgM3ZrTefs7jD5YXU1X/79MtwI5KP8VVVfV93zRNvhhB/P5X2QoF6fF2e9M0l8tltehDFAf5W1YFygFFNtoe59VTTkVh6bziXcvHxS9uu7+4PHSxevmj9NJ2u6DAlqAAv+x2UEil2pKVcKl0XC37lO+cx4u8RFzdT1aHXPkpSbe724Pc3gPH05SNbl/Y5TFB4aDDI8ON+pYj1ywBcEMa7h6d/FftDtF3L1XRrxopUlP9qotfXigcvEm6bdtlWaLTYvWxMSJtTwsx5OkhWl62kw/FhAKPGzNy0D4KAJxYKtjmed7e+K0a7FOSqOs6RgRs48XLyQf0rbZHlcPVfotpnuLcPM9xGWMwRfWn08Zeh3roGoXvKU24tJgvAeAYYnm/1az8IVJCfn+8GhoQG19xrERd15fLZRuMVrMMpx0iEuVnerW74rZK5p4ZF77gziaKV3G+M+Ln9FGAX1b6Ka4+dtOryh0YUy+06k//g1Xre16EVR/b+1edJY+j2vSluHpqeV+E/C2rKxO9PVZf8Yg+CpoeAHiCfNB7Lurko50+aojzRQHT2MgoEV+rv9rlcpnnOZ1Cvt5h0zRpe0xMvH1L13X5gga7T6xU8p/uEqeROoBAUtd/L8tft7cAO9rlp/jq3Lv5xhednLd02Numh0/f8hP3/xmdrVgVFNgSFOCX+Sk+vvv/jDQ9AABFggIAUCQoAABFggIAUCQoAABFZmYEAIpOWKNwYxYqAJ5rGIY6c+dSRqU5ieuP9l3j4AiGYYils1bnFbM7l843Xt3rGE4YFAA4pmEYxnFM8yvHqkipSPveCgVN08QqzDFXY1pC4gTquo7pF2NJyZQJYqWMeNp13fai5Tvv4M6pnl/F+c6In7PWA/yy0k9x9XFRhuXj2gRN0+RLIeRKL93+wJd2ewGItMpDTPa82q36uKjEVfdfJTUKAPye1Z3uMAxR1LVtm5ZLjpfatk0NCvd/fr7UZHxd3jCRPjmvdYh9vnMyD7Zabjt/KT+X6uP0z/svwH13uHkN5zsjfk6NAvyy0k9xKsOuru6YVxtEGRn3zfH4zhqF5eOKi9XHRSnj8fZm/dP77+fKD361Pa9pSO45o/uLy7MVq4ICW4IC/LIbP8V93+c3yqvloaN4W62nvGwWaM63b4vP9Dmr8jKvpa8+rmS9Xa/5OEo39ukybg9+36BwwuGRABxW3igQfRvHcVyuDUzLK9vzbPGpeZ5j/3meV9+Y7zZNU2qDOPISlMuyRGfG1TJO0x9d110ul8edgj4KAPySvGdADO2Lu/zdxzSmUnPVzJ+3faQBBfs36u+tbdurFyoFnYcOChUUAPgl4zhui7p73hh1A/fIOzNGShgyq92maZrn+ZjDKVeDRdPjmIjiN4/khEGhLnj2cQG8u6ZpVkP8o/y72qU/LxpvfOaUads25mlIb5znOX1dXdf5V8fBfKlR4zfFMJD8aZVVIazaUx7bdHJnX4ZXcb4z4ud0ZoRfduOneFswp5eigM/7MyY35lFY7bbq2bdqVti+dORujKtrlQ51dVI3unPecH9x+aFnxAms+npAVVV1/fey/HV7C7CjT3+KS/fBq+173S7/xm33Y+QTJFx96dsndX9xebZiVVBgS1CAX+an+Pju/zM6YR8FAGAvRwkKV/uq5INc79kOAOzrEEFhmqarY2a6rquqahzHfMxCrJpV/Zmd+2QrigKcQ2kAGgdx/x/lk2dmjJm5tgNkY2xraj6ps1XGx3FMU1BdXaUbgOd65w4K9en6Zzy/RuHqUlfDMOTDQpZliTSwmmvzasgAAPby5BqFNPNUTKWZzPPc931MoFGVZ9tIy2u+4qAXADi+59colIzjGBN4RR+FSAzTNH06i9ZDm2oA4K0cevXI1MwTHRvvbPU5WeMQADzRcWsU8o4LqWVhNff16062BQAv4aBBodS+sMoExjsAwEMdNCgMw5DPrDCOY0SHCAqpb2PaDgA8wkGDQtu2sQBodDZsmiaFhsvlEt0boxOiSgUAeJyjdGbc9kBMwxyqjy0ObdumaRX0TgCAhzpKUCgpRQERAQB+wUGbHgCAIxAUAICiozc9fENppkUTMQE8SF3/vdqyLH895UjY3QmDgkAA8PvyZLDNDbwuTQ8AQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAUCQoAQJGgAAAUnXDCJTMzAsBeThgUBAIA2IumBwCgSFAAAIoEBQCgSFAAAIoEBQCgSFAAAIoEBQCg6ITzKJhwCQD2csKgIBAAwF40PQAARYICAFAkKAAARUcJCsMw3HhpmqZ8yzRNbdveeAsAsItDBIVpmsZxXKWB0kvDMHRdFy/VdX31XQDALp4cFKJuIAr+q7YvjeN4uVymaZqmqWka9QoA8DjPr1Fo27bv+9JLTdM0TZO2RCxo2zY9nef54YcIAO/qyUEhuhpcrRWIEHC7ZSESg9YHAHiQ4064FE0Mq43R3HD7jaWZGW8wRxMAXHXQoFDXdd/3qYnhS5T6ALCX5/dR2IqmhOjn2LbtPM/jOEZoiKerPb+XJwCATx2xRmHVvXGe56ZpUlAYxzG9pHcCADzUEYNC9XH+pXx6pYgLqf/jOI6fdlkAAL7tiE0Pt10ul3Ec67qOTosqFQDgcY5So3CjB+IqCrRtuyyL3gmvoq7/Xm1Zlr+eciQAfMNRgsJXiQgvJE8G29wAwJG9alCAG8QRgL0ICpyTBg6AXZwwKJRmZjQREwB81QmDgkAAAHt5veGRAMCvERQAgCJBAQAoEhQAgCJBAQAoEhQAgCJBAQAoEhQAgKITTrhkZkYA2MsJg4JAAAB70fQAABQJCgBAkaAAABQJCgBAkaAAABQJCgBAkaAAABQJCgBA0QknXDIzIwDs5YRBQSAAgL2cMCjwRHX997MPAYA9CQrsbFn+evYhALCbo3RmHIbh6sa2bbcvTdN0dTsAsK9DBIVpmsZxnKYp31jX9TiOVVWN41jXdXp1GIau6+Jd+XYAYHdPDgpRNxAFf65t26qqlmWZpmlZlqZp0j7jOF4ul2mapmlqmka9AgA8zvNrFNq27ft+tXGe53xj5IbqTwtF/nSe58cfIwC8qSd3ZmzbNkr9aGVIVkMcS+0L8d6olnjI8QHAe3uBUQ9t287zfLlcqqqK5obb+5cmXLrB1AsAcNXzmx5uGIahrutICffXGSxf98iTAIAXdtyg0LbtOI593y/LklJC1C6kfaJJQrsDADzIQZseopfi9l4/0kN6amwkADzUQYPCOI5N06xyQOr5OAxDDH+I3Z5wfADwHg4aFKqqmud5Nb9CVDBcLpeu61K9gkoFAHicowSFVSvDjQ6GbdvGREyV3gkA8GBHCQpfJSIAwC847qgHAODpBAUAoOhVmx5uKM3MaGIlAPiqEwYFgQAA9qLpAQAoEhQAgCJBAQAoEhQAgCJBAQAoEhQAgCJBAQAoEhQAgKITTrhkZkYA2MsJg4JAAAB70fQAABQJCgBAkaAAABQJCgBAkaAAABQJCgBAkaAAABQJCgBA0QknXDIzIwDs5YRBQSAAgL1oegAAigQFAKBIUAAAio4SFIZhuLqxbdtpmlbbp2lq2/bqWwCAHR0iKEzTNI7jKhDUdT2OY1VVXde1bZu2D8PQdV28q67rbYwAAPby5KAQdQNR8OeitmBZlmmalmWZ5zkFgnEcL5fLNE3TNDVNo14BAB7n+TUKbdv2fb/aGCEgPU2BIP6fKhiGYZjn+VcOEwDe0ZODQnQ12NYKzPOcNze0bXs1EMQ+Wh8A4EFebMKlVU3DVaWZGW8wRxMAXPViQeEeSn0A2MuLBYW2bWMoRIhGh7yRAoAjqOu/V1uW5a+nHAk/dNCg0DRN3vMgtThcDQoAHMo2E2xzA6/i+aMerorhDJEDpmma5zkf75A6P47j+GmXBQDg2w5aoxBjJtP8Cn3fp/aFy+XSdV2qV1CpAACPc5SgsO2BGMMmY0amfHvbtjERU6V3AgA82FGCQkkpCogIAPALDtpHAQA4AkEBACg6etPDN5RmZjQREwB81QmDgkAAAHvR9AAAFAkKAECRoAAAFAkKAECRoAAAFAkKAECRoAAAFJ1wHgUTLgHAXk4YFAQCANiLpgcAoEhQAACKBAUAoEhQAACKBAUAoEhQAACKBAUAoEhQAACKTjjhkpkZAWAvJwwKAgEA7EXTAwBQJCgAAEWCAgBQdOigME3TMAxt207TtHqptB0A2NFxg8IwDF3XRRTouq5t2/RSXdfjOG63AwD7Om5QGMex7/tpmqZpulwu8zxHaBiGoaqqZVmmaVqWJW0HAHZ36OGRqbYgrzaYpqlpmvS0aZphGGSF31HXf283Lstfv38kAPyO49YoNE0TTQ/TNEVQiP/P85znhrZt53l+ziG+pWX5K//v2YcDwGMdt0Zhmqa6rruui6eXy+XON5ZmZrzBHE0AcNVxaxTqum6aZlmWZVn6vk8dGz+1fN2DTwUAXtVBg0JkgpQMogOjjggAr6uu/1799+wj4i7HbXpYaZomgkJ6EFZ9GwE4oG2XJkHhVRy0RiG6K0ZFQlVV0zTN8xxPh2FIQyLz7QDA7o5bo3C5XLqui4mVqqrq+z6NfYguC6vtAMDujhsU2raNWZWqj/MoVFU1DEPMnSAiAMBDHTcohBtRQEoAgEc7aB8FAOAIBAUAoOjoTQ/fUJqZ0cRKAPBVJwwKAgEA7EXTAwBQJCgAAEWCAgBQJCgAAEWCAgBQJCgAAEWCAgBQJCgAAEUnnHDJzIwAsJcTBgWBAAD2oukBACgSFACAIkEBACgSFACAIkEBACgSFACAohMOjzyluv57u3FZ/vr9IwHgrQgKL2MVC65GBwDY1wmDgpkZAR7KjcpbOWFQEAgAHk3T5/vQmREAKBIUAICioweFYRjath2G4er2aZqecEwA8DYO3UchuiU2TTOO4zRNKRak7V3XNU0jLgDAgxy3RqFt26ZplmWZpulyuczzHIEgahdi+7IsaTsAsLvj1ijM83y5XOJx27ZpLMM0TU3TpN2aphmGQVYAgEc4blCoqirvnZAezPPc932+zziOv35oAPAWDhoUooagruuoPJjneRzHOydIKE24dIOpFwDgquP2Uaiqqu/76MMYBfl27MNVy9c99jQA4GUdtEYh5MnA6IbdWWgKgE8dtEahbdvqTwNEmOc5HqwSw6pvI1+yLH/l/z37cAA4nOPWKOTDGfKBkcMwdF03TVNMuJQPjgDgEawC9c6OGxSmaarrOvVM7Ps+qhnatu37vuu61XYAHkel49s6blCo/syqVP1piUiGYYjKBhEBAB7q0EGh2kSEO18CAHZx0M6MAMARCAoAQNHRmx6+oTQz47tNrGSaBAB+7oRB4d0CwQ2rWGCAEwBfpekBACgSFACAohM2PfDLtGgAnJigwI/oHQlwbpoeAIAiQQEAKBIUAIAiQQEAKDphZ0YzMwLAXk4YFAQCANiLpgcAoEhQAACKBAUAoEhQAACKBAUAoEhQAACKBAUAoEhQAACKTjjhkpkZAWAvJwwKAgEA7EXTAwBQJCgAAEWCAgBQ9BpBoW3b1ZZhGNq2nabpCUcDAG/jBYJC27bzPOeZoK7rcRyrquq6bpshAIC9HH3UwzRN8zznW4ZhqLKhDXVdT9MkLgDAIxy9RqHrur7v8y3TNDVNk542TRPRAQDY3aGDQtu2fd+vcsA8z3n9QTRM/PKBAcCbOG7TwzAMq64JdyrNzHiDOZoA4KqDBoVpmsZx/F75rdQHgL0cNChEc0PexNB1XdM07zMesq7/fqGPBeCsjhsU8kwwz3PqtLiKC6u+jWeyLH8d/AMBOL2DBoW2bfPqhHEcY4alqqqGYei6LoZExuDJy+XyrOMEgHM7aFC4IYZCdF0XT/u+N4kCADzIawSFVf/EYRiibUJEeEXbfhLaRAAO6zWCwlVSwivaZgL9KwGO7NATLgEAzyUoAABFL9z0UFKamdFETADwVScMCgIBAOxF0wMAUCQoAABFggIAUCQoAABFggIAUCQoAABFggIAUHTCeRRMuAQAezlhUBAIAGAvmh4AgCJBAQAoEhQAgCJBAQAoOmFnRgBeQl3/nT9dlr+edSTcICgA8ASrWLAKDRyHpgcAoEhQAACKTtj0YGbGn9BkCGgFIHfCoCAQfJsmQyC4SSDR9AAAFAkKAECRoAAAFB06KEzTNAxD27bDMKxeiu3TND3hsADgbRw3KAzD0HVdRIFxHPOxDHVdj+NYVVXXdW3bPukAAeD8jhsUxnHs+36apmmaYiBD1CvE/5dlie3zPKtXAIAHOW5QqKoqry1omiYCwTRNTdPk27cNEwDALo4bFJZlyYPCPM/xND0IbdvO8/zbBwcA7+G4QSGZpik6KNxZc1B/3WNPAABe1tGDQtu2Xdc1TXP/fIvL1z30FADgdR16Cue4179cLoY2AMBTHDco1HWdOjDmVhtXfRsBgB0dNChEFNhOqRSTL8X8CvHqPM+Xy+UpB7mjXZZfsoYTALs7dFAYxzEmVgpRl9C2bd/3XdfFxr7vz9EwsctabRZ8A2BfBw0KwzDcGOMQr0Zo+L1jektqKQDe3EGDwj2khEdTPwHA0YdHAgBPJCgAAEUv3PRQUppp0cRKAPBVJwwKAgEA7EXTAwBQJCgAAEUnbHqAn7s6gYTxosAbEhTgulUsMPcU8J4EBd7UquBXWwBwlaDAO1JbAHAnnRkBgCJBAQAoOmHTg5kZAWAvJwwKAgEA7EXTAwBQJCgAAEUnbHrg5Xw6OtEkBwDPIijwZJ+GAJMcADyRoPDClKDBdQB4HEHhVamNz7kaAA8iKLykuq6NAr3hF67PSy8v6e/Pba7Pba7PuzlhUDjChEvbUsTiAufjzxQe6tMfUn7HCYPCQaJu/hf6pe8+AZ7l0x9SfsEJgwKn5/cC4NcICrwk9TEAv0NQ4AWoQoB34F/6Mb1qUBiGoaqqtm3btr2x2169c4/2OXs52nld/ZwnVh68xPV5oqOdl+vzO5+zl9XxfPtf+lmvz3G8XlCYpqnruqZpqqoax7Hv+wgN8HQ37od0aOXI3Mpzw+sFhWEYmqaZpikej+MoKHAEN0r9uv7f7at+mjkUsZWS1wsK8zxfLpd4HEFhmqbbDRCwi92L9ns+8JQ/368Skozjfy3q7R7kxYJCVCSsYoGgwC/Y/efmng98lQL1Gz6Oj//fJx7JbcbxvxbToD3Ci/XdiA4K+THXdZ1aIqrytIwA8Dteq2D91IvVKFyVVyec7I8HAJ7rv7lfLOIAAAaBSURBVJ59AADAcb1YUIjKg9TQkG8EAHb3YkGhqqqmabqui8dp2qUnHg8AnNjrBYWoTqjruq7rcRzTUMn3MQzDMAyrapW3Mk3TMAxt226n0Ijt24sTQ2PebcqNbYZ2fYK/Pzekf1/b6/Dm1+fqCX7jmrzeb/jymi6Xy+VyefZR/LZIRU3TxMSUfd8/+4ieoO/7/CLkf4fjaWxvmqb0ljf5m7M9WdcnlK6D67NsTtb1SeLnd3V2X70mL/ob/qpB4T3FX694HH8Rn3o4z7H615Weri5I/o8zf5xfwxNLNW3pxF2fkJ9g/tPv+oT835frEy6XS7otyYPCN67Ji/6Gv8ZRErY3iGcN7zeszjr9w1v9QqWnq3+N8dv3S8f6PPFzf+MX/G2vT+lfjesTSj8y73x9LpdL3/erf1DLt67J9vK+RKXC6/VReFulWSmfcSzPtCxLfhHmeY6n6UFo23ae5+3brw6cOZm2bbeLpbk+SbQch7TR9QlN03RdN01TmvHWv6/0F2a1/avX5Opv+Es4w4RL7+xk/yC/JKbprAo9jPLdUrXhOxiGYZ7n+/9ivNX1SV2h45TneR7Hcbk5S9tbXZ+qqqZpqus6jSz7tLf4u12fe9x/TZpsWuEjU6Pw2l4xnO6ibdtYbfz2r/y7mabp05KPvu/j9i4u1On76n9JpKiocO77PmoXnn1Qp3W1BuKABAVeT13XsYjoPT9hq/rA1639u0eaWSRUVdV13e2TfavrE/Jk8Okt3Vtdnzi7dEHiQrk+X/Wla/IS10pQeBlmpQzpjmd17qtf/FT791a9OoZh6Ps+DwpN08TPvetTXftHlH7QXZ+r0mVxfba+ek1e+Df8KV0o+Z74WxiPX2hozY6ixbTv+8tHy8ehXKsRz1XWtbj6ONz53KrNAG7XpzQ80vUJVWF4pOuzbMYsfOOavOhv+GscJUke8t5wbGT801pZjUsO+aCjVYespxz5U1TXhn27PvnJ5tfB9Vk2J+v65La/ut+4Jvn2V/kNrxf9nl6NVsDb0rCu7fbKdXN9qqq6ebKuT+X6fNFXr8nLXStBAQAo0pkRACgSFACAIkEBACgSFACAIkEBgKqqqrqu67reTprUtm1d1/mMlrEkRJJeGoahLvj5XEz3fEg6sNvH+XOrT37EVxyHoADA/9kWxqslCYZh6LouzXvW9/04jjHYbxiGNA1aVVVN06Snu4wG/DQoRDkdo/li3bh8/qh0nNWfTPPVA6jr+t///ndVVW3brqZ9O/PiWM+bwgGAA7laLqS5g/KpBvPJhZbCJIPV3rM0br93K595M38c8uP83sSIVVX961//Kr101lkp1SgA8B9RfOY37rGASP602iy5udqnZPgjNSJEo0bIqxzyiv2rVRFpNZPVxnme53lOL23rQiL3DMMwjmNVValSIT+SdHYxk1IcTKqB+Mc//rE92XjptEtdPDupAHAIVVX1fd80TX7jXlVVWmNluXabfuPTVntG5Xy0Ryx/7ulLayXEe/Ovzo+hVHht1/Ko/qwOs92z+li7sD2qFAjy2pR//vOfVz/qVeZj/gY1CgD8n3S3XV2bbHh1j/4NacLjaOZPj9PX5ctbt23b931+px7VBssdcwrH5zdNM45j13WryonV43Qk296Iy7KslibPXx2GIVa9f6Epmb/qv599AAAcSCqz27bdtin8sMte/vb0RSFtT+s1h7yQTgnmTtFwkL5lHMe6vrJwQX4kX/qK2D/Wdv/Sgb0WNQoAfNA0TRTPeXt/iBv67Vu+MTIwWv2vdnq4IYr5ewrm/DMj9MR7t9+VH8lq4cfbuq67+oEnIygA8EFUp0f5tw0K1abX3veWQ4x78WmahmG48d58bcbUgWCe5097Do7jeGcRfueRrMSe9zSCvDpBAYAPogiM4nP7UtM0XdflQwO6rmua5hvV76mwX42zyHNA3LVfPYa0/9VAEBM8rPo3VJ9VANyZLeIgI17k7nnv63lmT0oADqPKuvevevJXH+cwWPVUuDoOYrt9NWIiTyHxUrUZbpCXU6tjqLLREGmf1Vdse1SkM0pNDDeOZDvXQuzzz3/+s9RX45RTKVzp1gEAn/pei8OND8lbGfb6/Bufs/rqfJ/VkeT+/e9//8///M8PD+m1CAoAQJE+CgBAkaAAABQJCgBAkaAAABQJCgBAkaAAABQJCgBAkaAAABQJCgBAkaAAABQJCgBAkaAAABQJCgBA0f8HJkJZ2JXG4zkAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "myfile = ROOT.TFile(\"psi2s_jpsi2pi_jpsi_mumu_sim.root\")\n", "mytree = myfile.Get(\"pndsim\")\n", "can = ROOT.TCanvas()\n", "mytree.Draw(\"MCTrack.fStartZ\",\"\")\n", "can.Draw()" ] }, { "cell_type": "code", "execution_count": null, "id": "6990849b", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "@webio": { "lastCommId": null, "lastKernelId": null }, "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" }, "vscode": { "interpreter": { "hash": "aee8b7b246df8f9039afb4144a1f6fd8d2ca17a180786b69acc140d282b71a49" } } }, "nbformat": 4, "nbformat_minor": 5 }