Esempio n. 1
0
#configMgr.histBackupCacheFile =  "data/MyShapeFitExample_template.root" # the data file of your previous fit (= the backup cache file)

# Scaling calculated by outputLumi / inputLumi
configMgr.inputLumi = 0.001  # Luminosity of input TTree after weighting
configMgr.outputLumi = 4.713  # Luminosity required for output histograms
configMgr.setLumiUnits("fb-1")

# Set the files to read from
bgdFiles = []
if configMgr.readFromTree:
    bgdFiles.append("samples/tutorial/SusyFitterTree_OneSoftEle_BG_v3.root")
    bgdFiles.append("samples/tutorial/SusyFitterTree_OneSoftMuo_BG_v3.root")
else:
    bgdFiles = [configMgr.histCacheFile]
    pass
configMgr.setFileList(bgdFiles)

# Dictionnary of cuts for Tree->hist
configMgr.cutsDict[
    "SR"] = "((lep1Pt < 20 && lep2Pt<10 && met>250 && mt>100 && jet1Pt>130 && jet2Pt>25 && AnalysisType==7) || (lep1Pt < 25 && lep2Pt<10 && met>250 && mt>100 && jet1Pt>130 && jet2Pt>25 && AnalysisType==6))"

# Tuples of nominal weights without and with b-jet selection
configMgr.weights = [
    "genWeight", "eventWeight", "leptonWeight", "triggerWeight",
    "truthWptWeight", "bTagWeight2Jet"
]

# QCD weights without and with b-jet selection
# we turn the QCD background of for the tutorial as we do not want to use ATLAS data
#configMgr.weightsQCD = "qcdWeight"
#configMgr.weightsQCDWithB = "qcdBWeight"
configMgr.outputFileName = "results/" + configMgr.analysisName + "_Output.root"

# Scaling calculated by outputLumi / inputLumi
configMgr.inputLumi = 0.001  # Luminosity of input TTree after weighting
configMgr.outputLumi = 4.713  # Luminosity required for output histograms
configMgr.setLumiUnits("fb-1")

# Set the files to read from
bgdFiles = []
if configMgr.readFromTree:
    bgdFiles.append("samples/tutorial/SusyFitterTree_OneSoftEle_BG_v3.root")
    bgdFiles.append("samples/tutorial/SusyFitterTree_OneSoftMuo_BG_v3.root")
else:
    bgdFiles = [configMgr.histCacheFile]
    pass
configMgr.setFileList(bgdFiles)

# Dictionnary of cuts for Tree->hist
configMgr.cutsDict[
    "SR"
] = "((lep1Pt < 20 && lep2Pt<10 && met>250 && mt>100 && jet1Pt>130 && jet2Pt>25 && AnalysisType==7) || (lep1Pt < 25 && lep2Pt<10 && met>250 && mt>100 && jet1Pt>130 && jet2Pt>25 && AnalysisType==6))"


# Tuples of nominal weights without and with b-jet selection
configMgr.weights = ("genWeight", "eventWeight", "leptonWeight", "triggerWeight", "truthWptWeight", "bTagWeight2Jet")

# QCD weights without and with b-jet selection
configMgr.weightsQCD = "qcdWeight"
configMgr.weightsQCDWithB = "qcdBWeight"

# --------------------
Esempio n. 3
0
## Set the luminosity
################################################################################################
userPrint("Specifying Luminosity")

configMgr.inputLumi = userOpts.inputLumi
configMgr.outputLumi = userOpts.outputLumi
configMgr.setLumiUnits(userOpts.lumiUnits)

################################################################################################
## Specify the background input files
################################################################################################
userPrint("Specifying background input")

bkgFiles = []
bkgFiles.append(userOpts.bkgFile)
configMgr.setFileList(bkgFiles)

################################################################################################
## Configure the selection
################################################################################################
userPrint("Configuring selection")

configMgr = selectionConfig(configMgr)

################################################################################################
## Configure weight
################################################################################################
userPrint("Sepcifying Weight")

configMgr.weights = ["eventweight", "1.0"]