fileNames = cms.untracked.vstring('file:AOD_data.root')
)

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


process.TFileService=cms.Service(
    "TFileService",
    fileName=cms.string('testoutput_thrust.root')
)

process.ThrustAKPu4PF = cms.EDAnalyzer(
    'ThrustAnalyzer',
    JetType = cms.untracked.string('pf'),
    UEAlgo = cms.untracked.string('Pu'),
    radius = cms.int(4),
    src = cms.InputTag("akVs4PFJets"),
    #centralitycollection = cms.InputTag("hiCentrality"),
    #centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
    #JetCorrections = cms.string(""),
    #PFcands = cms.InputTag("particleFlowTmp"),
    RThreshold = cms.double(0.3),  
    mRecoJetPtThreshold = cms.double(10),        
    mReco_SubJetPtThreshold = cms.double(10)        
)


process.ThrustAKVs4PF = process.ThrustAKVs4PF.clone(
    UEAlgo = cms.untracked.string('Pu')
)
import FWCore.ParameterSet.Config as cms

phoJetNtuplizer = cms.EDAnalyzer('phoJetNtuplizer',
                                  debug                                = cms.bool(False),
                                  runYear                              = cms.int(2018),
                                  saveAll                              = cms.bool(False),
                                  is_Data                              = cms.bool(True),
                                  runEventInfo                         = cms.bool(True),
                                  runPhotons                           = cms.bool(True),
                                  runJets                              = cms.bool(True),
                                  runak8Jets                           = cms.bool(False),
                                  runJetWidthCalculator                = cms.bool(False),
                                  runEle                               = cms.bool(True),
                                  runMuon                              = cms.bool(True),
                                  runTaus                              = cms.bool(True),
                                  runMet                               = cms.bool(True),
				  runGenInfo                           = cms.bool(False),
				  addPSWeights                         = cms.bool(False),

				  rhoToken                             = cms.InputTag("fixedGridRhoFastjetAll"),
				  rhoCentralToken                      = cms.InputTag("fixedGridRhoFastjetCentralNeutral"),
				  vtxToken                             = cms.InputTag("offlineSlimmedPrimaryVertices"),

				  triggerResults                       = cms.InputTag("TriggerResults", "", "HLT"),
				  patTriggerResults                    = cms.InputTag("TriggerResults", "", "PAT"),
				  recoTriggerResults                   = cms.InputTag("TriggerResults", "", "RECO"),
				  trgFilterDeltaPtCut                  = cms.double(0.5),
				  trgFilterDeltaRCut                   = cms.double(0.3),
				  triggerEvent                         = cms.InputTag("slimmedPatTrigger", "", ""),

## Tau ID
from phoJetAnalysis.phoJetNtuplizer.runTauIdMVA import *
na = TauIDEmbedder(
    process,
    cms,  # pass tour process object
    debug=True,
    toKeep=[
        "2017v2"
    ]  # pick the one you need: ["2017v1", "2017v2", "newDM2017v2", "dR0p32017v2", "2016v1", "newDM2016v1"]
)
na.runTauID()

### Analyzer Related
process.load("phoJetAnalysis.phoJetNtuplizer.phoJetNtuplizer_cfi")
process.phoJetNtuplizer.debug = cms.bool(False)
process.phoJetNtuplizer.runYear = cms.int(2017)
process.phoJetNtuplizer.runak8Jets = cms.bool(False)
process.phoJetNtuplizer.runJetWidthCalculator = cms.bool(True)
# needed for monoZprime Analysis [Valid only if runJets is True]
process.phoJetNtuplizer.jetsAK4Token = cms.InputTag(
    "selectedUpdatedPatJetsUpdatedJECAK4")
#process.phoJetNtuplizer.jetsAK8Token = cms.InputTag("selectedUpdatedPatJetsUpdatedJECAK8")
#process.phoJetNtuplizer.jetsAK8TagToken = cms.InputTag("selectedUpdatedPatJetsAK8WithDeepTags")
process.phoJetNtuplizer.pfmetToken = cms.InputTag("slimmedMETsModifiedMET")

process.p = cms.Path(process.ecalBadCalibReducedMINIAODFilter *
                     process.fullPatMetSequenceModifiedMET *
                     process.egammaPostRecoSeq * process.prefiringweight *
                     process.rerunMvaIsolationSequence *
                     process.NewTauIDsEmbedded * process.phoJetNtuplizer)

#Array axis m
tracktriggerPhase2Params.MAxisMin=cms.vint32(-480, -440, -400, -360, -320, -280, -240, -200, 
-160, -120, -80, -40, 0, 40, 80, 120, 160, 200, 240, 280, 320, 360, 400, 440, 480)

tracktriggerPhase2Params.MAxisMax=cms.vint32(-480, -440, -400, -360, -320, -280, -240, -200, 
-160, -120, -80, -40, 0, 40, 80, 120, 160, 200, 240, 280, 320, 360, 400, 440, 480)

tracktriggerPhase2Params.MAxisAv=cms.vint32(-480, -440, -400, -360, -320, -280, -240, -200, 
-160, -120, -80, -40, 0, 40, 80, 120, 160, 200, 240, 280, 320, 360, 400, 440, 480)

#Array axis c
tracktriggerPhase2Params.CAxisMin=cms.vint32(0, 80000, 160000, 240000, 320000, 400000, 480000, 
560000, 640000, 720000, 800000, 880000, 960000, 1040000, 1120000, 1200000, 1280000,
1360000, 1440000, 1520000, 1600000, 1680000, 1760000, 1840000, 1920000)

tracktriggerPhase2Params.CAxisMax=cms.vint32(82000, 162000, 242000, 322000, 402000, 482000, 
562000, 642000, 722000, 802000, 962000, 1042000, 1122000, 1202000, 1282000, 1362000,
1442000, 1522000, 1602000, 1682000, 1762000, 1842000, 1922000, 2000000)

tracktriggerPhase2Params.CAxisAv=cms.vint32(82000, 162000, 242000, 322000, 402000, 482000, 
562000, 642000, 722000, 802000, 962000, 1042000, 1122000, 1202000, 1282000, 1362000,
1442000, 1522000, 1602000, 1682000, 1762000, 1842000, 1922000, 2000000)

#number of boards in parallel
tracktriggerPhase2Params.NBoards=cms.int(2)

#Radial origin for Hough transform (cm)
tracktriggerPhase2Params.T=cms.int(45)
    # replace 'myfile.root' with the source file you want to use                                                                                                                                                                              
    #fileNames = cms.untracked.vstring("/store/user/knash/JetHT/crab_Run2015D-PromptReco-v3_Oct9_v74x_V7_25ns_v2/151013_175415/0000/B2GEDMNtuple_1.root")
    #fileNames = cms.untracked.vstring("/store/user/jdolen/B2GAnaFW/JetHT/Run2015D-05Oct2015-v1/JetHT/B2GAnaFW_v74x_V8p2_25ns/151023_194322/0000/B2GEDMNtuple_1.root")
     fileNames = cms.untracked.vstring("/store/user/jpilot/ZprimeToTT_M-1000_W-100_TuneCUETP8M1_13TeV-madgraphMLM-pythia8/B2GAnaFW_v74x_V6p1_25ns_RunIISpring15DR74-Asympt25ns_MCRUN2_74_V9-v2/151012_125345/0000/B2GEDMNtuple_5.root")
    #fileNames = files
)

