PWGJE/macros/AddTaskPWG4HighPtSpectra.C

   1 //DEFINITION OF A FEW CONSTANTS
   2 #define round(x) ((x)>=0?(long)((x)+0.5):(long)((x)-0.5))
   3
   4 const Float_t phimin = 0.;
   5 const Float_t phimax = 2.*TMath::Pi();
   6 const Float_t etamin = -0.9;
   7 const Float_t etamax = 0.9;
   8
   9 const Int_t   mintrackrefsTPC = 0;
  10 const Int_t   mintrackrefsITS = 0;
  11
  12 void AddTaskPWG4HighPtSpectra(TString year = "2011", char *prodType = "LHC11h",Bool_t isPbPb=kTRUE, Int_t iAODanalysis = 0, UInt_t triggerMask = AliVEvent::kMB,Bool_t bSelectHijingParticles = kFALSE)
  13 {
  14   int cent = 10;
  15
  16   if(year.Contains("2011")) {
  17     AddTaskPWG4HighPtSpectra2011(prodType,isPbPb,bSelectHijingParticles);
  18   }
  19   else {
  20
  21     AliPWG4HighPtSpectra *taskSpectra00cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,0,triggerMask,bSelectHijingParticles);
  22     AliPWG4HighPtSpectra *taskSpectra01cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,1,triggerMask,bSelectHijingParticles);
  23     //  AliPWG4HighPtSpectra *taskSpectra02cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,2,triggerMask,bSelectHijingParticles);
  24     //  AliPWG4HighPtSpectra *taskSpectra10cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,1,0,triggerMask,bSelectHijingParticles);
  25     //  AliPWG4HighPtSpectra *taskSpectra20cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,triggerMask,bSelectHijingParticles);
  26     AliPWG4HighPtSpectra *taskSpectra70cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,0,triggerMask,bSelectHijingParticles);
  27     AliPWG4HighPtSpectra *taskSpectra71cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,1,triggerMask,bSelectHijingParticles);
  28     AliPWG4HighPtSpectra *taskSpectra72cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,2,triggerMask,bSelectHijingParticles);
  29
  30     if(isPbPb) {
  31       for(cent=0; cent<4; cent++) {
  32         AliPWG4HighPtSpectra *taskSpectra00 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,0,triggerMask,bSelectHijingParticles);
  33         AliPWG4HighPtSpectra *taskSpectra01 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,1,triggerMask,bSelectHijingParticles);
  34         //      AliPWG4HighPtSpectra *taskSpectra02 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,2,triggerMask,bSelectHijingParticles);
  35         //      AliPWG4HighPtSpectra *taskSpectra10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,1,0,triggerMask,bSelectHijingParticles);
  36         //      AliPWG4HighPtSpectra *taskSpectra20 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,triggerMask,bSelectHijingParticles);
  37         AliPWG4HighPtSpectra *taskSpectra70 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,0,triggerMask,bSelectHijingParticles);
  38         AliPWG4HighPtSpectra *taskSpectra71 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,1,triggerMask,bSelectHijingParticles);
  39         AliPWG4HighPtSpectra *taskSpectra72 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,2,triggerMask,bSelectHijingParticles);
  40       }
  41     }
  42   }
  43
  44 }
  45
  46 void AddTaskPWG4HighPtSpectra2011(char *prodType = "LHC10h", Bool_t isPbPb=kTRUE,Bool_t bSelHijingParticles = kFALSE)
  47 {
  48
  49   int cent = 10;
  50
  51   UInt_t iPhysicsSelectionFlag = AliVEvent::kMB|AliVEvent::kCentral|AliVEvent::kSemiCentral;
  52   UInt_t iPhysicsSelectionFlagEMCEJE = AliVEvent::kEMCEJE;
  53
  54   AliPWG4HighPtSpectra *taskSpectra00cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,0,iPhysicsSelectionFlag,bSelHijingParticles);
  55   AliPWG4HighPtSpectra *taskSpectra01cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,1,iPhysicsSelectionFlag,bSelHijingParticles);
  56   AliPWG4HighPtSpectra *taskSpectra70cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,0,iPhysicsSelectionFlag,bSelHijingParticles);
  57   AliPWG4HighPtSpectra *taskSpectra71cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,1,iPhysicsSelectionFlag,bSelHijingParticles);
  58   AliPWG4HighPtSpectra *taskSpectra72cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,2,iPhysicsSelectionFlag,bSelHijingParticles);
  59   AliPWG4HighPtSpectra *taskSpectra05cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,5,iPhysicsSelectionFlag,bSelHijingParticles);
  60   AliPWG4HighPtSpectra *taskSpectra74cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,4,iPhysicsSelectionFlag,bSelHijingParticles);
  61   AliPWG4HighPtSpectra *taskSpectra75cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,5,iPhysicsSelectionFlag,bSelHijingParticles);
  62   AliPWG4HighPtSpectra *taskSpectra20cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,iPhysicsSelectionFlag,bSelHijingParticles);
  63
  64   AliPWG4HighPtSpectra *taskSpectraEMCJE00cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  65   AliPWG4HighPtSpectra *taskSpectraEMCJE01cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,1,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  66   AliPWG4HighPtSpectra *taskSpectraEMCJE70cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  67   AliPWG4HighPtSpectra *taskSpectraEMCJE71cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,1,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  68   AliPWG4HighPtSpectra *taskSpectraEMCJE72cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,2,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  69   AliPWG4HighPtSpectra *taskSpectraEMCJE05cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,5,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  70   AliPWG4HighPtSpectra *taskSpectraEMCJE74cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,4,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  71   AliPWG4HighPtSpectra *taskSpectraEMCJE75cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,5,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  72   AliPWG4HighPtSpectra *taskSpectraEMCJE20cent10 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  73
  74
  75   if(isPbPb) {
  76     for(cent=0; cent<4; cent++) {
  77
  78       AliPWG4HighPtSpectra *taskSpectra00 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,0,iPhysicsSelectionFlag,bSelHijingParticles);
  79       AliPWG4HighPtSpectra *taskSpectra01 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,1,iPhysicsSelectionFlag,bSelHijingParticles);
  80       AliPWG4HighPtSpectra *taskSpectra70 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,0,iPhysicsSelectionFlag,bSelHijingParticles);
  81       AliPWG4HighPtSpectra *taskSpectra71 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,1,iPhysicsSelectionFlag,bSelHijingParticles);
  82       AliPWG4HighPtSpectra *taskSpectra72 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,2,iPhysicsSelectionFlag,bSelHijingParticles);
  83       AliPWG4HighPtSpectra *taskSpectra05 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,5,iPhysicsSelectionFlag,bSelHijingParticles);
  84       AliPWG4HighPtSpectra *taskSpectra74 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,4,iPhysicsSelectionFlag,bSelHijingParticles);
  85       AliPWG4HighPtSpectra *taskSpectra75 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,5,iPhysicsSelectionFlag,bSelHijingParticles);
  86       AliPWG4HighPtSpectra *taskSpectra20 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,iPhysicsSelectionFlag,bSelHijingParticles);
  87
  88       AliPWG4HighPtSpectra *taskSpectraEMCJE00 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  89       AliPWG4HighPtSpectra *taskSpectraEMCJE01 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,1,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  90       AliPWG4HighPtSpectra *taskSpectraEMCJE70 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  91       AliPWG4HighPtSpectra *taskSpectraEMCJE71 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,1,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  92       AliPWG4HighPtSpectra *taskSpectraEMCJE72 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,2,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  93       AliPWG4HighPtSpectra *taskSpectraEMCJE05 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,0,5,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  94       AliPWG4HighPtSpectra *taskSpectraEMCJE74 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,4,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  95       AliPWG4HighPtSpectra *taskSpectraEMCJE75 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,7,5,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  96       AliPWG4HighPtSpectra *taskSpectraEMCJE20 = ConfigureTaskPWG4HighPtSpectra(prodType,isPbPb,cent,2,0,iPhysicsSelectionFlagEMCEJE,bSelHijingParticles);
  97
  98     }
  99
 100   }
 101
 102 }
 103
 104
 105 AliPWG4HighPtSpectra* ConfigureTaskPWG4HighPtSpectra(char *prodType = "LHC10e14", Bool_t isPbPb=kTRUE,Int_t centClass = 0,
 106                                                      Int_t trackType = 0, Int_t cuts = 0, UInt_t triggerMask = AliVEvent::kMB,
 107                                                      Bool_t bSelectHijingParticles = kFALSE)
 108 {
 109
 110   /*
 111    trackType: 0 = global
 112                1 = TPC stand alone
 113                2 = TPC stand alone constrained to SPD vertex
 114     cuts:      0 (global) = standard ITSTPC2010 a la RAA analysis
 115                1 (global) = ITSrefit, no SPD requirements -> standard for jet analysis
 116                2 (global) = ITSrefit + no hits in SPD
 117                3 (global) = standard ITS tight cuts with nCrossed rows cut for hybrid tracks
 118                0 (TPC)    = standard TPC + NClusters>70
 119                1 (TPC)    = standard TPC + NClusters>0 --> to study new TPC QA recommendations
 120    */
 121
 122   // Creates HighPtSpectra analysis task and adds it to the analysis manager.
 123
 124   //Load common track cut class
 125   gROOT->LoadMacro("$ALICE_ROOT/PWGJE/macros/CreateTrackCutsPWGJE.C");
 126
 127   // A. Get the pointer to the existing analysis manager via the static access method.
 128   //==============================================================================
 129   AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
 130   if (!mgr) {
 131     Error("AddTaskPWG4HighPtSpectra", "No analysis manager to connect to.");
 132     return NULL;
 133   }
 134
 135   // B. Check the analysis type using the event handlers connected to the analysis
 136   //    manager. The availability of MC handler can also be checked here.
 137   //==============================================================================
 138   if (!mgr->GetInputEventHandler()) {
 139     ::Error("AddTaskPWG4HighPtSpectra", "This task requires an input event handler");
 140     return NULL;
 141   }
 142   TString type = mgr->GetInputEventHandler()->GetDataType(); // can be "ESD" or "AOD"
 143   const char *analysisType = "ESD";//"TPC"
 144
 145   // C. Create the task, add it to manager.
 146   //===========================================================================
 147  //CONTAINER DEFINITION
 148   Info("AliPWG4HighPtSpectra","SETUP CONTAINER");
 149   //the sensitive variables, their indices
 150   UInt_t ipt   = 0;
 151   UInt_t iphi  = 1;
 152   UInt_t ieta  = 2;
 153   UInt_t incls = 3;
 154
 155   //Setting up the container grid...
 156   UInt_t nstep = 4; //Steps/Modes for containers
 157   Int_t kStepReconstructed          = 0;
 158   Int_t kStepSecondaries            = 1;
 159   Int_t kStepReconstructedMC        = 2;
 160   Int_t kStepMCAcceptance           = 3;
 161
 162   //redefine pt ranges in case of Jet-Jet production
 163   Float_t ptBinEdges[4][4];
 164   Float_t ptmin =  0.0 ;
 165   Float_t ptmax =  100.0 ;
 166
 167   ptBinEdges[0][0] = 2.;
 168   ptBinEdges[0][1] = 0.2;
 169   ptBinEdges[1][0] = 6.;
 170   ptBinEdges[1][1] = 0.5;
 171   ptBinEdges[2][0] = 20.;
 172   ptBinEdges[2][1] = 2.;
 173   ptBinEdges[3][0] = 100.;
 174   ptBinEdges[3][1] = 5.;
 175
 176   const Int_t nvar   = 4; //number of variables on the grid: pt:phi:eta:NClsIter1
 177
 178   const Int_t nbin11 = round((ptBinEdges[0][0]-ptmin)/ptBinEdges[0][1]);
 179   const Int_t nbin12 = round((ptBinEdges[1][0]-ptBinEdges[0][0])/ptBinEdges[1][1])+nbin11;
 180   const Int_t nbin13 = round((ptBinEdges[2][0]-ptBinEdges[1][0])/ptBinEdges[2][1])+nbin12;
 181   const Int_t nbin14 = round((ptBinEdges[3][0]-ptBinEdges[2][0])/ptBinEdges[3][1])+nbin13;
 182
 183   const Int_t nbin1  = nbin14; //bins in pt
 184   const Int_t nbin2  =  36;    //bins in phi
 185   const Int_t nbin3  =  10;     //bins in eta
 186   const Int_t nbin4  =  1;//6;     //bins in NClsIter1: 0 70 80 90 100 120
 187
 188
 189   //arrays for the number of bins in each dimension
 190   Int_t iBin[nvar];
 191   iBin[0]=nbin1;
 192   iBin[1]=nbin2;
 193   iBin[2]=nbin3;
 194   iBin[3]=nbin4;
 195
 196   //arrays for lower bounds :
 197   Double_t *binLim1=new Double_t[nbin1+1];
 198   Double_t *binLim2=new Double_t[nbin2+1];
 199   Double_t *binLim3=new Double_t[nbin3+1];
 200   Double_t *binLim4=new Double_t[nbin4+1];
 201
 202   //values for bin lower bounds
 203   for(Int_t i=0; i<=nbin1; i++) {
 204     if(i<=nbin11) binLim1[i]=(Double_t)ptmin + (ptBinEdges[0][0]-ptmin)/nbin11*(Double_t)i ;
 205     if(i<=nbin12 && i>nbin11) binLim1[i]=(Double_t)ptBinEdges[0][0] + (ptBinEdges[1][0]-ptBinEdges[0][0])/(nbin12-nbin11)*((Double_t)i-(Double_t)nbin11) ;
 206     if(i<=nbin13 && i>nbin12) binLim1[i]=(Double_t)ptBinEdges[1][0] + (ptBinEdges[2][0]-ptBinEdges[1][0])/(nbin13-nbin12)*((Double_t)i-(Double_t)nbin12) ;
 207     if(i<=nbin14 && i>nbin13) binLim1[i]=(Double_t)ptBinEdges[2][0] + (ptBinEdges[3][0]-ptBinEdges[2][0])/(nbin14-nbin13)*((Double_t)i-(Double_t)nbin13) ;
 208   }
 209   for(Int_t i=0; i<=nbin2; i++) binLim2[i]=(Double_t)phimin + (phimax-phimin)/nbin2*(Double_t)i ;
 210   for(Int_t i=0; i<=nbin3; i++) binLim3[i]=(Double_t)etamin + (etamax-etamin)/nbin3*(Double_t)i ;
 211   binLim4[0] = 70.;//0.;
 212   binLim4[1] = 160.;//70.;
 213   // binLim4[2] = 80.;
 214   // binLim4[3] = 90.;
 215   // binLim4[4] = 100.;
 216   // binLim4[5] = 120.;
 217   // binLim4[6] = 160.;
 218
 219
 220   AliCFContainer* containerPos = new AliCFContainer("containerPos","container for positive tracks",nstep,nvar,iBin);
 221   //setting the bin limits
 222   containerPos -> SetBinLimits(ipt,binLim1);
 223   containerPos -> SetBinLimits(iphi,binLim2);
 224   containerPos -> SetBinLimits(ieta,binLim3);
 225   containerPos -> SetBinLimits(incls,binLim4);
 226
 227   AliCFContainer* containerNeg = new AliCFContainer("containerNeg","container for negative tracks",nstep,nvar,iBin);
 228   //setting the bin limits
 229   containerNeg -> SetBinLimits(ipt,binLim1);
 230   containerNeg -> SetBinLimits(iphi,binLim2);
 231   containerNeg -> SetBinLimits(ieta,binLim3);
 232   containerNeg -> SetBinLimits(incls,binLim4);
 233
 234   //CREATE THE  CUTS -----------------------------------------------
 235   //Use AliESDtrackCuts
 236   AliESDtrackCuts *trackCuts = new AliESDtrackCuts("AliESDtrackCuts","Standard Cuts");
 237   AliESDtrackCuts *trackCutsReject = 0x0;
 238   //Standard Cuts
 239   //Set track cuts for global tracks
 240   if(trackType==0 && cuts==0) {
 241     // tight global tracks - RAA analysis
 242     trackCuts = CreateTrackCutsPWGJE(1000);
 243   }
 244   if(trackType==0 && cuts==1) {
 245     //Cuts global tracks with ITSrefit requirement and SPDrequirement for jet analysis
 246     trackCuts = CreateTrackCutsPWGJE(10001006);
 247   }
 248   if(trackType==0 && cuts==5) {
 249     //Cuts global tracks with ITSrefit requirement and SPDrequirement for jet analysis + NCrossedRowsCut>120 recommended in 2011
 250     trackCuts = CreateTrackCutsPWGJE(10001008);
 251   }
 252
 253   if(trackType==0 && cuts==2) {
 254     //Cuts global tracks with ITSrefit requirement but without SPD
 255     if (!strcmp(prodType,"LHC12a15e"))
 256       trackCuts = CreateTrackCutsPWGJE(10011008);
 257     else
 258       trackCuts = CreateTrackCutsPWGJE(10011006);
 259   }
 260   if(trackType==7 && cuts==0) {
 261     // no requirements on SPD and ITSrefit failed
 262     trackCuts = CreateTrackCutsPWGJE(10041006);   //no ITSrefit requirement
 263     trackCutsReject = CreateTrackCutsPWGJE(1006); //ITSrefit requirement
 264     trackCutsReject->SetEtaRange(etamin,etamax);
 265     trackCutsReject->SetPtRange(0.15, 1e10);
 266   }
 267   if(trackType==7 && cuts==4) {
 268     // tight global tracks +  NCrossedRowsCut>120 recommended in 2011
 269     trackCuts = CreateTrackCutsPWGJE(10041008);
 270     trackCutsReject = CreateTrackCutsPWGJE(1008);
 271     trackCutsReject->SetEtaRange(-0.9,0.9);
 272     trackCutsReject->SetPtRange(0.15, 1e10);
 273   }
 274   if(trackType==7 && cuts==1) {
 275     //Cuts global tracks with ITSrefit requirement but without SPD
 276     trackCuts = CreateTrackCutsPWGJE(10011006);
 277   }
 278   if(trackType==7 && cuts==5) {
 279     // tight global tracks  + NCrossedRowsCut>120 recommended in 2011
 280     trackCuts = CreateTrackCutsPWGJE(10011008);
 281   }
 282
 283   if(trackType==7 && cuts==2) {
 284     // no requirements on SPD and ITSrefit failed
 285     trackCuts = CreateTrackCutsPWGJE(10041006);       //no ITSrefit requirement filter 256
 286     trackCutsReject = CreateTrackCutsPWGJE(10001006); //ITSrefit requirement filter 16
 287     trackCutsReject->SetEtaRange(etamin,etamax);
 288     trackCutsReject->SetPtRange(0.15, 1e10);
 289   }
 290   if(trackType==7 && cuts==6) {
 291     // no requirements on SPD and ITSrefit failed
 292     trackCuts = CreateTrackCutsPWGJE(10041008);       //no ITSrefit requirement filter 256
 293     trackCutsReject = CreateTrackCutsPWGJE(10001008); //ITSrefit requirement filter 16
 294     trackCutsReject->SetEtaRange(-0.9,0.9);
 295     trackCutsReject->SetPtRange(0.15, 1e10);
 296   }
 297
 298
 299
 300   if(trackType==1 && cuts==0) {
 301     //Set track cuts for TPConly tracks
 302     trackCuts = CreateTrackCutsPWGJE(2001);
 303   }
 304   if(trackType==2 && cuts==0) {
 305      //       Set track cuts for TPConly constrained tracks
 306     trackCuts = CreateTrackCutsPWGJE(2001);
 307   }
 308   trackCuts->SetEtaRange(etamin,etamax);
 309   trackCuts->SetPtRange(0.15, 1e10);
 310
 311   // Gen-Level kinematic cuts
 312   AliCFTrackKineCuts *mcKineCuts = new AliCFTrackKineCuts("mcKineCuts","MC-level kinematic cuts");
 313   mcKineCuts->SetPtRange(0.15,1e10);
 314   mcKineCuts->SetEtaRange(etamin,etamax);//-0.5,0.5);
 315   mcKineCuts->SetRequireIsCharged(kTRUE);
 316
 317   //Acceptance Cuts
 318   //AliCFAcceptanceCuts *mcAccCuts = new AliCFAcceptanceCuts("mcAccCuts","MC acceptance cuts");
 319   // mcAccCuts->SetMinNHitITS(mintrackrefsITS);
 320   //mcAccCuts->SetMinNHitTPC(mintrackrefsTPC);
 321
 322   TObjArray* recList = new TObjArray(0);
 323   TObjArray* secList = new TObjArray(0) ;
 324   TObjArray* recMCList = new TObjArray(0);
 325
 326   printf("CREATE MC KINE CUTS\n");
 327   TObjArray* mcList = new TObjArray(0) ;
 328   mcList->AddLast(mcKineCuts);
 329   //mcList->AddLast(mcAccCuts);
 330
 331   //CREATE THE INTERFACE TO CORRECTION FRAMEWORK USED IN THE TASK
 332   printf("CREATE INTERFACE AND CUTS\n");
 333   AliCFManager* manPos = new AliCFManager("manPos","Manager for Positive tracks") ;
 334   manPos->SetParticleContainer(containerPos);
 335   manPos->SetParticleCutsList(kStepReconstructed,recList);
 336   manPos->SetParticleCutsList(kStepSecondaries,secList);
 337   manPos->SetParticleCutsList(kStepReconstructedMC,recMCList);
 338   manPos->SetParticleCutsList(kStepMCAcceptance,mcList);
 339
 340   AliCFManager* manNeg = new AliCFManager("manNeg","Manager for Negative tracks") ;
 341   manNeg->SetParticleContainer(containerNeg);
 342   manNeg->SetParticleCutsList(kStepReconstructed,recList);
 343   manNeg->SetParticleCutsList(kStepSecondaries,secList);
 344   manNeg->SetParticleCutsList(kStepReconstructedMC,recMCList);
 345   manNeg->SetParticleCutsList(kStepMCAcceptance,mcList);
 346
 347
 348   AliPWG4HighPtSpectra *taskPWG4HighPtSpectra = new AliPWG4HighPtSpectra(Form("AliPWG4HighPtSpectraCent%dTrackType%dCuts%dPSF%d",centClass,trackType,cuts,triggerMask));
 349   taskPWG4HighPtSpectra->SetTrackType(trackType);
 350   taskPWG4HighPtSpectra->SetCuts(trackCuts);
 351   taskPWG4HighPtSpectra->SetCutsReject(trackCutsReject);
 352   taskPWG4HighPtSpectra->SetCFManagerPos(manPos); //here is set the CF manager +
 353   taskPWG4HighPtSpectra->SetCFManagerNeg(manNeg); //here is set the CF manager -
 354   taskPWG4HighPtSpectra->SetTriggerMask(triggerMask);
 355   taskPWG4HighPtSpectra->SelectHIJINGOnly(bSelectHijingParticles);
 356
 357   if(isPbPb) {
 358     taskPWG4HighPtSpectra->SetIsPbPb(kTRUE);
 359     taskPWG4HighPtSpectra->SetCentralityClass(centClass);
 360   }
 361   //  taskPWG4HighPtSpectra->SetSigmaConstrainedMax(5.);
 362
 363
 364   // E. Create ONLY the output containers for the data produced by the task.
 365   // Get and connect other common input/output containers via the manager as below
 366   //==============================================================================
 367
 368   //------ output containers ------
 369   TString outputfile = AliAnalysisManager::GetCommonFileName();
 370   if(triggerMask == AliVEvent::kMB)
 371     outputfile += Form(":PWG4_HighPtSpectraCent%dTrackType%dCuts%dkMB",centClass,trackType,cuts);
 372   else if(triggerMask == AliVEvent::kCentral)
 373     outputfile += Form(":PWG4_HighPtSpectraCent%dTrackType%dCuts%dkCentral",centClass,trackType,cuts);
 374   else if(triggerMask == AliVEvent::kSemiCentral)
 375     outputfile += Form(":PWG4_HighPtSpectraCent%dTrackType%dCuts%dkSemiCentral",centClass,trackType,cuts);
 376   else if(triggerMask == AliVEvent::kEMCEJE)
 377     outputfile += Form(":PWG4_HighPtSpectraCent%dTrackType%dCuts%dkEMCEJE",centClass,trackType,cuts);
 378   else if(triggerMask == AliVEvent::kMB|AliVEvent::kCentral|AliVEvent::kSemiCentral)
 379     outputfile += Form(":PWG4_HighPtSpectraCent%dTrackType%dCuts%dkMBkCentralkSemiCentral",centClass,trackType,cuts);
 380   else
 381     outputfile += Form(":PWG4_HighPtSpectraCent%dTrackType%dCuts%dPSF%d",centClass,trackType,cuts,triggerMask);
 382
 383
 384   AliAnalysisDataContainer *coutput0 = 0x0;
 385   AliAnalysisDataContainer *coutput1 = 0x0;
 386   AliAnalysisDataContainer *coutput2 = 0x0;
 387   AliAnalysisDataContainer *cout_cuts0 = 0x0;
 388
 389   if(triggerMask == AliVEvent::kMB) {
 390     coutput0 = mgr->CreateContainer(Form("chist0HighPtSpectraCent%dTrackType%dCuts%dkMB",centClass,trackType,cuts), TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 391     coutput1 = mgr->CreateContainer(Form("ccontainer0HighPtSpectraCent%dTrackType%dCuts%dkMB",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 392     coutput2 = mgr->CreateContainer(Form("ccontainer1HighPtSpectraCent%dTrackType%dCuts%dkMB",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 393     cout_cuts0 = mgr->CreateContainer(Form("qa_SpectraTrackCutsCent%dTrackType%dCuts%dkMB",centClass,trackType,cuts), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
 394   }
 395   else if(triggerMask == AliVEvent::kCentral) {
 396     coutput0 = mgr->CreateContainer(Form("chist0HighPtSpectraCent%dTrackType%dCuts%dkCentral",centClass,trackType,cuts), TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 397     coutput1 = mgr->CreateContainer(Form("ccontainer0HighPtSpectraCent%dTrackType%dCuts%dkCentral",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 398     coutput2 = mgr->CreateContainer(Form("ccontainer1HighPtSpectraCent%dTrackType%dCuts%dkCentral",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 399     cout_cuts0 = mgr->CreateContainer(Form("qa_SpectraTrackCutsCent%dTrackType%dCuts%dkCentral",centClass,trackType,cuts), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
 400   }
 401   else if(triggerMask == AliVEvent::kSemiCentral) {
 402     coutput0 = mgr->CreateContainer(Form("chist0HighPtSpectraCent%dTrackType%dCuts%dkSemiCentral",centClass,trackType,cuts), TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 403     coutput1 = mgr->CreateContainer(Form("ccontainer0HighPtSpectraCent%dTrackType%dCuts%dkSemiCentral",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 404     coutput2 = mgr->CreateContainer(Form("ccontainer1HighPtSpectraCent%dTrackType%dCuts%dkSemiCentral",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 405     cout_cuts0 = mgr->CreateContainer(Form("qa_SpectraTrackCutsCent%dTrackType%dCuts%dkSemiCentral",centClass,trackType,cuts), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
 406   }
 407   else if(triggerMask == AliVEvent::kEMCEJE) {
 408     coutput0 = mgr->CreateContainer(Form("chist0HighPtSpectraCent%dTrackType%dCuts%dkEMCEJE",centClass,trackType,cuts), TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 409     coutput1 = mgr->CreateContainer(Form("ccontainer0HighPtSpectraCent%dTrackType%dCuts%dkEMCEJE",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 410     coutput2 = mgr->CreateContainer(Form("ccontainer1HighPtSpectraCent%dTrackType%dCuts%dkEMCEJE",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 411     cout_cuts0 = mgr->CreateContainer(Form("qa_SpectraTrackCutsCent%dTrackType%dCuts%dkEMCEJE",centClass,trackType,cuts), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
 412   }
 413   else if(triggerMask == AliVEvent::kMB|AliVEvent::kCentral|AliVEvent::kSemiCentral) {
 414     coutput0 = mgr->CreateContainer(Form("chist0HighPtSpectraCent%dTrackType%dCuts%dkMBkCentralkSemiCentral",centClass,trackType,cuts), TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 415     coutput1 = mgr->CreateContainer(Form("ccontainer0HighPtSpectraCent%dTrackType%dCuts%dkMBkCentralkSemiCentral",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 416     coutput2 = mgr->CreateContainer(Form("ccontainer1HighPtSpectraCent%dTrackType%dCuts%dkMBkCentralkSemiCentral",centClass,trackType,cuts), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 417     cout_cuts0 = mgr->CreateContainer(Form("qa_SpectraTrackCutsCent%dTrackType%dCuts%dkMBkCentralkSemiCentral",centClass,trackType,cuts), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
 418   }
 419   else {
 420     coutput0 = mgr->CreateContainer(Form("chist0HighPtSpectraCent%dTrackType%dCuts%dPSF%d",centClass,trackType,cuts,triggerMask), TList::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 421     coutput1 = mgr->CreateContainer(Form("ccontainer0HighPtSpectraCent%dTrackType%dCuts%dPSF%d",centClass,trackType,cuts,triggerMask), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 422     coutput2 = mgr->CreateContainer(Form("ccontainer1HighPtSpectraCent%dTrackType%dCuts%dPSF%d",centClass,trackType,cuts,triggerMask), AliCFContainer::Class(),AliAnalysisManager::kOutputContainer,outputfile);
 423     cout_cuts0 = mgr->CreateContainer(Form("qa_SpectraTrackCutsCent%dTrackType%dCuts%dPSF%d",centClass,trackType,cuts,triggerMask), AliESDtrackCuts::Class(), AliAnalysisManager::kParamContainer,outputfile);
 424   }
 425
 426   mgr->AddTask(taskPWG4HighPtSpectra);
 427
 428   mgr->ConnectInput(taskPWG4HighPtSpectra,0,mgr->GetCommonInputContainer());
 429   mgr->ConnectOutput(taskPWG4HighPtSpectra,0,coutput0);
 430   mgr->ConnectOutput(taskPWG4HighPtSpectra,1,coutput1);
 431   mgr->ConnectOutput(taskPWG4HighPtSpectra,2,coutput2);
 432   mgr->ConnectOutput(taskPWG4HighPtSpectra,3,cout_cuts0);
 433
 434   // Return task pointer at the end
 435   return taskPWG4HighPtSpectra;
 436 }