AliRsnCutDaughterSigmaStar2010PP *cutPi = new AliRsnCutDaughterSigmaStar2010PP("cutPionForSigmaStar", AliPID::kPion);
cutPi->SetPIDCut(piPIDCut);
AliRsnCutTrackQuality *cutQuality = (AliRsnCutTrackQuality*) cutPi->CutQuality();
+ cutQuality->SetAODTestFilterBit(aodFilterBit);
cutQuality->SetDCARmax(0.05);
// cut set
// add to task
Int_t iCutPi = task->AddTrackCuts(cutSetPi);
- // added by EF
- AliRsnCutTrackQuality *cq = (AliRsnCutTrackQuality*) cutPi->CutQuality();
- cq->Print();
- cq->SetAODTestFilterBit(-1);
-
// quality cuts
AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("qualityDaughterLambda");
esdTrackCuts->SetAcceptKinkDaughters(0); // 0 = kFalse
- //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
+ esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
esdTrackCuts->SetMaxChi2PerClusterTPC(4);
esdTrackCuts->SetMinNClustersTPC(70);
esdTrackCuts->SetRequireTPCRefit();
cutLambda->SetMaxDCAVertex(0.3);
cutLambda->SetMinCosPointingAngle(0.99);
cutLambda->SetMaxDaughtersDCA(0.5);
- cutLambda->SetMaxRapidity(0.5);
+ cutLambda->SetMaxRapidity(0.8);
cutLambda->SetAODTestFilterBit(aodFilterBit);
cutLambda->SetPIDCut1(pPIDCut);
cutLambda->SetPIDCut2(piPIDCut);
cutAntiLambda->SetMaxDCAVertex(0.3);
cutAntiLambda->SetMinCosPointingAngle(0.99);
cutAntiLambda->SetMaxDaughtersDCA(0.5);
- cutAntiLambda->SetMaxRapidity(0.5);
+ cutAntiLambda->SetMaxRapidity(0.8);
cutAntiLambda->SetAODTestFilterBit(aodFilterBit);
cutAntiLambda->SetPIDCut1(pPIDCut);
cutAntiLambda->SetPIDCut2(piPIDCut);
AliRsnCutDaughterSigmaStar2010PP *cutPiLoose = new AliRsnCutDaughterSigmaStar2010PP("cutPionForSigmaStarLoose", AliPID::kPion);
cutPiLoose->SetPIDCut(piPIDCut);
AliRsnCutTrackQuality *cutQuality = (AliRsnCutTrackQuality*) cutPi->CutQuality();
+ cutQuality->SetAODTestFilterBit(aodFilterBit);
cutQuality->SetDCARmax(0.072);
// cut set
cutLambdaLoose->SetMaxDCAVertex(0.4);
cutLambdaLoose->SetMinCosPointingAngle(0.82);
cutLambdaLoose->SetMaxDaughtersDCA(0.67);
- cutLambdaLoose->SetMaxRapidity(0.5);
+ cutLambdaLoose->SetMaxRapidity(0.8);
cutLambdaLoose->SetAODTestFilterBit(aodFilterBit);
cutLambdaLoose->SetPIDCut1(pPIDCut);
cutLambdaLoose->SetPIDCut2(piPIDCut);
cutAntiLambdaLoose->SetMaxDCAVertex(0.4);
cutAntiLambdaLoose->SetMinCosPointingAngle(0.982);
cutAntiLambdaLoose->SetMaxDaughtersDCA(0.67);
- cutAntiLambdaLoose->SetMaxRapidity(0.5);
+ cutAntiLambdaLoose->SetMaxRapidity(0.8);
cutAntiLambdaLoose->SetAODTestFilterBit(aodFilterBit);
cutAntiLambdaLoose->SetPIDCut1(pPIDCut);
cutAntiLambdaLoose->SetPIDCut2(piPIDCut);
AliRsnCutDaughterSigmaStar2010PP *cutPiTight = new AliRsnCutDaughterSigmaStar2010PP("cutPionForSigmaStarTight", AliPID::kPion);
cutPiTight->SetPIDCut(piPIDCut);
AliRsnCutTrackQuality *cutQuality = (AliRsnCutTrackQuality*) cutPi->CutQuality();
+ cutQuality->SetAODTestFilterBit(aodFilterBit);
cutQuality->SetDCARmax(0.036);
// cut set
cutLambdaTight->SetMaxDCAVertex(0.22);
cutLambdaTight->SetMinCosPointingAngle(0.995);
cutLambdaTight->SetMaxDaughtersDCA(0.35);
- cutLambdaTight->SetMaxRapidity(0.5);
+ cutLambdaTight->SetMaxRapidity(0.8);
cutLambdaTight->SetAODTestFilterBit(aodFilterBit);
cutLambdaTight->SetPIDCut1(pPIDCut);
cutLambdaTight->SetPIDCut2(piPIDCut);
//AliRsnCutV0 *cutAntiLambdaTight = new AliRsnCutV0("cutAntiLambdaTight", kLambda0Bar);
AliRsnCutV0 *cutAntiLambdaTight = new AliRsnCutV0("cutAntiLambdaTight", kLambda0Bar, AliPID::kProton, AliPID::kPion);
cutAntiLambdaTight->SetESDtrackCuts(esdTrackCuts);
- cutAntiLambdaTight->SetTolerance(0.076);
+ cutAntiLambdaTight->SetTolerance(0.0076);
cutAntiLambdaTight->SetMaxDCAVertex(0.22);
cutAntiLambdaTight->SetMinCosPointingAngle(0.995);
cutAntiLambdaTight->SetMaxDaughtersDCA(0.35);
- cutAntiLambdaTight->SetMaxRapidity(0.5);
+ cutAntiLambdaTight->SetMaxRapidity(0.8);
cutAntiLambdaTight->SetAODTestFilterBit(aodFilterBit);
cutAntiLambdaTight->SetPIDCut1(pPIDCut);
cutAntiLambdaTight->SetPIDCut2(piPIDCut);
+++ /dev/null
-//
-// *** Configuration script for phi->KK analysis with 2010 runs ***
-//
-// A configuration script for RSN package needs to define the followings:
-//
-// (1) decay tree of each resonance to be studied, which is needed to select
-// true pairs and to assign the right mass to all candidate daughters
-// (2) cuts at all levels: single daughters, tracks, events
-// (3) output objects: histograms or trees
-//
-Bool_t ConfigSigmaStarMC
-(
- AliRsnMiniAnalysisTask *task,
- Bool_t isPP,
- const char *suffix,
- AliRsnCutSet *cutsPair
-)
-{
- // manage suffix
- if (strlen(suffix) > 0) suffix = Form("_%s", suffix);
-
- //
- // -- Define track cuts -------------------------------------------------------------------------
- //
-
- /*** EMPTY FOR TRUE PAIRS COMPUTATION ***/
-
- //
- // -- Values ------------------------------------------------------------------------------------
- //
-
- /* invariant mass */ Int_t imID = task->CreateValue(AliRsnMiniValue::kInvMass, kFALSE);
- /* transv. momentum */ Int_t ptID = task->CreateValue(AliRsnMiniValue::kPt, kFALSE);
- /* centrality */ Int_t centID = task->CreateValue(AliRsnMiniValue::kMult, kFALSE);
-
- //
- // -- Create all needed outputs -----------------------------------------------------------------
- //
-
- TString mode = "HIST";
- if (!isPP) mode = "SPARSE";
-
- // create output
- AliRsnMiniOutput *out = task->CreateOutput(Form("sigmastarP_TrueMC%s", suffix), mode.Data(), "MOTHER");
- // selection settings
- out->SetDaughter(0, AliRsnDaughter::kLambda);
- out->SetDaughter(1, AliRsnDaughter::kPion);
- out->SetMotherPDG(3224);
- out->SetMotherMass(1382.3);
- // pair cuts
- out->SetPairCuts(cutsPair);
- // binnings
- out->AddAxis(imID, 800, 1.2, 2.0);
- out->AddAxis(ptID, 100, 0.0, 10.0);
- if (!isPP) out->AddAxis(centID, 100, 0.0, 100.0);
-
- // create output
- AliRsnMiniOutput *out = task->CreateOutput(Form("sigmastarM_TrueMC%s", suffix), mode.Data(), "MOTHER");
- // selection settings
- out->SetDaughter(0, AliRsnDaughter::kLambda);
- out->SetDaughter(1, AliRsnDaughter::kPion);
- out->SetMotherPDG(3114);
- out->SetMotherMass(1387.4);
- // pair cuts
- out->SetPairCuts(cutsPair);
- // binnings
- out->AddAxis(imID, 800, 1.2, 2.0);
- out->AddAxis(ptID, 100, 0.0, 10.0);
- if (!isPP) out->AddAxis(centID, 100, 0.0, 100.0);
-
- return kTRUE;
-}