FitNtuple/python/FitNtuple_cfg.py

import FWCore.ParameterSet.Config as cms
import os

realData = False
getQscale = False
getAlternating = False

WHICH = int(os.getenv("PROCESS"))

# if WHICH not in (1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111): raise Exception

if WHICH == 0:
    getQscale = True
    sample = "ppMuX_CMSSW382_v6"
elif WHICH in (1, 2, 3, 4):
    getQscale = True
    sample = "InclusiveMu5_Pt30_CMSSW382_v6"
elif WHICH in (5, 6, 7, 8):
    getQscale = True
    sample = "InclusiveMu5_Pt50_CMSSW382_v6"
elif WHICH in (9, 10):
    getQscale = True
    sample = "InclusiveMu5_Pt150_CMSSW382_v6"
elif WHICH == 11:
    getQscale = True
    sample = "InclusiveMu5_Pt250_CMSSW382_v6"
elif WHICH == 12:
    getQscale = True
    sample = "InclusiveMu5_Pt350_CMSSW382_v6"

# elif WHICH == 13: sample = "DYmumu_M0_5_CMSSW382_v6"
# elif WHICH == 14: sample = "DYmumu_M5_10_CMSSW382_v6"
# elif WHICH in (15, 16): sample = "DYmumu_M10_20_CMSSW382_v6"
# elif WHICH in (17, 18, 19, 20, 21, 22, 23, 24, 25, 26): sample = "Zmumu_M20_CMSSW382_v6"

elif WHICH in (13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26): sample = "DrellYanPythia8_CMSSW382_v6"

elif WHICH == 27: sample = "JPsiToMuMu_CMSSW382_v6"
elif WHICH == 28: sample = "Psi2SToJpsiPiPi_CMSSW382_v6"
elif WHICH == 29: sample = "Psi2SToMuMu_CMSSW382_v6"

elif WHICH in (30, 31, 32, 33, 34, 35, 36, 37, 38, 39):
    getAlternating = True
    sample = "PairGun100_mass50_alternating_CMSSW382_v6"
elif WHICH in (40, 41):
    getAlternating = True
    sample = "PairGun100_alternating_CMSSW382_v6"
elif WHICH in (42, 43):
    getAlternating = True
    sample = "PairGun100_pileup_alternating_CMSSW382_v6"
elif WHICH in (44, 45):
    getAlternating = True
    sample = "PairGun100_pileup5_alternating_CMSSW382_v6"
elif WHICH in (46, 47):
    getAlternating = True
    sample = "PairGun100_dispvert100_alternating_CMSSW382_v6"
elif WHICH in (48, 49):
    getAlternating = True
    sample = "PairPairGun100_pileup_alternating_CMSSW382_v6"

elif WHICH in (50, 51, 52): sample = "NMSSMHiggs_h100_a2_pileup_CMSSW382_v6"
elif WHICH in (53, 54, 55): sample = "EXTRAU1_DMAA_pileup_CMSSW382_v6"
elif WHICH in (56, 57, 58): sample = "EXTRAU1_DMMIXED_pileup_CMSSW382_v6"
elif WHICH in (59, 60, 61): sample = "EXTRAU1_electrons_pileup_CMSSW382_v6"
elif WHICH in (62, 63, 64): sample = "EXTRAU1_pions_pileup_CMSSW382_v6"

elif WHICH == 65: sample = "WZ_CMSSW382_v6"
elif WHICH == 66: sample = "ZZ_CMSSW382_v6"
elif WHICH == 67: sample = "WW_CMSSW382_v6"

elif WHICH in (68, 69, 70, 71, 72, 73, 74, 75): sample = "inclBB2mu_CMSSW384_v6"
elif WHICH in (76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87): sample = "inclBB4mu_CMSSW384_v6"

elif WHICH == 88: sample = "inclBB4mu_Pt50_CMSSW384_v6"
elif WHICH == 89: sample = "inclBB4mu_Pt150_CMSSW384_v6"
elif WHICH == 90: sample = "inclBB4mu_Pt250_CMSSW384_v6"
elif WHICH == 91: sample = "inclBB4mu_Pt350_CMSSW384_v6"

elif WHICH in (92, 93, 94, 95, 96, 97, 98, 99, 100, 101):
    sample = "DataSep17_CMSSW382_v6"
    realData = True
elif WHICH in (102, 103, 104, 105, 106, 107, 108, 109, 110, 111):
    sample = "Data2010BPromptAll_CMSSW384_v6"
    realData = True

else: raise Exception

process = cms.Process("FITNTUPLE")
process.load("AnalysisDataFormats.MuJetAnalysis.filelists.%s_cff" % sample)
lenFileNames = len(process.source.fileNames)

if WHICH == 1: process.source.fileNames = process.source.fileNames[lenFileNames*0/4:lenFileNames*1/4]
if WHICH == 2: process.source.fileNames = process.source.fileNames[lenFileNames*1/4:lenFileNames*2/4]
if WHICH == 3: process.source.fileNames = process.source.fileNames[lenFileNames*2/4:lenFileNames*3/4]
if WHICH == 4: process.source.fileNames = process.source.fileNames[lenFileNames*3/4:lenFileNames*4/4]

