Skip to content

Instantly share code, notes, and snippets.

View jpata's full-sized avatar

Joosep Pata jpata

View GitHub Profile
@jpata
jpata / eoscp.sh
Created August 10, 2012 08:14
Copy all files in EOS directory
#!/bin/bash
INDIR=/eos/cms/store/caf/user/dkonst/T_t-channel_TuneZ2star_8TeV-powheg-tauola_CAF_EDM/T_t-channel_TuneZ2star_8TeV-powheg-tauola-AOD-FASTSIM_5_2_6/ c02e037c328b91752ebcd32fa1215b6f
EOSCMD="/afs/cern.ch/project/eos/installation/0.1.0-22d/bin/eos.select"
COPYCMD=xrdcp
FILES=$($EOSCMD ls $INDIR)
OUTDIR=/afs/cern.ch/user/j/jpata/work/public
for f in $FILES
do
@jpata
jpata / crab_parse.py
Created August 16, 2012 09:26
parse the output of crab and resubmit jobs
import re
import subprocess
import os
import sys
import datetime
import time
import signal
import glob
import pdb
import datetime
[MULTICRAB]
cfg=crab_ntuples_rew.cfg
[COMMON]
USER.ui_working_dir = /afs/cern.ch/user/j/jpata/private/stpol/CMSSW/CMSSW_4_2_8_patch7/src/TopQuarkAnalysis/SingleTop/test/crab_multicrab2
USER.publish_data_name = systematics_reweight_7to8TeV_v2
[T_t_default]
CMSSW.datasetpath=/T_TuneZ2_t-channel_7TeV-powheg-tauola/Fall11-PU_S6_START42_V14B-v1/AODSIM
@jpata
jpata / crab.cfg
Created August 21, 2012 09:49
crab.cfg for T2_EE_Estonia
[CMSSW]
total_number_of_events=-1
events_per_job=20000
#output_file=STPol_WJets_ntuple.root
#pycfg_params=outputFile=STPol_WJets_ntuple.root
pset=step1_cfg.py
datasetpath=/T_t-channel_TuneZ2star_8TeV-powheg-tauola_CAF_EDM/dkonst-T_t-channel_TuneZ2star_8TeV-powheg-tauola-AOD-FASTSIM_5_2_6_filter-f4dff057d52a2128cc15dd525ba19b60/USER
dbs_url= https://cmsdbsprod.cern.ch:8443/cms_dbs_ph_analysis_02_writer/servlet/DBSServlet
get_edm_output = 1
@jpata
jpata / gist:5495494
Last active December 16, 2015 20:49
int addr = 0
void loop()
{
if (digitalRead(D1) == 1) {
// need to divide by 4 because analog inputs range from
// 0 to 1023 and each byte of the EEPROM can only hold a
// value from 0 to 255.
int val = analogRead(0) / 4;
using namespace reco;
using namespace PhysicsTools;
CombinedMVAJetTagComputerETH::CombinedMVAJetTagComputerETH(
const edm::ParameterSet &params) :
GenericMVAJetTagComputer(params)
{
std::cout << "CombinedMVAJetTagComputerETH::constructor" << std::endl;
inputComputerNames = params.getParameter< std::vector<std::string>>(
using namespace reco;
using namespace PhysicsTools;
CombinedMVAJetTagComputerETH::CombinedMVAJetTagComputerETH(
const edm::ParameterSet &params) :
GenericMVAJetTagComputer(params)
{
std::cout << "CombinedMVAJetTagComputerETH::constructor" << std::endl;
inputComputerNames = params.getParameter< std::vector<std::string>>(
@jpata
jpata / dumpjets.py
Created January 20, 2016 10:31
cMVAv2 CMSSW 763 validation
import ROOT
from DataFormats.FWLite import Events, Handle
events = Events ('exerciseIPartII_histos.root')
handle = Handle ('std::vector<pat::Jet>')
label = ("selectedPatJets")
for event in events:
event.getByLabel(label, handle)
import ROOT
f = ROOT.TFile("/Users/joosep/Documents/tth/data/ntp/v16/Jan18_fixbdt/TT_TuneCUETP8M1_13TeV-powheg-pythia8.root")
t = f.Get("tree")
of = ROOT.TFile("x.root", "RECREATE")
of.cd()
t.CopyTree("is_sl==1")
of.Add(t)
MEM::get_support: E(reco) = 148.35699 ==> range at 0.98 CL is [68.356995, 430.85699] (stepping every 2.5 GeV)
MEM::get_support: E(reco) = 148.35699 ==> range at 0.98 CL is [70.856995, 650.85699] (stepping every 2.5 GeV)
MEM::get_support: E(reco) = 138.43381 ==> range at 0.98 CL is [63.433807, 403.43381] (stepping every 2.5 GeV)
MEM::get_support: E(reco) = 138.43381 ==> range at 0.98 CL is [65.933807, 610.93381] (stepping every 2.5 GeV)
MEM::get_support: E(reco) = 75.970924 ==> range at 0.98 CL is [23.470924, 243.47092] (stepping every 2.5 GeV)
MEM::get_support: E(reco) = 75.970924 ==> range at 0.98 CL is [30.970924, 365.97092] (stepping every 2.5 GeV)
MEM::get_support: E(reco) = 43.183689 ==> range at 0.98 CL is [0, 155.68369] (stepping every 2.5 GeV)
MEM::get_support: E(reco) = 43.183689 ==> range at 0.98 CL is [0, 183.18369] (stepping every 2.5 GeV)
MEM::get_support: E(reco) = 176.26241 ==> range at 0.98 CL is [98.762405, 498.76241] (stepping every 2.5 GeV)
MEM::get_support: E(reco) = 176.26241 ==> range at