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     }
 261
 262     cent = 10;
 263
 264     if(strRunPeriod.Contains("LHC11h")) {
 265       AliPWG4HighPtTrackQA *taskTrackQAC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagCentral);
 266       taskTrackQAC->SetFilterMask(filterBit);
 267
 268       AliPWG4HighPtTrackQA *taskTrackQASC = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagSemiCentral);
 269       taskTrackQASC->SetFilterMask(filterBit);
 270     }
 271   }
 272   else {
 273     cent = 10;
 274
 275     if(strRunPeriod.Contains("LHC13")) {
 276       AliPWG4HighPtTrackQA *taskTrackQAMB = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagINT7);
 277       taskTrackQAMB->SetFilterMask(768);
 278
 279       AliPWG4HighPtTrackQA *taskTrackQAMB1 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagINT7);
 280       taskTrackQAMB1->SetFilterMask(256);
 281
 282       AliPWG4HighPtTrackQA *taskTrackQAMB2 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagINT7);
 283       taskTrackQAMB2->SetFilterMask(512);
 284
 285       if(strRunPeriod.EqualTo("LHC13d") || strRunPeriod.EqualTo("LHC13e") || strRunPeriod.EqualTo("LHC13f")) {
 286         AliPWG4HighPtTrackQA *taskTrackQAEMCEJE = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagEMCEJE);
 287         taskTrackQAEMCEJE->SetFilterMask(768);
 288
 289         AliPWG4HighPtTrackQA *taskTrackQAEMCEJE1 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,5,iPhysicsSelectionFlagEMCEJE);
 290         taskTrackQAEMCEJE1->SetFilterMask(256);
 291
 292         AliPWG4HighPtTrackQA *taskTrackQAEMCEJE2 = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,7,5,iPhysicsSelectionFlagEMCEJE);
 293         taskTrackQAEMCEJE2->SetFilterMask(512);
 294       }
 295     }
 296     else {
 297       AliPWG4HighPtTrackQA *taskTrackQAMB = ConfigureTaskPWG4HighPtTrackQA(prodType,isPbPb,iAODanalysis,cent,0,0,iPhysicsSelectionFlagMB);
 298       taskTrackQAMB->SetFilterMask(filterBit);
 299     }
 300   }
 301 }
 302
 303 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)
 304 {
 305   /*
 306     trackType: 0 = global
 307                1 = TPC stand alone
 308                2 = TPC stand alone constrained to SPD vertex
 309                4 = TPC stand alone constrained to SPD vertex with QA track selection on global tracks
 310                5 = Hybrid tracks: constrained TPConly for which no tight ITS is available
 311                6 = Hybrid tracks: constrained loose global for which no tight ITS is available
 312     cuts:      0 (global) = standard ITSTPC2010 a la RAA analysis
 313                1 (global) = ITSrefit, no SPD requirements -> standard for jet analysis
 314                2 (global) = ITSrefit + no hits in SPD
 315                3 (global) = standard ITS tight cuts with nCrossed rows cut for hybrid tracks
 316                0 (TPC)    = standard TPC + NClusters>70
 317                1 (TPC)    = standard TPC + NClusters>0 --> to study new TPC QA recommendations
 318                0 (hybrid 5) = constrained TPConly for which no tight ITS is available
 319                0 (hybrid 6) = constrained loose global for which no tight ITS is available
 320    */
 321
 322   //Load common track cut class
 323   gROOT->LoadMacro("$ALICE_ROOT/PWGJE/macros/CreateTrackCutsPWGJE.C");
 324
 325   // Creates HighPtTrackQA analysis task and adds it to the analysis manager.
 326
 327   // A. Get the pointer to the existing analysis manager via the static access method.
 328   //==============================================================================
 329   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
 330   if (!mgr) {
 331     Error("AddTaskPWG4HighPtQMC", "No analysis manager to connect to.");
 332     return NULL;
 333   }
 334
 335   // B. Check the analysis type using the event handlers connected to the analysis
 336   //    manager. The availability of MC handler can also be checked here.
 337   //==============================================================================
 338   if (!mgr->GetInputEventHandler()) {
 339     ::Error("AddPWG4TaskHighPtTrackQA", "This task requires an input event handler");
 340     return NULL;
 341   }
 342
 343   // C. Create the task, add it to manager.
 344   //===========================================================================
 345
 346   //CREATE THE  CUTS -----------------------------------------------
 347   //Use AliESDtrackCuts
 348   AliESDtrackCuts *trackCuts = new AliESDtrackCuts("AliESDtrackCuts","Standard Cuts");
 349   AliESDtrackCuts *trackCutsReject = 0x0;
 350   AliESDtrackCuts *trackCutsTPConly = new AliESDtrackCuts("AliESDtrackCutsTPConly","TPC only Cuts");
 351
 352   //Standard Cuts
 353   //Set track cuts for global tracks
 354   if(trackType==0 && cuts==0) {
 355     // tight global tracks - RAA analysis
 356     trackCuts = CreateTrackCutsPWGJE(1000);
 357   }
 358   if(trackType==0 && cuts==1) {
 359     //Cuts global tracks with ITSrefit requirement and SPDrequirement for jet analysis
 360     trackCuts = CreateTrackCutsPWGJE(10001006);
 361   }
 362   if(trackType==0 && cuts==5) {
 363     //Cuts global tracks with ITSrefit requirement and SPDrequirement for jet analysis + NCrossedRowsCut>120 recommended in 2011
 364     trackCuts = CreateTrackCutsPWGJE(10001008);
 365   }
 366
 367   if(trackType==0 && cuts==2) {
 368     //Cuts global tracks with ITSrefit requirement but without SPD
 369     trackCuts = CreateTrackCutsPWGJE(10011006);
 370   }
 371   if(trackType==7 && cuts==0) {
 372     // tight global tracks
 373     trackCuts = CreateTrackCutsPWGJE(10041006);
 374     trackCutsReject = CreateTrackCutsPWGJE(1006);
 375     trackCutsReject->SetEtaRange(-0.9,0.9);
 376     trackCutsReject->SetPtRange(0.15, 1e10);
 377   }
 378   if(trackType==7 && cuts==4) {
 379     // tight global tracks +  NCrossedRowsCut>120 recommended in 2011
 380     trackCuts = CreateTrackCutsPWGJE(10041008);
 381     trackCutsReject = CreateTrackCutsPWGJE(1008);
 382     trackCutsReject->SetEtaRange(-0.9,0.9);
 383     trackCutsReject->SetPtRange(0.15, 1e10);
 384   }
 385   if(trackType==7 && cuts==1) {
 386     // tight global tracks
 387     trackCuts = CreateTrackCutsPWGJE(10011006);
 388   }
 389   if(trackType==7 && cuts==5) {
 390     // tight global tracks  + NCrossedRowsCut>120 recommended in 2011
 391     trackCuts = CreateTrackCutsPWGJE(10011008);
 392   }
 393   if(trackType==7 && cuts==2) {
 394     // no requirements on SPD and ITSrefit failed
 395     trackCuts = CreateTrackCutsPWGJE(10041006);       //no ITSrefit requirement filter 256
 396     trackCutsReject = CreateTrackCutsPWGJE(10001006); //ITSrefit requirement filter 16
 397     trackCutsReject->SetEtaRange(-0.9,0.9);
 398     trackCutsReject->SetPtRange(0.15, 1e10);
 399   }
 400   if(trackType==7 && cuts==6) {
 401     // no requirements on SPD and ITSrefit failed
 402     trackCuts = CreateTrackCutsPWGJE(10041008);       //no ITSrefit requirement filter 256
 403     trackCutsReject = CreateTrackCutsPWGJE(10001008); //ITSrefit requirement filter 16
 404     trackCutsReject->SetEtaRange(-0.9,0.9);
 405     trackCutsReject->SetPtRange(0.15, 1e10);
 406   }
 407
 408   if(trackType==1 && cuts==0) {
 409     //Set track cuts for TPConly tracks
 410     trackCuts = CreateTrackCutsPWGJE(2001);
 411   }
 412   if(trackType==1 && cuts==1) {
 413     //Set track cuts for TPConly tracks
 414     trackCuts = CreateTrackCutsPWGJE(10032001);
 415   }
 416
 417   if(trackType==2 && cuts==0) {
 418     //Set track cuts for TPConly constrained tracks
 419     trackCuts = CreateTrackCutsPWGJE(2001);
 420   }
 421   if(trackType==2 && cuts==1) {
 422     //Set track cuts for TPConly constrained tracks w/o cut on NClusters or NCrossedRows
 423     trackCuts = CreateTrackCutsPWGJE(10032001);
 424   }
 425
 426   if(trackType==4 && cuts==0) {
 427      //       Set track cuts for TPConly constrained tracks
 428     trackCuts = CreateTrackCutsPWGJE(2001);
 429   }
 430   if(trackType==4 && cuts==1) {
 431      //       Set track cuts for TPConly constrained tracks
 432     trackCuts = CreateTrackCutsPWGJE(10032001);
 433   }
 434   if(trackType==5 || trackType==6) {
 435     // tight global tracks
 436     trackCuts = CreateTrackCutsPWGJE(1003);
 437
 438     trackCutsReject = CreateTrackCutsPWGJE(10021003);
 439
 440     trackCutsTPConly = CreateTrackCutsPWGJE(2002);
 441
 442     trackCutsReject->SetEtaRange(-0.9,0.9);
 443     trackCutsReject->SetPtRange(0.15, 1e10);
 444
 445     trackCutsTPConly->SetEtaRange(-0.9,0.9);
 446     trackCutsTPConly->SetPtRange(0.15, 1e10);
 447
 448   }
 449
 450   trackCuts->SetEtaRange(-0.9,0.9);
 451   trackCuts->SetPtRange(0.15, 1e10);
 452
 453   TString trigName = "";
 454   if (iPhysicsSelectionFlag == AliVEvent::kAnyINT)
 455     trigName += "kAnyINT";
 456   else if (iPhysicsSelectionFlag == AliVEvent::kAny)
 457     trigName += "kAny";
 458   else if(iPhysicsSelectionFlag == AliVEvent::kINT7)
 459     trigName += "kINT7";
 460   else if(iPhysicsSelectionFlag == AliVEvent::kMB)
 461     trigName += "kMB";
 462   else if(iPhysicsSelectionFlag == AliVEvent::kEMC7)
 463     trigName += "kEMC7";
 464   else if(iPhysicsSelectionFlag == AliVEvent::kEMCEJE)
 465     trigName += "kEMCEJE";
 466   else if(iPhysicsSelectionFlag == AliVEvent::kEMCEGA)
 467     trigName += "kEMCEGA";
 468
 469
 470   //Create the task
 471   AliPWG4HighPtTrackQA *taskPWG4TrackQA = new AliPWG4HighPtTrackQA(Form("AliPWG4HighPtTrackQACent%dTrack%dCuts%d%s",centClass,trackType,cuts,trigName.Data()));
 472   taskPWG4TrackQA->SetTrackType(trackType);
 473   taskPWG4TrackQA->SetCuts(trackCuts);
 474   taskPWG4TrackQA->SetCutsITSLoose(trackCutsReject);
 475   taskPWG4TrackQA->SetCutsTPConly(trackCutsTPConly);
 476
 477   taskPWG4TrackQA->SetPtMax(100.);
 478
 479   if(iAODanalysis)
 480     taskPWG4TrackQA->SetDataType(AliPWG4HighPtTrackQA::kAOD);
 481   else
 482     taskPWG4TrackQA->SetDataType(AliPWG4HighPtTrackQA::kESD);
 483
 484   if(isPbPb) {
 485     taskPWG4TrackQA->SetIsPbPb(kTRUE);
 486     taskPWG4TrackQA->SetCentralityClass(centClass);
 487   }
 488   //  taskPWG4TrackQA->SetSigmaConstrainedMax(5.);
 489
 490   cout << "iPhysicsSelectionFlag: " << iPhysicsSelectionFlag << endl;
 491   taskPWG4TrackQA->SelectCollisionCandidates(iPhysicsSelectionFlag);
 492
 493
 494
 495
 496   // E. Create ONLY the output containers for the data produced by the task.
 497   // Get and connect other common input/output containers via the manager as below
 498   //==============================================================================
 499   TString outputfile = AliAnalysisManager::GetCommonFileName();
 500   outputfile += Form(":PWG4_HighPtTrackQACent%dTrackType%dCuts%d%s",centClass,trackType,cuts,trigName.Data());
 501
 502   AliAnalysisDataContainer *cout_histQAtrack = 0x0;
 503   TString contName = Form("qa_histsQAtrackCent%dType%dcuts%d%s",centClass,trackType,cuts,trigName.Data());
 504   cout_histQAtrack = mgr->CreateContainer(contName.Data(), TList::Class(), AliAnalysisManager::kOutputContainer,outputfile);
 505
 506   mgr->AddTask(taskPWG4TrackQA);
 507   mgr->ConnectInput(taskPWG4TrackQA,0,mgr->GetCommonInputContainer());
 508   mgr->ConnectOutput(taskPWG4TrackQA,1,cout_histQAtrack);
 509
 510   // Return task pointer at the end
 511   return taskPWG4TrackQA;
 512 }