+
+
+
+
+ ////////////////////////////////////////////////////////////////////////////////////////////////
+ // Cut set for mothers -- To be used only for the efficiency evaluation when using the charm enriched sample
+ /*
+ AliRsnCutDaughterD0 *cutPi_mother = new AliRsnCutDaughterD0("cutPionForD0_MC", AliPID::kPion);
+ cutPi_mother->SetIsCheckOnMother(kTRUE);
+ cutPi_mother->SetNoPID(kTRUE);
+ cutPi_mother->SetTPCPionPIDCut(1E20);
+ cutPi_mother->SetTOFPionPIDCut(1E20);
+ cutPi_mother->SetPtDependentPIDCut(kFALSE);
+ AliRsnCutTrackQuality *cutQuality = (AliRsnCutTrackQuality*) cutPi_mother->CutQuality();
+ cutQuality->SetCheckOnlyFilterBit(kFALSE);
+ cutQuality->SetAODTestFilterBit(-1);
+ if(maxDCAcutFixed)cutQuality->SetDCARmax(1E20);
+ else cutQuality->SetDCARPtFormula("");
+ if(minDCAcutFixed) cutQuality->SetDCARmin(0);
+ else cutQuality->SetDCARPtFormulaMin("");
+ if(!isPP)cutQuality->SetTPCminNClusters(0);
+ if(isPP)cutQuality->SetMinNCrossedRowsTPC(0,kTRUE);
+ if(isPP)cutQuality->SetMinNCrossedRowsOverFindableClsTPC(0,kTRUE);
+ cutQuality->SetPtRange(0,1E20);
+ cutQuality->SetEtaRange(-1E20, 1E20);
+ cutQuality->SetDCAZmax(1E20);
+ cutQuality->SetSPDminNClusters(0);
+ cutQuality->SetITSminNClusters(0);
+ cutQuality->SetITSmaxChi2(1E20);
+ cutQuality->SetTPCmaxChi2(1E20);
+ cutQuality->SetRejectKinkDaughters();
+ cutQuality->Print();
+
+
+ // cut set
+ AliRsnCutSet *cutSetPi_mother = new AliRsnCutSet("setPionD0_mother", AliRsnTarget::kDaughter);
+ cutSetPi_mother->AddCut(cutPi_mother);
+ cutSetPi_mother->SetCutScheme(cutPi_mother->GetName());
+ // add to task
+ Int_t iCutPi_mother = task->AddTrackCuts(cutSetPi_mother);
+
+
+ AliRsnCutDaughterD0 *cutK_mother = new AliRsnCutDaughterD0("cutKaonForD0_MC", AliPID::kKaon);
+ cutK_mother->SetIsCheckOnMother(kTRUE);
+ cutK_mother->SetNoPID(kTRUE);
+ cutK_mother->SetTPCKaonPIDCut(1E20);
+ cutK_mother->SetTOFKaonPIDCut(1E20);
+ cutK_mother->SetPtDependentPIDCut(kFALSE);
+ AliRsnCutTrackQuality *cutQuality = (AliRsnCutTrackQuality*) cutK_mother->CutQuality();
+ cutQuality->SetCheckOnlyFilterBit(kFALSE);
+ cutQuality->SetAODTestFilterBit(-1);
+ if(maxDCAcutFixed)cutQuality->SetDCARmax(1E20);
+ else cutQuality->SetDCARPtFormula("");
+ if(minDCAcutFixed) cutQuality->SetDCARmin(0);
+ else cutQuality->SetDCARPtFormulaMin("");
+ if(!isPP)cutQuality->SetTPCminNClusters(0);
+ if(isPP)cutQuality->SetMinNCrossedRowsTPC(0,kTRUE);
+ if(isPP)cutQuality->SetMinNCrossedRowsOverFindableClsTPC(0,kTRUE);
+ cutQuality->SetPtRange(0,1E20);
+ cutQuality->SetEtaRange(-1E20, 1E20);
+ cutQuality->SetDCAZmax(1E20);
+ cutQuality->SetSPDminNClusters(0);
+ cutQuality->SetITSminNClusters(0);
+ cutQuality->SetITSmaxChi2(1E20);
+ cutQuality->SetTPCmaxChi2(1E20);
+ cutQuality->SetRejectKinkDaughters();
+ cutQuality->Print();
+
+
+ // cut set
+ AliRsnCutSet *cutSetK_mother = new AliRsnCutSet("setKaonD0_mother", AliRsnTarget::kDaughter);
+ cutSetK_mother->AddCut(cutK_mother);
+ cutSetK_mother->SetCutScheme(cutK_mother->GetName());
+ // add to task
+ Int_t iCutK_mother = task->AddTrackCuts(cutSetK_mother);
+ */
+
+ ////////////////////////////////////////////////////////////////////////////////////////////////