Python Analysis Beispiele, icf.core.Analysis Python Beispiele

Beispiel #1

0

Datei anzeigen

Datei: McNoSmear.py Projekt: mengzx/CMS-SUSY-2012

vbtfElectronIdFilter = Electron_IDFilter(vbtfelectronidWP95ps.ps())
ra3PhotonIdFilter = Photon_IDFilter(ra3photonidps.ps())


def addCutFlowMC(b):
    #  b.AddWeightFilter("Weight", vertex_reweight)
    b.AddMuonFilter("PreCC", vbtfMuonId_cff)
    b.AddPhotonFilter("PreCC", ra3PhotonIdFilter)
    b.AddElectronFilter("PreCC", vbtfElectronIdFilter)
    b += cutTreeMC


#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC = Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

outDir = "../results_" + strftime("%d_%b_%H") + "//NoSmear/"
ensure_dir(outDir)

from montecarlo.Summer11.QCD_Summer11_madgraph_All import *

from montecarlo.Summer11.TTJets_TuneZ2_7TeV_madgraph_tauola_Summer11_PU_S4_START42_V11_v1_V15_03_19_jetCorrections_L1FastJet_L2Relative_L3Absolute_jetCollections_ak5calo_ak5pf_hbheNoiseFilterDefaultIsoReq_1 import *
anal_ak5_caloMC.Run(outDir, conf_ak5_caloMC, [
    TTJets_TuneZ2_7TeV_madgraph_tauola_Summer11_PU_S4_START42_V11_v1_V15_03_19_jetCorrections_L1FastJet_L2Relative_L3Absolute_jetCollections_ak5calo_ak5pf_hbheNoiseFilterDefaultIsoReq_1
])

Beispiel #2

0

Datei anzeigen

Datei: MuonNoSmear_MC_None_TreeThr86_Zinv.py Projekt: mengzx/CMS-SUSY-2012

    elif PUS == "PUS6" and DataSet == "MC":
        b.AddWeightFilter("Weight", vertex_reweight_PUS6)
    elif PUS != "PUS4" or PUS != "PUS6" or DataSet == "Data":
        print "Vertex reweighting is applied"

    if MuonForm == None:
        b.AddMuonFilter("PreCC", vbtfMuonId_cff)
    else:
        b.AddMuonFilter("PreCC", CustomMuID)

    b.AddPhotonFilter("PreCC", ra3PhotonIdFilter)
    b.AddElectronFilter("PreCC", vbtfElectronIdFilter)
    b += cutTreeMC


#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC = Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

#outDir = "../results_"+strftime("%d_%b")+"//MuonNoSmearTT/"
outDir = "../results_" + strftime("%d_%b") + (
    "//%s%s_%s%s/" % (RunScriptName, SampleIndicated, PUS, ExtraInstruction))
ensure_dir(outDir)

anal_ak5_caloMC.Run(outDir, conf_ak5_caloMC, sample)

Beispiel #3

0

Datei anzeigen

Datei: scaling_cfg.py Projekt: bainbrid/usercode

cut_flow.TAttach(numComPhotons, oddPhoton)
cut_flow.TAttach(oddPhoton, numComElectrons)
cut_flow.TAttach(numComElectrons, oddElectron)
cut_flow.TAttach(oddElectron, numComMuons)
cut_flow.TAttach(numComMuons, oddMuon)
cut_flow.TAttach(oddMuon, test)

# -----------------------------------------------------------------------------
# Dataset

tmp = PSet(
    Name="LM1",
    Events=100000,
    Format=("ICF", 2),
    File=[
        "/Users/bainbrid/Desktop/work/data/LM1.root",
    ],
    CrossSection=4.888,
    #FirstEntry = 1,
    LastEntry=1000,
)

# -----------------------------------------------------------------------------
# Analysis

anal = Analysis("Ratio")
anal.AddPhotonFilter("PreCC", phot_filter)
anal.AddElectronFilter("PreCC", ele_filter)
anal += cut_flow
anal.Run(".", conf_ak5_calo, [tmp])

Beispiel #4

0

Datei anzeigen

Datei: RA1_MuonData_2.py Projekt: DarrenBurton/RA1_Muon_Analysis


