Community contribution
This page contains CompHEP community contributions. Please contact us (comphep@theory.sinp.msu.ru) if you wish to put some extra modules for CompHEP on the site.
Please sign your contribution and indicate a contact information.
Program to plot CompHEP-distributions in ROOT
- imposes two or more curves on the same plot
- draws normalized distributions
- uses LATEX notations
- download: tab_plots.tgz
òÀÔ Maxim Perfilov 21/01/2007 05:35
CERN LCG MCDB Project
If you generate event samples for LHC or you are looking for some samples, probably, you will be interesting in CERN LCG MCDB project
Short Abstract about LCG MCDB: The LCG MCDB project is developed to store and document sophisticated event samples simulated for the LHC collaborations by experts. In many cases, the modern Monte-Carlo simulation of physics processes require expert knowledge in Monte-Carlo generators or significant amount of CPU time to produce an event sample. MCDB is a knowledge base mainly dedicated to accumulate such type of simulated events. The main motivation behind the LCG MCDB project is to make sophisticated MC event samples and their structured description available for various physics groups. All the data from MCDB is accessible in several convenient ways.
òÀÔ Lev Dudko 25/02/2007 05:35
Script to make Root ntuple for CompHep events
The script comphep_events2ntuple.pl allows to create Root ntuple for events, generated by CompHep, and coming in PEV type file, which name is a script argument.
Example of usage :
comphep_events2ntuple.pl Mixed.PEV
On output it gives two files :
comphep.dat -- data for ntuple doCompHepNTuple.C -- macro to be run in Root
File sizes for 11.5K ttbar events at LHC :
comphep.dat - 5.6M
CompHepData.root - 1.6M
To create ntuple run the macro in Root :
.x doCompHepNTuple.C
Documentation. Script is selfdocumented. To get help, run it with '-h' option :
comphep_events2ntuple.pl -h
Created ntuple file CompHepData.root contains tree CompHepNtuple with following structure [var_name (type) òÀÓ what] :
- IP (I) òÀÓ process ID
- CS (F) òÀÓ process CrosSec
- Qsq (F) òÀÓ Qsquared
- CStot (F) òÀÓ total CrosSec
òÀÔbeams:
- bm1KF (I) òÀÓ beam-1 KF
- bm1M (F) òÀÓ beam-1 Mass
- bm1E (F) òÀÓ beam-1 E
- bm1PDF (C) òÀÓ beam-1 PDF (combined)
- bm1PDFfile (C) òÀÓ beam-1 PDFfile
- bm2KF (I) òÀÓ beam-2 KF
- bm2M (F) òÀÓ beam-2 Mass
- bm2E (F) òÀÓ beam-2 E
- bm2PDF (C) òÀÓ beam-2 PDF (combined)
- bm2PDFfile (C) òÀÓ beam-2 PDFfile
òÀÔpartons:
- Np (I) òÀÓ process Nparton
- KF1 (I) òÀÓ KF for initial parton-1
- m1 (F) òÀÓ Mass for initial parton-1
- pz1 (F) òÀÓ Pz for initial parton-1
- KF2 (I) òÀÓ KF for initial parton-2
- m2 (F) òÀÓ Mass for initial parton-2
- pz2 (F) òÀÓ Pz for initial parton-2
òÀÔFor every outgoing parton / particle (3 - for the first) :
- KF3 (I) òÀÓ KF
- m3 (F) òÀÓ Mass
- px3 (F) òÀÓ Px
- py3 (F) òÀÓ Py
- pz3 (F) òÀÓ Pz
- pt3 (F) òÀÓ Pt
- e3 (F) òÀÓ E
- phi3 (F) òÀÓ Theta
- tet3 (F) òÀÓ Phi
- y3 (F) òÀÓ Y
- eta3 (F) òÀÓ Eta
download: comphep_events2ntuple.pl.gz (Version date: Fri Jun 5 15:28:33 MSD 2009)
Enjoy! òÀÔ Vyacheslav Rud 10/04/2008 19:16
Script to convert OLD pev format into NEW one
The script convert_comphep_old2new.pl allows to convert CompHep OLD format (cpyth-1.2.x) .pev file into NEW one (cpyth-2.0.x). Rarely, but sometime you need that. Obvious example is a needness to process CompHep generated files, stored in above mentioned CERN LCG MCDB in old format, in a framework, wich requires input CompHep events in new format.
The script accepts one argument (old.pev file) and writes out .pev file in new format (new.pev file) :
./convert_comphep_old2new.pl tq_tqb_lhc_e.pev
Output file tq_tqb_lhc_e.new.pev will be produced.
Documentation. Script is selfdocumented. To get help, run it with '-h' option :
convert_comphep_old2new.pl -h
Not all variables used in new format available in old one. So, they were predefined to some default values. Correcr them if they don't match your needs.
download: convert_comphep_old2new.pl.gz
òÀÔ Vyacheslav Rud 05/05/2008 16:07
Script for including cuts, kinematics, or regularizations into batch.dat file
When one need to calculate in CompHep complecate process with many subprocesses, he should define in CompHep GUI cuts, kinematics, and regularizations for every subprocess. These definitions, saved in results/session.dat for selected subprocess, on next step could be shared with other similar subprocesses running num_batch.pl script with -proc option. And that procedure (GUI - num_batch.pl calls) one should repeate many times, manually changing cuts, kinematics, or regularizations, because even when number of unique cuts, kinematics, and regularizations is small, they distributed between subprocesses differently and number of such combinations (and GUI - num_batch.pl calls) is not small. As a consequence, that is a long and tedious work.
That motivated me to write script tune_batch_dat.pl, which takes unique information blocks from external file and inserts them into results/batch.dat file into subprocesses, listed for each block. This external file preparation is easy - blocks format one takes from batch.dat file, cuts, kinematics, and regularizations for every subprocess you know anyway from physics.
This script could be used only after at least one CompHep GUI run with results/batch.dat file creating.
Short description.
Arguments :
1 -- batch.dat location 2 -- file with corrections
Example : ./tune_batch_dat.pl results/batch.dat batch.dat.tqb.corrections
Output files :
results/batch.dat -- corrected file results/batch.dat.orig -- saved original file
Documentation. Script is selfdocumented. To get help, run it with '-h' option :
./tune_batch_dat.pl -h
Along with some explanations this help shows sample batch.dat.tqb.corrections file.
Lines like
Cuts: 17-32 Kine: 1-2,7-8,11-14,17-18,23-24,27-30 Regs: 17-18,23-24,27-30
precede each block and contain obvious key words and list of subprocesses to which the block need to be assigned. Blocks include terminating lines (empty line for Kine, lines with equal signes for Cuts and Regs). Also used special terminators like EOKine, EOParms.
If you would like to replace for all subprocesses some parameters, which stored in batch.dat file in lines of following structure
<key> <value>
you may use Parms block (with any number of lines) :
Parms:
SQRT(S) 1.000000E+04
EOParms
This sample block sets process sqs energy to 10000 GeV. SQRT(S) will be considered as a key. Of course, one such replacement more easely could be done in text editor!
Obviously, that script works instantly.
Download: tune_batch_dat.pl.gz (Version date: Jun 30 12:42:18 MSD 2008)
òÀÔ Vyacheslav Rud 01/07/2008 15:35
Script for CompHep cross-sections energy dependence calculation
When you done CompHep calculations for the process under study at one energy, you may easely get cross-sections energy dependence running script cs_energy_dep.pl. The script also can produce plots for relative cross-sections for subprocesses.
Short description.
This script could be used only with fully prepared by CompHep results/batch.dat file.
For every energy point the script modifies results/batch.dat file and runs CompHep script num_batch.pl to calculate given process cross-sections. On output it creates Root ntuple and macros to produce appropriate cross-section plots.
Arguments :
1 -- min energy in CMS (in GeV)
2 -- max energy in CMS (in GeV)
3 -- step in energy (in GeV)
4 -- process type (string)
Example : ./cs_energy_dep.pl 10000 14000 4000 tqb
Output files :
results_vs_sqs/ -- working directory (copy of results/)
cs_vs_sqs.dat -- input file for ntuple creation
doCS_vs_sqsNTuple.C -- Root macro for ntuple creation
cs_vs_sqs_plots.C -- Root macro for CS-s plotting
cs_vs_sqs_plots.ps -- final CS plots
There is important for calculation accuracy parameter $slope which defines number of Vegas calls at every energy point - tune it for your needs (the higher Ecms the more dense grid needed for numerical calculations).
The macro cs_vs_sqs_plots.C has internal switches ifTeV, ifRelCS. To get relative cross-sections set ifRelCS to true.
Script is selfdocumented. To get help, run it with '-h' option :
./cs_energy_dep.pl -h
Enjoy!
Download: cs_energy_dep.pl.gz (Version date: Jul 8 15:07:45 MSD 2008)
Sample plots: cs_vs_sqs_plots_log_h1.ps.gz, cs_vs_sqs_plots_std_h2.ps.gz
òÀÔ Vyacheslav Rud 01/07/2008 15:55
