TOF examples using tails

Author: preghenella

examples/cards/propagate.2kG.tails.tcl

@@ -0,0 +1,62 @@
+set barrel_Bz 0.2
+set barrel_Radius 100.e-2
+set barrel_HalfLength 200.e-2
+set barrel_TimeResolution 0.020e-9
+set barrel_TailRight 1.0
+set barrel_TailLeft 1.0
+set barrel_Acceptance { 0.0 + 1.0 * fabs(eta) < 1.443 }
+
+set ExecutionPath {
+    ParticlePropagator
+    Merger
+    Acceptance
+    DecayFilter
+    TimeSmearingTail
+    TreeWriter
+}
+
+# module Module Name
+module ParticlePropagator ParticlePropagator {
+    set InputArray Delphes/stableParticles
+    set OutputArray stableParticles
+    set ChargedHadronOutputArray chargedHadrons
+    set ElectronOutputArray electrons
+    set MuonOutputArray muons
+
+    set Bz $barrel_Bz
+    set Radius $barrel_Radius
+    set HalfLength $barrel_HalfLength
+}
+
+module Merger Merger {
+    add InputArray ParticlePropagator/chargedHadrons
+    add InputArray ParticlePropagator/electrons
+    add InputArray ParticlePropagator/muons
+    set OutputArray tracks
+}
+
+module Efficiency Acceptance {
+    add InputArray Merger/tracks
+    add OutputArray tracks
+    set EfficiencyFormula $barrel_Acceptance
+}
+
+module DecayFilter DecayFilter {
+    set InputArray Acceptance/tracks
+    set OutputArray tracks
+}
+
+module TimeSmearingTail TimeSmearingTail {
+    add InputArray DecayFilter/tracks
+    add OutputArray tracks
+    set TimeResolution $barrel_TimeResolution
+    set TailRight $barrel_TailRight
+    set TailLeft $barrel_TailLeft
+}
+
+module TreeWriter TreeWriter {
+    # add Branch InputArray BranchName BranchClass
+    add Branch Delphes/allParticles Particle GenParticle
+    add Branch TimeSmearingTail/tracks Track Track
+}
+

examples/scripts/tof.sh

@@ -6,6 +6,8 @@ NEVENTS=10000  # number of events in a run
 
 BFIELD=5.      # magnetic field  [kG]
 SIGMAT=0.020   # time resolution [ns]
+TAILLX=1.0     # tail on left    [q]
+TAILRX=1.3     # tail on right   [q]
 TOFRAD=100.    # TOF radius      [cm]
 TOFLEN=200.    # TOF half length [cm]
 TOFETA=1.443   # TOF max pseudorapidity
@@ -15,13 +17,12 @@ TOFETA=`awk -v a=$TOFRAD -v b=$TOFLEN 'BEGIN {th=atan2(a,b)*0.5; sth=sin(th); ct
 echo "maxEta = $TOFETA"
 
 ### copy relevant files in the working directory
-cp $DELPHESO2_ROOT/examples/cards/propagate.2kG.tcl propagate.tcl
+cp $DELPHESO2_ROOT/examples/cards/propagate.2kG.tails.tcl propagate.tcl
 cp $DELPHESO2_ROOT/examples/pythia8/pythia8_inel.cfg .
 cp $DELPHESO2_ROOT/examples/smearing/tof.C .
 
 ### set magnetic field
 sed -i -e "s/set barrel_Bz .*$/set barrel_Bz ${BFIELD}e\-1/" propagate.tcl
-sed -i -e "s/double Bz = .*$/double Bz = ${BFIELD}e\-1\;/" tof.C
 ### set TOF radius
 sed -i -e "s/set barrel_Radius .*$/set barrel_Radius ${TOFRAD}e\-2/" propagate.tcl
 sed -i -e "s/double tof_radius = .*$/double tof_radius = ${TOFRAD}\;/" tof.C
@@ -30,15 +31,18 @@ sed -i -e "s/set barrel_HalfLength .*$/set barrel_HalfLength ${TOFLEN}e\-2/" pro
 sed -i -e "s/double tof_length = .*$/double tof_length = ${TOFLEN}\;/" tof.C
 ### set TOF acceptance
 sed -i -e "s/set barrel_Acceptance .*$/set barrel_Acceptance \{ 0.0 + 1.0 * fabs(eta) < ${TOFETA} \}/" propagate.tcl
-### set TOF time resolution
+### set TOF time resolution and tails
 sed -i -e "s/set barrel_TimeResolution .*$/set barrel_TimeResolution ${SIGMAT}e\-9/" propagate.tcl
+sed -i -e "s/set barrel_TailRight .*$/set barrel_TailRight ${TAILRX}/" propagate.tcl
+sed -i -e "s/set barrel_TailLeft  .*$/set barrel_TailLeft ${TAILLX}/" propagate.tcl
 sed -i -e "s/double tof_sigmat = .*$/double tof_sigmat = ${SIGMAT}\;/" tof.C
 
 ### create LUTs
 BFIELDT=`awk -v a=$BFIELD 'BEGIN {print a*0.1}'`
 $DELPHESO2_ROOT/examples/scripts/create_luts.sh werner $BFIELDT $TOFRAD
 
 ### loop over runs
+rm -f .running.* delphes.*.root
 for I in $(seq 1 $NRUNS); do
 
     ### wait for a free slot