Fork me on GitHub

source: git/README@ 9da65a5

ImprovedOutputFile Timing dual_readout llp
Last change on this file since 9da65a5 was 1d17f21, checked in by Michele Selvaggi <michele.selvaggi@…>, 9 years ago

added instructions for configuring Delphes on lxplus

  • Property mode set to 100644
File size: 3.6 KB
RevLine 
[d7d2da3]1Quick start with Delphes
2========================
3
4Commands to get the code:
5
[9d03c3b]6 wget http://cp3.irmp.ucl.ac.be/downloads/Delphes-3.2.0.tar.gz
[d7d2da3]7
[9d03c3b]8 tar -zxf Delphes-3.2.0.tar.gz
[d7d2da3]9
10Commands to compile the code:
11
[9d03c3b]12 cd Delphes-3.2.0
[d7d2da3]13
14 make
15
16Finally, we can run Delphes:
17
18 ./DelphesHepMC
19
20Command line parameters:
21
22 ./DelphesHepMC config_file output_file [input_file(s)]
23 config_file - configuration file in Tcl format
24 output_file - output file in ROOT format,
25 input_file(s) - input file(s) in HepMC format,
26 with no input_file, or when input_file is -, read standard input.
27
28Let's simulate some Z->ee events:
29
30 wget http://cp3.irmp.ucl.ac.be/downloads/z_ee.hep.gz
31 gunzip z_ee.hep.gz
[a8782e8]32 ./DelphesSTDHEP cards/delphes_card_CMS.tcl delphes_output.root z_ee.hep
[d7d2da3]33
34or
35
[a8782e8]36 curl -s http://cp3.irmp.ucl.ac.be/downloads/z_ee.hep.gz | gunzip | ./DelphesSTDHEP cards/delphes_card_CMS.tcl delphes_output.root
[d7d2da3]37
38For more detailed documentation, please visit
39
40https://cp3.irmp.ucl.ac.be/projects/delphes/wiki/WorkBook
41
42
[1d17f21]43Configure Delphes on lxplus.cern.ch
44====================================
45
46git clone git://github.com/delphes/delphes.git Delphes
47
48cd Delphes
49
50source /afs/cern.ch/sw/lcg/app/releases/ROOT/5.34.32/x86_64-slc6-gcc48-opt/root/bin/thisroot
51
52make
53
54
[d7d2da3]55Simple analysis using TTree::Draw
56=================================
57
58Now we can start ROOT and look at the data stored in the output ROOT file.
59
60Start ROOT and load Delphes shared library:
61
62 root -l
63 gSystem->Load("libDelphes");
64
65Open ROOT file and do some basic analysis using Draw or TBrowser:
66
67 TFile::Open("delphes_output.root");
68 Delphes->Draw("Electron.PT");
69 TBrowser browser;
70
71Note 1: Delphes - tree name, it can be learned e.g. from TBrowser
72
73Note 2: Electron - branch name; PT - variable (leaf) of this branch
74
75Complete description of all branches can be found in
76
77 doc/RootTreeDescription.html
78
79This information is also available at
80
81 https://cp3.irmp.ucl.ac.be/projects/delphes/wiki/WorkBook/RootTreeDescription
82
83
84Macro-based analysis
85====================
86
87Analysis macro consists of histogram booking, event loop (histogram filling),
88histogram display.
89
90Start ROOT and load Delphes shared library:
91
92 root -l
93 gSystem->Load("libDelphes");
94
95Basic analysis macro:
96
97{
98 // Create chain of root trees
99 TChain chain("Delphes");
100 chain.Add("delphes_output.root");
101
102 // Create object of class ExRootTreeReader
103 ExRootTreeReader *treeReader = new ExRootTreeReader(&chain);
104 Long64_t numberOfEntries = treeReader->GetEntries();
105
106 // Get pointers to branches used in this analysis
107 TClonesArray *branchElectron = treeReader->UseBranch("Electron");
108
109 // Book histograms
110 TH1 *histElectronPT = new TH1F("electron pt", "electron P_{T}", 50, 0.0, 100.0);
111
112 // Loop over all events
113 for(Int_t entry = 0; entry < numberOfEntries; ++entry)
114 {
115
116 // Load selected branches with data from specified event
117 treeReader->ReadEntry(entry);
118
119 // If event contains at least 1 electron
120 if(branchElectron->GetEntries() > 0)
121 {
122 // Take first electron
123 Electron *electron = (Electron*) branchElectron->At(0);
124
125 // Plot electron transverse momentum
126 histElectronPT->Fill(electron->PT);
127
128 // Print electron transverse momentum
129 cout << electron->PT << endl;
130 }
131
132 }
133
134 // Show resulting histograms
135 histElectronPT->Draw();
136}
137
138
139More advanced macro-based analysis
140==================================
141
142The 'examples' directory contains ROOT macros Example1.C, Example2.C and Example3.C.
143
144Here are the commands to run these ROOT macros:
145
146 root -l
147 .X examples/Example1.C("delphes_output.root");
148
149or
150
[21f3c04]151 root -l examples/Example1.C'("delphes_output.root")'
Note: See TracBrowser for help on using the repository browser.