D00 cuts in 1-2 + small updates in charm fraction task (AndreaR)

diff --git a/PWGHF/vertexingHF/AliAnalysisTaskSECharmFraction.cxx b/PWGHF/vertexingHF/AliAnalysisTaskSECharmFraction.cxx
index 02bc500c04f883342f8457675b148bed1398b210..926789ab62a1a83acf169cfbde20051dc7850064 100644 (file)
--- a/PWGHF/vertexingHF/AliAnalysisTaskSECharmFraction.cxx
+++ b/PWGHF/vertexingHF/AliAnalysisTaskSECharmFraction.cxx
@@ -5638,13 +5638,16 @@ void AliAnalysisTaskSECharmFraction::UserExec(Option_t */*option*/)
   
     signallevel=-1;
     
+   
 
     AliAODRecoDecayHF2Prong *d = (AliAODRecoDecayHF2Prong*)arrayD0toKpi->UncheckedAt(iD0toKpi);
-   //  Bool_t unsetvtx=kFALSE;
-//     if(!d->GetOwnPrimaryVtx()) {
-//       d->SetOwnPrimaryVtx(vtx1); // needed to compute all variables
-//       unsetvtx=kTRUE;
-//     }
+   
+ 
+    //  Bool_t unsetvtx=kFALSE;
+    //     if(!d->GetOwnPrimaryVtx()) {
+    //       d->SetOwnPrimaryVtx(vtx1); // needed to compute all variables
+    //       unsetvtx=kTRUE;
+    //     }
 
     //recalculate vertex w/o daughters
     AliAODVertex *origownvtx=0x0;
@@ -6505,6 +6508,11 @@ AliAODRecoDecayHF* AliAnalysisTaskSECharmFraction::GetD0toKPiSignalType(const Al
     return aodDMC;    
   }
 
+  if(labMum<0){
+    signaltype=-1;
+    return aodDMC;    
+  }
+
   // get daughter AOD tracks
   AliAODTrack *trk0 = (AliAODTrack*)d->GetDaughter(0);
   AliAODTrack *trk1 = (AliAODTrack*)d->GetDaughter(1);
@@ -6514,9 +6522,14 @@ AliAODRecoDecayHF* AliAnalysisTaskSECharmFraction::GetD0toKPiSignalType(const Al
   
     signaltype=-1;
     return aodDMC;
-   
+    
   }
-
+  if(trk0->GetLabel()<0||trk1->GetLabel()<0){
+    AliDebug(2,"Fake tracks? \n");
+    signaltype=-1;
+    return aodDMC;
+  }
+  
   b1=(AliAODMCParticle*)arrayMC->At(trk0->GetLabel());
   b2=(AliAODMCParticle*)arrayMC->At(trk1->GetLabel());
   mum1=(AliAODMCParticle*)arrayMC->At(labMum);  
@@ -6529,8 +6542,11 @@ AliAODRecoDecayHF* AliAnalysisTaskSECharmFraction::GetD0toKPiSignalType(const Al
   }
 
   Bool_t isfromDstar=kFALSE;
-  grandmoth1=(AliAODMCParticle*)arrayMC->At(mum1->GetMother());
-  if(TMath::Abs(grandmoth1->GetPdgCode())==413||TMath::Abs(grandmoth1->GetPdgCode())==423)isfromDstar=kTRUE;// D0 COMING FROM A D*
+ 
+  if(mum1->GetMother()>=0){
+     grandmoth1=(AliAODMCParticle*)arrayMC->At(mum1->GetMother());
+    if(TMath::Abs(grandmoth1->GetPdgCode())==413||TMath::Abs(grandmoth1->GetPdgCode())==423)isfromDstar=kTRUE;// D0 COMING FROM A D*
+  }
   
   Int_t origin=CheckOrigin(arrayMC,mum1);
   if(origin==4){

diff --git a/PWGHF/vertexingHF/AliRDHFCutsD0toKpi.cxx b/PWGHF/vertexingHF/AliRDHFCutsD0toKpi.cxx
index 2595282c2f380c4fb44225ee7b94336c7858f043..f234eef82765d70f33fd80455fca6e3bcd44977e 100644 (file)
--- a/PWGHF/vertexingHF/AliRDHFCutsD0toKpi.cxx
+++ b/PWGHF/vertexingHF/AliRDHFCutsD0toKpi.cxx
@@ -1463,7 +1463,7 @@ void AliRDHFCutsD0toKpi::SetStandardCutsPbPb2010() {
 
   Float_t cutsMatrixD0toKpiStand[nptbins][nvars]={{0.400,400.*1E-4,0.8,0.3,0.3,1000.*1E-4,1000.*1E-4,-10000.*1E-8,0.85,0.99,2.},/* pt<0.5*/
                                                  {0.400,400.*1E-4,0.8,0.3,0.3,1000.*1E-4,1000.*1E-4,-35000.*1E-8,0.9,0.99,2.},/* 0.5<pt<1*/
-                                                 {0.400,400.*1E-4,0.8,0.4,0.4,1000.*1E-4,1000.*1E-4,-43000.*1E-8,0.85,0.995,2.},/* 1<pt<2 */
+                                                 {0.400,400.*1E-4,0.8,0.4,0.4,1000.*1E-4,1000.*1E-4,-43000.*1E-8,0.85,0.995,8.},/* 1<pt<2 */
                                                  {0.400,250.*1E-4,0.8,0.7,0.7,1000.*1E-4,1000.*1E-4,-45000.*1E-8,0.95,0.998,7.},/* 2<pt<3 */
                                                  {0.400,250.*1E-4,0.8,0.7,0.7,1000.*1E-4,1000.*1E-4,-36000.*1E-8,0.95,0.998,5.},/* 3<pt<4 */
                                                  {0.400,250.*1E-4,0.8,0.7,0.7,1000.*1E-4,1000.*1E-4,-27000.*1E-8,0.95,0.998,5.},/* 4<pt<5 */