How to run Rivet :
- setupATLAS
- asetup 21.6,latest,AthGeneration
- source setupRivet.sh
- Copy local_j0.py from /raid01/users/veena/top_charm/scratch
- local_j0.py is the script that runs over the root files using the analysis code specified in the rivet.Analyses += [ ‘ ‘ ] line in the script.
- If the analysis code is to be written from the scratch, the skeleton of the analysis code can be created using the command : rivet-mkanalysis MY_ANALYSIS where MY_ANALYSIS name of the analysis code.
- export RIVET_ANALYSIS_PATH=$PWD : to set environment variable.
- To compile the code : rivet-build RivetMY _ANALYSIS.so MY_ANALYSIS.cc
- To run the code : athena local_j0.py
- This creates a yoda file
- gzip -d <name of the output file>
- Copy convert.py from /raid01/users/veena/top_charm/scratch
- python convert.py
More details given in :
https://gitlab.cern.ch/atlas/athena/-/tree/21.6/Generators/Rivet_i
Background files that rivet code runs over is situated at
/raid01/users/veena/top_charm/scratch/bg_test_files The background file name are :
- mc16_13TeV.410470.PhPy8EG_A14_ttbar_hdamp258p75_nonallhad.merge.EVNT.e6337_e5984
- mc16_13TeV.410644.PowhegPythia8EvtGen_A14_singletop_schan_lept_top.merge.EVNT.e6527_e5984
- mc16_13TeV.410646.PowhegPythia8EvtGen_A14_Wt_DR_inclusive_top.merge.EVNT.e6552_e598
- mc16_13TeV.410658.PhPy8EG_A14_tchan_BW50_lept_top.merge.EVNT.e6671_e5984
Signal files that rivet code runs over is
/raid01/users/veena/top_charm/EVNT/output_000001_2022-09-16-1663367381/runs/1/EVNT.pool.root
Signal and Background output files produced can be located at
/raid01/users/veena/top_charm/scratch/bg_test_files/output_rootfiles
The analysis code used to produce these can be located at
/raid01/users/veena/top_charm/scratch/bg_test_files/My_analysis.cc
- My_analysis.cc ran over the background files with DSID 410470,410658,410646,410644.
- The signal files are generated using MadGraph.
Cross section [signal] is 1.908 pb
How to get grid certificate ?
How to run on the grid ?
- setupATLAS
- asetup 21.6,latest,AthGeneration
- source setupRivet.sh
- lsetup rucio
- lsetup panda
- pathena --extOutFile=Rivet.yoda.gz --inDS=mc16_13TeV.410644.PowhegPythia8EvtGen_A14_singletop_schan_lept_top.merge.EVNT.e6527_e5984 --outDS=user.vbalakri.bg410644.v1 --extFile=RivetMy_analysis.so rivet_jobOptions_grid.py
- extOutFile : the name of the output yoda file
- inDS : the name of the mc file to run over
- outDS : the name of the ouput directory [it has to be in the format , user.<username>.<your choice of the outfile name>.<version>] [everytime you resubmit the same job, the version has to be different]
- extFile : .so script use to compile the analysis code
rivet_jobOptions_grid.py : the script to run over the code [/raid01/users/veena/top_charm/scratch/rivet_jobOptions_grid.py]
To download the output file :
rucio download <name of the file > [it will be of the format .EXTO]
How to run Madgraph+Pythia in Athena and generate EVNT files ?
- Copy the jobConfig directory which contains the configuration files needed to generate the process from /raid01/users/veena/top_charm/EVNT/000001
- 000001 is the jobConfig directory and it contains the file MadGraphControl_WWyAQGC.py where the process is specified.
- Copy generation_job.sh, generation_job1.sh submit_generation.sh ,condor.submit from /raid01/users/veena/top_charm/EVNT
- To run : ./submit_generation.sh 000001 1
- An output directory like output_000001_2023-03-21-1679375124 created.
- You can find the root files inside: /raid01/users/veena/top_charm/ EVNT/output_000001_2023-03-21-1679375124/runs/
- Copy the TRexfitter directory : cp -rf /raid01/users/veena/top_charm/TRExFitter .
- cd TRexFitter
- Remove the old build directory : rm -rf build
- Make a new build directory : mkdir build && cd build
- Run cmake : cmake ../
- Compile the code with : cmake --build ./
- To run the code, after compiling, use the command: trex-fitter <action (s)> <config file>
- Refer the attachment to see the details regarding "actions"
- An example configuration file which can be edited accordingly can be found here : /raid01/users/veena/top_charm/TRExFitter/config/topcharm.config
- The input histograms can be found here : /raid01/users/veena/top_charm/TRExFitter/histo_files
- trex-fitter h topcharm.config [read input histogram files and turns them into a form suitable for further use within framework and puts them in an output directory]
- trex-fitter d topcharm.config [to produce prefit plots]
- trex-fitter w topcharm.config [to produce a workspace containing the fit model]
- trex-fitter f topcharm.config [runs a fit]
- trex-fitter p topcharm.config [runs post fit plots]
- trex-fitter l topcharm.config [to calculate the expected and observed upper limits]
- trex-fitter s topcharm.config [to calculate the observed and expected significance]
- The output directory produced after performing all these actions can be found here : /raid01/users/veena/top_charm/TRExFitter/TopCharm
- More details : https://gitlab.cern.ch/TRExStats/TRExFitter#config-file
, https://trexfitter-docs.web.cern.ch/trexfitter-docs/