Changeset 1345 in svn for trunk/examples

Timestamp:

Dec 21, 2013, 3:00:11 PM (11 years ago)

Author:

Michele Selvaggi

Message:

timing implemented

Location:

trunk/examples

Files:

• delphes_card_ATLAS_PileUp.tcl (modified) (1 diff)
• delphes_card_CMS_PileUp.tcl (modified) (1 diff)
• delphes_card_CMS_PileUp_Timing.tcl (added)

trunk/examples/delphes_card_ATLAS_PileUp.tcl

@@ -65 +65 @@
   # average expected pile up
   set MeanPileUp 50
-  # spread in the beam direction in m (assumes gaussian)
-  set ZVertexSpread 0.07
+  
+  # maximum spread in the beam direction in m 
+  set ZVertexSpread 0.10
+  
+  # maximum spread in time in s
+  set TVertexSpread 1.5E-09
+
+  # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
+  
+  set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
+
+  #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \
+  #                                (abs(t) >  1.0e-09) * (abs(z) <= 0.15) * (0.00) + \
+  #                               (abs(t) <= 1.0e-09) * (abs(z) > 0.15)  * (0.00) + \
+  #                               (abs(t) >  1.0e-09) * (abs(z) > 0.15)  * (0.00)}
+
+
 }
 

trunk/examples/delphes_card_CMS_PileUp.tcl

@@ -64 +64 @@
 
   # average expected pile up
-  set MeanPileUp 50
-  # spread in the beam direction in m (assumes gaussian)
-  set ZVertexSpread 0.07
+  set MeanPileUp 10
+  
+  # maximum spread in the beam direction in m 
+  set ZVertexSpread 0.10
+  
+  # maximum spread in time in s
+  set TVertexSpread 1.5E-09
+
+  # vertex smearing formula f(z,t) (z,t need to be respectively given in m,s)
+  
+  set VertexDistributionFormula {exp(-(t^2/(2*(0.05/2.99792458E8*exp(-(z^2/(2*(0.05)^2))))^2)))}
+
+  #set VertexDistributionFormula { (abs(t) <= 1.0e-09) * (abs(z) <= 0.15) * (1.00) + \
+  #                                (abs(t) >  1.0e-09) * (abs(z) <= 0.15) * (0.00) + \
+  #                               (abs(t) <= 1.0e-09) * (abs(z) > 0.15)  * (0.00) + \
+  #                               (abs(t) >  1.0e-09) * (abs(z) > 0.15)  * (0.00)}
+
+
 }