#Karoly Makonyi 12/20/14 sim_radmap - full simulation of p-p-bar event with DPM event generator usage: [path-to-the-executable]/sim_rad.exe -o [name of the outout file [blablabla.root]]\ -n [number of events to simulate [2]]\ -t [used tracker [Geant4]] PndRadMapBoxMesh - The class that is used to make the maps. One can find examples for its usage in the plot_radmap.C file plot_radmap - The code that demonstrates the abilities of the Mesh class. It produces the 'RadMap_Out.root' file where the required histos are stored. It asssumes that the simulated file (RadMap_Out_Sim.root' is in the directory where the plot_radmap is run. histodivider.C - It calculates the quantities for the required running time nominal run and - optionally - removes the low-statistic bins. The output of this code is 'figures.pdf' and 'figures.root'. TODO: add the file name as a command line parameter One have to compile it and run the resulting executable In a separate textfile (the name must be "names.dat") some details of the histos can be given: #name of the histo/directory title Xlabel Ylabel norm - name of the directory - the histos are arranged in directories one must give the name of the directory where the 'Histo' and the 'StatHisto' histograms can be found - title - title of the histogram taht will appear in the figures.* files - Xlabel - the label on the X-axis - Ylabel - the label on the Y-axis - norm - the scaling factor on the histos (1 -> quantity per event) the fields should be separated by semi-comma (;) dirmaker.sh - produces 10x10x10 disrectories (this is the directory-structure that I used [see /hera/panda/makonyi/Work/FieldSimu/SimulationMacros/Geant[3,4].15GeV]) job_submit.sh - Sends the jobs (sinulation or plot_radmap.exe)to the cluster. Presently this script is pretty dumb and can only process the files in a directory structure which was created by the dirmaker.sh. But it - runs the simulation in every directory or - runs the plot_radmap.exe in every directory - 'hadd's the resulting files comman line parameters: - -b - [b]asename -- name of the the dir where the 'dirmaker' produced the directory structure - -s - [s]imulation -- the script will run simulation in every directories in the 'basename' directory - -n - [n]umber -- number of simulated events per file - -a - [a]nalysis -- the script will run analysis - -h - [h]elp -- print out a help screen job_submit_kronos.sh - same as the job_submit.sh but runs on the kronos cluster Workflow: ---------------------------- This have to be made once ---------------- You produce files that will be processed further. If the geometry doesn't change (or obvious bug was not found) these steps should be done only once. - Make a directory where the data will be stored, for example 'mkdir [Data]'. (the 'Data' name can be changed as one prefers it) - Go to this directory: 'cd Data'. - Copy the 'dirmaker.sh' to this directory: 'cp ..Scripts/dirmaker.sh .' - Run it: './dirmaker.sh' -> This will take some time, but finally the director structures (1000 directories) will be produced - go up 'cd ..' - run the simulation 'job_submit_kronos.sh -s -b Data -n 100' this will run _100_ events in one job, and will run a separated job for every directory (1000) in the 'Data' (produced by the 'dirmaked') parallel, using the cluster ... some output will be produced see the parametrisation of the job_submit_kronos.sh script via the './job_submit_kronos.sh -h' ---------------------------- This can be made several times ---------------- After having the simulated files, several different analysis can be run on the same set of simulated files. - run the 'plot_radmap' 'job_submit_kronos.sh -a -b Data' this will run 'plot_radmap.exe' on every simulated files in the 'Data' directory parallel, using the cluster ... some output will be produced. If every simulated files are processed the results will be added into one file (RadMap_Out_Final_Geant4.15GeV.root). - be sure that all the files are added properly together in the [Data]/RadMap_Out_Fineal.root - rescale the results to one nominal year operation run the histodivider ('./histodivider')