PWGCF/Correlations/macros/AddTaskLV0Correlation.C

   1 AliLeadingV0Correlation* AddTaskLV0Correlation(TString  fListName                = "LV0Correlation",
   2                                                                                            TString  fCollisiontype                       = "PP",
   3                                                                                            Bool_t       fAnalysisMC              = 1,
   4                                                                                            Int_t    fCase                    = 2,
   5                                                                                            Bool_t   fRemoveAutoCorr          = 0,
   6                                                                                            Double_t fPVzCut                  = 10,
   7                                                                                            Int_t    fFilterBit               = 256,
   8                                                                                            Double_t fRapidityCut             = 0.75,
   9                                                                                            Int_t        fmaxEventsinPool         = 2000,
  10                                                                                            Int_t        fminTracksinPool         = 100,
  11                                                                                            Int_t    fMinEventsToMix          = 5,
  12                                                                                            Double_t fV0radius                        = 0.3,
  13                                                                                            Double_t fV0PostoPVz                      = 0.060,
  14                                                                                            Double_t fV0NegtoPVz                      = 0.060,
  15                                                                                            Double_t fDCAV0Daughters                  = 1.00,
  16                                                                                            Double_t fCPAK0                                   = 0.960,
  17                                                                                            Double_t fCPALam                                  = 0.997,
  18                                                                                            Double_t fRejectLamK0                     = 0.005,
  19                                                                                            Double_t fRejectK0Lam                     = 0.010,
  20                                                                                            Double_t fSigmaPID                = 3.0,
  21                                                                                            Double_t fCutCTK0                             = 20.0,
  22                                                                                            Double_t fCutCTLa                             = 30.0,
  23                                                                                            Double_t fMassCutK0               = 0.0105,
  24                                                                                            Double_t fMassCutLa               = 0.0105,
  25                                                                                            Double_t fTriglow                         = 6.0,
  26                                                                                            Double_t fTrighigh                = 12.0,
  27                                                                                            Double_t fTPCClusters                         = 70,
  28                                                                                            Double_t fTPCfindratio            = 0.8)
  29 {
  30         // Get the current analysis manager.
  31     AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  32     if (!mgr) {Error("AddTaskLV0Correlation.C", "No Analysis Manager ");return 0;}
  33
  34         //PVz Binning for pool PP or PbPb
  35         Double_t pvzbinlimits[] = {-12,-11,-10,-9,-8,-7,-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12};
  36         Int_t pvzbinnumb = sizeof(pvzbinlimits)/sizeof(Double_t) - 1;
  37
  38         //Mult Binning for pool pp
  39         Double_t cent_mult_binlimitsPP[] = {  0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
  40                 100,110,120,130,140,150,160,170,180,190,
  41                 200,210,220,230,240,250,260,270,280,290,
  42                 300,500,1000,2000};
  43
  44         Int_t cent_mult_bin_numbPP = sizeof(cent_mult_binlimitsPP)/sizeof(Double_t) - 1;
  45
  46         //Cent Binning for pool PbPb
  47         Double_t cent_mult_binlimitsPbPb[] = { 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,
  48                 14,15,16,17,18,19,20,21,22,23,24,25,26,27,
  49                 28,29,30,31,32,33,34,35,36,37,38,39,40,42,
  50                 44,46,48,50,52,54,56,58,60,65,70,75,80,90};
  51
  52         Int_t cent_mult_bin_numbPbPb = sizeof(cent_mult_binlimitsPbPb)/sizeof(Double_t) - 1;
  53         //Correlation task
  54         AliLeadingV0Correlation *myTask = new AliLeadingV0Correlation(fListName.Data());
  55
  56         myTask->SetCollidingSystem(fCollisiontype);
  57         myTask->SetMCAnalysis(fAnalysisMC);
  58         myTask->SetCase(fCase);
  59         myTask->SetRemoveAutoCorr(fRemoveAutoCorr);
  60         myTask->SetMaxNEventsInPool(fmaxEventsinPool);
  61         myTask->SetMinNTracksInPool(fminTracksinPool);
  62         myTask->SetMinEventsToMix(fMinEventsToMix);
  63         myTask->SetPoolPVzBinLimits(pvzbinnumb,pvzbinlimits);
  64         if(fCollisiontype=="PP")myTask->SetPoolCentBinLimits(cent_mult_bin_numbPP,cent_mult_binlimitsPP);
  65         if(fCollisiontype=="PbPb")myTask->SetPoolCentBinLimits(cent_mult_bin_numbPbPb,cent_mult_binlimitsPbPb);
  66         myTask->SetPrimeryVertexCut(fPVzCut);
  67         myTask->SetFilterBit(fFilterBit);
  68         myTask->SetCutRap(fRapidityCut);
  69         myTask->SetV0Radius(fV0radius);
  70         myTask->SetV0PostoPVz(fV0PostoPVz);
  71         myTask->SetV0NegtoPVz(fV0NegtoPVz);
  72         myTask->SetDCAV0Daughters(fDCAV0Daughters);
  73         myTask->SetCPAK0(fCPAK0);
  74         myTask->SetCPALam(fCPALam);
  75         myTask->SetRejectLamK0(fRejectLamK0);
  76         myTask->SetRejectK0Lam(fRejectK0Lam);
  77         myTask->SetSigmaPID(fSigmaPID);
  78         myTask->SetCTK0(fCutCTK0);
  79         myTask->SetCTLa(fCutCTLa);
  80         myTask->SetMassCutK0(fMassCutK0);
  81         myTask->SetMassCutLa(fMassCutLa);
  82         myTask->SetTrigLow(fTriglow);
  83         myTask->SetTrigHigh(fTrighigh);
  84         myTask->SetTPCClusters(fTPCClusters);
  85         myTask->SetTPCFinables(fTPCfindratio);
  86
  87         mgr->AddTask(myTask);
  88
  89         // Create containers for input/output
  90         TString outputFileName = AliAnalysisManager::GetCommonFileName();
  91         AliAnalysisDataContainer *cinput  = mgr->GetCommonInputContainer();
  92         AliAnalysisDataContainer *coutput = mgr->CreateContainer(fListName.Data(), TList::Class(),AliAnalysisManager::kOutputContainer,outputFileName);
  93
  94         //____________________________________________//
  95         mgr->ConnectInput(myTask,0,cinput);
  96         mgr->ConnectOutput(myTask,1,coutput);
  97
  98 return myTask;
  99 }
 100
 101
 102