Exemplo n.º 1
0
def options() :
    options = VP('standard')
    options.output = "topTuple.root"
    options.maxEvents = 100
    
    options.register('isData', default = True, mytype = VP.varType.bool)
    options.register('isFastSim', default = False, mytype = VP.varType.bool)
    options.register('quiet', default = True, mytype = VP.varType.bool)
    options.register('requireLepton', default = True, mytype = VP.varType.bool)
    options.register('globalTag', mytype = VP.varType.string )
    options.register('postfix','TR', mytype = VP.varType.string )
    options.register('btags', mytype = VP.varType.string, mult = VP.multiplicity.list )
    options.register('doElectronEA', default = True, mytype = VP.varType.bool)
    options.parseArguments()
    options._tagOrder =[]

    defaultGT = ('FT53_V21A_AN6' if options.isData else 'START53_V15')

    sync53 = '/store/mc/Summer12_DR53X/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/AODSIM/PU_S10_START53_V7A-v1/0000/FE4C2F81-D0E1-E111-9080-0030487E0A2D.root'
    SingleElRun2012A = "root://cmsxrootd.fnal.gov///store/data/Run2012A/SingleElectron/AOD/13Jul2012-v1/0000/001A2EB8-47D4-E111-B527-003048679070.root" 
    defaultFiles = ( [sync53] if not options.isData else [SingleElRun2012A] )
    # pickRelValInputFiles( cmsswVersion = 'CMSSW_5_3_6',
    #                       dataTier     = 'RECO',
    #                       relVal       = 'SingleMu',
    #                       globalTag    = 'GR_R_53_V15_RelVal_mu2012A',
    #                       numberOfFiles = 10
    #                      )
    
    options.files = options.files if options.files else defaultFiles
    if not options.globalTag : options.globalTag = defaultGT

    if not options.quiet : print options
    options.btags = ['combinedSecondaryVertex','jetProbability']
    return options
Exemplo n.º 2
0
def options() :
    options = VP('standard')
    options.output = "topTuple.root"
    options.maxEvents = -1
    

    options.register('procName', default = "PAT2", mytype = VP.varType.string)
    options.register('isData', default = False, mytype = VP.varType.bool)
    options.register('skim', default = True, mytype = VP.varType.bool)
    options.register('btag', default = False, mytype = VP.varType.bool)
#    options.register('outputModule', default = True, mytype = VP.varType.bool)
    options.register('quiet', default = False, mytype = VP.varType.bool)
    options.register('requireLepton', default = True, mytype = VP.varType.bool)
    options.register('globalTag', mytype = VP.varType.string )
    options.register('postfix','TR', mytype = VP.varType.string )
    options.register('btags', mytype = VP.varType.string, mult = VP.multiplicity.list )
    options.register('doElectronEA', default = True, mytype = VP.varType.bool)
    options.register('noJetSmearing', default = False, mytype = VP.varType.bool)

    try:
        for key in options._register.keys():
            setattr(options, key, cms_var.get(key, getattr(options, key)))
    except NameError:
        print "cms_var is not in __builtin__"

    print sys.argv
    if not "-create" in sys.argv:
        options.parseArguments()
    else:
        print "Found '-create' in sys.argv! Omitting options.parseArguments()"
    options._tagOrder =[]



    defaultGT = ('GR_R_53_V18' if options.isData else 'START53_V15')

    # FT_53_V21_AN4
    # START53_V23

    sync53 = 'file:/user/kuessel/files/Synchfiles/syncExercise53.root'
    diffStep1="file:/user/kuessel/CMSSW/Synch/CMSSW_5_3_8_patch3/src/pickevents.root"
    subset="file:/user/kuessel/CMSSW/Synch/CMSSW_5_3_8_patch3/src/testSkim.root"
    skim="/store/user/kuessel/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/YK_MC_MARCH13_TTbar/accebdcbeca2810af478d8af2493d41f/SynchSelMuJets_103_1_lgS.root"
    defaultFiles = [skim]
#    defaultFiles = [subset]
#    defaultFiles = [diffStep1] 
    
    options.files = defaultFiles
    if not options.globalTag : 
        options.globalTag = defaultGT

    options.btags = ['combinedSecondaryVertex','jetProbability']

    if not options.quiet : print options
    return options
