Ejemplo n.º 1
0
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00039.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00040.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00041.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00042.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00043.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00044.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00045.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00046.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00047.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00048.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00049.pool.root",
    # # "/nfs/atlas/nwarrack/user.dhirsch.AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV.evgen/AcerMC39.110069.singletop_tchan_lept_scale_115_CTEQ6L1_8TeV._00050.pool.root",
]

from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()

from Rivet_i.Rivet_iConf import Rivet_i

rivet = Rivet_i()
rivet.AnalysisPath = os.environ['PWD']
#rivet.Analyses += [ 'ATLAS_2017_I1512776' ]
rivet.Analyses += ['RivetMC_tchan_particle']
rivet.RunName = ""
#rivet.HistoFile = "MG5_atNLO_Herwig"
rivet.HistoFile = "Acer"
#rivet.OutputLevel = DEBUG
#rivet.CrossSection = 26.587 #<- xSec from https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopMC12SingleTopSamples#Singletop_Samples
rivet.CrossSection = 28.4348  #<- with kFactor
job += rivet
Ejemplo n.º 2
0
# /unix/bsmxsec2/zgrout/Samples/mc15_13TeV.364250.Sherpa_222_NNPDF30NNLO_llll.evgen.EVNT.e5894/
#svcMgr.EventSelector.InputCollections = glob.glob("/unix/bsmxsec2/zgrout/Samples/mc16_13TeV.345706.Sherpa_222_NNPDF30NNLO_ggllll_130M4l.merge.EVNT.e6213_e5984/EVNT.12842743._000027.pool.root.1")

# svcMgr.EventSelector.InputCollections = glob.glob("/home/goblirsc/Datasets/H4l/EVNT/mc15_13TeV.361603.PowhegPy8EG_CT10nloME_AZNLOCTEQ6L1_ZZllll_mll4.evgen.EVNT.e4475/*")[0:2]

from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()

from GaudiSvc.GaudiSvcConf import THistSvc
ServiceMgr += THistSvc()
ServiceMgr.THistSvc.Output = ["Rivet DATAFILE='Rivet.root' OPT='RECREATE'"]

from Rivet_i.Rivet_iConf import Rivet_i

rivet = Rivet_i()
rivet.AnalysisPath = os.environ['PWD']
#rivet.Analyses += [ 'ATLAS_2019_VBS' ]
rivet.Analyses += ['ATLAS_2019_I00001']
#rivet.Analyses += [ 'ATLAS_2019_I1720442' ]
rivet.Analyses += ['MC_XS']
rivet.RunName = ""

#rivet.WeightName=" nnlops-nominal "
#rivet.WeightName= ' nominal '
#rivet.WeightName = ' PDF set = 90420 '

rivet.HistoFile = "myanalysis"
rivet.DoRootHistos = True
#rivet.CrossSection = 9.1185E+03
job += rivet
Ejemplo n.º 3
0
from AthenaCommon.AppMgr import ServiceMgr as svcMgr
theApp.EvtMax = -1

import AthenaPoolCnvSvc.ReadAthenaPool
svcMgr.EventSelector.InputCollections = ["EVNT.05915947._001376.pool.root.1"]

from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()

from Rivet_i.Rivet_iConf import Rivet_i
rivet = Rivet_i()
rivet.Analyses += [ 'MC_BOOSTEDVBB' ]
rivet.RunName = ""
rivet.HistoFile = "boostedvbb.yoda"
import os
rivet.AnalysisPath = os.environ['PWD']
rivet.CrossSection = 1 
job += rivet 

from GaudiSvc.GaudiSvcConf import THistSvc

svcMgr += THistSvc() 
svcMgr.THistSvc.Output = ["Rivet DATAFILE='boostedvbb.root' OPT='RECREATE'"] 

job += rivet
Ejemplo n.º 4
0
    #sys.exit()

from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()

