Example #1
0
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)
Example #2
0
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)
    #anal.Run(output_dir,conf_ak5_calo,[HT_Run2011A_AllReco_17June]) # 2579 files (and growing)
    #anal.Run(output_dir,conf_ak5_calo,[HT_Run2011A_L1OffSet]) # 2579 files (and growing)
elif ( dataset == 2 ) : # QCD Pythia
    anal.Run(output_dir,conf_ak5_calo,QCD_TuneZ2_7TeV_pythia6_Spring11_PU_START311_ALL) # 2462 files