"How to" for analysis code: Difference between revisions
No edit summary |
|||
Line 83: | Line 83: | ||
= How to simulate the kinematics? = | = How to simulate the kinematics? = | ||
If the code is git clone or download recently, please go to ~/digios/analysis/Cleopatra, and make. | ---------------------------------- | ||
CAUTION: If the code is git clone or download recently, please go to ~/digios/analysis/Cleopatra, and make. | |||
---------------------------------- | |||
edit | edit |
Revision as of 20:21, July 11, 2019
[last updated July 11, 2019]
Everything should be performed in digios/analysis/working directory,
- .txt file extension is editable for user.
- .dat file extension is not supposed to be edited, they are generated.
How to get option for a method or function ?
For most of the bash script, if you have no argument or the argument is "-help" it will display help massage
working>./process_run
working>../Cleopatra/Transfer -help
For most root macro, for example, RDTCutsCreator(), after typed the first bracket, use tab, it will display the arguments.
root>RDTCutsCreator(
then press tab. it will show
root>RDTCutsCreator(TString filename, TString saveFileName = "rdtCuts.root", int eRange = 4000, int eRange = 60000 )
How to process run?
CAUTION: If the code is git clone or download recently, please go to ~/digios/analysis/GEBSort, and
GEBSort>make offline
working>./process_run [RUNNUM]
The bash script will download run [RUNNUM], merge, event building, mapping, and Monitors.C
to see more options, ./process_run without any argument
working>./process_run
How to process multiple runs ?
edit ChainMonitors.C
you can see it use the TChain class.
to add files, for example, add the run 33 to 40, 42, 45, and 50
chain->Add("../root_data/gen_run03[3-9].root"); chain->Add("../root_data/gen_run04[0,2,5].root"); chain->Add("../root_data/gen_run050.root");
Then, in process_run
working>./process_run 50 1 1 1 2
Notice that it will download run 50 ONLY.
How to change or add a recoil Cuts ?
see #How to use RDTCutsCreator ?
How to change the array position ?
edit the file detectorGeo.txt
How to change the calculation for Ex and thetaCM ?
edit the file reactionConfig.txt, than recalculate the reaction.dat by running ../Cleopatra/Transfer
working>../Cleopatra/Transfer
How to change the kinematic E-Z line?
see #How to simulate the kinematics?
The E-Z lines are stored in fxList in transfer.root
In Monitors.C, at the end or Monitors::Terminate().
you can see the code call fxList.
How to simulate the kinematics?
CAUTION: If the code is git clone or download recently, please go to ~/digios/analysis/Cleopatra, and make.
edit
- reactionConfig.txt
- detectorGeo.txt
- Ex.txt
Then run
working>../Cleopatra/Transfer
after that, 2 files will be generated, 1) reaction.dat, 2) transfer.root
a tree is inside transfer.root.
you can see the tree branches by
root> tree->Print()
(I know some variable names are confusing, and I will update them soon)
How to use RDTCutsCreator ?
After ./process_run, the ../Armory/RDTCutsCreator.C will be loaded.
root>RDTCutsCreator( [filename], [save file name], [E range], [dE range] )
The default [save file name] is "rdtCuts.root". If you like to use other [save file name], please correct the TString rdtCutFile in Monitors.C
if something go wrong during the RDTCutsCreator, simply delete the "rdtCuts.root".
How to do DWBA calculation ?
If the code is git clone or download recently, please go to ~/digios/analysis/Cleopatra, and make.
ONLY IN LINUX, because Ptolemy only work in Linux.
edit the file example, then
working>./Cleopatra example 1 1 1 1 0 [minAng] [maxAng] [angStep]
This will generate the example.in, feed to Ptolemy, output example.out, extract the differential cross section into
- example.Xsec.txt
- example.root
And plot the d.s.c. using example.root
It will also generate example.Ex.txt
For arguments meaning, use Cleopatra without any argument
working>./Cleopatra
How to do kinematics simulation with DWBA differential cross section ?
see #How to do DWBA calculation ? first.
edit
- reactionConfig.txt
- detectorGeo.txt
- example
Make sure the reaction in the reactionConfig.txt is matching with that in example
working>./Cleopatra example 1 1 1 0 1 0 180 1
The excitation levels is using example.Ex.txt.
The bash script will calculate the DWBA d.s.c. using Cleopatra, then using ../Cleopatra/Transfer to do the simulation.