bugfix (dead code)

[u/mrichter/AliRoot.git] / PWG / Tools / AliHelperPID.cxx

diff --git a/PWG/Tools/AliHelperPID.cxx b/PWG/Tools/AliHelperPID.cxx
index 875a97dca7a43357e67bbadb16eb57e45c0241b3..3409cce79846e3138c43e3567b0f05cce8f4644f 100644 (file)
--- a/PWG/Tools/AliHelperPID.cxx
+++ b/PWG/Tools/AliHelperPID.cxx
@@ -44,7 +44,7 @@ const char * kParticleSpeciesName[]={"Pions","Kaons","Protons","Undefined"} ;
 
 ClassImp(AliHelperPID)
 
-AliHelperPID::AliHelperPID() : TNamed("HelperPID", "PID object"),fisMC(0), fPIDType(kNSigmaTPCTOF), fNSigmaPID(3), fBayesCut(0.8), fPIDResponse(0), fPIDCombined(0),fOutputList(0),fRequestTOFPID(1),fRemoveTracksT0Fill(0),fUseExclusiveNSigma(0),fPtTOFPID(.6),fHasTOFPID(0){
+AliHelperPID::AliHelperPID() : TNamed("HelperPID", "PID object"),fisMC(0), fPIDType(kNSigmaTPCTOF), fNSigmaPID(3), fBayesCut(0.8), fPIDResponse(0x0), fPIDCombined(0x0),fOutputList(0x0),fRequestTOFPID(1),fRemoveTracksT0Fill(0),fUseExclusiveNSigma(0),fPtTOFPID(.6),fHasTOFPID(0){
   
   for(Int_t ipart=0;ipart<kNSpecies;ipart++)
     for(Int_t ipid=0;ipid<=kNSigmaPIDType;ipid++)
@@ -133,6 +133,15 @@ AliHelperPID::AliHelperPID() : TNamed("HelperPID", "PID object"),fisMC(0), fPIDT
       fOutputList->Add(fHistoPID);
     }
   }
+  //PID signal plot, before PID cut
+  for(Int_t idet=0;idet<kNDetectors;idet++){
+    Double_t maxy=500;
+    if(idet==kTOF)maxy=1.1;
+    TH2F *fHistoPID=new TH2F(Form("PIDAll_%d",idet),Form("%s signal",kDetectorName[idet]),200,0,10,500,-maxy,maxy);
+    fHistoPID->GetXaxis()->SetTitle("P (GeV / c)");
+    fHistoPID->GetYaxis()->SetTitle(Form("%s signal",kDetectorName[idet]));
+    fOutputList->Add(fHistoPID);
+  }
 }
 
 //////////////////////////////////////////////////////////////////////////////////////////////////
@@ -167,7 +176,13 @@ Int_t AliHelperPID::GetParticleSpecies(AliVTrack * trk, Bool_t FIllQAHistos){
   if(fPIDType==kBayes){//use bayesianPID
     
     if(!fPIDCombined) {
-      AliFatal("PIDCombined object has to be set in the steering macro");
+      // ------- setup PIDCombined
+      fPIDCombined=new AliPIDCombined;
+      fPIDCombined->SetDefaultTPCPriors();
+      fPIDCombined->SetDetectorMask(AliPIDResponse::kDetTPC+AliPIDResponse::kDetTOF);  
+    }
+    if(!fPIDCombined) {
+      AliFatal("PIDCombined object not found");
     }
     
     ID = GetIDBayes(trk,FIllQAHistos);
@@ -194,6 +209,14 @@ Int_t AliHelperPID::GetParticleSpecies(AliVTrack * trk, Bool_t FIllQAHistos){
       if(idet==kTOF && fHasTOFPID)h->Fill(trk->P(),TOFBetaCalc(trk)*trk->Charge());
     }
   }
+  //Fill PID signal plot without cuts
+  for(Int_t idet=0;idet<kNDetectors;idet++){
+    TH2F *h=GetHistogram2D(Form("PIDAll_%d",idet));
+    if(idet==kITS)h->Fill(trk->P(),trk->GetITSsignal()*trk->Charge());
+    if(idet==kTPC)h->Fill(trk->P(),trk->GetTPCsignal()*trk->Charge());
+    if(idet==kTOF && fHasTOFPID)h->Fill(trk->P(),TOFBetaCalc(trk)*trk->Charge());
+  }
+  
   return ID;
 }