def addCutFlowBtagData(a) :
  a.AddPhotonFilter("PreCC",ra3PhotonIdFilter)
  a.AddElectronFilter("PreCC",vbtfElectronIdFilter)
  a.AddMuonFilter("PreCC",CustomMuID)
  a+=cutBtagTree

# AK5 Calo

conf_ak5_caloData = deepcopy(defaultConfig)
conf_ak5_caloData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloData.XCleaning = deepcopy(default_cc)
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData=Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)


conf_ak5_caloBtagData = deepcopy(defaultConfig)
conf_ak5_caloBtagData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloBtagData.XCleaning = deepcopy(default_cc)
conf_ak5_caloBtagData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloBtagData=Analysis("AK5Calo")
addCutFlowBtagData(anal_ak5_caloBtagData)


from data.Run2011.HT_Run2011AB import *
from data.Run2011.MuHad_Run2011AB import *
from HT_L1OffSet import *

Beispiel #5

0

Datei anzeigen

    a.AddWeightFilter("Weight", PreScaleWeights)
    # a.AddJetFilter("PreCC",JetCorrections)
    a.AddPhotonFilter("PreCC", ra3PhotonIdFilter)
    a.AddElectronFilter("PreCC", vbtfElectronIdFilter)
    a.AddMuonFilter("PreCC", vbtfMuonId_cff)
    a += cutTree


# AK5 Calo

conf_ak5_caloData = deepcopy(defaultConfig)
conf_ak5_caloData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloData.XCleaning = deepcopy(default_cc)
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData = Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)

# AK5 PF

conf_ak5_pfData = deepcopy(defaultConfig)
conf_ak5_pfData.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfData.XCleaning = deepcopy(default_cc)
conf_ak5_pfData.Common = deepcopy(default_common)
anal_ak5_pfData = Analysis("AK5PF")
addCutFlowData(anal_ak5_pfData)

tedSkim = PSet(Name="TedSkim",
               Format=("ICF", 3),
               Weight=1.0,
               File="../resultsWithSingleTop//Data/AK5Calo_Jets.root")

Beispiel #6

0

Datei anzeigen

cutTreeMC.TAttach(htCut350, TriggerPlots)

from ra1objectid.vbtfElectronId_cff import *
from ra1objectid.vbtfMuonId_cff import *
from ra1objectid.ra3PhotonId_cff import *

#JetSmear = JetSmear(0.1,30)
vbtfElectronIdFilter = Electron_IDFilter(vbtfelectronidWP95ps.ps())
ra3PhotonIdFilter = Photon_IDFilter(ra3photonidps.ps())


def addCutFlowMC(b):
    b.AddPhotonFilter("PreCC", ra3PhotonIdFilter)
    b.AddElectronFilter("PreCC", vbtfElectronIdFilter)
    b += cutTreeMC


#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC = Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

Triggersamples = [LM0, LM1, LM2, LM3, LM4]

ensure_dir("../results/Triggers/")
anal_ak5_caloMC.Run("../results/Triggers/", conf_ak5_caloMC, Triggersamples)

Beispiel #7

0

Datei anzeigen

Datei: validation_cfg.py Projekt: bainbrid/usercode

data = PSet(
    Name="Data",
    Format=("ICF", 3),
    File=[
        "dcap://gfe02.grid.hep.ph.ic.ac.uk:22128//pnfs/hep.ph.ic.ac.uk/data/cms/store/user/georgia//ICF/automated/2011_04_09_18_07_10///SusyCAF_Tree_100_1_OjA.root",
        #"/vols/cms02/bainbrid/qcd/data/HT_Run2011A_PromptReco_v1.root",
    ],
    Weight=1.,
    #FirstEntry = 1,
    #LastEntry = 10,
)

from montecarlo.LMx import LM0
from montecarlo.QCD_Pythia6_384patch3_V14_00_02_ALL import *
from data.Jet_35pb_WithTP_json051110 import Jet_35pb_WithTP_json051110

# -----------------------------------------------------------------------------
# Analysis

anal = Analysis("Validation")
anal.AddPhotonFilter("PreCC", phot_filter)
anal.AddElectronFilter("PreCC", ele_filter)
#anal.AddJetFilter("PreCC",JetCorrections)
anal += cut_flow

