// Origin: Panos Christakoglou | Panos.Christakoglou@cern.ch
//-----------------------------------------------------------------
#include <Riostream.h>
-#include <TFile.h>
-#include <TSystem.h>
-#include <TF1.h>
-#include <TH2D.h>
-#include <TH1D.h>
-#include <TH1I.h>
-#include <TParticle.h>
#include <TCanvas.h>
#include <TLatex.h>
-
-#include "AliProtonAnalysisBase.h"
+#include <TF1.h>
+#include <TList.h>
+#include <TH1F.h>
#include <AliExternalTrackParam.h>
#include <AliESDEvent.h>
-#include <AliLog.h>
#include <AliPID.h>
-#include <AliESDVertex.h>
#include <AliVertexerTracks.h>
+#include <AliESDpid.h>
+#include <AliTPCPIDResponse.h>
+class AliLog;
+class AliESDVertex;
+
+#include "AliProtonAnalysisBase.h"
ClassImp(AliProtonAnalysisBase)
//____________________________________________________________________//
AliProtonAnalysisBase::AliProtonAnalysisBase() :
- TObject(), fProtonAnalysisLevel("ESD"),
- fTriggerMode(kMB2), fProtonAnalysisMode(kTPC), fProtonPIDMode(kBayesian),
+ TObject(), fProtonAnalysisLevel("ESD"), fAnalysisMC(kFALSE),
+ fTriggerMode(kMB2), kUseOnlineTrigger(kFALSE), kUseOfflineTrigger(kFALSE),
+ fPhysicsSelection(0),
+ fProtonAnalysisMode(kTPC), fProtonPIDMode(kBayesian),
fAnalysisEtaMode(kFALSE),
- fVxMax(100.), fVyMax(100.), fVzMax(100.),
+ fVxMax(100.), fVyMax(100.), fVzMax(100.), fMinNumOfContributors(0),
fNBinsX(0), fMinX(0), fMaxX(0),
fNBinsY(0), fMinY(0), fMaxY(0),
fMinTPCClusters(0), fMinITSClusters(0),
fMaxDCA3DFlag(kFALSE), fMaxDCA3DTPCFlag(kFALSE),
fMaxConstrainChi2Flag(kFALSE),
fITSRefitFlag(kFALSE), fTPCRefitFlag(kFALSE),
- fESDpidFlag(kFALSE), fTPCpidFlag(kFALSE),
+ fESDpidFlag(kFALSE), fTPCpidFlag(kFALSE), fTOFpidFlag(kFALSE),
fPointOnITSLayer1Flag(0), fPointOnITSLayer2Flag(0),
fPointOnITSLayer3Flag(0), fPointOnITSLayer4Flag(0),
fPointOnITSLayer5Flag(0), fPointOnITSLayer6Flag(0),
fNSigma(0),
fElectronFunction(0), fMuonFunction(0),
fPionFunction(0), fKaonFunction(0), fProtonFunction(0),
- fDebugMode(kFALSE) {
+ fDebugMode(kFALSE), fListVertexQA(new TList()) {
//Default constructor
for(Int_t i = 0; i < 5; i++) fPartFrac[i] = 0.0;
- for(Int_t i = 0; i < 5; i++) {
+ for(Int_t i = 0; i < 24; i++) {
fdEdxMean[i] = 0.0;
fdEdxSigma[i] = 0.0;
}
+ fListVertexQA->SetName("fListVertexQA");
+ TH1F *gHistVx = new TH1F("gHistVx",
+ "Vx distribution;V_{x} [cm];Entries",
+ 100,-5.,5.);
+ gHistVx->SetFillColor(kRed-2);
+ fListVertexQA->Add(gHistVx);
+ TH1F *gHistVxAccepted = new TH1F("gHistVxaccepted",
+ "Vx distribution;V_{x} [cm];Entries",
+ 100,-5.,5.);
+ fListVertexQA->Add(gHistVxAccepted);
+ TH1F *gHistVy = new TH1F("gHistVy",
+ "Vy distribution;V_{y} [cm];Entries",
+ 100,-5.,5.);
+ gHistVy->SetFillColor(kRed-2);
+ fListVertexQA->Add(gHistVy);
+ TH1F *gHistVyAccepted = new TH1F("gHistVyaccepted",
+ "Vy distribution;V_{y} [cm];Entries",
+ 100,-5.,5.);
+ fListVertexQA->Add(gHistVyAccepted);
+ TH1F *gHistVz = new TH1F("gHistVz",
+ "Vz distribution;V_{z} [cm];Entries",
+ 100,-25.,25.);
+ gHistVz->SetFillColor(kRed-2);
+ fListVertexQA->Add(gHistVz);
+ TH1F *gHistVzAccepted = new TH1F("gHistVzaccepted",
+ "Vz distribution;V_{z} [cm];Entries",
+ 100,-25.,25.);
+ fListVertexQA->Add(gHistVzAccepted);
+
+ TH1F *gHistNumberOfContributors = new TH1F("gHistNumberOfContributors",
+ "Number of contributors;N_{contr.};Entries",
+ 100,0.,100.);
+ fListVertexQA->Add(gHistNumberOfContributors);
+
+
}
//____________________________________________________________________//
if(fPionFunction) delete fPionFunction;
if(fKaonFunction) delete fKaonFunction;
if(fProtonFunction) delete fProtonFunction;
+ if(fListVertexQA) delete fListVertexQA;
}
//____________________________________________________________________//
//____________________________________________________________________//
Bool_t AliProtonAnalysisBase::IsInPhaseSpace(AliESDtrack* const track) {
// Checks if the track is outside the analyzed y-Pt phase space
- Double_t gPt = 0.0, gPx = 0.0, gPy = 0.0, gPz = 0.0;
+ Double_t gP = 0.0, gPt = 0.0, gPx = 0.0, gPy = 0.0, gPz = 0.0;
Double_t eta = 0.0;
- if(kTPC) {
+ if((fProtonAnalysisMode == kTPC) || (fProtonAnalysisMode == kHybrid)) {
AliExternalTrackParam *tpcTrack = (AliExternalTrackParam *)track->GetTPCInnerParam();
if(!tpcTrack) {
- gPt = 0.0; gPx = 0.0; gPy = 0.0; gPz = 0.0; eta = -10.0;
+ gP = 0.0; gPt = 0.0; gPx = 0.0; gPy = 0.0; gPz = 0.0; eta = -10.0;
}
else {
+ gP = tpcTrack->P();
gPt = tpcTrack->Pt();
gPx = tpcTrack->Px();
gPy = tpcTrack->Py();
gPz = tpcTrack->Pz();
eta = tpcTrack->Eta();
}
- }
+ }//standalone TPC or Hybrid TPC approaches
else {
+ gP = track->P();
gPt = track->Pt();
gPx = track->Px();
gPy = track->Py();
Printf("IsInPhaseSpace: Track rejected because it has a Pt value of %lf (accepted interval: %lf - %lf)",gPt,fMinY,fMaxY);
return kFALSE;
}
+ if((gP < fMinY) || (gP > fMaxY)) {
+ if(fDebugMode)
+ Printf("IsInPhaseSpace: Track rejected because it has a P value of %lf (accepted interval: %lf - %lf)",gP,fMinY,fMaxY);
+ return kFALSE;
+ }
if(fAnalysisEtaMode) {
if((eta < fMinX) || (eta > fMaxX)) {
if(fDebugMode)
}
//____________________________________________________________________//
-Bool_t AliProtonAnalysisBase::IsAccepted(AliESDEvent *esd,
- const AliESDVertex *vertex,
- AliESDtrack* track) {
+Bool_t AliProtonAnalysisBase::IsAccepted(AliESDtrack* track) {
// Checks if the track is excluded from the cuts
- Double_t gPt = 0.0, gPx = 0.0, gPy = 0.0, gPz = 0.0;
- Double_t dca[2] = {0.0,0.0}, cov[3] = {0.0,0.0,0.0}; //The impact parameters and their covariance.
- Double_t dca3D = 0.0;
-
- if((kTPC)&&(!kHybrid)) {
- AliExternalTrackParam *tpcTrack = (AliExternalTrackParam *)track->GetTPCInnerParam();
- if(!tpcTrack) {
- gPt = 0.0; gPx = 0.0; gPy = 0.0; gPz = 0.0;
- dca[0] = -100.; dca[1] = -100.; dca3D = -100.;
- cov[0] = -100.; cov[1] = -100.; cov[2] = -100.;
- }
- else {
- gPt = tpcTrack->Pt();
- gPx = tpcTrack->Px();
- gPy = tpcTrack->Py();
- gPz = tpcTrack->Pz();
- tpcTrack->PropagateToDCA(vertex,
- esd->GetMagneticField(),
- 100.,dca,cov);
- }
- }
- else if(kHybrid) {
- AliExternalTrackParam *tpcTrack = (AliExternalTrackParam *)track->GetTPCInnerParam();
- if(!tpcTrack) {
- gPt = 0.0; gPx = 0.0; gPy = 0.0; gPz = 0.0;
- dca[0] = -100.; dca[1] = -100.;
- cov[0] = -100.; cov[1] = -100.; cov[2] = -100.;
- }
- else {
- gPt = tpcTrack->Pt();
- gPx = tpcTrack->Px();
- gPy = tpcTrack->Py();
- gPz = tpcTrack->Pz();
- tpcTrack->PropagateToDCA(vertex,
- esd->GetMagneticField(),
- 100.,dca,cov);
- }
- }
- else{
- gPt = track->Pt();
- gPx = track->Px();
- gPy = track->Py();
- gPz = track->Pz();
- track->PropagateToDCA(vertex,
- esd->GetMagneticField(),
- 100.,dca,cov);
- }
- dca3D = TMath::Sqrt(TMath::Power(dca[0],2) +
- TMath::Power(dca[1],2));
-
Int_t fIdxInt[200];
Int_t nClustersITS = track->GetITSclusters(fIdxInt);
Int_t nClustersTPC = track->GetTPCclusters(fIdxInt);
return kFALSE;
}
}
- if(fMaxSigmaToVertexFlag) {
- if(GetSigmaToVertex(track) > fMaxSigmaToVertex) {
+ if(fMinTPCdEdxPointsFlag) {
+ if(track->GetTPCsignalN() < fMinTPCdEdxPoints) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has a %lf sigmas to vertex (max. requested: %lf)",GetSigmaToVertex(track),fMaxSigmaToVertex);
+ Printf("IsAccepted: Track rejected because it has %d TPC points for the calculation of the energy loss (min. requested: %d)",track->GetTPCsignalN(),fMinTPCdEdxPoints);
return kFALSE;
}
}
- if(fMaxSigmaToVertexTPCFlag) {
- if(GetSigmaToVertex(track) > fMaxSigmaToVertexTPC) {
+ if(fITSRefitFlag) {
+ if ((track->GetStatus() & AliESDtrack::kITSrefit) == 0) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has a %lf sigmas to vertex TPC (max. requested: %lf)",GetSigmaToVertex(track),fMaxSigmaToVertexTPC);
+ Printf("IsAccepted: Track rejected because it has no ITS refit flag");
return kFALSE;
}
}
- if(fMaxDCAXYFlag) {
- if(TMath::Abs(dca[0]) > fMaxDCAXY) {
+ if(fTPCRefitFlag) {
+ if ((track->GetStatus() & AliESDtrack::kTPCrefit) == 0) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has a value of dca(xy) of %lf (max. requested: %lf)",TMath::Abs(dca[0]),fMaxDCAXY);
+ Printf("IsAccepted: Track rejected because it has no TPC refit flag");
return kFALSE;
}
}
- if(fMaxDCAXYTPCFlag) {
- if(TMath::Abs(dca[0]) > fMaxDCAXYTPC) {
+ if(fESDpidFlag) {
+ if ((track->GetStatus() & AliESDtrack::kESDpid) == 0) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has a value of dca(xy) (TPC) of %lf (max. requested: %lf)",TMath::Abs(dca[0]),fMaxDCAXYTPC);
+ Printf("IsAccepted: Track rejected because it has no ESD pid flag");
return kFALSE;
}
}
- if(fMaxDCAZFlag) {
- if(TMath::Abs(dca[1]) > fMaxDCAZ) {
+ if(fTPCpidFlag) {
+ if ((track->GetStatus() & AliESDtrack::kTPCpid) == 0) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has a value of dca(z) of %lf (max. requested: %lf)",TMath::Abs(dca[1]),fMaxDCAZ);
+ Printf("IsAccepted: Track rejected because it has no TPC pid flag");
return kFALSE;
}
}
- if(fMaxDCAZTPCFlag) {
- if(TMath::Abs(dca[1]) > fMaxDCAZTPC) {
+ if(fTOFpidFlag) {
+ if ((track->GetStatus() & AliESDtrack::kTOFpid) == 0) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has a value of dca(z) (TPC) of %lf (max. requested: %lf)",TMath::Abs(dca[1]),fMaxDCAZTPC);
+ Printf("IsAccepted: Track rejected because it has no TOF pid flag");
return kFALSE;
}
}
- if(fMaxDCA3DFlag) {
- if(TMath::Abs(dca3D) > fMaxDCA3D) {
+
+ return kTRUE;
+}
+
+//____________________________________________________________________//
+Bool_t AliProtonAnalysisBase::IsPrimary(AliESDEvent *esd,
+ const AliESDVertex *vertex,
+ AliESDtrack* track) {
+ // Checks if the track is a primary-like candidate
+ Double_t gPt = 0.0, gPx = 0.0, gPy = 0.0, gPz = 0.0;
+ Double_t dca[2] = {0.0,0.0}, cov[3] = {0.0,0.0,0.0}; //The impact parameters and their covariance.
+ Double_t dca3D = 0.0;
+
+ if((fProtonAnalysisMode == kTPC)||(fProtonAnalysisMode == kHybrid)) {
+ AliExternalTrackParam *tpcTrack = (AliExternalTrackParam *)track->GetTPCInnerParam();
+ if(!tpcTrack) {
+ gPt = 0.0; gPx = 0.0; gPy = 0.0; gPz = 0.0;
+ dca[0] = -100.; dca[1] = -100.; dca3D = -100.;
+ cov[0] = -100.; cov[1] = -100.; cov[2] = -100.;
+ }
+ else {
+ gPt = tpcTrack->Pt();
+ gPx = tpcTrack->Px();
+ gPy = tpcTrack->Py();
+ gPz = tpcTrack->Pz();
+ tpcTrack->PropagateToDCA(vertex,
+ esd->GetMagneticField(),
+ 100.,dca,cov);
+ }
+ }//standalone TPC or hybrid TPC approaches
+ else {
+ gPt = track->Pt();
+ gPx = track->Px();
+ gPy = track->Py();
+ gPz = track->Pz();
+ track->PropagateToDCA(vertex,
+ esd->GetMagneticField(),
+ 100.,dca,cov);
+ }
+ dca3D = TMath::Sqrt(TMath::Power(dca[0],2) +
+ TMath::Power(dca[1],2));
+
+ if(fMaxSigmaToVertexFlag) {
+ if(GetSigmaToVertex(track) > fMaxSigmaToVertex) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has a value of dca(3D) of %lf (max. requested: %lf)",TMath::Abs(dca3D),fMaxDCA3D);
+ Printf("IsPrimary: Track rejected because it has a %lf sigmas to vertex (max. requested: %lf)",GetSigmaToVertex(track),fMaxSigmaToVertex);
return kFALSE;
}
}
- if(fMaxDCA3DTPCFlag) {
- if(TMath::Abs(dca3D) > fMaxDCA3DTPC) {
+ if(fMaxSigmaToVertexTPCFlag) {
+ if(GetSigmaToVertex(track) > fMaxSigmaToVertexTPC) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has a value of dca(3D) (TPC) of %lf (max. requested: %lf)",TMath::Abs(dca3D),fMaxDCA3DTPC);
+ Printf("IsPrimary: Track rejected because it has a %lf sigmas to vertex TPC (max. requested: %lf)",GetSigmaToVertex(track),fMaxSigmaToVertexTPC);
return kFALSE;
}
}
- if(fMaxConstrainChi2Flag) {
- if(track->GetConstrainedChi2() > 0)
- if(TMath::Log(track->GetConstrainedChi2()) > fMaxConstrainChi2) {
+ if(fMaxDCAXYFlag) {
+ if(TMath::Abs(dca[0]) > fMaxDCAXY) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has a value of the constrained chi2 to the vertex of %lf (max. requested: %lf)",TMath::Log(track->GetConstrainedChi2()),fMaxConstrainChi2);
+ Printf("IsPrimary: Track rejected because it has a value of dca(xy) of %lf (max. requested: %lf)",TMath::Abs(dca[0]),fMaxDCAXY);
return kFALSE;
- }
+ }
}
- if(fMinTPCdEdxPointsFlag) {
- if(track->GetTPCsignalN() < fMinTPCdEdxPoints) {
+ if(fMaxDCAXYTPCFlag) {
+ if(TMath::Abs(dca[0]) > fMaxDCAXYTPC) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has %d TPC points for the calculation of the energy loss (min. requested: %d)",track->GetTPCsignalN(),fMinTPCdEdxPoints);
+ Printf("IsPrimary: Track rejected because it has a value of dca(xy) (TPC) of %lf (max. requested: %lf)",TMath::Abs(dca[0]),fMaxDCAXYTPC);
return kFALSE;
}
}
- if(fITSRefitFlag) {
- if ((track->GetStatus() & AliESDtrack::kITSrefit) == 0) {
+ if(fMaxDCAZFlag) {
+ if(TMath::Abs(dca[1]) > fMaxDCAZ) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has no ITS refit flag");
+ Printf("IsPrimary: Track rejected because it has a value of dca(z) of %lf (max. requested: %lf)",TMath::Abs(dca[1]),fMaxDCAZ);
return kFALSE;
}
}
- if(fTPCRefitFlag) {
- if ((track->GetStatus() & AliESDtrack::kTPCrefit) == 0) {
+ if(fMaxDCAZTPCFlag) {
+ if(TMath::Abs(dca[1]) > fMaxDCAZTPC) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has no TPC refit flag");
+ Printf("IsPrimary: Track rejected because it has a value of dca(z) (TPC) of %lf (max. requested: %lf)",TMath::Abs(dca[1]),fMaxDCAZTPC);
return kFALSE;
}
}
- if(fESDpidFlag) {
- if ((track->GetStatus() & AliESDtrack::kESDpid) == 0) {
+ if(fMaxDCA3DFlag) {
+ if(TMath::Abs(dca3D) > fMaxDCA3D) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has no ESD pid flag");
+ Printf("IsPrimary: Track rejected because it has a value of dca(3D) of %lf (max. requested: %lf)",TMath::Abs(dca3D),fMaxDCA3D);
return kFALSE;
}
}
- if(fTPCpidFlag) {
- if ((track->GetStatus() & AliESDtrack::kTPCpid) == 0) {
+ if(fMaxDCA3DTPCFlag) {
+ if(TMath::Abs(dca3D) > fMaxDCA3DTPC) {
if(fDebugMode)
- Printf("IsAccepted: Track rejected because it has no TPC pid flag");
+ Printf("IsPrimary: Track rejected because it has a value of dca(3D) (TPC) of %lf (max. requested: %lf)",TMath::Abs(dca3D),fMaxDCA3DTPC);
return kFALSE;
}
}
+ if(fMaxConstrainChi2Flag) {
+ if(track->GetConstrainedChi2() > 0)
+ if(TMath::Log(track->GetConstrainedChi2()) > fMaxConstrainChi2) {
+ if(fDebugMode)
+ Printf("IsPrimary: Track rejected because it has a value of the constrained chi2 to the vertex of %lf (max. requested: %lf)",TMath::Log(track->GetConstrainedChi2()),fMaxConstrainChi2);
+ return kFALSE;
+ }
+ }
return kTRUE;
}
//____________________________________________________________________//
Float_t AliProtonAnalysisBase::GetSigmaToVertex(AliESDtrack* esdTrack) const {
// Calculates the number of sigma to the vertex.
-
Float_t b[2];
Float_t bRes[2];
Float_t bCov[3];
- if((kTPC)&&(!kHybrid))
+ if((fProtonAnalysisMode == kTPC)&&(fProtonAnalysisMode != kHybrid))
esdTrack->GetImpactParametersTPC(b,bCov);
else
esdTrack->GetImpactParameters(b,bCov);
Printf("GetVertex: Event rejected because the value of the vertex resolution in z is 0");
return 0;
}
-
+ ((TH1F *)(fListVertexQA->At(0)))->Fill(vertex->GetXv());
+ ((TH1F *)(fListVertexQA->At(2)))->Fill(vertex->GetYv());
+ ((TH1F *)(fListVertexQA->At(4)))->Fill(vertex->GetZv());
+
//check position
if(TMath::Abs(vertex->GetXv()) > gVxMax) {
if(fDebugMode)
Printf("GetVertex: Event rejected because it has a Vz value of %lf cm (accepted interval: -%lf - %lf)",TMath::Abs(vertex->GetZv()),gVzMax,gVzMax);
return 0;
}
+ ((TH1F *)(fListVertexQA->At(1)))->Fill(vertex->GetXv());
+ ((TH1F *)(fListVertexQA->At(3)))->Fill(vertex->GetYv());
+ ((TH1F *)(fListVertexQA->At(5)))->Fill(vertex->GetZv());
+ ((TH1F *)(fListVertexQA->At(6)))->Fill(vertex->GetNContributors());
+
+ //check number of contributors
+ if(fMinNumOfContributors > 0) {
+ if(fMinNumOfContributors > vertex->GetNContributors()) {
+ if(fDebugMode)
+ Printf("GetVertex: Event rejected because it has %d number of contributors (requested minimum: %d)",vertex->GetNContributors(),fMinNumOfContributors);
+
+ return 0;
+ }
+ }
return vertex;
}
TriggerMode trigger) {
// check if the event was triggered
ULong64_t triggerMask = esd->GetTriggerMask();
-
- // definitions from p-p.cfg
- ULong64_t spdFO = (1 << 14);
- ULong64_t v0left = (1 << 11);
- ULong64_t v0right = (1 << 12);
-
- switch (trigger) {
- case kMB1: {
- if (triggerMask & spdFO || ((triggerMask & v0left) || (triggerMask & v0right)))
- return kTRUE;
- break;
- }
- case kMB2: {
- if (triggerMask & spdFO && ((triggerMask & v0left) || (triggerMask & v0right)))
- return kTRUE;
- break;
+ TString firedTriggerClass = esd->GetFiredTriggerClasses();
+
+ if(fAnalysisMC) {
+ // definitions from p-p.cfg
+ ULong64_t spdFO = (1 << 14);
+ ULong64_t v0left = (1 << 11);
+ ULong64_t v0right = (1 << 12);
+
+ switch (trigger) {
+ case kMB1: {
+ if (triggerMask & spdFO || ((triggerMask & v0left) || (triggerMask & v0right)))
+ return kTRUE;
+ break;
+ }
+ case kMB2: {
+ if (triggerMask & spdFO && ((triggerMask & v0left) || (triggerMask & v0right)))
+ return kTRUE;
+ break;
+ }
+ case kSPDFASTOR: {
+ if (triggerMask & spdFO)
+ return kTRUE;
+ break;
+ }
+ }//switch
}
- case kSPDFASTOR: {
- if (triggerMask & spdFO)
+ else {
+ if(kUseOnlineTrigger) {
+ if(firedTriggerClass.Contains("CINT1B-ABCE-NOPF-ALL"))
+ return kTRUE;
+ }
+ else if(!kUseOnlineTrigger)
return kTRUE;
- break;
}
- }//switch
return kFALSE;
}
if(fMaxChi2PerTPCClusterFlag) listOfCuts += fMaxChi2PerTPCCluster;
else listOfCuts += "Not used";
l.DrawLatex(0.1,0.5,listOfCuts.Data());
+ listOfCuts = "Minimum number of TPC points for the dE/dx: ";
+ if(fMinTPCdEdxPointsFlag) listOfCuts += fMinTPCdEdxPoints;
+ else listOfCuts += "Not used";
+ l.DrawLatex(0.1,0.3,listOfCuts.Data());
c->cd(4)->SetFillColor(4);
l.DrawLatex(0.3,0.9,"Tracking related cuts");
}
//Definition of an N-sigma area around the dE/dx vs P band
else if(fProtonPIDMode == kSigma1) {
- AliExternalTrackParam *tpcTrack = (AliExternalTrackParam *)track->GetTPCInnerParam();
- if(tpcTrack) {
- gPt = tpcTrack->Pt();
- gP = tpcTrack->P();
- gEta = tpcTrack->Eta();
+ Double_t fAlephParameters[5];
+ if(fAnalysisMC) {
+ fAlephParameters[0] = 2.15898e+00/50.;
+ fAlephParameters[1] = 1.75295e+01;
+ fAlephParameters[2] = 3.40030e-09;
+ fAlephParameters[3] = 1.96178e+00;
+ fAlephParameters[4] = 3.91720e+00;
}
- //We start the P slices at >0.3GeV/c with a bining of 50MeV/c ==> Int_t(0.3001/0.05) = 6
- Int_t nbinP = Int_t(gP/0.05) - 6;
- Double_t tpcSignal = track->GetTPCsignal();
- Double_t dEdxMean = fdEdxMean[nbinP];
- Double_t dEdxSigma = fdEdxSigma[nbinP];
- if((tpcSignal <= dEdxMean + fNSigma*dEdxSigma)&&
- (tpcSignal <= dEdxMean + fNSigma*dEdxSigma))
+ else {
+ fAlephParameters[0] = 0.0283086;
+ fAlephParameters[1] = 2.63394e+01;
+ fAlephParameters[2] = 5.04114e-11;
+ fAlephParameters[3] = 2.12543e+00;
+ fAlephParameters[4] = 4.88663e+00;
+ }
+
+ Double_t nsigma = 100.0;
+ AliESDpid *fESDpid = new AliESDpid();
+ fESDpid->GetTPCResponse().SetBetheBlochParameters(fAlephParameters[0],fAlephParameters[1],fAlephParameters[2],fAlephParameters[3],fAlephParameters[4]);
+
+ AliExternalTrackParam *tpcTrack = (AliExternalTrackParam *)track->GetTPCInnerParam();
+ if(tpcTrack)
+ nsigma = TMath::Abs(fESDpid->NumberOfSigmasTPC(track,AliPID::kProton));
+
+ if(nsigma <= fNSigma)
return kTRUE;
}//kSigma1 PID method
//Another definition of an N-sigma area around the dE/dx vs P band
gP = tpcTrack->P();
gEta = tpcTrack->Eta();
}
- //We start the P slices at >0.3GeV/c with a bining of 50MeV/c ==> Int_t(0.3001/0.05) = 6
- Int_t nbinP = Int_t(gP/0.05) - 6;
- Double_t tpcSignal = track->GetTPCsignal();
- Double_t dEdxTheory = Bethe(gP/9.38270000000000048e-01);
- Double_t dEdxSigma = fdEdxSigma[nbinP];
- Double_t nsigma = TMath::Abs(tpcSignal - dEdxTheory)/(tpcSignal*(dEdxSigma/TMath::Sqrt(track->GetTPCsignalN())));
- if(nsigma <= fNSigma)
+ Double_t fAlephParameters[5];
+ if(fAnalysisMC) {
+ fAlephParameters[0] = 2.15898e+00/50.;
+ fAlephParameters[1] = 1.75295e+01;
+ fAlephParameters[2] = 3.40030e-09;
+ fAlephParameters[3] = 1.96178e+00;
+ fAlephParameters[4] = 3.91720e+00;
+ }
+ else {
+ fAlephParameters[0] = 0.0283086;
+ fAlephParameters[1] = 2.63394e+01;
+ fAlephParameters[2] = 5.04114e-11;
+ fAlephParameters[3] = 2.12543e+00;
+ fAlephParameters[4] = 4.88663e+00;
+ }
+
+ AliESDpid *fESDpid = new AliESDpid();
+ AliTPCPIDResponse tpcResponse = fESDpid->GetTPCResponse();
+ tpcResponse.SetBetheBlochParameters(fAlephParameters[0],fAlephParameters[1],fAlephParameters[2],fAlephParameters[3],fAlephParameters[4]);
+ Double_t normalizeddEdx = TMath::Log(track->GetTPCsignal()/tpcResponse.GetExpectedSignal(gP,AliPID::kProton));
+
+ if(normalizeddEdx >= -0.15)
return kTRUE;
}
}
//________________________________________________________________________
-Double_t AliProtonAnalysisBase::Bethe(Double_t bg) {
+Double_t AliProtonAnalysisBase::Bethe(Double_t bg) const {
// This is the Bethe-Bloch function normalised to 1 at the minimum
// We renormalize it based on the MC information
// WARNING: To be revised soon!!!