Skip to content

Commit

Permalink
Merge pull request #2 from gechen/branch-14-10-2014
Browse files Browse the repository at this point in the history
Merged thanks!
  • Loading branch information
xshi committed Oct 14, 2014
2 parents c5a93fb + 75991a1 commit 3fef448
Show file tree
Hide file tree
Showing 11 changed files with 397 additions and 476 deletions.
65 changes: 65 additions & 0 deletions python/BuToJPSiK_8TeV.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,65 @@
import FWCore.ParameterSet.Config as cms
from btokmumu_2012_cfi import process


process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10) )
#process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )

process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(
#'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_1_1_96m.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_2_1_dfI.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_3_1_ud5.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_4_1_neu.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_5_1_8Kc.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarJPsi_7TeV_PYTHIA6_100.root',
# 'file:/afs/cern.ch/user/w/wguo/work/CMSSW_5_3_11/src/wguo_xin/BToKMuMu/python/2012testMC_numEvent5000.root',
# 'file:/afs/cern.ch/work/n/nsahoo/BPH-ANALYSIS/afb/SE/aod/mc/2012_test_BToKMuMu_500.root',
# 'file:/afs/cern.ch/work/n/nsahoo/BPH-ANALYSIS/afb/SE/aod/mc/2012_test_BToKMuMu_1000.root',
#'root://xrootd.unl.edu//store/mc/Summer12_DR53X/Bu2MuMuK_TuneZ2star_8TeV_Pythia6/AODSIM/PU_RD2_START53_V19F-v1/00000/000771A0-9408-E411-94C0-002590DB92A8.root',
'root://xrootd.unl.edu//store/mc/Summer12_DR53X/BuJpsiK_TuneZ2star_8TeV_Pythia6/AODSIM/PU_RD2_START53_V19F-v1/00000/00086026-C705-E411-B655-90E6BAE8CC0C.root',
#'root://xrootd.unl.edu//store/mc/Summer12_DR53X/BuToPsiK_KFilter_TuneZ2star_SVS_8TeV-pythia6-evtgen/AODSIM/PU_RD2_START53_V19F-v1/00000/001E2794-0711-E411-B5F0-02163E008EBC.root',
# '/store/mc/Summer12_DR53X/BuToMuMuK_EtaPtFilter_8TeV-pythia6-evtgen/AODSIM/PU_S10_START53_V7C-v1/20000/9C4AFE3E-BA65-E211-821F-AC162DACC3E8.root',

)
)
#process.GlobalTag.globaltag = cms.string('GR_R_42_V25::All')r
#process.GlobalTag.globaltag = cms.string('START53_V23::All')
process.GlobalTag.globaltag = cms.string('START53_V7G::All')

# do trigger matching for muons
triggerProcessName = 'HLT'


process.cleanMuonTriggerMatchHLT0 = cms.EDProducer(
# match by DeltaR only (best match by DeltaR)
'PATTriggerMatcherDRLessByR',
src = cms.InputTag('cleanPatMuons'),
# default producer label as defined in
# PhysicsTools/PatAlgos/python/triggerLayer1/triggerProducer_cfi.py
matched = cms.InputTag('patTrigger'),
matchedCuts = cms.string('path("HLT_DoubleMu3p5_LowMass_Displaced*",0,0)'),
maxDeltaR = cms.double(0.1),
# only one match per trigger object
resolveAmbiguities = cms.bool(True),
# take best match found per reco object (by DeltaR here, see above)
resolveByMatchQuality = cms.bool(False))


from PhysicsTools.PatAlgos.tools.trigTools import *
switchOnTriggerMatchEmbedding(process, triggerMatchers = ['cleanMuonTriggerMatchHLT0'],
hltProcess = triggerProcessName, outputModule = '')

g_TriggerNames_LastFilterNames = [
('HLT_DoubleMu3p5_LowMass_Displaced', 'hltDisplacedmumuFilterDoubleMu3p5LowMass')
]

g_TriggerNames = [i[0] for i in g_TriggerNames_LastFilterNames]
g_LastFilterNames = [i[1] for i in g_TriggerNames_LastFilterNames]