if WHICH == 5: process.source.fileNames = process.source.fileNames[lenFileNames*0/4:lenFileNames*1/4]
if WHICH == 6: process.source.fileNames = process.source.fileNames[lenFileNames*1/4:lenFileNames*2/4]
if WHICH == 7: process.source.fileNames = process.source.fileNames[lenFileNames*2/4:lenFileNames*3/4]
if WHICH == 8: process.source.fileNames = process.source.fileNames[lenFileNames*3/4:lenFileNames*4/4]

if WHICH == 9: process.source.fileNames = process.source.fileNames[lenFileNames*0/2:lenFileNames*1/2]
if WHICH == 10: process.source.fileNames = process.source.fileNames[lenFileNames*1/2:lenFileNames*2/2]

if WHICH == 13: process.source.fileNames = process.source.fileNames[lenFileNames*0/14:lenFileNames*1/14]
if WHICH == 14: process.source.fileNames = process.source.fileNames[lenFileNames*1/14:lenFileNames*2/14]
if WHICH == 15: process.source.fileNames = process.source.fileNames[lenFileNames*2/14:lenFileNames*3/14]
if WHICH == 16: process.source.fileNames = process.source.fileNames[lenFileNames*3/14:lenFileNames*4/14]
if WHICH == 17: process.source.fileNames = process.source.fileNames[lenFileNames*4/14:lenFileNames*5/14]
if WHICH == 18: process.source.fileNames = process.source.fileNames[lenFileNames*5/14:lenFileNames*6/14]
if WHICH == 19: process.source.fileNames = process.source.fileNames[lenFileNames*6/14:lenFileNames*7/14]
if WHICH == 20: process.source.fileNames = process.source.fileNames[lenFileNames*7/14:lenFileNames*8/14]
if WHICH == 21: process.source.fileNames = process.source.fileNames[lenFileNames*8/14:lenFileNames*9/14]
if WHICH == 22: process.source.fileNames = process.source.fileNames[lenFileNames*9/14:lenFileNames*10/14]
if WHICH == 23: process.source.fileNames = process.source.fileNames[lenFileNames*10/14:lenFileNames*11/14]
if WHICH == 24: process.source.fileNames = process.source.fileNames[lenFileNames*11/14:lenFileNames*12/14]
if WHICH == 25: process.source.fileNames = process.source.fileNames[lenFileNames*12/14:lenFileNames*13/14]
if WHICH == 26: process.source.fileNames = process.source.fileNames[lenFileNames*13/14:lenFileNames*14/14]

if WHICH == 30: process.source.fileNames = process.source.fileNames[lenFileNames*0/10:lenFileNames*1/10]
if WHICH == 31: process.source.fileNames = process.source.fileNames[lenFileNames*1/10:lenFileNames*2/10]
if WHICH == 32: process.source.fileNames = process.source.fileNames[lenFileNames*2/10:lenFileNames*3/10]
if WHICH == 33: process.source.fileNames = process.source.fileNames[lenFileNames*3/10:lenFileNames*4/10]
if WHICH == 34: process.source.fileNames = process.source.fileNames[lenFileNames*4/10:lenFileNames*5/10]
if WHICH == 35: process.source.fileNames = process.source.fileNames[lenFileNames*5/10:lenFileNames*6/10]
if WHICH == 36: process.source.fileNames = process.source.fileNames[lenFileNames*6/10:lenFileNames*7/10]
if WHICH == 37: process.source.fileNames = process.source.fileNames[lenFileNames*7/10:lenFileNames*8/10]
if WHICH == 38: process.source.fileNames = process.source.fileNames[lenFileNames*8/10:lenFileNames*9/10]
if WHICH == 39: process.source.fileNames = process.source.fileNames[lenFileNames*9/10:lenFileNames*10/10]

if WHICH == 40: process.source.fileNames = process.source.fileNames[lenFileNames*0/2:lenFileNames*1/2]
if WHICH == 41: process.source.fileNames = process.source.fileNames[lenFileNames*1/2:lenFileNames*2/2]

if WHICH == 42: process.source.fileNames = process.source.fileNames[lenFileNames*0/2:lenFileNames*1/2]
if WHICH == 43: process.source.fileNames = process.source.fileNames[lenFileNames*1/2:lenFileNames*2/2]

if WHICH == 44: process.source.fileNames = process.source.fileNames[lenFileNames*0/2:lenFileNames*1/2]
if WHICH == 45: process.source.fileNames = process.source.fileNames[lenFileNames*1/2:lenFileNames*2/2]

if WHICH == 46: process.source.fileNames = process.source.fileNames[lenFileNames*0/2:lenFileNames*1/2]
if WHICH == 47: process.source.fileNames = process.source.fileNames[lenFileNames*1/2:lenFileNames*2/2]

