Adding the AliAnalysisGUI class which is the main class that controls the GUI.
[u/mrichter/AliRoot.git] / ITS / AliITSpidESD1.cxx
1 /**************************************************************************
2  * Copyright(c) 2005-2007, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 //-----------------------------------------------------------------
17 // ITS PID method # 1
18 //           Implementation of the ITS PID class
19 // Very naive one... Should be made better by the detector experts...
20 //      Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
21 //-----------------------------------------------------------------
22 #include "AliITSpidESD.h"
23 #include "AliITSpidESD1.h"
24 #include "AliESD.h"
25 #include "AliESDtrack.h"
26
27 ClassImp(AliITSpidESD1)
28
29 AliITSpidESD1::AliITSpidESD1(): AliITSpidESD(),
30 fMIP(0),
31 fRes(0),
32 fRange(0) 
33 {
34   //Default constructor
35 }
36 //_________________________________________________________________________
37 AliITSpidESD1::AliITSpidESD1(Double_t *param): AliITSpidESD(),
38 fMIP(param[0]),
39 fRes(param[1]),
40 fRange(param[2])
41 {
42   //
43   //  The main constructor
44   //
45 }
46
47
48 //_________________________________________________________________________
49 Int_t AliITSpidESD1::MakePID(AliESD *event)
50 {
51   //
52   //  This function calculates the "detector response" PID probabilities 
53   //
54   Int_t ntrk=event->GetNumberOfTracks();
55   for (Int_t i=0; i<ntrk; i++) {
56     AliESDtrack *t=event->GetTrack(i);
57     if ((t->GetStatus()&AliESDtrack::kITSin )==0)
58       if ((t->GetStatus()&AliESDtrack::kITSout)==0) continue;
59     Double_t mom=t->GetP();
60     Double_t dedx=t->GetITSsignal()/fMIP;
61     Int_t ns=AliPID::kSPECIES;
62     Double_t p[10];
63     for (Int_t j=0; j<ns; j++) {
64       Double_t mass=AliPID::ParticleMass(j);
65       Double_t bethe=AliITSpidESD::Bethe(mom/mass); 
66       Double_t sigma=fRes*bethe;
67       if (TMath::Abs(dedx-bethe) > fRange*sigma) {
68         p[j]=TMath::Exp(-0.5*fRange*fRange)/sigma;
69         continue;
70       }
71       p[j]=TMath::Exp(-0.5*(dedx-bethe)*(dedx-bethe)/(sigma*sigma))/sigma;
72     }
73     t->SetITSpid(p);
74   }
75   return 0;
76 }