1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
16 /* $Id: AliAnalysisTaskPIDResponse.cxx 43811 2010-09-23 14:13:31Z wiechula $ */
20 #include <AliAnalysisManager.h>
21 #include <AliInputEventHandler.h>
22 #include <AliVEventHandler.h>
23 #include <AliVEvent.h>
24 #include <AliVParticle.h>
25 #include <AliVTrack.h>
27 #include <AliPIDResponse.h>
28 #include <AliESDpid.h>
30 #include "AliAnalysisTaskPIDResponse.h"
32 ClassImp(AliAnalysisTaskPIDResponse)
34 //______________________________________________________________________________
35 AliAnalysisTaskPIDResponse::AliAnalysisTaskPIDResponse():
40 fSpecialDetResponse(),
45 fIsTunedOnData(kFALSE),
53 //______________________________________________________________________________
54 AliAnalysisTaskPIDResponse::AliAnalysisTaskPIDResponse(const char* name):
55 AliAnalysisTaskSE(name),
59 fSpecialDetResponse(),
64 fIsTunedOnData(kFALSE),
68 // Default constructor
70 DefineInput(0,TChain::Class());
73 //______________________________________________________________________________
74 AliAnalysisTaskPIDResponse::~AliAnalysisTaskPIDResponse()
81 //______________________________________________________________________________
82 void AliAnalysisTaskPIDResponse::UserCreateOutputObjects()
85 // Create the output QA objects
88 AliLog::SetClassDebugLevel("AliAnalysisTaskPIDResponse",10);
91 AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
92 AliInputEventHandler *inputHandler=dynamic_cast<AliInputEventHandler*>(man->GetInputEventHandler());
93 if (!inputHandler) AliFatal("Input handler needed");
96 inputHandler->CreatePIDResponse(fIsMC);
97 fPIDResponse=inputHandler->GetPIDResponse();
98 if (!fPIDResponse) AliFatal("PIDResponse object was not created");
100 fPIDResponse->SetOADBPath(AliAnalysisManager::GetOADBPath());
101 // fPIDResponse->SetCachePID(fCachePID);
102 if (!fOADBPath.IsNull()) fPIDResponse->SetOADBPath(fOADBPath.Data());
104 if(fIsTunedOnData) fPIDResponse->SetTunedOnData(kTRUE,fRecoPassTuned);
106 if (!fSpecialDetResponse.IsNull()){
107 TObjArray *arr=fSpecialDetResponse.Tokenize("; ");
108 for (Int_t i=0; i<arr->GetEntriesFast();++i){
109 TString resp(arr->At(i)->GetName());
110 if (resp.BeginsWith("TPC:")){
111 resp.ReplaceAll("TPC:","");
112 fPIDResponse->SetCustomTPCpidResponse(resp.Data());
113 AliInfo(Form("Setting custom TPC response file: '%s'",resp.Data()));
119 //______________________________________________________________________________
120 void AliAnalysisTaskPIDResponse::UserExec(Option_t */*option*/)
122 // Setup the PID response functions and fill the QA histograms
124 AliVEvent *event=InputEvent();
126 fRun=event->GetRunNumber();
133 fPIDResponse->InitialiseEvent(event,fRecoPass);
134 AliESDpid *pidresp = dynamic_cast<AliESDpid*>(fPIDResponse);
135 if(pidresp && AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()){
136 pidresp->SetEventHandler(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler());
138 //create and attach transient PID object
140 fPIDResponse->FillTrackDetectorPID();
144 //______________________________________________________________________________
145 void AliAnalysisTaskPIDResponse::SetRecoInfo()
148 // Set reconstruction information
154 //Get the current file to check the reconstruction pass (UGLY, but not stored in ESD... )
155 AliAnalysisManager *mgr=AliAnalysisManager::GetAnalysisManager();
156 AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
157 if (!inputHandler) return;
159 TTree *tree= (TTree*)inputHandler->GetTree();
160 TFile *file= (TFile*)tree->GetCurrentFile();
163 AliError("Current file not found, cannot set reconstruction information");
167 //find pass from file name (UGLY, but not stored in ESD... )
168 TString fileName(file->GetName());
169 if (fileName.Contains("pass1") ) {
171 } else if (fileName.Contains("pass2") ) {
173 } else if (fileName.Contains("pass3") ) {
175 } else if (fileName.Contains("pass4") ) {
177 } else if (fileName.Contains("pass5") ) {
181 fPIDResponse->SetCurrentFile(fileName.Data());