cwb_setifar¶
Usage
This command adds a new ifar (inverse false alarm rate) leaf to the
selected entries in the merged simulation root file.
For each selected entry a new leaf is created “ifar” which value is
obtained from the background FAR vs RHO file.
It is a text file which contains a corrispondence between far and rho (from background estimation)
Syntax
cwb_setifar mlabel/list tsel farfile ilabel mode
- mlabel/lists same label used by cwb_merge, if the user writes lists it list the available labels
- tsel treen selections, only wave tree
- must contains only leaf of waveburst tree (not the leaf run) the syntax is the the same as the one used for the selection TTree
- ilabel label that will be added to the new merged files names (M1 -> M1.S_ilabel)
- mode inclusive/exclusive
- inclusive : if a previous ifar is already associated to the entry then it is updated if the new ifar is greatar exclusive : the ifar is updated regardless what is is the old one
Alternative Syntax
cwb_setifar mlabel/list options
- mlabel/lists same label used by cwb_merge
- options options (options must be declared inside ‘…’) : see previous definitions
- –label ilabel default=”“
- –tsel tsel default=”” : tsel is the name or a logical expression of names of TCut selections defined in user_pparameters.C
- –file file default=”” : file is the name of the far file defined in user_pparameters.C, name can be an element of an array (Ex: name[2])
- –mode mode default=”“
- –xtsel xtsel default=”” : the tsel are explicitly defined (must be used as alternative to –tsel)
- –xfile xfile default=”” : the far file path is declared explicitly (must be used as alternative to –file)
Note
label, tsel/xtsel, file/xfile, mode are mandatory
Note
farfile is produced by the background report, and it is saved under the report data dir and its name is far_rho.txt
Note
examples of allowed logical expressions are X&&Y, (X||Y), !X, ((X||Y)&&Z)
Examples
* cwb_setifar M1 'netcc[1]>0.8' /home/waveburst/work_dir/report/postprod/far_rho.txt cc1_gt_0d8 exclusive - add ifar leaf for all entries with ((netcc[1]>0.8)) in the wave file - far is read from /home/frank/work_dir/report/postprod/far_rho.txt : is a two columns file "rho far" Alternative Syntax : cwb_setifar M1 '--xtsel netcc[1]>0.8 --label cc1_gt_0d8 \ --xfile /home/waveburst/work_dir/report/postprod/far_rho.txt --mode exclusive' or : cwb_setifar M1 '--tsel cc1_tsel --label cc1_gt_0d8 --file farfile --mode exclusive' - cc1_tsel is the name of TCut selection defined in user_pparameters.C Ex : TCut cc1_tsel("cc1_tsel","netcc[1]>0.8"); - farfile is the name of the far file defined in user_pparameters.C Ex : TString farfile = "/home/waveburst/work_dir/report/postprod/far_rho.txt"; or : cwb_setifar M1 '--tsel cc0_cut&&cc2_cut --label cc0_cc2_cut --file farfile[2] --mode exclusive' - cc0_cut and cc2_cut are the names of TCut selections defined in user_pparameters.C Ex : TCut cc0_cut("cc0_cut","netcc[0]>0.8"); TCut cc2_cut("cc2_cut","netcc[2]>0.8"); - farfile[2] is the iname of the second element of the far file array defined in user_pparameters.C Ex : TString farfile[3]; farfile[0] = "/home/waveburst/work_dir_0/report/postprod/far_rho.txt"; farfile[1] = "/home/waveburst/work_dir_1/report/postprod/far_rho.txt"; farfile[2] = "/home/waveburst/work_dir_2/report/postprod/far_rho.txt";