PWGJE/macros/AddTaskPWG4HighPtTrackQA.C

   1 void AddTaskPWG4HighPtTrackQA(TString year = "2010", TString prodType = "LHC10h",Bool_t isPbPb = kTRUE, Int_t iAODanalysis = 0, Bool_t bReduced = kTRUE, Int_t filterBit = 272)
   2 {
   3
   4   if(iAODanalysis==0) { //run on ESDs
   5     if(year.Contains("2010")) {
   6       if(bReduced)
   7         AddTaskPWG4HighPtTrackQAAllReduced(prodType.Data(),isPbPb,iAODanalysis);
   8       else
   9         AddTaskPWG4HighPtTrackQAAll(prodType.Data(),isPbPb,iAODanalysis);
  10     }
  11     else if(year.Contains("2011")) {
  12       if(bReduced)
  13         AddTaskPWG4HighPtTrackQAAllReduced2011(prodType.Data(),isPbPb,iAODanalysis);
  14       else
  15         AddTaskPWG4HighPtTrackQAAll2011(prodType.Data(),isPbPb,iAODanalysis);
  16     }
  17     else if(year.Contains("2013")) {
  18       AddTaskPWG4HighPtTrackQApPb();
  19     }
  20
  21   }
  22   else if(iAODanalysis==1) {
  23     AddTaskPWG4HighPtTrackQAAOD(prodType.Data(),isPbPb,iAODanalysis,filterBit);
  24   }
  25
  26 }
  27
  28 void AddTaskPWG4HighPtTrackQApPb(char *prodType = "LHC13b") {
  29
  30   AliPWG4HighPtTrackQA *taskTrackQA05cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,kFALSE,0,10,0,5,AliVEvent::kINT7);
  31   AliPWG4HighPtTrackQA *taskTrackQA75cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,kFALSE,0,10,7,5,AliVEvent::kINT7);
  32
  33   if(!strcmp(prodType,"LHC13d") || !strcmp(prodType,"LHC13e") || !strcmp(prodType,"LHC13f")) {
  34     AliPWG4HighPtTrackQA *taskTrackQA05cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,kFALSE,0,10,0,5,AliVEvent::kEMCEJE);
  35     AliPWG4HighPtTrackQA *taskTrackQA75cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,kFALSE,0,10,7,5,AliVEvent::kEMCEJE);
  36   }
  37
  38 }
  39
  40
  41 void AddTaskPWG4HighPtTrackQAAll(char *prodType = "LHC10h",Bool_t isPbPb=kTRUE, Int_t iAODanalysis = 0)
  42 {
  43
  44   int cent = 10;
  45
  46   AliPWG4HighPtTrackQA *taskTrackQA00cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0);
  47   AliPWG4HighPtTrackQA *taskTrackQA01cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1);
  48   //  AliPWG4HighPtTrackQA *taskTrackQA02cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,2);
  49   // AliPWG4HighPtTrackQA *taskTrackQA10cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,1,0);
  50   // AliPWG4HighPtTrackQA *taskTrackQA11cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,1,1);
  51   //  AliPWG4HighPtTrackQA *taskTrackQA20cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,0);
  52   //  AliPWG4HighPtTrackQA *taskTrackQA21cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,1);
  53   //  AliPWG4HighPtTrackQA *taskTrackQA40cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0);
  54   //  AliPWG4HighPtTrackQA *taskTrackQA41cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,1);
  55   //  AliPWG4HighPtTrackQA *taskTrackQA50cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,5,0);
  56   //  AliPWG4HighPtTrackQA *taskTrackQA60cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,6,0);
  57   AliPWG4HighPtTrackQA *taskTrackQA70cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0);
  58   AliPWG4HighPtTrackQA *taskTrackQA71cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1);
  59   AliPWG4HighPtTrackQA *taskTrackQA72cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2);
  60
  61   if(isPbPb) {
  62     for(cent=0; cent<4; cent++) {
  63       AliPWG4HighPtTrackQA *taskTrackQA00 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0);
  64       AliPWG4HighPtTrackQA *taskTrackQA01 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1);
  65       //    AliPWG4HighPtTrackQA *taskTrackQA02 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,2);
  66       // AliPWG4HighPtTrackQA *taskTrackQA10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,1,0);
  67       // AliPWG4HighPtTrackQA *taskTrackQA11 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,1,1);
  68       //      AliPWG4HighPtTrackQA *taskTrackQA20 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,0);
  69       //      AliPWG4HighPtTrackQA *taskTrackQA21 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,1);
  70       //      AliPWG4HighPtTrackQA *taskTrackQA40 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0);
  71       //      AliPWG4HighPtTrackQA *taskTrackQA41 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,1);
  72       //      AliPWG4HighPtTrackQA *taskTrackQA50 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,5,0);
  73       //      AliPWG4HighPtTrackQA *taskTrackQA60 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,6,0);
  74       AliPWG4HighPtTrackQA *taskTrackQA70 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0);
  75       AliPWG4HighPtTrackQA *taskTrackQA71 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1);
  76       AliPWG4HighPtTrackQA *taskTrackQA72 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2);
  77     }
  78   }
  79
  80 }
  81
  82 void AddTaskPWG4HighPtTrackQAAll2011(char *prodType = "LHC11h",Bool_t isPbPb=kTRUE, Int_t iAODanalysis = 0)
  83 {
  84
  85   Int_t cent = 10;
  86   UInt_t iPhysicsSelectionFlag = AliVEvent::kMB;
  87   UInt_t iPhysicsSelectionFlagCentral = AliVEvent::kCentral;
  88   UInt_t iPhysicsSelectionFlagSemiCentral = AliVEvent::kSemiCentral;
  89
  90
  91   AliPWG4HighPtTrackQA *taskTrackQA00cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlag);
  92   AliPWG4HighPtTrackQA *taskTrackQA01cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlag);
  93   AliPWG4HighPtTrackQA *taskTrackQA70cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlag);
  94   AliPWG4HighPtTrackQA *taskTrackQA71cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlag);
  95   AliPWG4HighPtTrackQA *taskTrackQA72cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlag);
  96
  97   AliPWG4HighPtTrackQA *taskTrackQA05cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlag);
  98   AliPWG4HighPtTrackQA *taskTrackQA74cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlag);
  99   AliPWG4HighPtTrackQA *taskTrackQA75cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlag);
 100
 101   if(isPbPb) {
 102     for(cent=0; cent<4; cent++) {
 103       AliPWG4HighPtTrackQA *taskTrackQA00 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlag);
 104       AliPWG4HighPtTrackQA *taskTrackQA01 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlag);
 105       AliPWG4HighPtTrackQA *taskTrackQA70 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlag);
 106       AliPWG4HighPtTrackQA *taskTrackQA71 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlag);
 107       AliPWG4HighPtTrackQA *taskTrackQA72 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlag);
 108       AliPWG4HighPtTrackQA *taskTrackQA05 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlag);
 109       AliPWG4HighPtTrackQA *taskTrackQA74 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlag);
 110       AliPWG4HighPtTrackQA *taskTrackQA75 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlag);
 111
 112       if(cent==0) {
 113         AliPWG4HighPtTrackQA *taskTrackQA00 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagCentral);
 114         AliPWG4HighPtTrackQA *taskTrackQA01 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlagCentral);
 115         AliPWG4HighPtTrackQA *taskTrackQA70 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlagCentral);
 116         AliPWG4HighPtTrackQA *taskTrackQA71 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlagCentral);
 117         AliPWG4HighPtTrackQA *taskTrackQA72 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlagCentral);
 118         AliPWG4HighPtTrackQA *taskTrackQA05 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagCentral);
 119         AliPWG4HighPtTrackQA *taskTrackQA74 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlagCentral);
 120         AliPWG4HighPtTrackQA *taskTrackQA75 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagCentral);
 121       }
 122       else {
 123         AliPWG4HighPtTrackQA *taskTrackQA00 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagSemiCentral);
 124         AliPWG4HighPtTrackQA *taskTrackQA01 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlagSemiCentral);
 125         AliPWG4HighPtTrackQA *taskTrackQA70 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlagSemiCentral);
 126         AliPWG4HighPtTrackQA *taskTrackQA71 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlagSemiCentral);
 127         AliPWG4HighPtTrackQA *taskTrackQA72 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlagSemiCentral);
 128         AliPWG4HighPtTrackQA *taskTrackQA05 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagSemiCentral);
 129         AliPWG4HighPtTrackQA *taskTrackQA74 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlagSemiCentral);
 130         AliPWG4HighPtTrackQA *taskTrackQA75 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagSemiCentral);
 131       }
 132
 133     }
 134   }
 135
 136 }
 137
 138 void AddTaskPWG4HighPtTrackQAAllReduced(char *prodType = "LHC11h",Bool_t isPbPb=kTRUE, Int_t iAODanalysis = 0)
 139 {
 140
 141   int cent = 10;
 142
 143   if(isPbPb) {
 144     for(cent=0; cent<4; cent++) {
 145       AliPWG4HighPtTrackQA *taskTrackQA00 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0);
 146       AliPWG4HighPtTrackQA *taskTrackQA01 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1);
 147       AliPWG4HighPtTrackQA *taskTrackQA70 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0);
 148       AliPWG4HighPtTrackQA *taskTrackQA71 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1);
 149       AliPWG4HighPtTrackQA *taskTrackQA72 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2);
 150     }
 151   }
 152
 153 }
 154
 155
 156 void AddTaskPWG4HighPtTrackQALHC11hLTS(char *prodType = "LHC10h",Bool_t isPbPb=kTRUE, Int_t iAODanalysis = 0)
 157 {
 158
 159   Int_t cent = 10;
 160   UInt_t iPhysicsSelectionFlag = AliVEvent::kMB|AliVEvent::kCentral|AliVEvent::kSemiCentral;
 161   UInt_t iPhysicsSelectionFlagEMCEJE = AliVEvent::kEMCEJE;
 162
 163   AliPWG4HighPtTrackQA *taskTrackQA00cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlag);
 164   AliPWG4HighPtTrackQA *taskTrackQA01cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlag);
 165   AliPWG4HighPtTrackQA *taskTrackQA70cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlag);
 166   AliPWG4HighPtTrackQA *taskTrackQA71cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlag);
 167   AliPWG4HighPtTrackQA *taskTrackQA72cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlag);
 168   AliPWG4HighPtTrackQA *taskTrackQA05cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlag);
 169   AliPWG4HighPtTrackQA *taskTrackQA74cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlag);
 170   AliPWG4HighPtTrackQA *taskTrackQA75cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlag);
 171   AliPWG4HighPtTrackQA *taskTrackQA40cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0,iPhysicsSelectionFlag);
 172
 173   AliPWG4HighPtTrackQA *taskTrackQAEMCJE00cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagEMCEJE);
 174   AliPWG4HighPtTrackQA *taskTrackQAEMCJE01cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlagEMCEJE);
 175   AliPWG4HighPtTrackQA *taskTrackQAEMCJE70cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlagEMCEJE);
 176   AliPWG4HighPtTrackQA *taskTrackQAEMCJE71cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlagEMCEJE);
 177   AliPWG4HighPtTrackQA *taskTrackQAEMCJE72cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlagEMCEJE);
 178   AliPWG4HighPtTrackQA *taskTrackQAEMCJE05cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagEMCEJE);
 179   AliPWG4HighPtTrackQA *taskTrackQAEMCJE74cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlagEMCEJE);
 180   AliPWG4HighPtTrackQA *taskTrackQAEMCJE75cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagEMCEJE);
 181   AliPWG4HighPtTrackQA *taskTrackQAEMCJE40cent10 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0,iPhysicsSelectionFlagEMCEJE);
 182
 183   if(isPbPb) {
 184     for(cent=0; cent<4; cent++) {
 185
 186       AliPWG4HighPtTrackQA *taskTrackQA00 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlag);
 187       AliPWG4HighPtTrackQA *taskTrackQA01 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlag);
 188       AliPWG4HighPtTrackQA *taskTrackQA70 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlag);
 189       AliPWG4HighPtTrackQA *taskTrackQA71 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlag);
 190       AliPWG4HighPtTrackQA *taskTrackQA72 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlag);
 191       AliPWG4HighPtTrackQA *taskTrackQA05 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlag);
 192       AliPWG4HighPtTrackQA *taskTrackQA74 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlag);
 193       AliPWG4HighPtTrackQA *taskTrackQA75 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlag);
 194       AliPWG4HighPtTrackQA *taskTrackQA40 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0,iPhysicsSelectionFlag);
 195
 196       AliPWG4HighPtTrackQA *taskTrackQAEMCJE00 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagEMCEJE);
 197       AliPWG4HighPtTrackQA *taskTrackQAEMCJE01 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlagEMCEJE);
 198       AliPWG4HighPtTrackQA *taskTrackQAEMCJE70 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlagEMCEJE);
 199       AliPWG4HighPtTrackQA *taskTrackQAEMCJE71 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlagEMCEJE);
 200       AliPWG4HighPtTrackQA *taskTrackQAEMCJE72 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlagEMCEJE);
 201       AliPWG4HighPtTrackQA *taskTrackQAEMCJE05 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagEMCEJE);
 202       AliPWG4HighPtTrackQA *taskTrackQAEMCJE74 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlagEMCEJE);
 203       AliPWG4HighPtTrackQA *taskTrackQAEMCJE75 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagEMCEJE);
 204       AliPWG4HighPtTrackQA *taskTrackQAEMCJE40 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,4,0,iPhysicsSelectionFlagEMCEJE);
 205
 206     }
 207   }
 208
 209 }
 210
 211
 212
 213
 214 void AddTaskPWG4HighPtTrackQAAllReduced2011(char *prodType = "LHC10h",Bool_t isPbPb=kTRUE, Int_t iAODanalysis = 0)
 215 {
 216
 217   int cent = 10;
 218
 219   UInt_t iPhysicsSelectionFlagCentral = AliVEvent::kCentral;
 220   UInt_t iPhysicsSelectionFlagSemiCentral = AliVEvent::kSemiCentral;
 221
 222   AliPWG4HighPtTrackQA *taskTrackQA00C = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagCentral);
 223   AliPWG4HighPtTrackQA *taskTrackQA01C = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlagCentral);
 224   AliPWG4HighPtTrackQA *taskTrackQA21C = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,1,iPhysicsSelectionFlagCentral);
 225   AliPWG4HighPtTrackQA *taskTrackQA70C = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlagCentral);
 226   AliPWG4HighPtTrackQA *taskTrackQA71C = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlagCentral);
 227   AliPWG4HighPtTrackQA *taskTrackQA72C = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlagCentral);
 228   AliPWG4HighPtTrackQA *taskTrackQA05C = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagCentral);
 229   AliPWG4HighPtTrackQA *taskTrackQA74C = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlagCentral);
 230   AliPWG4HighPtTrackQA *taskTrackQA75C = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagCentral);
 231
 232   AliPWG4HighPtTrackQA *taskTrackQA00SC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagSemiCentral);
 233   AliPWG4HighPtTrackQA *taskTrackQA01SC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,1,iPhysicsSelectionFlagSemiCentral);
 234   AliPWG4HighPtTrackQA *taskTrackQA21SC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,2,1,iPhysicsSelectionFlagSemiCentral);
 235   AliPWG4HighPtTrackQA *taskTrackQA70SC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,0,iPhysicsSelectionFlagSemiCentral);
 236   AliPWG4HighPtTrackQA *taskTrackQA71SC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,1,iPhysicsSelectionFlagSemiCentral);
 237   AliPWG4HighPtTrackQA *taskTrackQA72SC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,2,iPhysicsSelectionFlagSemiCentral);
 238   AliPWG4HighPtTrackQA *taskTrackQA05SC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagSemiCentral);
 239   AliPWG4HighPtTrackQA *taskTrackQA74SC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,4,iPhysicsSelectionFlagSemiCentral);
 240   AliPWG4HighPtTrackQA *taskTrackQA75SC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagSemiCentral);
 241
 242 }
 243
 244 void AddTaskPWG4HighPtTrackQAAOD(char *prodType = "LHC10h",Bool_t isPbPb=kTRUE, Int_t iAODanalysis = 1, Int_t filterBit = 272)
 245 {
 246   UInt_t iPhysicsSelectionFlagMB = AliVEvent::kMB;
 247   UInt_t iPhysicsSelectionFlagCentral = AliVEvent::kCentral;
 248   UInt_t iPhysicsSelectionFlagSemiCentral = AliVEvent::kSemiCentral;
 249   UInt_t iPhysicsSelectionFlagINT7 = AliVEvent::kINT7;
 250   UInt_t iPhysicsSelectionFlagEMCEJE = AliVEvent::kEMCEJE;
 251
 252   Int_t cent = 10;
 253
 254   TString strRunPeriod = TString(prodType);
 255
 256   if(isPbPb) {
 257     for(cent=0; cent<4; cent++) {
 258       AliPWG4HighPtTrackQA *taskTrackQAMB = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagMB);
 259       taskTrackQAMB->SetFilterMask(filterBit);
 260       if(strRunPeriod.Contains("LHC10h"))
 261         taskTrackQAMB->SetIncludeNoITS(kTRUE);
 262       else
 263         taskTrackQAMB->SetIncludeNoITS(kFALSE);
 264     }
 265
 266     cent = 10;
 267
 268     if(strRunPeriod.Contains("LHC11h")) {
 269       AliPWG4HighPtTrackQA *taskTrackQAC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagCentral);
 270       taskTrackQAC->SetFilterMask(filterBit);
 271       taskTrackQAC->SetIncludeNoITS(kFALSE);
 272
 273       AliPWG4HighPtTrackQA *taskTrackQASC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagSemiCentral);
 274       taskTrackQASC->SetFilterMask(filterBit);
 275       taskTrackQASC->SetIncludeNoITS(kFALSE);
 276     }
 277   }
 278   else {
 279     cent = 10;
 280
 281     if(strRunPeriod.Contains("LHC13")) {
 282       AliPWG4HighPtTrackQA *taskTrackQAMB = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagINT7);
 283       taskTrackQAMB->SetFilterMask(768);
 284       taskTrackQAMB->SetIncludeNoITS(kFALSE);
 285
 286       AliPWG4HighPtTrackQA *taskTrackQAMB1 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagINT7);
 287       taskTrackQAMB1->SetFilterMask(256);
 288       taskTrackQAMB1->SetIncludeNoITS(kFALSE);
 289
 290       AliPWG4HighPtTrackQA *taskTrackQAMB2 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagINT7);
 291       taskTrackQAMB2->SetFilterMask(512);
 292       taskTrackQAMB2->SetIncludeNoITS(kFALSE);
 293
 294       if(strRunPeriod.EqualTo("LHC13d") || strRunPeriod.EqualTo("LHC13e") || strRunPeriod.EqualTo("LHC13f") || strRunPeriod.EqualTo("LHC13g")) {
 295         AliPWG4HighPtTrackQA *taskTrackQAEMCEJE = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagEMCEJE);
 296         taskTrackQAEMCEJE->SetFilterMask(768);
 297         taskTrackQAEMCEJE->SetIncludeNoITS(kFALSE);
 298
 299         AliPWG4HighPtTrackQA *taskTrackQAEMCEJE1 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagEMCEJE);
 300         taskTrackQAEMCEJE1->SetFilterMask(256);
 301         taskTrackQAEMCEJE1->SetIncludeNoITS(kFALSE);
 302
 303         AliPWG4HighPtTrackQA *taskTrackQAEMCEJE2 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagEMCEJE);
 304         taskTrackQAEMCEJE2->SetFilterMask(512);
 305         taskTrackQAEMCEJE2->SetIncludeNoITS(kFALSE);
 306       }
 307     }
 308     else {
 309       AliPWG4HighPtTrackQA *taskTrackQAMB = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagMB);
 310       taskTrackQAMB->SetFilterMask(filterBit);
 311       if(strRunPeriod.Contains("LHC10h")) taskTrackQAMB->SetIncludeNoITS(kTRUE);
 312     }
 313   }
 314 }
 315
 316 AliPWG4HighPtTrackQA* ConfigureTaskPWG4HighPtTrackQA(char *prodType = "LHC10e14",Bool_t isPbPb=kTRUE,Int_t iAODanalysis = 0, Int_t centClass = 0, Int_t trackType = 0, Int_t cuts = 0, UInt_t iPhysicsSelectionFlag = AliVEvent::kMB)
 317 {
 318   /*
 319     trackType: 0 = global
 320                1 = TPC stand alone
 321                2 = TPC stand alone constrained to SPD vertex
 322                4 = TPC stand alone constrained to SPD vertex with QA track selection on global tracks
 323                5 = Hybrid tracks: constrained TPConly for which no tight ITS is available
 324                6 = Hybrid tracks: constrained loose global for which no tight ITS is available
 325     cuts:      0 (global) = standard ITSTPC2010 a la RAA analysis
 326                1 (global) = ITSrefit, no SPD requirements -> standard for jet analysis
 327                2 (global) = ITSrefit + no hits in SPD
 328                3 (global) = standard ITS tight cuts with nCrossed rows cut for hybrid tracks
 329                0 (TPC)    = standard TPC + NClusters>70
 330                1 (TPC)    = standard TPC + NClusters>0 --> to study new TPC QA recommendations
 331                0 (hybrid 5) = constrained TPConly for which no tight ITS is available
 332                0 (hybrid 6) = constrained loose global for which no tight ITS is available
 333    */
 334
 335   //Load common track cut class
 336   gROOT->LoadMacro("$ALICE_ROOT/PWGJE/macros/CreateTrackCutsPWGJE.C");
 337
 338   // Creates HighPtTrackQA analysis task and adds it to the analysis manager.
 339
 340   // A. Get the pointer to the existing analysis manager via the static access method.
 341   //==============================================================================
 342   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
 343   if (!mgr) {
 344     Error("AddTaskPWG4HighPtQMC", "No analysis manager to connect to.");
 345     return NULL;
 346   }
 347
 348   // B. Check the analysis type using the event handlers connected to the analysis
 349   //    manager. The availability of MC handler can also be checked here.
 350   //==============================================================================
 351   if (!mgr->GetInputEventHandler()) {
 352     ::Error("AddPWG4TaskHighPtTrackQA", "This task requires an input event handler");
 353     return NULL;
 354   }
 355
 356   // C. Create the task, add it to manager.
 357   //===========================================================================
 358
 359   //CREATE THE  CUTS -----------------------------------------------
 360   //Use AliESDtrackCuts
 361   AliESDtrackCuts *trackCuts = new AliESDtrackCuts("AliESDtrackCuts","Standard Cuts");
 362   AliESDtrackCuts *trackCutsReject = 0x0;
 363   AliESDtrackCuts *trackCutsTPConly = new AliESDtrackCuts("AliESDtrackCutsTPConly","TPC only Cuts");
 364
 365   //Standard Cuts
 366   //Set track cuts for global tracks
 367   if(trackType==0 && cuts==0) {
 368     // tight global tracks - RAA analysis
 369     trackCuts = CreateTrackCutsPWGJE(1000);
 370   }
 371   if(trackType==0 && cuts==1) {
 372     //Cuts global tracks with ITSrefit requirement and SPDrequirement for jet analysis
 373     trackCuts = CreateTrackCutsPWGJE(10001006);
 374   }
 375   if(trackType==0 && cuts==5) {
 376     //Cuts global tracks with ITSrefit requirement and SPDrequirement for jet analysis + NCrossedRowsCut>120 recommended in 2011
 377     trackCuts = CreateTrackCutsPWGJE(10001008);
 378   }
 379
 380   if(trackType==0 && cuts==2) {
 381     //Cuts global tracks with ITSrefit requirement but without SPD
 382     trackCuts = CreateTrackCutsPWGJE(10011006);
 383   }
 384   if(trackType==7 && cuts==0) {
 385     // tight global tracks
 386     trackCuts = CreateTrackCutsPWGJE(10041006);
 387     trackCutsReject = CreateTrackCutsPWGJE(1006);
 388     trackCutsReject->SetEtaRange(-0.9,0.9);
 389     trackCutsReject->SetPtRange(0.15, 1e10);
 390   }
 391   if(trackType==7 && cuts==4) {
 392     // tight global tracks +  NCrossedRowsCut>120 recommended in 2011
 393     trackCuts = CreateTrackCutsPWGJE(10041008);
 394     trackCutsReject = CreateTrackCutsPWGJE(1008);
 395     trackCutsReject->SetEtaRange(-0.9,0.9);
 396     trackCutsReject->SetPtRange(0.15, 1e10);
 397   }
 398   if(trackType==7 && cuts==1) {
 399     // tight global tracks
 400     trackCuts = CreateTrackCutsPWGJE(10011006);
 401   }
 402   if(trackType==7 && cuts==5) {
 403     // tight global tracks  + NCrossedRowsCut>120 recommended in 2011
 404     trackCuts = CreateTrackCutsPWGJE(10011008);
 405   }
 406   if(trackType==7 && cuts==2) {
 407     // no requirements on SPD and ITSrefit failed
 408     trackCuts = CreateTrackCutsPWGJE(10041006);       //no ITSrefit requirement filter 256
 409     trackCutsReject = CreateTrackCutsPWGJE(10001006); //ITSrefit requirement filter 16
 410     trackCutsReject->SetEtaRange(-0.9,0.9);
 411     trackCutsReject->SetPtRange(0.15, 1e10);
 412   }
 413   if(trackType==7 && cuts==6) {
 414     // no requirements on SPD and ITSrefit failed
 415     trackCuts = CreateTrackCutsPWGJE(10041008);       //no ITSrefit requirement filter 256
 416     trackCutsReject = CreateTrackCutsPWGJE(10001008); //ITSrefit requirement filter 16
 417     trackCutsReject->SetEtaRange(-0.9,0.9);
 418     trackCutsReject->SetPtRange(0.15, 1e10);
 419   }
 420
 421   if(trackType==1 && cuts==0) {
 422     //Set track cuts for TPConly tracks
 423     trackCuts = CreateTrackCutsPWGJE(2001);
 424   }
 425   if(trackType==1 && cuts==1) {
 426     //Set track cuts for TPConly tracks
 427     trackCuts = CreateTrackCutsPWGJE(10032001);
 428   }
 429
 430   if(trackType==2 && cuts==0) {
 431     //Set track cuts for TPConly constrained tracks
 432     trackCuts = CreateTrackCutsPWGJE(2001);
 433   }
 434   if(trackType==2 && cuts==1) {
 435     //Set track cuts for TPConly constrained tracks w/o cut on NClusters or NCrossedRows
 436     trackCuts = CreateTrackCutsPWGJE(10032001);
 437   }
 438
 439   if(trackType==4 && cuts==0) {
 440      //       Set track cuts for TPConly constrained tracks
 441     trackCuts = CreateTrackCutsPWGJE(2001);
 442   }
 443   if(trackType==4 && cuts==1) {
 444      //       Set track cuts for TPConly constrained tracks
 445     trackCuts = CreateTrackCutsPWGJE(10032001);
 446   }
 447   if(trackType==5 || trackType==6) {
 448     // tight global tracks
 449     trackCuts = CreateTrackCutsPWGJE(1003);
 450
 451     trackCutsReject = CreateTrackCutsPWGJE(10021003);
 452
 453     trackCutsTPConly = CreateTrackCutsPWGJE(2002);
 454
 455     trackCutsReject->SetEtaRange(-0.9,0.9);
 456     trackCutsReject->SetPtRange(0.15, 1e10);
 457
 458     trackCutsTPConly->SetEtaRange(-0.9,0.9);
 459     trackCutsTPConly->SetPtRange(0.15, 1e10);
 460
 461   }
 462
 463   trackCuts->SetEtaRange(-0.9,0.9);
 464   trackCuts->SetPtRange(0.15, 1e10);
 465
 466   TString trigName = "";
 467   if (iPhysicsSelectionFlag == AliVEvent::kAnyINT)
 468     trigName += "kAnyINT";
 469   else if (iPhysicsSelectionFlag == AliVEvent::kAny)
 470     trigName += "kAny";
 471   else if(iPhysicsSelectionFlag == AliVEvent::kINT7)
 472     trigName += "kINT7";
 473   else if(iPhysicsSelectionFlag == AliVEvent::kMB)
 474     trigName += "kMB";
 475   else if(iPhysicsSelectionFlag == AliVEvent::kEMC7)
 476     trigName += "kEMC7";
 477   else if(iPhysicsSelectionFlag == AliVEvent::kEMCEJE)
 478     trigName += "kEMCEJE";
 479   else if(iPhysicsSelectionFlag == AliVEvent::kEMCEGA)
 480     trigName += "kEMCEGA";
 481
 482
 483   //Create the task
 484   AliPWG4HighPtTrackQA *taskPWG4TrackQA = new AliPWG4HighPtTrackQA(Form("AliPWG4HighPtTrackQACent%dTrack%dCuts%d%s",centClass,trackType,cuts,trigName.Data()));
 485   taskPWG4TrackQA->SetTrackType(trackType);
 486   taskPWG4TrackQA->SetCuts(trackCuts);
 487   taskPWG4TrackQA->SetCutsITSLoose(trackCutsReject);
 488   taskPWG4TrackQA->SetCutsTPConly(trackCutsTPConly);
 489
 490   taskPWG4TrackQA->SetPtMax(100.);
 491
 492   if(iAODanalysis)
 493     taskPWG4TrackQA->SetDataType(AliPWG4HighPtTrackQA::kAOD);
 494   else
 495     taskPWG4TrackQA->SetDataType(AliPWG4HighPtTrackQA::kESD);
 496
 497   if(isPbPb) {
 498     taskPWG4TrackQA->SetIsPbPb(kTRUE);
 499     taskPWG4TrackQA->SetCentralityClass(centClass);
 500   }
 501   //  taskPWG4TrackQA->SetSigmaConstrainedMax(5.);
 502
 503   cout << "iPhysicsSelectionFlag: " << iPhysicsSelectionFlag << endl;
 504   taskPWG4TrackQA->SelectCollisionCandidates(iPhysicsSelectionFlag);
 505
 506
 507
 508
 509   // E. Create ONLY the output containers for the data produced by the task.
 510   // Get and connect other common input/output containers via the manager as below
 511   //==============================================================================
 512   TString outputfile = AliAnalysisManager::GetCommonFileName();
 513   outputfile += Form(":PWG4_HighPtTrackQACent%dTrackType%dCuts%d%s",centClass,trackType,cuts,trigName.Data());
 514
 515   AliAnalysisDataContainer *cout_histQAtrack = 0x0;
 516   TString contName = Form("qa_histsQAtrackCent%dType%dcuts%d%s",centClass,trackType,cuts,trigName.Data());
 517   cout_histQAtrack = mgr->CreateContainer(contName.Data(), TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);
 518
 519   mgr->AddTask(taskPWG4TrackQA);
 520   mgr->ConnectInput(taskPWG4TrackQA,0,mgr->GetCommonInputContainer());
 521   mgr->ConnectOutput(taskPWG4TrackQA,1,cout_histQAtrack);
 522
 523   // Return task pointer at the end
 524   return taskPWG4TrackQA;
 525 }