if WHICH == 48: process.source.fileNames = process.source.fileNames[lenFileNames*0/2:lenFileNames*1/2]
if WHICH == 49: process.source.fileNames = process.source.fileNames[lenFileNames*1/2:lenFileNames*2/2]

if WHICH == 50: process.source.fileNames = process.source.fileNames[lenFileNames*0/3:lenFileNames*1/3]
if WHICH == 51: process.source.fileNames = process.source.fileNames[lenFileNames*1/3:lenFileNames*2/3]
if WHICH == 52: process.source.fileNames = process.source.fileNames[lenFileNames*2/3:lenFileNames*3/3]

if WHICH == 53: process.source.fileNames = process.source.fileNames[lenFileNames*0/3:lenFileNames*1/3]
if WHICH == 54: process.source.fileNames = process.source.fileNames[lenFileNames*1/3:lenFileNames*2/3]
if WHICH == 55: process.source.fileNames = process.source.fileNames[lenFileNames*2/3:lenFileNames*3/3]

if WHICH == 56: process.source.fileNames = process.source.fileNames[lenFileNames*0/3:lenFileNames*1/3]
if WHICH == 57: process.source.fileNames = process.source.fileNames[lenFileNames*1/3:lenFileNames*2/3]
if WHICH == 58: process.source.fileNames = process.source.fileNames[lenFileNames*2/3:lenFileNames*3/3]

if WHICH == 59: process.source.fileNames = process.source.fileNames[lenFileNames*0/3:lenFileNames*1/3]
if WHICH == 60: process.source.fileNames = process.source.fileNames[lenFileNames*1/3:lenFileNames*2/3]
if WHICH == 61: process.source.fileNames = process.source.fileNames[lenFileNames*2/3:lenFileNames*3/3]

if WHICH == 62: process.source.fileNames = process.source.fileNames[lenFileNames*0/3:lenFileNames*1/3]
if WHICH == 63: process.source.fileNames = process.source.fileNames[lenFileNames*1/3:lenFileNames*2/3]
if WHICH == 64: process.source.fileNames = process.source.fileNames[lenFileNames*2/3:lenFileNames*3/3]

if WHICH == 68: process.source.fileNames = process.source.fileNames[lenFileNames*0/8:lenFileNames*1/8]
if WHICH == 69: process.source.fileNames = process.source.fileNames[lenFileNames*1/8:lenFileNames*2/8]
if WHICH == 70: process.source.fileNames = process.source.fileNames[lenFileNames*2/8:lenFileNames*3/8]
if WHICH == 71: process.source.fileNames = process.source.fileNames[lenFileNames*3/8:lenFileNames*4/8]
if WHICH == 72: process.source.fileNames = process.source.fileNames[lenFileNames*4/8:lenFileNames*5/8]
if WHICH == 73: process.source.fileNames = process.source.fileNames[lenFileNames*5/8:lenFileNames*6/8]
if WHICH == 74: process.source.fileNames = process.source.fileNames[lenFileNames*6/8:lenFileNames*7/8]
if WHICH == 75: process.source.fileNames = process.source.fileNames[lenFileNames*7/8:lenFileNames*8/8]

if WHICH == 76: process.source.fileNames = process.source.fileNames[lenFileNames*0/12:lenFileNames*1/12]
if WHICH == 77: process.source.fileNames = process.source.fileNames[lenFileNames*1/12:lenFileNames*2/12]
if WHICH == 78: process.source.fileNames = process.source.fileNames[lenFileNames*2/12:lenFileNames*3/12]
if WHICH == 79: process.source.fileNames = process.source.fileNames[lenFileNames*3/12:lenFileNames*4/12]
if WHICH == 80: process.source.fileNames = process.source.fileNames[lenFileNames*4/12:lenFileNames*5/12]
if WHICH == 81: process.source.fileNames = process.source.fileNames[lenFileNames*5/12:lenFileNames*6/12]
if WHICH == 82: process.source.fileNames = process.source.fileNames[lenFileNames*6/12:lenFileNames*7/12]
if WHICH == 83: process.source.fileNames = process.source.fileNames[lenFileNames*7/12:lenFileNames*8/12]
if WHICH == 84: process.source.fileNames = process.source.fileNames[lenFileNames*8/12:lenFileNames*9/12]
if WHICH == 85: process.source.fileNames = process.source.fileNames[lenFileNames*9/12:lenFileNames*10/12]
if WHICH == 86: process.source.fileNames = process.source.fileNames[lenFileNames*10/12:lenFileNames*11/12]
if WHICH == 87: process.source.fileNames = process.source.fileNames[lenFileNames*11/12:lenFileNames*12/12]