Exemplo n.º 3
0
def options() :
    from FWCore.ParameterSet.VarParsing import VarParsing as VP
    options = VP('standard')
    options.output = "DJ_Tree.root"
    options.maxEvents = 100

    options.register('GlobalTag', mytype = VP.varType.string)
    options.register('quiet', default = True )
    options.register('isData', default = True )
    options.register('signal', default = False )
    options.register('jetCorrections', default = ['L1FastJet','L2Relative','L3Absolute','L2L3Residual'], #L2L3Residual removed from options for simulation (below)
                     info = "jet correction levels to apply", mult = VP.multiplicity.list, mytype = VP.varType.string)
    
    __hack_ListVarparsingBug__( options, 'jetCorrections')
    options.parseArguments()
    options._tagOrder =[] # weird, but something to do with options.output

    defaultGT,defaultFile = {
        "44X" : [('START44_V9B::All','/store/mc/Fall11/QCD_Pt-120to170_TuneZ2_7TeV_pythia6/AODSIM/PU_S6_START42_V14B-v1/0000/0A8E5F68-73F4-E011-AB1C-00E08178C0D5.root'),
                 ('START44_V9B::All','/store/mc/Fall11/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau-350_7TeV-pythia6/GEN-SIM-RECODEBUG/DEBUG-PU_S6_START44_V9B-v4/0000/E83B1DB8-158B-E111-9DB4-0017A477003C.root'),
                 ('FT_R_44_V11::All','/store/data/Run2011B/HT/RECO/19Nov2011-v1/10001/B2E314BC-4746-E111-9BC1-00261894390A.root')],

        "53X" : [
                 #('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau35To3500_8TeV-pythia6/AODSIM/DEBUG_PU_S10_START53_V7A-v1/0000/9C61CA50-0BF3-E111-91E1-003048D2BB22.root'),
                 #('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-400_MFF-50_CTau8To800_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/80A0705C-7518-E211-A17D-0030487E4D11.root'),
                 #('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-200_MFF-50_CTau20To2000_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/7A2CB2BB-6C18-E211-AB1B-003048C69416.root'),
		 ('START53_V7F::All','/store/mc/Summer12_DR53X/LongLivedChi0ToMuQQ_MSquark-1500_MChi-494_TuneZ2Star_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/B69AB21F-D50F-E211-B2A7-00215E21DCA2.root'),
		 ('START53_V7F::All','/store/mc/Summer12_DR53X/QCD_Pt-170to300_TuneZ2star_8TeV_pythia6/AODSIM/PU_S10_START53_V7A-v2/00000/0CECB336-B70E-E211-B22D-0018F3D096CE.root'),
		 ('FT_53_V6_AN2::All','/store/data/Run2012B/JetHT/AOD/13Jul2012-v1/00000/04989369-BFD2-E111-B190-002618943933.root'),
		 ('FT_53_V10_AN2::All','/store/data/Run2012C/JetHT/AOD/24Aug2012-v2/00000/22AFA398-D600-E211-AD2E-00261894385A.root'),
		 ('GR_P_V41_AN2::All','/store/data/Run2012C/JetHT/AOD/PromptReco-v2/000/198/955/144C3A0C-D8CF-E111-8A6C-001D09F27067.root'),
		 ('GR_P_V42_AN2::All','/store/data/Run2012D/JetHT/AOD/PromptReco-v1/000/204/506/F8F3A290-8C10-E211-B2D3-001D09F2525D.root'),],

        "53XCtrl" : [('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau35To3500_8TeV-pythia6/AODSIM/DEBUG_PU_S10_START53_V7A-v1/0000/9C61CA50-0BF3-E111-91E1-003048D2BB22.root'),
                 ('START53_V7F::All','/store/mc/Summer12_DR53X/QCD_Pt-170to300_TuneZ2star_8TeV_pythia6/AODSIM/PU_S10_START53_V7A-v2/00000/0CECB336-B70E-E211-B22D-0018F3D096CE.root'),
                 ('FT_53_V6_AN2::All','/store/data/Run2012B/SingleMu/AOD/13Jul2012-v1/0000/24DF4BD2-0DD4-E111-BB00-E0CB4E4408DF.root'),
                 ('FT_53_V10_AN2::All','/store/data/Run2012C/SingleMu/AOD/24Aug2012-v1/00000/10FD4F26-FAEE-E111-9AB5-00259073E4FC.root'),
                 ('GR_P_V41_AN2::All','/store/data/Run2012C/SingleMu/AOD/PromptReco-v2/000/198/955/247727CE-98CF-E111-B2EF-5404A63886A0.root'),
                 ('GR_P_V42_AN2::All','/store/data/Run2012D/SingleMu/AOD/PromptReco-v1/000/204/506/8078B6DD-9210-E211-80BE-BCAEC5329700.root'),]
        }["53X"][int(not options.signal)+ options.isData]
    options.files = options.files if options.files else defaultFile
    options.GlobalTag = options.GlobalTag if options.GlobalTag else defaultGT

    if not options.isData : #remove L2L3Residual correction from simulation options
        jecs = [jc for jc in options.jetCorrections if jc!='L2L3Residual']
        options.clearList('jetCorrections')
        options.jetCorrections = jecs

    return options
