#include "AliComparisonSelector.h"
+ClassImp(AliComparisonSelector)
AliComparisonSelector::AliComparisonSelector(TTree *) :
-AliComparisonDraw::AliComparisonDraw():
- TObject(),
- fPtResolLPT(0),
- fPtResolHPT(0)
-{
- InitHisto();
-}
-
-void AliComparisonDraw::InitHisto(){
- //
- //
- //
- fPtResolLPT = new TH2F("Pt resol","pt resol",10, 0.1,3,200,-0.2,0.2);
- fPtResolHPT = new TH2F("Pt resol","pt resol",10, 2,100,200,-0.3,0.3);
- //
- fPtPoolLPT = new TH2F("Pt pool","pt pool",10, 0.1,3,200,-6,6);
- fPtPoolHPT = new TH2F("Pt pool","pt pool",10, 2,100,200,-6,6);
-}
-
-void AliComparisonDraw::Process(AliMCInfo* infoMC, AliESDRecInfo *infoRC){
- //
- //
- //
- Float_t mcpt = infoMC->GetParticle().Pt();
-
- //
- //
- if (infoRC->GetStatus(1)==0) return;
- if (!infoRC->GetESDtrack()) return; //buggy line
-
-
- if (infoRC->GetESDtrack()->GetTPCNcls()<10) return;
-
- // printf("Pt\t%f\t%f\n",mcpt, infoRC->GetESDtrack()->Pt());
-
- Float_t deltaPt= (mcpt-infoRC->GetESDtrack()->Pt())/mcpt;
- Float_t poolPt= (1/mcpt-infoRC->GetESDtrack()->OneOverPt())/
- TMath::Sqrt(infoRC->GetESDtrack()->GetSigma1Pt2());
- fPtResolLPT->Fill(mcpt,deltaPt);
- fPtResolHPT->Fill(mcpt,deltaPt);
- fPtPoolLPT->Fill(mcpt,poolPt);
- fPtPoolHPT->Fill(mcpt,poolPt);
-}