// tighter cuts on primary particles for high pT tracks
// take the standard cuts, which include already
// ITSrefit and use only primaries...
- AliESDtrackCuts* esdTrackCutsH = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kTRUE);
- esdTrackCutsH->SetName("StandardFromAliESDTrackCuts2010");
+
+ AliESDtrackCuts* esdTrackCutsH = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kFALSE);
+ esdTrackCutsH->SetMaxDCAToVertexXY(2.4);
+ esdTrackCutsH->SetMaxDCAToVertexZ(3.2);
+ esdTrackCutsH->SetDCAToVertex2D(kTRUE);
+
+ // standard cuts with tight DCA cut
+ AliESDtrackCuts* esdTrackCutsH2 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010();
+
+ // standard cuts with tight DCA but with requiring the first SDD cluster instead of an SPD cluster
+ // tracks selected by this cut are exclusive to those selected by the previous cut
+ AliESDtrackCuts* esdTrackCutsH3 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010();
+ esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);
+ esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSDD, AliESDtrackCuts::kFirst);
+
// these are needed for older MC
esdTrackCutsH1->SetName("StandardFromAliESDTrackCutsMaxRelPt");
esdTrackCutsH1->SetMaxRel1PtUncertainty(0.4); // new
-
- AliESDtrackCuts* esdTrackCutsH2 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010();
- esdTrackCutsH2->SetName("StandardFromAliESDTrackCutsExtraDCA");
- esdTrackCutsH2->SetMaxDCAToVertexZ(2);
-
- AliESDtrackCuts* esdTrackCutsH3 = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010();
- esdTrackCutsH2->SetName("StandardFromAliESDTrackCutsExtraDCANoITS");
- esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);
- esdTrackCutsH3->SetClusterRequirementITS(AliESDtrackCuts::kSDD, AliESDtrackCuts::kFirst);
- esdTrackCutsH3->SetMaxDCAToVertexZ(2);
-
// Compose the filter
AliAnalysisFilter* trackFilter = new AliAnalysisFilter("trackFilter");
// 1
// 16
trackFilter->AddCuts(esdTrackCutsH);
// 32
- trackFilter->AddCuts(esdTrackCutsH1);
- // 64
trackFilter->AddCuts(esdTrackCutsH2);
+ // 64
+ trackFilter->AddCuts(esdTrackCutsH3);
// 128
- trackFilter->AddCuts(esdTrackCutsH2);
+ trackFilter->AddCuts(esdTrackCutsH1);
// Filter with cuts on V0s
AliESDv0Cuts* esdV0Cuts = new AliESDv0Cuts("Standard V0 Cuts pp", "ESD V0 Cuts");
-AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2(const char* bRec = "jets",const char* bGen = "jetsAODMC_UA104",const char* bBkg="jeteventbackground_clustersAOD_KT06",UInt_t filterMask = 16, Int_t iPhysicsSelection = 1,UInt_t iEventSelectionMask = 0,Bool_t kBackground=kTRUE);\r
+AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2(const char* bRec = "jets",const char* bGen = "jetsAODMC_UA104",const char* bBkg="jeteventbackground_clustersAOD_KT06",UInt_t filterMask = 32, Int_t iPhysicsSelection = 1,UInt_t iEventSelectionMask = 0,Bool_t kBackground=kTRUE);\r
\r
\r
-AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2Delta(UInt_t filterMask = 16,Bool_t kUseAODMC = kFALSE,Int_t iPhysicsSelection = 1,UInt_t iFlag = 0xfffffff, UInt_t iEventSelectionMask = 0,Bool_t kBackground = kFALSE){\r
+AliAnalysisTaskJetSpectrum2 *AddTaskJetSpectrum2Delta(UInt_t filterMask = 32,Bool_t kUseAODMC = kFALSE,Int_t iPhysicsSelection = 1,UInt_t iFlag = 0xfffffff, UInt_t iEventSelectionMask = 0,Bool_t kBackground = kFALSE){\r
AliAnalysisTaskJetSpectrum2 *js = 0;\r
\r
TString cBack = "";\r
jetana = new AliAnalysisTaskJets(Form("JetAnalysis%s_%s%s",jr,jf,cRadius));\r
\r
TString cAdd = "";\r
- if(filterMask==16){\r
+ if(filterMask==32){\r
// this is the standard mask do not add anything\r
// should be changed after current train over all data is finished \r
// now needed for merging\r
AliJetAODReaderHeader *jrh = new AliJetAODReaderHeader();\r
jrh->SetComment("AOD Reader");\r
jrh->SetPtCut(0.15); // set low p_T cut of to 150 MeV\r
- jrh->SetTestFilterMask(16); // Change this one for a different set of cuts\r
+ jrh->SetTestFilterMask(32); // Change this one for a different set of cuts\r
if(filterMask>0)jrh->SetTestFilterMask(filterMask); \r
// Define reader and set its header\r
er = new AliJetAODReader();\r