Exemplo n.º 4
0
def options() :
    from FWCore.ParameterSet.VarParsing import VarParsing as VP
    options = VP('standard')
    options.output = "SusyCAF_Tree.root"
    options.secondaryOutput = "" #switch PAT-tuple output off by default
    options.maxEvents = 100

    options.register('GlobalTag', mytype = VP.varType.string)
    options.register('quiet', default = True )
    options.register('isData', default = True )
    options.register('dropMore', default = False , info = "drop more for a smaller event size")
    options.register('leptonPtMin', default = 0 )
    options.register('patify', default = True,  info = "run SUSYPAT on the fly")
    options.register('AllTracks', default = False, info = "include all tracks")
    options.register('SourceName', info = "'S:stream' or 'DS:dataset' to store those HLT paths", mytype = VP.varType.string,)
    options.register('jetCollections', default = ['ak5calo','ak5pf'], # + ['ak7calo','ak7pf','ak5pf2pat','ak7pf2pat']
                     info = "jet types to store", mult = VP.multiplicity.list, mytype = VP.varType.string)
    options.register('jetCorrections', default = ['L1FastJet','L2Relative','L3Absolute','L2L3Residual'], #L2L3Residual removed from options for simulation (below)
                     info = "jet correction levels to apply", mult = VP.multiplicity.list, mytype = VP.varType.string)
    options.register('scan', default = "", info = "code for CMSSM or SMS scan", mytype = VP.varType.string)
    options.register('triggers', default = -1, info = "store trigger information in the ntuple") # defaults to isData (below)
    options.register('taus', default = True, info = "store tau information in the ntuple")
    options.register('dqm', default = False, info = "store dqm information in the ntuple")
    options.register('doTauReco', default = False, info = "redo tau recostruction (needed in order to read 51X files with 52X recipe)")
    options.register('doTypeIMetReco', default = False, info = "do type I MET recostruction")
    options.register('doTypeIMetPat', default = True, info = "do type I MET Pat (RA4 Recipe + also caloMET)")
    options.register('doPfMetPhiCorrections', default = False, info = "do phi corrections in pfMET")
    options.register('beamHaloVars', default = True, info = "store beam halo information in the ntuple")
    options.register('doPhotonPFIso', default = True, info = "include simple cut based Photon IDs 2012 in the ntuple")
    options.register('doPoolOutput', default = False, info = "Output edm file with full event content")
    options.register('pdfSets', default = "", info = "pdf Sets to use for pdf Uncertainty weights", mytype = VP.varType.string)
    options.register('defPdf', default = "cteq6ll.LHpdf", info = "pdf set used when producing sample", mytype = VP.varType.string)
    
    __hack_ListVarparsingBug__( options, 'jetCollections')
    __hack_ListVarparsingBug__( options, 'jetCorrections')
    options.parseArguments()
    options._tagOrder =[] # weird, but something to do with options.output

    defaultGT,defaultFile = {
        "53X" : [('START53_V10::All','/store/mc/Summer12_DR53X/ZZ_TuneZ2star_8TeV_pythia6_tauola/AODSIM/PU_S10_START53_V7A-v1/0000/FEEB0275-24FA-E111-B9A6-00266CF9B970.root'),
		('GR_P_V41_AN1::All' ,'/store/data/Run2012C/SingleMu/AOD/PromptReco-v2/000/201/624/F40EFD20-9CF0-E111-8F16-5404A63886C4.root')],
	"52X" : [('START52_V4::All','/store/relval/CMSSW_5_2_5_cand1/RelValTTbar/GEN-SIM-RECO/START52_V9-v1/0261/249F81B9-8C91-E111-AE72-003048679236.root'),               
		('GR_R_52_V4::All' ,'/store/relval/CMSSW_5_2_5_cand1/SingleMu/RECO/GR_R_52_V7_RelVal_mu2011B-v1/0262/7425B24B-9E91-E111-BA46-002618943800.root')],
        "51X" : [('START50_V15::All','/store/relval/CMSSW_5_1_2/RelValTTbar/GEN-SIM-RECO/PU_START50_V15A-v1/0003/0AB02E79-7C61-E111-8594-002481E0D73C.root'),
                 ('GR_R_50_V12::All' ,'/store/relval/CMSSW_5_1_2/Jet/RECO/GR_R_50_V12_RelVal_jet2010B-v1/0237/326E6DA4-B460-E111-AAFA-0026189438D6.root')],
        "44X" : [('START44_V12::All','/../user/b/bbetchar//CMSSW_4_4_2_patch10/RelValProdTTbar/AODSIM/START44_V7_special_120119-v1/0088/0ECC7FA5-7A43-E111-9B36-002618943861.root'),
                 ('GR_R_44_V13::All','/../user/b/bbetchar/data/Run2011A/SingleMu/AOD/08Nov2011-v1/0001/80D61727-200D-E111-B7F4-1CC1DE0571C8.root')],
        "42X" : [('START42_V17::All','/../user/b/bbetchar/CMSSW_4_2_8/RelValProdTTbar/GEN-SIM-RECO/MC_42_V12-v1/0026/9C0E8835-9ABB-E011-95B0-0026189438BA.root'),
                 ('GR_R_42_V24::All','/store/data/Run2011B/MultiJet/AOD/PromptReco-v1/000/175/832/485ABBC8-10DC-E011-980B-BCAEC518FF8E.root')]
        }["53X"][options.isData]
    options.files = options.files if options.files else defaultFile
    options.GlobalTag = options.GlobalTag if options.GlobalTag else defaultGT

    if options.triggers < 0 : options.triggers = int(options.isData)

    if not options.isData : #remove L2L3Residual correction from simulation options
        jecs = [jc for jc in options.jetCorrections if jc!='L2L3Residual']
        options.clearList('jetCorrections')
        options.jetCorrections = jecs

    return options
