i)Adding the file containing the 5 TF1 objects to calculate the momentum dependent...

diff --git a/PWG2/data/PriorProbabilities.root b/PWG2/data/PriorProbabilities.root
new file mode 100644 (file)
index 0000000..0a0e880
Binary files /dev/null and b/PWG2/data/PriorProbabilities.root differ

diff --git a/PWG2/data/README b/PWG2/data/README
new file mode 100644 (file)
index 0000000..169b9d2
--- /dev/null
+++ b/PWG2/data/README
@@ -0,0 +1,41 @@
+===========================================================
+=          Author: Panos.Christakoglou@cern.ch            =
+===========================================================
+
+
+PriorPropabilities.root
+_______________________
+
+The file contains 5 TF1 objects namely: 
+fitElectrons
+fitMuons
+fitPions
+fitKaons
+fitProtons
+
+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).
+
+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.
+
+To get the momentum dependent particle concentrations one has to do on the client side:
+
+TFile *f = TFile::Open("PriorProb/PriorProbabilities.root ");
+TF1 *fitElectrons = (TF1 *)f->Get("fitElectrons");
+TF1 *fitMuons = (TF1 *)f->Get("fitMuons");
+TF1 *fitPions = (TF1 *)f->Get("fitPions");
+TF1 *fitKaons = (TF1 *)f->Get("fitKaons");
+TF1 *fitProtons = (TF1 *)f->Get("fitProtons");
+
+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:
+
+fElectronFunction->Eval(p);
+fMuonFunction->Eval(p);
+fPionFunction->Eval(p);
+fKaonFunction->Eval(p);
+fProtonFunction->Eval(p);
+
+where p is the momentum of the track/particle.
+
+A complete aexample of the implementation of the usage of these functions can be found in:
+o) $ALICE_ROOT/PWG2/runProtonAnalysis.C
+o) $ALICE_ROOT/PWG2/SPECTRA/AliProtonAnalysis.cxx (.h)
\ No newline at end of file