ram1123

import ROOT

# Set ROOT to run in batch mode 
ROOT.gROOT.SetBatch(True)

def print_bins_noSort(hist):
    xaxis = hist.GetXaxis()
    for i in range(1, hist.GetNbinsX() + 1):
        bin_content = hist.GetBinContent(i)
        bin_label = xaxis.GetBinLabel(i)

ram1123

#!/usr/bin/python

print("Converting all of the .dav files in this current directory into .mp4 files using ffmpeg")
import os
from subprocess import call

files = [f for f in os.listdir('.') if os.path.isfile(f)]
for f in files:
    ext = f.split(".")[-1]
    if ext == "dav" or ext == "DAV":

ram1123

import ROOT as ROOT
import sys
ROOT.gROOT.SetBatch(ROOT.kTRUE)

def get_standard_deviation(histogram):
    # Compute the standard deviation of a histogram
    mean = histogram.GetMean()
    variance = histogram.GetStdDev() ** 2
    return variance ** 0.5

ram1123

import os

# Specify the directory containing the PDF files
pdf_dir = '/afs/cern.ch/user/r/rasharma/work/LearnCombine/CMSSW_11_3_4/src/2l2Q_limitSettingTool/ImpactPlotCollection/'

# Get a list of the PDF files in the directory
pdf_files = [f for f in os.listdir(pdf_dir) if f.endswith('.pdf')]

# Create a LaTeX file for the presentation
with open('presentation.tex', 'w') as f:

ram1123

import ROOT

def SignificanceScan(var="evalDNN_HH"):
    sigFile = "/eos/user/r/rasharma/post_doc_ihep/double-higgs/ntuples/January_2021_Production/DNN_Evaluate_HHWWyyDNN_NewModel6_E700_LR10em4_B150_VarMixV1MaxBRmvCorr_Trial1_BalanceYields/GluGluToHHTo2G4Q_node_cHHH1_2017.root"
    # bkgFiles = ["/eos/user/r/rasharma/post_doc_ihep/double-higgs/ntuples/January_2021_Production/DNN_Evaluate_HHWWyyDNN_NewModel6_E700_LR10em4_B150_VarMixV1MaxBRmvCorr_Trial1_BalanceYields/Data_2017.root"]
    bkgFiles = ["/eos/user/r/rasharma/post_doc_ihep/double-higgs/ntuples/January_2021_Production/DNN_Evaluate_HHWWyyDNN_NewModel6_E700_LR10em4_B150_VarMixV1MaxBRmvCorr_Trial1_BalanceYields/Data_2017.root", "/eos/user/r/rasharma/post_doc_ihep/double-higgs/ntuples/January_2021_Production/DNN_Evaluate_HHWWyyDNN_NewModel6_E700_LR10em4_B150_VarMixV1MaxBRmvCorr_Trial1_BalanceYields/Data_2018.root"]

    fileCount = 0

    # Open the input ROOT file

ram1123

import os
import sys


if sys.version_info[0] < 3:
    print("This script requires Python 3. Please use Python 3 to run this script.")
    sys.exit()

# Run crab status command and append in log file
ifRunCrabStatus = False

ram1123

import sys

if sys.version_info[0] < 3:
    print("This script requires Python 3. Please use Python 3 to run this script.")
    sys.exit()

ram1123

request_cpus    = 1
request_memory  = 4GB
request_disk    = 4GB
executable = /afs/cern.ch/user/r/rasharma/work/slc6_condor/generate_events.sh
output     = /afs/cern.ch/user/r/rasharma/work/slc6_condor/$(ClusterId).out
error      = /afs/cern.ch/user/r/rasharma/work/slc6_condor/$(ClusterId).out 
getenv = FALSE
should_transfer_files = NO 
+RequestRuntime  = 3600
+JobFlavour = "workday"

ram1123

import os

sh_file_template = '''#!/bin/bash

echo "Starting job on " `date`
echo "Running on: `uname -a`"
echo "System software: `cat /etc/redhat-release`"
source /cvmfs/cms.cern.ch/cmsset_default.sh
echo "###################################################"
echo "#    List of Input Arguments: "

ram1123

from ROOT import *

f_xBF = TFile('xs_125.0_2017/hzz4l_4muS_13TeV_xs_SM_125_mass4l_v3.Databin0.root','READ')
w_xBF = f_xBF.Get("w")
# w_xBF.Print()

f_LLR = TFile('xs_125.0_2017/hzz4l_4muS_13TeV_xs_SM_125_mass4l_v3.Databin0.root','READ')
w_LLR = f_LLR.Get("w")
# w_LLR.Print()