process.ntuple.TriggerNames = cms.vstring(g_TriggerNames)
process.ntuple.LastFilterNames = cms.vstring(g_LastFilterNames)
process.ntuple.IsMonteCarlo = cms.untracked.bool(True)
process.ntuple.KeepGENOnly = cms.untracked.bool(False)
1 change: 1 addition & 0 deletions python/BuToPsi2SK_8TeV.py
Original file line number Diff line number Diff line change
Expand Up @@ -61,3 +61,4 @@
process.ntuple.TriggerNames = cms.vstring(g_TriggerNames)
process.ntuple.LastFilterNames = cms.vstring(g_LastFilterNames)
process.ntuple.IsMonteCarlo = cms.untracked.bool(True)
process.ntuple.KeepGENOnly = cms.untracked.bool(False)
144 changes: 144 additions & 0 deletions python/PYTHIA6_Bu2MuMuK_TuneZ2star_8TeV_GENOnly.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,144 @@
# Auto generated configuration file
# using:
# Revision: 1.381.2.28
# Source: /local/reps/CMSSW/CMSSW/Configuration/PyReleaseValidation/python/ConfigBuilder.py,v
# with command line options: ana/BToKMuMu/python/PYTHIA6_Bu2MuMuK_TuneZ2star_8TeV_GEN_cff.py --step GEN --beamspot Realistic8TeVCollision --conditions START53_LV6A1::All --pileup NoPileUp --datamix NODATAMIXER --eventcontent RAWSIM --datatier GEN --no_exec -n 10 --fileout BToKMuMu_MC_OnlyGEN_8TeV.root --python_filename PYTHIA6_Bu2MuMuK_TuneZ2star_8TeV_GENOnly.py
import FWCore.ParameterSet.Config as cms

process = cms.Process('GEN')

# import of standard configurations
process.load('Configuration.StandardSequences.Services_cff')
process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
process.load('FWCore.MessageService.MessageLogger_cfi')
process.load('Configuration.EventContent.EventContent_cff')
process.load('SimGeneral.MixingModule.mixNoPU_cfi')
process.load('Configuration.StandardSequences.GeometryRecoDB_cff')
process.load('Configuration.StandardSequences.MagneticField_38T_cff')
process.load('Configuration.StandardSequences.Generator_cff')
process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic8TeVCollision_cfi')
process.load('GeneratorInterface.Core.genFilterSummary_cff')
process.load('Configuration.StandardSequences.EndOfProcess_cff')
process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')

process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(10)
)

# Input source
process.source = cms.Source("EmptySource")

process.options = cms.untracked.PSet(

)

# Production Info
process.configurationMetadata = cms.untracked.PSet(
version = cms.untracked.string('$Revision: 1.1 $'),
annotation = cms.untracked.string('B+ -> mu+ mu- K+ at 8TeV'),
name = cms.untracked.string('$Source: /local/reps/CMSSW/CMSSW/Configuration/GenProduction/python/EightTeV/PYTHIA6_Bu2MuMuK_TuneZ2star_8TeV_cff.py,v $')
)

# Output definition

process.RAWSIMoutput = cms.OutputModule("PoolOutputModule",
splitLevel = cms.untracked.int32(0),
eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
outputCommands = process.RAWSIMEventContent.outputCommands,
fileName = cms.untracked.string('BToKMuMu_MC_OnlyGEN_8TeV.root'),
dataset = cms.untracked.PSet(
filterName = cms.untracked.string(''),
dataTier = cms.untracked.string('GEN')
),
SelectEvents = cms.untracked.PSet(
SelectEvents = cms.vstring('generation_step')
)
)

# Additional output definition

# Other statements
process.genstepfilter.triggerConditions=cms.vstring("generation_step")
from Configuration.AlCa.GlobalTag import GlobalTag
process.GlobalTag = GlobalTag(process.GlobalTag, 'START53_LV6A1::All', '')

process.decayfilter = cms.EDFilter("PythiaDauVFilter",
MinPt = cms.untracked.vdouble(-1.0, 0, 0),
ParticleID = cms.untracked.int32(521),
MaxEta = cms.untracked.vdouble(9999.0, 9999, 9999),
MinEta = cms.untracked.vdouble(-9999.0, -9999, -9999),
DaughterIDs = cms.untracked.vint32(321, -13, 13),
NumberDaughters = cms.untracked.int32(3),
verbose = cms.untracked.int32(0)
)


