diff --git a/distribsForHNLandBG_byEvent.py b/distribsForHNLandBG_byEvent.py
index 37d3655..f1f3661 100755
--- a/distribsForHNLandBG_byEvent.py
+++ b/distribsForHNLandBG_byEvent.py
@@ -671,6 +671,17 @@
                          tc.ndf(25)[1], tc.doca(1.)[1], tc.redChi2(5.)[1], tc.daughtersP(1.)[1], tc.ip(250.)[1],
                          tc.notVetoedNoRPC(1000.)[1], tc.muon1(1)[1], tc.muon2(1)[1] ]
 
+        elif answers == 5:
+            # vetoes at the beginning
+            cuts =     [ tc.recoed()[0], tc.notVetoed(1000.)[0], tc.nRecoed(1)[0],
+                         tc.fiducial(zmin + 20., zmax - 20., 249., 499.)[0], tc.goodtracks()[0],
+                         tc.ndf(25)[0], tc.doca(1.)[0], tc.redChi2(5.)[0], tc.daughtersP(1.)[0], tc.ip(10.)[0],
+                         tc.muon1(1)[0], tc.muon2(1)[0] ]
+            cutnames = [ tc.recoed()[1], tc.notVetoed(1000.)[1], tc.nRecoed(1)[1],
+                         tc.fiducial(zmin + 20., zmax - 20., 249., 499.)[1], tc.goodtracks()[1],
+                         tc.ndf(25)[1], tc.doca(1.)[1], tc.redChi2(5.)[1], tc.daughtersP(1.)[1], tc.ip(10.)[1],
+                         tc.muon1(1)[1], tc.muon2(1)[1] ]
+
     #elif (study == 'nu'):
     #    cuts =     [ tc.recoed()[0], tc.notVetoed(1000.)[0], tc.redChi2(5.)[0], tc.ndf(25)[0], tc.fiducial(zmin, zmax, 250., 500.)[0], tc.goodtracks()[0], 
     #                 tc.ecal()[0], tc.muon1(1)[0], tc.muon2(1)[0], tc.doca(30.)[0], tc.ip(250.)[0] ]
@@ -763,6 +774,17 @@
                          tc.fiducial(zmin + 20., zmax - 20., 249., 499.)[1], tc.goodtracks()[1],
                          tc.ndf(25)[1], tc.doca(1.)[1], tc.redChi2(5.)[1], tc.daughtersP(1.)[1], tc.ip(250.)[1],
                          tc.notVetoedNoRPC(1000.)[1], tc.muon1(2)[1], tc.muon2(2)[1] ]
+        elif answers == 5:
+            # vetoes at the beginning
+            cuts =     [ tc.recoed()[0], tc.notVetoed(1000.)[0], tc.nRecoed(1)[0],
+                         tc.fiducial(zmin + 20., zmax - 20., 249., 499.)[0], tc.goodtracks()[0],
+                         tc.ndf(25)[0], tc.doca(1.)[0], tc.redChi2(5.)[0], tc.daughtersP(1.)[0], tc.ip(250.)[0],
+                         tc.muon1(2)[0], tc.muon2(2)[0] ]
+            cutnames = [ tc.recoed()[1], tc.notVetoed(1000.)[1], tc.nRecoed(1)[1],
+                         tc.fiducial(zmin + 20., zmax - 20., 249., 499.)[1], tc.goodtracks()[1],
+                         tc.ndf(25)[1], tc.doca(1.)[1], tc.redChi2(5.)[1], tc.daughtersP(1.)[1], tc.ip(250.)[1],
+                         tc.muon1(2)[1], tc.muon2(2)[1] ]
+
 
     elif (study == 'tiny'):
         cuts =     [ tc.recoed()[0], tc.notVetoed(1000.)[0], tc.redChi2(5.)[0], tc.ndf(25)[0], tc.fiducial(zmin, zmax, 250., 500.)[0], tc.goodtracks()[0], 
@@ -810,6 +832,15 @@
                          tc.fiducial(zmin + 20., zmax - 20., 249., 499.)[1], tc.goodtracks()[1],
                          tc.ndf(25)[1], tc.doca(1.)[1], tc.redChi2(5.)[1], tc.daughtersP(1.)[1], tc.ip(250.)[1],
                          tc.notVetoedNoRPC(1000.)[1], ]
+        elif answers == 5:
+            # vetoes at the beginning
+            cuts =     [ tc.recoed()[0], tc.notVetoed(1000.)[0], tc.nRecoed(1)[0],
+                         tc.fiducial(zmin + 20., zmax - 20., 249., 499.)[0], tc.goodtracks()[0],
+                         tc.ndf(25)[0], tc.doca(1.)[0], tc.redChi2(5.)[0], tc.daughtersP(1.)[0], tc.ip(250.)[0],]
+            cutnames = [ tc.recoed()[1], tc.notVetoed(1000.)[1], tc.nRecoed(1)[1],
+                         tc.fiducial(zmin + 20., zmax - 20., 249., 499.)[1], tc.goodtracks()[1],
+                         tc.ndf(25)[1], tc.doca(1.)[1], tc.redChi2(5.)[1], tc.daughtersP(1.)[1], tc.ip(250.)[1],]
+
 
 
     # if not cutnames:
@@ -856,7 +887,8 @@
 
 def tableWithDraw(answers=False):
     results = {}
-    for study in ['nu', 'pimu', 'fakedp']:  # ['nu', 'pimu', 'mumunu', 'tiny', 'fakedp']:
+    #for study in ['nu', 'pimu', 'fakedp']:  # ['nu', 'pimu', 'mumunu', 'tiny', 'fakedp']:
+    for study in ['nu', 'pimu', 'mumunu', 'tiny', 'fakedp']:
         results[study] = {}
         numbers, numbersWithWeights, cuts = countWithDraw(study, answers)
         results[study]['numbers'] = numbers