3 gSystem->SetIncludePath("-I. -I$ROOTSYS/include -I$ALICE_ROOT -I$ALICE_ROOT/build/include -I$ALICE_ROOT/include -I$ALICE_ROOT/ITS -I$ALICE_ROOT/TPC -I$ALICE_ROOT/CONTAINERS -I$ALICE_ROOT/STEER -I$ALICE_ROOT/TRD -I$ALICE_ROOT/macros -I$ALICE_ROOT/ANALYSIS -I$ALICE_ROOT/PWG2/FLOW/AliFlowCommon -I$ALICE_ROOT/PWG2/FLOW/AliFlowTasks -I$ALICE_ROOT/PWG3/dielectron/ -g");
5 gSystem->Load("libCore");// no
6 gSystem->Load("libTree");
7 gSystem->Load("libGeom");
8 gSystem->Load("libVMC");
9 gSystem->Load("libXMLIO");// no
10 gSystem->Load("libPhysics");
11 gSystem->Load("libXMLParser");
12 gSystem->Load("libProof");
13 gSystem->Load("libMinuit");
15 gSystem->Load("libSTEERBase");
16 gSystem->Load("libCDB");
17 gSystem->Load("libRAWDatabase");
18 gSystem->Load("libRAWDatarec");
19 gSystem->Load("libESD");
20 gSystem->Load("libAOD");
21 gSystem->Load("libSTEER");
22 gSystem->Load("libANALYSIS");
23 gSystem->Load("libANALYSISalice");
24 gSystem->Load("libTOFbase");
25 gSystem->Load("libTOFrec");
26 gSystem->Load("libT0base");
27 gSystem->Load("libT0rec");
28 gSystem->Load("libPWG2flowCommon");
29 gSystem->Load("libPWG2flowTasks");
31 gSystem->Load("libTender");
32 gSystem->Load("libTenderSupplies");
36 gSystem->Load("libCORRFW.so");
37 gSystem->Load("libPWG3base.so");
38 gSystem->Load("libPWG3dielectron.so");
39 gSystem->Load("libPWG3hfe.so");
42 TChain *chain = new TChain("esdTree");
43 for(int i=1; i!=4; ++i)
44 chain->Add( Form("/home/gunji/softwares/dielectron/data/esd137549035/%d0/AliESDs.root",i) );
46 AliAnalysisManager *mgr = new AliAnalysisManager("DielectronAnalysisManager");
47 AliESDInputHandler *esdH = new AliESDInputHandler();
48 mgr->SetInputEventHandler(esdH);
49 gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
50 AddTaskPhysicsSelection(kFALSE);
51 gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskCentrality.C");
52 AliCentralitySelectionTask *taskCentrality =AddTaskCentrality();
53 // taskCentrality->SetPass(2);
55 /* gROOT->LoadMacro("AliDielectronDebugTreeTaku.cxx++");
56 gROOT->LoadMacro("AliDielectronHistosTaku.cxx++");
57 gROOT->LoadMacro("AliDielectronTaku.cxx++");
58 gROOT->LoadMacro("AliAnalysisTaskMultiDielectronNewTaku.cxx++");
61 gSystem->Load("./AliDielectronHistosTaku_cxx.so");
62 gSystem->Load("./AliDielectronDebugTreeTaku_cxx.so");
63 gSystem->Load("./AliDielectronTaku_cxx.so");
72 AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
73 AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler());
74 Bool_t isESD=man->GetInputEventHandler()->IsA()==AliESDInputHandler::Class();
76 if ( inputHandler->GetPIDResponse() ){
77 AliDielectronVarManager::SetPIDResponse( inputHandler->GetPIDResponse() );
80 if (!AliDielectronVarManager::GetESDpid()){
81 if (AliDielectronMC::Instance()->HasMC()) {
82 AliDielectronVarManager::InitESDpid();
84 cout<<" set pid as 1"<<endl;
85 AliDielectronVarManager::InitESDpid(1);
91 AliPIDResponse * f= (AliPIDResponse*) AliDielectronVarManager::GetESDpid() ;
92 AliTPCPIDResponse *ff = (AliTPCPIDResponse*)f->GetTPCResponse();
94 for (Int_t j=0; j<AliPID::kSPECIES; j++) {
95 //AliPID::EParticleType type=AliPID::EParticleType(j);
96 for(int ii=0;ii<100;ii++){
97 float p = 0.05+0.1*ii;
98 //Double_t bethe=fTPCResponse.GetExpectedSignal(mom,type);
99 Double_t bethe=ff->GetExpectedSignal(p,j);
103 //cout<<j<<" "<<p<<" "<<bethe<<endl;
108 Double_t fAlephParam[5]={2.11543/57,
116 Double_t fAlephParam[5]={2.11543/122,
124 AliESDpid *fESDpid = new AliESDpid();
125 fESDpid->GetTPCResponse().SetBetheBlochParameters(fAlephParam[0],
132 AliTPCPIDResponse *ff = (AliTPCPIDResponse*)fESDpid->GetTPCResponse();
134 for (Int_t j=0; j<AliPID::kSPECIES; j++) {
135 for(int ii=0;ii<200;ii++){
136 float p = 0.025+0.05*ii;
137 Double_t bethe=ff->GetExpectedSignal(p,j);
140 cout<<j<<" "<<p<<" "<<bethe<<endl;
147 for(int i=0;i<5;i++){
148 g[i] = new TGraph(200, mom[i], dedx[i]);
149 sprintf(name,"g%d",i);
151 g[i]->SetLineColor(i+2);
154 TFile *fin=new TFile("ana/tmp.root");
156 hdedx_pt->SetAxisRange(0,4);
163 TFile *fout = new TFile("dedx.root","recreate");
165 for(int i=0;i<5;i++){