#OutputRoot="Rivet_"+fitt+".root"
from GaudiSvc.GaudiSvcConf import THistSvc
ServiceMgr += THistSvc()
ServiceMgr.THistSvc.Output = ["Rivet DATAFILE='Rivet.root' OPT='RECREATE'"]
ServiceMgr.MessageSvc.OutputLevel = 4 #INFO #DEBUG 

from Rivet_i.Rivet_iConf import Rivet_i

import os
rivet = Rivet_i()
rivet.AnalysisPath = os.environ['PWD']

#rivet.Analyses +=['ttw_ttH']
rivet.Analyses +=['ttwqq']
rivet.RunName = ""
# https://twiki.cern.ch/twiki/bin/view/AtlasProtected/RivetForAtlas#Use_with_Multiple_LHE_Weights
#rivet.WeightName="MUR0.5_MUF0.5_PDF261000"
#rivet.WeightName="MUR2_MUF2_PDF261000"
#rivet.HistoFile = "ttw_ttH"
rivet.HistoFile = "ttw_qq"
rivet.DoRootHistos=True
rivet.CrossSection = xs # xs to be read from athena or pathena commandline
print "CrossSection",xs
job += rivet
Ejemplo n.º 5
0
#M4--Best Limits
#svcMgr.EventSelector.InputCollections = ["/afs/cern.ch/user/n/ngubbi/work/AnalysisFiles-aQGC-WWZ/NewObservables-samples/DownloadedSamples/user.ngubbi.aQGC-M4-bestlimits-40-10pow12.10k-3lepevents.v1_EXT0/user.ngubbi.16195971.EXT0._000001.EVNT-aQGC-01.pool-M4-bestlimits-40-10pow12.10k-3lepevents.v1.root"]
#M2-Best Limits                                                                                                                                                                                
#svcMgr.EventSelector.InputCollections = ["/afs/cern.ch/user/n/ngubbi/work/AnalysisFiles-aQGC-WWZ/NewObservables-samples/DownloadedSamples/user.ngubbi.aQGC-M2-bestlimits-26-10pow12.10k-3lepevents.v1_EXT0/user.ngubbi.16192461.EXT0._000001.EVNT-aQGC-01.pool-M2-bestlimits-26-10pow12.10k-3lepevents.v1.root"]
#M3 Best Limits
#svcMgr.EventSelector.InputCollections = ["/afs/cern.ch/user/n/ngubbi/work/AnalysisFiles-aQGC-WWZ/NewObservables-samples/DownloadedSamples/user.ngubbi.aQGC-M3-bestlimits-43.5-10pow12.10k-3lepevents.v1_EXT0/user.ngubbi.16193004.EXT0._000001.EVNT-aQGC-01.pool-M3-bestlimits-43.5-10pow12.10k-3lepevents.v1.root"]
#T5 Best Limits
svcMgr.EventSelector.InputCollections = ["/afs/cern.ch/user/n/ngubbi/work/AnalysisFiles-aQGC-WWZ/NewObservables-samples/DownloadedSamples/user.ngubbi.aQGC-T5-bestlimits-3.8-10pow12.10k-3lepevents.v1_EXT0/user.ngubbi.16179036.EXT0._000001.EVNT-aQGC-01.pool-T5-bestlimits-3.8-10pow12.10k-3lepevents.v1.root"]

from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()

from Rivet_i.Rivet_iConf import Rivet_i

rivet = Rivet_i()
rivet.AnalysisPath = os.environ['PWD']
rivet.Analyses += [ 'WWZ' ]#The name of analysis here must be the same as the name of the class in thr rivet routine, else it will not know which analysis to run and also the respective rivet routine needs to be available in the same folder as this JO so that you can give PWD in the previous line for the analysis path. Else specify the path where rivet routine is located. 
rivet.RunName = ""
#rivet.HistoFile = "aQGC_rivetAnalysis-addnlObservables.yoda"
#rivet.HistoFile = OutputYoda
rivet.HistoFile = "WWZ"
#OutputYoda =  'WWZ-v2'
#OutputRoot  = 'WWZ-v2'

