- histos with expected TPC signal added

diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODPID.cxx b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODPID.cxx
index 5bb19cd58247ac88b598274db2fcd006f5f5273b..a546eccc0e35ff98f718b125c32a12fec693f09e 100644 (file)
--- a/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODPID.cxx
+++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/AliSpectraAODPID.cxx
@@ -38,10 +38,12 @@ void AliSpectraAODPID::FillQAHistos(AliSpectraAODHistoManager * hman, AliAODTrac
   
   //Response
   AliVParticle *inEvHMain = dynamic_cast<AliVParticle *>(track);
-
+  
   hman->GetPIDHistogram(kHistPIDTPC)->Fill(track->GetTPCmomentum(), track->GetTPCsignal()*track->Charge()); // PID histo
-
-
+  hman->GetPIDHistogram(kHistPIDTPCPion)->Fill(track->GetTPCmomentum(),fPIDResponse->GetTPCResponse().GetExpectedSignal(track->GetTPCmomentum(),AliPID::kPion)*track->Charge()); // Expected PIDPion histo
+  hman->GetPIDHistogram(kHistPIDTPCKaon)->Fill(track->GetTPCmomentum(),fPIDResponse->GetTPCResponse().GetExpectedSignal(track->GetTPCmomentum(),AliPID::kKaon)*track->Charge()); // Expected PIDKaon histo
+  hman->GetPIDHistogram(kHistPIDTPCProton)->Fill(track->GetTPCmomentum(),fPIDResponse->GetTPCResponse().GetExpectedSignal(track->GetTPCmomentum(),AliPID::kProton)*track->Charge()); // Expected PIDProton histo
+  
   Double_t nsigmaTPCkProton = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(inEvHMain, AliPID::kProton));
   Double_t nsigmaTPCkKaon = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(inEvHMain, AliPID::kKaon)); 
   Double_t nsigmaTPCkPion = TMath::Abs(fPIDResponse->NumberOfSigmasTPC(inEvHMain, AliPID::kPion)); 
@@ -49,6 +51,7 @@ void AliSpectraAODPID::FillQAHistos(AliSpectraAODHistoManager * hman, AliAODTrac
 
   if(track->Pt()>trackCuts->GetPtTOFMatching()){
     hman->GetPIDHistogram(kHistPIDTOF)->Fill(track->P(),(track->GetTOFsignal()/100)*track->Charge()); // PID histo
+    
     nsigmaTOFkProton = TMath::Abs(fPIDResponse->NumberOfSigmasTOF(inEvHMain, AliPID::kProton));
     nsigmaTOFkKaon = TMath::Abs(fPIDResponse->NumberOfSigmasTOF(inEvHMain, AliPID::kKaon)); 
     nsigmaTOFkPion = TMath::Abs(fPIDResponse->NumberOfSigmasTOF(inEvHMain, AliPID::kPion)); 

diff --git a/PWGLF/SPECTRA/PiKaPr/TestAOD/QAPlots.C b/PWGLF/SPECTRA/PiKaPr/TestAOD/QAPlots.C
index d6a4d0e7b25f65243aa7288ce3753ea693c14b16..0dee9bba87a0df88448ccb38e16ef8026e7a606b 100644 (file)
--- a/PWGLF/SPECTRA/PiKaPr/TestAOD/QAPlots.C
+++ b/PWGLF/SPECTRA/PiKaPr/TestAOD/QAPlots.C
@@ -149,6 +149,12 @@ void QAPlots( AliSpectraAODHistoManager* hman_data, AliSpectraAODHistoManager* h
   gPad->SetGridy();
   gPad->SetGridx();
   PIDSig_data->DrawClone("colz");
+  for(Int_t ipart=0;ipart<3;ipart++){
+    TH2F *PIDSig_dataPart = (TH2F*)((TH2F*)hman_data->GetPIDHistogram(Form("hHistPIDTPC%s",Particle[ipart].Data())))->Clone();
+    PIDSig_dataPart->SetMarkerStyle(22);
+    PIDSig_dataPart->SetMarkerSize(.5);
+    PIDSig_dataPart->DrawClone("same");
+  }
   cPIDSig->cd(2);
   TH2F *PIDSig_mc = (TH2F*)((TH2F*)hman_mc->GetPIDHistogram("hHistPIDTPC"))->Clone();
   PIDSig_mc->SetYTitle("TPC signal");
@@ -156,6 +162,12 @@ void QAPlots( AliSpectraAODHistoManager* hman_data, AliSpectraAODHistoManager* h
   gPad->SetGridy();
   gPad->SetGridx();
   PIDSig_mc->DrawClone("colz");
+  for(Int_t ipart=0;ipart<3;ipart++){
+    TH2F *PIDSig_mcPart = (TH2F*)((TH2F*)hman_mc->GetPIDHistogram(Form("hHistPIDTPC%s",Particle[ipart].Data())))->Clone();
+    PIDSig_mcPart->SetMarkerStyle(22);
+    PIDSig_mcPart->SetMarkerSize(.5);
+    PIDSig_mcPart->DrawClone("same");
+  }
   cPIDSig->cd(3);
   TH2F *PIDSig_data = (TH2F*)((TH2F*)hman_data->GetPIDHistogram("hHistPIDTOF"))->Clone();
   PIDSig_data->SetYTitle("TOF signal");
@@ -337,8 +349,7 @@ void QAPlots( AliSpectraAODHistoManager* hman_data, AliSpectraAODHistoManager* h
     nsig_mc_Proj2->DrawNormalized("lhistsame");
     gPad->BuildLegend();
   }
-  return;
-
+  
   //Muon over Pion Ratio
   Printf("\n\n-> Muon Over Pion");
   TCanvas *cMu=new TCanvas("cMu","cMu");