cwb_setcuts


Usage

This command applies thresholds on tree leaves to select events with certain characteristics defined by the user. We remind that for each merge label we have two files type: the wave and the mdc/live, with different structure but some leafs are in common. So when choosing selection cuts, you should distinguish between the cuts that should be applied only to wave tree from the ones that will be applied also to mdc/live. One example of the former cut is a selection on rho or cc, on the contray a selection on the run is an example of the latter.

Syntax

  • cwb_setcuts mlabel/list wcuts clabel jcuts(opt) mcuts(opt) lcuts(opt)
  • mlabel/lists same label used by cwb_merge, if the user writes lists it list the available labels
  • clabel label that will be added to the new merged files names (M1 -> M1.C_clabel)
  • wcuts selections for 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
  • jcuts job selection for wave & mdc trees must contains only the leaf run the syntax is the the same as the one used for the selection TTree
  • mcuts selections for only mdc tree must contains only leaf of mdc tree (not the leaf run) the syntax is the the same as the one used for the selection TTree
  • lcuts selections of lag/slag : is applied only to wave and livetime tree must contains only lag/slag leafs the syntax is the the same as the one used for the selection TTree

Alternative Syntax

  • cwb_setcuts mlabel/list options
  • mlabel/lists same label used by cwb_merge
  • options options (options must be declared inside ‘…’) : see previous definitions
  • cwb_setcuts mlabel/list wcuts clabel jcuts(opt) mcuts(opt) lcuts(opt)

  • –label clabel default=”“

  • –tcuts tcuts default=”” : tcuts is the name of TCut selection defined in user_pparameters.C, it is added to wcuts

  • –wcuts wcuts default=”“

  • –jcuts jcuts default=”“

  • –mcuts mcuts default=”“

  • –lcuts lcuts default=”“

  • –unique true/false select from the root file a unique reconstructed event for each injected signal.

    The signal is selected from the reconstructed events which have the same injected gps time, the selected one has the maximum rho. When this option is used the label “.C_U” will be added to the new merged files names. (default=false) A new leaf is added to the waveburst tree : fsize=number of fragments found in the iwindow NOTE : this option can be applied only when simulation>0.

Note

not use !expression but !(expression)
EX : ot use !veto_cat3_H1 but !(veto_cat3_H1)
use !expression instead of “expression”
EX : use !veto_cat3_H1 instead of “!(veto_cat3_H1)”

Examples

* cwb_setcuts M1 'netcc[1]>0.8'  cc1_gt_0d8_run_lt_100 'run<100'
  - apply the cut ((netcc[1]>0.8)&&(run<100)) to wave file
  - apply the cut (run<100) to the live file

  if the input file is with name: merge/wave_S6A_BKG_L1H1_run1.M1.root
  the output file assumes the name: merge/wave_S6A_BKG_L1H1_run1.M1.C_cc1_gt_0d8_run_lt_100.root
  and the same for live file.

  Alternative Syntax :
  cwb_setcuts M1 '--wcuts netcc[1]>0.8  --label cc1_gt_0d8_run_lt_100 --jcuts run<100'

* cwb_setcuts M1.C_rho0_gt6.V_cat2LH_hvetoLH_cat3LH  '!(veto_cat3_H1) && !(veto_cat3_L1)' nohwinj

* cwb_setcuts M1 'netcc[2]>0.45'  cc2_gt_0d45_lag1_eq_0 '' '' 'lag[1]==0'
  - apply the cut ((netcc[2]>0.45)&&(lag[1]==0)) to wave file
  - apply the cut (lag[1]==0) to the live file

  Alternative Syntax :
  cwb_setcuts M1 '--wcuts netcc[2]>0.45 --label cc2_gt_0d45_lag1_eq_0 --lcuts lag[1]==0'

* cwb_setcuts M1 '--unique true'
  - apply unique selection of detected events

  if the input file is with name: merge/wave_S6A_BKG_L1H1_run1.M1.root
  the output file assumes the name: merge/wave_S6A_BKG_L1H1_run1.M1.C_U.root

* cwb_setcuts M1 '--unique true --wcuts netcc[2]>0.45  --label cc2_gt_0d45'
  - apply unique selection of detected events and tree cuts

  if the input file is with name: merge/wave_S6A_BKG_L1H1_run1.M1.root
  the output file assumes the name: merge/wave_S6A_BKG_L1H1_run1.M1.C_U.C_cc2_gt_0d45.root

* cwb_setcuts M1 '--tcuts MyTCut --wcuts netcc[2]>0.45 --label cc2_gt_0d45_mytcut'
  - select events with (MyTCut.GetTitle()) && (netcc[2]>0.45)

  The MyTCut must be defined in user_pparameters.C
  Ex : TCut MyTCut('MyTCut','rho[0]>6'); // MyTCut.GetTitle() is 'rho[0]>6'
  or it can be defined in a file.hh and included in user_pparameters.C inside the brackets {}
  Ex : #include "file.hh"