from AthenaCommon.AppMgr import ServiceMgr as svcMgr
from GaudiSvc.GaudiSvcConf import THistSvc
svcMgr += THistSvc()
svcMgr.THistSvc.Output = ["Rivet DATAFILE='"+OutputRoot+"' OPT='RECREATE'"]
 
job += rivet
                                         loc+head+'4332_00/EVNT.01404332._002820.pool.root.1',
                                         loc+head+'4328_00/EVNT.01404328._002641.pool.root.1',
                                         loc+head+'4328_00/EVNT.01404328._002604.pool.root.1',
                                         loc+head+'4327_00/EVNT.01404327._002576.pool.root.1',
                                         loc+head+'4325_00/EVNT.01404325._002476.pool.root.1',
                                         loc+head+'4325_00/EVNT.01404325._002461.pool.root.2',
                                         loc+head+'4325_00/EVNT.01404325._002459.pool.root.1',
                                         loc+head+'4324_00/EVNT.01404324._002427.pool.root.1',
                                         loc+head+'4324_00/EVNT.01404324._002419.pool.root.1',
                                         loc+head+'4225_00/EVNT.01404225._001488.pool.root.1',
                                         loc+head+'4225_00/EVNT.01404225._001474.pool.root.1',
                                         loc+head+'4218_00/EVNT.01404218._001111.pool.root.1',
                                         loc+head+'4202_00/EVNT.01404202._000337.pool.root.1',
                                         loc+head+'4197_00/EVNT.01404197._000082.pool.root.1',
                                         loc+head+'4197_00/EVNT.01404197._000080.pool.root.1']


from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()
from Rivet_i.Rivet_iConf import Rivet_i
rivet = Rivet_i()

# Choose which Rivet analyses to run (just the name of the analysis)
rivet.Analyses += ['ANALYSIS_NAME'] # You can add to this list, ['analysis1','analysis2'] as necessary

# You need to include the path of your own analysis (location of RivetANALYSIS_NAME.so) here so Rivet will find it:
rivet.AnalysisPath = '/path/to/analysis' # just use '.' for thise directory or maybe '..' or '../some/other/directory' etc
rivet.HistoFile = 'Rivet.yoda'

job += rivet
Ejemplo n.º 7
0
theApp.EvtMax = -1

import AthenaPoolCnvSvc.ReadAthenaPool
svcMgr.EventSelector.InputCollections = [
    'Samples/MoreSamples/mc15_13TeV.312484.Sh_227_NN30NNLO_Znunu_PTV100_140_MJJ0_500_KtMerging.evgen.EVNT.e7728/EVNT.18992357._000007.pool.root.1'
]

from AthenaCommon.AlgSequence import AlgSequence
job = AlgSequence()

from Rivet_i.Rivet_iConf import Rivet_i
rivet = Rivet_i()
import os
rivet.AnalysisPath = os.environ['PWD']

rivet.Analyses += ['Re_Vjj_vbf_nn_asym']
rivet.RunName = ''
rivet.HistoFile = 'File312484v2.yoda.gz'
rivet.CrossSection = 105.869824
#rivet.IgnoreBeamCheck = True
#rivet.SkipWeights=True
job += rivet
Ejemplo n.º 8
0
theApp.EvtMax = -1

import AthenaPoolCnvSvc.ReadAthenaPool

from AthenaCommon.AlgSequence import AlgSequence

job = AlgSequence()

from Rivet_i.Rivet_iConf import Rivet_i

rivet = Rivet_i()
import os

rivet.AnalysisPath = os.environ['PWD']

rivet.Analyses += ['Re_Vjj_vbf_nn_asym']
rivet.RunName = ''
rivet.HistoFile = '364222.yoda.gz'
rivet.CrossSection = 0.29612032
#rivet.IgnoreBeamCheck = True
#rivet.SkipWeights = True
job += rivet