ANALYSIS/macros/AddTaskESDFilter.C

   1 \r
   2 Bool_t AddTrackCutsLHC10bcde(AliAnalysisTaskESDfilter* esdFilter);\r
   3 Bool_t AddTrackCutsLHC10h(AliAnalysisTaskESDfilter* esdFilter);\r
   4 Bool_t AddTrackCutsLHC11h(AliAnalysisTaskESDfilter* esdFilter);\r
   5 Bool_t enableTPCOnlyAODTracksLocalFlag=kFALSE;\r
   6 \r
   7 \r
   8 AliAnalysisTaskESDfilter *AddTaskESDFilter(Bool_t useKineFilter=kTRUE, \r
   9                                            Bool_t writeMuonAOD=kFALSE,\r
  10                                            Bool_t writeDimuonAOD=kFALSE,\r
  11                                            Bool_t usePhysicsSelection=kFALSE,\r
  12                                            Bool_t useCentralityTask=kFALSE, /*obsolete*/\r
  13                                            Bool_t enableTPCOnlyAODTracks=kFALSE,\r
  14                                            Bool_t disableCascades=kFALSE,\r
  15                                            Bool_t disableKinks=kFALSE, \r
  16                                            Int_t runFlag = 1100,\r
  17                                            Int_t  muonMCMode = 2, \r
  18                                            Bool_t useV0Filter=kTRUE)\r
  19 {\r
  20   // Creates a filter task and adds it to the analysis manager.\r
  21    // Get the pointer to the existing analysis manager via the static access method.\r
  22    //==============================================================================\r
  23    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();\r
  24    if (!mgr) {\r
  25       ::Error("AddTaskESDFilter", "No analysis manager to connect to.");\r
  26      return NULL;\r
  27    }   \r
  28    \r
  29    // This task requires an ESD input handler and an AOD output handler.\r
  30    // Check this using the analysis manager.\r
  31    //===============================================================================\r
  32    TString type = mgr->GetInputEventHandler()->GetDataType();\r
  33    if (!type.Contains("ESD")) {\r
  34       ::Error("AddTaskESDFilter", "ESD filtering task needs the manager to have an ESD input handler.");\r
  35       return NULL;\r
  36    }   \r
  37    // Check if AOD output handler exist.\r
  38    AliAODHandler *aod_h = (AliAODHandler*)mgr->GetOutputEventHandler();\r
  39    if (!aod_h) {\r
  40       ::Error("AddTaskESDFilter", "ESD filtering task needs the manager to have an AOD output handler.");\r
  41       return NULL;\r
  42    }\r
  43    // Check if MC handler is connected in case kine filter requested\r
  44    AliMCEventHandler *mcH = (AliMCEventHandler*)mgr->GetMCtruthEventHandler();\r
  45    if (!mcH && useKineFilter) {\r
  46       ::Error("AddTaskESDFilter", "No MC handler connected while kine filtering requested");\r
  47       return NULL;\r
  48    }   \r
  49    \r
  50    // Create the task, add it to the manager and configure it.\r
  51    //===========================================================================   \r
  52    // Barrel tracks filter\r
  53    AliAnalysisTaskESDfilter *esdfilter = new AliAnalysisTaskESDfilter("ESD Filter");\r
  54    if (disableCascades) esdfilter->DisableCascades();\r
  55    if  (disableKinks) esdfilter->DisableKinks();\r
  56   \r
  57    mgr->AddTask(esdfilter);\r
  58   \r
  59    // Muons\r
  60    Bool_t onlyMuon=kTRUE;\r
  61    Bool_t keepAllEvents=kTRUE;\r
  62    Int_t mcMode= useKineFilter ? muonMCMode : 0; // use 1 instead of 2 to get all MC information instead of just ancestors of mu tracks\r
  63    AliAnalysisTaskESDMuonFilter *esdmuonfilter = new AliAnalysisTaskESDMuonFilter("ESD Muon Filter",onlyMuon,keepAllEvents,mcMode);\r
  64    mgr->AddTask(esdmuonfilter);\r
  65    if(usePhysicsSelection){\r
  66      esdfilter->SelectCollisionCandidates(AliVEvent::kAny);\r
  67      esdmuonfilter->SelectCollisionCandidates(AliVEvent::kAny);\r
  68    }  \r
  69 \r
  70    // Filtering of MC particles (decays conversions etc)\r
  71    // this task has to go AFTER all other filter tasks\r
  72    // since it fills the AODMC array with all\r
  73    // selected MC Particles, only this way we have the \r
  74    // AODMCparticle information available for following tasks\r
  75    AliAnalysisTaskMCParticleFilter *kinefilter = 0;\r
  76    if (useKineFilter) {\r
  77       kinefilter = new AliAnalysisTaskMCParticleFilter("Particle Kine Filter");\r
  78       mgr->AddTask(kinefilter);\r
  79    }   \r
  80 \r
  81    enableTPCOnlyAODTracksLocalFlag = enableTPCOnlyAODTracks;\r
  82    if((runFlag/100)==10){\r
  83      if((runFlag%100)==0) AddTrackCutsLHC10bcde(esdfilter);\r
  84      else AddTrackCutsLHC10h(esdfilter);\r
  85    }\r
  86    else {\r
  87      // default 11h\r
  88      AddTrackCutsLHC11h(esdfilter);\r
  89    }\r
  90 \r
  91    // Filter with cuts on V0s\r
  92    if (useV0Filter) {\r
  93      AliESDv0Cuts*   esdV0Cuts = new AliESDv0Cuts("Standard V0 Cuts pp", "ESD V0 Cuts");\r
  94      esdV0Cuts->SetMinRadius(0.2);\r
  95      esdV0Cuts->SetMaxRadius(200);\r
  96      esdV0Cuts->SetMinDcaPosToVertex(0.05);\r
  97      esdV0Cuts->SetMinDcaNegToVertex(0.05);\r
  98      esdV0Cuts->SetMaxDcaV0Daughters(1.5);\r
  99      esdV0Cuts->SetMinCosinePointingAngle(0.99);\r
 100      AliAnalysisFilter* v0Filter = new AliAnalysisFilter("v0Filter");\r
 101      v0Filter->AddCuts(esdV0Cuts);\r
 102 \r
 103      esdfilter->SetV0Filter(v0Filter);\r
 104    }  \r
 105 \r
 106    // Enable writing of Muon AODs\r
 107    esdmuonfilter->SetWriteMuonAOD(writeMuonAOD);\r
 108    \r
 109    // Enable writing of Dimuon AODs\r
 110    esdmuonfilter->SetWriteDimuonAOD(writeDimuonAOD);\r
 111  \r
 112    // Create ONLY the output containers for the data produced by the task.\r
 113    // Get and connect other common input/output containers via the manager as below\r
 114    //==============================================================================\r
 115    mgr->ConnectInput  (esdfilter,  0, mgr->GetCommonInputContainer());\r
 116    mgr->ConnectOutput (esdfilter,  0, mgr->GetCommonOutputContainer());\r
 117    mgr->ConnectInput  (esdmuonfilter, 0, mgr->GetCommonInputContainer());\r
 118    if (useKineFilter) {\r
 119       mgr->ConnectInput  (kinefilter,  0, mgr->GetCommonInputContainer());\r
 120       mgr->ConnectOutput (kinefilter,  0, mgr->GetCommonOutputContainer());\r
 121       AliAnalysisDataContainer *coutputEx = mgr->CreateContainer("cFilterList", TList::Class(),\r
 122                                                                    AliAnalysisManager::kOutputContainer,"pyxsec_hists.root");\r
 123       mgr->ConnectOutput (kinefilter,  1,coutputEx);\r
 124    }   \r
 125    return esdfilter;\r
 126 }\r
 127  \r
 128 \r
 129 \r
 130 \r
 131 Bool_t AddTrackCutsLHC10h(AliAnalysisTaskESDfilter* esdfilter){\r
 132 \r
 133   Printf("%s%d: Creating Track Cuts for LHC10h",(char*)__FILE__,__LINE__);\r
 134 \r
 135   // Cuts on primary tracks\r
 136   AliESDtrackCuts* esdTrackCutsL = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();\r
 137   \r
 138   // ITS stand-alone tracks\r
 139   AliESDtrackCuts* esdTrackCutsITSsa = new AliESDtrackCuts("ITS stand-alone Track Cuts", "ESD Track Cuts");\r
 140   esdTrackCutsITSsa->SetRequireITSStandAlone(kTRUE);\r
 141   \r
 142   // Pixel OR necessary for the electrons\r
 143   AliESDtrackCuts *itsStrong = new AliESDtrackCuts("ITSorSPD", "pixel requirement for ITS");\r
 144   itsStrong->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);\r
 145   \r
 146   \r
 147   // PID for the electrons\r
 148   AliESDpidCuts *electronID = new AliESDpidCuts("Electrons", "Electron PID cuts");\r
 149   electronID->SetTPCnSigmaCut(AliPID::kElectron, 3.5);\r
 150   \r
 151   // tighter cuts on primary particles for high pT tracks\r
 152   // take the standard cuts, which include already \r
 153   // ITSrefit and use only primaries...\r
 154   \r
 155   // ITS cuts for new jet analysis \r
 156   //  gROOT->LoadMacro("$ALICE_ROOT/PWGJE/macros/CreateTrackCutsPWGJE.C");\r
 157   //  AliESDtrackCuts* esdTrackCutsHG0 = CreateTrackCutsPWGJE(10001006);\r
 158 \r
 159   AliESDtrackCuts *jetCuts1006 = new AliESDtrackCuts("AliESDtrackCuts"); \r
 160 \r
 161   TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");\r
 162   jetCuts1006->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep,20.);\r
 163   jetCuts1006->SetMinNClustersTPC(70);\r
 164   jetCuts1006->SetMaxChi2PerClusterTPC(4);\r
 165   jetCuts1006->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1\r
 166   jetCuts1006->SetAcceptKinkDaughters(kFALSE);\r
 167   jetCuts1006->SetRequireTPCRefit(kTRUE);\r
 168   jetCuts1006->SetMaxFractionSharedTPCClusters(0.4);\r
 169   // ITS\r
 170   jetCuts1006->SetRequireITSRefit(kTRUE);\r
 171   //accept secondaries\r
 172   jetCuts1006->SetMaxDCAToVertexXY(2.4);\r
 173   jetCuts1006->SetMaxDCAToVertexZ(3.2);\r
 174   jetCuts1006->SetDCAToVertex2D(kTRUE);\r
 175   //reject fakes\r
 176   jetCuts1006->SetMaxChi2PerClusterITS(36);\r
 177   jetCuts1006->SetMaxChi2TPCConstrainedGlobal(36);\r
 178 \r
 179   jetCuts1006->SetRequireSigmaToVertex(kFALSE);\r
 180 \r
 181   jetCuts1006->SetEtaRange(-0.9,0.9);\r
 182   jetCuts1006->SetPtRange(0.15, 1E+15.);\r
 183 \r
 184   AliESDtrackCuts* esdTrackCutsHG0 = jetCuts1006->Clone("JetCuts10001006");\r
 185   esdTrackCutsHG0->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);\r
 186 \r
 187 \r
 188   // throw out tracks with too low number of clusters in\r
 189   // the first pass (be consistent with TPC only tracks)\r
 190   // N.B. the number off crossed rows still acts on the tracks after\r
 191   // all iterations if we require tpc standalone, number of clusters\r
 192   // and chi2 TPC cuts act on track after the first iteration\r
 193   //   esdTrackCutsH0->SetRequireTPCStandAlone(kTRUE);\r
 194   //   esdTrackCutsH0->SetMinNClustersTPC(80); // <--- first pass\r
 195   \r
 196   \r
 197   // the complement to the one with SPD requirement\r
 198   //  AliESDtrackCuts* esdTrackCutsHG1 = CreateTrackCutsPWGJE(10011006);\r
 199   AliESDtrackCuts* esdTrackCutsHG1 = jetCuts1006->Clone("JetCuts10011006");\r
 200   esdTrackCutsHG1->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);\r
 201 \r
 202   // the tracks that must not be taken pass this cut and\r
 203   // non HGC1 and HG\r
 204   //  AliESDtrackCuts* esdTrackCutsHG2 = CreateTrackCutsPWGJE(10021006);\r
 205   AliESDtrackCuts* esdTrackCutsHG2 = jetCuts1006->Clone("JetCuts10021006");\r
 206   esdTrackCutsHG2->SetMaxChi2PerClusterITS(1E10);\r
 207 \r
 208 \r
 209   // standard cuts also used in R_AA analysis\r
 210   //   "Global track RAA analysis QM2011 + Chi2ITS<36";\r
 211   //  AliESDtrackCuts* esdTrackCutsH2 = CreateTrackCutsPWGJE(1000);\r
 212   AliESDtrackCuts* esdTrackCutsH2 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kTRUE,1);\r
 213   esdTrackCutsH2->SetMinNCrossedRowsTPC(120);\r
 214   esdTrackCutsH2->SetMinRatioCrossedRowsOverFindableClustersTPC(0.8);\r
 215   esdTrackCutsH2->SetMaxChi2PerClusterITS(36);\r
 216   esdTrackCutsH2->SetMaxFractionSharedTPCClusters(0.4);\r
 217   esdTrackCutsH2->SetMaxChi2TPCConstrainedGlobal(36);\r
 218 \r
 219   esdTrackCutsH2->SetEtaRange(-0.9,0.9);\r
 220   esdTrackCutsH2->SetPtRange(0.15, 1e10);\r
 221 \r
 222 \r
 223   //  AliESDtrackCuts* esdTrackCutsGCOnly = CreateTrackCutsPWGJE(10041006);\r
 224   AliESDtrackCuts* esdTrackCutsGCOnly = jetCuts1006->Clone("JetCuts10041006");\r
 225   esdTrackCutsGCOnly->SetRequireITSRefit(kFALSE);\r
 226 \r
 227   // TPC only tracks\r
 228   AliESDtrackCuts* esdTrackCutsTPCCOnly = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();\r
 229   esdTrackCutsTPCCOnly->SetMinNClustersTPC(70);\r
 230   \r
 231   // Compose the filter\r
 232   AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");\r
 233   // 1, 1<<0\r
 234   trackFilter->AddCuts(esdTrackCutsL);\r
 235   // 2 1<<1\r
 236   trackFilter->AddCuts(esdTrackCutsITSsa);\r
 237   // 4 1<<2\r
 238   trackFilter->AddCuts(itsStrong);\r
 239   itsStrong->SetFilterMask(1);        // AND with Standard track cuts \r
 240   // 8 1<<3\r
 241   trackFilter->AddCuts(electronID);\r
 242   electronID->SetFilterMask(4);       // AND with Pixel Cuts\r
 243   // 16 1<<4\r
 244   trackFilter->AddCuts(esdTrackCutsHG0);\r
 245   // 32 1<<5\r
 246   trackFilter->AddCuts(esdTrackCutsHG1);\r
 247   // 64 1<<6\r
 248   trackFilter->AddCuts(esdTrackCutsHG2);\r
 249   // 128 1<<7\r
 250   trackFilter->AddCuts(esdTrackCutsTPCCOnly); // add QM TPC only track cuts\r
 251   if(enableTPCOnlyAODTracksLocalFlag)esdfilter->SetTPCOnlyFilterMask(128);\r
 252   // 256 1<<8\r
 253   trackFilter->AddCuts(esdTrackCutsGCOnly);\r
 254   // 512 1<<9                         \r
 255   AliESDtrackCuts* esdTrackCutsHG1_tmp = new AliESDtrackCuts(*esdTrackCutsHG1); // avoid double delete\r
 256   trackFilter->AddCuts(esdTrackCutsHG1_tmp); // add once more for tpc only tracks\r
 257   // 1024 1<<10                        \r
 258   trackFilter->AddCuts(esdTrackCutsH2); // add r_aa cuts\r
 259 \r
 260   \r
 261   \r
 262   esdfilter->SetGlobalConstrainedFilterMask(1<<8|1<<9); // these tracks are written out as global constrained tracks\r
 263   esdfilter->SetHybridFilterMaskGlobalConstrainedGlobal((1<<4)); // these normal global tracks will be marked as hybrid\r
 264   esdfilter->SetWriteHybridGlobalConstrainedOnly(kTRUE); // write only the complement\r
 265   //     esdfilter->SetTPCConstrainedFilterMask(1<<11); // these tracks are written out as tpc constrained tracks\r
 266 \r
 267   esdfilter->SetTrackFilter(trackFilter);\r
 268   return kTRUE;\r
 269   \r
 270 }\r
 271 \r
 272 \r
 273 Bool_t AddTrackCutsLHC11h(AliAnalysisTaskESDfilter* esdfilter){\r
 274 \r
 275 \r
 276   Printf("%s%d: Creating Track Cuts LHC11h",(char*)__FILE__,__LINE__);\r
 277 \r
 278   // Cuts on primary tracks\r
 279    AliESDtrackCuts* esdTrackCutsL = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();\r
 280 \r
 281    // ITS stand-alone tracks\r
 282    AliESDtrackCuts* esdTrackCutsITSsa = new AliESDtrackCuts("ITS stand-alone Track Cuts", "ESD Track Cuts");\r
 283    esdTrackCutsITSsa->SetRequireITSStandAlone(kTRUE);\r
 284 \r
 285    // Pixel OR necessary for the electrons\r
 286    AliESDtrackCuts *itsStrong = new AliESDtrackCuts("ITSorSPD", "pixel requirement for ITS");\r
 287    itsStrong->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);\r
 288 \r
 289 \r
 290    // PID for the electrons\r
 291    AliESDpidCuts *electronID = new AliESDpidCuts("Electrons", "Electron PID cuts");\r
 292    electronID->SetTPCnSigmaCut(AliPID::kElectron, 3.5);\r
 293 \r
 294    // standard cuts with very loose DCA\r
 295    AliESDtrackCuts* esdTrackCutsH = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kFALSE); \r
 296    esdTrackCutsH->SetMaxDCAToVertexXY(2.4);\r
 297    esdTrackCutsH->SetMaxDCAToVertexZ(3.2);\r
 298    esdTrackCutsH->SetDCAToVertex2D(kTRUE);\r
 299 \r
 300    // standard cuts with tight DCA cut\r
 301    AliESDtrackCuts* esdTrackCutsH2 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011();\r
 302 \r
 303    // standard cuts with tight DCA but with requiring the first SDD cluster instead of an SPD cluster\r
 304    // tracks selected by this cut are exclusive to those selected by the previous cut\r
 305    AliESDtrackCuts* esdTrackCutsH3 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(); \r
 306    esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);\r
 307    esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSDD, AliESDtrackCuts::kFirst);\r
 308  \r
 309    // TPC only tracks: Optionally enable the writing of TPConly information\r
 310    // constrained to SPD vertex in the filter below\r
 311    AliESDtrackCuts* esdTrackCutsTPCOnly = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();\r
 312    // The following line is needed for 2010 PbPb reprocessing and pp, but not for 2011 PbPb\r
 313    //esdTrackCutsTPCOnly->SetMinNClustersTPC(70);\r
 314 \r
 315    // Extra cuts for hybrids\r
 316    // first the global tracks we want to take\r
 317    AliESDtrackCuts* esdTrackCutsHTG = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kFALSE); \r
 318    esdTrackCutsHTG->SetName("Global Hybrid tracks, loose DCA");\r
 319    esdTrackCutsHTG->SetMaxDCAToVertexXY(2.4);\r
 320    esdTrackCutsHTG->SetMaxDCAToVertexZ(3.2);\r
 321    esdTrackCutsHTG->SetDCAToVertex2D(kTRUE);\r
 322    esdTrackCutsHTG->SetMaxChi2TPCConstrainedGlobal(36);\r
 323    esdTrackCutsHTG->SetMaxFractionSharedTPCClusters(0.4);\r
 324    \r
 325 \r
 326    // Than the complementary tracks which will be stored as global\r
 327    // constraint, complement is done in the ESDFilter task\r
 328    AliESDtrackCuts* esdTrackCutsHTGC = new AliESDtrackCuts(*esdTrackCutsHTG);\r
 329    esdTrackCutsHTGC->SetName("Global Constraint Hybrid tracks, loose DCA no it requirement");\r
 330    esdTrackCutsHTGC->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kOff);\r
 331    esdTrackCutsHTGC->SetRequireITSRefit(kTRUE);\r
 332 \r
 333    // standard cuts with tight DCA cut, using cluster cut instead of crossed rows (a la 2010 default)\r
 334    AliESDtrackCuts* esdTrackCutsH2Cluster = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011(kTRUE, 0);\r
 335 \r
 336    // Compose the filter\r
 337    AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");\r
 338    // 1, 1<<0\r
 339    trackFilter->AddCuts(esdTrackCutsL);\r
 340    // 2, 1<<1\r
 341    trackFilter->AddCuts(esdTrackCutsITSsa);\r
 342    // 4, 1<<2\r
 343    trackFilter->AddCuts(itsStrong);\r
 344    itsStrong->SetFilterMask(1);        // AND with Standard track cuts \r
 345    // 8, 1<<3\r
 346    trackFilter->AddCuts(electronID);\r
 347    electronID->SetFilterMask(4);       // AND with Pixel Cuts\r
 348    // 16, 1<<4\r
 349    trackFilter->AddCuts(esdTrackCutsH);\r
 350    // 32, 1<<5\r
 351    trackFilter->AddCuts(esdTrackCutsH2);\r
 352    // 64, 1<<6\r
 353    trackFilter->AddCuts(esdTrackCutsH3);\r
 354    // 128 , 1 << 7\r
 355    trackFilter->AddCuts(esdTrackCutsTPCOnly);\r
 356    if(enableTPCOnlyAODTracksLocalFlag)esdfilter->SetTPCOnlyFilterMask(128);\r
 357    // 256, 1 << 8 Global Hybrids\r
 358    trackFilter->AddCuts(esdTrackCutsHTG);\r
 359    esdfilter->SetHybridFilterMaskGlobalConstrainedGlobal((1<<8)); // these normal global tracks will be marked as hybrid    \r
 360    // 512, 1<< 9 GlobalConstraint Hybrids\r
 361    trackFilter->AddCuts(esdTrackCutsHTGC);\r
 362    esdfilter->SetGlobalConstrainedFilterMask(1<<9); // these tracks are written out as global constrained tracks \r
 363    esdfilter->SetWriteHybridGlobalConstrainedOnly(kTRUE); // write only the complement\r
 364    // 1024, 1<< 10 // tight DCA cuts\r
 365    trackFilter->AddCuts(esdTrackCutsH2Cluster);\r
 366    esdfilter->SetTrackFilter(trackFilter);\r
 367 \r
 368    return kTRUE;\r
 369 \r
 370 }\r
 371 \r
 372 Bool_t AddTrackCutsLHC10bcde(AliAnalysisTaskESDfilter* esdfilter){\r
 373 \r
 374 \r
 375   Printf("%s%d: Creating Track Cuts LHC10bcde",(char*)__FILE__,__LINE__);\r
 376 \r
 377   // Cuts on primary tracks\r
 378    AliESDtrackCuts* esdTrackCutsL = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();\r
 379 \r
 380    // ITS stand-alone tracks\r
 381    AliESDtrackCuts* esdTrackCutsITSsa = new AliESDtrackCuts("ITS stand-alone Track Cuts", "ESD Track Cuts");\r
 382    esdTrackCutsITSsa->SetRequireITSStandAlone(kTRUE);\r
 383 \r
 384    // Pixel OR necessary for the electrons\r
 385    AliESDtrackCuts *itsStrong = new AliESDtrackCuts("ITSorSPD", "pixel requirement for ITS");\r
 386    itsStrong->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);\r
 387 \r
 388 \r
 389    // PID for the electrons\r
 390    AliESDpidCuts *electronID = new AliESDpidCuts("Electrons", "Electron PID cuts");\r
 391    electronID->SetTPCnSigmaCut(AliPID::kElectron, 3.5);\r
 392 \r
 393    // standard cuts with very loose DCA\r
 394    AliESDtrackCuts* esdTrackCutsH = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kFALSE); \r
 395    esdTrackCutsH->SetMaxDCAToVertexXY(2.4);\r
 396    esdTrackCutsH->SetMaxDCAToVertexZ(3.2);\r
 397    esdTrackCutsH->SetDCAToVertex2D(kTRUE);\r
 398 \r
 399    // standard cuts with tight DCA cut\r
 400    AliESDtrackCuts* esdTrackCutsH2 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010();\r
 401 \r
 402    // standard cuts with tight DCA but with requiring the first SDD cluster instead of an SPD cluster\r
 403    // tracks selected by this cut are exclusive to those selected by the previous cut\r
 404    AliESDtrackCuts* esdTrackCutsH3 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(); \r
 405    esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);\r
 406    esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSDD, AliESDtrackCuts::kFirst);\r
 407  \r
 408    // TPC only tracks: Optionally enable the writing of TPConly information\r
 409    // constrained to SPD vertex in the filter below\r
 410    AliESDtrackCuts* esdTrackCutsTPCOnly = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();\r
 411    // The following line is needed for 2010 PbPb reprocessing and pp, but not for 2011 PbPb\r
 412    esdTrackCutsTPCOnly->SetMinNClustersTPC(70);\r
 413 \r
 414    // Extra cuts for hybrids\r
 415    // first the global tracks we want to take\r
 416    // take the HTGs from 10h\r
 417 \r
 418    //  gROOT->LoadMacro("$ALICE_ROOT/PWGJE/macros/CreateTrackCutsPWGJE.C");\r
 419    //  AliESDtrackCuts* esdTrackCutsHG0 = CreateTrackCutsPWGJE(10001006);\r
 420 \r
 421    AliESDtrackCuts *jetCuts1006 = new AliESDtrackCuts("AliESDtrackCuts"); \r
 422 \r
 423    TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");\r
 424    jetCuts1006->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep,20.);\r
 425    jetCuts1006->SetMinNClustersTPC(70);\r
 426    jetCuts1006->SetMaxChi2PerClusterTPC(4);\r
 427    jetCuts1006->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1\r
 428    jetCuts1006->SetAcceptKinkDaughters(kFALSE);\r
 429    jetCuts1006->SetRequireTPCRefit(kTRUE);\r
 430    jetCuts1006->SetMaxFractionSharedTPCClusters(0.4);\r
 431    // ITS\r
 432    jetCuts1006->SetRequireITSRefit(kTRUE);\r
 433    //accept secondaries\r
 434    jetCuts1006->SetMaxDCAToVertexXY(2.4);\r
 435    jetCuts1006->SetMaxDCAToVertexZ(3.2);\r
 436    jetCuts1006->SetDCAToVertex2D(kTRUE);\r
 437    //reject fakes\r
 438    jetCuts1006->SetMaxChi2PerClusterITS(36);\r
 439    jetCuts1006->SetMaxChi2TPCConstrainedGlobal(36);\r
 440 \r
 441    jetCuts1006->SetRequireSigmaToVertex(kFALSE);\r
 442    \r
 443    jetCuts1006->SetEtaRange(-0.9,0.9);\r
 444    jetCuts1006->SetPtRange(0.15, 1E+15.);\r
 445    \r
 446    AliESDtrackCuts* esdTrackCutsHTG = jetCuts1006->Clone("JetCuts10001006");\r
 447    esdTrackCutsHTG->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);\r
 448 \r
 449    // Than the complementary tracks which will be stored as global\r
 450    // constraint, complement is done in the ESDFilter task\r
 451    // HGC from 10h\r
 452 \r
 453    AliESDtrackCuts* esdTrackCutsHTGC = jetCuts1006->Clone("JetCuts10041006");\r
 454    esdTrackCutsHTGC->SetRequireITSRefit(kFALSE);\r
 455 \r
 456    // standard cuts with tight DCA cut, using cluster cut instead of crossed rows (a la 2010 default)\r
 457    AliESDtrackCuts* esdTrackCutsH2Cluster = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kTRUE, 0);\r
 458 \r
 459    // Compose the filter\r
 460    AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");\r
 461    // 1, 1<<0\r
 462    trackFilter->AddCuts(esdTrackCutsL);\r
 463    // 2, 1<<1\r
 464    trackFilter->AddCuts(esdTrackCutsITSsa);\r
 465    // 4, 1<<2\r
 466    trackFilter->AddCuts(itsStrong);\r
 467    itsStrong->SetFilterMask(1);        // AND with Standard track cuts \r
 468    // 8, 1<<3\r
 469    trackFilter->AddCuts(electronID);\r
 470    electronID->SetFilterMask(4);       // AND with Pixel Cuts\r
 471    // 16, 1<<4\r
 472    trackFilter->AddCuts(esdTrackCutsH);\r
 473    // 32, 1<<5\r
 474    trackFilter->AddCuts(esdTrackCutsH2);\r
 475    // 64, 1<<6\r
 476    trackFilter->AddCuts(esdTrackCutsH3);\r
 477    // 128 , 1 << 7\r
 478    trackFilter->AddCuts(esdTrackCutsTPCOnly);\r
 479    if(enableTPCOnlyAODTracksLocalFlag)esdfilter->SetTPCOnlyFilterMask(128);\r
 480    // 256, 1 << 8 Global Hybrids\r
 481    trackFilter->AddCuts(esdTrackCutsHTG);\r
 482    esdfilter->SetHybridFilterMaskGlobalConstrainedGlobal((1<<8)); // these normal global tracks will be marked as hybrid    \r
 483    // 512, 1<< 9 GlobalConstraint Hybrids\r
 484    trackFilter->AddCuts(esdTrackCutsHTGC);\r
 485    esdfilter->SetGlobalConstrainedFilterMask(1<<9); // these tracks are written out as global constrained tracks \r
 486    esdfilter->SetWriteHybridGlobalConstrainedOnly(kTRUE); // write only the complement\r
 487    // 1024, 1<< 10 // tight DCA cuts\r
 488    trackFilter->AddCuts(esdTrackCutsH2Cluster);\r
 489    esdfilter->SetTrackFilter(trackFilter);\r
 490 \r
 491    return kTRUE;\r
 492 \r
 493 }\r