{"id":5987,"date":"2018-05-05T11:44:12","date_gmt":"2018-05-05T11:44:12","guid":{"rendered":"https:\/\/particlebites.com\/?p=5987"},"modified":"2018-05-05T11:44:12","modified_gmt":"2018-05-05T11:44:12","slug":"going-rogue-the-search-for-anything-and-everything-with-atlas","status":"publish","type":"post","link":"https:\/\/www.particlebites.com\/?p=5987","title":{"rendered":"Going Rogue: The Search for Anything (and Everything) with ATLAS"},"content":{"rendered":"<p><strong>Title<\/strong>: \u201cA model-independent general search for new phenomena with the ATLAS detector at\u00a0\u221as=13 TeV\u201d<\/p>\n<p><strong>Author:<\/strong> The ATLAS Collaboration<\/p>\n<p><strong>Reference:<\/strong>\u00a0<a href=\"http:\/\/inspirehep.net\/record\/1511481\">ATLAS-PHYS-CONF-2017-001<\/a><\/p>\n<p>&nbsp;<\/p>\n<p class=\"p1\"><span class=\"s1\">When a single experimental collaboration has a few thousand contributors (and even more opinions), there are a lot of rules. These rules dictate everything from how you get authorship rights to how you get chosen to give a conference talk. In fact, this rulebook is so thorough that it could be the topic of a whole other post. But for now, I want to focus on one rule in particular, a rule that has only been around for a few decades in particle physics but is considered one of the most important practices of good science: <strong>blinding<\/strong>. <\/span><\/p>\n<p class=\"p1\"><span class=\"s1\">In brief, blinding is the notion that it\u2019s experimentally compromising for a scientist to look at the data before finalizing the analysis. As much as we like to think of ourselves as perfectly objective observers, the truth is, when we really <i>really <\/i>want a particular result (let\u2019s say a SUSY discovery), that desire can bias our work. For instance, imagine you were looking at actual collision data while you were designing a signal region. You might <a href=\"https:\/\/www.nature.com\/news\/blind-analysis-hide-results-to-seek-the-truth-1.18510\">unconsciously craft your selection<\/a> in such a way to force an excess of data over background prediction. To avoid such human influences, particle physics experiments \u201cblind\u201d their analyses while they are under construction, and only look at the data once everything else is in place and validated. <\/span><\/p>\n<figure id=\"attachment_5988\" aria-describedby=\"caption-attachment-5988\" style=\"width: 512px\" class=\"wp-caption alignnone\"><a href=\"https:\/\/particlebites.com\/wp-content\/uploads\/2018\/05\/blindAnalysis.jpg\"><img loading=\"lazy\" decoding=\"async\" class=\" wp-image-5988\" src=\"https:\/\/particlebites.com\/wp-content\/uploads\/2018\/05\/blindAnalysis-300x157.jpg\" alt=\"\" width=\"512\" height=\"268\" srcset=\"https:\/\/www.particlebites.com\/wp-content\/uploads\/2018\/05\/blindAnalysis-300x157.jpg 300w, https:\/\/www.particlebites.com\/wp-content\/uploads\/2018\/05\/blindAnalysis.jpg 630w\" sizes=\"auto, (max-width: 512px) 100vw, 512px\" \/><\/a><figcaption id=\"caption-attachment-5988\" class=\"wp-caption-text\">Figure 1: &#8220;Blind analysis: Hide results to seek the truth&#8221;, R. MacCounor &amp; S. Perlmutter for Nature.com<\/figcaption><\/figure>\n<p class=\"p1\"><span class=\"s1\">This technique has kept the field of particle physics in rigorous shape for quite a while. But there\u2019s always been a subtle downside to this practice. If we only ever look at the data after we finalize an analysis, we are trapped within the confines of theoretically motivated signatures. In this blinding paradigm, we\u2019ll look at all the places that theory has shone a spotlight on, but we won\u2019t look <i>everywhere.<\/i> Our whole game is to search for new physics. But what if amongst all our signal regions and hypothesis testing and neural net classifications\u2026 <strong>we\u2019ve simply missed something?<\/strong><\/span><\/p>\n<p class=\"p1\"><span class=\"s1\">It is this nagging question that motivates a specific method of combing the LHC datasets for new physics, one that the authors of this paper call a <strong>\u201cstructured, global and automated way to search for new physics.\u201d<\/strong> With this proposal, we can let the data itself tell us where to look and throw unblinding caution to the winds. <\/span><\/p>\n<p class=\"p1\"><span class=\"s1\">The idea is simple: scan the whole ATLAS dataset for discrepancies, setting a threshold for what defines a feature as \u201cinteresting\u201d. If this preliminary scan stumbles upon a mysterious excess of data over Standard Model background, don\u2019t just run straight to Stockholm proclaiming a discovery. Instead, simply remember to look at this area again once more data is collected. If your feature of interest is a fluctuation, it will wash out and go away. If not, you can keep watching it until you collect enough statistics to do the running to Stockholm bit. Essentially, you let a first scan of the data rather than theory define your signal regions of interest. In fact, all the cool kids are doing it: <\/span><span class=\"s2\"><a href=\"https:\/\/arxiv.org\/abs\/hep-ex\/0408044\">H1<\/a><\/span><span class=\"s2\">, <\/span><span class=\"s2\"><a href=\"https:\/\/arxiv.org\/abs\/0809.3781\">CDF<\/a><\/span><span class=\"s2\">, <\/span><span class=\"s2\"><a href=\"https:\/\/arxiv.org\/abs\/hep-ex\/0011067\">D0<\/a><\/span><span class=\"s2\">, and even<\/span> <span class=\"s2\"><a href=\"http:\/\/inspirehep.net\/record\/1204310\">ATLAS<\/a><\/span> <span class=\"s2\">and<\/span><span class=\"s2\"> <a href=\"http:\/\/cdsweb.cern.ch\/record\/1360173\">CMS<\/a><\/span><span class=\"s2\"> have perform<\/span><span class=\"s2\">ed<\/span><strong><span class=\"s2\"> earlier versions <\/span><\/strong><span class=\"s2\">of this general search<\/span><span class=\"s2\">.<\/span><\/p>\n<p class=\"p1\"><span class=\"s1\">The nuts and bolts of this particular paper include 3.2 fb<\/span><span class=\"s2\"><sup>-1<\/sup><\/span><span class=\"s1\"> of 2015 13 TeV LHC data to try out. Since the whole goal of this strategy is to be as general as possible, we might as well go big or go home with potential topologies. To that end, the authors comb through all the data and select any event \u201cinvolving high pT isolated leptons (electrons and muons), photons, jets, b-tagged jets and missing transverse momentum\u201d.\u00a0All of the backgrounds are simply modeled with Monte Carlo simulation.<\/span><\/p>\n<p class=\"p1\"><span class=\"s1\">Once we have all these events, we need to sort them. Here, \u201cthe classification includes all possible final state configurations and object multiplicities, e.g. if a data event with seven reconstructed muons is found it is classified in a \u20187- muon\u2019 event class (7\u03bc).\u201d When you add up all the possible permutations of objects and multiplicities, you come up with a cool <strong>639 event classes<\/strong> with at least 1 data event and a Standard Model expectation of at least 0.1. <\/span><\/p>\n<p class=\"p1\"><span class=\"s1\">From here, it\u2019s just a matter of checking data vs. MC agreement and the <a href=\"http:\/\/physics.rockefeller.edu\/luc\/technical_reports\/cdf5776_pulls.pdf\">pulls<\/a> for each event class. The authors also apply some measures to weed out the low stat or otherwise sketchy regions; for instance, 1 electron + many jets is more likely to be multijet faking a lepton and shouldn\u2019t necessarily be considered as a good event category. Once this logic applied, you can plot all of your SRs together grouped by category; Figure 2 shows an example for the multijet events. The paper includes 10 of these plots in total, with regions ranging in complexity from nothing but 1<\/span><span class=\"s3\">\u03bc1j <\/span><span class=\"s1\">to more complicated final states like <i>E<\/i><\/span><span class=\"s2\"><i><sub>T<\/sub><\/i><\/span><span class=\"s4\">miss<\/span><span class=\"s1\">2\u03bc1\u03b34j (say that five times fast.)<\/span><\/p>\n<figure id=\"attachment_5991\" aria-describedby=\"caption-attachment-5991\" style=\"width: 607px\" class=\"wp-caption aligncenter\"><a href=\"https:\/\/particlebites.com\/wp-content\/uploads\/2018\/05\/multijets-1.png\"><img loading=\"lazy\" decoding=\"async\" class=\" wp-image-5991\" src=\"https:\/\/particlebites.com\/wp-content\/uploads\/2018\/05\/multijets-1-300x125.png\" alt=\"\" width=\"607\" height=\"253\" srcset=\"https:\/\/www.particlebites.com\/wp-content\/uploads\/2018\/05\/multijets-1-300x125.png 300w, https:\/\/www.particlebites.com\/wp-content\/uploads\/2018\/05\/multijets-1-768x319.png 768w, https:\/\/www.particlebites.com\/wp-content\/uploads\/2018\/05\/multijets-1-1024x426.png 1024w, https:\/\/www.particlebites.com\/wp-content\/uploads\/2018\/05\/multijets-1.png 1508w\" sizes=\"auto, (max-width: 607px) 100vw, 607px\" \/><\/a><figcaption id=\"caption-attachment-5991\" class=\"wp-caption-text\">Figure 2: The number of events in data and for the different SM background predictions considered. The classes are labeled according to the multiplicity and type (e, \u03bc, \u03b3, j, b, ETmiss) of the reconstructed objects for this event class. The hatched bands indicate the total uncertainty of the SM prediction.<\/figcaption><\/figure>\n<p>&nbsp;<\/p>\n<p class=\"p1\"><span class=\"s1\">Once we can see data next to Standard Model prediction for all these categories, it\u2019s necessary to have a way to measure just how unusual an excess may be. The authors of this paper implement an algorithm that searches for the region of largest deviation in the distributions of two variables that are good at discriminating background from new physics. These are the <\/span><span class=\"s2\"><b>effective mass<i>,\u00a0 <\/i><\/b><i><\/i>the sum of all jet and missing momenta,\u00a0<\/span><span class=\"s1\">and the <strong>invariant mass<\/strong>, computed with all visible objects and no missing energy. <\/span><\/p>\n<p class=\"p1\"><span class=\"s1\">For each deviation found, a simple likelihood function is built as the convolution of probability density functions (pdfs): one<a href=\"http:\/\/hyperphysics.phy-astr.gsu.edu\/hbase\/Math\/poiex.html\"> Poissonian<\/a> pdf to describe the event yields, and Gaussian pdfs for each systematic uncertainty. The integral of this function, <strong>p0<\/strong>, is the probability that the Standard Model expectation fluctuated to the observed yield. This p0 value is an industry standard in particle physics: a value of p0 &lt; 3e-7 is our threshold for discovery. <\/span><\/p>\n<p class=\"p1\"><span class=\"s1\">Sadly (or reassuringly), the <strong>smallest p0 value found in this scan is 3e-04<\/strong> (in the 1<i>m<\/i>1<i>e<\/i>4<i>b<\/i>2<i>j <\/i>event class). To figure out precisely how significant this value is, the authors ran a series of pseudoexperiments for each event class and applied the same scanning algorithm to them, to determine how often such a deviation would occur in a wholly different fake dataset. In fact, a p0 of 3e-04 was expected 70% of the pseudoexperiments. <\/span><\/p>\n<p class=\"p1\"><span class=\"s1\">So the excesses that were observed are not (so far) significant enough to focus on. But the beauty of this analysis strategy is that this deviation can be easily followed up with the addition of a newer dataset. Think of these general searches as the sidekick of the superheros that are our flagship SUSY, exotics, and dark matter searches. They can help us dot i\u2019s and cross t\u2019s, make sure nothing falls through the cracks\u2014 and eventually, just maybe, make a discovery.<\/span><\/p>\n<p class=\"p1\">\n","protected":false},"excerpt":{"rendered":"<p>Title: \u201cA model-independent general search for new phenomena with the ATLAS detector at\u00a0\u221as=13 TeV\u201d Author: The ATLAS Collaboration Reference:\u00a0ATLAS-PHYS-CONF-2017-001 &nbsp; When a single experimental collaboration has a few thousand contributors (and even more opinions), there are a lot of rules. These rules dictate everything from how you get authorship rights to how you get chosen &hellip; <\/p>\n<p class=\"link-more\"><a href=\"https:\/\/www.particlebites.com\/?p=5987\" class=\"more-link\">Continue reading<span class=\"screen-reader-text\"> &#8220;Going Rogue: The Search for Anything (and Everything) with ATLAS&#8221;<\/span><\/a><\/p>\n","protected":false},"author":8,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"aside","meta":{"footnotes":""},"categories":[47,57],"tags":[],"class_list":["post-5987","post","type-post","status-publish","format-aside","hentry","category-beyond-standard-model","category-experimental-techniques","post_format-post-format-aside"],"_links":{"self":[{"href":"https:\/\/www.particlebites.com\/index.php?rest_route=\/wp\/v2\/posts\/5987","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.particlebites.com\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.particlebites.com\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.particlebites.com\/index.php?rest_route=\/wp\/v2\/users\/8"}],"replies":[{"embeddable":true,"href":"https:\/\/www.particlebites.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=5987"}],"version-history":[{"count":11,"href":"https:\/\/www.particlebites.com\/index.php?rest_route=\/wp\/v2\/posts\/5987\/revisions"}],"predecessor-version":[{"id":6002,"href":"https:\/\/www.particlebites.com\/index.php?rest_route=\/wp\/v2\/posts\/5987\/revisions\/6002"}],"wp:attachment":[{"href":"https:\/\/www.particlebites.com\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=5987"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.particlebites.com\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=5987"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.particlebites.com\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=5987"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}