PWG1/macros/AddTaskImpParRes.C

   1 AliAnalysisTaskSEImpParRes *AddTaskImpParRes(Bool_t readMC=kFALSE,
   2                                              Int_t selPdg=-1,
   3                                              Bool_t diamond=kTRUE,
   4                                              Bool_t skipTrack=kTRUE,
   5                                              Int_t minmult=0,
   6                                              Int_t maxmult=1000000)
   7 {
   8   //
   9   // Configuration for the study of the impact parameter resolution
  10   //
  11   // xianbao.yuan@pd.infn.it, andrea.dainese@pd.infn.it
  12   //
  13
  14
  15   // Get the pointer to the existing analysis manager via the static access method.
  16   //==============================================================================
  17   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  18   if (!mgr) {
  19     ::Error("AddTaskImpParRes", "No analysis manager to connect to.");
  20     return NULL;
  21   }
  22
  23   // Aanalysis task
  24   AliAnalysisTaskSEImpParRes *d0ResTask = new AliAnalysisTaskSEImpParRes("d0ResAnalysis");
  25   d0ResTask->SetDebugLevel(2);
  26   d0ResTask->SetReadMC(readMC);
  27   d0ResTask->SetSelectedPdg(selPdg);
  28   d0ResTask->SetUseDiamond(diamond);
  29   d0ResTask->SetSkipTrack(skipTrack);
  30   d0ResTask->SetMultiplicityRange(minmult,maxmult);
  31   mgr->AddTask(d0ResTask);
  32
  33   TString fname=Form("%s:ImpParRes_Performance",mgr->GetCommonFileName());
  34   if(selPdg>0) {fname+=selPdg;}
  35
  36
  37   //
  38   // Create containers for input/output
  39   AliAnalysisDataContainer *cinputd0distr = mgr->CreateContainer("cinputd0distr",TChain::Class(),
  40                                                                  AliAnalysisManager::kInputContainer);
  41
  42   AliAnalysisDataContainer *coutputd0ITSpureSARec = mgr->CreateContainer(Form("coutputd0ITSpureSARec_%d_%d",minmult,maxmult),TList::Class(),
  43                                                                    AliAnalysisManager::kOutputContainer,
  44                                                                    fname.Data());
  45
  46   AliAnalysisDataContainer *coutputd0ITSpureSASkip = mgr->CreateContainer(Form("coutputd0ITSpureSASkip_%d_%d",minmult,maxmult),TList::Class(),
  47                                                                    AliAnalysisManager::kOutputContainer,
  48                                                                    fname.Data());
  49
  50   AliAnalysisDataContainer *coutputd0allPointRec = mgr->CreateContainer(Form("coutputd0allPointRec_%d_%d",minmult,maxmult),TList::Class(),
  51                                                                    AliAnalysisManager::kOutputContainer,
  52                                                                    fname.Data());
  53
  54   AliAnalysisDataContainer *coutputd0allPointSkip = mgr->CreateContainer(Form("coutputd0allPointSkip_%d_%d",minmult,maxmult),TList::Class(),
  55                                                                    AliAnalysisManager::kOutputContainer,
  56                                                                    fname.Data());
  57
  58   AliAnalysisDataContainer *coutputd0partPointRec = mgr->CreateContainer(Form("coutputd0partPointRec_%d_%d",minmult,maxmult),TList::Class(),
  59                                                                         AliAnalysisManager::kOutputContainer,
  60                                                                         fname.Data());
  61
  62   AliAnalysisDataContainer *coutputd0partPointSkip = mgr->CreateContainer(Form("coutputd0partPointSkip_%d_%d",minmult,maxmult),TList::Class(),
  63                                                                         AliAnalysisManager::kOutputContainer,
  64                                                                         fname.Data());
  65
  66   AliAnalysisDataContainer *coutputd0onepointSPDRec = mgr->CreateContainer(Form("coutputd0onepointSPDRec_%d_%d",minmult,maxmult),TList::Class(),
  67                                                                         AliAnalysisManager::kOutputContainer,
  68                                                                         fname.Data());
  69
  70   AliAnalysisDataContainer *coutputd0onepointSPDSkip = mgr->CreateContainer(Form("coutputd0onepointSPDSkip_%d_%d",minmult,maxmult),TList::Class(),
  71                                                                         AliAnalysisManager::kOutputContainer,
  72                                                                         fname.Data());
  73
  74   AliAnalysisDataContainer *coutputd0postvTracRec = mgr->CreateContainer(Form("coutputd0postvTracRec_%d_%d",minmult,maxmult),TList::Class(),
  75                                                                         AliAnalysisManager::kOutputContainer,
  76                                                                         fname.Data());
  77
  78
  79   AliAnalysisDataContainer *coutputd0postvTracSkip = mgr->CreateContainer(Form("coutputd0postvTracSkip_%d_%d",minmult,maxmult),TList::Class(),
  80                                                                         AliAnalysisManager::kOutputContainer,
  81                                                                         fname.Data());
  82
  83   AliAnalysisDataContainer *coutputd0negtvTracRec = mgr->CreateContainer(Form("coutputd0negtvTracRec_%d_%d",minmult,maxmult),TList::Class(),
  84                                                                         AliAnalysisManager::kOutputContainer,
  85                                                                         fname.Data());
  86
  87   AliAnalysisDataContainer *coutputd0negtvTracSkip = mgr->CreateContainer(Form("coutputd0negtvTracSkip_%d_%d",minmult,maxmult),TList::Class(),
  88                                                                         AliAnalysisManager::kOutputContainer,
  89                                                                         fname.Data());
  90
  91   AliAnalysisDataContainer *coutputd0pullAllpointRec = mgr->CreateContainer(Form("coutputd0pullAllpointRec_%d_%d",minmult,maxmult),TList::Class(),
  92                                                                         AliAnalysisManager::kOutputContainer,
  93                                                                         fname.Data());
  94
  95
  96   AliAnalysisDataContainer *coutputd0pullAllpointSkip = mgr->CreateContainer(Form("coutputd0pullAllpointSkip_%d_%d",minmult,maxmult),TList::Class(),
  97                                                                         AliAnalysisManager::kOutputContainer,
  98                                                                         fname.Data());
  99
 100   AliAnalysisDataContainer *coutputd0onlyRefitRec = mgr->CreateContainer(Form("coutputd0onlyRefitRec_%d_%d",minmult,maxmult),TList::Class(),
 101                                                                         AliAnalysisManager::kOutputContainer,
 102                                                                         fname.Data());
 103
 104   AliAnalysisDataContainer *coutputd0onlyRefitSkip = mgr->CreateContainer(Form("coutputd0onlyRefitSkip_%d_%d",minmult,maxmult),TList::Class(),
 105                                                                         AliAnalysisManager::kOutputContainer,
 106                                                                         fname.Data());
 107
 108
 109   AliAnalysisDataContainer *coutputd0sinThetaRec = mgr->CreateContainer(Form("coutputd0sinThetaRec_%d_%d",minmult,maxmult),TList::Class(),
 110                                                                         AliAnalysisManager::kOutputContainer,
 111                                                                         fname.Data());
 112
 113
 114   AliAnalysisDataContainer *coutputd0sinThetaSkip = mgr->CreateContainer(Form("coutputd0sinThetaSkip_%d_%d",minmult,maxmult),TList::Class(),
 115                                                                         AliAnalysisManager::kOutputContainer,
 116                                                                         fname.Data());
 117
 118
 119  AliAnalysisDataContainer *coutputd0allPointTrue = mgr->CreateContainer(Form("coutputd0allPointTrue_%d_%d",minmult,maxmult),TList::Class(),
 120                                                                         AliAnalysisManager::kOutputContainer,
 121                                                                         fname.Data());
 122
 123  AliAnalysisDataContainer *coutputd0postvTracTrue = mgr->CreateContainer(Form("coutputd0postvTracTrue_%d_%d",minmult,maxmult),TList::Class(),
 124                                                                         AliAnalysisManager::kOutputContainer,
 125                                                                         fname.Data());
 126
 127  AliAnalysisDataContainer *coutputd0negtvTracTrue = mgr->CreateContainer(Form("coutputd0negtvTracTrue_%d_%d",minmult,maxmult),TList::Class(),
 128                                                                         AliAnalysisManager::kOutputContainer,
 129                                                                         fname.Data());
 130
 131  AliAnalysisDataContainer *coutputd0pullAllpointTrue = mgr->CreateContainer(Form("coutputd0pullAllpointTrue_%d_%d",minmult,maxmult),TList::Class(),
 132                                                                         AliAnalysisManager::kOutputContainer,
 133                                                                         fname.Data());
 134
 135
 136  AliAnalysisDataContainer *coutputd0phiAllpointSkip = mgr->CreateContainer(Form("coutputd0phiAllpointSkip_%d_%d",minmult,maxmult),TList::Class(),
 137                                                                         AliAnalysisManager::kOutputContainer,
 138                                                                         fname.Data());
 139
 140  AliAnalysisDataContainer *coutputd0phiPostvtracSkip = mgr->CreateContainer(Form("coutputd0phiPostvtracSkip_%d_%d",minmult,maxmult),TList::Class(),
 141                                                                         AliAnalysisManager::kOutputContainer,
 142                                                                         fname.Data());
 143
 144  AliAnalysisDataContainer *coutputd0phiNegtvtracSkip = mgr->CreateContainer(Form("coutputd0phiNegtvtracSkip_%d_%d",minmult,maxmult),TList::Class(),
 145                                                                             AliAnalysisManager::kOutputContainer,
 146                                                                           fname.Data());
 147
 148
 149  AliAnalysisDataContainer *coutputd0clusterTypeSPD01Skip = mgr->CreateContainer(Form("coutputd0clusterTypeSPD01Skip_%d_%d",minmult,maxmult),TList::Class(),
 150                                                                         AliAnalysisManager::kOutputContainer,
 151                                                                         fname.Data());
 152
 153  AliAnalysisDataContainer *coutputd0clusterTypeSPD02Skip = mgr->CreateContainer(Form("coutputd0clusterTypeSPD02Skip_%d_%d",minmult,maxmult),TList::Class(),
 154                                                                         AliAnalysisManager::kOutputContainer,
 155                                                                         fname.Data());
 156
 157  AliAnalysisDataContainer *coutputd0clusterTypeSPD03Skip = mgr->CreateContainer(Form("coutputd0clusterTypeSPD03Skip_%d_%d",minmult,maxmult),TList::Class(),
 158                                                                         AliAnalysisManager::kOutputContainer,
 159                                                                         fname.Data());
 160
 161  AliAnalysisDataContainer *coutputd0clusterTypeSPD11Skip = mgr->CreateContainer(Form("coutputd0clusterTypeSPD11Skip_%d_%d",minmult,maxmult),TList::Class(),
 162                                                                         AliAnalysisManager::kOutputContainer,
 163                                                                         fname.Data());
 164
 165  AliAnalysisDataContainer *coutputd0clusterTypeSPD12Skip = mgr->CreateContainer(Form("coutputd0clusterTypeSPD12Skip_%d_%d",minmult,maxmult),TList::Class(),
 166                                                                         AliAnalysisManager::kOutputContainer,
 167                                                                         fname.Data());
 168
 169  AliAnalysisDataContainer *coutputd0clusterTypeSPD13Skip = mgr->CreateContainer(Form("coutputd0clusterTypeSPD13Skip_%d_%d",minmult,maxmult),TList::Class(),
 170                                                                         AliAnalysisManager::kOutputContainer,
 171                                                                         fname.Data());
 172
 173   AliAnalysisDataContainer *coutputd0PID = mgr->CreateContainer(Form("coutputd0PID_%d_%d",minmult,maxmult),TList::Class(),
 174                                                                         AliAnalysisManager::kOutputContainer,
 175                                                                         fname.Data());
 176
 177   AliAnalysisDataContainer *coutputd0Pt = mgr->CreateContainer(Form("coutputd0Pt_%d_%d",minmult,maxmult),TList::Class(),
 178                                                                         AliAnalysisManager::kOutputContainer,
 179                                                                         fname.Data());
 180
 181   AliAnalysisDataContainer *coutputNentries = mgr->CreateContainer(Form("coutputNentries_%d_%d",minmult,maxmult),TH1F::Class(),
 182                                                                      AliAnalysisManager::kOutputContainer,
 183                                                                    fname.Data());
 184
 185   AliAnalysisDataContainer *coutputEstimVtx = mgr->CreateContainer(Form("coutputEstimVtx_%d_%d",minmult,maxmult),TH1F::Class(),
 186                                                                      AliAnalysisManager::kOutputContainer,
 187                                                                    fname.Data());
 188
 189   // Attach input
 190   mgr->ConnectInput(d0ResTask,0,mgr->GetCommonInputContainer());
 191   // Attack output
 192   mgr->ConnectOutput(d0ResTask,1,coutputd0ITSpureSARec);
 193   mgr->ConnectOutput(d0ResTask,2,coutputd0ITSpureSASkip);
 194   mgr->ConnectOutput(d0ResTask,3,coutputd0allPointRec);
 195   mgr->ConnectOutput(d0ResTask,4,coutputd0allPointSkip);
 196   mgr->ConnectOutput(d0ResTask,5,coutputd0partPointRec);
 197   mgr->ConnectOutput(d0ResTask,6,coutputd0partPointSkip);
 198   mgr->ConnectOutput(d0ResTask,7,coutputd0onepointSPDRec);
 199   mgr->ConnectOutput(d0ResTask,8,coutputd0onepointSPDSkip);
 200   mgr->ConnectOutput(d0ResTask,9,coutputd0postvTracRec);
 201   mgr->ConnectOutput(d0ResTask,10,coutputd0postvTracSkip);
 202   mgr->ConnectOutput(d0ResTask,11,coutputd0negtvTracRec);
 203   mgr->ConnectOutput(d0ResTask,12,coutputd0negtvTracSkip);
 204   mgr->ConnectOutput(d0ResTask,13,coutputd0pullAllpointRec);
 205   mgr->ConnectOutput(d0ResTask,14,coutputd0pullAllpointSkip);
 206   mgr->ConnectOutput(d0ResTask,15,coutputd0onlyRefitRec);
 207   mgr->ConnectOutput(d0ResTask,16,coutputd0onlyRefitSkip);
 208   mgr->ConnectOutput(d0ResTask,17,coutputd0sinThetaRec);
 209   mgr->ConnectOutput(d0ResTask,18,coutputd0sinThetaSkip);
 210   mgr->ConnectOutput(d0ResTask,19,coutputd0allPointTrue);
 211   mgr->ConnectOutput(d0ResTask,20,coutputd0postvTracTrue);
 212   mgr->ConnectOutput(d0ResTask,21,coutputd0negtvTracTrue);
 213   mgr->ConnectOutput(d0ResTask,22,coutputd0pullAllpointTrue);
 214   mgr->ConnectOutput(d0ResTask,23,coutputd0phiAllpointSkip);
 215   mgr->ConnectOutput(d0ResTask,24,coutputd0phiPostvtracSkip);
 216   mgr->ConnectOutput(d0ResTask,25,coutputd0phiNegtvtracSkip);
 217   mgr->ConnectOutput(d0ResTask,26,coutputd0clusterTypeSPD01Skip);
 218   mgr->ConnectOutput(d0ResTask,27,coutputd0clusterTypeSPD02Skip);
 219   mgr->ConnectOutput(d0ResTask,28,coutputd0clusterTypeSPD03Skip);
 220   mgr->ConnectOutput(d0ResTask,29,coutputd0clusterTypeSPD11Skip);
 221   mgr->ConnectOutput(d0ResTask,30,coutputd0clusterTypeSPD12Skip);
 222   mgr->ConnectOutput(d0ResTask,31,coutputd0clusterTypeSPD13Skip);
 223   mgr->ConnectOutput(d0ResTask,32,coutputd0PID);
 224   mgr->ConnectOutput(d0ResTask,33,coutputd0Pt);
 225   mgr->ConnectOutput(d0ResTask,34,coutputNentries);
 226   mgr->ConnectOutput(d0ResTask,35,coutputEstimVtx);
 227
 228   return d0ResTask;
 229 }