anal.Run(".", conf_ak5_calo, [data])
#anal.Run(".",conf_ak5_calo,[lm1])
#anal.Run(".",conf_ak5_calo,[LM0])
#anal.Run(".",conf_ak5_calo,[Jet_35pb_WithTP_json051110])
#anal.Run(".",conf_ak5_calo,QCD_Pythia6_384patch3_V14_00_02_ALL)

Beispiel #8

0

Datei anzeigen

Datei: McNoSmear.py Projekt: CMSRA1/bryn

#JetSmear = JetSmear(0.1,30)
vbtfElectronIdFilter = Electron_IDFilter( vbtfelectronidWP95ps.ps() )
ra3PhotonIdFilter    = Photon_IDFilter( ra3photonidps.ps() )
def addCutFlowMC(b) :
  b.AddPhotonFilter("PreCC",ra3PhotonIdFilter)
  b.AddElectronFilter("PreCC",vbtfElectronIdFilter)
  b+=cutTreeMC

#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC=Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

# AK5 PF

conf_ak5_pfMC = deepcopy(defaultConfig)
conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfMC.XCleaning = deepcopy(default_cc)
conf_ak5_pfMC.Common = deepcopy(default_common)
anal_ak5_pfMC=Analysis("AK5PF")
addCutFlowMC(anal_ak5_pfMC)

# AK5 JPT

conf_ak5_jptMC = deepcopy(defaultConfig)
conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)

Beispiel #9

0

Datei anzeigen

Datei: myKinCorrZ.py Projekt: taroni/usercode

conf = defaultConfig.copy()
conf.Ntuple.Jets.Prefix = "ak5JetPF"
conf.Common.Jets.EtCut = 20.0
conf.Common.Jets.EtaCut = 5.0
conf.Common.Muons.EtaCut = 2.4
conf.Common.Muons.PtCut = 15.0
conf.Common.Muons.TrkIsoCut = -1.
conf.Common.Muons.CombIsoCut = 0.1
conf.Common.Electrons.PtCut = 15.0
conf.Common.Electrons.EtaCut = 2.4
conf.Common.Photons.EtCut = 30.
conf.Common.Electrons.ApplyID = False

# Create the analysis
a = Analysis("KinCorrection")

tree = Tree("Main")
Trigger = PSet(McAnal=False, MSugraScan=0., TriggerBits=PSet(bit1=""))
ssdlTrigg = SSDLTrigger("Trigger", Trigger.ps())

ZeroMuons = OP_NumComMuons("==", 0)

KinCorrParPlus = PSet(
    zMassMin=71.,
    zMassMax=111.,
    BarrelEtaMax=1.4442,
    EndCapEtaMin=1.56,
    EndCapEtaMax=2.5,
    MinElecRescaledPt=0.,
    c_ErNu_pt=0.,

Beispiel #10

0

Datei anzeigen

  b+=cutTreeMC


def addCutFlowData(a):
  a.AddPhotonFilter("PreCC",ra3PhotonIdFilter)
  a.AddElectronFilter("PreCC",vbtfElectronIdFilter)
  a+=cutTreeData


#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
# conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC=Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

#AK5 Calo
conf_ak5_caloData = deepcopy(defaultConfig)
conf_ak5_caloData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloData.XCleaning = deepcopy(default_cc)
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_caloMC.Common.print_out()
anal_ak5_caloData=Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)



testFile = PSet(
Name = "TestFile",

Beispiel #11

0

Datei anzeigen

Datei: McJESetaMinus.py Projekt: CMSRA1/bryn

# JESUncert(-0.02,true)   means -2% per unit eta
JESUncert = JESUncert(-0.02, True)


def addCutFlowMC(b):
    b.AddJetFilter("PreCC", JESUncert)
    b += cutTreeMC


#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC = Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

# AK5 PF

conf_ak5_pfMC = deepcopy(defaultConfig)
conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfMC.XCleaning = deepcopy(default_cc)
conf_ak5_pfMC.Common = deepcopy(default_common)
anal_ak5_pfMC = Analysis("AK5PF")
addCutFlowMC(anal_ak5_pfMC)

# AK5 JPT

conf_ak5_jptMC = deepcopy(defaultConfig)
conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)