###############################################################################
#PRODUCER
process.b2gntuple = cms.EDAnalyzer('ZprimeB2Ganalyzer',
                                   #testSrc = cms.string('hola')
                                   negativeWeights = cms.bool(True),
                                   isMC = cms.bool(True),
                                   isFlat = cms.bool(False),
                                   applyFilters = cms.bool(True),
                                   JERshift = cms.int(0),
                                   JECshift = cms.int(-1),
                                   reweightTopPt = cms.bool(False),
                                   puFile = cms.string('/uscms/home/camclean/nobackup/CMSSW_7_4_1/src/B2GTTbar/test/pileup_reweight.root'),
                                   )

###############################################################################
#OUTPUT
 
#process.TFileService = cms.Service("TFileService", fileName = cms.string("histo_DATA.root"))
process.TFileService = cms.Service("TFileService", fileName = cms.string("output/B2Gtree_B2Gtree_ZprimeToTT_M-1000_W-100_JECdown.root") )

process.out = cms.OutputModule("PoolOutputModule",
                               fileName = cms.untracked.string("ZprimeB2G_out.root"),
                               outputCommands = cms.untracked.vstring('drop *'))
Exemple #6
0
    "PoolSource",
    fileNames=cms.untracked.vstring(
        'file:/afs/cern.ch/work/r/rkunnawa/Run2_Analysis/CMSSW_7_5_5_patch4/src/ZplusJet/RootFiles/ZMM-PromptReco-v1-test.root'
    ))

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

process.TFileService = cms.Service(
    "TFileService", fileName=cms.string('testoutput_ZplusJet.root'))

process.ZplusAKPu4PF = cms.EDAnalyzer(
    'ZplusJetAnalyzer',
    src=cms.InputTag("ak4PFJets"),
    JetType=cms.untracked.string('PF'),
    UEAlgo=cms.untracked.string('pp'),
    radius=cms.int(4),
    muons=cms.InputTag("muons"),
    vertices=cms.InputTag("offlinePrimaryVerticesWithBS"),
    #centralitycollection = cms.InputTag("hiCentrality"),
    #centralitybincollection = cms.InputTag("centralityBin","HFtowers"),
    #JetCorrections = cms.string(""),
    #PFcands = cms.InputTag("particleFlowTmp"),
    RThreshold=cms.double(0.3),
    mRecoJetPtThreshold=cms.double(10))

process.ZplusAKVs4PF = process.ZplusAKVs4PF.clone(
    UEAlgo=cms.untracked.string('Pu'))

process.p = cms.Path(process.ZplusAKPu4PF * process.ZplusAKVs4PF)

process.MessageLogger.cerr.FwkReport.reportEvery = 10000