process.bfilter = cms.EDFilter("PythiaFilter",
MaxEta = cms.untracked.double(9999.0),
MinEta = cms.untracked.double(-9999.0),
ParticleID = cms.untracked.int32(521)
)


process.generator = cms.EDFilter("Pythia6GeneratorFilter",
ExternalDecays = cms.PSet(
EvtGen = cms.untracked.PSet(
use_default_decay = cms.untracked.bool(False),
decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY_NOLONGLIFE.DEC'),
particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
user_decay_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/Bu_mumuK.dec'),
list_forced_decays = cms.vstring('MyB+',
'MyB-'),
operates_on_particles = cms.vint32(0)
),
parameterSets = cms.vstring('EvtGen')
),
pythiaPylistVerbosity = cms.untracked.int32(0),
filterEfficiency = cms.untracked.double(0.003),
pythiaHepMCVerbosity = cms.untracked.bool(False),
comEnergy = cms.double(8000.0),
crossSection = cms.untracked.double(48440000000.0),
maxEventsToPrint = cms.untracked.int32(0),
PythiaParameters = cms.PSet(
pythiaUESettings = cms.vstring('MSTU(21)=1 ! Check on possible errors during program execution',
'MSTJ(22)=2 ! Decay those unstable particles',
'PARJ(71)=10 . ! for which ctau 10 mm',
'MSTP(33)=0 ! no K factors in hard cross sections',
'MSTP(2)=1 ! which order running alphaS',
'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)',
'MSTP(52)=2 ! work with LHAPDF',
'PARP(82)=1.921 ! pt cutoff for multiparton interactions',
'PARP(89)=1800. ! sqrts for which PARP82 is set',
'PARP(90)=0.227 ! Multiple interactions: rescaling power',
'MSTP(95)=6 ! CR (color reconnection parameters)',
'PARP(77)=1.016 ! CR',
'PARP(78)=0.538 ! CR',
'PARP(80)=0.1 ! Prob. colored parton from BBR',
'PARP(83)=0.356 ! Multiple interactions: matter distribution parameter',
'PARP(84)=0.651 ! Multiple interactions: matter distribution parameter',
'PARP(62)=1.025 ! ISR cutoff',
'MSTP(91)=1 ! Gaussian primordial kT',
'PARP(93)=10.0 ! primordial kT-max',
'MSTP(81)=21 ! multiple parton interactions 1 is Pythia default',
'MSTP(82)=4 ! Defines the multi-parton model'),
bbbarSettings = cms.vstring('MSEL = 1'),
parameterSets = cms.vstring('pythiaUESettings',
'bbbarSettings')
)
)


process.ProductionFilterSequence = cms.Sequence(process.generator+process.bfilter+process.decayfilter)

# Path and EndPath definitions
process.generation_step = cms.Path(process.pgen)
process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
process.endjob_step = cms.EndPath(process.endOfProcess)
process.RAWSIMoutput_step = cms.EndPath(process.RAWSIMoutput)

# Schedule definition
process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.endjob_step,process.RAWSIMoutput_step)
# filter all path with the production filter sequence
for path in process.paths:
getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq

9 changes: 5 additions & 4 deletions python/btokmumu_2012_cfi.py
Original file line number Diff line number Diff line change
Expand Up @@ -151,7 +151,7 @@
MuonMass = cms.untracked.double(0.10565837),
MuonMassErr = cms.untracked.double(3.5*1e-9),
KaonMass = cms.untracked.double(0.493677),
KaonMassErr = cms.untracked.double(1.6*1e-5),
KaonMassErr = cms.untracked.double(1.6*1e-5),
# PionMass = cms.untracked.double(0.13957018),
# PionMassErr = cms.untracked.double(0.13957018*1e-6),
#KshortMass = cms.untracked.double(0.497614),
Expand All @@ -173,7 +173,8 @@

