X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PWGPP%2Fpid%2FAliAnalysisTaskPhiBayes.cxx;h=eeaa0a0a28bd64c9bf83384a47b9125246d5a321;hb=e72663d5ea704f30a9a1df4852c98d2cb07d29d7;hp=b3b3b428141f63163588d3458b4b750b5638fb01;hpb=30235f51e268837f42111dd7d2e434efcde9e00a;p=u%2Fmrichter%2FAliRoot.git

diff --git a/PWGPP/pid/AliAnalysisTaskPhiBayes.cxx b/PWGPP/pid/AliAnalysisTaskPhiBayes.cxx
index b3b3b428141..eeaa0a0a28b 100644
--- a/PWGPP/pid/AliAnalysisTaskPhiBayes.cxx
+++ b/PWGPP/pid/AliAnalysisTaskPhiBayes.cxx
@@ -497,7 +497,8 @@ void AliAnalysisTaskPhiBayes::Analyze(AliAODEvent* aodEvent)
   Double_t Qx2 = 0, Qy2 = 0;
   Double_t Qx3 = 0, Qy3 = 0;
   for(Int_t iT = 0; iT < ntrack; iT++) {
-    AliAODTrack* aodTrack = aodEvent->GetTrack(iT);
+    AliAODTrack* aodTrack = dynamic_cast<AliAODTrack*>(aodEvent->GetTrack(iT));
+    if(!aodTrack) AliFatal("Not a standard AOD");
     
     if (!aodTrack){
       continue;
@@ -528,6 +529,8 @@ void AliAnalysisTaskPhiBayes::Analyze(AliAODEvent* aodEvent)
   AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
   AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler());
   AliPIDResponse *PIDResponse=inputHandler->GetPIDResponse();
+  PIDResponse->SetTOFResponse(aodEvent,AliPIDResponse::kTOF_T0);
+  PIDResponse->GetTOFResponse().SetTOFtailAllPara(-23,1.1);
 
 //   PIDResponse->GetTOFResponse().SetTrackParameter(0,0.);
 //   PIDResponse->GetTOFResponse().SetTrackParameter(1,0.);
@@ -551,7 +554,8 @@ void AliAnalysisTaskPhiBayes::Analyze(AliAODEvent* aodEvent)
     nmc = mcArray->GetEntries();
 
   for(Int_t i=0;i < ntrack;i++){ // loop on tracks
-    AliAODTrack* track = aodEvent->GetTrack(i);
+    AliAODTrack* track = dynamic_cast<AliAODTrack*>(aodEvent->GetTrack(i));
+    if(!track) AliFatal("Not a standard AOD");
     
     AliAODMCParticle *mcp = NULL;
     Int_t pdg = 0;
@@ -616,7 +620,8 @@ void AliAnalysisTaskPhiBayes::Analyze(AliAODEvent* aodEvent)
 
   // start analysis phi
   for(Int_t i=0;i < ntrack;i++){ // loop on positive tracks
-    AliAODTrack* KpTrack = aodEvent->GetTrack(i);
+    AliAODTrack* KpTrack = dynamic_cast<AliAODTrack*>(aodEvent->GetTrack(i));
+    if(!KpTrack) AliFatal("Not a standard AOD");
         
     if (!KpTrack){
       continue;
@@ -739,7 +744,8 @@ void AliAnalysisTaskPhiBayes::Analyze(AliAODEvent* aodEvent)
     if(fPtKp > 4.299) fPtKp = 4.299;
 
     for(Int_t j=0;j < ntrack;j++){ // loop on negative tracks
-      AliAODTrack* KnTrack = aodEvent->GetTrack(j);
+      AliAODTrack* KnTrack = dynamic_cast<AliAODTrack*>(aodEvent->GetTrack(j));
+      if(!KnTrack) AliFatal("Not a standard AOD");
       
       if (!KnTrack){
 	continue;
@@ -897,8 +903,8 @@ void AliAnalysisTaskPhiBayes::Analyze(AliAODEvent* aodEvent)
 
       if(fPtKn > 4.299) fPtKn = 4.299;
 
-      Float_t xTOfill[] = {fPtPhi,KpTrack->Eta(),fPtKp,fPtKn,(fPidKp%128)*0.01,(fPidKn%128)*0.01,tofMatch1,tofMatch2,isTrue,nSigmaComb,nSigmaComb2,deltaphi1,deltaphi2,fPsi};
-      Float_t xTOfill2[] = {fPtPhi,KpTrack->Eta(),fPtKp,fPtKn,(fPidKp%128)*0.01,(fPidKn%128)*0.01,tofMatch1,tofMatch2,isTrue,nSigmaComb,nSigmaComb2,deltaphi1,deltaphi2,fPsi};
+      Float_t xTOfill[] = {static_cast<Float_t>(fPtPhi),static_cast<Float_t>(KpTrack->Eta()),static_cast<Float_t>(fPtKp),static_cast<Float_t>(fPtKn),static_cast<Float_t>((fPidKp%128)*0.01),static_cast<Float_t>((fPidKn%128)*0.01),static_cast<Float_t>(tofMatch1),static_cast<Float_t>(tofMatch2),static_cast<Float_t>(isTrue),static_cast<Float_t>(nSigmaComb),static_cast<Float_t>(nSigmaComb2),static_cast<Float_t>(deltaphi1),static_cast<Float_t>(deltaphi2),static_cast<Float_t>(fPsi)};
+      Float_t xTOfill2[] = {static_cast<Float_t>(fPtPhi),static_cast<Float_t>(KpTrack->Eta()),static_cast<Float_t>(fPtKp),static_cast<Float_t>(fPtKn),static_cast<Float_t>((fPidKp%128)*0.01),static_cast<Float_t>((fPidKn%128)*0.01),static_cast<Float_t>(tofMatch1),static_cast<Float_t>(tofMatch2),static_cast<Float_t>(isTrue),static_cast<Float_t>(nSigmaComb),static_cast<Float_t>(nSigmaComb2),static_cast<Float_t>(deltaphi1),static_cast<Float_t>(deltaphi2),static_cast<Float_t>(fPsi)};
       
       Int_t ipt = 0;
       while(fPtPhiMin[ipt] < fPtPhi && ipt < nPtBin){
@@ -922,8 +928,8 @@ void AliAnalysisTaskPhiBayes::Analyze(AliAODEvent* aodEvent)
 	if((fPidKp%128) > 80) fContPid2->Fill(0,fMassV0,fCentrality,xTOfill2);// use tagging on positive track
 
 	if(fPIDuserCut){
-	  Float_t xUser[] = {KpTrack->Eta(),fPtKp,isTrue,0,deltaphi1,fPsi};
-	  Float_t xUser2[] = {KnTrack->Eta(),fPtKn,isTrue,0,deltaphi2,fPsi};
+	  Float_t xUser[] = {static_cast<Float_t>(KpTrack->Eta()),static_cast<Float_t>(fPtKp),static_cast<Float_t>(isTrue),0,static_cast<Float_t>(deltaphi1),static_cast<Float_t>(fPsi)};
+	  Float_t xUser2[] = {static_cast<Float_t>(KnTrack->Eta()),static_cast<Float_t>(fPtKn),static_cast<Float_t>(isTrue),0,static_cast<Float_t>(deltaphi2),static_cast<Float_t>(fPsi)};
 	  
 	  if(fPIDuserCut->IsSelected(KpTrack,AliPID::kPion)){ // to be filled for positive
 	    xUser[3] = 1;