Question: Goals: All projects are intended to make use what we learned over the course; in particular all projects will require you to implement an "event-like" ROOT Class and usage of container classes of your own ROOT object classes (depends on the problem). In particular the resonance problems (2-3 groups), since you should be able to analyze the resonance of the other group(s), requires that you discuss beforehand on a common event class interface. In principle you could implement it as you wish, but the interface should be in a generic way the you can use the resonance simulation data of the other group(s).
Relativistic Resonance Decay: Implement a relativistic resonance decay Monte Carlo Simulation: - Pick a resonance (will be assigned) - Implement decay in resonance rest frame (isotropic decay) - Boost in the center-of-mass frame of the collision use for example sqrt(sNN)=200 GeV - Save this event (of course lots of them) - Write an analysis macro: Invariant Mass calculation (to identify/confirm the resonance) - Additions: i) You can implement gaussian smearing of uncertainties in momentum (mimic detector resolution) ii) Add uncorrelated background PHY6860 Computational Physics, J. Putschke, Winter 2015 - Final Project Discussion - Relativistic Resonance Decay Cont. 5 Preparations: What is needed and how do you want to store/generate? Track class, Event class, containers/trees ... ? Prepare a short outline of what functionalities the classes need and a quick flow chart of the project itself! We will discuss it in class accordingly (April 13th)! Hints: You clearly need a 4 mom. vector (look at TLorentzVector), but it would be beneficial to extend this by creating your own track class allowing to store further information like particle type ... A good starting point for an event class we did in Lecture 10 ...
Project Tasks: Write a code solving the Pred-Prey ODE system - Create a root output file saving a TGraph (TH1) for the pred-prey population as function of time (allow for meaning file names and TGraph names depending on what parameters you used; in particular important ones we include population fluctuations) - Allow for pred-prey gaussian population fluctuations for each time step (Use ROOT TF1 in order create these fluctuations) - Write a separate plotting macro (reading in the different root files you created) and show several figures comparing the different solutions by varying the amount of population fluctuations). The figures should contain axis labeling a figure legend (use ROOT TLegend) ... (make it look pretty and professional!) - Additional: Think about other ways to characterize the distributions, mean, rms, ratios, 2d plots ... !? - Prepare a short presentation discussing your findings and an explanation of your code ...
Remark: This project is intended for you to learn how to define and outline a more complex computing program and in addition how to organize and work with other people on such a problem (so using OO programing you can think about that you can, once an interface is defined, indeed split work). It would be of course great if all projects will indeed be successful, but even if not, in the presentation you can discuss/point out the obstacles/problems and this is certainly part of the intended learning experience!