# gen particle
IsMonteCarlo = cms.untracked.bool(False),
TruthMatchMuonMaxR = cms.untracked.double(0.004), # [eta-phi]
KeepGENOnly = cms.untracked.bool(False),
TruthMatchMuonMaxR = cms.untracked.double(0.004), # [eta-phi]
TruthMatchKaonMaxR = cms.untracked.double(0.3), # [eta-phi]
# TruthMatchKsMaxVtx = cms.untracked.double(10.0),

Expand All @@ -192,8 +193,8 @@
MuMuMinCosAlphaBs = cms.untracked.double(0.9),

# pre-selection cuts
TrkMinPt = cms.untracked.double(0.4), # [GeV/c]
TrkMinDcaSigBs = cms.untracked.double(1.2), # hadron DCA/sigma w/respect to BS [1.2]
TrkMinPt = cms.untracked.double(0.2), # [GeV/c]
TrkMinDcaSigBs = cms.untracked.double(0.1), # hadron DCA/sigma w/respect to BS [1.2]
TrkMaxR = cms.untracked.double(110.0), # [cm]
TrkMaxZ = cms.untracked.double(280.0), # [cm]

Expand Down
10 changes: 5 additions & 5 deletions python/btokmumu_Data_8TeV.py
Original file line number Diff line number Diff line change
@@ -1,13 +1,13 @@
import FWCore.ParameterSet.Config as cms
from btokmumu_2012_cfi import process


process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10) )
#process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
#process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10) )
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )

process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(
'root://xrootd.unl.edu//store/data/Run2012A/DoubleMuParked/AOD/22Jan2013-v1/20000/002557DC-16D4-E211-9415-20CF3027A5AB.root',
#'root://xrootd.unl.edu//store/data/Run2012A/DoubleMuParked/AOD/22Jan2013-v1/20000/002557DC-16D4-E211-9415-20CF3027A5AB.root',
'root://xrootd.unl.edu//store/data/Run2012B/MuOniaParked/AOD/22Jan2013-v1/20000/00054BD0-5668-E211-8091-00215E21DC7E.root',
#'file:/afs/cern.ch/work/n/nsahoo/BPH-ANALYSIS/afb/SE/aod/data/2012_test_data_Run2012B_Jan2013_1000.root',
#'file:/afs/cern.ch/work/n/nsahoo/BPH-ANALYSIS/afb/SE/aod/data/2012_test_data_Run2012B_15000.root'
#'/store/data/Run2012A/DoubleMuParked/AOD/22Jan2013-v1/20000/00A96C5F-0ED4-E211-A18C-485B39800C15.root'
Expand Down Expand Up @@ -52,4 +52,4 @@

process.ntuple.TriggerNames = cms.vstring(g_TriggerNames)
process.ntuple.LastFilterNames = cms.vstring(g_LastFilterNames)

process.ntuple.KeepGENOnly = cms.untracked.bool(False)
1 change: 1 addition & 0 deletions python/btokmumu_MC_8TeV.py
Original file line number Diff line number Diff line change
Expand Up @@ -59,3 +59,4 @@
process.ntuple.TriggerNames = cms.vstring(g_TriggerNames)
process.ntuple.LastFilterNames = cms.vstring(g_LastFilterNames)
process.ntuple.IsMonteCarlo = cms.untracked.bool(True)
process.ntuple.KeepGENOnly = cms.untracked.bool(False)
29 changes: 29 additions & 0 deletions python/btokmumu_MC_GENOnly.py
Original file line number Diff line number Diff line change
@@ -0,0 +1,29 @@
import FWCore.ParameterSet.Config as cms
from btokmumu_2012_cfi import process

#process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(100) )
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )

