"/vols/cms02/gouskos/QCD_Pythia_Pt80_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt170_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt300_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt470_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt800_Jun2010/SusyCAF_Tree_*.root"
),
Weights = PSet(
CrossSection = [8.762e+08,6.041e+07,9.238e+05,2.547e+04,1.256e+03 ,
8.798e+01,2.186,0.01122 ],
Events = [6095857,5069664,2065792,3171950,2976108,2159497,2181700,1185024],
PtBin = [15.,30.,80.,170.,300., 470.,800.,1400.],
)
)
utils.reweight([LM0,LM1, LM01, LM01p, LMParis, QCD_AllPtBins_7TeV_Pythia,
zjets_madgraph, Wj,ttbarTauola
],10.)
#Operations
#trigger = OP_TriggerCut("HLT_Jet70U")
#trigger = OP_TwoTriggerCut("HLT_HT100U","HLT_HT120U")
#trigger = OP_TriggerCut("HLT_Mu9")
selection = OP_GoodEventSelection()
HBHE = OP_HadronicHBHEnoiseFilter()
MuonTreeReq=OP_NumComMuons(">=",1)
LeadingJetEta = OP_FirstJetEta(2.4)
unCorLeadJetCut = OP_UnCorLeadJetCut(70.)
PrintOut = PrintAlphaTEvents()
OneMuon = OP_NumComMuons("==",1)
MaxPtCut = OP_UpperMuPtCut(25.)
Weights = PSet(
CrossSection = [8.762e+08,6.041e+07,9.238e+05,2.547e+04,1.256e+03 ,
8.798e+01,2.186,0.01122 ],
Events = [6095857,5069664,2065792,3171950,2976108,2159497,2181700,1185024],
PtBin = [15.,30.,80.,170.,300., 470.,800.,1400.],
)
)
QCD170=PSet(
Name="QCD170",
Format=("ICF",2),
File=["/vols/cms02/gouskos/QCD_Pythia_Pt170_Jun2010/SusyCAF_Tree_*.root"],
CrossSection=25470.
)
utils.reweight([LM0,LM1, LM01, LM01p, LMParis, QCD_AllPtBins_7TeV_Pythia,
zjets_madgraph, Wj,ttbarTauola,em_20to30,em_30to80,em_80to170,bce_20to30,bce_30to80,bce_80to170,QCD170
],10.)
#Operations
trigger = OP_TriggerCut("HLT_Jet50")
#trigger = OP_TriggerCut("HLT_Mu9")
selection = OP_GoodEventSelection()
HBHE = OP_HadronicHBHEnoiseFilter()
EleTreeReq=OP_NumComElectrons(">=",1)
LeadingJetEta = OP_FirstJetEta(2.4)
unCorLeadJetCut = OP_UnCorLeadJetCut(70.)
PrintOut = PrintAlphaTEvents()
OneMuon = OP_NumComMuons("==",1)
MaxPtCut = OP_UpperElePtCut(25.)
ZeroElectrons = OP_NumComElectrons("<=",0)
samples=[
# icf.qcd_pt80,
# icf.QCD_EMenriched_20to30,
icf.qcd_pt80,
# icf.QCD_EMenriched_20to30,
# icf.QCD_EMenriched_30to80,
# icf.QCD_EMenriched_80to170,
# icf.qcd_bctoe_20to30,
# icf.QCD_BCtoE_30to80,
# icf.QCD_BCtoE_80to170,
# icf.qcd_em_30to80,
# icf.qcd_em_80to170,
# icf.Wjets,
# icf.Zjets,
# icf.gammajets_pt30
]
icf.Wjets.LastEntry=2000000
# NLO CrossSection!!!
icf.Wjets.CrossSection = 24170.
#local_ntuple(samples, my_ntuple_path, True)
#utils.switch_ntuple(samples, "/home/alex/Research/Samples/NtupleMap.txt", True)
utils.reweight(samples, lumi)
sig.Run("../results", conf, samples)
#control.Run("../results",conf,samples)
T1TauNu = [SimpModel_t1TauNu]
SMS_test = PSet(
Name = 'test',
Format = ('ICF',3),
File = [
"root://xrootd.grid.hep.ph.ic.ac.uk//store/user/arlogb/SimplifiedModel_T1TauNu_Second/0.5_100.0_50.0to325.0_50.0.lhe_1.root"
],
Weight = 1.0,
# LastEntry = 10
)
SMS_test2 = PSet(
Name = 'test2',
Format = ('ICF',3),
File = [
'root://xrootd.grid.hep.ph.ic.ac.uk//store/user/arlogb/ICF/automated/2011_08_23_16_46_43/SusyCAF_Tree_79_3_rar.root'
],
Weight = 1.0
)
test = [SMS_test]
test2 = [SMS_test2]
reweight(TanB10,976.)
if mode == 'mSugra':
a.Run(outpath, conf, TanB10)
elif mode == 'SMS':
a.Run(outpath, conf, T1TauNu)
"/vols/cms02/gouskos/QCD_Pythia_Pt300_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt470_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt800_Jun2010/SusyCAF_Tree_*.root"
),
Weights = PSet(
CrossSection = [8.762e+08,6.041e+07,9.238e+05,2.547e+04,1.256e+03 ,
8.798e+01,2.186,0.01122 ],
Events = [6095857,5069664,2065792,3171950,2976108,2159497,2181700,1185024],
PtBin = [15.,30.,80.,170.,300., 470.,800.,1400.]
)
)
utils.reweight([LM0,LM1, LM01, LM01p, LMParis, QCD_AllPtBins_7TeV_Pythia,
zjets_madgraph, Wj,ttbarTauola,QCD_Mu_Pt80to120_Summer10_START36_V9_S09_v1,QCD_Mu_Pt50to80_Summer10_START36_V9_S09_v1,QCD_Mu_Pt30to50_Summer10_START36_V9_S09_v1,QCD_Mu_Pt20to30_Summer10_START36_V9_S09_v1,QCD_Mu_Pt170toInf_Summer10_START36_V9_S09_v1,QCD_Mu_Pt120to170_Summer10_START36_V9_S09_v1
],10.)
#Operations
trigger = OP_TriggerCut("HLT_Jet50")
#trigger = OP_TriggerCut("HLT_Mu9")
selection = OP_GoodEventSelection()
HBHE = OP_HadronicHBHEnoiseFilter()
MuonTreeReq=OP_NumComMuons(">=",1)
LeadingJetEta = OP_FirstJetEta(2.4)
unCorLeadJetCut = OP_UnCorLeadJetCut(70.)
PrintOut = PrintAlphaTEvents()
OneMuon = OP_NumComMuons("==",1)
MaxPtCut = OP_UpperMuPtCut(25.)
ZeroElectrons = OP_NumComElectrons("<=",0)
#tree.TAttach(HBHE,OneMuon)
#tree.TAttach(OneMuon, ZeroElectrons)
#tree.TAttach(ZeroElectrons,BasicPlotsPreMHT)
#Run on MC
tree.Attach(OneMuon)
tree.TAttach(OneMuon, ZeroElectrons)
tree.TAttach(ZeroElectrons,BasicPlotsPreMHT)
#sig=Analysis("muAlphaT_TMOneStLoose")
sig=Analysis("muAlphaT")
#sig.AddElectronFilter("PreCC",eleS)
#sig.AddJetFilter("PreCC",jet)
sig.AddMuonFilter("PreCC",customMuID)
sig+=tree
#samples=[sd_jmt_skmd,mc_qcd15p8_skmd,sd_eg_skmd]
#samples=[QCD_Pt30_7TeV_Pythia]
#samples=[sd_jmt_skmd]
#samples=[sd_jmt_20100608]
#samples=[QCD_Pt30_7TeV_Pythia_skim]
#samples=[QCD_Pt30and80_7TeV_Pythia_skim]
#samples=[mc_qcd15p8_skmd]
samples=[Wjets_June2010_vols]
#scale to the desired luminosity
utils.reweight(samples,1.)
sig.Run("muTrees", conf, samples)
QCD_Pythia_AllPtBins_June2010,
QCD_38_Pythia6_Pt30to50,
QCD_38_Pythia6_Pt50to80,
QCD_38_Pythia6_Pt80to120,
QCD_38_Pythia6_Pt120to170,
QCD_38_Pythia6_Pt170to300,
QCD_38_Pythia6_Pt300to470,
QCD_38_Pythia6_Pt470to600,
QCD_38_Pythia6_Pt600to800,
QCD_38_Pythia6_Pt800to1000,
QCD_38_Pythia6_Pt1000to1400,
QCD_38_Pythia6_Pt1400to1800,
QCD_38_Pythia6_Pt1800toInf
]
samples = [
mu38X_A_1.Mu_Run2010A_CMSSW38X_1
# mu38X_B_1.Mu_Run2010B_CMSSW38X_1,
# mu38X_B_2.Mu_Run2010B_CMSSW38X_2,
# mu38X_B_3.Mu_Run2010B_CMSSW38X_3,
# mu38X_B_4.Mu_Run2010B_CMSSW38X_4,
# mu38X_B_5.Mu_Run2010B_CMSSW38X_5
]
utils.reweight(simplesMC, 1.)
sig.Run("/vols/cms02/gouskos/onelepton/20110120", conf, simplesMC)
#BatchRun(samples,"sig","conf","/vols/cms02/gouskos/onelepton/Muons_MuDS_20101026",50000)
# Global configuration
data = False
mc_int_lumi = 10.
conf = ewpol_conf.copy()
e = eWPol(data, loose_id, tight_id)
a = Analysis("eWPol", )
e.installFilters(a)
e.buildTree()
skim = fwk.SkimOp( PSet(
SkimName = "eWPol_MHT50Skim",
DropBranches = False,
Branches = types.StringVector([])
).ps() )
e.tree.TAttach(e.ops[len(e.ops)-1][1], skim)
a += e.tree
lm0 = PSet(
Name = "LM0",
Format = ("ICF", 2),
File = glob.glob("/home/alex/Research/Samples/ICF/LM0/SusyCAF*.root")
)
samples = [lm0]
utils.reweight(samples, mc_int_lumi)
a.Run("../results", conf, samples)
"/vols/cms02/gouskos/QCD_Pythia_Pt15_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt30_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt80_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt170_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt300_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt470_Jun2010/SusyCAF_Tree_*.root",
"/vols/cms02/gouskos/QCD_Pythia_Pt800_Jun2010/SusyCAF_Tree_*.root",
),
Weights=PSet(
CrossSection=[8.762e08, 6.041e07, 9.238e05, 2.547e04, 1.256e03, 8.798e01, 2.186, 0.01122],
Events=[6095857, 5069664, 2065792, 3171950, 2976108, 2159497, 2181700, 1185024],
PtBin=[15.0, 30.0, 80.0, 170.0, 300.0, 470.0, 800.0, 1400.0],
),
)
utils.reweight([LM0, LM1, LM01, LM01p, LMParis, QCD_AllPtBins_7TeV_Pythia, zjets_madgraph, bce_80to170], 10.0)
# Operations
trigger = OP_TwoTriggerCut("HLT_HT100U", "HLT_HT120U")
# trigger = OP_TriggerCut("HLT_Mu9")
selection = OP_GoodEventSelection()
HBHE = OP_HadronicHBHEnoiseFilter()
EleTreeReq = OP_NumComElectrons(">=", 1)
LeadingJetEta = OP_FirstJetEta(2.4)
unCorLeadJetCut = OP_UnCorLeadJetCut(70.0)
PrintOut = PrintAlphaTEvents()
OneMuon = OP_NumComMuons("==", 1)
MaxPtCut = OP_UpperElePtCut(25.0)
ZeroElectrons = OP_NumComElectrons("<=", 0)
OneElectron = OP_NumComElectrons("==", 1)