PWGCF/FLOW/macros/AddTaskPidVZEROSP.C

   1 AddTaskPidVZEROSP(Int_t centralityselection=AliVEvent::kAny,Float_t etamin=-0.8,Float_t etamax=0.8,Int_t side=0,Int_t filterbit=1,Bool_t TOFbeta=kFALSE){
   2   gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityPIDSP.C");
   3
   4   const Int_t ncentr = 5;
   5   Int_t cmin[ncentr]={0,5,10,20,40};
   6   Int_t cmax[ncentr]={5,10,20,30,50};
   7
   8   for(Int_t i=0;i < ncentr;i++){
   9     if(!TOFbeta){
  10       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side,filterbit); // no pid
  11       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side,filterbit);
  12       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side,filterbit);
  13       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax,"",side,filterbit);
  14       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax,"",side,filterbit);
  15     }
  16     else{
  17       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbeta,0,2,0,etamin,etamax,"",side,filterbit); // no pid
  18       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbeta,0,2,0,etamin,etamax,"",side,filterbit);
  19       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbeta,0,2,0,etamin,etamax,"",side,filterbit);
  20       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbeta,-1,2,0,etamin,etamax,"",side,filterbit);
  21       AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbeta,0,2,0,etamin,etamax,"",side,filterbit);
  22     }
  23     AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax,"",side,filterbit);
  24     AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax,"",side,filterbit);
  25     AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax,"",side,filterbit);
  26     AddTaskFlowCentralityPIDSP(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax,"",side,filterbit);
  27   }
  28 }
  29
  30 createSPres(){
  31   gSystem->Load("libVMC");
  32   gSystem->Load("libPhysics");
  33   gSystem->Load("libTree");
  34   gSystem->Load("libMinuit");
  35   gSystem->Load("libSTEERBase");
  36   gSystem->Load("libANALYSIS");
  37   gSystem->Load("libAOD");
  38   gSystem->Load("libESD");
  39   gSystem->Load("libANALYSIS");
  40   gSystem->Load("libANALYSISalice");
  41   gSystem->Load("libCORRFW");
  42   gSystem->Load("libNetx");
  43   gSystem->Load("libPWGflowBase");
  44
  45   char name[200];
  46   char *spe[4]={"pion","kaon","antipr","proton"};
  47   char *tech[2]={"TOF","TPC"};
  48   const Int_t ncentr = 5;
  49   Int_t cmin[ncentr]={0,5,10,20,40};
  50   Int_t cmax[ncentr]={5,10,20,30,50};
  51
  52   TFile *f = new TFile("AnalysisResults.root");
  53   TFile *fo = new TFile("results.root","RECREATE");
  54   TDirectory* directory = dynamic_cast<TDirectory*>(f->Get("outputSPanalysisTPCstandalone"));
  55   TList* listTemp = directory->GetListOfKeys();
  56   for(Int_t i=0;i < ncentr;i++){
  57     TList* list2 = dynamic_cast<TList*>(directory->Get(listTemp->At(i*ncentr)->GetName()));
  58     AliFlowAnalysisWithScalarProduct* sp2 = new AliFlowAnalysisWithScalarProduct();
  59     sp2->GetOutputHistograms(list2);
  60     sp2->Finish();
  61     AliFlowCommonHistResults* res2=sp2->GetCommonHistsRes();
  62     TH1D *h2=res2->GetHistDiffFlowPtPOI();
  63     sprintf(name,"v2SP_%s_%i_%i","AllCharged",cmin[i],cmax[i]);
  64     h2->SetName(name);
  65     fo->cd();
  66     h2->Write();
  67     for(Int_t j=0;j < 2;j++){
  68       for(Int_t k=0;k < 4;k++){
  69         TList* list = dynamic_cast<TList*>(directory->Get(listTemp->At(i*ncentr+j*4+k+1)->GetName()));
  70         AliFlowAnalysisWithScalarProduct* sp = new AliFlowAnalysisWithScalarProduct();
  71         sp->GetOutputHistograms(list);
  72         sp->Finish();
  73         AliFlowCommonHistResults* res=sp->GetCommonHistsRes();
  74         TH1D *h=res->GetHistDiffFlowPtPOI();
  75         sprintf(name,"v2SP_%s_%i_%i%s",spe[k],cmin[i],cmax[i],tech[j]);
  76         h->SetName(name);
  77         fo->cd();
  78         h->Write();
  79       }
  80     }
  81   }
  82   fo->Close();
  83 }
  84
  85
  86
  87