예제 #1
0
파일: datasets.py 프로젝트: mayoub/hhdb
 def xsec_kfact_effic(self):
     global XSEC_CACHE_MODIFIED
     year = self.year % 1E3
     if self.datatype == DATA:
         return 1., 1., 1.
     if year in XSEC_CACHE and self.name in XSEC_CACHE[year]:
         log.warning("using cached cross section for dataset %s" % self.ds)
         return XSEC_CACHE[year][self.name]
     try:
         return xsec.xsec_kfact_effic(self.year, self.id)
     except KeyError:
         log.warning("cross section of dataset %s not available locally. "
                     "Looking it up in AMI instead. AMI cross sections can be very"
                     "wrong! You have been warned!"
                     % self.ds)
     if USE_PYAMI:
         if self.ds in DS_NOPROV:
             xs, effic = get_dataset_xsec_effic(amiclient, DS_NOPROV[self.ds])
         else:
             xs, effic = get_dataset_xsec_effic(amiclient, self.ds)
         if year not in XSEC_CACHE:
             XSEC_CACHE[year] = {}
         XSEC_CACHE[year][self.name] = (xs, 1., effic)
         XSEC_CACHE_MODIFIED = True
         return xs, 1., effic
     raise Exception("cross section of dataset %s is not known!" % self.ds)
예제 #2
0
파일: datasets.py 프로젝트: mayoub/hhdb
 def xsec_kfact_effic(self):
     global XSEC_CACHE_MODIFIED
     year = self.year % 1E3
     if self.datatype == DATA:
         return 1., 1., 1.
     if year in XSEC_CACHE and self.name in XSEC_CACHE[year]:
         log.warning("using cached cross section for dataset %s" % self.ds)
         return XSEC_CACHE[year][self.name]
     try:
         return xsec.xsec_kfact_effic(self.year, self.id)
     except KeyError:
         log.warning(
             "cross section of dataset %s not available locally. "
             "Looking it up in AMI instead. AMI cross sections can be very"
             "wrong! You have been warned!" % self.ds)
     if USE_PYAMI:
         if self.ds in DS_NOPROV:
             xs, effic = get_dataset_xsec_effic(amiclient,
                                                DS_NOPROV[self.ds])
         else:
             xs, effic = get_dataset_xsec_effic(amiclient, self.ds)
         if year not in XSEC_CACHE:
             XSEC_CACHE[year] = {}
         XSEC_CACHE[year][self.name] = (xs, 1., effic)
         XSEC_CACHE_MODIFIED = True
         return xs, 1., effic
     raise Exception("cross section of dataset %s is not known!" % self.ds)
예제 #3
0
	sys.exit(1)

treeline = "TREENAME=%s\n" % treename
outfile.write(treeline)
for ds in dslist:

	isData = True
	if not datare.search(ds):
		isData = False

	xs = 0
	eff = 0
	slash = False
	if not isData:
		try:
			xs, eff = get_dataset_xsec_effic(client, ds+"/")
			slash = True
		except:
			try:
				xs, eff = get_dataset_xsec_effic(client, ds)
			except:	
				try:
					Slash = True
					orig = ds
					ds = re.sub(r'(e[0-9]*_)',r'merge.NTUP_COMMON.\1', ds) 
					print "trying to add NTUP_COMMON, like "+ds
					xs, eff = get_dataset_xsec_effic(client, ds+"/")
					#print xs
				except:
					print "Dataset", ds+"/", "not found in AMI, skipping."
					
예제 #4
0
파일: dataset_xsec.py 프로젝트: ndawe/pyAMI
from pyAMI.client import AMIClient
from pyAMI.query import get_dataset_xsec_effic

client = AMIClient()

dataset = 'mc11_7TeV.125206.PowHegPythia_VBFH130_tautauhh.merge.NTUP_TAUMEDIUM.e893_s1310_s1300_r2730_r2780_p795'
xsec, effic = get_dataset_xsec_effic(client, dataset)
print 'xsec: %.5E nb' % xsec
print 'effic: %.2f' % effic