Beispiel #12

0

Datei anzeigen

Datei: McJESminus.py Projekt: CMSRA1/bryn


def addCutFlowMC(b):
    b.AddPhotonFilter("PreCC", ra3PhotonIdFilter)
    b.AddElectronFilter("PreCC", vbtfElectronIdFilter)
    b.AddJetFilter("PreCC", JESUncert)
    b += cutTreeMC


#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC = Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

# AK5 PF

conf_ak5_pfMC = deepcopy(defaultConfig)
conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfMC.XCleaning = deepcopy(default_cc)
conf_ak5_pfMC.Common = deepcopy(default_common)
anal_ak5_pfMC = Analysis("AK5PF")
addCutFlowMC(anal_ak5_pfMC)

# AK5 JPT

conf_ak5_jptMC = deepcopy(defaultConfig)
conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)

Beispiel #13

0

Datei anzeigen

    a += kinsuitecomplot
    a += myplots
    a += HadAlphaT


# -----------------------------------------------------------------------------
# Analyses

# AK5 Calo

conf_ak5_calo = deepcopy(defaultConfig)
conf_ak5_calo.Ntuple = deepcopy(ak5_calo)
conf_ak5_calo.XCleaning = deepcopy(default_cc)
conf_ak5_calo.Common = deepcopy(default_common)

anal_ak5_calo = Analysis("AK5Calo")
addCutFlow(anal_ak5_calo)

# AK5 JPT

conf_ak5_jpt = deepcopy(defaultConfig)
conf_ak5_jpt.Ntuple = deepcopy(ak5_jpt)
conf_ak5_jpt.XCleaning = deepcopy(default_cc)
conf_ak5_jpt.Common = deepcopy(default_common)

anal_ak5_jpt = Analysis("AK5JPT")
addCutFlow(anal_ak5_jpt)

# Run analyses

#anal_ak5_calo.Run("results",conf_ak5_calo,[lm0])

Beispiel #14

0

Datei anzeigen

Datei: turnon.py Projekt: bainbrid/usercode

lm1 = PSet(
    Name="lm1",
    Events=100000,
    Format=("ICF", 2),
    File=[
        #"dcap://gfe02.grid.hep.ph.ic.ac.uk:22128//pnfs/hep.ph.ic.ac.uk/data/cms/store/user/bainbrid/ICF/automated/2010_07_12_17_52_54/LM1.Spring10-START3X_V26_S09-v1.GEN-SIM-RECO/SusyCAF_Tree_1_1.root",
        "/Users/bainbrid/Desktop/work/data/LM1.root",
    ],
    CrossSection=4.888,
    #FirstEntry = 1,
    #LastEntry = 10,
)

from data.Run2011.HT_Run2011A_PromptReco_v1 import HT_Run2011A_PromptReco_v1
from data.Run2011.HT_Run2011_promptReco_DCS import HT_Run2011_promptReco_DCS

# -----------------------------------------------------------------------------
# Analysis

anal = Analysis("TurnOn")
anal.AddMuonFilter("PreCC", vbtf_muon)
anal.AddElectronFilter("PreCC", vbtf_electron)
anal.AddPhotonFilter("PreCC", ra3_photon)
anal += cut_flow

#anal.Run(".",conf,[data])
#anal.Run(".",conf,[tmp])
#anal.Run(".",conf,[lm1])
#anal.Run(".",conf_ak5_calo,[HT_Run2011A_PromptReco_v1])
anal.Run(".", conf_ak5_calo, [HT_Run2011_promptReco_DCS])

Beispiel #15

0

Datei anzeigen

Datei: original_cfg_v81.py Projekt: bainbrid/usercode

from montecarlo.Spring11.TToBLNu_TuneZ2_tW_channel_7TeV_madgraph_Spring11_PU_S1_START311_V1G1_v1 import *
from montecarlo.Spring11.TToBLNu_TuneZ2_s_channel_7TeV_madgraph_Spring11_PU_S1_START311_V1G1_v1 import *
from montecarlo.Spring11.LMx_SUSY_sftsht_7TeV_pythia6_Spring11_PU_S1_START311_V1G1_v1 import *

