Logo ROOT  
Reference Guide
No Matches
pythia_display.C File Reference

Detailed Description

Demo showing H -> ZZ -> 4 mu generated by Pythia.

Requires libPythia6.

void pythia_display()
TString dir = gSystem->UnixPathName(__FILE__);
gROOT->LoadMacro(dir +"MultiView.C+");
#ifndef R__WIN32 // libPythia6 is a static library on Windoze
if (gSystem->Load("libPythia6") < 0)
"Could not load 'libPythia6', make sure it is available!");
gROOT->ProcessLine("#define __RUN_PYTHIA_DISPLAY__ 1");
gROOT->ProcessLine("#include \"pythia_display.C\"");
gROOT->ProcessLine("#undef __RUN_PYTHIA_DISPLAY__");
// Constants.
const Double_t kR_min = 240;
const Double_t kR_max = 250;
const Double_t kZ_d = 300;
// Solenoid field along z, in Tesla.
const Double_t kMagField = 4;
// Color for Higgs, Zs and muons
const Color_t kColors[3] = { kRed, kGreen, kYellow };
// Global variables.
#include "TEveTrack.h"
#include "TEveElement.h"
#include "TEveGeoShape.h"
#include "TPythia6.h"
#include "TGeoTube.h"
#include "TMCParticle.h"
TPythia6 *g_pythia = 0;
// Implemented in MultiView.C
class MultiView;
MultiView* gMultiView = 0;
TEveTrackList *gTrackList = 0;
// Forward decalarations of CINT functions.
void pythia_next_event();
void pythia_make_gui();
// Main - pythia_display()
void run_pythia_display()
if (g_pythia != 0)
Warning("pythia_display()", "Already initialized.");
// Create an instance of the Pythia event generator ...
g_pythia = new TPythia6;
TPythia6& P = * g_pythia;
P.SetMSEL(0); // full user controll;
P.SetMSUB(102, 1); // g + g -> H0
//P.SetMSUB(123, 1); // f + f' -> f + f' + H0
//P.SetMSUB(124, 1); // f + f' -> f" + f"' + H0
P.SetPMAS(6, 1, 175); // mass of TOP
P.SetPMAS(25, 1, 180); // mass of Higgs
P.SetCKIN(1, 170.0); // range of allowed mass
P.SetCKIN(2, 190.0);
P.SetMSTP(61, 0); // switch off ISR
P.SetMSTP(71, 0); // switch off FSR
P.SetMSTP(81, 0); // switch off multiple interactions
P.SetMSTP(111, 0); // Switch off fragmentation
// Force h0 -> ZZ
for (Int_t i = 210; i <= 288; ++i)
P.SetMDME(i, 1, 0);
P.SetMDME(225, 1, 1);
// Force Z -> mumu
for (Int_t i = 174; i <= 189; ++i)
P.SetMDME(i, 1, 0);
P.SetMDME(184, 1, 1);
P.Initialize("cms", "p", "p", 14000);
// Create views and containers.
TEveElementList *fake_geom = new TEveElementList("Geometry");
b = new TEveGeoShape("Barell 1");
b->SetShape(new TGeoTube(kR_min, kR_max, kZ_d));
b = new TEveGeoShape("Barell 2");
b->SetShape(new TGeoTube(2*kR_min, 2*kR_max, 2*kZ_d));
gMultiView = new MultiView;
gTrackList = new TEveTrackList("Pythia Tracks");
TEveTrackPropagator* trkProp = gTrackList->GetPropagator();
// Next event
void pythia_next_event()
TPythia6& P = * g_pythia;
int nh = P.GetMSTU(72);
// printf("N = %d, Nhard = %d :: NumSec = %d, separators (%d,%d,%d,%d)\n",
// P.GetN(), nh, P.GetMSTU(70), P.GetMSTU(71), P.GetMSTU(72), P.GetMSTU(73), P.GetMSTU(74));
// 2->2 hard postfrag final
TEveTrackPropagator *trkProp = gTrackList->GetPropagator();
for (Int_t i = 0; i < 7; ++i)
TMCParticle& p = (TMCParticle&)*MC[nh+i];
TParticle pb(p.GetKF(), p.GetKS(), 0, 0,
p.GetFirstChild()-nh-1, p.GetLastChild()-nh-1,
p.GetPx(), p.GetPy(), p.GetPz(), p.GetEnergy(),
p.GetVx(), p.GetVy(), p.GetVz(), p.GetTime());
TEveTrack* track = new TEveTrack(&pb, i, trkProp);
track->SetName(Form("%s [%d]", pb.GetName(), i));
if (i == 0)
else if (i <= 2)
printf("%d - %d %d %d %d %d %d\n", i,
p.GetKF(), p.GetKS(), 0, 0,
p.GetFirstChild()-nh-1, p.GetLastChild()-nh-1);
printf("%d - %f %f %f %f\n", i,
p.GetPx(), p.GetPy(), p.GetPz(), p.GetEnergy(),
printf("%d - %f %f %f %f\n", i,
p.GetVx(), p.GetVy(), p.GetVz(), p.GetTime());
TEveElement* top = static_cast<TEveElement *>(gEve->GetCurrentEvent());
// GUI stuff
class EvNavHandler
void Fwd()
void Bck()
void pythia_make_gui()
// Create minimal GUI for event navigation.
TEveBrowser* browser = gEve->GetBrowser();
TGMainFrame* frmMain = new TGMainFrame(gClient->GetRoot(), 1000, 600);
frmMain->SetWindowName("XX GUI");
TString icondir( Form("%s/icons/", gSystem->Getenv("ROOTSYS")) );
EvNavHandler *fh = new EvNavHandler;
b = new TGPictureButton(hf, gClient->GetPicture(icondir+"GoBack.gif"));
b->SetToolTipText("Go to previous event - not supported.");
b->Connect("Clicked()", "EvNavHandler", fh, "Bck()");
b = new TGPictureButton(hf, gClient->GetPicture(icondir+"GoForward.gif"));
b->SetToolTipText("Generate new event.");
b->Connect("Clicked()", "EvNavHandler", fh, "Fwd()");
browser->SetTabTitle("Event Control", 0);
#define b(i)
Definition RSha256.hxx:100
int Int_t
Definition RtypesCore.h:45
short Color_t
Definition RtypesCore.h:92
constexpr Bool_t kFALSE
Definition RtypesCore.h:101
double Double_t
Definition RtypesCore.h:59
constexpr Bool_t kTRUE
Definition RtypesCore.h:100
@ kPink
Definition Rtypes.h:67
@ kRed
Definition Rtypes.h:66
@ kGreen
Definition Rtypes.h:66
@ kCyan
Definition Rtypes.h:66
@ kYellow
Definition Rtypes.h:66
void Error(const char *location, const char *msgfmt,...)
Use this function in case an error occurred.
Definition TError.cxx:185
void Warning(const char *location, const char *msgfmt,...)
Use this function in warning situations.
Definition TError.cxx:229
R__EXTERN TEveManager * gEve
#define gClient
Definition TGClient.h:157
@ kDeepCleanup
Definition TGFrame.h:42
winID h TVirtualViewer3D TVirtualGLPainter p
#define gROOT
Definition TROOT.h:407
char * Form(const char *fmt,...)
Formats a string in a circular formatting buffer.
Definition TString.cxx:2467
R__EXTERN TSystem * gSystem
Definition TSystem.h:560
An array of clone (identical) objects.
Specialization of TRootBrowser for Eve.
A list of TEveElements.
Base class for TEveUtil visualization elements, providing hierarchy management, rendering control and...
Definition TEveElement.h:36
virtual void AddElement(TEveElement *el)
Add el to the list of children.
virtual void DestroyElements()
Destroy all children of this element.
Wrapper for TGeoShape with absolute positioning and color attributes allowing display of extracted TG...
void SetLineColor(Color_t col) override
Set the line color.
Definition TEveLine.h:48
void AddElement(TEveElement *element, TEveElement *parent=nullptr)
Add an element.
void AddGlobalElement(TEveElement *element, TEveElement *parent=nullptr)
Add a global element, i.e.
TEveBrowser * GetBrowser() const
static TEveManager * Create(Bool_t map_window=kTRUE, Option_t *opt="FIV")
If global TEveManager* gEve is not set initialize it.
void Redraw3D(Bool_t resetCameras=kFALSE, Bool_t dropLogicals=kFALSE)
TEveEventManager * GetCurrentEvent() const
A list of tracks supporting change of common attributes and selection based on track parameters.
Definition TEveTrack.h:140
void SetMarkerStyle(Style_t s) override
Set marker style for the list and the elements.
void SetMarkerColor(Color_t c) override
Set marker color for the list and the elements.
void SetMainColor(Color_t c) override
Set main (line) color for the list and the elements.
void SetMarkerSize(Size_t s) override
Set marker size for the list and the elements.
void MakeTracks(Bool_t recurse=kTRUE)
Regenerate the visual representations of tracks.
TEveTrackPropagator * GetPropagator()
Definition TEveTrack.h:175
Holding structure for a number of track rendering parameters.
void SetMaxR(Double_t x)
Set maximum radius and rebuild tracks.
void SetMaxZ(Double_t x)
Set maximum z and rebuild tracks.
void SetMagField(Double_t bX, Double_t bY, Double_t bZ)
Set constant magnetic field and rebuild tracks.
Visual representation of a track.
Definition TEveTrack.h:33
void SetAttLineAttMarker(TEveTrackList *tl)
Set line and marker attributes from TEveTrackList.
virtual void SetStdTitle()
Set standard track title based on most data-member values.
virtual void AddFrame(TGFrame *f, TGLayoutHints *l=nullptr)
Add frame to the composite frame using the specified layout hints.
Definition TGFrame.cxx:1117
void MapSubwindows() override
Map all sub windows that are part of the composite frame.
Definition TGFrame.cxx:1164
void SetCleanup(Int_t mode=kLocalCleanup) override
Turn on automatic cleanup of child frames in dtor.
Definition TGFrame.cxx:1072
void Resize(UInt_t w=0, UInt_t h=0) override
Resize the frame.
Definition TGFrame.cxx:605
void MapWindow() override
map window
Definition TGFrame.h:204
A composite frame that layout their children in horizontal way.
Definition TGFrame.h:385
Defines top level windows that interact with the system Window Manager.
Definition TGFrame.h:397
void SetWindowName(const char *name=nullptr) override
Set window name. This is typically done via the window manager.
Definition TGFrame.cxx:1788
Yield an action as soon as it is clicked.
Definition TGButton.h:228
virtual Bool_t SetTab(Int_t tabIndex, Bool_t emit=kTRUE)
Brings the composite frame with the index tabIndex to the front and generate the following event if t...
Definition TGTab.cxx:558
virtual TObjArray * GetListOfParticles() const
Definition TGenerator.h:176
Cylindrical tube class.
Definition TGeoTube.h:17
This class serves as a data storage for description of one particle.
Definition TMCParticle.h:20
Description of the dynamic properties of a particle.
Definition TParticle.h:26
virtual void SetName(const char *name)
Change (i.e.
TPythia is an interface class to F77 version of Pythia 6.2
Definition TPythia6.h:84
int GetMSTU(int i)
Definition TPythia6.h:152
void SetMDME(int i, int j, int m)
Definition TPythia6.h:184
void Initialize(const char *frame, const char *beam, const char *target, float win)
Calls PyInit with the same parameters after performing some checking, sets correct title.
Definition TPythia6.cxx:438
void SetMSEL(int m)
Definition TPythia6.h:210
void SetCKIN(int i, double c)
Definition TPythia6.h:213
void GenerateEvent() override
generate event and copy the information from /HEPEVT/ to fPrimaries
Definition TPythia6.cxx:299
void SetPMAS(int ip, int i, double m)
Definition TPythia6.h:171
void SetMSUB(int i, int m)
Definition TPythia6.h:212
void SetMSTP(int i, int m)
Definition TPythia6.h:224
void StartEmbedding(Int_t pos=kRight, Int_t subpos=-1) override
Start embedding external frame in the tab "pos" and tab element "subpos".
void SetTabTitle(const char *title, Int_t pos=kRight, Int_t subpos=-1)
Set text "title" of Tab "subpos" in TGTab "pos".
void StopEmbedding(const char *name=nullptr) override
TGTab * GetTabRight() const
Basic string class.
Definition TString.h:139
TString & ReplaceAll(const TString &s1, const TString &s2)
Definition TString.h:704
virtual const char * Getenv(const char *env)
Get environment variable.
Definition TSystem.cxx:1650
virtual int Load(const char *module, const char *entry="", Bool_t system=kFALSE)
Load a shared library.
Definition TSystem.cxx:1842
virtual const char * UnixPathName(const char *unixpathname)
Convert from a local pathname to a Unix pathname.
Definition TSystem.cxx:1050
const Double_t kR_max
const Double_t kZ_d
const Double_t kR_min
Matevz Tadel

Definition in file pythia_display.C.