if WHICH == 92: process.source.fileNames = process.source.fileNames[lenFileNames*0/10:lenFileNames*1/10]
if WHICH == 93: process.source.fileNames = process.source.fileNames[lenFileNames*1/10:lenFileNames*2/10]
if WHICH == 94: process.source.fileNames = process.source.fileNames[lenFileNames*2/10:lenFileNames*3/10]
if WHICH == 95: process.source.fileNames = process.source.fileNames[lenFileNames*3/10:lenFileNames*4/10]
if WHICH == 96: process.source.fileNames = process.source.fileNames[lenFileNames*4/10:lenFileNames*5/10]
if WHICH == 97: process.source.fileNames = process.source.fileNames[lenFileNames*5/10:lenFileNames*6/10]
if WHICH == 98: process.source.fileNames = process.source.fileNames[lenFileNames*6/10:lenFileNames*7/10]
if WHICH == 99: process.source.fileNames = process.source.fileNames[lenFileNames*7/10:lenFileNames*8/10]
if WHICH == 100: process.source.fileNames = process.source.fileNames[lenFileNames*8/10:lenFileNames*9/10]
if WHICH == 101: process.source.fileNames = process.source.fileNames[lenFileNames*9/10:lenFileNames*10/10]

if WHICH == 102: process.source.fileNames = process.source.fileNames[lenFileNames*0/10:lenFileNames*1/10]
if WHICH == 103: process.source.fileNames = process.source.fileNames[lenFileNames*1/10:lenFileNames*2/10]
if WHICH == 104: process.source.fileNames = process.source.fileNames[lenFileNames*2/10:lenFileNames*3/10]
if WHICH == 105: process.source.fileNames = process.source.fileNames[lenFileNames*3/10:lenFileNames*4/10]
if WHICH == 106: process.source.fileNames = process.source.fileNames[lenFileNames*4/10:lenFileNames*5/10]
if WHICH == 107: process.source.fileNames = process.source.fileNames[lenFileNames*5/10:lenFileNames*6/10]
if WHICH == 108: process.source.fileNames = process.source.fileNames[lenFileNames*6/10:lenFileNames*7/10]
if WHICH == 109: process.source.fileNames = process.source.fileNames[lenFileNames*7/10:lenFileNames*8/10]
if WHICH == 110: process.source.fileNames = process.source.fileNames[lenFileNames*8/10:lenFileNames*9/10]
if WHICH == 111: process.source.fileNames = process.source.fileNames[lenFileNames*9/10:lenFileNames*10/10]

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

process.MessageLogger = cms.Service("MessageLogger", destinations = cms.untracked.vstring("cout"), cout = cms.untracked.PSet(threshold = cms.untracked.string("ERROR")))
process.load("Configuration/StandardSequences/FrontierConditions_GlobalTag_cff")
process.GlobalTag.globaltag = "MC_38Y_V12::All"

process.NoScraping = cms.EDFilter("FilterOutScraping",
                                  applyfilter = cms.untracked.bool(True),
                                  debugOn = cms.untracked.bool(False),
                                  numtrack = cms.untracked.uint32(10),
                                  thresh = cms.untracked.double(0.25)
                                  )

process.load("RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi")
process.OneGoodVertexFilter = cms.EDFilter("VertexSelector",
                                           src = cms.InputTag("offlinePrimaryVertices"),
                                           cut = cms.string("!isFake && tracksSize > 3 && abs(z) <= 24"),
                                           filter = cms.bool(True),
                                           )

process.load("AnalysisAlgos.MuJetProducer.MuJetProducer_cff")
process.MuJetProducer.muons = cms.InputTag("cleanPatMuonsTriggerMatch")

process.FitNtuple = cms.EDAnalyzer("FitNtuple",
                                   src = cms.InputTag("MuJetProducer"),
                                   srcOrphans = cms.InputTag("MuJetProducer", "Orphans"),
                                   getQscale = cms.bool(getQscale),
                                   getAlternating = cms.bool(getAlternating),
                                   )

if getAlternating:
    process.Path = cms.Path(process.MuJetProducer * process.FitNtuple)
else:
    process.Path = cms.Path(process.NoScraping * process.OneGoodVertexFilter * process.MuJetProducer * process.FitNtuple)

