#!/usr/bin/env python from sys import argv import argparse parser = argparse.ArgumentParser(description='plot stuff') parser.add_argument('--outdir', metavar='outdir', type=str) parser.add_argument('--label', metavar='label', type=str) args = parser.parse_args() argv = [] import ROOT as root from PandaCore.Tools.Load import Load Load('Tools', 'EventSyncher') Load('Drawers', 'CanvasDrawer') es = root.EventSyncher() if args.label == 'MET': f1 = '/mnt/hadoop/scratch/bmaier/panda/v7/flat/MET.root' f2 = '/mnt/hadoop/scratch/bmaier/panda/v9/flat/MET.root' else: f1 = '/mnt/hadoop/scratch/bmaier/panda/v7/flat/TTbar_CUETP8M2T4.root' f2 = '/mnt/hadoop/scratch/bmaier/panda/v9/flat/TTbar.root' es.preselection = '(( metFilter==1 ) && ( (( nFatjet==1 && fj1Pt>250 && fabs(fj1Eta)<2.4 ) && ( pfUWmag>250 && dphipfUW>0.5 && nLoosePhoton==0 && nTau==0 && nLooseLep==1 && looseLep1IsTight==1 && abs(looseLep1PdgId)==13 )) ))' formulae = [ ('fj1MSD', 0, 300, 'm_{SD}'), ('fj1Pt', 250, 500, 'p_{T}'), ('fj1RawPt', 250, 500, 'raw p_{T}'),
from argparse import ArgumentParser sname = argv[0] parser = ArgumentParser() parser.add_argument('--silent', action='store_true') parser.add_argument('arguments', type=str, nargs='+') args = parser.parse_args() arguments = args.arguments argv = [] import ROOT as root from PandaCore.Tools.process import * from PandaCore.Tools.Misc import * from PandaCore.Tools.Load import Load Load('Normalizer') pds = {} for k, v in processes.iteritems(): if v[1] == 'MC': pds[v[0]] = (k, v[2]) else: pds[v[0]] = (k, -1) VERBOSE = not args.silent user = environ['USER'] system('mkdir -p /tmp/%s/split' % user) # tmp dir system('mkdir -p /tmp/%s/merged' % user) # tmp dir inbase = environ['SUBMIT_OUTDIR']
infile = basedir + '/WJets.root' parser = argparse.ArgumentParser(description='plot stuff') parser.add_argument('--outdir', metavar='outdir', type=str, default=None) parser.add_argument('--cut', metavar='cut', type=str, default=None) args = parser.parse_args() figsdir = args.outdir argv = [] import ROOT as root from PandaCore.Tools.Load import Load from PandaCore.Tools.Misc import * from PandaCore.Tools.root_interface import Selector Load('HistogramDrawer') lumi = 35900 import PandaAnalysis.VBF.PandaSelection as sel base_cut = sel.cuts['signal'] label = 'inc' if args.cut == 'cnc': base_cut = tAND(sel.cnc, sel.cuts['signal']) label = 'cnc' elif args.cut == 'mjj': base_cut = tAND(sel.mjj, sel.cuts['signal']) label = 'mjj' weight = sel.weights['signal'] % lumi plot = root.HistogramDrawer()
#!/usr/bin/env python from PandaCore.Tools.Load import Load from PandaCore.Tools.Misc import PInfo, PError import ROOT as root Load('PandaAnalyzer') def _dump(a): PInfo('PandaAnalysis.Flat.analysis','Summary of analysis %s:'%(a.name)) print a.monojet for k in dir(a): if k[0] == '_': continue if type(getattr(a, k)) != int: continue PInfo('PandaAnalysis.Flat.analysis',' %20s = %s'%(k, 'True' if getattr(a, k) else 'False')) def _analysis(name, verbose, **kwargs): a = root.Analysis(name) for k,v in kwargs.iteritems(): if not hasattr(a, k): PError('PandaAnalysis.Flat.analysis','Could not set property %s'%k) return None setattr(a, k, bool(v)) setattr(a, 'dump', lambda : _dump(a)) if verbose: a.dump() return a
parser.add_argument('--infile', metavar='infile', type=str, default=None) parser.add_argument('--variable', metavar='variable', type=str, default='barrelHTMiss') args = parser.parse_args() figsdir = args.outdir argv = [] import ROOT as root from PandaCore.Tools.Load import Load from PandaCore.Tools.Misc import * from math import sqrt Load('GraphAsymmErrDrawer') root.gROOT.SetBatch() lumi = 36000 plot = root.GraphAsymmErrDrawer() plot.SetTDRStyle() plot.SetLumi(lumi / 1000.) plot.DrawEmpty(True) plot.InitLegend(0.55, 0.76, 0.88, 0.9) #plot.SetAutoRange(False) root.gStyle.SetOptStat(0) fin = root.TFile(args.infile) fout = root.TFile(args.infile.replace('param', 'fit'), 'RECREATE') #ffit = root.TF1("sigmoidFit", "[2]*pow(1+[3]*TMath::Exp(-[0]*(x-[1])),-1./[4])", 10, 1000) #ffit.SetParNames("k", "E_{c}","#varepsilon_{obs}","Q","#nu");
import argparse parser = argparse.ArgumentParser(description='interpolate') parser.add_argument('--mass', type=str) parser.add_argument('--mass_reco', type=str, default=None) args = parser.parse_args() sname = argv[0].split('/')[-1] argv = [] from PandaCore.Tools.Misc import * from PandaCore.Tools.Load import Load from PandaCore.Tools.root_interface import read_tree, draw_hist import ROOT as root from math import sqrt Load('BranchAdder') ba = root.BranchAdder() ba.verbose = False ba.formula = 'genBosonPt' bins = array('f', [250, 280, 310, 350, 400, 450, 600, 1000]) N = len(bins) - 1 h_base = root.TH1F('hbase', '', N, bins) systs = [ 'btagUp', 'btagDown', 'mistagUp', 'mistagDown', 'sjbtagUp', 'sjbtagDown', 'sjmistagUp', 'sjmistagDown' ] basedir = getenv('PANDA_FITTING')
parser = argparse.ArgumentParser(description='plot stuff') parser.add_argument('--outdir',metavar='outdir',type=str) parser.add_argument('--infile',metavar='infile',type=str,default=None) args = parser.parse_args() if not args.infile: args.infile = getenv('PANDA_FITTING') + '/fittingForest.root' argv = [] import ROOT as root from math import sqrt from collections import namedtuple from array import array from PandaCore.Tools.Load import Load from PandaCore.Tools.root_interface import draw_hist, read_tree Load('PandaCoreTools') Load('PandaCoreDrawers') # create some global variables plot = root.HistogramDrawer() plot.SetTDRStyle() plot.InitLegend() plot.SetLumi(35.8) f_input = root.TFile(args.infile) fztow = root.TFile(getenv('PANDA')+'/data/theory_uncs/wtoz_unc.root') fztoa = root.TFile(getenv('PANDA')+'/data/theory_uncs/atoz_unc.root') f_input.cd() recoil_bins = array('f',[250,280,310,350,400,450,600,1000]) hbase = root.TH1D('dummy','',len(recoil_bins)-1,recoil_bins)
#!/usr/bin/env python import ROOT as root import numpy as np import root_numpy as rnp import PandaAnalysis.Monotop.TestSelection as sel from PandaCore.Tools.Misc import * from PandaCore.Tools.Load import Load ''' Load('Drawers','HistogramDrawer') class testclass(): def __init__(self): self.plot = root.HistogramDrawer() members = [x for x in dir(self.plot) if ('__' not in x)] for m in members: if callable(getattr(self.plot,m)): self.__getPlotMethod(m) else: self.__getPlotMember(m) return def __getPlotMethod(self,x): method = getattr(self.plot,x) setattr(self,x, lambda *args : method(*args)) def __getPlotMember(self,x): member = getattr(self.plot,x) setattr(self,x, member) test = testclass() print dir(test) test.SetLumi(10) '''
import argparse from collections import namedtuple import numpy as np parser = argparse.ArgumentParser(description='plot stuff') parser.add_argument('--outdir',metavar='outdir',type=str) parser.add_argument('--sample',metavar='sample',type=str) args = parser.parse_args() argv = [] basedir = getenv('PANDA_FLATDIR') import root_numpy as rnp import ROOT as root from PandaCore.Tools.Load import Load Load('CanvasDrawer') plot = root.CanvasDrawer() plot.SetTDRStyle() c = root.TCanvas() if args.sample=='QCD': weightname = 'ptweight_qcd' fin = root.TFile.Open(basedir+'/QCD_evt8.root') cut = None else: weightname = 'ptweight_binned' fin = root.TFile.Open(basedir+'/ZpTT.root') cut = 'fj1IsMatched==1 && fj1GenSize<1.44' #fin = root.TFile.Open('/home/snarayan/home000/store/panda/v_8024_4_3//Diboson.root') #weightname = 'normalizedWeight'