> For the complete documentation index, see [llms.txt](https://antonio-iuliano1993.gitbook.io/software-notes/llms.txt). Markdown versions of documentation pages are available by appending `.md` to page URLs; this page is available as [Markdown](https://antonio-iuliano1993.gitbook.io/software-notes/other-software/useful-code-structures.md). # Useful code structures I will report here both the C++ and the python equivalent, when different ## ROOT ### Quickly converting numpy arrays to ROOT TTrees When working with a python numpy script, often we want to store the numbers in a ROOT TTree for easy sharing. Before, it was done with external libraries, such as uproot, rootpy, etc. Now ROOT itself supports an user-friendly conversion, by passing through a RDataframe: `energydata = {"Erec": Ehad_rec, "Etrue": y_test, "HitsTreeEntry": entrylist}` `energydf = ROOT.RDF.MakeNumpyDataFrame(energydata) energydf.Snapshot("energytree","energyfile.root")` ### Plotting numpy arrays in ROOT histograms Another common conversion problem is plotting scipy variables into ROOT TH1, TGraph, etc. usually for fitting, or sharing issues. The most straightforward way to do this is to use the **root\_numpy** library, which implements general functions, to fill histograms, profile histograms, and graphs. See for example {% embed url="" %} Note: for more than 1D plots, variables need to have the proper dimension. Remember how to increase dimension in python: ``` x = x[:, np.newsize] y = y[:, np.newsize] xy = np.concatenate([x,y], axis=1) root_numpy.fill_hist(h,xy) ``` ### Formatting variables to string Commonly used for nice printouts. We pass variables, with information about decimals and scientific notations. %d for integer, %.2f for floats with 2 decimals, %.2e etc * Python works in a fantastic way with `format` : here you can just write your string and a format, like `'testing with {:.2f}'.format(a)`. Now also directly with fstrings (but more strange to use) * C++ ROOT has a `TString::Format` which works similarly, but with % symbols, like TString::Format("testing with %.2f",a). Remember that a TString needs to be re-converted to chars with Data(), while for a std::string it is data() ### Numerical containers Used to perform easily selections and operations: * RVec containers are ROOT equivalent of numpy collections: `ROOT::VecOps::RVec v = {2,2,3}` and can be used like v\[v<2] * RDataFrame as equivalent of pandas dataframe ### TTree Indexing It is often useful to sort a TTree according to a particular branch, instead of using the default entries. Luckily, there is a builtin function to do that: **BuildIndex** mytree->BuildIndex("mybranch") creates an index according to mybranch, which identifies in a unique way an entry (if that is not a case, a minor second branch can be added). To the get the entry, then use mytree->GetEntryWithIndex("mybranch") ### TEfficiency The recommended way to produce a plot with the acceptances from the ratio of two histograms. Instead of the more intuitive TH1::Divide, this class automatically computes the errors in the expected binomial way. It is built with using: ``` TEfficiency *peff = new TEfficiency(&hpassed, &hall) ``` ROOT reference guide recommends to check consistency between histograms when using this class. All details can be found there (I want to have this page as clear as possible, just as a reminder of class names to look them in the ROOT reference). ### Shape of TCanvas Especially in 2D plots, it may be required to have exactly the same shape of both axes lengths. To do that, not only the range should be the same, but also the canvas size. So, you should set the TCanvas size as such: r.TCanvas("name","Title",800,800) ## BASH ### Loops Bash uses a list loop, like for example: ``` for VARIABLE in file1 file2 file3 do command1 on $VARIABLE command2 commandN done ``` a numerical loop can be done with seq, which uses the syntax FIRST INCREMENT LAST (see man seq): ``` for alpha in $(seq 0.02 0.01 0.08) do echo $alpha done ``` --- # Agent Instructions This documentation is published with GitBook. GitBook is the documentation platform designed so that both humans and AI agents can read, navigate, and reason over technical content effectively. Learn more at gitbook.com. ## Querying This Documentation If you need additional information that is not directly available in this page, you can query the documentation dynamically by asking a question. Perform an HTTP GET request on the current page URL with the `ask` query parameter, and the optional `goal` query parameter: ``` GET https://antonio-iuliano1993.gitbook.io/software-notes/other-software/useful-code-structures.md?ask=&goal= ``` `ask` is the immediate question: it should be specific, self-contained, and written in natural language. `goal` is optional and describes the broader end goal you are ultimately trying to accomplish on behalf of the user. GitBook uses it to tailor the answer towards what is most useful for that goal. The response will contain a direct answer to the question and relevant excerpts and sources from the documentation. Use this mechanism when the answer is not explicitly present in the current page, you need clarification or additional context, or you want to retrieve related documentation sections.