Exemplo n.º 5
0
def options():
    from FWCore.ParameterSet.VarParsing import VarParsing as VP
    options = VP('standard')
    options.output = "DJ_Tree.root"
    options.maxEvents = 100

    options.register('GlobalTag', mytype=VP.varType.string)
    options.register('quiet', default=True)
    options.register('isData', default=True)
    options.register('signal', default=False)
    options.register(
        'jetCorrections',
        default=['L1FastJet', 'L2Relative', 'L3Absolute', 'L2L3Residual'
                 ],  #L2L3Residual removed from options for simulation (below)
        info="jet correction levels to apply",
        mult=VP.multiplicity.list,
        mytype=VP.varType.string)

    __hack_ListVarparsingBug__(options, 'jetCorrections')
    options.parseArguments()
    options._tagOrder = []  # weird, but something to do with options.output

    defaultGT, defaultFile = {
        "44X":
        [('START44_V9B::All',
          '/store/mc/Fall11/QCD_Pt-120to170_TuneZ2_7TeV_pythia6/AODSIM/PU_S6_START42_V14B-v1/0000/0A8E5F68-73F4-E011-AB1C-00E08178C0D5.root'
          ),
         ('START44_V9B::All',
          '/store/mc/Fall11/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau-350_7TeV-pythia6/GEN-SIM-RECODEBUG/DEBUG-PU_S6_START44_V9B-v4/0000/E83B1DB8-158B-E111-9DB4-0017A477003C.root'
          ),
         ('FT_R_44_V11::All',
          '/store/data/Run2011B/HT/RECO/19Nov2011-v1/10001/B2E314BC-4746-E111-9BC1-00261894390A.root'
          )],
        "53X": [
            #('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau35To3500_8TeV-pythia6/AODSIM/DEBUG_PU_S10_START53_V7A-v1/0000/9C61CA50-0BF3-E111-91E1-003048D2BB22.root'),
            #('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-400_MFF-50_CTau8To800_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/80A0705C-7518-E211-A17D-0030487E4D11.root'),
            #('START53_V7F::All','/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-200_MFF-50_CTau20To2000_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/7A2CB2BB-6C18-E211-AB1B-003048C69416.root'),
            ('START53_V7F::All',
             '/store/mc/Summer12_DR53X/LongLivedChi0ToMuQQ_MSquark-1500_MChi-494_TuneZ2Star_8TeV-pythia6/AODSIM/PU_S10_START53_V7A-v1/00000/B69AB21F-D50F-E211-B2A7-00215E21DCA2.root'
             ),
            ('START53_V7F::All',
             '/store/mc/Summer12_DR53X/QCD_Pt-170to300_TuneZ2star_8TeV_pythia6/AODSIM/PU_S10_START53_V7A-v2/00000/0CECB336-B70E-E211-B22D-0018F3D096CE.root'
             ),
            ('FT_53_V6_AN2::All',
             '/store/data/Run2012B/JetHT/AOD/13Jul2012-v1/00000/04989369-BFD2-E111-B190-002618943933.root'
             ),
            ('FT_53_V10_AN2::All',
             '/store/data/Run2012C/JetHT/AOD/24Aug2012-v2/00000/22AFA398-D600-E211-AD2E-00261894385A.root'
             ),
            ('GR_P_V41_AN2::All',
             '/store/data/Run2012C/JetHT/AOD/PromptReco-v2/000/198/955/144C3A0C-D8CF-E111-8A6C-001D09F27067.root'
             ),
            ('GR_P_V42_AN2::All',
             '/store/data/Run2012D/JetHT/AOD/PromptReco-v1/000/204/506/F8F3A290-8C10-E211-B2D3-001D09F2525D.root'
             ),
        ],
        "53XCtrl": [
            ('START53_V7F::All',
             '/store/mc/Summer12_DR53X/HTo2LongLivedTo4F_MH-1000_MFF-350_CTau35To3500_8TeV-pythia6/AODSIM/DEBUG_PU_S10_START53_V7A-v1/0000/9C61CA50-0BF3-E111-91E1-003048D2BB22.root'
             ),
            ('START53_V7F::All',
             '/store/mc/Summer12_DR53X/QCD_Pt-170to300_TuneZ2star_8TeV_pythia6/AODSIM/PU_S10_START53_V7A-v2/00000/0CECB336-B70E-E211-B22D-0018F3D096CE.root'
             ),
            ('FT_53_V6_AN2::All',
             '/store/data/Run2012B/SingleMu/AOD/13Jul2012-v1/0000/24DF4BD2-0DD4-E111-BB00-E0CB4E4408DF.root'
             ),
            ('FT_53_V10_AN2::All',
             '/store/data/Run2012C/SingleMu/AOD/24Aug2012-v1/00000/10FD4F26-FAEE-E111-9AB5-00259073E4FC.root'
             ),
            ('GR_P_V41_AN2::All',
             '/store/data/Run2012C/SingleMu/AOD/PromptReco-v2/000/198/955/247727CE-98CF-E111-B2EF-5404A63886A0.root'
             ),
            ('GR_P_V42_AN2::All',
             '/store/data/Run2012D/SingleMu/AOD/PromptReco-v1/000/204/506/8078B6DD-9210-E211-80BE-BCAEC5329700.root'
             ),
        ]
    }["53X"][int(not options.signal) + options.isData]
    options.files = options.files if options.files else defaultFile
    options.GlobalTag = options.GlobalTag if options.GlobalTag else defaultGT

    if not options.isData:  #remove L2L3Residual correction from simulation options
        jecs = [jc for jc in options.jetCorrections if jc != 'L2L3Residual']
        options.clearList('jetCorrections')
        options.jetCorrections = jecs

    return options