WEBVTT 1 00:00:00.930 --> 00:00:03.030 Edgar Fernando Carrera Jarrin: lesson for analysis flow. 2 00:00:04.710 --> 00:00:06.990 Edgar Fernando Carrera Jarrin: And we have to go through the solution. 3 00:00:08.519 --> 00:00:10.590 Edgar Fernando Carrera Jarrin: I don't know if you had the chance to. 4 00:00:12.360 --> 00:00:18.540 Edgar Fernando Carrera Jarrin: To work on the on the on the final assignment on the final challenge, which was to. 5 00:00:19.650 --> 00:00:26.880 Edgar Fernando Carrera Jarrin: derive this Cassidy contribution to the analysis, we were performing yesterday. 6 00:00:29.130 --> 00:00:53.880 Edgar Fernando Carrera Jarrin: take note of the time, so this 10 past the hour Okay, and so just let me remind you what we were up to last time or during the analysis flow lesson and essentially what we were trying to do is to go through all the different steps of the regular physics analysis. 7 00:00:55.980 --> 00:00:57.300 Edgar Fernando Carrera Jarrin: So, if you remember. 8 00:00:58.800 --> 00:01:00.960 Edgar Fernando Carrera Jarrin: This blood that I had here somewhere. 9 00:01:02.490 --> 00:01:03.810 Edgar Fernando Carrera Jarrin: So many lessons. 10 00:01:05.550 --> 00:01:06.120 You know. 11 00:01:07.620 --> 00:01:17.670 Edgar Fernando Carrera Jarrin: Here, so this is the blood that I had, and we have started from our a God samples add data sets they are quite heavy. 12 00:01:18.450 --> 00:01:40.560 Edgar Fernando Carrera Jarrin: contain a lot of information now, not all of it is good for your analysis, once you pick a final process final signal, then you would like to extract only the information needed in order to perform your study, and for that we use these poet, which is a very modular and very simple. 13 00:01:42.480 --> 00:01:53.460 Edgar Fernando Carrera Jarrin: put together a code that you can use to filter essentially of those events, and if you have enough computing resources that can be done. 14 00:01:54.660 --> 00:01:59.370 Edgar Fernando Carrera Jarrin: Depending on what you use can be done in in a few hours, perhaps. 15 00:02:00.750 --> 00:02:14.940 Edgar Fernando Carrera Jarrin: And then you end up with the either route file, which is flat that's what we are using in this tutorial or you can choose to dump your data information in your own form okay. 16 00:02:15.420 --> 00:02:26.670 Edgar Fernando Carrera Jarrin: So just reference, the group at MIT had their own MIT open data format and so on, and they have actually their their code somewhere. 17 00:02:27.840 --> 00:02:29.820 Edgar Fernando Carrera Jarrin: On the Web available. 18 00:02:31.320 --> 00:02:34.530 Edgar Fernando Carrera Jarrin: But you can choose whatever you you feel more comfortable with. 19 00:02:35.730 --> 00:02:37.980 Edgar Fernando Carrera Jarrin: And then the idea was that. 20 00:02:39.480 --> 00:02:41.070 Edgar Fernando Carrera Jarrin: Either at the poet. 21 00:02:42.090 --> 00:02:58.380 Edgar Fernando Carrera Jarrin: stage, you can apply calibrations and corrections as scotty i'm sorry i'm Julie, was a teaching you how to do that, so you can apply those directly at the skimmer at the poet using cms w or you can. 22 00:02:59.400 --> 00:03:03.390 Edgar Fernando Carrera Jarrin: apply that after that's also possible you would need some. 23 00:03:05.040 --> 00:03:09.510 Edgar Fernando Carrera Jarrin: Perhaps some files, with the corrections and you need to apply them afterwards. 24 00:03:10.500 --> 00:03:22.560 Edgar Fernando Carrera Jarrin: And then, once you have those enter balls goes flat enter both a one thing that you you'd like to do is to go over them again over all the events and apply. 25 00:03:23.520 --> 00:03:43.740 Edgar Fernando Carrera Jarrin: more sophisticated selection criteria and then you would obtain file with histograms and then you can just keep this histogram file and do all sorts of plots whatever you need in order to interpret your your data, something which we haven't done yet, but. 26 00:03:45.480 --> 00:03:47.220 Edgar Fernando Carrera Jarrin: This is sort of the full. 27 00:03:48.390 --> 00:03:53.370 Edgar Fernando Carrera Jarrin: Analysis flow scheme that I had we had been going through. 28 00:03:55.110 --> 00:04:08.580 Edgar Fernando Carrera Jarrin: Okay, so it's time for the also just to remind you what the task was about so we went through all this and one important thing is that we had very different. 29 00:04:10.170 --> 00:04:13.290 Edgar Fernando Carrera Jarrin: backgrounds coming back to this. 30 00:04:17.190 --> 00:04:18.240 Edgar Fernando Carrera Jarrin: Right and. 31 00:04:21.300 --> 00:04:22.800 we're on my backgrounds. 32 00:04:24.210 --> 00:04:26.550 Edgar Fernando Carrera Jarrin: Here, so we have different backgrounds. 33 00:04:30.450 --> 00:04:44.400 Edgar Fernando Carrera Jarrin: draw young to the bar w plus jets those were our main backgrounds, but something that we didn't consider explicitly here excuse ED that's a kind of background, which is very hard to deal with. 34 00:04:45.420 --> 00:04:46.230 Edgar Fernando Carrera Jarrin: And so. 35 00:04:48.330 --> 00:04:59.250 Edgar Fernando Carrera Jarrin: Many times experimental is derived this day, it is a contribution directly from data and that's that was that was the idea of the of the project. 36 00:05:01.500 --> 00:05:04.920 Edgar Fernando Carrera Jarrin: Because we had a little problem, remember, we. 37 00:05:06.090 --> 00:05:11.640 Edgar Fernando Carrera Jarrin: When we finished our our analysis flow we ended up with some. 38 00:05:12.810 --> 00:05:14.130 Edgar Fernando Carrera Jarrin: Sorry, the one. 39 00:05:15.270 --> 00:05:22.140 Edgar Fernando Carrera Jarrin: Here, with we ended up with some plots which didn't really look that correct. 40 00:05:23.400 --> 00:05:36.090 Edgar Fernando Carrera Jarrin: You know, we have initially we were missing some contribution, there is not a CD here yeah you have to blame it on the multijet contribution, and so the challenge was to derive that. 41 00:05:37.350 --> 00:05:43.560 Edgar Fernando Carrera Jarrin: So the I gave you a few hints I don't know if they were that useful. 42 00:05:45.660 --> 00:05:51.720 Edgar Fernando Carrera Jarrin: If you were in really in touch with this analysis or this way of. 43 00:05:53.910 --> 00:06:00.990 Edgar Fernando Carrera Jarrin: deriving things it's probably too difficult, but if you are familiar with these tools, maybe we're so busy I don't know. 44 00:06:02.100 --> 00:06:05.550 Edgar Fernando Carrera Jarrin: In any case, the idea was the following. 45 00:06:07.650 --> 00:06:20.940 Edgar Fernando Carrera Jarrin: I was trying to explain during this tutorial and it is analysis this this simplified analysis was already performed and that's where all these blocks come from right so. 46 00:06:23.190 --> 00:06:31.230 Edgar Fernando Carrera Jarrin: I gave you the links actually to the to the code, this was done with this art our data frame. 47 00:06:32.370 --> 00:06:35.640 Edgar Fernando Carrera Jarrin: Class from from route, which is a very powerful tool. 48 00:06:36.870 --> 00:06:51.510 Edgar Fernando Carrera Jarrin: And one of those scripts actually the script that is used in order to produce these plots was pointed here, and particularly this comment which said that you could. 49 00:06:53.820 --> 00:06:54.630 Edgar Fernando Carrera Jarrin: That you could. 50 00:06:55.680 --> 00:07:13.290 Edgar Fernando Carrera Jarrin: estimate that you see the background, by doing simple trick, essentially, which is is done here, this is Python, so it is probably easy enough to to understand it just by reading it, and what they do is essentially. 51 00:07:14.400 --> 00:07:30.060 Edgar Fernando Carrera Jarrin: estimating the CCD background, by taking into account all the different the other contributors to the background processes and then subtracting so they subtract the. 52 00:07:31.860 --> 00:07:42.300 Edgar Fernando Carrera Jarrin: The contribution from from the data itself okay that's the idea and so that's what they do here in this piece this block of code. 53 00:07:43.560 --> 00:07:49.230 Edgar Fernando Carrera Jarrin: But for that So the idea was I don't know if it was too hard and has that, as I said, are too easy. 54 00:07:49.470 --> 00:07:58.440 Edgar Fernando Carrera Jarrin: But the idea is that if you were trying to understand this, is that they get the histograms but they use a different a different tack for those histograms. 55 00:07:58.710 --> 00:08:12.300 Edgar Fernando Carrera Jarrin: So it meant as you go through this function of getting all the histograms from our route file, they were getting not the same histograms but but histograms from from with a different tack. 56 00:08:13.500 --> 00:08:16.560 Edgar Fernando Carrera Jarrin: Because, as it was mentioned in the comment. 57 00:08:18.240 --> 00:08:29.460 Edgar Fernando Carrera Jarrin: You had to define a control region in order to extract this contribution right, so you, you are not working, on the other way, is my camera turned on. 58 00:08:30.690 --> 00:08:35.160 Edgar Fernando Carrera Jarrin: So you're not you're not working on the signal region because. 59 00:08:36.870 --> 00:08:46.950 Edgar Fernando Carrera Jarrin: You don't want to submit your background directly from from this region you're taking a different phases space or different world of. 60 00:08:49.050 --> 00:08:59.250 Edgar Fernando Carrera Jarrin: Where the particles do their things and then you extract the contribution from that region so that's what that was sort of the idea and and. 61 00:09:00.480 --> 00:09:04.590 Edgar Fernando Carrera Jarrin: And I also I also pointed to you that that was done. 62 00:09:06.930 --> 00:09:10.440 Edgar Fernando Carrera Jarrin: So here's the whole package right it's called hicks toto. 63 00:09:12.630 --> 00:09:25.830 Edgar Fernando Carrera Jarrin: Now, known as the outreach analysis and it was done in three steps basically your skin, then you plug a you obtain the histograms and then you plot we simplify that to two steps. 64 00:09:26.730 --> 00:09:40.650 Edgar Fernando Carrera Jarrin: And so you could try to look around and figured out how to what is needed for those histograms to be in your histogram file, but at the end of the day. 65 00:09:41.370 --> 00:09:51.270 Edgar Fernando Carrera Jarrin: And i'm going to present the solution now, this is what you could what you could do so I have posted a solution here's the link, you can you can take a look, this is. 66 00:09:52.170 --> 00:10:05.520 Edgar Fernando Carrera Jarrin: Remember that we were working with the event loop analysis different versions of it, the idea was that we presented at first, a very simple version of it with just one variable. 67 00:10:06.360 --> 00:10:14.370 Edgar Fernando Carrera Jarrin: And even one data set, and then we expanded it we added more backgrounds, we made it is a monitor so so it can loop. 68 00:10:16.290 --> 00:10:30.150 Edgar Fernando Carrera Jarrin: Over all the data samples that we had already skimmed, of course, and it added more variables more histograms and then it was they were automatically or mostly automatically failed. 69 00:10:31.170 --> 00:10:32.670 Edgar Fernando Carrera Jarrin: So, here it is. 70 00:10:33.930 --> 00:10:45.540 Edgar Fernando Carrera Jarrin: also remember that this was, this is a very archived in some sense, a way of dealing with analysis actually does make glass is a very old glass. 71 00:10:47.640 --> 00:10:48.540 Edgar Fernando Carrera Jarrin: Of root. 72 00:10:50.340 --> 00:11:09.960 Edgar Fernando Carrera Jarrin: But you know it's very instructive I think I believe, if this is your first time trying to do analysis or particle physics analysis day and then, if you want to use the our data frame from the original approach or simply simplified approach, and it is probably easier. 73 00:11:10.110 --> 00:11:11.220 teeba: that that was the whole idea. 74 00:11:12.870 --> 00:11:13.650 Edgar Fernando Carrera Jarrin: Okay, so. 75 00:11:15.300 --> 00:11:36.390 Edgar Fernando Carrera Jarrin: here's the answer so initially we had all these histograms for the signal region defined okay and notice that I define here the histograms explicitly So you can see how they are constructed this, by the way, is the name of the histogram. 76 00:11:36.660 --> 00:11:39.000 Edgar Fernando Carrera Jarrin: This is just a title of the histogram. 77 00:11:39.750 --> 00:11:47.460 Edgar Fernando Carrera Jarrin: This is the real name of the histogram, this is the the tag by which it is called within c++. 78 00:11:49.230 --> 00:11:56.790 Edgar Fernando Carrera Jarrin: And then the name of the histogram and then you have some beans skin, which was the same as for the other simplified analysis. 79 00:11:58.020 --> 00:12:12.060 Edgar Fernando Carrera Jarrin: But if you if you look at the other version of the analysis and this was kind of hidden it is hidden in the our data frame because it's just a method in and you never is not very intuitive how. 80 00:12:12.510 --> 00:12:24.360 Edgar Fernando Carrera Jarrin: How they build this histogram so I just wanted to to explicitly put them out here, so you can see what's going on, but certainly one can find a way to to just. 81 00:12:25.740 --> 00:12:41.910 Edgar Fernando Carrera Jarrin: Do this definition in a more automatic way wish Okay, so this, these are the signal region, the signal histograms but then you define the control region we're going to see how that is done. 82 00:12:43.140 --> 00:12:54.120 Edgar Fernando Carrera Jarrin: And you repeat all of those histograms but for a different region Okay, this is our control region, and this is the ritual you're using in order to estimate the tct background. 83 00:12:55.140 --> 00:12:57.780 Edgar Fernando Carrera Jarrin: All right, and then the idea was that. 84 00:12:58.920 --> 00:13:00.960 Edgar Fernando Carrera Jarrin: You create another one of these. 85 00:13:03.360 --> 00:13:09.210 Edgar Fernando Carrera Jarrin: arrays of histograms so you can easily look over them and start filling them out. 86 00:13:10.260 --> 00:13:13.950 Edgar Fernando Carrera Jarrin: So I added this that should have done that well. 87 00:13:15.480 --> 00:13:21.180 Edgar Fernando Carrera Jarrin: The variables are the same because you're you're reading the same information from the root file. 88 00:13:22.950 --> 00:13:23.460 Edgar Fernando Carrera Jarrin: and 89 00:13:25.320 --> 00:13:32.160 Edgar Fernando Carrera Jarrin: What else any and here you have to load them say you know here at this stage what you're doing is is. 90 00:13:33.000 --> 00:13:47.100 Edgar Fernando Carrera Jarrin: letting the object, because the histograms are defined outside that is analysis object so you're letting the object know about the existence of the systems which were already defined so you're kind of copying them. 91 00:13:48.570 --> 00:14:00.270 Edgar Fernando Carrera Jarrin: Remember i'm not a c++ expert, so I probably do things in a weird way, like most of us so feel free to to change this at some point have. 92 00:14:02.280 --> 00:14:08.430 Edgar Fernando Carrera Jarrin: You guys, am I hearing you say yes right i'm afraid old was with this headphones that. 93 00:14:08.640 --> 00:14:10.230 Matthew Bellis: Can you can you hear me. 94 00:14:10.890 --> 00:14:11.700 Edgar Fernando Carrera Jarrin: Yes, yes, I can. 95 00:14:11.790 --> 00:14:13.050 Matthew Bellis: Okay, great so they're working. 96 00:14:14.040 --> 00:14:15.030 Edgar Fernando Carrera Jarrin: that's good to check. 97 00:14:15.960 --> 00:14:18.750 Edgar Fernando Carrera Jarrin: sorry about yesterday, yet so. 98 00:14:20.070 --> 00:14:21.390 Edgar Fernando Carrera Jarrin: yeah that's what i'm doing here. 99 00:14:22.410 --> 00:14:23.820 Edgar Fernando Carrera Jarrin: Okay, and. 100 00:14:25.470 --> 00:14:37.470 Edgar Fernando Carrera Jarrin: And yeah you just pass a copy of those of all those histograms and, of course, if you did that, for the signal you have to do that for the control region as well, so you have to define, you have to. 101 00:14:39.150 --> 00:14:56.040 Edgar Fernando Carrera Jarrin: populate this race as well, so Okay, this is this, this can also be done automatically, I suppose, once you do this automatically you can do the other thing automatically to but, again, I just wanted to expose this, so it is, it is more understandable. 102 00:14:57.270 --> 00:14:57.780 Probably. 103 00:15:00.300 --> 00:15:16.380 Edgar Fernando Carrera Jarrin: don't forget don't forget I didn't probably mentioned this, while it was in the instructions, but now they were going through the solution it's probably worth repeating don't forget that the advantage, perhaps, of this poet. 104 00:15:17.700 --> 00:15:21.660 Edgar Fernando Carrera Jarrin: Technology is that I like that main technology. 105 00:15:23.580 --> 00:15:24.270 Edgar Fernando Carrera Jarrin: is to. 106 00:15:25.290 --> 00:15:42.510 Edgar Fernando Carrera Jarrin: is to be able to perform this study modular early So if you don't care about meals, you can just throw them away you don't put them in the in the path in the cms sw path and and yeah they they don't they don't pull it your your sample your your final histogram. 107 00:15:43.980 --> 00:15:46.260 Edgar Fernando Carrera Jarrin: And the idea here is that. 108 00:15:47.400 --> 00:16:03.840 Edgar Fernando Carrera Jarrin: If you do some change to it, it looks like the the most complicated tree should go first always so that's why the trigger goes first, because the trigger has this map which has as its first element, the. 109 00:16:04.890 --> 00:16:15.060 Edgar Fernando Carrera Jarrin: The name of the trigger and as its second element in acceptance bit i'm kind of like a mixture of acceptance and and freescale. 110 00:16:16.290 --> 00:16:30.690 Edgar Fernando Carrera Jarrin: And then you, you can add the rest of the histograms which live in different sub directories in your word file and then you have to make them France so friendship is performed here, and then you can read all the variables, you have. 111 00:16:32.490 --> 00:16:39.960 Edgar Fernando Carrera Jarrin: Okay, so and then and then everything has to be in agreement with with those changes. 112 00:16:41.940 --> 00:16:58.980 Edgar Fernando Carrera Jarrin: And, but the whole the the bulk of the of the exercises here, of course, in the analysis part and we we define a few cuts remember Okay, the original selection was the most. 113 00:17:00.420 --> 00:17:10.440 Edgar Fernando Carrera Jarrin: harder district their selection was already done at the skimmer you remember We selected at least one noon and at least one tower the past certain requirement. 114 00:17:11.190 --> 00:17:25.920 Edgar Fernando Carrera Jarrin: But here we refined that we find neon tauber right, so we try to to to check or make sure we have the best 60 K into. 115 00:17:26.910 --> 00:17:43.830 Edgar Fernando Carrera Jarrin: detail but it went I went to a meal, on the case sonically in the other one hydroponically that's what we do here all right, we find the best PR and then we pose as Julia was mentioning and she was already talking about the transverse mass and how that. 116 00:17:46.680 --> 00:17:49.590 Edgar Fernando Carrera Jarrin: can help suppress this w plus jets. 117 00:17:51.150 --> 00:18:08.100 Edgar Fernando Carrera Jarrin: background and then yeah you guys came out with this isolated new on selection as well, for some reason, the analysis uses this Cone 03 but, if I remember correctly Julie mentioned that. 118 00:18:09.300 --> 00:18:19.290 Edgar Fernando Carrera Jarrin: Point four is using the actual analysis but anyway, we wanted to replicate this exercise so that's fine, for example. 119 00:18:21.240 --> 00:18:34.800 Edgar Fernando Carrera Jarrin: And here's where you define your control region right, and this is the addition that you should have done, if you understood this, so what we were doing is just fill in the histograms for the signal region. 120 00:18:36.030 --> 00:18:49.050 Edgar Fernando Carrera Jarrin: And we impose this additional requirement for the charge essentially we were saying Okay, these two towns have to have a different sign for the charge of course it is a hicks is. 121 00:18:50.220 --> 00:18:56.760 Edgar Fernando Carrera Jarrin: An intro particles that make sense, and it is very simple now to define. 122 00:18:58.590 --> 00:19:12.870 Edgar Fernando Carrera Jarrin: Control region because we are essentially requiring requiring the the opposite yeah because the top is identification is is independent of that, and then we we filled out called the. 123 00:19:14.100 --> 00:19:19.050 Edgar Fernando Carrera Jarrin: Now the underscore control region histograms. 124 00:19:20.100 --> 00:19:22.140 Edgar Fernando Carrera Jarrin: With the same selection, except for this one. 125 00:19:23.760 --> 00:19:39.240 Edgar Fernando Carrera Jarrin: And that that does the trick so at the end we're going to have the same set of histograms with that with these additional tax, which is what you see in the original platter So it goes for the. 126 00:19:40.890 --> 00:19:42.150 Edgar Fernando Carrera Jarrin: It goes for. 127 00:19:44.400 --> 00:19:46.200 Edgar Fernando Carrera Jarrin: For trying to read those histograms. 128 00:19:47.400 --> 00:19:54.150 Edgar Fernando Carrera Jarrin: Okay that's that's what it's done and it works, certainly, so we can try to test it. 129 00:19:55.560 --> 00:19:58.560 Edgar Fernando Carrera Jarrin: So here's my area, and when I go to the test. 130 00:20:01.470 --> 00:20:18.960 Edgar Fernando Carrera Jarrin: area, we were working on last time, and I could I could simply go to the to the solution and download my my file I should put it in the raw mode, because otherwise you you download a bunch of. 131 00:20:21.570 --> 00:20:41.910 Edgar Fernando Carrera Jarrin: html, and so let me copy this and you can exercise this if you have your wish, of course, you're very welcome to try it and yeah I have it, and I have to compile it remember it's always this line, by the way, I don't know if you noticed but. 132 00:20:43.440 --> 00:20:49.170 Edgar Fernando Carrera Jarrin: In the end, the latest version or the most sophisticated version of this. 133 00:20:50.190 --> 00:21:05.790 Edgar Fernando Carrera Jarrin: Of this program we had to add these other tech for the for the compiler was the same for the first few ones that were in that complicated, but after we had this. 134 00:21:07.290 --> 00:21:16.260 Edgar Fernando Carrera Jarrin: library which no i'm not sure if I needed it but anyway it didn't combine unless you you put this extra time. 135 00:21:18.630 --> 00:21:20.190 Edgar Fernando Carrera Jarrin: So that's what i'm going to use. 136 00:21:23.850 --> 00:21:26.580 Edgar Fernando Carrera Jarrin: That library had something some function, there was needed. 137 00:21:28.200 --> 00:21:30.090 Edgar Fernando Carrera Jarrin: And that's what why we included it. 138 00:21:31.290 --> 00:21:34.050 Edgar Fernando Carrera Jarrin: Okay, it combined fine and then you run this. 139 00:21:41.700 --> 00:21:44.910 yeah but I switch switch this. 140 00:21:55.110 --> 00:21:56.250 You know it's picking up. 141 00:22:01.020 --> 00:22:08.010 Edgar Fernando Carrera Jarrin: Again, the problem that I don't remember her name but she made us realize that this was. 142 00:22:10.620 --> 00:22:12.420 Edgar Fernando Carrera Jarrin: needed to be changed at some point. 143 00:22:13.980 --> 00:22:15.180 Edgar Fernando Carrera Jarrin: Finding my solution. 144 00:22:20.790 --> 00:22:26.190 Edgar Fernando Carrera Jarrin: I came back to because I I love down when was the trick again. 145 00:22:28.830 --> 00:22:31.680 Edgar Fernando Carrera Jarrin: Let me grab that trick serious strict. 146 00:22:32.880 --> 00:22:33.510 Edgar Fernando Carrera Jarrin: The quake. 147 00:22:34.950 --> 00:22:36.660 Edgar Fernando Carrera Jarrin: can see the core thing. 148 00:22:37.950 --> 00:22:40.980 Edgar Fernando Carrera Jarrin: was very simple right just export. 149 00:22:46.230 --> 00:22:49.440 Edgar Fernando Carrera Jarrin: All right, and i'm gonna I can't read it. 150 00:22:58.560 --> 00:23:00.270 Edgar Fernando Carrera Jarrin: Alright, so let me change this. 151 00:23:03.150 --> 00:23:03.540 yeah. 152 00:23:09.930 --> 00:23:12.540 Edgar Fernando Carrera Jarrin: Now I can execute this. 153 00:23:13.590 --> 00:23:14.520 Edgar Fernando Carrera Jarrin: and 154 00:23:16.500 --> 00:23:19.680 Edgar Fernando Carrera Jarrin: For some reason i'm not working in the same. 155 00:23:22.620 --> 00:23:23.610 Edgar Fernando Carrera Jarrin: I have so many. 156 00:23:25.980 --> 00:23:27.450 Edgar Fernando Carrera Jarrin: know, I have some many. 157 00:23:28.620 --> 00:23:30.990 Edgar Fernando Carrera Jarrin: containers that I didn't confuse. 158 00:23:32.040 --> 00:23:35.550 Edgar Fernando Carrera Jarrin: But that not a problem, I think I have them. 159 00:23:36.630 --> 00:23:37.260 Edgar Fernando Carrera Jarrin: here. 160 00:23:38.790 --> 00:23:39.960 Edgar Fernando Carrera Jarrin: In my test. 161 00:23:42.690 --> 00:23:45.360 Edgar Fernando Carrera Jarrin: copy of the files again. 162 00:23:46.650 --> 00:23:48.780 Edgar Fernando Carrera Jarrin: Because I don't know what I did. 163 00:23:50.070 --> 00:23:51.120 Edgar Fernando Carrera Jarrin: I move them somewhere. 164 00:23:52.440 --> 00:24:00.210 Cecilia Uribe Estrada: yeah I think that the a new you put the route files in the home directory and that's the reason why. 165 00:24:00.780 --> 00:24:01.890 Cecilia Uribe Estrada: yeah yeah. 166 00:24:02.460 --> 00:24:02.790 Maybe. 167 00:24:09.960 --> 00:24:10.920 I shouldn't take long. 168 00:24:12.690 --> 00:24:25.830 Edgar Fernando Carrera Jarrin: So when I don't know if you realize when you executed this event analysis program and any went through all the data sets I don't know if it gives me the chills to. 169 00:24:26.250 --> 00:24:38.040 Edgar Fernando Carrera Jarrin: think that you're actually going through the whole data set from 2012 and four, which many of us didn't sleep and there's so much for then collecting those data, and you know. 170 00:24:38.610 --> 00:24:47.130 Edgar Fernando Carrera Jarrin: You do this and less than a minute, and then you realize, everything is it's really I mean, at least for me is kind of. 171 00:24:49.590 --> 00:24:50.280 Edgar Fernando Carrera Jarrin: Oh sorry. 172 00:24:51.330 --> 00:24:55.800 Edgar Fernando Carrera Jarrin: yeah this I should read there right so. 173 00:24:57.690 --> 00:25:11.220 Edgar Fernando Carrera Jarrin: Not everything is perfect, you know I left out here this directory if you see here, this is Kim five So if you if you read this carefully, or you try to understand when the files are formed. 174 00:25:12.690 --> 00:25:13.590 Edgar Fernando Carrera Jarrin: They are for. 175 00:25:14.850 --> 00:25:20.970 Edgar Fernando Carrera Jarrin: Using samples like a base back right basically means, what is the. 176 00:25:22.560 --> 00:25:35.190 Edgar Fernando Carrera Jarrin: district, that you have to put behind the name of the file and since I was working in a different the container I was always always adding this directory that I had before, and so. 177 00:25:36.480 --> 00:25:37.560 Edgar Fernando Carrera Jarrin: I have to take that out. 178 00:25:38.910 --> 00:25:39.450 So. 179 00:25:43.080 --> 00:25:47.640 Edgar Fernando Carrera Jarrin: Sometimes we don't read our own errors that is. 180 00:25:58.620 --> 00:26:01.350 Edgar Fernando Carrera Jarrin: And I take this opportunity to mention that. 181 00:26:01.920 --> 00:26:13.770 Edgar Fernando Carrera Jarrin: If you didn't have the chance to as it was explained in the in the tutorial if you didn't have the chance to copy those files over you could also access them if you put if you enable this. 182 00:26:14.190 --> 00:26:26.280 Edgar Fernando Carrera Jarrin: descriptor here Okay, because these files are also in a public server at CERN and you can access them if you execute this program you can access them remotely through the X. 183 00:26:27.810 --> 00:26:48.810 Edgar Fernando Carrera Jarrin: X Ray D protocol, I found that in our container with disorder route version it takes a while like really take more than an hour so, but if you if you go to a more updated container with the with a newer route version, then things go go faster. 184 00:26:50.460 --> 00:26:53.190 Edgar Fernando Carrera Jarrin: Okay, so let me save this. 185 00:26:54.390 --> 00:26:56.190 Edgar Fernando Carrera Jarrin: and compile again. 186 00:27:03.690 --> 00:27:14.520 Edgar Fernando Carrera Jarrin: And now it goes yeah, this is the fun part, you are going through the whole nother whole dataset because we have released only half of the 2012 data set but still I mean. 187 00:27:15.780 --> 00:27:21.240 Edgar Fernando Carrera Jarrin: If you do this for the geeks different leptons which we have the example also. 188 00:27:22.770 --> 00:27:26.880 Edgar Fernando Carrera Jarrin: Then, then you are rediscovering the higgs boson in like a minute. 189 00:27:28.410 --> 00:27:33.630 Edgar Fernando Carrera Jarrin: Very good, all right, and then, yes, you know you can. 190 00:27:35.310 --> 00:27:46.020 Edgar Fernando Carrera Jarrin: You I gave you actually I gave you the final version Just to give you a head start right, I gave you the final version of the of the clutter that I call. 191 00:27:49.650 --> 00:27:52.260 Edgar Fernando Carrera Jarrin: I called it, the same but it's a little bit different from the. 192 00:27:54.270 --> 00:28:01.860 Edgar Fernando Carrera Jarrin: From the original version, because of their boots version differences but it works and if I. 193 00:28:02.910 --> 00:28:03.750 Edgar Fernando Carrera Jarrin: apply this. 194 00:28:06.750 --> 00:28:19.800 Edgar Fernando Carrera Jarrin: Then, my my histograms are created for the variables that I that we have and the variables that we had are the three that are displayed in the original reference. 195 00:28:20.670 --> 00:28:44.100 Edgar Fernando Carrera Jarrin: We had the Ada of the towel, we have just as an example of the variables that are stored for the Tower particle and then we have the primary vertex which is somehow independent of these two objects and then the visible mass you know Julie was talking about the. 196 00:28:45.540 --> 00:29:00.540 Edgar Fernando Carrera Jarrin: The possibility for us to reconstruct the full mass of the hicks because of the neutrinos that the escape, and so we don't have that much information about them as a transverse information, and so we when we do is reconstruct this. 197 00:29:01.770 --> 00:29:02.790 Edgar Fernando Carrera Jarrin: festival mass. 198 00:29:03.870 --> 00:29:18.540 Edgar Fernando Carrera Jarrin: And thanks to senior have out that display to show you these these one of these plots or maybe I can show you just to make sure that you see that if we run without actually the blood. 199 00:29:19.560 --> 00:29:20.520 that were mentioned. 200 00:29:21.540 --> 00:29:23.700 Edgar Fernando Carrera Jarrin: But you have them all here in the solution. 201 00:29:26.130 --> 00:29:26.610 Edgar Fernando Carrera Jarrin: Okay. 202 00:29:28.710 --> 00:29:30.570 Edgar Fernando Carrera Jarrin: right they look like this. 203 00:29:35.370 --> 00:29:55.590 Edgar Fernando Carrera Jarrin: Smaller they look like this, so you, you have successfully determine the QC the background or estimated it more or less okay we're not showing here or we're not doing any systematics study, which is perhaps the most difficult part in physics analysis. 204 00:29:57.510 --> 00:30:01.590 Edgar Fernando Carrera Jarrin: But but it's already pretty convincing that the script the. 205 00:30:03.000 --> 00:30:05.220 Edgar Fernando Carrera Jarrin: script is going on in a correct way. 206 00:30:06.450 --> 00:30:11.550 Edgar Fernando Carrera Jarrin: Right, then we have the detail data data for the towel. 207 00:30:12.780 --> 00:30:14.430 Edgar Fernando Carrera Jarrin: Right and. 208 00:30:15.450 --> 00:30:18.720 Edgar Fernando Carrera Jarrin: For the hydroponic decay of town essentially. 209 00:30:20.040 --> 00:30:24.090 Edgar Fernando Carrera Jarrin: And then, this is the visible visible mass. 210 00:30:25.980 --> 00:30:28.860 Edgar Fernando Carrera Jarrin: Of course you don't see the peak. 211 00:30:29.880 --> 00:30:32.610 Edgar Fernando Carrera Jarrin: In the the mass that you. 212 00:30:33.750 --> 00:30:39.090 Edgar Fernando Carrera Jarrin: You know 125 just because we don't have that much information so it's lower. 213 00:30:41.160 --> 00:30:54.870 Edgar Fernando Carrera Jarrin: And, and also, we have this other two samples for the signal, so if we get it, you know if we were able to observe this we couldn't with this data set by the way. 214 00:30:55.980 --> 00:31:07.470 Edgar Fernando Carrera Jarrin: I guess they hicks to tau tau was observed in 2008, I think, or something like that, but anyway, this is a good exercise. 215 00:31:08.640 --> 00:31:12.540 Edgar Fernando Carrera Jarrin: Okay, and so I was asking you, at some point. 216 00:31:15.570 --> 00:31:28.500 Edgar Fernando Carrera Jarrin: Oh, and here's the comparison Okay, so I copied their reference of the original simplified analysis which is on on the on the right and our version is on the left. 217 00:31:30.960 --> 00:31:47.970 Edgar Fernando Carrera Jarrin: And yeah I have a comparison for all three variables, the original histogram actually makes up a lot more variables that you can you can try to reproduce those as well if you're interested in learning about a different approach to this exercise. 218 00:31:49.200 --> 00:31:58.560 Edgar Fernando Carrera Jarrin: But I asked you to notice 111 difference, or at least one difference between our analysis, what what are you, what do you say. 219 00:31:59.670 --> 00:32:02.070 Edgar Fernando Carrera Jarrin: What is the What is the difference. 220 00:32:10.050 --> 00:32:10.770 Dillon Fitzgerald: distinguish. 221 00:32:14.190 --> 00:32:19.620 Edgar Fernando Carrera Jarrin: Exactly that's the difference right, so they are both correct, I believe. 222 00:32:20.730 --> 00:32:25.080 Edgar Fernando Carrera Jarrin: But Stefan is able to distinguish between the. 223 00:32:26.160 --> 00:32:39.450 Edgar Fernando Carrera Jarrin: The the true TAO right the irreducible contribution from the from the towel from the towels from the CD game two truths house and the contribution which is a smaller, of course. 224 00:32:40.710 --> 00:32:47.970 Edgar Fernando Carrera Jarrin: The contribution from the fact doubts, whereas leptons wherever they are our other leptons or they're. 225 00:32:49.290 --> 00:32:54.240 Edgar Fernando Carrera Jarrin: charged electrons are not identified or are mistaken bye bye bye house. 226 00:32:56.100 --> 00:33:08.730 Edgar Fernando Carrera Jarrin: And he can do that because he had well first a different approach, but he had more information in there in his original scheme and i'm gonna. 227 00:33:09.420 --> 00:33:20.040 Edgar Fernando Carrera Jarrin: Just as an exercise i'm gonna i'm going to show you why we can't do this, why we can't distinguish them and it's a good opportunity to. 228 00:33:21.000 --> 00:33:32.880 Edgar Fernando Carrera Jarrin: To talk about particle particle matching so Okay, the first thing to notice is that these are simulations Okay, so these backgrounds come from simulations. 229 00:33:33.270 --> 00:33:45.960 Edgar Fernando Carrera Jarrin: And, and that means that we know the truth, the Monte Carlo truth, as we, as we say and and we know exactly which ones were towers in which was which ones were not towers, and so we can distinguish them. 230 00:33:46.980 --> 00:33:48.930 Edgar Fernando Carrera Jarrin: And this is essentially what. 231 00:33:50.550 --> 00:33:52.680 Edgar Fernando Carrera Jarrin: Stephen does because. 232 00:33:53.970 --> 00:33:57.210 Edgar Fernando Carrera Jarrin: At some point, you know, in his analysis. 233 00:33:58.440 --> 00:34:10.140 Edgar Fernando Carrera Jarrin: He checks for security checks for generator House so he makes that that check so that don't mind the different language here it doesn't really matter. 234 00:34:11.040 --> 00:34:26.430 Edgar Fernando Carrera Jarrin: Of the matters is what i'm going to show you, but he checks whether the talus rejects the region of this Dallas other generator level that's what they what he does so we can search for this for this function. 235 00:34:31.290 --> 00:34:36.690 Edgar Fernando Carrera Jarrin: And we if we do that we arrived here at the this check. 236 00:34:38.790 --> 00:34:39.240 Edgar Fernando Carrera Jarrin: and 237 00:34:48.210 --> 00:35:00.870 Edgar Fernando Carrera Jarrin: Right and okay it's kind of complicated how this works, but he checks, essentially in his in his original scheme, he had this variable okay. 238 00:35:01.560 --> 00:35:15.480 Edgar Fernando Carrera Jarrin: Which is something called gen part PG ID basically he's asking for the generator particle pgd to be of the index of the right index right through the analysis. 239 00:35:17.520 --> 00:35:32.640 Edgar Fernando Carrera Jarrin: And so he's had this information here and i'm going to show you how he obtained that but how he obtains that is in a different package is in his schemer and here's the skimmer the the analogy of the poet schemer. 240 00:35:34.740 --> 00:35:36.570 Edgar Fernando Carrera Jarrin: And that is. 241 00:35:37.650 --> 00:35:38.970 Edgar Fernando Carrera Jarrin: That is going to be. 242 00:35:40.140 --> 00:35:40.620 here. 243 00:35:43.170 --> 00:35:48.090 Edgar Fernando Carrera Jarrin: i'm showing you this because the code is very understandable, I believe. 244 00:35:50.430 --> 00:35:52.470 Edgar Fernando Carrera Jarrin: And so we can learn something about this. 245 00:36:19.980 --> 00:36:33.660 Edgar Fernando Carrera Jarrin: A section here in his schemer as a matter of fact, where he matches, for example, neon with the JET the generator particles Okay, so he has information about the best match, and so on. 246 00:36:34.890 --> 00:36:39.060 Edgar Fernando Carrera Jarrin: And you know, and you notice that he has dysfunction. 247 00:36:40.380 --> 00:36:59.310 Edgar Fernando Carrera Jarrin: find best visible match, and he passes all the interesting generator particles and the for momentum of of the moon, and he does this for all the other means you know and so compares all the meals, with all the generator particles essentially and find the. 248 00:37:00.630 --> 00:37:06.480 Edgar Fernando Carrera Jarrin: The best visible match so i'm going to show you how this is done. 249 00:37:09.270 --> 00:37:10.350 Edgar Fernando Carrera Jarrin: Because the function. 250 00:37:11.640 --> 00:37:13.110 Edgar Fernando Carrera Jarrin: If those functions are. 251 00:37:14.460 --> 00:37:15.180 So we're. 252 00:37:30.300 --> 00:37:30.660 here. 253 00:37:32.460 --> 00:37:32.940 Okay. 254 00:37:34.500 --> 00:37:45.360 Edgar Fernando Carrera Jarrin: So what he does is, as I said, take the new one right there for momentum of me on take all the generator particles. 255 00:37:46.680 --> 00:38:03.060 Edgar Fernando Carrera Jarrin: And he looks over all the generator particles and try to compare that or to make a geometrical match, you know, on delta are, I think you went through this with with Julie, you know this typical. 256 00:38:04.200 --> 00:38:13.260 Edgar Fernando Carrera Jarrin: Distance in this are kind of weird our space in fi and this, so the rapidity variable. 257 00:38:14.190 --> 00:38:30.900 Edgar Fernando Carrera Jarrin: And and yeah you have you find the best of the best match for that, but before he does this part of matching this geometrically right you look at the direction of this every constructed particle and you see all the particles that were generated that go through. 258 00:38:31.470 --> 00:38:34.470 Edgar Fernando Carrera Jarrin: The same region and find the best that the closest. 259 00:38:35.700 --> 00:38:48.930 Edgar Fernando Carrera Jarrin: Of those particles, but before doing that he's abstracts the invisible, of course, right because the momentum of the new trainers is not known only the transverse momentum is known and so. 260 00:38:49.500 --> 00:39:01.530 Edgar Fernando Carrera Jarrin: That that messes up the whole approach of matching geometrically and so he takes care of that actually because he has a function for subtracting visible. 261 00:39:02.040 --> 00:39:17.430 Edgar Fernando Carrera Jarrin: But in order to subtract investable you see what he does here, he looks over all the generator particles and he finds the daughters right in order to find those neutrinos and subtract they're there for momentum. 262 00:39:19.980 --> 00:39:39.240 Edgar Fernando Carrera Jarrin: And this is the piece of information, the information about the daughters that we have we haven't worked on in the end, the point is schema and that's the reason why we couldn't do this okay so Okay, this project is still in progress and we are certainly. 263 00:39:40.410 --> 00:39:50.970 Edgar Fernando Carrera Jarrin: find a way to store this information, but that was just the reason for it, I thought it was interesting I don't know if you enjoyed this exercise, but that's the reason why we kept this. 264 00:39:51.900 --> 00:40:05.550 Edgar Fernando Carrera Jarrin: Is not discrepancy, but a different way of doing yeah you know, so there is no way for us to distinguish them, because we cannot do the matching at this stage, but it can be done, certainly. 265 00:40:06.780 --> 00:40:07.380 Edgar Fernando Carrera Jarrin: All right. 266 00:40:09.600 --> 00:40:15.060 Edgar Fernando Carrera Jarrin: yeah and I think I went over my time apologize apologies. 267 00:40:16.680 --> 00:40:21.450 Edgar Fernando Carrera Jarrin: We should go for a break, I know we have to start our hands on tutorial right yeah. 268 00:40:23.460 --> 00:40:23.850 Edgar Fernando Carrera Jarrin: Okay. 269 00:40:25.020 --> 00:40:25.620 Edgar Fernando Carrera Jarrin: and 270 00:40:27.810 --> 00:40:28.380 let's go. 271 00:40:31.470 --> 00:40:35.130 Edgar Fernando Carrera Jarrin: Do you have any questions anything any. 272 00:40:40.620 --> 00:40:42.600 Edgar Fernando Carrera Jarrin: Any comments anything that's. 273 00:40:45.690 --> 00:40:49.110 Edgar Fernando Carrera Jarrin: that's concerning you about this analysis flow. 274 00:40:59.820 --> 00:41:04.740 Dillon Fitzgerald: Just to clarify I need this cms software for this given part yeah and then. 275 00:41:05.760 --> 00:41:08.550 Dillon Fitzgerald: import everything we've been through my local mission. 276 00:41:10.860 --> 00:41:13.260 Edgar Fernando Carrera Jarrin: that's correct that's correct, if you are one. 277 00:41:14.520 --> 00:41:16.200 Edgar Fernando Carrera Jarrin: route hater because. 278 00:41:16.440 --> 00:41:21.570 Edgar Fernando Carrera Jarrin: We are all half and half of my hater and I read lover. 279 00:41:23.070 --> 00:41:40.050 Edgar Fernando Carrera Jarrin: But you're a complete route hater, then yes you're free to get rid of cms is w, which is the most complicated part, perhaps, and we try to make it make it easier for you, with this poet tool that we keep updating of course still be improved. 280 00:41:41.490 --> 00:41:48.240 Edgar Fernando Carrera Jarrin: But then you're you're free, as I said, you can dump this information in whatever format you want and be done with it. 281 00:41:52.590 --> 00:41:53.250 Edgar Fernando Carrera Jarrin: that's right. 282 00:41:55.440 --> 00:41:55.950 Dillon Fitzgerald: Okay. 283 00:41:57.450 --> 00:41:58.140 Edgar Fernando Carrera Jarrin: Very good. 284 00:41:59.490 --> 00:42:05.430 Edgar Fernando Carrera Jarrin: So let's jump to the hands on tutorial for the cloud computing part. 285 00:42:06.510 --> 00:42:10.260 Matthew Bellis: Edgar i'm gonna stop this session and then i'll get the next session starting to record. 286 00:42:10.860 --> 00:42:12.930 Edgar Fernando Carrera Jarrin: yeah that sounds like a good idea, thank you man.