Index: cards/MuonCollider/MuonColliderDet_BTag_50.tcl
===================================================================
--- cards/MuonCollider/MuonColliderDet_BTag_50.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/MuonCollider/MuonColliderDet_BTag_50.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,81 @@
+    # efficiency formula for b-jets
+    add EfficiencyFormula {5} {0.5}
+
+    # default efficiency formula (misidentification rate)
+add EfficiencyFormula {0} {
+	(energy >= 500 )* ( abs(eta) <=2.66 && abs(eta) > 2.09  ) * (  6e-3 )+ \
+		(energy >= 500 )* ( abs(eta) <=2.09 && abs(eta) > 1.53  ) * (  3e-3 )+ \
+		(energy >= 500 )* ( abs(eta) <=1.53 && abs(eta)> 1.165  ) * (  2e-3 )+ \
+		(energy >= 500 )* ( abs(eta) <=1.165 && abs(eta)>0.885  ) * (  2e-3 )+ \
+		(energy >= 500 )* ( abs(eta) <=0.885 && abs(eta)>0.655  ) * (  2e-3 )+ \
+		(energy >= 500 )* ( abs(eta) <=0.655 && abs(eta)> 0.455 ) * (  2e-3 )+ \
+		(energy >= 500 )* ( abs(eta) <=0.455 && abs(eta)>0.27   ) * (  3e-3 )+ \
+		(energy >= 500 )* ( abs(eta) <=0.27 && abs(eta)> 0.09   ) * (  3e-3 )+ \
+		(energy >= 500 )* ( abs(eta) <=0.09                     ) * (  3e-3 )+ \
+		(energy < 500 && energy >= 250)* ( abs(eta) <=2.66 && abs(eta) > 2.09  ) * (  5e-3 )+ \
+		(energy < 500 && energy >= 250)* ( abs(eta) <=2.09 && abs(eta) > 1.53  ) * (  2e-3 )+ \
+		(energy < 500 && energy >= 250)* ( abs(eta) <=1.53 && abs(eta)> 1.165  ) * (  1e-3 )+ \
+		(energy < 500 && energy >= 250)* ( abs(eta) <=1.165 && abs(eta)>0.885  ) * (  1e-3 )+ \
+		(energy < 500 && energy >= 250)* ( abs(eta) <=0.885 && abs(eta)>0.655  ) * (  1e-3 )+ \
+		(energy < 500 && energy >= 250)* ( abs(eta) <=0.655 && abs(eta)> 0.455 ) * (  1e-3 )+ \
+		(energy < 500 && energy >= 250)* ( abs(eta) <=0.455 && abs(eta)>0.27   ) * (  1e-3 )+ \
+		(energy < 500 && energy >= 250)* ( abs(eta) <=0.27 && abs(eta)> 0.09   ) * (  1e-3 )+ \
+		(energy < 500 && energy >= 250)* ( abs(eta) <=0.09                     ) * (  1e-3 )+ \
+		(energy < 250 && energy >= 100)* ( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 3e-3)+ \
+		(energy < 250 && energy >= 100)* ( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 1e-3)+ \
+		(energy < 250 && energy >= 100)* ( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 9e-4)+ \
+		(energy < 250 && energy >= 100)* ( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 9e-4)+ \
+		(energy < 250 && energy >= 100)* ( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 8e-4)+ \
+		(energy < 250 && energy >= 100)* ( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 7e-4)+ \
+		(energy < 250 && energy >= 100)* ( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 7e-4)+ \
+		(energy < 250 && energy >= 100)* ( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 7e-4)+ \
+		(energy < 250 && energy >= 100)* ( abs(eta) <=0.09                     ) * ( 7e-4)+ \
+		(energy < 100 )* ( abs(eta) <=2.66 && abs(eta) > 2.09  ) * (   2e-3 )+ \
+		(energy < 100 )* ( abs(eta) <=2.09 && abs(eta) > 1.53  ) * (   7e-4 )+ \
+		(energy < 100 )* ( abs(eta) <=1.53 && abs(eta)> 1.165  ) * (   5e-4 )+ \
+		(energy < 100 )* ( abs(eta) <=1.165 && abs(eta)>0.885  ) * (   4e-4 )+ \
+		(energy < 100 )* ( abs(eta) <=0.885 && abs(eta)>0.655  ) * (   3e-4 )+ \
+		(energy < 100 )* ( abs(eta) <=0.655 && abs(eta)> 0.455 ) * (   2e-4 )+ \
+		(energy < 100 )* ( abs(eta) <=0.455 && abs(eta)>0.27   ) * (   2e-4 )+ \
+		(energy < 100 )* ( abs(eta) <=0.27 && abs(eta)> 0.09   ) * (   2e-4 )+ \
+		(energy < 100 )* ( abs(eta) <=0.09                     ) * (   3e-4 )}
+
+    # efficiency formula for c-jets (misidentification rate)
+    add EfficiencyFormula {4} {
+		(energy >= 500 )*(abs(eta) <=2.66 && abs(eta) > 2.09  ) * (  2e-2)+ \
+			(energy >= 500 )*(abs(eta) <=2.09 && abs(eta) > 1.53  ) * (  1e-2)+ \
+			(energy >= 500 )*(abs(eta) <=1.53 && abs(eta)> 1.165  ) * (  8e-3)+ \
+			(energy >= 500 )*(abs(eta) <=1.165 && abs(eta)>0.885  ) * (  8e-3)+ \
+			(energy >= 500 )*(abs(eta) <=0.885 && abs(eta)>0.655  ) * (  9e-3)+ \
+			(energy >= 500 )*(abs(eta) <=0.655 && abs(eta)> 0.455 ) * (  1e-2)+ \
+			(energy >= 500 )*(abs(eta) <=0.455 && abs(eta)>0.27   ) * (  1e-2)+ \
+			(energy >= 500 )*(abs(eta) <=0.27 && abs(eta)> 0.09   ) * (  1e-2)+ \
+			(energy >= 500 )*(abs(eta) <=0.09                     ) * (  1e-2)+ \
+			(energy < 500 && energy >= 250)*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * (  2e-2)+ \
+			(energy < 500 && energy >= 250)*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * (  6e-3)+ \
+			(energy < 500 && energy >= 250)*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * (  5e-3)+ \
+			(energy < 500 && energy >= 250)*( abs(eta) <=1.165 && abs(eta)>0.885  ) * (  5e-3)+ \
+			(energy < 500 && energy >= 250)*( abs(eta) <=0.885 && abs(eta)>0.655  ) * (  4e-3)+ \
+			(energy < 500 && energy >= 250)*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * (  4e-3)+ \
+			(energy < 500 && energy >= 250)*( abs(eta) <=0.455 && abs(eta)>0.27   ) * (  3e-3)+ \
+			(energy < 500 && energy >= 250)*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * (  4e-3)+ \
+			(energy < 500 && energy >= 250)*( abs(eta) <=0.09                     ) * (    3e-3)+ \
+			(energy < 250 && energy >= 100)*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * (    3e-2)+ \
+			(energy < 250 && energy >= 100)*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * (    4e-3)+ \
+			(energy < 250 && energy >= 100)*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * (    3e-3)+ \
+			(energy < 250 && energy >= 100)*( abs(eta) <=1.165 && abs(eta)>0.885  ) * (    2e-3)+ \
+			(energy < 250 && energy >= 100)*( abs(eta) <=0.885 && abs(eta)>0.655  ) * (    1e-3)+ \
+			(energy < 250 && energy >= 100)*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * (    1e-3)+ \
+			(energy < 250 && energy >= 100)*( abs(eta) <=0.455 && abs(eta)>0.27   ) * (    1e-3)+ \
+			(energy < 250 && energy >= 100)*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * (    1e-3)+ \
+			(energy < 250 && energy >= 100)*( abs(eta) <=0.09                     ) * (    1e-3)+ \
+			(energy < 100 )*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * (    3e-2 )+ \
+			(energy < 100 )*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * (    7e-3 )+ \
+			(energy < 100 )*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * (    3e-3 )+ \
+			(energy < 100 )*( abs(eta) <=1.165 && abs(eta)>0.885  ) * (    2e-3 )+ \
+			(energy < 100 )*( abs(eta) <=0.885 && abs(eta)>0.655  ) * (    1e-3 )+ \
+			(energy < 100 )*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * (    1e-3 )+ \
+			(energy < 100 )*( abs(eta) <=0.455 && abs(eta)>0.27   ) * (    9e-4 )+ \
+			(energy < 100 )*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * (    9e-4 )+ \
+			(energy < 100 )*( abs(eta) <=0.09                     ) * (    7e-4 )}
+
Index: cards/MuonCollider/MuonColliderDet_BTag_70.tcl
===================================================================
--- cards/MuonCollider/MuonColliderDet_BTag_70.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/MuonCollider/MuonColliderDet_BTag_70.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,81 @@
+# efficiency formula for b-jets
+add EfficiencyFormula {5} {0.7}
+# default efficiency formula (misidentification rate)
+add EfficiencyFormula {0} {(energy >= 500 )*
+    (energy >= 500 )*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 2e-2 )+ \
+			       (energy >= 500 )*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 9e-3 )+ \
+			       (energy >= 500 )*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 8e-3 )+ \
+			       (energy >= 500 )*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 8e-3 )+ \
+			       (energy >= 500 )*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 1e-2 )+ \
+			       (energy >= 500 )*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 9e-3 )+ \
+			       (energy >= 500 )*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 1e-2 )+ \
+			       (energy >= 500 )*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 1e-2 )+ \
+			       (energy >= 500 )*( abs(eta) <=0.09                     ) * ( 1e-2 )+ \
+			       (energy < 500 && energy >= 250)*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 2e-2 )+ \
+			       (energy < 500 && energy >= 250)*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 7e-3 )+ \
+			       (energy < 500 && energy >= 250)*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 6e-3 )+ \
+			       (energy < 500 && energy >= 250)*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 6e-3 )+ \
+			       (energy < 500 && energy >= 250)*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 5e-3 )+ \
+			       (energy < 500 && energy >= 250)*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 5e-3 )+ \
+			       (energy < 500 && energy >= 250)*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 5e-3 )+ \
+			       (energy < 500 && energy >= 250)*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 5e-3 )+ \
+			       (energy < 500 && energy >= 250)*( abs(eta) <=0.09                     ) * ( 5e-3 )+ \
+			       (energy < 250 && energy >= 100)*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 2e-2 )+ \
+			       (energy < 250 && energy >= 100)*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 4e-3 )+ \
+			       (energy < 250 && energy >= 100)*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 3e-3 )+ \
+			       (energy < 250 && energy >= 100)*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 3e-3 )+ \
+			       (energy < 250 && energy >= 100)*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 2e-3 )+ \
+			       (energy < 250 && energy >= 100)*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 2e-3 )+ \
+			       (energy < 250 && energy >= 100)*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 2e-3 )+ \
+			       (energy < 250 && energy >= 100)*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 2e-3 )+ \
+			       (energy < 250 && energy >= 100)*( abs(eta) <=0.09                     ) * ( 2e-3 )+ \
+			       (energy < 100 )*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 2e-2 )+ \
+			       (energy < 100 )*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 6e-3 )+ \
+			       (energy < 100 )*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 3e-3 )+ \
+			       (energy < 100 )*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 2e-3 )+ \
+			       (energy < 100 )*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 1e-3 )+ \
+			       (energy < 100 )*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 1e-3 )+ \
+			       (energy < 100 )*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 1e-3 )+ \
+			       (energy < 100 )*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 1e-3 )+ \
+			       (energy < 100 )*( abs(eta) <=0.09                     ) * ( 1e-3 )
+}
+
+# efficiency formula for c-jets (misidentification rate)
+add EfficiencyFormula {4} {
+    (energy >= 500 )*(abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 1e-1 )+ \
+	(energy >= 500 )*(abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 4e-2 )+ \
+	(energy >= 500 )*(abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 3e-2 )+ \
+	(energy >= 500 )*(abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 3e-2 )+ \
+	(energy >= 500 )*(abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 4e-2 )+ \
+	(energy >= 500 )*(abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 4e-2 )+ \
+	(energy >= 500 )*(abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 5e-2 )+ \
+	(energy >= 500 )*(abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 5e-2 )+ \
+	(energy >= 500 )*(abs(eta) <=0.09                     ) * ( 6e-2 )+ \
+	(energy < 500 && energy >= 250)*(abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 1e-1 )+ \
+	(energy < 500 && energy >= 250)*(abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 4e-2 )+ \
+	(energy < 500 && energy >= 250)*(abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 3e-2 )+ \
+	(energy < 500 && energy >= 250)*(abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 2e-2 )+ \
+	(energy < 500 && energy >= 250)*(abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 1e-2 )+ \
+	(energy < 500 && energy >= 250)*(abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 1e-2 )+ \
+	(energy < 500 && energy >= 250)*(abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 1e-2 )+ \
+	(energy < 500 && energy >= 250)*(abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 1e-2 )+ \
+	(energy < 500 && energy >= 250)*(abs(eta) <=0.09                     ) * ( 1e-2 )+ \
+	(energy < 250 && energy >= 100)*(abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 1e-1 )+ \
+	(energy < 250 && energy >= 100)*(abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 4e-2 )+ \
+	(energy < 250 && energy >= 100)*(abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 2e-2 )+ \
+	(energy < 250 && energy >= 100)*(abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 1e-2 )+ \
+	(energy < 250 && energy >= 100)*(abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 1e-2 )+ \
+	(energy < 250 && energy >= 100)*(abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 8e-3 )+ \
+	(energy < 250 && energy >= 100)*(abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 7e-3 )+ \
+	(energy < 250 && energy >= 100)*(abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 7e-3 )+ \
+	(energy < 250 && energy >= 100)*(abs(eta) <=0.09                     ) * ( 7e-3 )+ \
+	(energy < 100 )* (abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 1e-1 )+ \
+	(energy < 100 )* (abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 7e-2 )+ \
+	(energy < 100 )* (abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 4e-2 )+ \
+	(energy < 100 )* (abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 2e-2 )+ \
+	(energy < 100 )* (abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 1e-2 )+ \
+	(energy < 100 )* (abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 1e-2 )+ \
+	(energy < 100 )* (abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 1e-2 )+ \
+	(energy < 100 )* (abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 9e-3 )+ \
+	(energy < 100 )* (abs(eta) <=0.09                     ) * ( 9e-3 )
+}
Index: cards/MuonCollider/MuonColliderDet_BTag_90.tcl
===================================================================
--- cards/MuonCollider/MuonColliderDet_BTag_90.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/MuonCollider/MuonColliderDet_BTag_90.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,81 @@
+# efficiency formula for b-jets
+add EfficiencyFormula {5} {0.9}
+# default efficiency formula (misidentification rate)
+add EfficiencyFormula {0} {
+	(energy >= 500 )*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 2e-1 )+ \
+		(energy >= 500 )*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 5e-2 )+ \
+		(energy >= 500 )*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 4e-2 )+ \
+		(energy >= 500 )*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 5e-2 )+ \
+		(energy >= 500 )*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 6e-2 )+ \
+		(energy >= 500 )*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 7e-2 )+ \
+		(energy >= 500 )*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 7e-2 )+ \
+		(energy >= 500 )*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 8e-2 )+ \
+		(energy >= 500 )*( abs(eta) <=0.09                     ) * ( 8e-2 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 2e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 4e-2 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 3e-2 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 3e-2 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 2e-2 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 2e-2 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 2e-2 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 2e-2 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.09                     ) * ( 2e-2 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 2e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 5e-2 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 3e-2 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 2e-2 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 2e-2 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 1e-2 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 1e-2 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 1e-2 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.09                     ) * ( 1e-2 )+ \
+		(energy < 100 )*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 2e-1 )+ \
+		(energy < 100 )*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 8e-2 )+ \
+		(energy < 100 )*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 4e-2 )+ \
+		(energy < 100 )*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 3e-2 )+ \
+		(energy < 100 )*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 2e-2 )+ \
+		(energy < 100 )*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 1e-2 )+ \
+		(energy < 100 )*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 1e-2 )+ \
+		(energy < 100 )*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 1e-2 )+ \
+		(energy < 100 )*( abs(eta) <=0.09                     ) * ( 1e-2 )
+}
+
+# efficiency formula for c-jets (misidentification rate)
+add EfficiencyFormula {4} {
+	(energy >= 500 )*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 8e-1 )+ \
+		(energy >= 500 )*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 4e-1 )+ \
+		(energy >= 500 )*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 3e-1 )+ \
+		(energy >= 500 )*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 3e-1 )+ \
+		(energy >= 500 )*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 3e-1 )+ \
+		(energy >= 500 )*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 3e-1 )+ \
+		(energy >= 500 )*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 4e-1 )+ \
+		(energy >= 500 )*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 4e-1 )+ \
+		(energy >= 500 )*( abs(eta) <=0.09                     ) * ( 4e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 6e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 5e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 4e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 3e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 2e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 2e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 2e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 2e-1 )+ \
+		(energy < 500 && energy >= 250)*( abs(eta) <=0.09                     ) * ( 2e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 5e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 4e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 3e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 2e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 2e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 2e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 2e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 2e-1 )+ \
+		(energy < 250 && energy >= 100)*( abs(eta) <=0.09                     ) * ( 2e-1 )+ \
+		(energy < 100 )*( abs(eta) <=2.66 && abs(eta) > 2.09  ) * ( 5e-1 )+ \
+		(energy < 100 )*( abs(eta) <=2.09 && abs(eta) > 1.53  ) * ( 5e-1 )+ \
+		(energy < 100 )*( abs(eta) <=1.53 && abs(eta)> 1.165  ) * ( 4e-1 )+ \
+		(energy < 100 )*( abs(eta) <=1.165 && abs(eta)>0.885  ) * ( 3e-1 )+ \
+		(energy < 100 )*( abs(eta) <=0.885 && abs(eta)>0.655  ) * ( 3e-1 )+ \
+		(energy < 100 )*( abs(eta) <=0.655 && abs(eta)> 0.455 ) * ( 2e-1 )+ \
+		(energy < 100 )*( abs(eta) <=0.455 && abs(eta)>0.27   ) * ( 2e-1 )+ \
+		(energy < 100 )*( abs(eta) <=0.27 && abs(eta)> 0.09   ) * ( 2e-1 )+ \
+		(energy < 100 )*( abs(eta) <=0.09                     ) * ( 2e-1 )
+}
Index: cards/MuonCollider/MuonColliderDet_BTagging.tcl
===================================================================
--- cards/MuonCollider/MuonColliderDet_BTagging.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/MuonCollider/MuonColliderDet_BTagging.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,569 @@
+module BTagging BTagging_WP50_R02N2 {
+	set JetInputArray JetMomentumSmearing_VLCR02N2/JER_VLCjetsR02N2
+	set BitNumber 0
+
+	# 50% efficiency working point
+    # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
+    # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
+    # gluon's PDG code has the lowest priority
+
+	# based on CLICdp-Note-2014-002
+
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R02N2 {
+	set JetInputArray JetMomentumSmearing_VLCR02N2/JER_VLCjetsR02N2
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R02N2 {
+	set JetInputArray JetMomentumSmearing_VLCR02N2/JER_VLCjetsR02N2
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R02N3 {
+	set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R02N3 {
+	set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R02N3 {
+	set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R02N4 {
+	set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R02N4 {
+	set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R02N4 {
+	set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R02N5 {
+	set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R02N5 {
+	set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R02N5 {
+	set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R02N6 {
+	set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R02N6 {
+	set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R02N6 {
+	set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R07N2 {
+	set JetInputArray JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+
+module BTagging BTagging_WP50_R05N2 {
+	set JetInputArray JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2
+	set BitNumber 0
+
+	# 50% efficiency working point
+    # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
+    # PDG code = the highest PDG code of a quark or gluon inside DeltaR cone around jet axis
+    # gluon's PDG code has the lowest priority
+
+	# based on CLICdp-Note-2014-002
+
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R05N2 {
+	set JetInputArray JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R05N2 {
+	set JetInputArray JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R05N3 {
+	set JetInputArray JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R05N3 {
+	set JetInputArray JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R05N3 {
+	set JetInputArray JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R05N4 {
+	set JetInputArray JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R05N4 {
+	set JetInputArray JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R05N4 {
+	set JetInputArray JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R05N5 {
+	set JetInputArray JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R05N5 {
+	set JetInputArray JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R05N5 {
+	set JetInputArray JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R05N6 {
+	set JetInputArray JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R05N6 {
+	set JetInputArray JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R05N6 {
+	set JetInputArray JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R07N2 {
+	set JetInputArray JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+
+
+module BTagging BTagging_WP70_R07N2 {
+	set JetInputArray JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R07N2 {
+	set JetInputArray JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R07N3 {
+	set JetInputArray JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R07N3 {
+	set JetInputArray JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R07N3 {
+	set JetInputArray JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R07N4 {
+	set JetInputArray JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R07N4 {
+	set JetInputArray JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R07N4 {
+	set JetInputArray JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R07N5 {
+	set JetInputArray JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R07N5 {
+	set JetInputArray JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R07N5 {
+	set JetInputArray JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R07N6 {
+	set JetInputArray JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R07N6 {
+	set JetInputArray JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R07N6 {
+	set JetInputArray JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R10N2 {
+	set JetInputArray JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R10N2 {
+	set JetInputArray JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R10N2 {
+	set JetInputArray JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R10N3 {
+	set JetInputArray JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R10N3 {
+	set JetInputArray JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R10N3 {
+	set JetInputArray JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R10N4 {
+	set JetInputArray JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R10N4 {
+	set JetInputArray JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R10N4 {
+	set JetInputArray JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R10N5 {
+	set JetInputArray JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R10N5 {
+	set JetInputArray JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R10N5 {
+	set JetInputArray JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R10N6 {
+	set JetInputArray JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R10N6 {
+	set JetInputArray JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R10N6 {
+	set JetInputArray JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R12N2 {
+	set JetInputArray JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R12N2 {
+	set JetInputArray JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R12N2 {
+	set JetInputArray JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R12N3 {
+	set JetInputArray JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R12N3 {
+	set JetInputArray JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R12N3 {
+	set JetInputArray JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R12N4 {
+	set JetInputArray JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R12N4 {
+	set JetInputArray JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R12N4 {
+	set JetInputArray JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R12N5 {
+	set JetInputArray JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R12N5 {
+	set JetInputArray JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R12N5 {
+	set JetInputArray JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R12N6 {
+	set JetInputArray JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R12N6 {
+	set JetInputArray JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R12N6 {
+	set JetInputArray JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R15N2 {
+	set JetInputArray JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R15N2 {
+	set JetInputArray JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R15N2 {
+	set JetInputArray JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R15N3 {
+	set JetInputArray JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R15N3 {
+	set JetInputArray JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R15N3 {
+	set JetInputArray JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R15N4 {
+	set JetInputArray JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R15N4 {
+	set JetInputArray JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R15N4 {
+	set JetInputArray JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R15N5 {
+	set JetInputArray JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R15N5 {
+	set JetInputArray JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R15N5 {
+	set JetInputArray JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R15N6 {
+	set JetInputArray JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6
+	set BitNumber 0
+	source MuonCollider/MuonColliderDet_BTag_50.tcl
+}
+module BTagging BTagging_WP70_R15N6 {
+	set JetInputArray JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6
+	set BitNumber 1
+	source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R15N6 {
+	set JetInputArray JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6
+	set BitNumber 2
+	source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+
+########################
+# inclusive clustering
+########################
+
+module BTagging BTagging_WP50_R02_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR02_inclusive/JER_VLCjetsR02_inclusive
+ set BitNumber 0
+ source MuonCollider/MuonColliderDet_BTag_50.tcl
+ }
+module BTagging BTagging_WP70_R02_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR02_inclusive/JER_VLCjetsR02_inclusive
+ set BitNumber 1
+ source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R02_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR02_inclusive/JER_VLCjetsR02_inclusive
+ set BitNumber 2
+ source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R05_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive
+ set BitNumber 0
+ source MuonCollider/MuonColliderDet_BTag_50.tcl
+ }
+module BTagging BTagging_WP70_R05_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive
+ set BitNumber 1
+ source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R05_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive
+ set BitNumber 2
+ source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R07_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive
+ set BitNumber 0
+ source MuonCollider/MuonColliderDet_BTag_50.tcl
+ }
+module BTagging BTagging_WP70_R07_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive
+ set BitNumber 1
+ source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R07_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive
+ set BitNumber 2
+ source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R10_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive
+ set BitNumber 0
+ source MuonCollider/MuonColliderDet_BTag_50.tcl
+ }
+module BTagging BTagging_WP70_R10_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive
+ set BitNumber 1
+ source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R10_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive
+ set BitNumber 2
+ source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R12_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive
+ set BitNumber 0
+ source MuonCollider/MuonColliderDet_BTag_50.tcl
+ }
+module BTagging BTagging_WP70_R12_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive
+ set BitNumber 1
+ source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R12_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive
+ set BitNumber 2
+ source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
+module BTagging BTagging_WP50_R15_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive
+ set BitNumber 0
+ source MuonCollider/MuonColliderDet_BTag_50.tcl
+ }
+module BTagging BTagging_WP70_R15_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive
+ set BitNumber 1
+ source MuonCollider/MuonColliderDet_BTag_70.tcl
+}
+module BTagging BTagging_WP90_R15_inclusive {
+ set JetInputArray JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive
+ set BitNumber 2
+ source MuonCollider/MuonColliderDet_BTag_90.tcl
+}
Index: cards/MuonCollider/MuonColliderDet_JetFlavorAssociation.tcl
===================================================================
--- cards/MuonCollider/MuonColliderDet_JetFlavorAssociation.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/MuonCollider/MuonColliderDet_JetFlavorAssociation.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,444 @@
+module JetFlavorAssociation JetFlavorAssociation_R02N2 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR02N2/JER_VLCjetsR02N2
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R02N3 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R02N4 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R02N5 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R02N6 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+
+module JetFlavorAssociation JetFlavorAssociation_R05N2 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R05N3 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R05N4 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R05N5 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R05N6 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+
+#R07
+module JetFlavorAssociation JetFlavorAssociation_R07N2 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R07N3 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R07N4 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R07N5 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R07N6 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+
+
+#R10
+module JetFlavorAssociation JetFlavorAssociation_R10N2 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R10N3 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R10N4 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R10N5 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R10N6 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+
+
+# R 12
+module JetFlavorAssociation JetFlavorAssociation_R12N2 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R12N3 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R12N4 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R12N5 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R12N6 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+
+# R15
+module JetFlavorAssociation JetFlavorAssociation_R15N2 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R15N3 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R15N4 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R15N5 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R15N6 {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+
+module JetFlavorAssociation JetFlavorAssociation_R02_inclusive {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR02_inclusive/JER_VLCjetsR02_inclusive
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R05_inclusive {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R07_inclusive {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R10_inclusive {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R12_inclusive {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
+module JetFlavorAssociation JetFlavorAssociation_R15_inclusive {
+
+    set PartonInputArray Delphes/partons
+    set ParticleInputArray Delphes/allParticles
+    set ParticleLHEFInputArray Delphes/allParticlesLHEF
+    set JetInputArray JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive
+
+    set DeltaR 0.5
+    set PartonPTMin 1.0
+    set PartonEtaMax 2.5
+
+}
Index: cards/MuonCollider/MuonColliderDet_JetReco.tcl
===================================================================
--- cards/MuonCollider/MuonColliderDet_JetReco.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/MuonCollider/MuonColliderDet_JetReco.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,448 @@
+############
+# Jet finder VLC
+############
+#R02 N2
+module FastJetFinder FastJetFinderVLC_R02_N2 {
+    #  set InputArray Calorimeter/towers
+    set InputArray EFlowFilter/eflow
+
+    set OutputArray VLCjetsR02N2
+
+    # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
+    set NJets 2
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.2
+    set Beta 1.0
+    set Gamma 1.0
+
+    set JetPTMin 20.0
+}
+#R02 N3
+module FastJetFinder FastJetFinderVLC_R02_N3 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR02N3
+    set NJets 3
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R02 N4
+module FastJetFinder FastJetFinderVLC_R02_N4 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR02N4
+    set NJets 4
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R02 N5
+module FastJetFinder FastJetFinderVLC_R02_N5 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR02N5
+    set NJets 5
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R02 N6
+module FastJetFinder FastJetFinderVLC_R02_N6 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR02N6
+    set NJets 6
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R05 N2
+module FastJetFinder FastJetFinderVLC_R05_N2 {
+    #  set InputArray Calorimeter/towers
+    set InputArray EFlowFilter/eflow
+
+    set OutputArray VLCjetsR05N2
+
+    # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
+    set NJets 2
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.5
+    set Beta 1.0
+    set Gamma 1.0
+
+    set JetPTMin 20.0
+}
+#R05 N3
+module FastJetFinder FastJetFinderVLC_R05_N3 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR05N3
+    set NJets 3
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R05 N4
+module FastJetFinder FastJetFinderVLC_R05_N4 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR05N4
+    set NJets 4
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R05 N5
+module FastJetFinder FastJetFinderVLC_R05_N5 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR05N5
+    set NJets 5
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R05 N6
+module FastJetFinder FastJetFinderVLC_R05_N6 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR05N6
+    set NJets 6
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R07 N2
+module FastJetFinder FastJetFinderVLC_R07_N2 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR07N2
+    set NJets 2
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.7
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R07 N3
+module FastJetFinder FastJetFinderVLC_R07_N3 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR07N3
+    set NJets 3
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.7
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R07 N4
+module FastJetFinder FastJetFinderVLC_R07_N4 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR07N4
+    set NJets 4
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.7
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R07 N5
+module FastJetFinder FastJetFinderVLC_R07_N5 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR07N5
+    set NJets 5
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.7
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R07 N6
+module FastJetFinder FastJetFinderVLC_R07_N6 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR07N6
+    set NJets 6
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 0.7
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+
+#R10N2
+module FastJetFinder FastJetFinderVLC_R10_N2 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR10N2
+    set NJets 2
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.0
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R10 N3
+module FastJetFinder FastJetFinderVLC_R10_N3 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR10N3
+    set NJets 3
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.0
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R10 N4
+module FastJetFinder FastJetFinderVLC_R10_N4 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR10N4
+    set NJets 4
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.0
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R10 N5
+module FastJetFinder FastJetFinderVLC_R10_N5 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR10N5
+    set NJets 5
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.0
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R10 N6
+module FastJetFinder FastJetFinderVLC_R10_N6 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR10N6
+    set NJets 6
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.0
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+
+#R12 N2
+module FastJetFinder FastJetFinderVLC_R12_N2 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR12N2
+    set NJets 2
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R12 N3
+module FastJetFinder FastJetFinderVLC_R12_N3 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR12N3
+    set NJets 3
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R12 N4
+module FastJetFinder FastJetFinderVLC_R12_N4 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR12N4
+    set NJets 4
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R12 N5
+module FastJetFinder FastJetFinderVLC_R12_N5 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR12N5
+    set NJets 5
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R12 N6
+module FastJetFinder FastJetFinderVLC_R12_N6 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR12N6
+    set NJets 6
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+
+
+#R15 N2
+module FastJetFinder FastJetFinderVLC_R15_N2 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR15N2
+    set NJets 2
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R15 N3
+module FastJetFinder FastJetFinderVLC_R15_N3 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR15N3
+    set NJets 3
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R15 N4
+module FastJetFinder FastJetFinderVLC_R15_N4 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR15N4
+    set NJets 4
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R15 N5
+module FastJetFinder FastJetFinderVLC_R15_N5 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR15N5
+    set NJets 5
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R15 N6
+module FastJetFinder FastJetFinderVLC_R15_N6 {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR15N6
+    set NJets 6
+    set ExclusiveClustering true
+    set JetAlgorithm 9
+    set ParameterR 1.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+
+########################
+# inclusive clustering
+# as means of comparison
+########################
+#R02
+module FastJetFinder FastJetFinderVLC_R02_inclusive {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR02_inclusive
+    set ExclusiveClustering false
+    set JetAlgorithm 9
+    set ParameterR 0.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R05
+module FastJetFinder FastJetFinderVLC_R05_inclusive {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR05_inclusive
+    set ExclusiveClustering false
+    set JetAlgorithm 9
+    set ParameterR 0.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R07
+module FastJetFinder FastJetFinderVLC_R07_inclusive {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR07_inclusive
+    set ExclusiveClustering false
+    set JetAlgorithm 9
+    set ParameterR 0.7
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R10
+module FastJetFinder FastJetFinderVLC_R10_inclusive {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR10_inclusive
+    set ExclusiveClustering false
+    set JetAlgorithm 9
+    set ParameterR 1.0
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R12
+module FastJetFinder FastJetFinderVLC_R12_inclusive {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR12_inclusive
+    set ExclusiveClustering false
+    set JetAlgorithm 9
+    set ParameterR 1.2
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
+#R15
+module FastJetFinder FastJetFinderVLC_R15_inclusive {
+    set InputArray EFlowFilter/eflow
+    set OutputArray VLCjetsR15_inclusive
+    set ExclusiveClustering false
+    set JetAlgorithm 9
+    set ParameterR 1.5
+    set Beta 1.0
+    set Gamma 1.0
+    set JetPTMin 20.0
+}
Index: cards/MuonCollider/MuonColliderDet_JetSmearing.tcl
===================================================================
--- cards/MuonCollider/MuonColliderDet_JetSmearing.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/MuonCollider/MuonColliderDet_JetSmearing.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,184 @@
+module MomentumSmearing JetMomentumSmearing_VLCR02N2 {
+    set InputArray  FastJetFinderVLC_R02_N2/VLCjetsR02N2
+    set OutputArray JER_VLCjetsR02N2
+
+     # Resolution given in dpT/pT.
+    set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR02N3 {
+ set InputArray  FastJetFinderVLC_R02_N3/VLCjetsR02N3
+  set OutputArray JER_VLCjetsR02N3
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR02N4 {
+ set InputArray  FastJetFinderVLC_R02_N4/VLCjetsR02N4
+  set OutputArray JER_VLCjetsR02N4
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR02N5 {
+ set InputArray  FastJetFinderVLC_R02_N5/VLCjetsR02N5
+  set OutputArray JER_VLCjetsR02N5
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR02N6 {
+ set InputArray  FastJetFinderVLC_R02_N6/VLCjetsR02N6
+  set OutputArray JER_VLCjetsR02N6
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR02_inclusive {
+ set InputArray  FastJetFinderVLC_R02_inclusive/VLCjetsR02_inclusive
+  set OutputArray JER_VLCjetsR02_inclusive
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR05N2 {
+    set InputArray  FastJetFinderVLC_R05_N2/VLCjetsR05N2
+    set OutputArray JER_VLCjetsR05N2
+
+     # Resolution given in dpT/pT.
+    set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR05N3 {
+ set InputArray  FastJetFinderVLC_R05_N3/VLCjetsR05N3
+  set OutputArray JER_VLCjetsR05N3
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR05N4 {
+ set InputArray  FastJetFinderVLC_R05_N4/VLCjetsR05N4
+  set OutputArray JER_VLCjetsR05N4
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR05N5 {
+ set InputArray  FastJetFinderVLC_R05_N5/VLCjetsR05N5
+  set OutputArray JER_VLCjetsR05N5
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR05N6 {
+ set InputArray  FastJetFinderVLC_R05_N6/VLCjetsR05N6
+  set OutputArray JER_VLCjetsR05N6
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR05_inclusive {
+ set InputArray  FastJetFinderVLC_R05_inclusive/VLCjetsR05_inclusive
+  set OutputArray JER_VLCjetsR05_inclusive
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR07N2 {
+ set InputArray  FastJetFinderVLC_R07_N2/VLCjetsR07N2
+  set OutputArray JER_VLCjetsR07N2
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR07N3 {
+ set InputArray  FastJetFinderVLC_R07_N3/VLCjetsR07N3
+  set OutputArray JER_VLCjetsR07N3
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR07N4 {
+ set InputArray  FastJetFinderVLC_R07_N4/VLCjetsR07N4
+  set OutputArray JER_VLCjetsR07N4
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR07N5 {
+ set InputArray  FastJetFinderVLC_R07_N5/VLCjetsR07N5
+  set OutputArray JER_VLCjetsR07N5
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR07N6 {
+ set InputArray  FastJetFinderVLC_R07_N6/VLCjetsR07N6
+  set OutputArray JER_VLCjetsR07N6
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR07_inclusive {
+ set InputArray  FastJetFinderVLC_R07_inclusive/VLCjetsR07_inclusive
+  set OutputArray JER_VLCjetsR07_inclusive
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR10N2 {
+ set InputArray  FastJetFinderVLC_R10_N2/VLCjetsR10N2
+  set OutputArray JER_VLCjetsR10N2
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR10N3 {
+ set InputArray  FastJetFinderVLC_R10_N3/VLCjetsR10N3
+  set OutputArray JER_VLCjetsR10N3
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR10N4 {
+ set InputArray  FastJetFinderVLC_R10_N4/VLCjetsR10N4
+  set OutputArray JER_VLCjetsR10N4
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR10N5 {
+ set InputArray  FastJetFinderVLC_R10_N5/VLCjetsR10N5
+  set OutputArray JER_VLCjetsR10N5
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR10N6 {
+ set InputArray  FastJetFinderVLC_R10_N6/VLCjetsR10N6
+  set OutputArray JER_VLCjetsR10N6
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR10_inclusive {
+ set InputArray  FastJetFinderVLC_R10_inclusive/VLCjetsR10_inclusive
+  set OutputArray JER_VLCjetsR10_inclusive
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR12N2 {
+ set InputArray  FastJetFinderVLC_R12_N2/VLCjetsR12N2
+  set OutputArray JER_VLCjetsR12N2
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR12N3 {
+ set InputArray  FastJetFinderVLC_R12_N3/VLCjetsR12N3
+  set OutputArray JER_VLCjetsR12N3
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR12N4 {
+ set InputArray  FastJetFinderVLC_R12_N4/VLCjetsR12N4
+  set OutputArray JER_VLCjetsR12N4
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR12N5 {
+ set InputArray  FastJetFinderVLC_R12_N5/VLCjetsR12N5
+  set OutputArray JER_VLCjetsR12N5
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR12N6 {
+ set InputArray  FastJetFinderVLC_R12_N6/VLCjetsR12N6
+  set OutputArray JER_VLCjetsR12N6
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR12_inclusive {
+ set InputArray  FastJetFinderVLC_R12_inclusive/VLCjetsR12_inclusive
+  set OutputArray JER_VLCjetsR12_inclusive
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR15N2 {
+ set InputArray  FastJetFinderVLC_R15_N2/VLCjetsR15N2
+  set OutputArray JER_VLCjetsR15N2
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR15N3 {
+ set InputArray  FastJetFinderVLC_R15_N3/VLCjetsR15N3
+  set OutputArray JER_VLCjetsR15N3
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR15N4 {
+ set InputArray  FastJetFinderVLC_R15_N4/VLCjetsR15N4
+  set OutputArray JER_VLCjetsR15N4
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR15N5 {
+ set InputArray  FastJetFinderVLC_R15_N5/VLCjetsR15N5
+  set OutputArray JER_VLCjetsR15N5
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR15N6 {
+ set InputArray  FastJetFinderVLC_R15_N6/VLCjetsR15N6
+  set OutputArray JER_VLCjetsR15N6
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR15_inclusive {
+ set InputArray  FastJetFinderVLC_R15_inclusive/VLCjetsR15_inclusive
+  set OutputArray JER_VLCjetsR15_inclusive
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.02 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
Index: cards/MuonCollider/MuonColliderDet_JetSmearing_1500.tcl
===================================================================
--- cards/MuonCollider/MuonColliderDet_JetSmearing_1500.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/MuonCollider/MuonColliderDet_JetSmearing_1500.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,152 @@
+module MomentumSmearing JetMomentumSmearing_VLCR05N2 {
+    set InputArray  FastJetFinderVLC_R05_N2/VLCjetsR05N2
+    set OutputArray JER_VLCjetsR05N2
+
+     # Resolution given in dpT/pT.
+    set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }
+}
+module MomentumSmearing JetMomentumSmearing_VLCR05N3 {                                       
+ set InputArray  FastJetFinderVLC_R05_N3/VLCjetsR05N3                                        
+  set OutputArray JER_VLCjetsR05N3                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR05N4 {                                       
+ set InputArray  FastJetFinderVLC_R05_N4/VLCjetsR05N4                                        
+  set OutputArray JER_VLCjetsR05N4                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR05N5 {                                       
+ set InputArray  FastJetFinderVLC_R05_N5/VLCjetsR05N5                                        
+  set OutputArray JER_VLCjetsR05N5                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR05N6 {                                       
+ set InputArray  FastJetFinderVLC_R05_N6/VLCjetsR05N6                                        
+  set OutputArray JER_VLCjetsR05N6                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR05_inclusive {                               
+ set InputArray  FastJetFinderVLC_R05_inclusive/VLCjetsR05_inclusive                        
+  set OutputArray JER_VLCjetsR05_inclusive                                                   
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR07N2 {                                       
+ set InputArray  FastJetFinderVLC_R07_N2/VLCjetsR07N2                                        
+  set OutputArray JER_VLCjetsR07N2                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR07N3 {                                       
+ set InputArray  FastJetFinderVLC_R07_N3/VLCjetsR07N3                                        
+  set OutputArray JER_VLCjetsR07N3                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR07N4 {                                       
+ set InputArray  FastJetFinderVLC_R07_N4/VLCjetsR07N4                                        
+  set OutputArray JER_VLCjetsR07N4                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR07N5 {                                       
+ set InputArray  FastJetFinderVLC_R07_N5/VLCjetsR07N5                                        
+  set OutputArray JER_VLCjetsR07N5                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR07N6 {                                       
+ set InputArray  FastJetFinderVLC_R07_N6/VLCjetsR07N6                                        
+  set OutputArray JER_VLCjetsR07N6                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR07_inclusive {                               
+ set InputArray  FastJetFinderVLC_R07_inclusive/VLCjetsR07_inclusive                        
+  set OutputArray JER_VLCjetsR07_inclusive                                                   
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR10N2 {                                       
+ set InputArray  FastJetFinderVLC_R10_N2/VLCjetsR10N2                                        
+  set OutputArray JER_VLCjetsR10N2                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR10N3 {                                       
+ set InputArray  FastJetFinderVLC_R10_N3/VLCjetsR10N3                                        
+  set OutputArray JER_VLCjetsR10N3                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR10N4 {                                       
+ set InputArray  FastJetFinderVLC_R10_N4/VLCjetsR10N4                                        
+  set OutputArray JER_VLCjetsR10N4                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR10N5 {                                       
+ set InputArray  FastJetFinderVLC_R10_N5/VLCjetsR10N5                                        
+  set OutputArray JER_VLCjetsR10N5                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR10N6 {                                       
+ set InputArray  FastJetFinderVLC_R10_N6/VLCjetsR10N6                                        
+  set OutputArray JER_VLCjetsR10N6                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR10_inclusive {                               
+ set InputArray  FastJetFinderVLC_R10_inclusive/VLCjetsR10_inclusive                        
+  set OutputArray JER_VLCjetsR10_inclusive                                                   
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR12N2 {                                       
+ set InputArray  FastJetFinderVLC_R12_N2/VLCjetsR12N2                                        
+  set OutputArray JER_VLCjetsR12N2                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR12N3 {                                       
+ set InputArray  FastJetFinderVLC_R12_N3/VLCjetsR12N3                                        
+  set OutputArray JER_VLCjetsR12N3                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR12N4 {                                       
+ set InputArray  FastJetFinderVLC_R12_N4/VLCjetsR12N4                                        
+  set OutputArray JER_VLCjetsR12N4                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR12N5 {                                       
+ set InputArray  FastJetFinderVLC_R12_N5/VLCjetsR12N5                                        
+  set OutputArray JER_VLCjetsR12N5                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR12N6 {                                       
+ set InputArray  FastJetFinderVLC_R12_N6/VLCjetsR12N6                                        
+  set OutputArray JER_VLCjetsR12N6                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR12_inclusive {                               
+ set InputArray  FastJetFinderVLC_R12_inclusive/VLCjetsR12_inclusive                        
+  set OutputArray JER_VLCjetsR12_inclusive                                                   
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR15N2 {                                       
+ set InputArray  FastJetFinderVLC_R15_N2/VLCjetsR15N2                                        
+  set OutputArray JER_VLCjetsR15N2                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR15N3 {                                       
+ set InputArray  FastJetFinderVLC_R15_N3/VLCjetsR15N3                                        
+  set OutputArray JER_VLCjetsR15N3                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR15N4 {                                       
+ set InputArray  FastJetFinderVLC_R15_N4/VLCjetsR15N4                                        
+  set OutputArray JER_VLCjetsR15N4                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR15N5 {                                       
+ set InputArray  FastJetFinderVLC_R15_N5/VLCjetsR15N5                                        
+  set OutputArray JER_VLCjetsR15N5                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR15N6 {                                       
+ set InputArray  FastJetFinderVLC_R15_N6/VLCjetsR15N6                                        
+  set OutputArray JER_VLCjetsR15N6                                                           
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                                                            
+module MomentumSmearing JetMomentumSmearing_VLCR15_inclusive {                               
+ set InputArray  FastJetFinderVLC_R15_inclusive/VLCjetsR15_inclusive                        
+  set OutputArray JER_VLCjetsR15_inclusive                                                   
+  set ResolutionFormula { (abs(eta) < 0.76) * ( 0.01 ) + (abs(eta) >= 0.76 ) * (0.05) }                                                              
+}                                                 
Index: cards/MuonCollider/MuonColliderDet_TauTagging.tcl
===================================================================
--- cards/MuonCollider/MuonColliderDet_TauTagging.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/MuonCollider/MuonColliderDet_TauTagging.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,553 @@
+module TauTagging TauTagging_R02N2 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+    set JetInputArray JetMomentumSmearing_VLCR02N2/JER_VLCjetsR02N2
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+    # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
+    # default efficiency formula (misidentification rate)
+    add EfficiencyFormula {0} {0.02}
+    add EfficiencyFormula {11} {0.001}
+    # efficiency formula for tau-jets
+
+    add EfficiencyFormula {15} {
+	(pt < 5) * (0.0) +
+	(pt >=5) * (0.80)
+    }
+}
+
+module TauTagging TauTagging_R02N3 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR02N3/JER_VLCjetsR02N3
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R02N4 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR02N4/JER_VLCjetsR02N4
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R02N5 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR02N5/JER_VLCjetsR02N5
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R02N6 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR02N6/JER_VLCjetsR02N6
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+
+module TauTagging TauTagging_R05N2 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+    set JetInputArray JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+    # add EfficiencyFormula {abs(PDG code)} {efficiency formula as a function of eta and pt}
+    # default efficiency formula (misidentification rate)
+    add EfficiencyFormula {0} {0.02}
+    add EfficiencyFormula {11} {0.001}
+    # efficiency formula for tau-jets
+
+    add EfficiencyFormula {15} {
+	(pt < 5) * (0.0) +
+	(pt >=5) * (0.80)
+    }
+}
+
+module TauTagging TauTagging_R05N3 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R05N4 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R05N5 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R05N6 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R07N2 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R07N3 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R07N4 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R07N5 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R07N6 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R10N2 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R10N3 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R10N4 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R10N5 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R10N6 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R12N2 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R12N3 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R12N4 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R12N5 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R12N6 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R15N2 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R15N3 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R15N4 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R15N5 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+module TauTagging TauTagging_R15N6 {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+
+ module TauTagging TauTagging_R02_inclusive {
+  set ParticleInputArray Delphes/allParticles
+  set PartonInputArray Delphes/partons
+  set JetInputArray JetMomentumSmearing_VLCR02_inclusive/JER_VLCjetsR02_inclusive
+  set DeltaR 0.5
+  set TauPTMin 1.0
+  set TauEtaMax 2.5
+  add EfficiencyFormula {0} {0.02}
+  add EfficiencyFormula {11} {0.001}
+  add EfficiencyFormula {15} {
+  (pt < 10) * (0.0) +
+  (pt >=10) * (0.80)
+  }
+  }
+
+  module TauTagging TauTagging_R05_inclusive {
+   set ParticleInputArray Delphes/allParticles
+   set PartonInputArray Delphes/partons
+   set JetInputArray JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive
+   set DeltaR 0.5
+   set TauPTMin 1.0
+   set TauEtaMax 2.5
+   add EfficiencyFormula {0} {0.02}
+   add EfficiencyFormula {11} {0.001}
+   add EfficiencyFormula {15} {
+   (pt < 10) * (0.0) +
+   (pt >=10) * (0.80)
+   }
+   }
+
+
+module TauTagging TauTagging_R07_inclusive {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+
+module TauTagging TauTagging_R10_inclusive {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+
+module TauTagging TauTagging_R12_inclusive {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
+
+
+module TauTagging TauTagging_R15_inclusive {
+ set ParticleInputArray Delphes/allParticles
+ set PartonInputArray Delphes/partons
+ set JetInputArray JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive
+ set DeltaR 0.5
+ set TauPTMin 1.0
+ set TauEtaMax 2.5
+ add EfficiencyFormula {0} {0.02}
+ add EfficiencyFormula {11} {0.001}
+ add EfficiencyFormula {15} {
+ (pt < 10) * (0.0) +
+ (pt >=10) * (0.80)
+ }
+ }
Index: cards/delphes_card_MuonColliderDet.tcl
===================================================================
--- cards/delphes_card_MuonColliderDet.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
+++ cards/delphes_card_MuonColliderDet.tcl	(revision a02a49e880f7fdc400f6a689cdce26889ecc510a)
@@ -0,0 +1,1201 @@
+#######################################
+# Muon Collider Detector TARGET model
+#
+# Michele Selvaggi michele.selvaggi@cern.ch
+# Ulrike Schnoor ulrike.schnoor@cern.ch
+#
+#
+# !!! DISCLAIMER !!!
+#
+# The parameterisation of the Muon Collider
+# has to be intended as a target performance.
+# This has not been validated by full simulation.
+# Very similar to CLIC performance, added DenseTrackFilter
+#
+# Comments:
+#   - fix: angle param in DenseTrackFilter
+#   - fix: add FWD muon collection
+#   - what to do with Eta acceptance?
+#       - for now everything reduced at 2.5
+#   - add jet substructure to valencia jets ?
+#   - added R02 jets
+#   - added electron misId for taus, and make flat eff
+#
+#######################################
+
+#######################################
+# Order of execution of various modules
+#######################################
+
+set ExecutionPath {
+    ParticlePropagator
+    TrackMergerProp
+
+    DenseProp
+    DenseMergeTracks
+    DenseTrackFilter
+
+    ChargedHadronTrackingEfficiency
+    ElectronTrackingEfficiency
+    MuonTrackingEfficiency
+
+    ChargedHadronMomentumSmearing
+    ElectronMomentumSmearing
+    MuonMomentumSmearing
+
+    TrackMerger
+
+    ECal
+    HCal
+
+    Calorimeter
+    EFlowMerger
+    EFlowFilter
+
+    PhotonEfficiency
+    PhotonIsolation
+
+    ElectronFilter
+    ElectronEfficiency
+    ElectronIsolation
+
+    ChargedHadronFilter
+
+    MuonEfficiency
+    MuonIsolation
+
+    UniqueObjectFinder
+
+    NeutrinoFilter
+    GenJetFinder
+
+
+    FastJetFinderKt
+    FastJetFinderVLC_R02_N2
+    FastJetFinderVLC_R02_N3
+    FastJetFinderVLC_R02_N4
+    FastJetFinderVLC_R02_N5
+    FastJetFinderVLC_R02_N6
+    FastJetFinderVLC_R05_N2
+    FastJetFinderVLC_R05_N3
+    FastJetFinderVLC_R05_N4
+    FastJetFinderVLC_R05_N5
+    FastJetFinderVLC_R05_N6
+    FastJetFinderVLC_R07_N2
+    FastJetFinderVLC_R07_N3
+    FastJetFinderVLC_R07_N4
+    FastJetFinderVLC_R07_N5
+    FastJetFinderVLC_R07_N6
+    FastJetFinderVLC_R10_N2
+    FastJetFinderVLC_R10_N3
+    FastJetFinderVLC_R10_N4
+    FastJetFinderVLC_R10_N5
+    FastJetFinderVLC_R10_N6
+    FastJetFinderVLC_R12_N2
+    FastJetFinderVLC_R12_N3
+    FastJetFinderVLC_R12_N4
+    FastJetFinderVLC_R12_N5
+    FastJetFinderVLC_R12_N6
+    FastJetFinderVLC_R15_N2
+    FastJetFinderVLC_R15_N3
+    FastJetFinderVLC_R15_N4
+    FastJetFinderVLC_R15_N5
+    FastJetFinderVLC_R15_N6
+
+    FastJetFinderVLC_R05_inclusive
+    FastJetFinderVLC_R02_inclusive
+    FastJetFinderVLC_R07_inclusive
+    FastJetFinderVLC_R10_inclusive
+    FastJetFinderVLC_R12_inclusive
+    FastJetFinderVLC_R15_inclusive
+
+    MissingET
+    GenMissingET
+
+    JetMomentumSmearing_VLCR02N2
+    JetMomentumSmearing_VLCR02N3
+    JetMomentumSmearing_VLCR02N4
+    JetMomentumSmearing_VLCR02N5
+    JetMomentumSmearing_VLCR02N6
+    JetMomentumSmearing_VLCR02_inclusive
+    JetMomentumSmearing_VLCR05N2
+    JetMomentumSmearing_VLCR05N3
+    JetMomentumSmearing_VLCR05N4
+    JetMomentumSmearing_VLCR05N5
+    JetMomentumSmearing_VLCR05N6
+    JetMomentumSmearing_VLCR05_inclusive
+    JetMomentumSmearing_VLCR07N2
+    JetMomentumSmearing_VLCR07N3
+    JetMomentumSmearing_VLCR07N4
+    JetMomentumSmearing_VLCR07N5
+    JetMomentumSmearing_VLCR07N6
+    JetMomentumSmearing_VLCR07_inclusive
+    JetMomentumSmearing_VLCR10N2
+    JetMomentumSmearing_VLCR10N3
+    JetMomentumSmearing_VLCR10N4
+    JetMomentumSmearing_VLCR10N5
+    JetMomentumSmearing_VLCR10N6
+    JetMomentumSmearing_VLCR10_inclusive
+    JetMomentumSmearing_VLCR12N2
+    JetMomentumSmearing_VLCR12N3
+    JetMomentumSmearing_VLCR12N4
+    JetMomentumSmearing_VLCR12N5
+    JetMomentumSmearing_VLCR12N6
+    JetMomentumSmearing_VLCR12_inclusive
+    JetMomentumSmearing_VLCR15N2
+    JetMomentumSmearing_VLCR15N3
+    JetMomentumSmearing_VLCR15N4
+    JetMomentumSmearing_VLCR15N5
+    JetMomentumSmearing_VLCR15N6
+    JetMomentumSmearing_VLCR15_inclusive
+
+
+    JetFlavorAssociation_R02N2
+    JetFlavorAssociation_R02N3
+    JetFlavorAssociation_R02N4
+    JetFlavorAssociation_R02N5
+    JetFlavorAssociation_R02N6
+    JetFlavorAssociation_R05N2
+    JetFlavorAssociation_R05N3
+    JetFlavorAssociation_R05N4
+    JetFlavorAssociation_R05N5
+    JetFlavorAssociation_R05N6
+
+    JetFlavorAssociation_R07N2
+    JetFlavorAssociation_R07N3
+    JetFlavorAssociation_R07N4
+    JetFlavorAssociation_R07N5
+    JetFlavorAssociation_R07N6
+
+    JetFlavorAssociation_R10N2
+    JetFlavorAssociation_R10N3
+    JetFlavorAssociation_R10N4
+    JetFlavorAssociation_R10N5
+    JetFlavorAssociation_R10N6
+
+    JetFlavorAssociation_R12N2
+    JetFlavorAssociation_R12N3
+    JetFlavorAssociation_R12N4
+    JetFlavorAssociation_R12N5
+    JetFlavorAssociation_R12N6
+
+    JetFlavorAssociation_R15N2
+    JetFlavorAssociation_R15N3
+    JetFlavorAssociation_R15N4
+    JetFlavorAssociation_R15N5
+    JetFlavorAssociation_R15N6
+
+    JetFlavorAssociation_R02_inclusive
+    JetFlavorAssociation_R05_inclusive
+    JetFlavorAssociation_R07_inclusive
+    JetFlavorAssociation_R10_inclusive
+    JetFlavorAssociation_R12_inclusive
+    JetFlavorAssociation_R15_inclusive
+
+
+    BTagging_WP50_R02N2
+    BTagging_WP70_R02N2
+    BTagging_WP90_R02N2
+    BTagging_WP50_R02N3
+    BTagging_WP70_R02N3
+    BTagging_WP90_R02N3
+    BTagging_WP50_R02N4
+    BTagging_WP70_R02N4
+    BTagging_WP90_R02N4
+    BTagging_WP50_R02N5
+    BTagging_WP70_R02N5
+    BTagging_WP90_R02N5
+    BTagging_WP50_R02N6
+    BTagging_WP70_R02N6
+    BTagging_WP90_R02N6
+    BTagging_WP50_R05N2
+    BTagging_WP70_R05N2
+    BTagging_WP90_R05N2
+    BTagging_WP50_R05N3
+    BTagging_WP70_R05N3
+    BTagging_WP90_R05N3
+    BTagging_WP50_R05N4
+    BTagging_WP70_R05N4
+    BTagging_WP90_R05N4
+    BTagging_WP50_R05N5
+    BTagging_WP70_R05N5
+    BTagging_WP90_R05N5
+    BTagging_WP50_R05N6
+    BTagging_WP70_R05N6
+    BTagging_WP90_R05N6
+    BTagging_WP50_R07N2
+    BTagging_WP70_R07N2
+    BTagging_WP90_R07N2
+    BTagging_WP50_R07N3
+    BTagging_WP70_R07N3
+    BTagging_WP90_R07N3
+    BTagging_WP50_R07N4
+    BTagging_WP70_R07N4
+    BTagging_WP90_R07N4
+    BTagging_WP50_R07N5
+    BTagging_WP70_R07N5
+    BTagging_WP90_R07N5
+    BTagging_WP50_R07N6
+    BTagging_WP70_R07N6
+    BTagging_WP90_R07N6
+    BTagging_WP50_R10N2
+    BTagging_WP70_R10N2
+    BTagging_WP90_R10N2
+    BTagging_WP50_R10N3
+    BTagging_WP70_R10N3
+    BTagging_WP90_R10N3
+    BTagging_WP50_R10N4
+    BTagging_WP70_R10N4
+    BTagging_WP90_R10N4
+    BTagging_WP50_R10N5
+    BTagging_WP70_R10N5
+    BTagging_WP90_R10N5
+    BTagging_WP50_R10N6
+    BTagging_WP70_R10N6
+    BTagging_WP90_R10N6
+    BTagging_WP50_R12N2
+    BTagging_WP70_R12N2
+    BTagging_WP90_R12N2
+    BTagging_WP50_R12N3
+    BTagging_WP70_R12N3
+    BTagging_WP90_R12N3
+    BTagging_WP50_R12N4
+    BTagging_WP70_R12N4
+    BTagging_WP90_R12N4
+    BTagging_WP50_R12N5
+    BTagging_WP70_R12N5
+    BTagging_WP90_R12N5
+    BTagging_WP50_R12N6
+    BTagging_WP70_R12N6
+    BTagging_WP90_R12N6
+    BTagging_WP50_R15N2
+    BTagging_WP70_R15N2
+    BTagging_WP90_R15N2
+    BTagging_WP50_R15N3
+    BTagging_WP70_R15N3
+    BTagging_WP90_R15N3
+    BTagging_WP50_R15N4
+    BTagging_WP70_R15N4
+    BTagging_WP90_R15N4
+    BTagging_WP50_R15N5
+    BTagging_WP70_R15N5
+    BTagging_WP90_R15N5
+    BTagging_WP50_R15N6
+    BTagging_WP70_R15N6
+    BTagging_WP90_R15N6
+    BTagging_WP50_R02_inclusive
+    BTagging_WP70_R02_inclusive
+    BTagging_WP90_R02_inclusive
+    BTagging_WP50_R05_inclusive
+    BTagging_WP70_R05_inclusive
+    BTagging_WP90_R05_inclusive
+    BTagging_WP50_R07_inclusive
+    BTagging_WP70_R07_inclusive
+    BTagging_WP90_R07_inclusive
+    BTagging_WP50_R10_inclusive
+    BTagging_WP70_R10_inclusive
+    BTagging_WP90_R10_inclusive
+    BTagging_WP50_R12_inclusive
+    BTagging_WP70_R12_inclusive
+    BTagging_WP90_R12_inclusive
+    BTagging_WP50_R15_inclusive
+    BTagging_WP70_R15_inclusive
+    BTagging_WP90_R15_inclusive
+
+
+    TauTagging_R02N2
+    TauTagging_R02N3
+    TauTagging_R02N4
+    TauTagging_R02N5
+    TauTagging_R02N6
+    TauTagging_R05N2
+    TauTagging_R05N3
+    TauTagging_R05N4
+    TauTagging_R05N5
+    TauTagging_R05N6
+    TauTagging_R07N2
+    TauTagging_R07N3
+    TauTagging_R07N4
+    TauTagging_R07N5
+    TauTagging_R07N6
+    TauTagging_R10N2
+    TauTagging_R10N3
+    TauTagging_R10N4
+    TauTagging_R10N5
+    TauTagging_R10N6
+    TauTagging_R12N2
+    TauTagging_R12N3
+    TauTagging_R12N4
+    TauTagging_R12N5
+    TauTagging_R12N6
+    TauTagging_R15N2
+    TauTagging_R15N3
+    TauTagging_R15N4
+    TauTagging_R15N5
+    TauTagging_R15N6
+    TauTagging_R02_inclusive
+    TauTagging_R05_inclusive
+    TauTagging_R07_inclusive
+    TauTagging_R10_inclusive
+    TauTagging_R12_inclusive
+    TauTagging_R15_inclusive
+
+    ScalarHT
+
+    TreeWriter
+}
+
+#################################
+# Propagate particles in cylinder
+#################################
+
+module ParticlePropagator ParticlePropagator {
+    set InputArray Delphes/stableParticles
+
+    set OutputArray stableParticles
+    set ChargedHadronOutputArray chargedHadrons
+    set ElectronOutputArray electrons
+    set MuonOutputArray muons
+
+    # radius of the magnetic field coverage in the calorimeter, in m
+    set Radius 1.5
+    # half-length of the magnetic field coverage in the calorimeter, in m
+    set HalfLength 2.31
+
+    # magnetic field, in T
+    set Bz 4.0
+}
+
+
+##############
+# Track merger
+##############
+
+module Merger TrackMergerProp {
+# add InputArray InputArray
+  add InputArray ParticlePropagator/chargedHadrons
+  add InputArray ParticlePropagator/electrons
+  add InputArray ParticlePropagator/muons
+  set OutputArray tracks
+}
+
+
+
+####################################
+# Track propagation to pseudo-pixel
+####################################
+
+module ParticlePropagator DenseProp {
+
+  set InputArray TrackMergerProp/tracks
+
+  # radius of the magnetic field coverage, in m
+  set Radius 0.45
+  set RadiusMax 1.5
+  # half-length of the magnetic field coverage, in m
+  set HalfLength 0.8
+  set HalfLengthMax 2.31
+
+  # magnetic field
+  set Bz 4.0
+}
+
+#####################
+# Dense Track merger
+#####################
+
+module Merger DenseMergeTracks {
+# add InputArray InputArray
+  add InputArray DenseProp/chargedHadrons
+  add InputArray DenseProp/electrons
+  add InputArray DenseProp/muons
+  set OutputArray tracks
+}
+
+
+######################
+#   Dense Track Filter
+######################
+
+module DenseTrackFilter DenseTrackFilter {
+
+  set TrackInputArray DenseMergeTracks/tracks
+
+  set TrackOutputArray tracks
+  set ChargedHadronOutputArray chargedHadrons
+  set ElectronOutputArray electrons
+  set MuonOutputArray muons
+
+  set EtaPhiRes 0.01
+  set EtaMax 2.5
+
+  set pi [expr {acos(-1)}]
+
+  set nbins_phi [expr {$pi/$EtaPhiRes} ]
+  set nbins_phi [expr {int($nbins_phi)} ]
+
+  set PhiBins {}
+  for {set i -$nbins_phi} {$i <= $nbins_phi} {incr i} {
+    add PhiBins [expr {$i * $pi/$nbins_phi}]
+  }
+
+  set nbins_eta [expr {$EtaMax/$EtaPhiRes} ]
+  set nbins_eta [expr {int($nbins_eta)} ]
+
+  for {set i -$nbins_eta} {$i <= $nbins_eta} {incr i} {
+    set eta [expr {$i * $EtaPhiRes}]
+    add EtaPhiBins $eta $PhiBins
+  }
+}
+
+
+
+
+####################################
+# Charged hadron tracking efficiency
+####################################
+
+module Efficiency ChargedHadronTrackingEfficiency {
+    set InputArray DenseTrackFilter/chargedHadrons
+    set OutputArray chargedHadrons
+    # tracking efficiency formula for charged hadrons
+
+    set EfficiencyFormula { (pt <= 0.5) * (0.00) +
+    (abs(eta) <= 2.0) * (pt > 0.5 && pt <= 1) * (0.90) +
+    (abs(eta) <= 2.0) * (pt > 1) * (0.95) +
+    (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 0.5 && pt <= 1) * (0.80) +
+    (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 1.0) * (0.85) +
+    (abs(eta) > 2.5 ) * (0.00)
+   }
+}
+
+##############################
+# Electron tracking efficiency
+##############################
+
+module Efficiency ElectronTrackingEfficiency {
+    set InputArray DenseTrackFilter/electrons
+    set OutputArray electrons
+
+    set EfficiencyFormula { (pt <= 0.5) * (0.00) +
+    (abs(eta) <= 2.0) * (pt > 0.5 && pt <= 1) * (0.90) +
+    (abs(eta) <= 2.0) * (pt > 1) * (0.95) +
+    (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 0.5 && pt <= 1) * (0.80) +
+    (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 1.0) * (0.85) +
+    (abs(eta) > 2.5 ) * (0.00)
+   }
+}
+
+##########################
+# Muon tracking efficiency
+##########################
+
+module Efficiency MuonTrackingEfficiency {
+    set InputArray DenseTrackFilter/muons
+    set OutputArray muons
+
+    set EfficiencyFormula { (pt <= 0.5) * (0.00) +
+    (abs(eta) <= 2.0) * (pt > 0.5 && pt <= 1) * (0.95) +
+    (abs(eta) <= 2.0) * (pt > 1) * (0.99) +
+    (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 0.5 && pt <= 1) * (0.90) +
+    (abs(eta) > 2.0 && abs(eta) < 2.5) * (pt > 1.0) * (0.95) +
+    (abs(eta) > 2.5 ) * (0.00)
+   }
+}
+
+########################################
+# Momentum resolution for charged tracks
+########################################
+
+module MomentumSmearing ChargedHadronMomentumSmearing {
+    set InputArray ChargedHadronTrackingEfficiency/chargedHadrons
+    set OutputArray chargedHadrons
+
+
+    # Resolution given in dpT/pT.
+    # CLICdet internal studies
+    set ResolutionFormula {
+	(abs(eta) < 2.50 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2    ) +
+      	(abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.56036e-05^2 * pt^2 +0.0148987^2    ) +
+	(abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.12382e-05^2 * pt^2 +0.00391722^2   ) +
+	(abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.16768e-05^2 * pt^2 +0.00255204^2    ) +
+	(abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.28327e-05^2 * pt^2 +0.00220587^2   ) +
+	(abs(eta) < 0.18)                      * sqrt( 1.32845e-05^2 * pt^2 +0.00209325^2   )
+
+    }
+}
+
+###################################
+# Momentum resolution for electrons
+###################################
+
+module MomentumSmearing ElectronMomentumSmearing {
+    set InputArray ElectronTrackingEfficiency/electrons
+    set OutputArray electrons
+
+    # Resolution given in dpT/pT.
+    # CLICdet internal studies
+    set ResolutionFormula {
+	(abs(eta) < 2.50 && abs(eta) >= 1.74 ) * 2 * sqrt( 8.62283e-05^2 * pt^2  + 0.0177556^2   ) +
+	(abs(eta) < 1.74 && abs(eta) >= 1.32 ) * 1.5 * sqrt( 8.62283e-05^2 * pt^2  + 0.0177556^2   ) +
+	(abs(eta) < 1.32 && abs(eta) >= 0.76 ) * sqrt( 1.0915e-05 ^2 * pt^2  + 0.00663766^2  ) +
+	(abs(eta) < 0.76 && abs(eta) >= 0.36 ) * sqrt( 1.15518e-05^2 * pt^2  + 0.00398644^2  ) +
+	(abs(eta) < 0.36 && abs(eta) >= 0.18 ) * sqrt( 1.3307e-05 ^2 * pt^2  + 0.00317807^2  ) +
+	(abs(eta) < 0.18)                      * sqrt( 1.40722e-05^2 * pt^2 + 0.00292138^2  )
+
+    }
+}
+
+###############################
+# Momentum resolution for muons
+###############################
+
+module MomentumSmearing MuonMomentumSmearing {
+    set InputArray MuonTrackingEfficiency/muons
+    set OutputArray muons
+
+    # Resolution given in dpT/pT.
+
+    # CLICdet internal studies
+    set ResolutionFormula {
+
+	(abs(eta) < 2.50 && abs(eta) >= 1.74 ) 	 * 2 * sqrt(4.57439e-05^2 * pt^2   + 0.0149328^2	   ) +
+	(abs(eta) < 1.74 && abs(eta) >= 1.32 )   * 1.5 * sqrt(4.57439e-05^2 * pt^2   + 0.0149328^2	   ) +
+	(abs(eta) < 1.32 && abs(eta) >= 0.76 )	 * sqrt(9.81626e-06^2 * pt^2   + 0.00379895^2  ) +
+	(abs(eta) < 0.76 && abs(eta) >= 0.36 )	 * sqrt(1.1959e-05^2 * pt^2   +  0.00242417^2 ) +
+	(abs(eta) < 0.36 && abs(eta) >= 0.18 )	 * sqrt(1.20149e-05^2 * pt^2  + 0.00219291^2  ) +
+	(abs(eta) < 0.18)                     	 * sqrt(1.29686e-05^2 * pt^2  + 0.0020392^2      )
+
+    }
+}
+
+##############
+# Track merger
+##############
+
+module Merger TrackMerger {
+    # add InputArray InputArray
+    add InputArray ChargedHadronMomentumSmearing/chargedHadrons
+    add InputArray ElectronMomentumSmearing/electrons
+    add InputArray MuonMomentumSmearing/muons
+    set OutputArray tracks
+}
+
+#############
+#   ECAL
+#############
+
+module SimpleCalorimeter ECal {
+    set ParticleInputArray ParticlePropagator/stableParticles
+    set TrackInputArray TrackMerger/tracks
+
+    set TowerOutputArray ecalTowers
+    set EFlowTrackOutputArray eflowTracks
+    set EFlowTowerOutputArray eflowPhotons
+
+    set IsEcal true
+
+    set EnergyMin 0.5
+    set EnergySignificanceMin 1.0
+
+    set SmearTowerCenter true
+
+    set pi [expr {acos(-1)}]
+
+    # lists of the edges of each tower in eta and phi
+    # each list starts with the lower edge of the first tower
+    # the list ends with the higher edged of the last tower
+
+    #ECAL barrel: dphi = 0.2 degree, deta=0.003 towers up to |eta| <=1.2
+    #ECAL endcaps: dphi = 0.8 degree, deta=0.02 towers up to |eta| <=2.5
+
+    #barrel:
+    #dphi = 0.2 degree towers up to eta <=1.2
+    set PhiBins {}
+    for {set i -900} {$i <= 900} {incr i} {
+	add PhiBins [expr {$i * $pi/900.0 }]
+    }
+    # 0.003 unit (5x5 mm^2) in eta up to eta <=1.2
+    for {set i -400} {$i <=400} {incr i} {
+	set eta [expr {$i * 0.003}]
+	add EtaPhiBins $eta $PhiBins
+    }
+
+    #endcaps:
+    #dphi = 0.8 degree towers for 1.2 < eta <=2.5
+    set PhiBins {}
+    for {set i -225} {$i <= 225} {incr i} {
+	add PhiBins [expr {$i * $pi/225.}]
+    }
+    #deta=0.02 units for 1.2 < |eta| <=2.5
+    #first, from -2.5 to -1.2, there will be (1.3/0.02=)65 segments
+    for {set i 1} {$i <=66} {incr i} {
+	set eta [expr {-2.52 + $i * 0.02}]
+	add EtaPhiBins $eta $PhiBins
+    }
+    #same for 1.2 to 2.5
+    for  {set i 1} {$i <=66} {incr i} {
+	set eta [expr {1.18 + $i*0.02}]
+	add EtaPhiBins $eta $PhiBins
+    }
+
+
+    # default energy fractions {abs(PDG code)} {fraction of energy deposited in ECAL}
+
+    add EnergyFraction {0} {0.0}
+    # energy fractions for e, gamma and pi0
+    add EnergyFraction {11} {1.0}
+    add EnergyFraction {22} {1.0}
+    add EnergyFraction {111} {1.0}
+    # energy fractions for muon, neutrinos and neutralinos
+    add EnergyFraction {12} {0.0}
+    add EnergyFraction {13} {0.0}
+    add EnergyFraction {14} {0.0}
+    add EnergyFraction {16} {0.0}
+    add EnergyFraction {1000022} {0.0}
+    add EnergyFraction {1000023} {0.0}
+    add EnergyFraction {1000025} {0.0}
+    add EnergyFraction {1000035} {0.0}
+    add EnergyFraction {1000045} {0.0}
+    # energy fractions for K0short and Lambda
+    add EnergyFraction {310} {0.3}
+    add EnergyFraction {3122} {0.3}
+
+    # set ECalResolutionFormula {resolution formula as a function of eta and energy}
+    set ResolutionFormula {
+	(abs(eta) <= 0.78 )  * sqrt(energy^2*0.01^2 + energy*0.156^2)+
+	(abs(eta) > 0.78 && abs(eta) <=0.83 ) * sqrt( energy^0.01^2 + energy*0.175^2  ) +
+	(abs(eta) <= 2.5 && abs(eta) > 0.83) * sqrt( energy^2*0.01^2 + energy*0.151^2  )}
+}
+
+#############
+#   HCAL
+#############
+
+module SimpleCalorimeter HCal {
+    set ParticleInputArray ParticlePropagator/stableParticles
+    set TrackInputArray ECal/eflowTracks
+
+    set TowerOutputArray hcalTowers
+    set EFlowTrackOutputArray eflowTracks
+    set EFlowTowerOutputArray eflowNeutralHadrons
+
+    set IsEcal false
+
+    set EnergyMin 1.0
+    set EnergySignificanceMin 1.0
+
+    set SmearTowerCenter true
+
+    set pi [expr {acos(-1)}]
+
+    # lists of the edges of each tower in eta and phi
+    # each list starts with the lower edge of the first tower
+    # the list ends with the higher edged of the last tower
+
+
+    #HCAL barrel: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.8
+    #HCAL ring: dphi = 1 degree, deta= 0.02 towers up to |eta| <=0.9
+    #HCAL endcaps: dphi = 6 degree, deta = 0.1  up to |eta| <=2.5
+    #HCAL cell sizes always 30x30 mm^2
+
+    #barrel and ring:
+    #dphi = 1 degree up to |eta| <=0.9
+    set PhiBins {}
+    for {set i -180} {$i <=180} {incr i} {
+	add PhiBins [expr {$i * $pi/180.0}]
+    }
+    #deta= 0.02 towers up to |eta| <=0.9
+    for {set i -45} {$i <=45} {incr i} {
+	set eta [expr {$i * 0.02}]
+	add EtaPhiBins $eta $PhiBins
+    }
+
+    #endcaps:
+    # dphi = 6 degree
+    set PhiBins {}
+    for {set i -30} {$i <=30} {incr i} {
+	add PhiBins [expr {$i * $pi/30.0}]
+    }
+    # deta =0.1 for 0.9 < |eta| <=2.5
+    #for -2.5 to -0.9, 21 segments
+    for {set i 1} {$i <=17} {incr i} {
+	set eta [expr {-2.5 + $i * 0.1}]
+	add EtaPhiBins $eta $PhiBins
+    }
+    #same for 0.9 to 2.5
+    for {set i 1} {$i <=17} {incr i} {
+	set eta [expr {0.8 + $i * 0.1 }]
+	add EtaPhiBins $eta $PhiBins
+    }
+
+    # default energy fractions {abs(PDG code)} {Fecal Fhcal}
+    add EnergyFraction {0} {1.0}
+    # energy fractions for e, gamma and pi0
+    add EnergyFraction {11} {0.0}
+    add EnergyFraction {22} {0.0}
+    add EnergyFraction {111} {0.0}
+    # energy fractions for muon, neutrinos and neutralinos
+    add EnergyFraction {12} {0.0}
+    add EnergyFraction {13} {0.0}
+    add EnergyFraction {14} {0.0}
+    add EnergyFraction {16} {0.0}
+    add EnergyFraction {1000022} {0.0}
+    add EnergyFraction {1000023} {0.0}
+    add EnergyFraction {1000025} {0.0}
+    add EnergyFraction {1000035} {0.0}
+    add EnergyFraction {1000045} {0.0}
+    # energy fractions for K0short and Lambda
+    add EnergyFraction {310} {0.7}
+    add EnergyFraction {3122} {0.7}
+
+    # set HCalResolutionFormula {resolution formula as a function of eta and energy}
+    #CLICdet internal studies
+    set ResolutionFormula {
+	(abs(eta)<= 0.3) * sqrt(1.38^2  + energy*0.308^2  + energy^2*0.050^2) +
+	 (abs(eta)<= 0.78 && abs(eta) > 0.3) * sqrt(1.25^2  + energy*0.322^2  + energy^2*0.048^2) +
+	 (abs(eta)<=1.099 && abs(eta) > 0.78) * sqrt(  1.159^2 + energy*0.341^2 + energy^2*0.049^2 ) +
+	 (abs(eta)<=2.5 && abs(eta)> 1.099) * sqrt( 1.09^2 + energy*0.319^2 + energy^2*0.052^2  )
+	}
+
+}
+
+#################
+# Electron filter
+#################
+
+module PdgCodeFilter ElectronFilter {
+    set InputArray HCal/eflowTracks
+    set OutputArray electrons
+    set Invert true
+    add PdgCode {11}
+    add PdgCode {-11}
+}
+
+######################
+# ChargedHadronFilter
+######################
+
+module PdgCodeFilter ChargedHadronFilter {
+    set InputArray HCal/eflowTracks
+    set OutputArray chargedHadrons
+
+    add PdgCode {11}
+    add PdgCode {-11}
+    add PdgCode {13}
+    add PdgCode {-13}
+}
+
+
+
+###################################################
+# Tower Merger (in case not using e-flow algorithm)
+###################################################
+
+module Merger Calorimeter {
+    # add InputArray InputArray
+    add InputArray ECal/ecalTowers
+    add InputArray HCal/hcalTowers
+    add InputArray MuonMomentumSmearing/muons
+    set OutputArray towers
+}
+
+
+####################
+# Energy flow merger
+####################
+
+module Merger EFlowMerger {
+    # add InputArray InputArray
+    add InputArray HCal/eflowTracks
+    add InputArray ECal/eflowPhotons
+    add InputArray HCal/eflowNeutralHadrons
+    set OutputArray eflow
+}
+
+######################
+# EFlowFilter
+######################
+
+module PdgCodeFilter EFlowFilter {
+  set InputArray EFlowMerger/eflow
+  set OutputArray eflow
+
+  add PdgCode {11}
+  add PdgCode {-11}
+  add PdgCode {13}
+  add PdgCode {-13}
+}
+
+
+###################
+# Photon efficiency
+###################
+
+module Efficiency PhotonEfficiency {
+    set InputArray ECal/eflowPhotons
+    set OutputArray photons
+
+    # set EfficiencyFormula {efficiency formula as a function of eta and pt}
+
+    # efficiency formula for photons
+    set EfficiencyFormula {
+	(energy < 2.0 ) * (0.000) +
+	(energy >= 2.0) * (abs(eta) < 0.7)*(0.94) +
+	(energy >= 2.0) * (abs(eta) >=0.7 && abs(eta) <=2.5) * (0.9)	}
+
+}
+
+
+##################
+# Photon isolation
+##################
+
+module Isolation PhotonIsolation {
+    set CandidateInputArray PhotonEfficiency/photons
+    set IsolationInputArray EFlowMerger/eflow
+
+    set OutputArray photons
+
+    set DeltaRMax 0.1
+
+    set PTMin 0.5
+
+    set PTRatioMax 0.2
+}
+
+
+#####################
+# Electron efficiency
+#####################
+
+module Efficiency ElectronEfficiency {
+    set InputArray ElectronFilter/electrons
+    set OutputArray electrons
+
+    # set EfficiencyFormula {efficiency formula as a function of eta and pt}
+
+    set EfficiencyFormula {
+	(energy < 3.0 ) * ( 0.00 ) +
+  (abs(eta) > 2.50) * ( 0.00 ) +
+	( energy >=3 && energy < 8  ) * (abs(eta) > 1.95 && abs(eta) < 2.50)  * (0.58 ) +
+	( energy >=3 && energy < 8  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.7 ) +
+	( energy >=3 && energy < 8  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.6 ) +
+	( energy >=3 && energy < 8  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.7 ) +
+	( energy >=3 && energy < 8  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.8 ) +
+	( energy >=3 && energy < 8  ) * (abs(eta) <= 0.69)                    * (0.84 ) +
+	( energy >=8 && energy < 13  ) * (abs(eta) > 1.95 && abs(eta) < 2.50)  * (  0.6 ) +
+	( energy >=8 && energy < 13  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.76 ) +
+	( energy >=8 && energy < 13  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.67 ) +
+	( energy >=8 && energy < 13  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.78 ) +
+	( energy >=8 && energy < 13  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.86 ) +
+	( energy >=8 && energy < 13  ) * (abs(eta) <= 0.69)                    * ( 0.88 ) +
+	( energy >=13 && energy < 18  ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (  0.6 ) +
+	( energy >=13 && energy < 18  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (  0.8 ) +
+	( energy >=13 && energy < 18  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.68 ) +
+	( energy >=13 && energy < 18  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.84 ) +
+	( energy >=13 && energy < 18  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.88 ) +
+	( energy >=13 && energy < 18  ) * (abs(eta) <= 0.69)                    * (  0.9 ) +
+	( energy >=18 && energy < 23  ) * (abs(eta) > 1.95 && abs(eta) < 2.50)  * (0.64 ) +
+	( energy >=18 && energy < 23  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.82 ) +
+	( energy >=18 && energy < 23  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.7 ) +
+	( energy >=18 && energy < 23  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.84 ) +
+	( energy >=18 && energy < 23  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.9 ) +
+	( energy >=18 && energy < 23  ) * (abs(eta) <= 0.69)                    * (0.92 ) +
+	( energy >= 23 && energy < 28  ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.64 ) +
+	( energy >= 23 && energy < 28  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.86 ) +
+	( energy >= 23 && energy < 28  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.74 ) +
+	( energy >= 23 && energy < 28  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.87 ) +
+	( energy >= 23 && energy < 28  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.91 ) +
+	( energy >= 23 && energy < 28  ) * (abs(eta) <= 0.69)                    * (0.94 ) +
+	( energy >=28 && energy < 35  ) * (abs(eta) > 1.95 && abs(eta) < 2.50)   * (0.67 ) +
+	( energy >=28 && energy < 35  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.88 ) +
+	( energy >=28 && energy < 35  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.78 ) +
+	( energy >=28 && energy < 35  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.9 ) +
+	( energy >=28 && energy < 35  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
+	( energy >=28 && energy < 35  ) * (abs(eta) <= 0.69)                    * (0.94 ) +
+	( energy >=35 && energy < 45  ) * (abs(eta) > 1.95 && abs(eta) < 2.50)  * (0.68 ) +
+	( energy >=35 && energy < 45  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.9 ) +
+	( energy >=35 && energy < 45  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.86 ) +
+	( energy >=35 && energy < 45  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.92 ) +
+	( energy >=35 && energy < 45  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.94 ) +
+	( energy >=35 && energy < 45  ) * (abs(eta) <= 0.69)                    * (0.96 ) +
+	( energy >=45 && energy < 80  ) * (abs(eta) > 1.95 && abs(eta) < 2.50)  * (  0.7 ) +
+	( energy >=45 && energy < 80  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.92 ) +
+	( energy >=45 && energy < 80  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (  0.8 ) +
+	( energy >=45 && energy < 80  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.94 ) +
+	( energy >=45 && energy < 80  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.96 ) +
+	( energy >=45 && energy < 80  ) * (abs(eta) <= 0.69)                    * ( 0.97 ) +
+	( energy >=80 && energy < 200  ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * (0.68 ) +
+	( energy >=80 && energy < 200  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
+	( energy >=80 && energy < 200  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.84 ) +
+	( energy >=80 && energy < 200  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.94 ) +
+	( energy >=80 && energy < 200  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
+	( energy >=80 && energy < 200  ) * (abs(eta) <= 0.69)                    * (0.98 ) +
+	( energy >=200 && energy < 400  ) * (abs(eta) > 1.95 && abs(eta) < 2.50) * ( 0.68 ) +
+	( energy >=200 && energy < 400  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * ( 0.97 ) +
+	( energy >=200 && energy < 400  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * ( 0.86 ) +
+	( energy >=200 && energy < 400  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * ( 0.96 ) +
+	( energy >=200 && energy < 400  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * ( 0.98 ) +
+	( energy >=200 && energy < 400  ) * (abs(eta) <= 0.69)                    * ( 0.98 ) +
+	( energy >=400  ) * (abs(eta) > 1.95 && abs(eta) < 2.50)  * (0.68 ) +
+	( energy >=400  ) * (abs(eta) <= 1.95 && abs(eta) > 1.22) * (0.96 ) +
+	( energy >=400  ) * (abs(eta) <= 1.22 && abs(eta) > 1.1 ) * (0.82 ) +
+	( energy >=400  ) * (abs(eta) <= 1.1 && abs(eta) > 0.91 ) * (0.96 ) +
+	( energy >=400  ) * (abs(eta) <= 0.91 && abs(eta) > 0.69) * (0.98 ) +
+	( energy >=400  ) * (abs(eta) <= 0.69)                    * (0.98 )
+    }
+}
+
+####################
+# Electron isolation
+####################
+
+module Isolation ElectronIsolation {
+    set CandidateInputArray ElectronEfficiency/electrons
+    set IsolationInputArray EFlowMerger/eflow
+
+    set OutputArray electrons
+
+    set DeltaRMax 0.1
+
+    set PTMin 0.5
+
+    set PTRatioMax 0.2
+}
+
+#################
+# Muon efficiency
+#################
+
+module Efficiency MuonEfficiency {
+    set InputArray MuonMomentumSmearing/muons
+    set OutputArray muons
+
+    # set EfficiencyFormula {efficiency as a function of eta and pt}
+
+    set EfficiencyFormula {
+	(energy < 2.5 )     * (0.00) +
+	(energy>=2.5  )     * (0.999)
+    }
+}
+
+################
+# Muon isolation
+################
+
+module Isolation MuonIsolation {
+    set CandidateInputArray MuonEfficiency/muons
+    set IsolationInputArray EFlowMerger/eflow
+
+    set OutputArray muons
+
+    set DeltaRMax 0.1
+
+    set PTMin 0.5
+
+    set PTRatioMax 0.2
+}
+
+
+
+###################
+# Missing ET merger
+###################
+
+module Merger MissingET {
+    # add InputArray InputArray
+    add InputArray EFlowMerger/eflow
+    set MomentumOutputArray momentum
+}
+
+
+##################
+# Scalar HT merger
+##################
+
+module Merger ScalarHT {
+    # add InputArray InputArray
+    add InputArray EFlowMerger/eflow
+    set EnergyOutputArray energy
+}
+######################
+# EFlowFilter (UniqueObjectFinder)
+######################
+module UniqueObjectFinder UniqueObjectFinder {
+    add InputArray PhotonIsolation/photons photons
+    add InputArray ElectronIsolation/electrons electrons
+    add InputArray MuonIsolation/muons muons
+    add InputArray EFlowMerger/eflow eflow
+}
+
+#################
+# Neutrino Filter
+#################
+
+module PdgCodeFilter NeutrinoFilter {
+
+    set InputArray Delphes/stableParticles
+    set OutputArray filteredParticles
+
+    set PTMin 0.0
+
+    add PdgCode {12}
+    add PdgCode {14}
+    add PdgCode {16}
+    add PdgCode {-12}
+    add PdgCode {-14}
+    add PdgCode {-16}
+
+}
+
+
+#####################
+# MC truth jet finder
+#####################
+
+module FastJetFinder GenJetFinder {
+    set InputArray NeutrinoFilter/filteredParticles
+
+    set OutputArray jets
+
+    # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
+    set JetAlgorithm 9
+    set ParameterR 0.5
+
+    set JetPTMin 20.0
+}
+
+#########################
+# Gen Missing ET merger
+########################
+
+module Merger GenMissingET {
+    # add InputArray InputArray
+    add InputArray NeutrinoFilter/filteredParticles
+    set MomentumOutputArray momentum
+}
+
+
+
+############
+# Jet finder
+############
+
+module FastJetFinder FastJetFinderKt {
+    #  set InputArray Calorimeter/towers
+    set InputArray EFlowMerger/eflow
+
+    set OutputArray KTjets
+
+    # algorithm: 1 CDFJetClu, 2 MidPoint, 3 SIScone, 4 kt, 5 Cambridge/Aachen, 6 antikt, 7 anti-kt with winner-take-all axis (for N-subjettiness), 8 N-jettiness, 9 Valencia
+    set JetAlgorithm 4
+    set ParameterR 0.5
+
+    set JetPTMin 20.0
+}
+
+################
+# Jet finder VLC
+################
+
+
+source MuonCollider/MuonColliderDet_JetReco.tcl
+
+
+#########################################
+# Jet Momentum Smearing to mimick overlay
+#########################################
+
+
+source MuonCollider/MuonColliderDet_JetSmearing.tcl
+
+
+
+########################
+# Jet Flavor Association
+########################
+
+source  MuonCollider/MuonColliderDet_JetFlavorAssociation.tcl
+
+###########
+# b-tagging
+###########
+# based on CLICdp-Note-2014-002
+
+source  MuonCollider/MuonColliderDet_BTagging.tcl
+
+
+#############
+# tau-tagging
+#############
+
+source MuonCollider/MuonColliderDet_TauTagging.tcl
+
+##################
+# ROOT tree writer
+##################
+
+module TreeWriter TreeWriter {
+    # add Branch InputArray BranchName BranchClass
+    add Branch Delphes/allParticles Particle GenParticle
+
+    add Branch GenJetFinder/jets GenJet Jet
+
+    add Branch FastJetFinderKt/KTjets KTjet Jet
+
+    add Branch JetMomentumSmearing_VLCR05N2/JER_VLCjetsR05N2 VLCjetR05N2 Jet
+    add Branch JetMomentumSmearing_VLCR05N3/JER_VLCjetsR05N3 VLCjetR05N3 Jet
+    add Branch JetMomentumSmearing_VLCR05N4/JER_VLCjetsR05N4 VLCjetR05N4 Jet
+    add Branch JetMomentumSmearing_VLCR05N5/JER_VLCjetsR05N5 VLCjetR05N5 Jet
+    add Branch JetMomentumSmearing_VLCR05N6/JER_VLCjetsR05N6 VLCjetR05N6 Jet
+
+    add Branch JetMomentumSmearing_VLCR07N2/JER_VLCjetsR07N2 VLCjetR07N2 Jet
+    add Branch JetMomentumSmearing_VLCR07N3/JER_VLCjetsR07N3 VLCjetR07N3 Jet
+    add Branch JetMomentumSmearing_VLCR07N4/JER_VLCjetsR07N4 VLCjetR07N4 Jet
+    add Branch JetMomentumSmearing_VLCR07N5/JER_VLCjetsR07N5 VLCjetR07N5 Jet
+    add Branch JetMomentumSmearing_VLCR07N6/JER_VLCjetsR07N6 VLCjetR07N6 Jet
+
+    add Branch JetMomentumSmearing_VLCR10N2/JER_VLCjetsR10N2 VLCjetR10N2 Jet
+    add Branch JetMomentumSmearing_VLCR10N3/JER_VLCjetsR10N3 VLCjetR10N3 Jet
+    add Branch JetMomentumSmearing_VLCR10N4/JER_VLCjetsR10N4 VLCjetR10N4 Jet
+    add Branch JetMomentumSmearing_VLCR10N5/JER_VLCjetsR10N5 VLCjetR10N5 Jet
+    add Branch JetMomentumSmearing_VLCR10N6/JER_VLCjetsR10N6 VLCjetR10N6 Jet
+
+    add Branch JetMomentumSmearing_VLCR12N2/JER_VLCjetsR12N2 VLCjetR12N2 Jet
+    add Branch JetMomentumSmearing_VLCR12N3/JER_VLCjetsR12N3 VLCjetR12N3 Jet
+    add Branch JetMomentumSmearing_VLCR12N4/JER_VLCjetsR12N4 VLCjetR12N4 Jet
+    add Branch JetMomentumSmearing_VLCR12N5/JER_VLCjetsR12N5 VLCjetR12N5 Jet
+    add Branch JetMomentumSmearing_VLCR12N6/JER_VLCjetsR12N6 VLCjetR12N6 Jet
+
+    add Branch JetMomentumSmearing_VLCR15N2/JER_VLCjetsR15N2 VLCjetR15N2 Jet
+    add Branch JetMomentumSmearing_VLCR15N3/JER_VLCjetsR15N3 VLCjetR15N3 Jet
+    add Branch JetMomentumSmearing_VLCR15N4/JER_VLCjetsR15N4 VLCjetR15N4 Jet
+    add Branch JetMomentumSmearing_VLCR15N5/JER_VLCjetsR15N5 VLCjetR15N5 Jet
+    add Branch JetMomentumSmearing_VLCR15N6/JER_VLCjetsR15N6 VLCjetR15N6 Jet
+
+    add Branch JetMomentumSmearing_VLCR05_inclusive/JER_VLCjetsR05_inclusive VLCjetR05_inclusive Jet
+    add Branch JetMomentumSmearing_VLCR07_inclusive/JER_VLCjetsR07_inclusive VLCjetR07_inclusive Jet
+    add Branch JetMomentumSmearing_VLCR10_inclusive/JER_VLCjetsR10_inclusive VLCjetR10_inclusive Jet
+    add Branch JetMomentumSmearing_VLCR12_inclusive/JER_VLCjetsR12_inclusive VLCjetR12_inclusive Jet
+    add Branch JetMomentumSmearing_VLCR15_inclusive/JER_VLCjetsR15_inclusive VLCjetR15_inclusive Jet
+
+
+    ####
+
+    add Branch GenMissingET/momentum GenMissingET MissingET
+
+    add Branch TrackMerger/tracks Track Track
+    add Branch Calorimeter/towers Tower Tower
+
+    add Branch HCal/eflowTracks EFlowTrack Track
+    add Branch ECal/eflowPhotons EFlowPhoton Tower
+    add Branch HCal/eflowNeutralHadrons EFlowNeutralHadron Tower
+
+    add Branch UniqueObjectFinder/photons Photon Photon
+    add Branch UniqueObjectFinder/electrons Electron Electron
+    add Branch UniqueObjectFinder/muons Muon Muon
+
+    add Branch MissingET/momentum MissingET MissingET
+    add Branch ScalarHT/energy ScalarHT ScalarHT
+}
