]> git.uio.no Git - u/mrichter/AliRoot.git/blob - PWG2/data/README
Last modifictaions for PWG2 -> PWGLF
[u/mrichter/AliRoot.git] / PWG2 / data / README
1 ===========================================================
2 =          Author: Panos.Christakoglou@cern.ch            =
3 ===========================================================
4
5
6 PriorPropabilities.root
7 _______________________
8
9 The file contains 5 TF1 objects namely: 
10 fitElectrons
11 fitMuons
12 fitPions
13 fitKaons
14 fitProtons
15
16 These objects represent the functions (pol9) used to fit the the momentum spectra of electrons, muons, pions, kaons and protons. The spectra were generated by looping over the primaries of the kine trees (generator PYTHIA - production PDC07).
17
18 This file can be used as input for any PID related study if one needs to have a momentum dependence of the a priori particle concentrations. This method gives better results (higher efficiency) with increasing momenta.
19
20 To get the momentum dependent particle concentrations one has to do on the client side:
21
22 TFile *f = TFile::Open("PriorProb/PriorProbabilities.root ");
23 TF1 *fitElectrons = (TF1 *)f->Get("fitElectrons");
24 TF1 *fitMuons = (TF1 *)f->Get("fitMuons");
25 TF1 *fitPions = (TF1 *)f->Get("fitPions");
26 TF1 *fitKaons = (TF1 *)f->Get("fitKaons");
27 TF1 *fitProtons = (TF1 *)f->Get("fitProtons");
28
29 If you want to pass these function in your task then you have to implement a function/setter that will do this. Then inside your task you do:
30
31 fElectronFunction->Eval(p);
32 fMuonFunction->Eval(p);
33 fPionFunction->Eval(p);
34 fKaonFunction->Eval(p);
35 fProtonFunction->Eval(p);
36
37 where p is the momentum of the track/particle.
38
39 A complete aexample of the implementation of the usage of these functions can be found in:
40 o) $ALICE_ROOT/PWG2/runProtonAnalysis.C
41 o) $ALICE_ROOT/PWG2/SPECTRA/AliProtonAnalysis.cxx (.h)