Notebook

It may take a few seconds

In [ ]:

!wget -q --show-progress --progress=dot:giga https://www.dropbox.com/s/zza28peyjy8qgg6/T2tt_700_50.root
!wget -q --show-progress --progress=dot:giga http://www.hepforge.org/archive/cutlang/atla_opendata_had_ttbar.root

CutLang parameters¶

file = Root file name (required)
filetype = Root file type (required) | ROOT file type can be:

LHCO FCC LVL0 DELPHES ATLASVLL ATLMIN ATLASOD ATLASODR2 CMSOD CMSNANO VLLBG3 VLLSIGNAL

adl | adlfile = ADL file name (default=cutlang-adl-name)
events = Numbers of events (number)
start = Starting event (number)
help = Help (bool)
deps = DEPS (string)
verbose = Verbose frequency (number)

Output file¶

Output file name will be histoOut-<adlfile>-<file>.root

In [ ]:

%%cutlang file=T2tt_700_50.root filetype=DELPHES adlfile=exHistos events=20000
# ADL file for ROOT101 tutorial: opening files and drawing histograms
# Runs with binder/ROOTintroPython.ipynb, binder/ROOTintroCpp.ipynb

# Object selection
object jets
  take jet
  select pT(jet) > 30
  select abs(eta(jet)) < 2.4

# Event selection
region presel
  select ALL
  select MET > 300
  select size(jets) >= 4
  histo hjet1pt , "jet 1 pT (GeV)", 40, 0, 1000, pT(jets[0])
  histo hjet2pt , "jet 2 pT (GeV)", 40, 0, 1000, pT(jets[1])
  histo hjet3pt , "jet 3 pT (GeV)", 40, 0, 1000, pT(jets[2])
  histo hjet4pt , "jet 4 pT (GeV)", 40, 0, 1000, pT(jets[3])
  histo hmetjet1pt, "MET vs. jet 1 pT", 20, 300, 1300, 20, 0, 1000, MET, pT(jets[0])

You can run the same analysis with different root files (separating each with ;) and add different parameters (separating each with ;) for each (optional)

In [ ]:

%%cutlang file=T2tt_700_50.root;atla_opendata_had_ttbar.root filetype=DELPHES;ATLASOD adlfile=exHistos events=5000;10000 start=0
# ADL file for ROOT101 tutorial: opening files and drawing histograms
# Runs with binder/ROOTintroPython.ipynb, binder/ROOTintroCpp.ipynb

# Object selection
object jets
  take jet
  select pT(jet) > 30
  select abs(eta(jet)) < 2.4

# Event selection
region presel
  select ALL
  select MET > 300
  select size(jets) >= 4
  histo hjet1pt , "jet 1 pT (GeV)", 40, 0, 1000, pT(jets[0])
  histo hjet2pt , "jet 2 pT (GeV)", 40, 0, 1000, pT(jets[1])
  histo hjet3pt , "jet 3 pT (GeV)", 40, 0, 1000, pT(jets[2])
  histo hjet4pt , "jet 4 pT (GeV)", 40, 0, 1000, pT(jets[3])
  histo hmetjet1pt, "MET vs. jet 1 pT", 20, 300, 1300, 20, 0, 1000, MET, pT(jets[0])