process.source = cms.Source("PoolSource",
fileNames = cms.untracked.vstring(
#'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_1_1_96m.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_2_1_dfI.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_3_1_ud5.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_4_1_neu.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarMuMu_7TeV_PYTHIA6_5_1_8Kc.root',
# 'file:/afs/cern.ch/user/x/xshi/work/cms/afb/dat/aod/mc/BuToKstarJPsi_7TeV_PYTHIA6_100.root',
# 'file:/afs/cern.ch/user/w/wguo/work/CMSSW_5_3_11/src/wguo_xin/BToKMuMu/python/2012testMC_numEvent5000.root',
# 'file:/afs/cern.ch/work/n/nsahoo/BPH-ANALYSIS/afb/SE/aod/mc/2012_test_BToKMuMu_500.root',
'file:/afs/cern.ch/user/g/gechen/gechen/work/BToKMuMu/Ana/RootFiles/MC_GENOnly/BToKMuMu_MC_OnlyGEN_8TeV_9.root',
# '/store/mc/Summer12_DR53X/BuToMuMuK_EtaPtFilter_8TeV-pythia6-evtgen/AODSIM/PU_S10_START53_V7C-v1/20000/9C4AFE3E-BA65-E211-821F-AC162DACC3E8.root',

)
)
#process.GlobalTag.globaltag = cms.string('GR_R_42_V25::All')r
#process.GlobalTag.globaltag = cms.string('START53_V23::All')
process.GlobalTag.globaltag = cms.string('START53_V7G::All')

process.ntuple.IsMonteCarlo = cms.untracked.bool(True)
process.ntuple.KeepGENOnly = cms.untracked.bool(True)
process.p = cms.Path(process.ntuple)
34 changes: 34 additions & 0 deletions python/cb_MC_GENOnly_8TeV.cfg
Original file line number Diff line number Diff line change
@@ -0,0 +1,34 @@
[CRAB]
jobtype = cmssw
scheduler = remoteGlidein
use_server = 0


[CMSSW]
get_edm_output = 1
generator = pythia
total_number_of_events = 500000000
number_of_jobs = 2000
pset = ./PYTHIA6_Bu2MuMuK_TuneZ2star_8TeV_GENOnly.py
datasetpath = None
output_file = BToKMuMu_MC_OnlyGEN_8TeV.root
#dbs_url = phys03

[USER]
ui_working_dir = BuToKMuMu_MC_GENOnly_8TeV_v3
user_remote_dir = BuToKMuMu_MC_GENOnly_8TeV_v3
#local_stage_out = 1
return_data = 0
storage_element = T2_CN_Beijing
#storage_element = T2_UK_SGrid_Bristol
copy_data = 1
publish_data = 1
publish_data_name = PYTHIA6_Bu2MuMuK_TuneZ2star_GENOnly_8TeV_v3
dbs_url_for_publication = phys03

[GRID]
rb = CERN
proxy_server = myproxy.cern.ch



35 changes: 35 additions & 0 deletions python/cb_MC_GENOnly_8TeV_Ntuple.cfg
Original file line number Diff line number Diff line change
@@ -0,0 +1,35 @@
[CRAB]
jobtype = cmssw
scheduler = remoteGlidein
use_server = 0


[CMSSW]
get_edm_output = 1
#generator = pythia
total_number_of_events = -1
number_of_jobs = 2000
pset = ./btokmumu_MC_GENOnly.py
datasetpath = /PYTHIA6_Bu2MuMuK_TuneZ2star_GENOnly_8TeV_v3/gechen-PYTHIA6_Bu2MuMuK_TuneZ2star_GENOnly_8TeV_v3-afd5e29ce33c42baca35170346e3548c/USER
#datasetpath = /PYTHIA6_Bu2MuMuK_TuneZ2star_GENOnly_8TeV_v2/gechen-PYTHIA6_Bu2MuMuK_TuneZ2star_GENOnly_8TeV_v2-afd5e29ce33c42baca35170346e3548c/USER
output_file = BToKMuMu_GENOnly_8TeV_Ntuple.root
dbs_url = phys03

[USER]
ui_working_dir = BuToKMuMu_MC_GENOnly_8TeV_Ntuples_v3
user_remote_dir = BuToKMuMu_MC_GENOnly_8TeV_Ntuples_v3
#local_stage_out = 1
return_data = 0
storage_element = T2_CN_Beijing
#storage_element = T2_UK_SGrid_Bristol
copy_data = 1
#publish_data = 1
#publish_data_name = PYTHIA6_Bu2MuMuK_TuneZ2star_GENOnly_8TeV_v2
#dbs_url_for_publication = phys03

[GRID]
rb = CERN
proxy_server = myproxy.cern.ch



Loading

0 comments on commit 3fef448

Please sign in to comment.