* provided "as is" without express or implied warranty. *
**************************************************************************/
+/* $Id$ */
+
//-----------------------------------------------------------------------//
// ITS PID class --- method # 2 //
// //
//-----------------------------------------------------------------------//
#include "AliITSpidESD2.h"
-#include "AliESD.h"
+#include "AliESDEvent.h"
#include "AliESDtrack.h"
#include "AliITStrackV2.h"
#include "AliITSRecPoint.h"
//_________________________________________________________________________
AliITSpidESD2::AliITSpidESD2():AliITSpidESD(),
fTracker(0),
-fLoader(0),
fSp(0)
{ //
// The main constructor
}
//_________________________________________________________________________
-AliITSpidESD2::AliITSpidESD2(AliITStrackerMI* tracker,AliITSLoader* loader):AliITSpidESD(),
-fTracker(0),
-fLoader(0),
+AliITSpidESD2::AliITSpidESD2(AliITStrackerMI* tracker):AliITSpidESD(),
+fTracker(tracker),
fSp(0)
{ //
// The main constructor
- fTracker=tracker;
- fLoader=loader;
fSp=new AliITSSteerPid();
fSp->InitLayer();
}
//______________________________________________________________________
AliITSpidESD2::AliITSpidESD2(const AliITSpidESD2 &ob) :AliITSpidESD(ob),
fTracker(ob.fTracker),
-fLoader(ob.fLoader),
fSp(ob.fSp)
{
// Copy constructor
}
//_________________________________________________________________________
-Int_t AliITSpidESD2::MakePID(AliESD *event)
+Int_t AliITSpidESD2::MakePID(TTree *clustersTree, AliESDEvent *event)
{
//
//
Double_t xr,par[5];
AliITStrackV2* track=0;
- fLoader->LoadRecPoints();
- TTree *cTree=fLoader->TreeR();
- fTracker->LoadClusters(cTree);
+ fTracker->LoadClusters(clustersTree);
printf("==== Landau Fit PID ITS ====== \n");
Int_t ntrk=event->GetNumberOfTracks();
Double_t momits;
}
AliITSRecPoint* cluarr[12];
Float_t qclu[8],qclucorr[8],nlay[8];
- for(Int_t i=0;i<8;i++){
- qclu[i]=0;
- qclucorr[i]=0;
- nlay[i]=0;
+ for(Int_t ii=0;ii<8;ii++){
+ qclu[ii]=0;
+ qclucorr[ii]=0;
+ nlay[ii]=0;
}
Int_t jj=0;
- for(Int_t i=0;i<12;i++){
- cluind[i]=track->GetClusterIndex(i);
- if(cluind[i]>0){
- cluarr[i]=(AliITSRecPoint*)fTracker->GetCluster(cluind[i]);
- Int_t lay=cluarr[i]->GetLayer();
+ for(Int_t ij=0;ij<12;ij++){
+ cluind[ij]=track->GetClusterIndex(ij);
+ if(cluind[ij]>0){
+ cluarr[ij]=(AliITSRecPoint*)fTracker->GetCluster(cluind[ij]);
+ Int_t lay=cluarr[ij]->GetLayer();
if(lay>1){//sdd+ssd only
- qclu[jj]=cluarr[i]->GetQ();
+ qclu[jj]=cluarr[ij]->GetQ();
qclucorr[jj]=qclu[jj]*TMath::Sqrt((1-snp*snp)/(1+tgl*tgl));
nlay[jj]=lay;
jj++;
esdtr->SetITSpid(condprobfun);
delete track;
- }
+ }
fTracker->UnloadClusters();
- fLoader->UnloadRecPoints();
return 0;
}