from data.Jet_35pb_WithTP_json051110 import Jet_35pb_WithTP_json051110
from data.Run2011.HT42_incomplete import HT42_incomplete
from data.Run2011.HT_Run2011_promptReco_DCS import HT_Run2011_promptReco_DCS
from data.Run2011.HT_Run2011A import HT_Run2011A
from data.Run2011.HT_Run2011A_PromptReco_v1 import HT_Run2011A_PromptReco_v1
from data.Run2011.HT_Run2011A_AllReco_17June import HT_Run2011A_AllReco_17June
from data.Run2011.HT_Run2011A_L1OffSet import HT_Run2011A_L1OffSet

# -----------------------------------------------------------------------------
# Analysis

anal=Analysis("Ratio")
if ( not DataOption ) :
    anal.AddWeightFilter("Weight", vertex_reweight)
anal.AddMuonFilter("PreCC",vbtf_muon)
anal.AddElectronFilter("PreCC",vbtf_electron)
anal.AddPhotonFilter("PreCC",ra3_photon)
#anal.AddJetFilter("PreCC",JetCorrections)
 
anal+=cut_flow

if   ( dataset == 0 ) : # Data 2010
    anal.Run(output_dir,conf_ak5_calo,[Jet_35pb_WithTP_json051110]) # 160 files
elif ( dataset == 1 ) : # Data 2011
    #anal.Run(output_dir,conf_ak5_calo,[HT_Run2011_promptReco_DCS]) # 1543 files (and growing)
    #anal.Run(output_dir,conf_ak5_calo,[HT42_incomplete]) # 2579 files (and growing)
    anal.Run(output_dir,conf_ak5_calo,[HT_Run2011A]) # 2579 files (and growing)

Beispiel #16

0

Datei anzeigen

  a+=plots_hadronic
  a+=plots_comjet
  a+=plots_objkin
  a+=plots_kinsuite
  # -----------------------------------------------------------------------------
# Definition of analyses
# Analyses

# IC5 Calo

conf_ic5_calo = deepcopy(defaultConfig)
conf_ic5_calo.Ntuple = deepcopy(ic5_calo)
conf_ic5_calo.XCleaning = deepcopy(default_cc)
conf_ic5_calo.Common = deepcopy(default_common)
# conf_ic5_calo.Common.print_out()
anal_ic5_calo=Analysis("IC5Calo")
addCutFlow(anal_ic5_calo)

# AK5 Calo

conf_ak5_calo = deepcopy(defaultConfig)
conf_ak5_calo.Ntuple = deepcopy(ak5_calo)
conf_ak5_calo.XCleaning = deepcopy(default_cc)
conf_ak5_calo.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_calo=Analysis("AK5Calo")
addCutFlow(anal_ak5_calo)

# AK5 JPT

conf_ak5_jpt = deepcopy(defaultConfig)

Beispiel #17

0

Datei anzeigen

Datei: McJESminus.py Projekt: bainbrid/usercode

# JESUncert(0.1,false)    means +10%, independent of eta
# JESUncert(-0.1,false)   means -10%, independent of eta
# JESUncert(0.02,true)    means +2% per unit eta
# JESUncert(-0.02,true)   means -2% per unit eta
JESUncert = JESUncert(-0.1,False)
def addCutFlowMC(b) :
  b.AddJetFilter("PreCC",JESUncert)
  b+=cutTreeMC

#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC=Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

# AK5 PF

conf_ak5_pfMC = deepcopy(defaultConfig)
conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfMC.XCleaning = deepcopy(default_cc)
conf_ak5_pfMC.Common = deepcopy(default_common)
anal_ak5_pfMC=Analysis("AK5PF")
addCutFlowMC(anal_ak5_pfMC)

# AK5 JPT

conf_ak5_jptMC = deepcopy(defaultConfig)
conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)

Beispiel #18

0

Datei anzeigen

Datei: trigger_cfg.py Projekt: bainbrid/usercode

        #"/Users/bainbrid/Desktop/work/data/HT_Run2011A_PromptReco_v1.root",
    ],
    Weight=1.,
    LastEntry=10,
)

