PWGCF/FLOW/macros/AddTaskPidQC.C

   1 AddTaskPidQC(Int_t centralityselection=AliVEvent::kAny,Float_t etamin=-0.8,Float_t etamax=0.8){
   2   gROOT->LoadMacro("$ALICE_ROOT/PWGCF/FLOW/macros/AddTaskFlowCentralityPIDQC.C");
   3
   4   const Int_t ncentr = 9;
   5   Int_t cmin[ncentr]={0,5,10,20,30,40,50,60,70};
   6   Int_t cmax[ncentr]={5,10,20,30,40,50,60,70,80};
   7
   8   for(Int_t i=0;i < ncentr;i++){
   9     AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kFALSE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,1,etamin,etamax); // no pid + QA
  10     AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);
  11     AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);
  12     AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,-1,2,0,etamin,etamax);
  13     AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTOFbayesian,0,2,0,etamin,etamax);
  14     AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kPion,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);
  15     AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kKaon,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);
  16     AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,-1,2,0,etamin,etamax);
  17     AddTaskFlowCentralityPIDQC(centralityselection,cmin[i],cmax[i],"AnalysisResults",kTRUE,AliPID::kProton,AliFlowTrackCuts::kTPCbayesian,0,2,0,etamin,etamax);
  18   }
  19 }
  20
  21 createQC4res(){
  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   Int_t cmin[9]={0,5,10,20,30,40,50,60,70};
  40   Int_t cmax[9]={5,10,20,30,40,50,60,70,80};
  41
  42   TFile *f = new TFile("AnalysisResults.root");
  43   TFile *fo = new TFile("resultsQC4.root","RECREATE");
  44   TDirectory* directory = dynamic_cast<TDirectory*>(f->Get("outputQCanalysisTPCstandalone"));
  45   TList* listTemp = directory->GetListOfKeys();
  46   for(Int_t i=0;i < 9;i++){
  47     TList* list2 = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9)->GetName()));
  48     AliFlowAnalysisWithQCumulants* sp2 = new AliFlowAnalysisWithQCumulants();
  49     sp2->GetOutputHistograms(list2);
  50     sp2->Finish();
  51     AliFlowCommonHistResults* res2=sp2->GetCommonHistsResults4th();
  52     TH1D *h2=res2->GetHistDiffFlowPtPOI();
  53     sprintf(name,"v2QC4_%s_%i_%i","AllCharged",cmin[i],cmax[i]);
  54     h2->SetName(name);
  55     fo->cd();
  56     h2->Write();
  57     for(Int_t j=0;j < 2;j++){
  58       for(Int_t k=0;k < 4;k++){
  59         TList* list = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9+j*4+k+1)->GetName()));
  60         AliFlowAnalysisWithQCumulants* sp = new AliFlowAnalysisWithQCumulants();
  61         sp->GetOutputHistograms(list);
  62         sp->Finish();
  63         AliFlowCommonHistResults* res=sp->GetCommonHistsResults4th();
  64         TH1D *h=res->GetHistDiffFlowPtPOI();
  65         sprintf(name,"v2QC4_%s_%i_%i%s",spe[k],cmin[i],cmax[i],tech[j]);
  66         h->SetName(name);
  67         fo->cd();
  68         h->Write();
  69       }
  70     }
  71   }
  72   fo->Close();
  73 }
  74
  75
  76
  77 createQC2res(){
  78   gSystem->Load("libVMC.so");
  79   gSystem->Load("libPhysics.so");
  80   gSystem->Load("libTree.so");
  81   gSystem->Load("libMinuit.so");
  82   gSystem->Load("libSTEERBase.so");
  83   gSystem->Load("libANALYSIS.so");
  84   gSystem->Load("libAOD.so");
  85   gSystem->Load("libESD.so");
  86   gSystem->Load("libANALYSIS.so");
  87   gSystem->Load("libANALYSISalice.so");
  88   gSystem->Load("libCORRFW.so");
  89   gSystem->Load("libNetx.so");
  90   gSystem->Load("libPWGflowBase.so");
  91
  92   char name[200];
  93   char *spe[4]={"pion","kaon","antipr","proton"};
  94   char *tech[2]={"TOF","TPC"};
  95   Int_t cmin[9]={0,5,10,20,30,40,50,60,70};
  96   Int_t cmax[9]={5,10,20,30,40,50,60,70,80};
  97
  98   TFile *f = new TFile("AnalysisResults.root");
  99   TFile *fo = new TFile("resultsQC2.root","RECREATE");
 100   TDirectory* directory = dynamic_cast<TDirectory*>(f->Get("outputQCanalysisTPCstandalone"));
 101   TList* listTemp = directory->GetListOfKeys();
 102   for(Int_t i=0;i < 9;i++){
 103     TList* list2 = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9)->GetName()));
 104     AliFlowAnalysisWithQCumulants* sp2 = new AliFlowAnalysisWithQCumulants();
 105     sp2->GetOutputHistograms(list2);
 106     sp2->Finish();
 107     AliFlowCommonHistResults* res2=sp2->GetCommonHistsResults2nd();
 108     TH1D *h2=res2->GetHistDiffFlowPtPOI();
 109     sprintf(name,"v2QC2_%s_%i_%i","AllCharged",cmin[i],cmax[i]);
 110     h2->SetName(name);
 111     fo->cd();
 112     h2->Write();
 113     for(Int_t j=0;j < 2;j++){
 114       for(Int_t k=0;k < 4;k++){
 115         TList* list = dynamic_cast<TList*>(directory->Get(listTemp->At(i*9+j*4+k+1)->GetName()));
 116         AliFlowAnalysisWithQCumulants* sp = new AliFlowAnalysisWithQCumulants();
 117         sp->GetOutputHistograms(list);
 118         sp->Finish();
 119         AliFlowCommonHistResults* res=sp->GetCommonHistsResults2nd();
 120         TH1D *h=res->GetHistDiffFlowPtPOI();
 121         sprintf(name,"v2QC2_%s_%i_%i%s",spe[k],cmin[i],cmax[i],tech[j]);
 122         h->SetName(name);
 123         fo->cd();
 124         h->Write();
 125       }
 126     }
 127   }
 128   fo->Close();
 129 }
 130
 131
 132