예제 #1
0
    "/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.)
예제 #2
0
    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)
예제 #3
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)

예제 #4
0
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)
    
예제 #5
0
    "/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)
예제 #6
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)
예제 #7
0
    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)
예제 #8
0
# 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)
예제 #9
0
        "/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)