Adding the normalized dE/dx vs P plot

diff --git a/PWG2/SPECTRA/AliProtonAnalysis.cxx b/PWG2/SPECTRA/AliProtonAnalysis.cxx
index 38a98509747f8355db3bef341893137aee4d2eb3..9867a5c3605ed428995179cb43c3c215d5c5a779 100644 (file)
--- a/PWG2/SPECTRA/AliProtonAnalysis.cxx
+++ b/PWG2/SPECTRA/AliProtonAnalysis.cxx
@@ -38,7 +38,8 @@
 #include <AliCFContainer.h>
 #include <AliCFEffGrid.h>
 #include <AliCFDataGrid.h>
-//#include <AliESDVertex.h>
+#include <AliTPCPIDResponse.h>
+#include <AliESDpid.h>
 class AliLog;
 class AliESDVertex;
 
@@ -690,6 +691,28 @@ void AliProtonAnalysis::Analyze(AliESDEvent* esd,
   Int_t nTracks = 0;
   Int_t nIdentifiedProtons = 0, nIdentifiedAntiProtons = 0;
   Int_t nSurvivedProtons = 0, nSurvivedAntiProtons = 0;
+  
+  //=========================================//
+  //Aleph parametrization
+  Double_t fAlephParameters[5];
+  if(fProtonAnalysisBase->GetMCAnalysisMode()) {
+    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]);
+  //=========================================//  
 
   fHistEvents->Fill(1); //number of analyzed events
   Double_t containerInput[2] ;