from data.Run2011.HT_Run2011A_PromptReco_v1 import HT_Run2011A_PromptReco_v1
from data.Run2011.HT_Run2011_promptReco_DCS import HT_Run2011_promptReco_DCS
from data.Run2011.HT42_incomplete import HT42_incomplete
from data.Run2011.HT_Run2011A_L1OffSet import HT_Run2011A_L1OffSet
from data.Run2011.HT_Run2011A import HT_Run2011A
from data.Run2011.HT_Run2011A_AllReco_17June import HT_Run2011A_AllReco_17June

# -----------------------------------------------------------------------------
# Analysis

anal = Analysis("Trigger")
#anal.AddPhotonFilter("PreCC",phot_filter)
#anal.AddElectronFilter("PreCC",ele_filter)
#anal.AddJetFilter("PreCC",JetCorrections)
anal += cut_flow

#anal.Run(".",conf,[data])
#anal.Run(".",conf,[Jet_35pb_WithTP_json051110])
#anal.Run(".",conf_ak5_calo,[HT_Run2011A_PromptReco_v1])
#anal.Run(".",conf_ak5_calo,[HT_Run2011_promptReco_DCS])
#anal.Run(output_dir,conf_ak5_calo,[HT42_incomplete])
#anal.Run(output_dir,conf_ak5_calo,[HT_Run2011A_L1OffSet])
anal.Run(output_dir, conf_ak5_calo, [HT_Run2011A])
#anal.Run(output_dir,conf_ak5_calo,[HT_Run2011A_AllReco_17June])

Beispiel #19

0

Datei anzeigen

def addCutFlowMC(b):
    b.AddJetFilter("PreCC", JetSmear)
    b += cutTreeMC


# -----------------------------------------------------------------------------
# Definition of analyses
# Analyses

#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC = Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

# AK5 Calo

conf_ak5_caloData = deepcopy(defaultConfig)
conf_ak5_caloData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloData.XCleaning = deepcopy(default_cc)
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData = Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)

# AK5 PF

conf_ak5_pfData = deepcopy(defaultConfig)

Beispiel #20

0

Datei anzeigen

    a += cutTreeData

    # -----------------------------------------------------------------------------


# Definition of analyses
# Analyses

# IC5 Calo

conf_ic5_calo = deepcopy(defaultConfig)
conf_ic5_calo.Ntuple = deepcopy(ic5_calo)
conf_ic5_calo.XCleaning = deepcopy(default_cc)
conf_ic5_calo.Common = deepcopy(default_common)
# conf_ic5_calo.Common.print_out()
anal_ic5_calo = Analysis("IC5Calo")
addCutDataFlow(anal_ic5_calo)

# AK5 Calo

# conf_ak5_caloMC = deepcopy(defaultCon# fig)
# conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
# conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
# conf_ak5_caloMC.Common = deepcopy(default_common)
# # conf_ak5_calo.Common.print_out()
# anal_ak5_caloMC=Analysis("AK5Calo")
# addCutMCFlow(anal_ak5_caloMC)
#
# # AK5 Calo
#
# conf_ak5_caloData = deepcopy(defaultConfig)

Beispiel #21

0

Datei anzeigen

Datei: RA1_MuonNoSmear43.py Projekt: DarrenBurton/RA1_Muon_Analysis

cutTreeMC, junkVar, l = MakeMCTree(100. * (325. / 375.), Muon=True)

vbtfElectronIdFilter = Electron_IDFilter(vbtfelectronidWP95ps.ps())
ra3PhotonIdFilter = Photon_IDFilter(ra3photonidps.ps())
CustomMuID = OL_CustomVBTFMuID(mu_id_lower.ps())


def addCutFlowMC(b):
    b.AddWeightFilter("Weight", vertex_reweight)
    b.AddMuonFilter("PreCC", CustomMuID)
    b.AddPhotonFilter("PreCC", ra3PhotonIdFilter)
    b.AddElectronFilter("PreCC", vbtfElectronIdFilter)
    b += cutTreeMC


#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC = Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

outDir = "../results_" + strftime("%d_%b") + "//MuonNoSmear43/"
ensure_dir(outDir)

anal_ak5_caloMC.Run(outDir, conf_ak5_caloMC, L1OffSet_MC_Lower_Bins)

