def theConfig(Nj, mcdir = '', initFile = '', btag = False, toydataFile=""):
    fitterPars = DibosonConfig.theConfig(Nj, mcdir, initFile, btag, toydataFile)
    fitterPars.includeElectrons = True
    fitterPars.includeMuons = False

    fitterPars.cuts += '&& (W_mt>50.) ' + \
                       '&& (event_met_pfmet>30) ' + \
                       '&& (W_electron_pt>35.) '

    fitterPars.QCDcuts = fitterPars.cuts
    fitterPars.cuts += '&& ( (abs(W_electron_eta)<1.5 && abs(W_electron_deltaphi_in)<0.03 && abs(W_electron_deltaeta_in)<0.004) || (abs(W_electron_eta)>1.5 && abs(W_electron_deltaphi_in)<0.02 && abs(W_electron_deltaeta_in)<0.005) ) '

    fitterPars.cuts += '&& ((W_electron_trackiso+W_electron_hcaliso+W_electron_ecaliso-event_RhoForLeptonIsolation*3.141592653589*0.09)/W_electron_pt<0.05) '


    print 'cuts:',fitterPars.cuts
    print 'QCD cuts:', fitterPars.QCDcuts
    
    return fitterPars
Example #2
0
def theConfig(Nj, mcdir='', initFile='', btag=False, toydataFile=""):
    fitterPars = DibosonConfig.theConfig(Nj, mcdir, initFile, btag,
                                         toydataFile)
    fitterPars.includeElectrons = True
    fitterPars.includeMuons = False

    fitterPars.cuts += '&& (W_mt>50.) ' + \
                       '&& (event_met_pfmet>30) ' + \
                       '&& (W_electron_pt>35.) '

    fitterPars.QCDcuts = fitterPars.cuts
    fitterPars.cuts += '&& ( (abs(W_electron_eta)<1.5 && abs(W_electron_deltaphi_in)<0.03 && abs(W_electron_deltaeta_in)<0.004) || (abs(W_electron_eta)>1.5 && abs(W_electron_deltaphi_in)<0.02 && abs(W_electron_deltaeta_in)<0.005) ) '

    fitterPars.cuts += '&& ((W_electron_trackiso+W_electron_hcaliso+W_electron_ecaliso-event_RhoForLeptonIsolation*3.141592653589*0.09)/W_electron_pt<0.05) '

    print 'cuts:', fitterPars.cuts
    print 'QCD cuts:', fitterPars.QCDcuts

    return fitterPars
def theConfig(Nj, mcdir = '', initFile = '', btag = False, toydataFile=""):
    fitterPars = DibosonConfig.theConfig(Nj, mcdir, initFile, btag, toydataFile)
    fitterPars.includeElectrons = False
    fitterPars.includeMuons = True

    fitterPars.cuts += '&& (W_mt>30.) ' + \
                       '&& (event_met_pfmet>25) ' + \
                       '&& (W_muon_pt>25.) ' + \
                       '&& (abs(W_muon_eta)<2.1) ' + \
                       '&& (abs(W_muon_d0bsp)<0.02) ' 

    fitterPars.QCDcuts = fitterPars.cuts
    
    fitterPars.cuts += '&& ((W_muon_trackiso+W_muon_hcaliso+W_muon_ecaliso-event_RhoForLeptonIsolation*3.141592653589*0.09)/W_muon_pt<0.1) '

    print 'cuts:',fitterPars.cuts
    print 'QCD cuts:', fitterPars.QCDcuts
    
    return fitterPars
def theConfig(Nj, mcdir='', initFile='', btag=False, toydataFile=""):
    fitterPars = DibosonConfig.theConfig(Nj, mcdir, initFile, btag,
                                         toydataFile)
    fitterPars.includeElectrons = False
    fitterPars.includeMuons = True

    fitterPars.cuts += '&& (W_mt>30.) ' + \
                       '&& (event_met_pfmet>25) ' + \
                       '&& (W_muon_pt>25.) ' + \
                       '&& (abs(W_muon_eta)<2.1) ' + \
                       '&& (abs(W_muon_d0bsp)<0.02) '

    fitterPars.QCDcuts = fitterPars.cuts

    fitterPars.cuts += '&& ((W_muon_trackiso+W_muon_hcaliso+W_muon_ecaliso-event_RhoForLeptonIsolation*3.141592653589*0.09)/W_muon_pt<0.1) '

    print 'cuts:', fitterPars.cuts
    print 'QCD cuts:', fitterPars.QCDcuts

    return fitterPars