"How to" for analysis code

How to Guide for analysis code

[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 process run?

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 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?

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".