#anal_ak5_caloMC.Run(outDir,conf_ak5_caloMC,Summer11_MC_Lower_Bins)

Beispiel #22

0

Datei anzeigen

Datei: TriggerTurnOns.py Projekt: bainbrid/usercode

    # a.AddWeightFilter("Weight",PreScaleWeights)
    # a.AddJetFilter("PreCC",JetCorrections)
    a.AddPhotonFilter("PreCC", ra3PhotonIdFilter)
    a.AddElectronFilter("PreCC", vbtfElectronIdFilter)
    a.AddMuonFilter("PreCC", vbtfMuonId_cff)
    a += cutTreeData


# AK5 Calo

conf_ak5_caloData = deepcopy(defaultConfig)
conf_ak5_caloData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloData.XCleaning = deepcopy(default_cc)
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData = Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)

# AK5 PF

conf_ak5_pfData = deepcopy(defaultConfig)
conf_ak5_pfData.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfData.XCleaning = deepcopy(default_cc)
conf_ak5_pfData.Common = deepcopy(default_common)
anal_ak5_pfData = Analysis("AK5PF")
addCutFlowData(anal_ak5_pfData)

from data.Run2011.HT_Run2011_promptReco_DCS import *
outDir = "../results/TriggerTurnOns/"
ensure_dir(outDir)
anal_ak5_caloData.Run(outDir, conf_ak5_caloData, [HT_Run2011_promptReco_DCS])

Beispiel #23

0

Datei anzeigen

def addCutFlowData(a):
    a.AddPhotonFilter("PreCC", ra3PhotonIdFilter)
    a.AddElectronFilter("PreCC", vbtfElectronIdFilter)
    a.AddMuonFilter("PreCC", vbtfMuonId_cff)
    a += cutTree


# AK5 Calo

conf_ak5_caloData = deepcopy(defaultConfig)
conf_ak5_caloData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloData.XCleaning = deepcopy(default_cc)
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData = Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)

from HT_L1OffSet import *
from data.Run2011.HT_Run2011AB import *
#from RA2b_Skims import *

#outDir = "../Split_Jsons_"+strftime("%d_%b")+"/Data43/"
#outDir = "../RA2B_"+strftime("%d_%b")+"/Data43/"
outDir = "../results_" + strftime("%d_%b") + "//Data43"
ensure_dir(outDir)

#anal_ak5_caloData.Run(outDir,conf_ak5_caloData,[RA2B_Skim])#HT_Run2011A])
#anal_ak5_caloData.Run(outDir,conf_ak5_caloData,[HT_Run2011AB])
anal_ak5_caloData.Run(outDir, conf_ak5_caloData, [HT_L1OffSet])  #HT_Run2011A])

Beispiel #24

0

Datei anzeigen

Datei: MakeStandardPlots.py Projekt: CMSRA1/bryn

def addCutFlowMC(b) :
  b+=cutTreeMC


# -----------------------------------------------------------------------------
# Definition of analyses
# Analyses

# AK5 Calo

conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloMC=Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

# AK5 Calo

conf_ak5_caloData = deepcopy(defaultConfig)
conf_ak5_caloData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloData.XCleaning = deepcopy(default_cc)
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData=Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)

# -----------------------------------------------------------------------------
# Run analyses
#IMPERIAL COLLEGE

Beispiel #25

0

Datei anzeigen

    a.AddWeightFilter("Weight", PreScaleWeights)
    # a.AddJetFilter("PreCC",JetCorrections)
    a.AddPhotonFilter("PreCC", ra3PhotonIdFilter)
    a.AddElectronFilter("PreCC", vbtfElectronIdFilter)
    a.AddMuonFilter("PreCC", vbtfMuonId_cff)
    a += cutTree


# AK5 Calo

conf_ak5_caloData = deepcopy(defaultConfig)
conf_ak5_caloData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloData.XCleaning = deepcopy(default_cc)
conf_ak5_caloData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloData = Analysis("AK5Calo")
addCutFlowData(anal_ak5_caloData)

# AK5 PF