process.TFileService = cms.Service("TFileService", fileName = cms.string("/uscmst1b_scratch/lpc1/3DayLifetime/pivarski/FITNTUPLE05/FitNtuple_%s_%d.root" % (sample, WHICH)))
Revision:	1.6
Committed:	Wed Dec 1 00:32:28 2010 UTC (14 years, 5 months ago) by pivarski
Content type:	text/x-python
Branch:	MAIN
CVS Tags:	JP-2010-11-30-a
Changes since 1.5:	+3 -2 lines
Log Message:	dimuorphan to get a mass template for un-triggered muons
#	User	Rev	Content
1	pivarski	1.1	import FWCore.ParameterSet.Config as cms
2			import os
3
4			realData = False
5			getQscale = False
6			getAlternating = False
7
8			WHICH = int(os.getenv("PROCESS"))
9
10	pivarski	1.6	# if WHICH not in (1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111): raise Exception
11	pivarski	1.4
12	pivarski	1.1	if WHICH == 0:
13			getQscale = True
14			sample = "ppMuX_CMSSW382_v6"
15			elif WHICH in (1, 2, 3, 4):
16			getQscale = True
17			sample = "InclusiveMu5_Pt30_CMSSW382_v6"
18			elif WHICH in (5, 6, 7, 8):
19			getQscale = True
20			sample = "InclusiveMu5_Pt50_CMSSW382_v6"
21			elif WHICH in (9, 10):
22			getQscale = True
23			sample = "InclusiveMu5_Pt150_CMSSW382_v6"
24			elif WHICH == 11:
25			getQscale = True
26			sample = "InclusiveMu5_Pt250_CMSSW382_v6"
27			elif WHICH == 12:
28			getQscale = True
29			sample = "InclusiveMu5_Pt350_CMSSW382_v6"
30
31	pivarski	1.4	# elif WHICH == 13: sample = "DYmumu_M0_5_CMSSW382_v6"
32			# elif WHICH == 14: sample = "DYmumu_M5_10_CMSSW382_v6"
33			# elif WHICH in (15, 16): sample = "DYmumu_M10_20_CMSSW382_v6"
34			# elif WHICH in (17, 18, 19, 20, 21, 22, 23, 24, 25, 26): sample = "Zmumu_M20_CMSSW382_v6"
35
36			elif WHICH in (13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26): sample = "DrellYanPythia8_CMSSW382_v6"
37	pivarski	1.1
38			elif WHICH == 27: sample = "JPsiToMuMu_CMSSW382_v6"
39			elif WHICH == 28: sample = "Psi2SToJpsiPiPi_CMSSW382_v6"
40			elif WHICH == 29: sample = "Psi2SToMuMu_CMSSW382_v6"
41
42			elif WHICH in (30, 31, 32, 33, 34, 35, 36, 37, 38, 39):
43			getAlternating = True
44			sample = "PairGun100_mass50_alternating_CMSSW382_v6"
45			elif WHICH in (40, 41):
46			getAlternating = True
47			sample = "PairGun100_alternating_CMSSW382_v6"
48			elif WHICH in (42, 43):
49			getAlternating = True
50			sample = "PairGun100_pileup_alternating_CMSSW382_v6"
51			elif WHICH in (44, 45):
52			getAlternating = True
53			sample = "PairGun100_pileup5_alternating_CMSSW382_v6"
54			elif WHICH in (46, 47):
55			getAlternating = True
56			sample = "PairGun100_dispvert100_alternating_CMSSW382_v6"
57			elif WHICH in (48, 49):
58			getAlternating = True
59			sample = "PairPairGun100_pileup_alternating_CMSSW382_v6"
60
61			elif WHICH in (50, 51, 52): sample = "NMSSMHiggs_h100_a2_pileup_CMSSW382_v6"
62			elif WHICH in (53, 54, 55): sample = "EXTRAU1_DMAA_pileup_CMSSW382_v6"
63			elif WHICH in (56, 57, 58): sample = "EXTRAU1_DMMIXED_pileup_CMSSW382_v6"
64			elif WHICH in (59, 60, 61): sample = "EXTRAU1_electrons_pileup_CMSSW382_v6"
65			elif WHICH in (62, 63, 64): sample = "EXTRAU1_pions_pileup_CMSSW382_v6"
66
67			elif WHICH == 65: sample = "WZ_CMSSW382_v6"
68			elif WHICH == 66: sample = "ZZ_CMSSW382_v6"
69			elif WHICH == 67: sample = "WW_CMSSW382_v6"
70
71			elif WHICH in (68, 69, 70, 71, 72, 73, 74, 75): sample = "inclBB2mu_CMSSW384_v6"
72			elif WHICH in (76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87): sample = "inclBB4mu_CMSSW384_v6"
73
74			elif WHICH == 88: sample = "inclBB4mu_Pt50_CMSSW384_v6"
75			elif WHICH == 89: sample = "inclBB4mu_Pt150_CMSSW384_v6"
76			elif WHICH == 90: sample = "inclBB4mu_Pt250_CMSSW384_v6"
77			elif WHICH == 91: sample = "inclBB4mu_Pt350_CMSSW384_v6"
78
79			elif WHICH in (92, 93, 94, 95, 96, 97, 98, 99, 100, 101):
80			sample = "DataSep17_CMSSW382_v6"
81			realData = True
82			elif WHICH in (102, 103, 104, 105, 106, 107, 108, 109, 110, 111):
83			sample = "Data2010BPromptAll_CMSSW384_v6"
84			realData = True
85
86			else: raise Exception
87
88			process = cms.Process("FITNTUPLE")
89			process.load("AnalysisDataFormats.MuJetAnalysis.filelists.%s_cff" % sample)
90			lenFileNames = len(process.source.fileNames)
91
92			if WHICH == 1: process.source.fileNames = process.source.fileNames[lenFileNames0/4:lenFileNames1/4]
93			if WHICH == 2: process.source.fileNames = process.source.fileNames[lenFileNames1/4:lenFileNames2/4]
94			if WHICH == 3: process.source.fileNames = process.source.fileNames[lenFileNames2/4:lenFileNames3/4]
95			if WHICH == 4: process.source.fileNames = process.source.fileNames[lenFileNames3/4:lenFileNames4/4]
96
97			if WHICH == 5: process.source.fileNames = process.source.fileNames[lenFileNames0/4:lenFileNames1/4]
98			if WHICH == 6: process.source.fileNames = process.source.fileNames[lenFileNames1/4:lenFileNames2/4]
99			if WHICH == 7: process.source.fileNames = process.source.fileNames[lenFileNames2/4:lenFileNames3/4]
100			if WHICH == 8: process.source.fileNames = process.source.fileNames[lenFileNames3/4:lenFileNames4/4]
101
102			if WHICH == 9: process.source.fileNames = process.source.fileNames[lenFileNames0/2:lenFileNames1/2]
103			if WHICH == 10: process.source.fileNames = process.source.fileNames[lenFileNames1/2:lenFileNames2/2]
104
105	pivarski	1.4	if WHICH == 13: process.source.fileNames = process.source.fileNames[lenFileNames0/14:lenFileNames1/14]
106			if WHICH == 14: process.source.fileNames = process.source.fileNames[lenFileNames1/14:lenFileNames2/14]
107			if WHICH == 15: process.source.fileNames = process.source.fileNames[lenFileNames2/14:lenFileNames3/14]
108			if WHICH == 16: process.source.fileNames = process.source.fileNames[lenFileNames3/14:lenFileNames4/14]
109			if WHICH == 17: process.source.fileNames = process.source.fileNames[lenFileNames4/14:lenFileNames5/14]
110			if WHICH == 18: process.source.fileNames = process.source.fileNames[lenFileNames5/14:lenFileNames6/14]
111			if WHICH == 19: process.source.fileNames = process.source.fileNames[lenFileNames6/14:lenFileNames7/14]
112			if WHICH == 20: process.source.fileNames = process.source.fileNames[lenFileNames7/14:lenFileNames8/14]
113			if WHICH == 21: process.source.fileNames = process.source.fileNames[lenFileNames8/14:lenFileNames9/14]
114			if WHICH == 22: process.source.fileNames = process.source.fileNames[lenFileNames9/14:lenFileNames10/14]
115			if WHICH == 23: process.source.fileNames = process.source.fileNames[lenFileNames10/14:lenFileNames11/14]
116			if WHICH == 24: process.source.fileNames = process.source.fileNames[lenFileNames11/14:lenFileNames12/14]
117			if WHICH == 25: process.source.fileNames = process.source.fileNames[lenFileNames12/14:lenFileNames13/14]
118			if WHICH == 26: process.source.fileNames = process.source.fileNames[lenFileNames13/14:lenFileNames14/14]
119	pivarski	1.1
120			if WHICH == 30: process.source.fileNames = process.source.fileNames[lenFileNames0/10:lenFileNames1/10]
121			if WHICH == 31: process.source.fileNames = process.source.fileNames[lenFileNames1/10:lenFileNames2/10]
122			if WHICH == 32: process.source.fileNames = process.source.fileNames[lenFileNames2/10:lenFileNames3/10]
123			if WHICH == 33: process.source.fileNames = process.source.fileNames[lenFileNames3/10:lenFileNames4/10]
124			if WHICH == 34: process.source.fileNames = process.source.fileNames[lenFileNames4/10:lenFileNames5/10]
125			if WHICH == 35: process.source.fileNames = process.source.fileNames[lenFileNames5/10:lenFileNames6/10]
126			if WHICH == 36: process.source.fileNames = process.source.fileNames[lenFileNames6/10:lenFileNames7/10]
127			if WHICH == 37: process.source.fileNames = process.source.fileNames[lenFileNames7/10:lenFileNames8/10]
128			if WHICH == 38: process.source.fileNames = process.source.fileNames[lenFileNames8/10:lenFileNames9/10]
129			if WHICH == 39: process.source.fileNames = process.source.fileNames[lenFileNames9/10:lenFileNames10/10]
130
131			if WHICH == 40: process.source.fileNames = process.source.fileNames[lenFileNames0/2:lenFileNames1/2]
132			if WHICH == 41: process.source.fileNames = process.source.fileNames[lenFileNames1/2:lenFileNames2/2]
133
134			if WHICH == 42: process.source.fileNames = process.source.fileNames[lenFileNames0/2:lenFileNames1/2]
135			if WHICH == 43: process.source.fileNames = process.source.fileNames[lenFileNames1/2:lenFileNames2/2]
136
137			if WHICH == 44: process.source.fileNames = process.source.fileNames[lenFileNames0/2:lenFileNames1/2]
138			if WHICH == 45: process.source.fileNames = process.source.fileNames[lenFileNames1/2:lenFileNames2/2]
139
140			if WHICH == 46: process.source.fileNames = process.source.fileNames[lenFileNames0/2:lenFileNames1/2]
141			if WHICH == 47: process.source.fileNames = process.source.fileNames[lenFileNames1/2:lenFileNames2/2]
142
143			if WHICH == 48: process.source.fileNames = process.source.fileNames[lenFileNames0/2:lenFileNames1/2]
144			if WHICH == 49: process.source.fileNames = process.source.fileNames[lenFileNames1/2:lenFileNames2/2]
145
146			if WHICH == 50: process.source.fileNames = process.source.fileNames[lenFileNames0/3:lenFileNames1/3]
147			if WHICH == 51: process.source.fileNames = process.source.fileNames[lenFileNames1/3:lenFileNames2/3]
148			if WHICH == 52: process.source.fileNames = process.source.fileNames[lenFileNames2/3:lenFileNames3/3]
149
150			if WHICH == 53: process.source.fileNames = process.source.fileNames[lenFileNames0/3:lenFileNames1/3]
151			if WHICH == 54: process.source.fileNames = process.source.fileNames[lenFileNames1/3:lenFileNames2/3]
152			if WHICH == 55: process.source.fileNames = process.source.fileNames[lenFileNames2/3:lenFileNames3/3]
153
154			if WHICH == 56: process.source.fileNames = process.source.fileNames[lenFileNames0/3:lenFileNames1/3]
155			if WHICH == 57: process.source.fileNames = process.source.fileNames[lenFileNames1/3:lenFileNames2/3]
156			if WHICH == 58: process.source.fileNames = process.source.fileNames[lenFileNames2/3:lenFileNames3/3]
157
158			if WHICH == 59: process.source.fileNames = process.source.fileNames[lenFileNames0/3:lenFileNames1/3]
159			if WHICH == 60: process.source.fileNames = process.source.fileNames[lenFileNames1/3:lenFileNames2/3]
160			if WHICH == 61: process.source.fileNames = process.source.fileNames[lenFileNames2/3:lenFileNames3/3]
161
162			if WHICH == 62: process.source.fileNames = process.source.fileNames[lenFileNames0/3:lenFileNames1/3]
163			if WHICH == 63: process.source.fileNames = process.source.fileNames[lenFileNames1/3:lenFileNames2/3]
164			if WHICH == 64: process.source.fileNames = process.source.fileNames[lenFileNames2/3:lenFileNames3/3]
165
166			if WHICH == 68: process.source.fileNames = process.source.fileNames[lenFileNames0/8:lenFileNames1/8]
167			if WHICH == 69: process.source.fileNames = process.source.fileNames[lenFileNames1/8:lenFileNames2/8]
168			if WHICH == 70: process.source.fileNames = process.source.fileNames[lenFileNames2/8:lenFileNames3/8]
169			if WHICH == 71: process.source.fileNames = process.source.fileNames[lenFileNames3/8:lenFileNames4/8]
170			if WHICH == 72: process.source.fileNames = process.source.fileNames[lenFileNames4/8:lenFileNames5/8]
171			if WHICH == 73: process.source.fileNames = process.source.fileNames[lenFileNames5/8:lenFileNames6/8]
172			if WHICH == 74: process.source.fileNames = process.source.fileNames[lenFileNames6/8:lenFileNames7/8]
173			if WHICH == 75: process.source.fileNames = process.source.fileNames[lenFileNames7/8:lenFileNames8/8]
174
175			if WHICH == 76: process.source.fileNames = process.source.fileNames[lenFileNames0/12:lenFileNames1/12]
176			if WHICH == 77: process.source.fileNames = process.source.fileNames[lenFileNames1/12:lenFileNames2/12]
177			if WHICH == 78: process.source.fileNames = process.source.fileNames[lenFileNames2/12:lenFileNames3/12]
178			if WHICH == 79: process.source.fileNames = process.source.fileNames[lenFileNames3/12:lenFileNames4/12]
179			if WHICH == 80: process.source.fileNames = process.source.fileNames[lenFileNames4/12:lenFileNames5/12]
180			if WHICH == 81: process.source.fileNames = process.source.fileNames[lenFileNames5/12:lenFileNames6/12]
181			if WHICH == 82: process.source.fileNames = process.source.fileNames[lenFileNames6/12:lenFileNames7/12]
182			if WHICH == 83: process.source.fileNames = process.source.fileNames[lenFileNames7/12:lenFileNames8/12]
183			if WHICH == 84: process.source.fileNames = process.source.fileNames[lenFileNames8/12:lenFileNames9/12]
184			if WHICH == 85: process.source.fileNames = process.source.fileNames[lenFileNames9/12:lenFileNames10/12]
185			if WHICH == 86: process.source.fileNames = process.source.fileNames[lenFileNames10/12:lenFileNames11/12]
186			if WHICH == 87: process.source.fileNames = process.source.fileNames[lenFileNames11/12:lenFileNames12/12]
187
188			if WHICH == 92: process.source.fileNames = process.source.fileNames[lenFileNames0/10:lenFileNames1/10]
189			if WHICH == 93: process.source.fileNames = process.source.fileNames[lenFileNames1/10:lenFileNames2/10]
190			if WHICH == 94: process.source.fileNames = process.source.fileNames[lenFileNames2/10:lenFileNames3/10]
191			if WHICH == 95: process.source.fileNames = process.source.fileNames[lenFileNames3/10:lenFileNames4/10]
192			if WHICH == 96: process.source.fileNames = process.source.fileNames[lenFileNames4/10:lenFileNames5/10]
193			if WHICH == 97: process.source.fileNames = process.source.fileNames[lenFileNames5/10:lenFileNames6/10]
194			if WHICH == 98: process.source.fileNames = process.source.fileNames[lenFileNames6/10:lenFileNames7/10]
195			if WHICH == 99: process.source.fileNames = process.source.fileNames[lenFileNames7/10:lenFileNames8/10]
196			if WHICH == 100: process.source.fileNames = process.source.fileNames[lenFileNames8/10:lenFileNames9/10]
197			if WHICH == 101: process.source.fileNames = process.source.fileNames[lenFileNames9/10:lenFileNames10/10]
198
199			if WHICH == 102: process.source.fileNames = process.source.fileNames[lenFileNames0/10:lenFileNames1/10]
200			if WHICH == 103: process.source.fileNames = process.source.fileNames[lenFileNames1/10:lenFileNames2/10]
201			if WHICH == 104: process.source.fileNames = process.source.fileNames[lenFileNames2/10:lenFileNames3/10]
202			if WHICH == 105: process.source.fileNames = process.source.fileNames[lenFileNames3/10:lenFileNames4/10]
203			if WHICH == 106: process.source.fileNames = process.source.fileNames[lenFileNames4/10:lenFileNames5/10]
204			if WHICH == 107: process.source.fileNames = process.source.fileNames[lenFileNames5/10:lenFileNames6/10]
205			if WHICH == 108: process.source.fileNames = process.source.fileNames[lenFileNames6/10:lenFileNames7/10]
206			if WHICH == 109: process.source.fileNames = process.source.fileNames[lenFileNames7/10:lenFileNames8/10]
207			if WHICH == 110: process.source.fileNames = process.source.fileNames[lenFileNames8/10:lenFileNames9/10]
208			if WHICH == 111: process.source.fileNames = process.source.fileNames[lenFileNames9/10:lenFileNames10/10]
209
210			process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(-1))
211
212	pivarski	1.2	process.MessageLogger = cms.Service("MessageLogger", destinations = cms.untracked.vstring("cout"), cout = cms.untracked.PSet(threshold = cms.untracked.string("ERROR")))
213	pivarski	1.1	process.load("Configuration/StandardSequences/FrontierConditions_GlobalTag_cff")
214			process.GlobalTag.globaltag = "MC_38Y_V12::All"
215
216	pivarski	1.2	process.NoScraping = cms.EDFilter("FilterOutScraping",
217			applyfilter = cms.untracked.bool(True),
218			debugOn = cms.untracked.bool(False),
219			numtrack = cms.untracked.uint32(10),
220			thresh = cms.untracked.double(0.25)
221			)
222	pivarski	1.1
223	pivarski	1.2	process.load("RecoVertex.PrimaryVertexProducer.OfflinePrimaryVertices_cfi")
224			process.OneGoodVertexFilter = cms.EDFilter("VertexSelector",
225			src = cms.InputTag("offlinePrimaryVertices"),
226			cut = cms.string("!isFake && tracksSize > 3 && abs(z) <= 24"),
227			filter = cms.bool(True),
228			)
229	pivarski	1.1
230			process.load("AnalysisAlgos.MuJetProducer.MuJetProducer_cff")
231			process.MuJetProducer.muons = cms.InputTag("cleanPatMuonsTriggerMatch")
232
233			process.FitNtuple = cms.EDAnalyzer("FitNtuple",
234	pivarski	1.2	src = cms.InputTag("MuJetProducer"),
235	pivarski	1.6	srcOrphans = cms.InputTag("MuJetProducer", "Orphans"),
236	pivarski	1.1	getQscale = cms.bool(getQscale),
237			getAlternating = cms.bool(getAlternating),
238			)
239
240	pivarski	1.4	if getAlternating:
241			process.Path = cms.Path(process.MuJetProducer * process.FitNtuple)
242			else:
243			process.Path = cms.Path(process.NoScraping * process.OneGoodVertexFilter * process.MuJetProducer * process.FitNtuple)
244
245	pivarski	1.6	process.TFileService = cms.Service("TFileService", fileName = cms.string("/uscmst1b_scratch/lpc1/3DayLifetime/pivarski/FITNTUPLE05/FitNtuple_%s_%d.root" % (sample, WHICH)))