@@ -729,16 +752,17 @@ void AliProtonAnalysis::Analyze(AliESDEvent* esd,
       if(fProtonAnalysisBase->IsPrimary(esd,vertex,track)) {
        if(fProtonAnalysisBase->IsAccepted(track)) {
          ((TH2F *)(fQA2DList->At(0)))->Fill(gP,track->GetTPCsignal());
-         ((TH3F *)(fQA2DList->At(2)))->Fill(tpcTrack->Eta(),
+         ((TH2F *)(fQA2DList->At(2)))->Fill(gP,TMath::Log(track->GetTPCsignal()/tpcResponse.GetExpectedSignal(gP,AliPID::kProton)));
+         ((TH3F *)(fQA2DList->At(4)))->Fill(tpcTrack->Eta(),
                                             tpcTrack->Phi()*180./TMath::Pi(),
                                             npointsTPCdEdx);
-         ((TH3F *)(fQA2DList->At(4)))->Fill(tpcTrack->Eta(),
+         ((TH3F *)(fQA2DList->At(6)))->Fill(tpcTrack->Eta(),
                                             tpcTrack->Phi()*180./TMath::Pi(),
                                             nClustersTPC);
-         ((TH3F *)(fQA2DList->At(6)))->Fill(gPt,
+         ((TH3F *)(fQA2DList->At(8)))->Fill(gPt,
                                             tpcTrack->Phi()*180./TMath::Pi(),
                                             npointsTPCdEdx);
-         ((TH3F *)(fQA2DList->At(8)))->Fill(gPt,
+         ((TH3F *)(fQA2DList->At(10)))->Fill(gPt,
                                             tpcTrack->Phi()*180./TMath::Pi(),
                                             nClustersTPC);     
        }//quality cuts
@@ -774,35 +798,36 @@ void AliProtonAnalysis::Analyze(AliESDEvent* esd,
          //Step: kStepSurvived
          if(fProtonAnalysisBase->IsAccepted(track)) {
            ((TH2F *)(fQA2DList->At(1)))->Fill(gP,track->GetTPCsignal());
-           ((TH3F *)(fQA2DList->At(3)))->Fill(tpcTrack->Eta(),
+           ((TH2F *)(fQA2DList->At(3)))->Fill(gP,TMath::Log(track->GetTPCsignal()/tpcResponse.GetExpectedSignal(gP,AliPID::kProton)));
+           ((TH3F *)(fQA2DList->At(5)))->Fill(tpcTrack->Eta(),
                                               tpcTrack->Phi()*180./TMath::Pi(),
                                               npointsTPCdEdx);
-           ((TH3F *)(fQA2DList->At(5)))->Fill(tpcTrack->Eta(),
+           ((TH3F *)(fQA2DList->At(7)))->Fill(tpcTrack->Eta(),
                                               tpcTrack->Phi()*180./TMath::Pi(),
                                               nClustersTPC);
-           ((TH3F *)(fQA2DList->At(7)))->Fill(gPt,
+           ((TH3F *)(fQA2DList->At(9)))->Fill(gPt,
                                               tpcTrack->Phi()*180./TMath::Pi(),
                                               npointsTPCdEdx);
-           ((TH3F *)(fQA2DList->At(9)))->Fill(gPt,
+           ((TH3F *)(fQA2DList->At(11)))->Fill(gPt,
                                               tpcTrack->Phi()*180./TMath::Pi(),
                                               nClustersTPC);
            
            if(tpcTrack->Charge() > 0) {
              fProtonContainer->Fill(containerInput,kStepSurvived);   
-             ((TH2F *)(fQA2DList->At(10)))->Fill(tpcTrack->Eta(),
+             ((TH2F *)(fQA2DList->At(12)))->Fill(tpcTrack->Eta(),
                                                  tpcTrack->Phi()*180./TMath::Pi());
-             ((TH2F *)(fQA2DList->At(12)))->Fill(tpcTrack->Pt(),
+             ((TH2F *)(fQA2DList->At(14)))->Fill(tpcTrack->Pt(),
                                                  TMath::Abs(dca[0]));
-             ((TH2F *)(fQA2DList->At(13)))->Fill(tpcTrack->Pt(),
+             ((TH2F *)(fQA2DList->At(15)))->Fill(tpcTrack->Pt(),
                                                  TMath::Abs(dca[1]));
            }//protons
            else if(tpcTrack->Charge() < 0) {
              fAntiProtonContainer->Fill(containerInput,kStepSurvived);   
-             ((TH2F *)(fQA2DList->At(11)))->Fill(tpcTrack->Eta(),
+             ((TH2F *)(fQA2DList->At(13)))->Fill(tpcTrack->Eta(),
                                                  tpcTrack->Phi()*180./TMath::Pi());
-             ((TH2F *)(fQA2DList->At(14)))->Fill(tpcTrack->Pt(),
+             ((TH2F *)(fQA2DList->At(16)))->Fill(tpcTrack->Pt(),
                                                  TMath::Abs(dca[0]));
-             ((TH2F *)(fQA2DList->At(15)))->Fill(tpcTrack->Pt(),
+             ((TH2F *)(fQA2DList->At(17)))->Fill(tpcTrack->Pt(),
                                                  TMath::Abs(dca[1]));
            }//antiprotons
            
@@ -835,16 +860,17 @@ void AliProtonAnalysis::Analyze(AliESDEvent* esd,
       if(fProtonAnalysisBase->IsPrimary(esd,vertex,track)) {
        if(fProtonAnalysisBase->IsAccepted(track)) {
          ((TH2F *)(fQA2DList->At(0)))->Fill(gP,track->GetTPCsignal());
-         ((TH3F *)(fQA2DList->At(2)))->Fill(track->Eta(),
+         ((TH2F *)(fQA2DList->At(2)))->Fill(gP,TMath::Log(track->GetTPCsignal()/tpcResponse.GetExpectedSignal(gP,AliPID::kProton)));
+         ((TH3F *)(fQA2DList->At(4)))->Fill(track->Eta(),
                                             track->Phi()*180./TMath::Pi(),
                                             npointsTPCdEdx);
-         ((TH3F *)(fQA2DList->At(4)))->Fill(track->Eta(),
+         ((TH3F *)(fQA2DList->At(6)))->Fill(track->Eta(),
                                             track->Phi()*180./TMath::Pi(),
                                             nClustersTPC);
-         ((TH3F *)(fQA2DList->At(6)))->Fill(gPt,
+         ((TH3F *)(fQA2DList->At(8)))->Fill(gPt,
                                             track->Phi()*180./TMath::Pi(),
                                             npointsTPCdEdx);
-         ((TH3F *)(fQA2DList->At(8)))->Fill(gPt,
+         ((TH3F *)(fQA2DList->At(10)))->Fill(gPt,
                                             track->Phi()*180./TMath::Pi(),
                                             nClustersTPC);     
        }//quality cuts
@@ -880,35 +906,36 @@ void AliProtonAnalysis::Analyze(AliESDEvent* esd,
          //Step: kStepSurvived
          if(fProtonAnalysisBase->IsAccepted(track)) {
            ((TH2F *)(fQA2DList->At(1)))->Fill(gP,track->GetTPCsignal());
-           ((TH3F *)(fQA2DList->At(3)))->Fill(track->Eta(),
+           ((TH2F *)(fQA2DList->At(3)))->Fill(gP,TMath::Log(track->GetTPCsignal()/tpcResponse.GetExpectedSignal(gP,AliPID::kProton)));
+           ((TH3F *)(fQA2DList->At(5)))->Fill(track->Eta(),
                                               track->Phi()*180./TMath::Pi(),
                                               npointsTPCdEdx);
-           ((TH3F *)(fQA2DList->At(5)))->Fill(track->Eta(),
+           ((TH3F *)(fQA2DList->At(7)))->Fill(track->Eta(),
                                               track->Phi()*180./TMath::Pi(),
                                               nClustersTPC);
-           ((TH3F *)(fQA2DList->At(7)))->Fill(gPt,
+           ((TH3F *)(fQA2DList->At(9)))->Fill(gPt,
                                               track->Phi()*180./TMath::Pi(),
                                               npointsTPCdEdx);
-           ((TH3F *)(fQA2DList->At(9)))->Fill(gPt,
+           ((TH3F *)(fQA2DList->At(11)))->Fill(gPt,
                                               track->Phi()*180./TMath::Pi(),
                                               nClustersTPC);
            
            if(track->Charge() > 0) {
              fProtonContainer->Fill(containerInput,kStepSurvived);   
-             ((TH2F *)(fQA2DList->At(10)))->Fill(track->Eta(),
+             ((TH2F *)(fQA2DList->At(12)))->Fill(track->Eta(),
                                                  track->Phi()*180./TMath::Pi());
-             ((TH2F *)(fQA2DList->At(12)))->Fill(track->Pt(),
+             ((TH2F *)(fQA2DList->At(14)))->Fill(track->Pt(),
                                                  TMath::Abs(dca[0]));
-             ((TH2F *)(fQA2DList->At(13)))->Fill(track->Pt(),
+             ((TH2F *)(fQA2DList->At(15)))->Fill(track->Pt(),
                                                  TMath::Abs(dca[1]));
            }//protons
            else if(track->Charge() < 0) {
              fAntiProtonContainer->Fill(containerInput,kStepSurvived);   
-             ((TH2F *)(fQA2DList->At(11)))->Fill(track->Eta(),
+             ((TH2F *)(fQA2DList->At(13)))->Fill(track->Eta(),
                                                  track->Phi()*180./TMath::Pi());
-             ((TH2F *)(fQA2DList->At(14)))->Fill(track->Pt(),
+             ((TH2F *)(fQA2DList->At(16)))->Fill(track->Pt(),
                                                  TMath::Abs(dca[0]));
-             ((TH2F *)(fQA2DList->At(15)))->Fill(track->Pt(),
+             ((TH2F *)(fQA2DList->At(17)))->Fill(track->Pt(),
                                                  TMath::Abs(dca[1]));
            }//antiprotons
            
@@ -1321,6 +1348,12 @@ void AliProtonAnalysis::InitQA() {
   TH2F *gHistProtonsdEdxP = new TH2F("gHistProtonsdEdxP","Accepted protons dE/dx (TPC); P [GeV/c]; dE/dx [a.u]",1000,0.05,20.05,600,0,600);
   fQA2DList->Add(gHistProtonsdEdxP);
 
+  //normalized dEdx plots
+  TH2F *gHistZP = new TH2F("gHistZP","Normalized dE/dx (TPC); P [GeV/c]; ln[(dE/dx)_{exp.}/(dE/dx)_{BB}] ",1000,0.05,20.05,100,-2.0,2.0);
+  fQA2DList->Add(gHistZP);
+  TH2F *gHistProtonsZP = new TH2F("gHistProtonsZP","Normalized dE/dx (TPC); P [GeV/c]; ln[(dE/dx)_{exp.}/(dE/dx)_{BB}] ",1000,0.05,20.05,100,-2.0,2.0);
+  fQA2DList->Add(gHistProtonsZP);
+
   //eta-phi-Npoints(dEdx)
   TH3F *gHistEtaPhiTPCdEdxNPoints = new TH3F("gHistEtaPhiTPCdEdxNPoints",
                                             ";#eta;#phi;N_{points}(TPC)",

diff --git a/PWG2/SPECTRA/AliProtonAnalysis.h b/PWG2/SPECTRA/AliProtonAnalysis.h
index 2e7c14d33b7ac7af3a3fc3a71e4b0e97cdff7e44..45bc75851428525b2a9c049dcc3cb2e638537221 100644 (file)
--- a/PWG2/SPECTRA/AliProtonAnalysis.h
+++ b/PWG2/SPECTRA/AliProtonAnalysis.h
@@ -75,12 +75,8 @@ class AliProtonAnalysis : public TObject {
   AliCFContainer *GetProtonContainer() const {return fProtonContainer;}
   AliCFContainer *GetAntiProtonContainer() const {return fAntiProtonContainer;}
 
-  //TH2D *GetProtonYPtHistogram() const {return fHistYPtProtons;}
-  //TH2D *GetAntiProtonYPtHistogram() const {return fHistYPtAntiProtons;}
-  TH2D *GetProtonYPtHistogram() const {
-    return fProtonContainer->ShowProjection(0,1,kStepInPhaseSpace);}
-  TH2D *GetAntiProtonYPtHistogram() const {
-    return fAntiProtonContainer->ShowProjection(0,1,kStepInPhaseSpace);}
+  TH2D *GetProtonYPtHistogram() const {return fHistYPtProtons;}
+  TH2D *GetAntiProtonYPtHistogram() const {return fHistYPtAntiProtons;}
   TH1D *GetProtonYHistogram();
   TH1D *GetAntiProtonYHistogram();
   TH1D *GetProtonPtHistogram();

diff --git a/PWG2/SPECTRA/macros/drawProtonResults.C b/PWG2/SPECTRA/macros/drawProtonResults.C
index 29401e9240a78d3953592afc14e3ba99041c1d00..563bbeab40144a21e78b7f4330ba00d260386d79 100644 (file)
--- a/PWG2/SPECTRA/macros/drawProtonResults.C
+++ b/PWG2/SPECTRA/macros/drawProtonResults.C
@@ -173,8 +173,15 @@ void drawQAPlots(const char* analysisOutput,
   grProtons->SetLineColor(4); grProtons->SetLineWidth(2);
   grProtons->SetName("grProtons");
 
+  //2D de/dx vs P
+  TH2F *gHistZP = dynamic_cast<TH2F *>(fQA2DList->At(2));
+  gHistZP->SetStats(kFALSE);
+  //drawdEdx(gHistdEdxP);
+  TH2F *gHistProtonsZP = dynamic_cast<TH2F *>(fQA2DList->At(3));
+  gHistProtonsZP->SetStats(kFALSE);
+
   //3D eta-phi-NPoints(dEdx)
-  TH3F *gHistEtaPhiTPCdEdxNPoints = dynamic_cast<TH3F *>(fQA2DList->At(2));
+  TH3F *gHistEtaPhiTPCdEdxNPoints = dynamic_cast<TH3F *>(fQA2DList->At(4));
   TH2D *gHistEtaPhi = dynamic_cast<TH2D *>gHistEtaPhiTPCdEdxNPoints->Project3D("yx");
   gHistEtaPhi->SetStats(kFALSE);
   TH2D *gHistEtaTPCdEdxNPoints = dynamic_cast<TH2D *>gHistEtaPhiTPCdEdxNPoints->Project3D("zx");
@@ -183,7 +190,7 @@ void drawQAPlots(const char* analysisOutput,
   gHistPhiTPCdEdxNPoints->SetStats(kFALSE);
 
   //3D eta-phi-NPoints(dEdx): protons
-  TH3F *gHistProtonsEtaPhiTPCdEdxNPoints = dynamic_cast<TH3F *>(fQA2DList->At(3));
+  TH3F *gHistProtonsEtaPhiTPCdEdxNPoints = dynamic_cast<TH3F *>(fQA2DList->At(5));
   TH2D *gHistProtonsEtaPhi = dynamic_cast<TH2D *>gHistProtonsEtaPhiTPCdEdxNPoints->Project3D("yx");
   gHistProtonsEtaPhi->SetStats(kFALSE);
   TH2D *gHistProtonsEtaTPCdEdxNPoints = dynamic_cast<TH2D *>gHistProtonsEtaPhiTPCdEdxNPoints->Project3D("zx");
@@ -192,7 +199,7 @@ void drawQAPlots(const char* analysisOutput,
   gHistProtonsPhiTPCdEdxNPoints->SetStats(kFALSE);
 
   //3D eta-phi-NPoints
-  TH3F *gHistEtaPhiTPCNPoints = dynamic_cast<TH3F *>(fQA2DList->At(4));
+  TH3F *gHistEtaPhiTPCNPoints = dynamic_cast<TH3F *>(fQA2DList->At(6));
   TH2D *gHistEtaPhi = dynamic_cast<TH2D *>gHistEtaPhiTPCNPoints->Project3D("yx");
   gHistEtaPhi->SetStats(kFALSE);
   TH2D *gHistEtaTPCNPoints = dynamic_cast<TH2D *>gHistEtaPhiTPCNPoints->Project3D("zx");
@@ -201,7 +208,7 @@ void drawQAPlots(const char* analysisOutput,
   gHistPhiTPCNPoints->SetStats(kFALSE);
 
   //3D eta-phi-NPoints: protons
-  TH3F *gHistProtonsEtaPhiTPCNPoints = dynamic_cast<TH3F *>(fQA2DList->At(5));
+  TH3F *gHistProtonsEtaPhiTPCNPoints = dynamic_cast<TH3F *>(fQA2DList->At(7));
   TH2D *gHistProtonsEtaPhi = dynamic_cast<TH2D *>gHistProtonsEtaPhiTPCNPoints->Project3D("yx");
   gHistProtonsEtaPhi->SetStats(kFALSE);
   TH2D *gHistProtonsEtaTPCNPoints = dynamic_cast<TH2D *>gHistProtonsEtaPhiTPCNPoints->Project3D("zx");
@@ -210,7 +217,7 @@ void drawQAPlots(const char* analysisOutput,
   gHistProtonsPhiTPCNPoints->SetStats(kFALSE);
 
   //3D pt-phi-NPoints(dEdx)
-  TH3F *gHistPtPhiTPCdEdxNPoints = dynamic_cast<TH3F *>(fQA2DList->At(6));
+  TH3F *gHistPtPhiTPCdEdxNPoints = dynamic_cast<TH3F *>(fQA2DList->At(8));
   TH2D *gHistPtPhi = dynamic_cast<TH2D *>gHistPtPhiTPCdEdxNPoints->Project3D("yx");
   gHistPtPhi->SetStats(kFALSE);
   TH2D *gHistPtTPCdEdxNPoints = dynamic_cast<TH2D *>gHistPtPhiTPCdEdxNPoints->Project3D("zx");
@@ -219,7 +226,7 @@ void drawQAPlots(const char* analysisOutput,
   gHistPhiTPCdEdxNPoints->SetStats(kFALSE);
 
   //3D pt-phi-NPoints(dEdx): protons
-  TH3F *gHistProtonsPtPhiTPCdEdxNPoints = dynamic_cast<TH3F *>(fQA2DList->At(7));
+  TH3F *gHistProtonsPtPhiTPCdEdxNPoints = dynamic_cast<TH3F *>(fQA2DList->At(9));
   TH2D *gHistProtonsPtPhi = dynamic_cast<TH2D *>gHistProtonsPtPhiTPCdEdxNPoints->Project3D("yx");
   gHistProtonsPtPhi->SetStats(kFALSE);
   TH2D *gHistProtonsPtTPCdEdxNPoints = dynamic_cast<TH2D *>gHistProtonsPtPhiTPCdEdxNPoints->Project3D("zx");
@@ -228,7 +235,7 @@ void drawQAPlots(const char* analysisOutput,
   gHistProtonsPhiTPCdEdxNPoints->SetStats(kFALSE);
 
   //3D pt-phi-NPoints
-  TH3F *gHistPtPhiTPCNPoints = dynamic_cast<TH3F *>(fQA2DList->At(8));
+  TH3F *gHistPtPhiTPCNPoints = dynamic_cast<TH3F *>(fQA2DList->At(10));
   TH2D *gHistPtPhi = dynamic_cast<TH2D *>gHistPtPhiTPCNPoints->Project3D("yx");
   gHistPtPhi->SetStats(kFALSE);
   TH2D *gHistPtTPCNPoints = dynamic_cast<TH2D *>gHistPtPhiTPCNPoints->Project3D("zx");
@@ -237,7 +244,7 @@ void drawQAPlots(const char* analysisOutput,
   gHistPhiTPCNPoints->SetStats(kFALSE);
 
   //3D pt-phi-NPoints: protons
-  TH3F *gHistProtonsPtPhiTPCNPoints = dynamic_cast<TH3F *>(fQA2DList->At(9));
+  TH3F *gHistProtonsPtPhiTPCNPoints = dynamic_cast<TH3F *>(fQA2DList->At(11));
   TH2D *gHistProtonsPtPhi = dynamic_cast<TH2D *>gHistProtonsPtPhiTPCNPoints->Project3D("yx");
   gHistProtonsPtPhi->SetStats(kFALSE);
   TH2D *gHistProtonsPtTPCNPoints = dynamic_cast<TH2D *>gHistProtonsPtPhiTPCNPoints->Project3D("zx");
@@ -246,19 +253,19 @@ void drawQAPlots(const char* analysisOutput,
   gHistProtonsPhiTPCNPoints->SetStats(kFALSE);
 
   //2D eta-phi- accepted protons & antiprotons
-  TH2F *gHistEtaPhiProtons = dynamic_cast<TH2F *>(fQA2DList->At(10));
+  TH2F *gHistEtaPhiProtons = dynamic_cast<TH2F *>(fQA2DList->At(12));
   gHistEtaPhiProtons->SetStats(kFALSE);
-  TH2F *gHistEtaPhiAntiProtons = dynamic_cast<TH2F *>(fQA2DList->At(11));
+  TH2F *gHistEtaPhiAntiProtons = dynamic_cast<TH2F *>(fQA2DList->At(13));
   gHistEtaPhiAntiProtons->SetStats(kFALSE);
 
   //2D dca vs pT - accepted protons & antiprotons
-  TH2F *gHistDCAxyPtProtons = dynamic_cast<TH2F *>(fQA2DList->At(12));
+  TH2F *gHistDCAxyPtProtons = dynamic_cast<TH2F *>(fQA2DList->At(14));
   gHistDCAxyPtProtons->SetStats(kFALSE);
-  TH2F *gHistDCAzPtProtons = dynamic_cast<TH2F *>(fQA2DList->At(13));
+  TH2F *gHistDCAzPtProtons = dynamic_cast<TH2F *>(fQA2DList->At(15));
   gHistDCAzPtProtons->SetStats(kFALSE);
-  TH2F *gHistDCAxyPtAntiProtons = dynamic_cast<TH2F *>(fQA2DList->At(14));
+  TH2F *gHistDCAxyPtAntiProtons = dynamic_cast<TH2F *>(fQA2DList->At(16));
   gHistDCAxyPtAntiProtons->SetStats(kFALSE);
-  TH2F *gHistDCAzPtAntiProtons = dynamic_cast<TH2F *>(fQA2DList->At(15));
+  TH2F *gHistDCAzPtAntiProtons = dynamic_cast<TH2F *>(fQA2DList->At(17));
   gHistDCAzPtAntiProtons->SetStats(kFALSE);
 
   //__________________________________________________//
@@ -1003,67 +1010,67 @@ void drawMCvsData(const char* analysisOutputMC,
   gDataAntiProtonsITSClustersPass->Add(gDataAntiProtonsITSClustersReject);
   gDataAntiProtonsITSClustersPass->Sumw2();
   gDataAntiProtonsITSClustersPass->Scale(1./gDataAntiProtonsITSClustersPass->Integral(1,gDataAntiProtonsITSClustersPass->GetNbinsX()));
-  gDataAntiProtonsITSClustersPass->SetMarkerStyle(24);
+  gDataAntiProtonsITSClustersPass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsTPCClustersPass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(2));
   gDataAntiProtonsTPCClustersPass->SetStats(kFALSE);
   gDataAntiProtonsTPCClustersPass->Add(gDataAntiProtonsTPCClustersReject);
   gDataAntiProtonsTPCClustersPass->Sumw2();
   gDataAntiProtonsTPCClustersPass->Scale(1./gDataAntiProtonsTPCClustersPass->Integral(1,gDataAntiProtonsTPCClustersPass->GetNbinsX()));
-  gDataAntiProtonsTPCClustersPass->SetMarkerStyle(24);
+  gDataAntiProtonsTPCClustersPass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsChi2PerClusterTPCPass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(3));
   gDataAntiProtonsChi2PerClusterTPCPass->SetStats(kFALSE);
   gDataAntiProtonsChi2PerClusterTPCPass->Add(gDataAntiProtonsChi2PerClusterTPCReject);
   gDataAntiProtonsChi2PerClusterTPCPass->Sumw2();
   gDataAntiProtonsChi2PerClusterTPCPass->Scale(1./gDataAntiProtonsChi2PerClusterTPCPass->Integral(1,gDataAntiProtonsChi2PerClusterTPCPass->GetNbinsX()));
-  gDataAntiProtonsChi2PerClusterTPCPass->SetMarkerStyle(24);
+  gDataAntiProtonsChi2PerClusterTPCPass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsExtCov11Pass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(4));
   gDataAntiProtonsExtCov11Pass->SetStats(kFALSE);
   gDataAntiProtonsExtCov11Pass->Add(gDataAntiProtonsExtCov11Reject);
   gDataAntiProtonsExtCov11Pass->Sumw2();
   gDataAntiProtonsExtCov11Pass->Scale(1./gDataAntiProtonsExtCov11Pass->Integral(1,gDataAntiProtonsExtCov11Pass->GetNbinsX()));
-  gDataAntiProtonsExtCov11Pass->SetMarkerStyle(24);
+  gDataAntiProtonsExtCov11Pass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsExtCov22Pass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(5));
   gDataAntiProtonsExtCov22Pass->SetStats(kFALSE);
   gDataAntiProtonsExtCov22Pass->Add(gDataAntiProtonsExtCov22Reject);
   gDataAntiProtonsExtCov22Pass->Sumw2();
   gDataAntiProtonsExtCov22Pass->Scale(1./gDataAntiProtonsExtCov22Pass->Integral(1,gDataAntiProtonsExtCov22Pass->GetNbinsX()));
-  gDataAntiProtonsExtCov22Pass->SetMarkerStyle(24);
+  gDataAntiProtonsExtCov22Pass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsExtCov33Pass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(6));
   gDataAntiProtonsExtCov33Pass->SetStats(kFALSE);
   gDataAntiProtonsExtCov33Pass->Add(gDataAntiProtonsExtCov33Reject);
   gDataAntiProtonsExtCov33Pass->Sumw2();
   gDataAntiProtonsExtCov33Pass->Scale(1./gDataAntiProtonsExtCov33Pass->Integral(1,gDataAntiProtonsExtCov33Pass->GetNbinsX()));
-  gDataAntiProtonsExtCov33Pass->SetMarkerStyle(24);
+  gDataAntiProtonsExtCov33Pass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsExtCov44Pass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(7));
   gDataAntiProtonsExtCov44Pass->SetStats(kFALSE);
   gDataAntiProtonsExtCov44Pass->Add(gDataAntiProtonsExtCov44Reject);
   gDataAntiProtonsExtCov44Pass->Sumw2();
   gDataAntiProtonsExtCov44Pass->Scale(1./gDataAntiProtonsExtCov44Pass->Integral(1,gDataAntiProtonsExtCov44Pass->GetNbinsX()));
-  gDataAntiProtonsExtCov44Pass->SetMarkerStyle(24);
+  gDataAntiProtonsExtCov44Pass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsExtCov55Pass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(8));
   gDataAntiProtonsExtCov55Pass->SetStats(kFALSE);
   gDataAntiProtonsExtCov55Pass->Add(gDataAntiProtonsExtCov55Reject);
   gDataAntiProtonsExtCov55Pass->Sumw2();
   gDataAntiProtonsExtCov55Pass->Scale(1./gDataAntiProtonsExtCov55Pass->Integral(1,gDataAntiProtonsExtCov55Pass->GetNbinsX()));
-  gDataAntiProtonsExtCov55Pass->SetMarkerStyle(24);
+  gDataAntiProtonsExtCov55Pass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsDCAXYPass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(11));
   gDataAntiProtonsDCAXYPass->SetStats(kFALSE);  
   gDataAntiProtonsDCAXYPass->Add(gDataAntiProtonsDCAXYReject);
   gDataAntiProtonsDCAXYPass->Sumw2();
   gDataAntiProtonsDCAXYPass->Scale(1./gDataAntiProtonsDCAXYPass->Integral(1,gDataAntiProtonsDCAXYPass->GetNbinsX()));
-  gDataAntiProtonsDCAXYPass->SetMarkerStyle(24);
+  gDataAntiProtonsDCAXYPass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsDCAZPass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(13));
   gDataAntiProtonsDCAZPass->SetStats(kFALSE);  
   gDataAntiProtonsDCAZPass->Add(gDataAntiProtonsDCAZReject);
   gDataAntiProtonsDCAZPass->Sumw2();
   gDataAntiProtonsDCAZPass->Scale(1./gDataAntiProtonsDCAZPass->Integral(1,gDataAntiProtonsDCAZPass->GetNbinsX()));
-  gDataAntiProtonsDCAZPass->SetMarkerStyle(24);
+  gDataAntiProtonsDCAZPass->SetMarkerStyle(1);
   TH1F *gDataAntiProtonsNumberOfTPCdEdxPointsPass = dynamic_cast<TH1F *>(fQADataAntiProtonsAcceptedList->At(26));
   gDataAntiProtonsNumberOfTPCdEdxPointsPass->SetStats(kFALSE);  
   gDataAntiProtonsNumberOfTPCdEdxPointsPass->Add(gDataAntiProtonsNumberOfTPCdEdxPointsReject);
   gDataAntiProtonsNumberOfTPCdEdxPointsPass->Sumw2();
   gDataAntiProtonsNumberOfTPCdEdxPointsPass->Scale(1./gDataAntiProtonsNumberOfTPCdEdxPointsPass->Integral(1,gDataAntiProtonsNumberOfTPCdEdxPointsPass->GetNbinsX()));
-  gDataAntiProtonsNumberOfTPCdEdxPointsPass->SetMarkerStyle(24);
+  gDataAntiProtonsNumberOfTPCdEdxPointsPass->SetMarkerStyle(1);
 
 
   //__________________________________________________//