Merge pull request cms-sw#4310 from giamman/new-mixing-option-with-oot

new pileup mixing with OOT PU
cbernet · Jun 27, 2014 · 13caed1 · 13caed1
2 parents 45e27b5 + 60f8ee2
commit 13caed1
Show file tree

Hide file tree

Showing 2 changed files with 81 additions and 0 deletions.
diff --git a/Configuration/StandardSequences/python/Mixing.py b/Configuration/StandardSequences/python/Mixing.py
@@ -114,6 +114,7 @@ def addMixingScenario(label,dict):
 addMixingScenario("FS_2012_Summer_inTimeOnly",{'file': 'FastSimulation.PileUpProducer.PileUpSimulator_2012_Summer_inTimeOnly_cff'})
 addMixingScenario("FS_mix_2012_Startup_inTimeOnly",{'file': 'FastSimulation.PileUpProducer.mix_2012_Startup_inTimeOnly_cff'})
 addMixingScenario("FS_mix_2012_Summer_inTimeOnly",{'file': 'FastSimulation.PileUpProducer.mix_2012_Summer_inTimeOnly_cff'})
+addMixingScenario("FS_mix_2012_Summer_50ns_PoissonOOTPU",{'file': 'FastSimulation.PileUpProducer.mix_2012_Summer_50ns_PoissonOOTPU_cfi'}) # note: eventually we want absolute synchronization between FastSim and FullSim PU scenarios, and files like this one should be moved to SimGeneral/MixingModule ; to be done after full validation
 addMixingScenario("FS_CSA14_inTimeOnly",{'file': 'FastSimulation.PileUpProducer.PileUpSimulator_CSA14_inTimeOnly_cff'})
 addMixingScenario("FS_E13TeV_AVE_10_inTimeOnly",{'file': 'FastSimulation.PileUpProducer.PileUpSimulator_E13TeV_AVE_10_inTimeOnly_cff'})
 addMixingScenario("FS_E13TeV_AVE_20_inTimeOnly",{'file': 'FastSimulation.PileUpProducer.PileUpSimulator_E13TeV_AVE_20_inTimeOnly_cff'})

diff --git a/FastSimulation/PileUpProducer/python/mix_2012_Summer_50ns_PoissonOOTPU_cfi.py b/FastSimulation/PileUpProducer/python/mix_2012_Summer_50ns_PoissonOOTPU_cfi.py
@@ -0,0 +1,80 @@
+import FWCore.ParameterSet.Config as cms
+
+from FastSimulation.Configuration.CommonInputs_cff import MixingMode # CaloMode and MixingMode are defined here
+
+npu = cms.vint32(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59)
+prob = cms.vdouble(
+                         2.560E-06,
+                         5.239E-06,
+                         1.420E-05,
+                         5.005E-05,
+                         1.001E-04,
+                         2.705E-04,
+                         1.999E-03,
+                         6.097E-03,
+                         1.046E-02,
+                         1.383E-02,
+                         1.685E-02,
+                         2.055E-02,
+                         2.572E-02,
+                         3.262E-02,
+                         4.121E-02,
+                         4.977E-02,
+                         5.539E-02,
+                         5.725E-02,
+                         5.607E-02,
+                         5.312E-02,
+                         5.008E-02,
+                         4.763E-02,
+                         4.558E-02,
+                         4.363E-02,
+                         4.159E-02,
+                         3.933E-02,
+                         3.681E-02,
+                         3.406E-02,
+                         3.116E-02,
+                         2.818E-02,
+                         2.519E-02,
+                         2.226E-02,
+                         1.946E-02,
+                         1.682E-02,
+                         1.437E-02,
+                         1.215E-02,
+                         1.016E-02,
+                         8.400E-03,
+                         6.873E-03,
+                         5.564E-03,
+                         4.457E-03,
+                         3.533E-03,
+                         2.772E-03,
+                         2.154E-03,
+                         1.656E-03,
+                         1.261E-03,
+                         9.513E-04,
+                         7.107E-04,
+                         5.259E-04,
+                         3.856E-04,
+                         2.801E-04,
+                         2.017E-04,
+                         1.439E-04,
+                         1.017E-04,
+                         7.126E-05,
+                         4.948E-05,
+                         3.405E-05,
+                         2.322E-05,
+                         1.570E-05,
+                         5.005E-06)
+
+
+if (MixingMode=='DigiRecoMixing'):
+    # mix at SIM and RECO levels:
+    from FastSimulation.Configuration.MixingHitsAndTracks_cff import *
+    mix.input.nbPileupEvents.probFunctionVariable = npu
+    mix.input.nbPileupEvents.probValue = prob
+else:
+    print '#####################################################################################################################'
+    print '#####################################################################################################################'
+    print 'You are requesting a PU scenario that requires the DigiRecoMixing option, but you are actually using a different one!'
+    print '#####################################################################################################################'
+    print '#####################################################################################################################'
+