PWGHF/hfe/macros/AddTaskHFEtpctofv2.C

   1 AliAnalysisTask *AddTaskHFEtpctofv2(UInt_t trigger=131073,Int_t aodfilter=16,Int_t tpcCls=110, Double_t tpcClsr=50, Int_t tpcClspid=60, Double_t tpcsharedfraction=10, Bool_t rejectkinkmother, Int_t itsCls=4, Double_t chi2peritscl=36, Int_t pixellayer=2, Double_t dcaxy=100,Double_t dcaz=200, Double_t tofsig=30., Double_t tpceff=50., Int_t vzero=1,Int_t debuglevel=5,Double_t etarange=80,Double_t ITSclustersback=0,Double_t minTPCback=-2.0,Double_t maxTPCback=5.0){
   2
   3   //
   4   // Define TPC cut for 2011 data
   5   //
   6   Double_t tpcdedx[8] = {0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0}; //0-5,5-10,10-20,20-30,30-40,40-50,50-60,60-100
   7   // -0.2 0-5%
   8   // -0.15 5-10%
   9   // -0.1 10-20%
  10   // -0.0 20-30%, sigma=1.17
  11   // 0.156 30-40%, sigma=1.2
  12   // 0.19 40-50%, sigma=1.2
  13   // 0.2 50-60%
  14   // 0.2 60-80%
  15   tpcdedx[0]=-0.2;
  16   tpcdedx[1]=-0.15;
  17   tpcdedx[2]=-0.1;
  18   tpcdedx[3]=0.0;
  19   tpcdedx[4]=0.156;
  20   tpcdedx[5]=0.19;
  21   tpcdedx[6]=0.2;
  22   tpcdedx[7]=0.2;
  23   if(TMath::Abs(tpceff-55)<0.01) {
  24     tpcdedx[0]=-0.365;
  25     tpcdedx[1]=-0.314;
  26     tpcdedx[2]=-0.267;
  27     tpcdedx[3]=-0.165;
  28     tpcdedx[4]=-0.022;
  29     tpcdedx[5]= 0.01;
  30     tpcdedx[6]= 0.018;
  31     tpcdedx[7]= 0.018;
  32   }
  33   if(TMath::Abs(tpceff-45)<0.01) {
  34     tpcdedx[0]=-0.062;
  35     tpcdedx[1]=-0.015;
  36     tpcdedx[2]=0.035;
  37     tpcdedx[3]=0.131;
  38     tpcdedx[4]=0.278;
  39     tpcdedx[5]=0.32;
  40     tpcdedx[6]=0.32;
  41     tpcdedx[7]=0.32;
  42   }
  43   if(TMath::Abs(tpceff-60)<0.01) {
  44     tpcdedx[0]=-0.518;
  45     tpcdedx[1]=-0.47;
  46     tpcdedx[2]=-0.42;
  47     tpcdedx[3]=-0.315;
  48     tpcdedx[4]=-0.178;
  49     tpcdedx[5]=-0.145;
  50     tpcdedx[6]=-0.135;
  51     tpcdedx[7]=-0.135;
  52   }
  53   if(TMath::Abs(tpceff-40)<0.01) {
  54     tpcdedx[0]=0.09;
  55     tpcdedx[1]=0.14;
  56     tpcdedx[2]=0.188;
  57     tpcdedx[3]=0.28;
  58     tpcdedx[4]=0.43;
  59     tpcdedx[5]=0.462;
  60     tpcdedx[6]=0.473;
  61     tpcdedx[7]=0.473;
  62   }
  63
  64   // Name
  65   TString appendixx(TString::Format("t%df%dTPC%dr%dp%ds%dkm%dITS%dC%dPi%dDCAr%dz%dTOF%dTPCe%dV%dD%der%di%dt%dt%d",(Int_t)trigger,aodfilter,tpcCls,(Int_t)tpcClsr,tpcClspid,(Int_t) tpcsharedfraction,(Int_t)rejectkinkmother,itsCls,(Int_t) chi2peritscl,(Int_t) pixellayer,(Int_t) dcaxy,(Int_t)dcaz,(Int_t) tofsig,(Int_t)tpceff,vzero,debuglevel,(Int_t)(etarange*0.1),(Int_t)ITSclustersback,(Int_t)(minTPCback*10.0),(Int_t)(maxTPCback*10.0)));
  66   //TString appendixx("tpctofv2");
  67
  68
  69   //set config file name
  70   TString configFile("$ALICE_ROOT/PWGHF/hfe/macros/configs/PbPb/ConfigHFE_FLOW_TOFTPC.C");
  71   //TString configFile("/d/alice12/bailhache/AliRootInstallations/07_10_2012/AliRoot/PWGHF/hfe/macros/configs/PbPb/ConfigHFE_FLOW_TOFTPC.C");
  72   TString checkconfig="ConfigHFE_FLOW_TOFTPC";
  73   if (!gROOT->GetListOfGlobalFunctions()->FindObject(checkconfig.Data()))
  74     gROOT->LoadMacro(configFile.Data());
  75
  76   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  77   AliAnalysisDataContainer *cinput  = mgr->GetCommonInputContainer();
  78   AliAnalysisTaskHFEFlow *task = ConfigHFE_FLOW_TOFTPC(kFALSE,appendixx,trigger,aodfilter,tpcCls, tpcClsr, tpcClspid, tpcsharedfraction,rejectkinkmother,itsCls, chi2peritscl, pixellayer, dcaxy, dcaz,tofsig,&tpcdedx[0],vzero,debuglevel,etarange,kFALSE,ITSclustersback,minTPCback,maxTPCback);
  79
  80
  81   task->SetNbBinsCentralityQCumulant(4);
  82   //task->SetBinCentralityLess(0,0.0);
  83   task->SetBinCentralityLess(0,0.0);
  84   task->SetBinCentralityLess(1,10.0);
  85   task->SetBinCentralityLess(2,20.0);
  86   task->SetBinCentralityLess(3,40.0);
  87   task->SetBinCentralityLess(4,50.0);
  88   //task->SetBinCentralityLess(5,60.0);
  89   //task->SetBinCentralityLess(7,80.0);
  90
  91   task->SetHFEVZEROEventPlane(0x0);
  92   task->SelectCollisionCandidates(AliVEvent::kSemiCentral);
  93   //AliLog::SetClassDebugLevel("AliAnalysisTaskHFEFlow",3);
  94
  95   mgr->AddTask(task);
  96
  97   TString containerName = mgr->GetCommonFileName();
  98   containerName += ":";
  99   containerName += appendixx.Data();
 100
 101   AliAnalysisDataContainer *cinput  = mgr->GetCommonInputContainer();
 102   mgr->ConnectOutput(task,1, mgr->CreateContainer(Form("list_%s",appendixx.Data()), TList::Class(),AliAnalysisManager::kOutputContainer,containerName.Data()));
 103   mgr->ConnectInput(task,0, cinput );
 104
 105   return NULL;
 106
 107
 108 }