PWGCF/FLOW/macros/AddTaskPidVZEROSP_2.C

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