['user.damazio.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_0_etaBin_2.t0002_tunedDiscrXYZ.tgz/',
'user.damazio.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_1_etaBin_2.t0002_tunedDiscrXYZ.tgz/',
'user.damazio.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_2_etaBin_2.t0002_tunedDiscrXYZ.tgz/'],

['user.nbullacr.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_0_etaBin_3.t0001_tunedDiscrXYZ.tgz/',
'user.nbullacr.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_1_etaBin_3.t0001_tunedDiscrXYZ.tgz/',
'user.nbullacr.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_2_etaBin_3.t0001_tunedDiscrXYZ.tgz/'],
]

summaryList = [
'mc14_13TeV.147406.129160.sgn.offCutID.bkg.truth.trig.e24_lhmedium_L1EM20VH_summary.pic.gz',
'mc14_13TeV.147406.129160.sgn.offCutID.bkg.truth.trig.e24_lhtight_L1EM20VH_summary.pic.gz',
]

lhmediumSummary = load(summaryList[0])
lhtightSummary  = load(summaryList[1])



for nEta in [3]:

  for nEt in [0, 1, 2]:


    outputName = ('crossValStat_mc14_13TeV.147406.129160.sgn.offCutID.bkg.truth.trig.e24.etaBin_%d_etBin_%d')%(nEta,nEt)

    stat = CrossValidStatAnalysis( basepath+'/'+fileList[nEta][nEt], level = LoggingLevel.DEBUG )

    Medium_LH_L2Calo_Pd = ReferenceBenchmark( "Medium_LH_L2Calo_Pd", "Pd", refVal = lhmediumSummary['sgn'][nEt][nEta][0].efficiency()/float(100)  )
    Medium_LH_L2Calo_Pf = ReferenceBenchmark( "Medium_LH_L2Calo_Pf", "Pf", refVal = lhmediumSummary['bkg'][nEt][nEta][0].efficiency()/float(100)  )
예제 #2
0
        'user.damazio.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_1_etaBin_2.t0002_tunedDiscrXYZ.tgz/',
        'user.damazio.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_2_etaBin_2.t0002_tunedDiscrXYZ.tgz/'
    ],
    [
        'user.nbullacr.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_0_etaBin_3.t0001_tunedDiscrXYZ.tgz/',
        'user.nbullacr.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_1_etaBin_3.t0001_tunedDiscrXYZ.tgz/',
        'user.nbullacr.tuned.mc14.sgn.offLH.bkg.truth.trig.l1cluscut_20.l2etcut_19.e24_medium_etBin_2_etaBin_3.t0001_tunedDiscrXYZ.tgz/'
    ],
]

summaryList = [
    'mc14_13TeV.147406.129160.sgn.offCutID.bkg.truth.trig.e24_lhmedium_L1EM20VH_summary.pic.gz',
    'mc14_13TeV.147406.129160.sgn.offCutID.bkg.truth.trig.e24_lhtight_L1EM20VH_summary.pic.gz',
]

lhmediumSummary = load(summaryList[0])
lhtightSummary = load(summaryList[1])

for nEta in [3]:

    for nEt in [0, 1, 2]:

        outputName = (
            'crossValStat_mc14_13TeV.147406.129160.sgn.offCutID.bkg.truth.trig.e24.etaBin_%d_etBin_%d'
        ) % (nEta, nEt)

        stat = CrossValidStatAnalysis(basepath + '/' + fileList[nEta][nEt],
                                      level=LoggingLevel.DEBUG)

        Medium_LH_L2Calo_Pd = ReferenceBenchmark(
            "Medium_LH_L2Calo_Pd",
예제 #3
0
triggerList=[
              'e24_lhmedium_L1EM20VH_L2EFCalo_ringer_pd',
              'e24_lhmedium_L1EM20VH_L2Calo_ringer_pd',
              'e24_lhmedium_L1EM20VH_ringer_sp',
              'e24_lhtight_L1EM20VH_L2EFCalo_ringer_pd',
              'e24_lhtight_L1EM20VH_L2Calo_ringer_pd',
            ]
tc=0

for key in ['Medium_LH_EFCalo_Pd', 'Medium_LH_L2Calo_Pd','Medium_LH_MaxSP','Tight_LH_EFCalo_Pd','Tight_LH_L2Calo_Pd']:

  netBins=dict()
  for nEta in [0,1,2,3]:  
    for nEt in [0,1,2]:

      crossValSummary = load((crossValBasepath+'/'+crossValName)%(nEta,nEt))
      stat = CrossValidStatAnalysis( basepath+'/'+fileList[nEta][nEt], level = LoggingLevel.DEBUG )
      # Add them to a list:
      data  = CrossValidStatAnalysis.exportDiscrFiles(crossValSummary, RingerOperation.L2,
                                                      refBenchmarkNameList=[key],
                                                      configList=[config[key][nEta][nEt]])

      network = data['tunedDiscr']
      network['etaBin']=[etaBins[nEta],etaBins[nEta+1]] 
      network['etBin'] =[etBins[nEt]  ,etBins[nEt+1]  ] 
      netBins[('eta%d_et%d')%(nEta,nEt)]=network

  output.write(('  signatures["%s"]=%s\n')%(triggerList[tc],str(netBins)))
  tc+=1

output.write('  return signatures')