]> git.uio.no Git - u/mrichter/AliRoot.git/commitdiff
updated bit 16 to 2010 cuts --> default cuts with wide (2-3 cm) DCA cut
authorjgrosseo <jgrosseo@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 20 Oct 2010 13:08:09 +0000 (13:08 +0000)
committerjgrosseo <jgrosseo@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 20 Oct 2010 13:08:09 +0000 (13:08 +0000)
adding bit 32 --> default cuts with strict DCA cut (for primaries)
adding bit 64 --> default cuts with strict DCA cut but requiring the first SDD layer instead of any of the SPD layers

ANALYSIS/macros/AddTaskESDFilter.C

index e11a8a5197e4ef6b6d0c47d440f83ac614d9b1c4..4d1e93bb1f1faf311b6c9fec899ed5bddd629c2e 100644 (file)
@@ -74,19 +74,20 @@ AliAnalysisTaskESDfilter *AddTaskESDFilter(Bool_t useKineFilter=kTRUE,
    AliESDpidCuts *electronID = new AliESDpidCuts("Electrons", "Electron PID cuts");
    electronID->SetTPCnSigmaCut(AliPID::kElectron, 3.);
 
-   // tighter cuts on primary particles for high pT tracks
-   // needed as input for jetfinder 
-   AliESDtrackCuts* esdTrackCutsH = new AliESDtrackCuts("Standard Track Cuts + ITSRefit", "High pT ESD Track Cuts");
-   esdTrackCutsH->SetMinNClustersTPC(50);
-   esdTrackCutsH->SetMaxChi2PerClusterTPC(3.5);
-   esdTrackCutsH->SetRequireTPCRefit(kTRUE);
+   // standard cuts with very loose DCA
+   AliESDtrackCuts* esdTrackCutsH = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kFALSE); 
    esdTrackCutsH->SetMaxDCAToVertexXY(2.4);
    esdTrackCutsH->SetMaxDCAToVertexZ(3.2);
    esdTrackCutsH->SetDCAToVertex2D(kTRUE);
-   esdTrackCutsH->SetRequireSigmaToVertex(kFALSE);
-   esdTrackCutsH->SetAcceptKinkDaughters(kFALSE);
-   esdTrackCutsH->SetRequireITSRefit(kTRUE); // additional cut 
 
+   // 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);
  
 
    // Compose the filter
@@ -103,6 +104,10 @@ AliAnalysisTaskESDfilter *AddTaskESDFilter(Bool_t useKineFilter=kTRUE,
    electronID->SetFilterMask(4);       // AND with Pixel Cuts
    // 16
    trackFilter->AddCuts(esdTrackCutsH);
+   // 32
+   trackFilter->AddCuts(esdTrackCutsH2);
+   // 64
+   trackFilter->AddCuts(esdTrackCutsH3);
  
    // Filter with cuts on V0s
    AliESDv0Cuts*   esdV0Cuts = new AliESDv0Cuts("Standard V0 Cuts pp", "ESD V0 Cuts");