conf_ak5_pfData = deepcopy(defaultConfig)
conf_ak5_pfData.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfData.XCleaning = deepcopy(default_cc)
conf_ak5_pfData.Common = deepcopy(default_common)
anal_ak5_pfData = Analysis("AK5PF")
addCutFlowData(anal_ak5_pfData)

tedSkim = PSet(Name="TedSkim",
               Format=("ICF", 3),
               Weight=1.0,
               File="../results_" + strftime("%d_%b_%H") +

Beispiel #26

0

Datei anzeigen

Datei: RA1_MuonNoSmear.py Projekt: DarrenBurton/RA1_Muon_Analysis

  b+=cutTreeMC

def addCutFlowBtagMC(a) :
  a.AddWeightFilter("Weight",vertex_reweight)
  a.AddPhotonFilter("PreCC",ra3PhotonIdFilter)
  a.AddElectronFilter("PreCC",vbtfElectronIdFilter)
  a.AddMuonFilter("PreCC",CustomMuID)
  a+=cutTreeBtagMC

#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
#conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC=Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

conf_ak5_caloBtagData = deepcopy(defaultConfig)
conf_ak5_caloBtagData.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloBtagData.XCleaning = deepcopy(default_cc)
conf_ak5_caloBtagData.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_caloBtagData=Analysis("AK5Calo")
addCutFlowBtagMC(anal_ak5_caloBtagData)


outDir = "../results_"+strftime("%d_%b")+"//MuonNoSmear_Btag/"

#outDir = "../results_"+strftime("%d_%b")+"//MuonNoSmear/"
ensure_dir(outDir)

Beispiel #27

0

Datei anzeigen

def addCutFlowMC(b) :
  b.AddJetFilter("PreCC",JetAdd)
#  b.AddWeightFilter("Weight", vertex_reweight)
  b.AddMuonFilter("PreCC",vbtfMuonId_cff)
  b.AddPhotonFilter("PreCC",ra3PhotonIdFilter)
  b.AddElectronFilter("PreCC",vbtfElectronIdFilter)
  b.AddJetFilter("PreCC",JESUncert)
  b+=cutTreeMC

#AK5 Calo
conf_ak5_caloMC = deepcopy(defaultConfig)
conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo)
conf_ak5_caloMC.XCleaning = deepcopy(default_cc)
conf_ak5_caloMC.Common = deepcopy(default_common)
conf_ak5_caloMC.Common.print_out()
anal_ak5_caloMC=Analysis("AK5Calo")
addCutFlowMC(anal_ak5_caloMC)

# AK5 PF

conf_ak5_pfMC = deepcopy(defaultConfig)
conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf)
conf_ak5_pfMC.XCleaning = deepcopy(default_cc)
conf_ak5_pfMC.Common = deepcopy(default_common)
anal_ak5_pfMC=Analysis("AK5PF")
addCutFlowMC(anal_ak5_pfMC)

# AK5 JPT

conf_ak5_jptMC = deepcopy(defaultConfig)
conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)

Beispiel #28

0

Datei anzeigen

Datei: met_cfg.py Projekt: CMSRA1/bryn

  a+=greaterThreeJet
  a+=PlotsMetG3
  a+=plots_comjet
  a+=plots_kinsuite
  # -----------------------------------------------------------------------------
# Definition of analyses
# Analyses

# IC5 Calo

conf_ic5_calo = deepcopy(defaultConfig)
conf_ic5_calo.Ntuple = deepcopy(ic5_calo)
conf_ic5_calo.XCleaning = deepcopy(default_cc)
conf_ic5_calo.Common = deepcopy(default_common)
# conf_ic5_calo.Common.print_out()
anal_ic5_calo=Analysis("IC5Calo")
addCutFlow(anal_ic5_calo)

# AK5 Calo

conf_ak5_calo = deepcopy(defaultConfig)
conf_ak5_calo.Ntuple = deepcopy(ak5_calo)
conf_ak5_calo.XCleaning = deepcopy(default_cc)
conf_ak5_calo.Common = deepcopy(default_common)
# conf_ak5_calo.Common.print_out()
anal_ak5_calo=Analysis("AK5Calo")
addCutFlow(anal_ak5_calo)

# AK5 JPT

conf_ak5_jpt = deepcopy(defaultConfig)