brynmathias

import FWCore.ParameterSet.Config as cms

from UserCode.L1TriggerDPG.l1Ntuple_cfg import *

# global tag

process.GlobalTag.globaltag = cms.string('GR_R_42_V14::All')

# load standard RECO for MET cleaning
process.load('Configuration.StandardSequences.Reconstruction_cff')

brynmathias

void L1Analysis::L1AnalysisRecoJet::SetCaloJet(const edm::Event& event,
                 const edm::EventSetup& setup,
                 edm::Handle<reco::CaloJetCollection> caloJets,
                 edm::Handle<edm::ValueMap<reco::JetID> > jetsID,
                 const JetCorrector* caloJetCorrector,
                 unsigned maxJet)
{

  recoJet_.nJets=0;

brynmathias

import FWCore.ParameterSet.Config as cms

l1RecoTreeProducer = cms.EDAnalyzer("L1RecoTreeProducer",

  jetTag                  = cms.untracked.InputTag("ak5CaloJets"),
  jetIdTag                = cms.untracked.InputTag("ak5JetID"),
  metTag                  = cms.untracked.InputTag("met"),
  ebRecHitsTag            = cms.untracked.InputTag("ecalRecHit:EcalRecHitsEB"),
  eeRecHitsTag            = cms.untracked.InputTag("ecalRecHit:EcalRecHitsEE"),
  superClustersBarrelTag  = cms.untracked.InputTag("hybridSuperClusters"),

brynmathias

#include "AlphaT.hh"
#include "CommonOps.hh"
#include "EventData.hh"
#include "GenMatrixBin.hh"
#include "Jet.hh"
#include "KinSuite.hh"
#include "Lepton.hh"
#include "Math/VectorUtil.h"
#include "Photon.hh"
#include "TH1D.h"

brynmathias

//////////////////////////////////////
// Example root macro for l1 ntuples
//////////////////////////////////////

#ifndef L1Ntuple_h
#define L1Ntuple_h

#include <vector>
#include <string>
#include <iostream>

brynmathias

import FWCore.ParameterSet.Config as cms

# make L1 ntuples from RAW+RECO

process = cms.Process("L1NTUPLE")

# import of standard configurations
process.load('Configuration/StandardSequences/Services_cff')
process.load('FWCore/MessageService/MessageLogger_cfi')
#process.load('Configuration/StandardSequences/GeometryIdeal_cff')

brynmathias

  asymErrorCumu = r.TGraphAsymmErrors()
  asymErrorCumu.Divide(cumuNomTotal,cumuDenomTotal)

  for bin in range(CumuNomList[0].GetNbinsX()):
    binLumi = 0
    BinEff = 0.
    point = bin + 1
    errPlus = 0.
    errMinus = 0.
    yval = 0.

brynmathias

#!/bin/bash

pwd=$PWD

# prompt analysis directory
export CAF_TRIGGER=/vols/cms03/bm409/
export L1CMSSW=CMSSW_4_2_5

# 64bit architecture
export SCRAM_ARCH=slc5_amd64_gcc434

brynmathias

/** \class HLTAlphaTFilter
*
*
*  \author Bryn Mathias
*
*/

#include "HLTrigger/JetMET/interface/HLTAlphaTFilter.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/HLTReco/interface/TriggerFilterObjectWithRefs.h"

brynmathias

#ifndef HLTAlphaTFilter_h
#define HLTAlphaTFilter_h

/** \class HLTAlphaTFilter
 *
 *  \author Bryn Mathias
 *
 */

#include "HLTrigger/HLTcore/interface/HLTFilter.h"