ROOT logo
/**************************************************************************
 * Copyright(c) 1998-2010, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

/////////////////////////////////////////////////////////////
//
// Class for cuts on AOD reconstructed Lc->pKpi
//
// Author: A.Dainese, andrea.dainese@pd.infn.it
/////////////////////////////////////////////////////////////

#include <TDatabasePDG.h>
#include <Riostream.h>
#include <AliAnalysisManager.h>
#include <AliInputEventHandler.h>
#include <AliPIDResponse.h>

#include "AliRDHFCutsLctopKpi.h"
#include "AliAODRecoDecayHF3Prong.h"
#include "AliRDHFCuts.h"
#include "AliAODTrack.h"
#include "AliESDtrack.h"
#include "AliKFParticle.h"
#include "AliESDVertex.h"

using std::cout;
using std::endl;

ClassImp(AliRDHFCutsLctopKpi)

//--------------------------------------------------------------------------
AliRDHFCutsLctopKpi::AliRDHFCutsLctopKpi(const char* name) : 
AliRDHFCuts(name),
fPidObjprot(0),
fPidObjpion(0),
fUseImpParProdCorrCut(kFALSE),
fPIDStrategy(kNSigma),
fCutsStrategy(kStandard),
fUseSpecialCut(kFALSE)
{
  //
  // Default Constructor
  //
  Int_t nvars=13;
  SetNVars(nvars);
  TString varNames[13]={"inv. mass [GeV]",
			"pTK [GeV/c]",
			"pTP [GeV/c]",
			"d0K [cm]   lower limit!",
			"d0Pi [cm]  lower limit!",
			"dist12 (cm)",
			"sigmavert (cm)",
			"dist prim-sec (cm)",
			"pM=Max{pT1,pT2,pT3} (GeV/c)",
			"cosThetaPoint",
			"Sum d0^2 (cm^2)",
			"dca cut (cm)",
			"cut on pTpion [GeV/c]"};
  Bool_t isUpperCut[13]={kTRUE,
			 kFALSE,
			 kFALSE,
			 kFALSE,
			 kFALSE,
			 kFALSE,
			 kTRUE,
			 kFALSE,
			 kFALSE,
			 kFALSE,
			 kFALSE,
			 kTRUE,
			 kFALSE
			 };
  SetVarNames(nvars,varNames,isUpperCut);
  Bool_t forOpt[13]={kFALSE,
		     kTRUE,
		     kTRUE,
		     kFALSE,
		     kFALSE,
		     kFALSE,
		     kFALSE,
		     kTRUE,
		     kFALSE,
		     kFALSE,
		     kFALSE,
		     kFALSE,
		     kTRUE};
  SetVarsForOpt(4,forOpt);
  Float_t limits[2]={0,999999999.};
  SetPtBins(2,limits);
  for (Int_t ispecies=0;ispecies<AliPID::kSPECIES;++ispecies)
      fPIDThreshold[ispecies]=0.;
}
//--------------------------------------------------------------------------
AliRDHFCutsLctopKpi::AliRDHFCutsLctopKpi(const AliRDHFCutsLctopKpi &source) :
  AliRDHFCuts(source),
  fPidObjprot(0x0),
  fPidObjpion(0x0),
  fUseImpParProdCorrCut(source.fUseImpParProdCorrCut),
  fPIDStrategy(source.fPIDStrategy),
  fCutsStrategy(source.fCutsStrategy),
  fUseSpecialCut(source.fUseSpecialCut)
{
  //
  // Copy constructor
  //
  if (source.fPidObjprot) fPidObjprot = new AliAODPidHF(*(source.fPidObjprot));
  else fPidObjprot = new AliAODPidHF();
  if (source.fPidObjpion) fPidObjpion = new AliAODPidHF(*(source.fPidObjpion));
  else fPidObjpion = new AliAODPidHF();
  memcpy(fPIDThreshold,source.fPIDThreshold,AliPID::kSPECIES*sizeof(Double_t));
}
//--------------------------------------------------------------------------
AliRDHFCutsLctopKpi &AliRDHFCutsLctopKpi::operator=(const AliRDHFCutsLctopKpi &source)
{
  //
  // assignment operator
  //
  if(this != &source) {
    
    AliRDHFCuts::operator=(source);
    delete fPidObjprot;
    fPidObjprot = new AliAODPidHF(*(source.fPidObjprot));
    delete fPidObjpion;
    fPidObjpion = new AliAODPidHF(*(source.fPidObjpion));
    fPIDStrategy=source.fPIDStrategy;
    fCutsStrategy=source.fCutsStrategy;
    memcpy(fPIDThreshold,source.fPIDThreshold,AliPID::kSPECIES*sizeof(Double_t));
  }
    
  return *this;
}
//---------------------------------------------------------------------------
AliRDHFCutsLctopKpi::~AliRDHFCutsLctopKpi() {
 //
 //  // Default Destructor
 //   
 if(fPidObjpion){
  delete fPidObjpion;
  fPidObjpion=0;
 }
 if(fPidObjprot){
  delete fPidObjprot;
  fPidObjprot=0;
 }

}

//---------------------------------------------------------------------------
void AliRDHFCutsLctopKpi::GetCutVarsForOpt(AliAODRecoDecayHF *d,Float_t *vars,Int_t nvars,Int_t *pdgdaughters, AliAODEvent *aod) {
  // 
  // Fills in vars the values of the variables 
  //

  if(nvars!=fnVarsForOpt) {
    printf("AliRDHFCutsLctopKpi::GetCutsVarsForOpt: wrong number of variables\n");
    return;
  }

  AliAODRecoDecayHF3Prong *dd = (AliAODRecoDecayHF3Prong*)d;

  Int_t iter=-1;
  if(fVarsForOpt[0]){
    iter++;
    vars[iter]=dd->InvMassLcpKpi();
  }
  if(fVarsForOpt[1]){
    iter++;
    for(Int_t iprong=0;iprong<3;iprong++){
      if(TMath::Abs(pdgdaughters[iprong])==321) {
	vars[iter]=dd->PtProng(iprong);
      }
    }
  }
  if(fVarsForOpt[2]){
    iter++;
    for(Int_t iprong=0;iprong<3;iprong++){
      if(TMath::Abs(pdgdaughters[iprong])==2212) {
	vars[iter]=dd->PtProng(iprong);
      }
    }
  }
  if(fVarsForOpt[3]){
    iter++;
    for(Int_t iprong=0;iprong<3;iprong++){
      if(TMath::Abs(pdgdaughters[iprong])==2212) {
	vars[iter]=dd->Getd0Prong(iprong);
      }
    }
  }
  if(fVarsForOpt[4]){
    iter++;
    for(Int_t iprong=0;iprong<3;iprong++){
      if(TMath::Abs(pdgdaughters[iprong])==211) {
	vars[iter]=dd->Getd0Prong(iprong);
      }
    }
  }
  if(fVarsForOpt[5]){
    iter++;
    vars[iter]=dd->GetDist12toPrim();
  }
  if(fVarsForOpt[6]){
    iter++;
    vars[iter]=dd->GetSigmaVert(aod);
  }
  if(fVarsForOpt[7]){
    iter++;
    vars[iter] = dd->DecayLength();
  }
  if(fVarsForOpt[8]){
    iter++;
    Float_t ptmax=0;
    for(Int_t i=0;i<3;i++){
      if(dd->PtProng(i)>ptmax)ptmax=dd->PtProng(i);
    }
    vars[iter]=ptmax;
  }
  if(fVarsForOpt[9]){
    iter++;
    vars[iter]=dd->CosPointingAngle();
  }
  if(fVarsForOpt[10]){
    iter++;
    vars[iter]=dd->Getd0Prong(0)*dd->Getd0Prong(0)+dd->Getd0Prong(1)*dd->Getd0Prong(1)+dd->Getd0Prong(2)*dd->Getd0Prong(2);
  }
  if(fVarsForOpt[11]){
    iter++;
    vars[iter]=dd->GetDCA();
  }
  if(fVarsForOpt[12]){
    iter++;
    for(Int_t iprong=0;iprong<3;iprong++){
      if(TMath::Abs(pdgdaughters[iprong])==211) {
	vars[iter]=dd->PtProng(iprong);
      }
    }
  }

  return;
}
//---------------------------------------------------------------------------
Int_t AliRDHFCutsLctopKpi::IsSelected(TObject* obj,Int_t selectionLevel,AliAODEvent *aod) {
  //
  // Apply selection
  //

  if(!fCutsRD){
    AliError("Cut matrice not inizialized. Exit...\n");
    return 0;
  }
  //PrintAll();
  AliAODRecoDecayHF3Prong* d=(AliAODRecoDecayHF3Prong*)obj;

  if(!d){
    AliError("AliAODRecoDecayHF3Prong null \n");
    return 0;
  }


  if(fKeepSignalMC) if(IsSignalMC(d,aod,4122)) return 3;

  Int_t returnvalue=3;
  Int_t returnvaluePID=3;

  if(d->Pt()<fMinPtCand) return 0;
  if(d->Pt()>fMaxPtCand) return 0;

  if(fUseTrackSelectionWithFilterBits && d->HasBadDaughters()) return 0;


  // selection on daughter tracks 
  if(selectionLevel==AliRDHFCuts::kAll || 
     selectionLevel==AliRDHFCuts::kTracks) {
    if(!AreDaughtersSelected(d)) return 0;
  }


  // PID selection
  if(selectionLevel==AliRDHFCuts::kAll ||
     selectionLevel==AliRDHFCuts::kCandidate|| 
     selectionLevel==AliRDHFCuts::kPID) {
    switch (fPIDStrategy) {
    case kNSigma:
      returnvaluePID = IsSelectedPID(d);
      break;
     case kNSigmaPbPb:
      returnvaluePID = IsSelectedNSigmaPbPb(d);
      break;
    case kCombined:
      returnvaluePID = IsSelectedCombinedPID(d);
      break;
    case kCombinedSoft:
      returnvaluePID = IsSelectedCombinedPIDSoft(d);
      break;
    case kNSigmaStrong:
      returnvaluePID = IsSelectedPIDStrong(d);
      break;
    case kCombinedpPb:
      returnvaluePID = IsSelectedCombinedPIDpPb(d);
      break;
    case kCombinedpPb2:
      returnvaluePID = IsSelectedCombinedPIDpPb2(d);
      break;
    }
    fIsSelectedPID=returnvaluePID;
  }
  //  if(fUsePID || selectionLevel==AliRDHFCuts::kPID) returnvaluePID = IsSelectedCombinedPID(d);   // to test!!
  if(returnvaluePID==0) return 0;




  // selection on candidate
  if(selectionLevel==AliRDHFCuts::kAll || 
     selectionLevel==AliRDHFCuts::kCandidate) {

    Double_t pt=d->Pt();
    
    Int_t ptbin=PtBin(pt);
    
    Double_t mLcpKpi=0.,mLcpiKp=0.;
    Int_t okLcpKpi=1,okLcpiKp=1;

    Double_t mLcPDG = TDatabasePDG::Instance()->GetParticle(4122)->Mass();

    mLcpKpi=d->InvMassLcpKpi();
    mLcpiKp=d->InvMassLcpiKp();

    if(TMath::Abs(mLcpKpi-mLcPDG)>fCutsRD[GetGlobalIndex(0,ptbin)]) okLcpKpi = 0;
    if(TMath::Abs(mLcpiKp-mLcPDG)>fCutsRD[GetGlobalIndex(0,ptbin)]) okLcpiKp = 0;
    if(!okLcpKpi && !okLcpiKp) return 0;

  switch (fCutsStrategy) {

    case kStandard:
    if(TMath::Abs(d->PtProng(1)) < fCutsRD[GetGlobalIndex(1,ptbin)] || TMath::Abs(d->Getd0Prong(1))<fCutsRD[GetGlobalIndex(3,ptbin)]) return 0;//Kaon
    if(d->Pt()>=3. && d->PProng(1)<0.55) return 0;
    if(fUseSpecialCut) {
      if(TMath::Abs(d->PtProng(0)) < TMath::Abs(d->PtProng(2)) )okLcpKpi=0;
      if(TMath::Abs(d->PtProng(2)) < TMath::Abs(d->PtProng(0)) )okLcpiKp=0;
    }
    if((TMath::Abs(d->PtProng(0)) < fCutsRD[GetGlobalIndex(2,ptbin)]) || (TMath::Abs(d->PtProng(2)) < fCutsRD[GetGlobalIndex(12,ptbin)])) okLcpKpi=0;
    if((TMath::Abs(d->PtProng(2)) < fCutsRD[GetGlobalIndex(2,ptbin)]) || (TMath::Abs(d->PtProng(0)) < fCutsRD[GetGlobalIndex(12,ptbin)]))okLcpiKp=0;
    if(!okLcpKpi && !okLcpiKp) return 0;
    //2track cuts
    if(d->GetDist12toPrim()<fCutsRD[GetGlobalIndex(5,ptbin)]|| d->GetDist23toPrim()<fCutsRD[GetGlobalIndex(5,ptbin)]) return 0;
    if(d->GetDist12toPrim()>0.5) return 0;
    if(d->GetDist23toPrim()>0.5) return 0;
    if(fUseImpParProdCorrCut){
      if(d->Getd0Prong(0)*d->Getd0Prong(1)<0. && d->Getd0Prong(2)*d->Getd0Prong(1)<0.) return 0;
    }
    //sec vert
    if(d->DecayLength()<fCutsRD[GetGlobalIndex(7,ptbin)]) return 0;
    if(d->DecayLength()>0.5) return 0;

  //  Double_t sumd0s=d->Getd0Prong(0)*d->Getd0Prong(0)+d->Getd0Prong(1)*d->Getd0Prong(1)+d->Getd0Prong(2)*d->Getd0Prong(2);
  //  if(sumd0s<fCutsRD[GetGlobalIndex(10,ptbin)]) return 0;
    if((d->Getd0Prong(0)*d->Getd0Prong(0)+d->Getd0Prong(1)*d->Getd0Prong(1)+d->Getd0Prong(2)*d->Getd0Prong(2))<fCutsRD[GetGlobalIndex(10,ptbin)]) return 0;
    
    if(TMath::Abs(d->PtProng(0))<fCutsRD[GetGlobalIndex(8,ptbin)] && TMath::Abs(d->PtProng(1))<fCutsRD[GetGlobalIndex(8,ptbin)] && TMath::Abs(d->PtProng(2))<fCutsRD[GetGlobalIndex(8,ptbin)]) return 0;
    if(d->CosPointingAngle()< fCutsRD[GetGlobalIndex(9,ptbin)]) return 0;
    if(d->GetSigmaVert(aod)>fCutsRD[GetGlobalIndex(6,ptbin)]) return 0;
    
    //DCA
    for(Int_t i=0;i<3;i++) if(d->GetDCA(i)>fCutsRD[GetGlobalIndex(11,ptbin)]) return 0;

    break;

   case kKF:
    Int_t pdgs[3]={0,321,0};
    Bool_t constraint=kFALSE;
    if(fCutsRD[GetGlobalIndex(1,ptbin)]>0.) constraint=kTRUE;
    Double_t field=aod->GetMagneticField();
    if (returnvaluePID==1 || returnvaluePID==3){

      pdgs[0]=2122;pdgs[2]=211;
      AliKFParticle *lc1=ReconstructKF(d,pdgs,field,constraint);
      if(!lc1){
	okLcpKpi=0;
      }else{
	if(lc1->GetChi2()/lc1->GetNDF()>fCutsRD[GetGlobalIndex(2,ptbin)]) okLcpKpi=0;
      }
    } else if(returnvaluePID>=2){

      pdgs[0]=211;pdgs[2]=2212;
      AliKFParticle *lc2=ReconstructKF(d,pdgs,field,constraint);
      if(!lc2){ 
	okLcpiKp=0;
      }else{
	if(lc2->GetChi2()/lc2->GetNDF()>fCutsRD[GetGlobalIndex(2,ptbin)])okLcpiKp=0; 
      }
    }
    break;

   }

    if(okLcpKpi) returnvalue=1; //cuts passed as Lc->pKpi
    if(okLcpiKp) returnvalue=2; //cuts passed as Lc->piKp
    if(okLcpKpi && okLcpiKp) returnvalue=3; //cuts passed as both pKpi and piKp
   
  }


  Int_t returnvalueTot=CombinePIDCuts(returnvalue,returnvaluePID);
  return returnvalueTot;
}
//---------------------------------------------------------------------------
Int_t AliRDHFCutsLctopKpi::IsSelectedPID(AliAODRecoDecayHF* obj) {


    if(!fUsePID || !obj) return 3;
    Int_t okLcpKpi=0,okLcpiKp=0;
    Int_t returnvalue=0;
    Bool_t isPeriodd=fPidHF->GetOnePad();
    Bool_t isMC=fPidHF->GetMC();
    Bool_t ispion0=kTRUE,ispion2=kTRUE;
    Bool_t isproton0=kFALSE,isproton2=kFALSE;
    Bool_t iskaon1=kFALSE;
    if(isPeriodd) {
     fPidObjprot->SetOnePad(kTRUE);
     fPidObjpion->SetOnePad(kTRUE);
    }
    if(isMC) {
     fPidObjprot->SetMC(kTRUE);
     fPidObjpion->SetMC(kTRUE);
    }

   if(fPidObjprot->GetPidResponse()==0x0){
      AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
      AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
      AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
      fPidObjprot->SetPidResponse(pidResp);
    }
    if(fPidObjpion->GetPidResponse()==0x0){
      AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
      AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
      AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
      fPidObjpion->SetPidResponse(pidResp);
    }
    if(fPidHF->GetPidResponse()==0x0){
      AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
      AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
      AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
      fPidHF->SetPidResponse(pidResp);
    }

    for(Int_t i=0;i<3;i++){
     AliAODTrack *track=(AliAODTrack*)obj->GetDaughter(i);
     if(!track) return 0;
     // identify kaon
     if(track->P()<0.55){
      fPidHF->SetTOF(kFALSE);
      fPidHF->SetTOFdecide(kFALSE);
     }
     if(i==1) {
      Int_t isKaon=fPidHF->MakeRawPid(track,3); 
      if(isKaon>=1) iskaon1=kTRUE;
      if(track->P()<0.55){
       fPidHF->SetTOF(kTRUE);
       fPidHF->SetTOFdecide(kTRUE);
      }
      
      if(!iskaon1) return 0;
     
     }else{
     //pion or proton
    if(track->P()<1.){
      fPidObjprot->SetTOF(kFALSE);
      fPidObjprot->SetTOFdecide(kFALSE);
     }
      
     Int_t isProton=fPidObjprot->MakeRawPid(track,4);
     

     Int_t isPion=fPidObjpion->MakeRawPid(track,2);
     
     if(track->P()<1.){
      fPidObjprot->SetTOF(kTRUE);
      fPidObjprot->SetTOFdecide(kTRUE);
     }
     

     if(i==0) {
      if(isPion<0) ispion0=kFALSE;
      if(isProton>=1) isproton0=kTRUE;

     }
      if(!ispion0 && !isproton0) return 0;
     if(i==2) {
      if(isPion<0) ispion2=kFALSE;
      if(isProton>=1) isproton2=kTRUE;
     }

    }
   }

    if(ispion2 && isproton0 && iskaon1) okLcpKpi=1;
    if(ispion0 && isproton2 && iskaon1) okLcpiKp=1;
    if(okLcpKpi) returnvalue=1; //cuts passed as Lc->pKpi
    if(okLcpiKp) returnvalue=2; //cuts passed as Lc->piKp
    if(okLcpKpi && okLcpiKp) returnvalue=3; //cuts passed as both pKpi and piKp

 return returnvalue;
}

//--------------------------------------------------------------------------

Int_t AliRDHFCutsLctopKpi::IsSelectedNSigmaPbPb(AliAODRecoDecayHF* obj) {


  if(!fUsePID || !obj) return 3;
  Int_t okLcpKpi=0,okLcpiKp=0;
  Int_t returnvalue=0;
  Bool_t isPeriodd=fPidHF->GetOnePad();
  Bool_t isMC=fPidHF->GetMC();
  Bool_t ispion0=kTRUE,ispion2=kTRUE;
  Bool_t isproton0=kFALSE,isproton2=kFALSE;
  Bool_t iskaon1=kFALSE;
  if(isPeriodd) {
    fPidObjprot->SetOnePad(kTRUE);
    fPidObjpion->SetOnePad(kTRUE);
  }
  if(isMC) {
    fPidObjprot->SetMC(kTRUE);
    fPidObjpion->SetMC(kTRUE);
  }

  if(fPidObjprot->GetPidResponse()==0x0){
    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
    AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
    AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
    fPidObjprot->SetPidResponse(pidResp);
  }
  if(fPidObjpion->GetPidResponse()==0x0){
    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
    AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
    AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
    fPidObjpion->SetPidResponse(pidResp);
  }
  if(fPidHF->GetPidResponse()==0x0){
    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
    AliInputEventHandler *inputHandler=(AliInputEventHandler*)mgr->GetInputEventHandler();
    AliPIDResponse *pidResp=inputHandler->GetPIDResponse();
    fPidHF->SetPidResponse(pidResp);
  }

  for(Int_t i=0;i<3;i++){
    AliAODTrack *track=(AliAODTrack*)obj->GetDaughter(i);
    if(!track) return 0;

    if(i==1) {
      //kaon
      Int_t isKaon=fPidHF->MakeRawPid(track,3); 
      if(isKaon>=1) iskaon1=kTRUE;
      if(!iskaon1) return 0;
    }

    else {
      //pion or proton
      Int_t isProton=fPidObjprot->MakeRawPid(track,4);
      Int_t isPion=fPidObjpion->MakeRawPid(track,2);

      if(i==0) {
	if(isPion<0) ispion0=kFALSE;
	if(isProton>=1) isproton0=kTRUE;
      }
      if(!ispion0 && !isproton0) return 0;

      if(i==2) {
	if(isPion<0) ispion2=kFALSE;
	if(isProton>=1) isproton2=kTRUE;
      }
    }
  }

  if(ispion2 && isproton0 && iskaon1) okLcpKpi=1;
  if(ispion0 && isproton2 && iskaon1) okLcpiKp=1;
  if(okLcpKpi) returnvalue=1; //cuts passed as Lc->pKpi
  if(okLcpiKp) returnvalue=2; //cuts passed as Lc->piKp
  if(okLcpKpi && okLcpiKp) returnvalue=3; //cuts passed as both pKpi and piKp

  return returnvalue;
}




//---------------------------------------------------------------------------
Int_t AliRDHFCutsLctopKpi::IsSelectedCombinedPID(AliAODRecoDecayHF* obj) {
    
    if(!fUsePID || !obj) {return 3;}
    Int_t okLcpKpi=0,okLcpiKp=0;
    Int_t returnvalue=0;
    Bool_t isPeriodd=fPidHF->GetOnePad();
    Bool_t isMC=fPidHF->GetMC();

    if(isPeriodd) {
	    fPidObjprot->SetOnePad(kTRUE);
	    fPidObjpion->SetOnePad(kTRUE);
    }
    if(isMC) {
	    fPidObjprot->SetMC(kTRUE);
	    fPidObjpion->SetMC(kTRUE);
    }

    AliVTrack *track0=dynamic_cast<AliVTrack*>(obj->GetDaughter(0));
    AliVTrack *track1=dynamic_cast<AliVTrack*>(obj->GetDaughter(1));
    AliVTrack *track2=dynamic_cast<AliVTrack*>(obj->GetDaughter(2));
    if (!track0 || !track1 || !track2) return 0;
    Double_t prob0[AliPID::kSPECIES];
    Double_t prob1[AliPID::kSPECIES];
    Double_t prob2[AliPID::kSPECIES];
    if(obj->Pt()<3. && track0->P()<1.) fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC);
    fPidHF->GetPidCombined()->ComputeProbabilities(track0,fPidHF->GetPidResponse(),prob0);
    if(obj->Pt()<3. && track0->P()<1.) fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC|AliPIDResponse::kDetTOF);

   if(obj->Pt()<3. && track1->P()<0.55) fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC);
    fPidHF->GetPidCombined()->ComputeProbabilities(track1,fPidHF->GetPidResponse(),prob1);
   if(obj->Pt()<3. && track1->P()<0.55) fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC|AliPIDResponse::kDetTOF);

    if(obj->Pt()<3. && track2->P()<1.) fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC);
    fPidHF->GetPidCombined()->ComputeProbabilities(track2,fPidHF->GetPidResponse(),prob2);
   if(obj->Pt()<3. && track2->P()<1.) fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC|AliPIDResponse::kDetTOF);

    if(fPIDThreshold[AliPID::kPion]>0. && fPIDThreshold[AliPID::kKaon]>0. && fPIDThreshold[AliPID::kProton]>0.){
    okLcpiKp=  (prob0[AliPID::kPion  ]>fPIDThreshold[AliPID::kPion  ])
             &&(prob1[AliPID::kKaon  ]>fPIDThreshold[AliPID::kKaon  ])
             &&(prob2[AliPID::kProton]>fPIDThreshold[AliPID::kProton]);
    okLcpKpi=  (prob0[AliPID::kProton]>fPIDThreshold[AliPID::kProton])
             &&(prob1[AliPID::kKaon  ]>fPIDThreshold[AliPID::kKaon  ])
             &&(prob2[AliPID::kPion  ]>fPIDThreshold[AliPID::kPion  ]);
   }else{ 
		    //pion or proton
		    
		    
    if(TMath::MaxElement(AliPID::kSPECIES,prob1) == prob1[AliPID::kKaon]){
    if(TMath::MaxElement(AliPID::kSPECIES,prob0) == prob0[AliPID::kProton] && TMath::MaxElement(AliPID::kSPECIES,prob2) == prob2[AliPID::kPion]) okLcpKpi = 1;  
    if(TMath::MaxElement(AliPID::kSPECIES,prob2) == prob2[AliPID::kProton] && TMath::MaxElement(AliPID::kSPECIES,prob0) == prob0[AliPID::kPion]) okLcpiKp = 1; 
	    }
	   }
    
    if(okLcpKpi) returnvalue=1; //cuts passed as Lc->pKpi
    if(okLcpiKp) returnvalue=2; //cuts passed as Lc->piKp
    if(okLcpKpi && okLcpiKp) returnvalue=3; //cuts passed as both pKpi and piKp
    
    return returnvalue;
}
//-----------------------
Int_t AliRDHFCutsLctopKpi::CombinePIDCuts(Int_t returnvalue, Int_t returnvaluePID) const {

 Int_t returnvalueTot=0;
 Int_t okLcpKpi=0,okLcpiKp=0;
 if(returnvaluePID==1){
   if(returnvalue==1 || returnvalue==3) okLcpKpi=1;
 }
 if(returnvaluePID==2){
   if(returnvalue>=2) okLcpiKp=1;
 }
 if(returnvaluePID==3 && returnvalue>0){
  if(returnvalue==1 || returnvalue==3) okLcpKpi=1;
  if(returnvalue>=2) okLcpiKp=1;
 } 

 if(okLcpKpi) returnvalueTot=1; //cuts passed as Lc->pKpi
 if(okLcpiKp) returnvalueTot=2; //cuts passed as Lc->piKp
 if(okLcpKpi && okLcpiKp) returnvalueTot=3; //cuts passed as both pKpi and piKp
 return returnvalueTot;
}
//----------------------------------
void AliRDHFCutsLctopKpi::SetStandardCutsPP2010() {

 SetName("LctopKpiProdCuts");
 SetTitle("Production cuts for Lc analysis");

 AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default");
 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
 esdTrackCuts->SetRequireTPCRefit(kTRUE);
 esdTrackCuts->SetMinNClustersTPC(70);
 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
                                          AliESDtrackCuts::kAny);
 esdTrackCuts->SetRequireITSRefit(kTRUE);
 esdTrackCuts->SetMinNClustersITS(4);
 esdTrackCuts->SetMinDCAToVertexXY(0.);
 esdTrackCuts->SetEtaRange(-0.8,0.8);
 esdTrackCuts->SetPtRange(0.3,1.e10);
 AddTrackCuts(esdTrackCuts);
 delete esdTrackCuts;
 esdTrackCuts=NULL;

 const Int_t nptbins=8;
 const Int_t nvars=13;
 Float_t* ptbins;
 ptbins=new Float_t[nptbins+1];
 
 ptbins[0]=1.;
 ptbins[1]=2.;
 ptbins[2]=3.;
 ptbins[3]=4.;
 ptbins[4]=5.;
 ptbins[5]=6.;
 ptbins[6]=8.;
 ptbins[7]=10.;
 ptbins[8]=20.;

 SetGlobalIndex(nvars,nptbins);
 SetPtBins(nptbins+1,ptbins);

 Float_t** prodcutsval;
 prodcutsval=new Float_t*[nvars];
 for(Int_t iv=0;iv<nvars;iv++){
  prodcutsval[iv]=new Float_t[nptbins];
 }

 for(Int_t ipt=0;ipt<nptbins;ipt++){
  prodcutsval[0][ipt]=0.13;
  prodcutsval[1][ipt]=0.4;
  prodcutsval[2][ipt]=0.4;
  prodcutsval[3][ipt]=0.;
  prodcutsval[4][ipt]=0.;
  prodcutsval[5][ipt]=0.;
  prodcutsval[6][ipt]=0.06;
  prodcutsval[7][ipt]=0.005;
  prodcutsval[8][ipt]=0.;
  prodcutsval[9][ipt]=0.;
  prodcutsval[10][ipt]=0.;
  prodcutsval[11][ipt]=0.05;
  prodcutsval[12][ipt]=0.4;
 }
 SetCuts(nvars,nptbins,prodcutsval);

 AliAODPidHF* pidObjK=new AliAODPidHF();
 Double_t sigmasK[5]={3.,1.,1.,3.,2.};
 pidObjK->SetSigma(sigmasK);
 pidObjK->SetAsym(kTRUE);
 pidObjK->SetMatch(1);
 pidObjK->SetTPC(kTRUE);
 pidObjK->SetTOF(kTRUE);
 pidObjK->SetITS(kTRUE);
 Double_t plimK[2]={0.5,0.8};
 pidObjK->SetPLimit(plimK,2);
 pidObjK->SetTOFdecide(kTRUE);

 SetPidHF(pidObjK);

 AliAODPidHF* pidObjpi=new AliAODPidHF();
 pidObjpi->SetTPC(kTRUE);
 Double_t sigmaspi[5]={3.,0.,0.,0.,0.};
 pidObjpi->SetSigma(sigmaspi);
 pidObjpi->SetTOFdecide(kTRUE);
 SetPidpion(pidObjpi);

 AliAODPidHF* pidObjp=new AliAODPidHF();
 Double_t sigmasp[5]={3.,1.,1.,3.,2.};
 pidObjp->SetSigma(sigmasp);
 pidObjp->SetAsym(kTRUE);
 pidObjp->SetMatch(1);
 pidObjp->SetTPC(kTRUE);
 pidObjp->SetTOF(kTRUE);
 pidObjp->SetITS(kTRUE);
 Double_t plimp[2]={1.,2.};
 pidObjp->SetPLimit(plimp,2);
 pidObjp->SetTOFdecide(kTRUE);

 SetPidprot(pidObjp);

 SetUsePID(kTRUE);
 SetOptPileup(kTRUE);
 
 PrintAll();

 for(Int_t iiv=0;iiv<nvars;iiv++){
  delete [] prodcutsval[iiv];
 }
 delete [] prodcutsval;
 prodcutsval=NULL;
 delete [] ptbins;
 ptbins=NULL;

 delete pidObjK;
 pidObjK=NULL;
 delete pidObjpi;
 pidObjpi=NULL;
 delete pidObjp;
 pidObjp=NULL;

 return;
}
//------------------
void AliRDHFCutsLctopKpi::SetStandardCutsPbPb2010() {

 SetName("LctopKpiProdCuts");
 SetTitle("Production cuts for Lc analysis");

 AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default");

 esdTrackCuts->SetRequireTPCRefit(kTRUE);
 esdTrackCuts->SetMinNClustersTPC(70);
 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
                                          AliESDtrackCuts::kAny);
 esdTrackCuts->SetRequireITSRefit(kTRUE);
 esdTrackCuts->SetMinNClustersITS(4);
 esdTrackCuts->SetMinDCAToVertexXYPtDep("0.0100*TMath::Max(0.,(1-TMath::Floor(TMath::Abs(pt)/2.)))");
 esdTrackCuts->SetEtaRange(-0.8,0.8);
 esdTrackCuts->SetMaxDCAToVertexXY(1.);
 esdTrackCuts->SetMaxDCAToVertexZ(1.);
 esdTrackCuts->SetPtRange(0.49,1.e10);
 AddTrackCuts(esdTrackCuts);
 delete esdTrackCuts;
 esdTrackCuts=NULL;

 const Int_t nptbins=8;
 const Int_t nvars=13;
 Float_t* ptbins;
 ptbins=new Float_t[nptbins+1];
 
 ptbins[0]=1.;
 ptbins[1]=2.;
 ptbins[2]=3.;
 ptbins[3]=4.;
 ptbins[4]=5.;
 ptbins[5]=6.;
 ptbins[6]=8.;
 ptbins[7]=10.;
 ptbins[8]=99999.;


 SetGlobalIndex(nvars,nptbins);
 SetPtBins(nptbins+1,ptbins);

 Float_t** prodcutsval;
 prodcutsval=new Float_t*[nvars];
 for(Int_t iv=0;iv<nvars;iv++){
  prodcutsval[iv]=new Float_t[nptbins];
 }

 for(Int_t ipt=0;ipt<nptbins;ipt++){
  prodcutsval[0][ipt]=0.13;
  prodcutsval[1][ipt]=0.5;
  prodcutsval[2][ipt]=0.6;
  prodcutsval[3][ipt]=0.;
  prodcutsval[4][ipt]=0.;
  prodcutsval[5][ipt]=0.01;
  prodcutsval[6][ipt]=0.04;
  prodcutsval[7][ipt]=0.006;
  prodcutsval[8][ipt]=0.8;
  prodcutsval[9][ipt]=0.3;
  prodcutsval[10][ipt]=0.;
  prodcutsval[11][ipt]=0.05;
  prodcutsval[12][ipt]=0.4;
 }
 SetCuts(nvars,nptbins,prodcutsval);

 AliAODPidHF* pidObj=new AliAODPidHF();
 pidObj->SetTPC(kTRUE);
 pidObj->SetTOF(kTRUE);
 SetPidHF(pidObj);
 SetPidpion(pidObj);
 SetPidprot(pidObj);


 // bayesian pid
 GetPidHF()->SetUseCombined(kTRUE);
 GetPidHF()->SetUseDefaultPriors(kTRUE);
 GetPidHF()->SetCombDetectors(AliAODPidHF::kTPCTOF);
 for (Int_t ispecies=0;ispecies<AliPID::kSPECIES;++ispecies){
   SetPIDThreshold(static_cast<AliPID::EParticleType>(ispecies),0);
 }
 SetPIDStrategy(AliRDHFCutsLctopKpi::kCombinedpPb);
 SetUsePID(kTRUE);


 PrintAll();

 for(Int_t iiv=0;iiv<nvars;iiv++){
  delete [] prodcutsval[iiv];
 }
 delete [] prodcutsval;
 prodcutsval=NULL;
 delete [] ptbins;
 ptbins=NULL;

 delete pidObj;
 pidObj=NULL;

 return;
}
//------------------
AliKFParticle* AliRDHFCutsLctopKpi::ReconstructKF(AliAODRecoDecayHF3Prong *d,Int_t *pdgs,Double_t field,Bool_t constraint) const{
  // Method to construct the KF particle from the candidate

 const Int_t nprongs=d->GetNProngs();
 if(nprongs<=0) return 0x0;

 Int_t iprongs[nprongs];
 for(Int_t i=0;i<nprongs;i++) iprongs[i]=i;

 Double_t mass[2]={0.,0.};

 AliKFParticle *decay=d->ApplyVertexingKF(iprongs,nprongs,pdgs,constraint,field,mass);
 if(!decay) return 0x0;
 AliESDVertex *vertexESD = new AliESDVertex(decay->Parameters(),
                                            decay->CovarianceMatrix(),
					    decay->GetChi2(),
					    nprongs);
 Double_t pos[3],cov[6],chi2perNDF;
 vertexESD->GetXYZ(pos);
 vertexESD->GetCovMatrix(cov);
 chi2perNDF = vertexESD->GetChi2toNDF();
 delete vertexESD; vertexESD=NULL;
 AliAODVertex *vertexAOD = new AliAODVertex(pos,cov,chi2perNDF,0x0,-1,AliAODVertex::kUndef,nprongs);
 d->SetSecondaryVtx(vertexAOD);
 return decay;
}

//------------------
void AliRDHFCutsLctopKpi::SetStandardCutsPbPb2011() {

  // Default 2010 PbPb cut object
  SetStandardCutsPbPb2010();

  //
  // Enable all 2011 PbPb run triggers
  //  
  SetTriggerClass("");
  ResetMaskAndEnableMBTrigger();
  EnableCentralTrigger();
  EnableSemiCentralTrigger();
}
//-----------------

Bool_t AliRDHFCutsLctopKpi::IsInFiducialAcceptance(Double_t pt, Double_t y) const
{
  //
  //  // Checking if Lc is in fiducial acceptance region 
  //    //
  //
  if(fMaxRapidityCand>-998.){
    if(TMath::Abs(y) > fMaxRapidityCand) return kFALSE;
    else return kTRUE;
  }

  if(pt > 5.) {
    // applying cut for pt > 5 GeV
   AliDebug(2,Form("pt of Lc = %f (> 5), cutting at |y| < 0.8",pt));
   if (TMath::Abs(y) > 0.8) return kFALSE;
  
  } else {
   // appliying smooth cut for pt < 5 GeV
   Double_t maxFiducialY = -0.2/15*pt*pt+1.9/15*pt+0.5;
   Double_t minFiducialY = 0.2/15*pt*pt-1.9/15*pt-0.5;
  AliDebug(2,Form("pt of Lc = %f (< 5), cutting  according to the fiducial zone [%f, %f]\n",pt,minFiducialY,maxFiducialY));
  if (y < minFiducialY || y > maxFiducialY) return kFALSE;
 }
  //
  return kTRUE;
}
//--------------------------------------------------------
Int_t AliRDHFCutsLctopKpi::IsSelectedCombinedPIDSoft(AliAODRecoDecayHF* obj) {
 if(!fUsePID || !obj) {return 3;}
 Int_t okLcpKpi=0,okLcpiKp=0;
 Int_t returnvalue=0;
 
 AliVTrack *track0=dynamic_cast<AliVTrack*>(obj->GetDaughter(0));
 AliVTrack *track1=dynamic_cast<AliVTrack*>(obj->GetDaughter(1));
 AliVTrack *track2=dynamic_cast<AliVTrack*>(obj->GetDaughter(2));
 if (!track0 || !track1 || !track2) return 0;
 Double_t prob0[AliPID::kSPECIES];
 Double_t prob1[AliPID::kSPECIES];
 Double_t prob2[AliPID::kSPECIES];

 Bool_t isTOF0=fPidHF->CheckTOFPIDStatus((AliAODTrack*)obj->GetDaughter(0));
 Bool_t isTOF1=fPidHF->CheckTOFPIDStatus((AliAODTrack*)obj->GetDaughter(1));
 Bool_t isTOF2=fPidHF->CheckTOFPIDStatus((AliAODTrack*)obj->GetDaughter(2));

Bool_t isK1=kFALSE;
 if(isTOF1){ //kaon
  if(track1->P()<1.8) {
    fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC|AliPIDResponse::kDetTOF);
    if(obj->Pt()<3. && track1->P()<0.55) fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC);
    fPidHF->GetPidCombined()->ComputeProbabilities(track1,fPidHF->GetPidResponse(),prob1);
   
  }else{
    AliAODTrack *trackaod1=(AliAODTrack*)(obj->GetDaughter(1));
    if(trackaod1->P()<0.55){
      fPidHF->SetTOF(kFALSE);
      fPidHF->SetTOFdecide(kFALSE);
     }
    Int_t isKaon=fPidHF->MakeRawPid(trackaod1,3);
    if(isKaon>=1) isK1=kTRUE;
    if(trackaod1->P()<0.55){
      fPidHF->SetTOF(kTRUE);
      fPidHF->SetTOFdecide(kTRUE);
     }
  }
 }else{
  if(track1->P()<0.8){
    fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC);
    fPidHF->GetPidCombined()->ComputeProbabilities(track1,fPidHF->GetPidResponse(),prob0);
  }else{
    AliAODTrack *trackaod1=(AliAODTrack*)(obj->GetDaughter(1));
     if(trackaod1->P()<0.55){
      fPidHF->SetTOF(kFALSE);
      fPidHF->SetTOFdecide(kFALSE);
     }
    Int_t isKaon=fPidHF->MakeRawPid(trackaod1,3);
    if(isKaon>=1) isK1=kTRUE;
     if(trackaod1->P()<0.55){
      fPidHF->SetTOF(kTRUE);
      fPidHF->SetTOFdecide(kTRUE);
     }
  }
 }

 Bool_t ispi0=kFALSE;
 Bool_t isp0=kFALSE;
 Bool_t ispi2=kFALSE;
 Bool_t isp2=kFALSE;

 if(isTOF0){ //proton
  if(track0->P()<2.2) {
    fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC|AliPIDResponse::kDetTOF);
    if(obj->Pt()<3. && track0->P()<1.) fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC);
    fPidHF->GetPidCombined()->ComputeProbabilities(track0,fPidHF->GetPidResponse(),prob0);
  }else{
   AliAODTrack *trackaod0=(AliAODTrack*)(obj->GetDaughter(0));
   if(trackaod0->P()<1.){
      fPidObjprot->SetTOF(kFALSE);
      fPidObjprot->SetTOFdecide(kFALSE);
   }
   Int_t isProton=fPidObjprot->MakeRawPid(trackaod0,4);
   if(isProton>=1) isp0=kTRUE;
   if(trackaod0->P()<1.){
      fPidObjprot->SetTOF(kTRUE);
      fPidObjprot->SetTOFdecide(kTRUE);
   }
  }
 }else{
   if(track0->P()<1.2){
    fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC);
    fPidHF->GetPidCombined()->ComputeProbabilities(track0,fPidHF->GetPidResponse(),prob0);
   }else{
    AliAODTrack *trackaod0=(AliAODTrack*)(obj->GetDaughter(0));
    if(trackaod0->P()<1.){
      fPidObjprot->SetTOF(kFALSE);
      fPidObjprot->SetTOFdecide(kFALSE);
    }
    Int_t isProton=fPidObjprot->MakeRawPid(trackaod0,4);
    if(isProton>=1) isp0=kTRUE;
    if(trackaod0->P()<1.){
      fPidObjprot->SetTOF(kTRUE);
      fPidObjprot->SetTOFdecide(kTRUE);
    }
   }
 }

 if(isTOF2){ //proton
  if(track2->P()<2.2) {
    fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC|AliPIDResponse::kDetTOF);
    if(obj->Pt()<3. && track2->P()<1.) fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC);
    fPidHF->GetPidCombined()->ComputeProbabilities(track2,fPidHF->GetPidResponse(),prob2);
  }else{
    AliAODTrack *trackaod2=(AliAODTrack*)(obj->GetDaughter(2));
    if(trackaod2->P()<1.){
      fPidObjprot->SetTOF(kFALSE);
      fPidObjprot->SetTOFdecide(kFALSE);
    }
   Int_t isProton=fPidObjprot->MakeRawPid(trackaod2,4);
   if(isProton>=1) isp2=kTRUE;
   if(trackaod2->P()<1.){
      fPidObjprot->SetTOF(kTRUE);
      fPidObjprot->SetTOFdecide(kTRUE);
    }
  }
 }else{
   if(track2->P()<1.2){
    fPidHF->GetPidCombined()->SetDetectorMask(AliPIDResponse::kDetTPC);
    fPidHF->GetPidCombined()->ComputeProbabilities(track2,fPidHF->GetPidResponse(),prob2);
   }else{
    AliAODTrack *trackaod2=(AliAODTrack*)(obj->GetDaughter(2));
     if(trackaod2->P()<1.){
      fPidObjprot->SetTOF(kFALSE);
      fPidObjprot->SetTOFdecide(kFALSE);
    }
    Int_t isProton=fPidObjprot->MakeRawPid(trackaod2,4);
    if(isProton>=1) isp2=kTRUE;
    if(trackaod2->P()<1.){
      fPidObjprot->SetTOF(kTRUE);
      fPidObjprot->SetTOFdecide(kTRUE);
    }
   }
 }
  AliAODTrack *trackaod2=(AliAODTrack*)(obj->GetDaughter(2));
  if(fPidObjpion->MakeRawPid(trackaod2,2)>=1)ispi2=kTRUE;
  AliAODTrack *trackaod0=(AliAODTrack*)(obj->GetDaughter(2));
  if(fPidObjpion->MakeRawPid(trackaod0,2)>=1)ispi0=kTRUE;

  if(TMath::MaxElement(AliPID::kSPECIES,prob1) == prob1[AliPID::kKaon]){

    if(TMath::MaxElement(AliPID::kSPECIES,prob0) == prob0[AliPID::kProton] && TMath::MaxElement(AliPID::kSPECIES,prob2) == prob2[AliPID::kPion]) okLcpKpi = 1;
    if(TMath::MaxElement(AliPID::kSPECIES,prob2) == prob2[AliPID::kProton] && TMath::MaxElement(AliPID::kSPECIES,prob0) == prob0[AliPID::kPion]) okLcpiKp = 1;
    }

   if(!isK1 && TMath::MaxElement(AliPID::kSPECIES,prob1) == prob1[AliPID::kKaon]) isK1=kTRUE;
    if(!ispi0 && TMath::MaxElement(AliPID::kSPECIES,prob0) == prob0[AliPID::kPion]) ispi0=kTRUE;
    if(!ispi2 && TMath::MaxElement(AliPID::kSPECIES,prob2) == prob2[AliPID::kPion]) ispi2=kTRUE;
    if(!isp0 && TMath::MaxElement(AliPID::kSPECIES,prob0) == prob0[AliPID::kProton]) isp0=kTRUE;
    if(!isp2 && TMath::MaxElement(AliPID::kSPECIES,prob2) == prob2[AliPID::kProton]) isp2=kTRUE;
    if(isK1 && ispi0 && isp2) okLcpiKp = 1;
    if(isK1 && isp0 && ispi2) okLcpKpi = 1;

   if(okLcpKpi) returnvalue=1; //cuts passed as Lc->pKpi
    if(okLcpiKp) returnvalue=2; //cuts passed as Lc->piKp
    if(okLcpKpi && okLcpiKp) returnvalue=3; //cuts passed as both pKpi and piKp

    return returnvalue;
 

}
//----------------------------------------------------------
Int_t AliRDHFCutsLctopKpi::IsSelectedPIDStrong(AliAODRecoDecayHF* obj) {


    if(!fUsePID || !obj) return 3;
    Int_t okLcpKpi=0,okLcpiKp=0;
    Int_t returnvalue=0;
    Bool_t isPeriodd=fPidHF->GetOnePad();
    Bool_t isMC=fPidHF->GetMC();
    Bool_t ispion0=kTRUE,ispion2=kTRUE;
    Bool_t isproton0=kFALSE,isproton2=kFALSE;
    Bool_t iskaon1=kFALSE;
    if(isPeriodd) {
     fPidObjprot->SetOnePad(kTRUE);
     fPidObjpion->SetOnePad(kTRUE);
    }
    if(isMC) {
     fPidObjprot->SetMC(kTRUE);
     fPidObjpion->SetMC(kTRUE);
    }

    for(Int_t i=0;i<3;i++){
     AliAODTrack *track=(AliAODTrack*)obj->GetDaughter(i);
     if(!track) return 0;
     // identify kaon
     if(i==1) {
      Int_t isKaon=fPidHF->MakeRawPid(track,3);
      if(isKaon>=1) {
       iskaon1=kTRUE;
      if(fPidHF->MakeRawPid(track,2)>=1) iskaon1=kFALSE;
      }
      if(!iskaon1) return 0;
     
     }else{
     //pion or proton
     
     Int_t isProton=fPidObjprot->MakeRawPid(track,4);
     if(isProton>=1){
      if(fPidHF->MakeRawPid(track,2)>=1) isProton=-1;
      if(fPidHF->MakeRawPid(track,3)>=1) isProton=-1;
     }

     Int_t isPion=fPidObjpion->MakeRawPid(track,2);
     if(fPidHF->MakeRawPid(track,3)>=1) isPion=-1;
     if(fPidObjprot->MakeRawPid(track,4)>=1) isPion=-1;


     if(i==0) {
      if(isPion<0) ispion0=kFALSE;
      if(isProton>=1) isproton0=kTRUE;

     }
      if(!ispion0 && !isproton0) return 0;
     if(i==2) {
      if(isPion<0) ispion2=kFALSE;
      if(isProton>=1) isproton2=kTRUE;
     }

    }
   }

    if(ispion2 && isproton0 && iskaon1) okLcpKpi=1;
    if(ispion0 && isproton2 && iskaon1) okLcpiKp=1;
    if(okLcpKpi) returnvalue=1; //cuts passed as Lc->pKpi
    if(okLcpiKp) returnvalue=2; //cuts passed as Lc->piKp
    if(okLcpKpi && okLcpiKp) returnvalue=3; //cuts passed as both pKpi and piKp

 return returnvalue;
}
//--------------------
Int_t AliRDHFCutsLctopKpi::IsSelectedCombinedPIDpPb(AliAODRecoDecayHF* obj) {
    
    if(!fUsePID || !obj) {return 3;}
    Int_t okLcpKpi=0,okLcpiKp=0;
    Int_t returnvalue=0;
   
    Bool_t isMC=fPidHF->GetMC();

    
    if(isMC) {
	    fPidObjprot->SetMC(kTRUE);
	    fPidObjpion->SetMC(kTRUE);
    }

    AliVTrack *track0=dynamic_cast<AliVTrack*>(obj->GetDaughter(0));
    AliVTrack *track1=dynamic_cast<AliVTrack*>(obj->GetDaughter(1));
    AliVTrack *track2=dynamic_cast<AliVTrack*>(obj->GetDaughter(2));
    if (!track0 || !track1 || !track2) return 0;
    Double_t prob0[AliPID::kSPECIES];
    Double_t prob1[AliPID::kSPECIES];
    Double_t prob2[AliPID::kSPECIES];
    
    fPidHF->GetPidCombined()->ComputeProbabilities(track0,fPidHF->GetPidResponse(),prob0);
    fPidHF->GetPidCombined()->ComputeProbabilities(track1,fPidHF->GetPidResponse(),prob1);
    fPidHF->GetPidCombined()->ComputeProbabilities(track2,fPidHF->GetPidResponse(),prob2);
  

    if(fPIDThreshold[AliPID::kPion]>0. && fPIDThreshold[AliPID::kKaon]>0. && fPIDThreshold[AliPID::kProton]>0.){
    okLcpiKp=  (prob0[AliPID::kPion  ]>fPIDThreshold[AliPID::kPion  ])
             &&(prob1[AliPID::kKaon  ]>fPIDThreshold[AliPID::kKaon  ])
             &&(prob2[AliPID::kProton]>fPIDThreshold[AliPID::kProton]);
    okLcpKpi=  (prob0[AliPID::kProton]>fPIDThreshold[AliPID::kProton])
             &&(prob1[AliPID::kKaon  ]>fPIDThreshold[AliPID::kKaon  ])
             &&(prob2[AliPID::kPion  ]>fPIDThreshold[AliPID::kPion  ]);
   }else{ 
		    //pion or proton
		    
		    
    if(TMath::MaxElement(AliPID::kSPECIES,prob1) == prob1[AliPID::kKaon]){
    if(TMath::MaxElement(AliPID::kSPECIES,prob0) == prob0[AliPID::kProton] && TMath::MaxElement(AliPID::kSPECIES,prob2) == prob2[AliPID::kPion]) okLcpKpi = 1;  
    if(TMath::MaxElement(AliPID::kSPECIES,prob2) == prob2[AliPID::kProton] && TMath::MaxElement(AliPID::kSPECIES,prob0) == prob0[AliPID::kPion]) okLcpiKp = 1; 
	    }
	   }
    
    if(okLcpKpi) returnvalue=1; //cuts passed as Lc->pKpi
    if(okLcpiKp) returnvalue=2; //cuts passed as Lc->piKp
    if(okLcpKpi && okLcpiKp) returnvalue=3; //cuts passed as both pKpi and piKp
    
    return returnvalue;
}
//-----------------------
Int_t AliRDHFCutsLctopKpi::IsSelectedCombinedPIDpPb2(AliAODRecoDecayHF* obj) {

  Int_t returnvalue =0;
  Double_t thresholdK =0.7;
  Double_t thresholdPi =0.3; //!
  Double_t thresholdPr =0.7;

  AliVTrack *track0=dynamic_cast<AliVTrack*>(obj->GetDaughter(0));
  AliVTrack *track1=dynamic_cast<AliVTrack*>(obj->GetDaughter(1));
  AliVTrack *track2=dynamic_cast<AliVTrack*>(obj->GetDaughter(2));

  if (!track0 || !track1 || !track2) return 0;
  Double_t prob0[AliPID::kSPECIES];
  Double_t prob1[AliPID::kSPECIES];
  Double_t prob2[AliPID::kSPECIES];

  fPidHF->GetPidCombined()->ComputeProbabilities(track0,fPidHF->GetPidResponse(),prob0);
  fPidHF->GetPidCombined()->ComputeProbabilities(track1,fPidHF->GetPidResponse(),prob1);
  fPidHF->GetPidCombined()->ComputeProbabilities(track2,fPidHF->GetPidResponse(),prob2);

  if(prob1[AliPID::kKaon]>thresholdK){
    if(TMath::MaxElement(AliPID::kSPECIES,prob0)>TMath::MaxElement(AliPID::kSPECIES,prob2)){
      if(((prob0[AliPID::kPion  ]>prob0[AliPID::kProton  ])&& prob0[AliPID::kPion]>thresholdPi) && prob2[AliPID::kProton]>thresholdPr) returnvalue=2;//piKp	
      else if(((prob0[AliPID::kProton  ]>prob0[AliPID::kPion  ])&& prob0[AliPID::kProton]>thresholdPr) && prob2[AliPID::kPion]>thresholdPi)returnvalue =1;//pKpi
    }

    else if(TMath::MaxElement(AliPID::kSPECIES,prob0)<TMath::MaxElement(AliPID::kSPECIES,prob2)){
      if(((prob2[AliPID::kPion  ]>prob2[AliPID::kProton  ])&& prob2[AliPID::kPion]>thresholdPi) && prob0[AliPID::kProton]>thresholdPr) returnvalue=1; //pKpi	
      else if(((prob2[AliPID::kProton  ]>prob2[AliPID::kPion  ])&& prob2[AliPID::kProton]>thresholdPr) && prob0[AliPID::kPion]>thresholdPi)returnvalue =2;	//piKp
    }

  }
  return returnvalue;

}
//------------------------------------------------------
void AliRDHFCutsLctopKpi::SetStandardCutsPPb2013() {

 SetName("LctopKpiProdCuts");
 SetTitle("Production cuts for Lc analysis");

 AliESDtrackCuts *esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts","default");
 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
 esdTrackCuts->SetRequireTPCRefit(kTRUE);
 esdTrackCuts->SetMinNClustersTPC(70);
 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,
                                          AliESDtrackCuts::kAny);
 esdTrackCuts->SetRequireITSRefit(kTRUE);
 esdTrackCuts->SetMinNClustersITS(4);
 esdTrackCuts->SetMinDCAToVertexXY(0.);
 esdTrackCuts->SetEtaRange(-0.8,0.8);
 esdTrackCuts->SetPtRange(0.3,1.e10);
 AddTrackCuts(esdTrackCuts);
 delete esdTrackCuts;
 esdTrackCuts=NULL;
 
 const Int_t nvars=13;
 const Int_t nptbins=9;
 Float_t* ptbins;
 ptbins=new Float_t[nptbins+1];
 ptbins[0]=0.;
 ptbins[1]=1.;
 ptbins[2]=2.;
 ptbins[3]=3.;
 ptbins[4]=4.;
 ptbins[5]=5.;
 ptbins[6]=6.;
 ptbins[7]=8.;
 ptbins[8]=10.;  
 ptbins[9]=99999.;

 SetGlobalIndex(nvars,nptbins);
 SetPtBins(nptbins+1,ptbins);

 Float_t** prodcutsval;
 prodcutsval=new Float_t*[nvars];
 for(Int_t iv=0;iv<nvars;iv++){
  prodcutsval[iv]=new Float_t[nptbins];
 }

 for(Int_t ipt=0;ipt<nptbins;ipt++){
    prodcutsval[0][ipt]=0.13;
    prodcutsval[1][ipt]=0.4;
    prodcutsval[2][ipt]=0.4;
    prodcutsval[3][ipt]=0.;
    prodcutsval[4][ipt]=0.;
    prodcutsval[5][ipt]=0.;
    prodcutsval[6][ipt]=0.06;
    prodcutsval[7][ipt]=0.;
    prodcutsval[8][ipt]=0.;
    prodcutsval[9][ipt]=0.;
    prodcutsval[10][ipt]=0.;
    prodcutsval[11][ipt]=0.05;
    prodcutsval[12][ipt]=0.4;
 }
 SetCuts(nvars,nptbins,prodcutsval);

 AliAODPidHF* pidObjK=new AliAODPidHF();
 Double_t sigmasK[5]={3.,1.,1.,3.,2.};
 pidObjK->SetSigma(sigmasK);
 pidObjK->SetAsym(kTRUE);
 pidObjK->SetMatch(1);
 pidObjK->SetTPC(kTRUE);
 pidObjK->SetTOF(kTRUE);
 pidObjK->SetITS(kTRUE);
 Double_t plimK[2]={0.5,0.8};
 pidObjK->SetPLimit(plimK,2);
 pidObjK->SetTOFdecide(kTRUE);
 SetPidHF(pidObjK);

 AliAODPidHF* pidObjpi=new AliAODPidHF();
 pidObjpi->SetTPC(kTRUE);
 Double_t sigmaspi[5]={3.,0.,0.,0.,0.};
 pidObjpi->SetSigma(sigmaspi);
 pidObjpi->SetTOFdecide(kTRUE);
 SetPidpion(pidObjpi);

 AliAODPidHF* pidObjp=new AliAODPidHF();
 Double_t sigmasp[5]={3.,1.,1.,3.,2.};
 pidObjp->SetSigma(sigmasp);
 pidObjp->SetAsym(kTRUE);
 pidObjp->SetMatch(1);
 pidObjp->SetTPC(kTRUE);
 pidObjp->SetTOF(kTRUE);
 pidObjp->SetITS(kTRUE);
 Double_t plimp[2]={1.,2.};
 pidObjp->SetPLimit(plimp,2);
 pidObjp->SetTOFdecide(kTRUE);

 SetUsePID(kTRUE);

 PrintAll();

 for(Int_t iiv=0;iiv<nvars;iiv++){
  delete [] prodcutsval[iiv];
 }
 delete [] prodcutsval;
 prodcutsval=NULL;
 delete [] ptbins;
 ptbins=NULL;

 delete pidObjK;
 pidObjK=NULL;
 delete pidObjpi;
 pidObjpi=NULL;
 delete pidObjp;
 pidObjp=NULL;

 return;
}

 AliRDHFCutsLctopKpi.cxx:1
 AliRDHFCutsLctopKpi.cxx:2
 AliRDHFCutsLctopKpi.cxx:3
 AliRDHFCutsLctopKpi.cxx:4
 AliRDHFCutsLctopKpi.cxx:5
 AliRDHFCutsLctopKpi.cxx:6
 AliRDHFCutsLctopKpi.cxx:7
 AliRDHFCutsLctopKpi.cxx:8
 AliRDHFCutsLctopKpi.cxx:9
 AliRDHFCutsLctopKpi.cxx:10
 AliRDHFCutsLctopKpi.cxx:11
 AliRDHFCutsLctopKpi.cxx:12
 AliRDHFCutsLctopKpi.cxx:13
 AliRDHFCutsLctopKpi.cxx:14
 AliRDHFCutsLctopKpi.cxx:15
 AliRDHFCutsLctopKpi.cxx:16
 AliRDHFCutsLctopKpi.cxx:17
 AliRDHFCutsLctopKpi.cxx:18
 AliRDHFCutsLctopKpi.cxx:19
 AliRDHFCutsLctopKpi.cxx:20
 AliRDHFCutsLctopKpi.cxx:21
 AliRDHFCutsLctopKpi.cxx:22
 AliRDHFCutsLctopKpi.cxx:23
 AliRDHFCutsLctopKpi.cxx:24
 AliRDHFCutsLctopKpi.cxx:25
 AliRDHFCutsLctopKpi.cxx:26
 AliRDHFCutsLctopKpi.cxx:27
 AliRDHFCutsLctopKpi.cxx:28
 AliRDHFCutsLctopKpi.cxx:29
 AliRDHFCutsLctopKpi.cxx:30
 AliRDHFCutsLctopKpi.cxx:31
 AliRDHFCutsLctopKpi.cxx:32
 AliRDHFCutsLctopKpi.cxx:33
 AliRDHFCutsLctopKpi.cxx:34
 AliRDHFCutsLctopKpi.cxx:35
 AliRDHFCutsLctopKpi.cxx:36
 AliRDHFCutsLctopKpi.cxx:37
 AliRDHFCutsLctopKpi.cxx:38
 AliRDHFCutsLctopKpi.cxx:39
 AliRDHFCutsLctopKpi.cxx:40
 AliRDHFCutsLctopKpi.cxx:41
 AliRDHFCutsLctopKpi.cxx:42
 AliRDHFCutsLctopKpi.cxx:43
 AliRDHFCutsLctopKpi.cxx:44
 AliRDHFCutsLctopKpi.cxx:45
 AliRDHFCutsLctopKpi.cxx:46
 AliRDHFCutsLctopKpi.cxx:47
 AliRDHFCutsLctopKpi.cxx:48
 AliRDHFCutsLctopKpi.cxx:49
 AliRDHFCutsLctopKpi.cxx:50
 AliRDHFCutsLctopKpi.cxx:51
 AliRDHFCutsLctopKpi.cxx:52
 AliRDHFCutsLctopKpi.cxx:53
 AliRDHFCutsLctopKpi.cxx:54
 AliRDHFCutsLctopKpi.cxx:55
 AliRDHFCutsLctopKpi.cxx:56
 AliRDHFCutsLctopKpi.cxx:57
 AliRDHFCutsLctopKpi.cxx:58
 AliRDHFCutsLctopKpi.cxx:59
 AliRDHFCutsLctopKpi.cxx:60
 AliRDHFCutsLctopKpi.cxx:61
 AliRDHFCutsLctopKpi.cxx:62
 AliRDHFCutsLctopKpi.cxx:63
 AliRDHFCutsLctopKpi.cxx:64
 AliRDHFCutsLctopKpi.cxx:65
 AliRDHFCutsLctopKpi.cxx:66
 AliRDHFCutsLctopKpi.cxx:67
 AliRDHFCutsLctopKpi.cxx:68
 AliRDHFCutsLctopKpi.cxx:69
 AliRDHFCutsLctopKpi.cxx:70
 AliRDHFCutsLctopKpi.cxx:71
 AliRDHFCutsLctopKpi.cxx:72
 AliRDHFCutsLctopKpi.cxx:73
 AliRDHFCutsLctopKpi.cxx:74
 AliRDHFCutsLctopKpi.cxx:75
 AliRDHFCutsLctopKpi.cxx:76
 AliRDHFCutsLctopKpi.cxx:77
 AliRDHFCutsLctopKpi.cxx:78
 AliRDHFCutsLctopKpi.cxx:79
 AliRDHFCutsLctopKpi.cxx:80
 AliRDHFCutsLctopKpi.cxx:81
 AliRDHFCutsLctopKpi.cxx:82
 AliRDHFCutsLctopKpi.cxx:83
 AliRDHFCutsLctopKpi.cxx:84
 AliRDHFCutsLctopKpi.cxx:85
 AliRDHFCutsLctopKpi.cxx:86
 AliRDHFCutsLctopKpi.cxx:87
 AliRDHFCutsLctopKpi.cxx:88
 AliRDHFCutsLctopKpi.cxx:89
 AliRDHFCutsLctopKpi.cxx:90
 AliRDHFCutsLctopKpi.cxx:91
 AliRDHFCutsLctopKpi.cxx:92
 AliRDHFCutsLctopKpi.cxx:93
 AliRDHFCutsLctopKpi.cxx:94
 AliRDHFCutsLctopKpi.cxx:95
 AliRDHFCutsLctopKpi.cxx:96
 AliRDHFCutsLctopKpi.cxx:97
 AliRDHFCutsLctopKpi.cxx:98
 AliRDHFCutsLctopKpi.cxx:99
 AliRDHFCutsLctopKpi.cxx:100
 AliRDHFCutsLctopKpi.cxx:101
 AliRDHFCutsLctopKpi.cxx:102
 AliRDHFCutsLctopKpi.cxx:103
 AliRDHFCutsLctopKpi.cxx:104
 AliRDHFCutsLctopKpi.cxx:105
 AliRDHFCutsLctopKpi.cxx:106
 AliRDHFCutsLctopKpi.cxx:107
 AliRDHFCutsLctopKpi.cxx:108
 AliRDHFCutsLctopKpi.cxx:109
 AliRDHFCutsLctopKpi.cxx:110
 AliRDHFCutsLctopKpi.cxx:111
 AliRDHFCutsLctopKpi.cxx:112
 AliRDHFCutsLctopKpi.cxx:113
 AliRDHFCutsLctopKpi.cxx:114
 AliRDHFCutsLctopKpi.cxx:115
 AliRDHFCutsLctopKpi.cxx:116
 AliRDHFCutsLctopKpi.cxx:117
 AliRDHFCutsLctopKpi.cxx:118
 AliRDHFCutsLctopKpi.cxx:119
 AliRDHFCutsLctopKpi.cxx:120
 AliRDHFCutsLctopKpi.cxx:121
 AliRDHFCutsLctopKpi.cxx:122
 AliRDHFCutsLctopKpi.cxx:123
 AliRDHFCutsLctopKpi.cxx:124
 AliRDHFCutsLctopKpi.cxx:125
 AliRDHFCutsLctopKpi.cxx:126
 AliRDHFCutsLctopKpi.cxx:127
 AliRDHFCutsLctopKpi.cxx:128
 AliRDHFCutsLctopKpi.cxx:129
 AliRDHFCutsLctopKpi.cxx:130
 AliRDHFCutsLctopKpi.cxx:131
 AliRDHFCutsLctopKpi.cxx:132
 AliRDHFCutsLctopKpi.cxx:133
 AliRDHFCutsLctopKpi.cxx:134
 AliRDHFCutsLctopKpi.cxx:135
 AliRDHFCutsLctopKpi.cxx:136
 AliRDHFCutsLctopKpi.cxx:137
 AliRDHFCutsLctopKpi.cxx:138
 AliRDHFCutsLctopKpi.cxx:139
 AliRDHFCutsLctopKpi.cxx:140
 AliRDHFCutsLctopKpi.cxx:141
 AliRDHFCutsLctopKpi.cxx:142
 AliRDHFCutsLctopKpi.cxx:143
 AliRDHFCutsLctopKpi.cxx:144
 AliRDHFCutsLctopKpi.cxx:145
 AliRDHFCutsLctopKpi.cxx:146
 AliRDHFCutsLctopKpi.cxx:147
 AliRDHFCutsLctopKpi.cxx:148
 AliRDHFCutsLctopKpi.cxx:149
 AliRDHFCutsLctopKpi.cxx:150
 AliRDHFCutsLctopKpi.cxx:151
 AliRDHFCutsLctopKpi.cxx:152
 AliRDHFCutsLctopKpi.cxx:153
 AliRDHFCutsLctopKpi.cxx:154
 AliRDHFCutsLctopKpi.cxx:155
 AliRDHFCutsLctopKpi.cxx:156
 AliRDHFCutsLctopKpi.cxx:157
 AliRDHFCutsLctopKpi.cxx:158
 AliRDHFCutsLctopKpi.cxx:159
 AliRDHFCutsLctopKpi.cxx:160
 AliRDHFCutsLctopKpi.cxx:161
 AliRDHFCutsLctopKpi.cxx:162
 AliRDHFCutsLctopKpi.cxx:163
 AliRDHFCutsLctopKpi.cxx:164
 AliRDHFCutsLctopKpi.cxx:165
 AliRDHFCutsLctopKpi.cxx:166
 AliRDHFCutsLctopKpi.cxx:167
 AliRDHFCutsLctopKpi.cxx:168
 AliRDHFCutsLctopKpi.cxx:169
 AliRDHFCutsLctopKpi.cxx:170
 AliRDHFCutsLctopKpi.cxx:171
 AliRDHFCutsLctopKpi.cxx:172
 AliRDHFCutsLctopKpi.cxx:173
 AliRDHFCutsLctopKpi.cxx:174
 AliRDHFCutsLctopKpi.cxx:175
 AliRDHFCutsLctopKpi.cxx:176
 AliRDHFCutsLctopKpi.cxx:177
 AliRDHFCutsLctopKpi.cxx:178
 AliRDHFCutsLctopKpi.cxx:179
 AliRDHFCutsLctopKpi.cxx:180
 AliRDHFCutsLctopKpi.cxx:181
 AliRDHFCutsLctopKpi.cxx:182
 AliRDHFCutsLctopKpi.cxx:183
 AliRDHFCutsLctopKpi.cxx:184
 AliRDHFCutsLctopKpi.cxx:185
 AliRDHFCutsLctopKpi.cxx:186
 AliRDHFCutsLctopKpi.cxx:187
 AliRDHFCutsLctopKpi.cxx:188
 AliRDHFCutsLctopKpi.cxx:189
 AliRDHFCutsLctopKpi.cxx:190
 AliRDHFCutsLctopKpi.cxx:191
 AliRDHFCutsLctopKpi.cxx:192
 AliRDHFCutsLctopKpi.cxx:193
 AliRDHFCutsLctopKpi.cxx:194
 AliRDHFCutsLctopKpi.cxx:195
 AliRDHFCutsLctopKpi.cxx:196
 AliRDHFCutsLctopKpi.cxx:197
 AliRDHFCutsLctopKpi.cxx:198
 AliRDHFCutsLctopKpi.cxx:199
 AliRDHFCutsLctopKpi.cxx:200
 AliRDHFCutsLctopKpi.cxx:201
 AliRDHFCutsLctopKpi.cxx:202
 AliRDHFCutsLctopKpi.cxx:203
 AliRDHFCutsLctopKpi.cxx:204
 AliRDHFCutsLctopKpi.cxx:205
 AliRDHFCutsLctopKpi.cxx:206
 AliRDHFCutsLctopKpi.cxx:207
 AliRDHFCutsLctopKpi.cxx:208
 AliRDHFCutsLctopKpi.cxx:209
 AliRDHFCutsLctopKpi.cxx:210
 AliRDHFCutsLctopKpi.cxx:211
 AliRDHFCutsLctopKpi.cxx:212
 AliRDHFCutsLctopKpi.cxx:213
 AliRDHFCutsLctopKpi.cxx:214
 AliRDHFCutsLctopKpi.cxx:215
 AliRDHFCutsLctopKpi.cxx:216
 AliRDHFCutsLctopKpi.cxx:217
 AliRDHFCutsLctopKpi.cxx:218
 AliRDHFCutsLctopKpi.cxx:219
 AliRDHFCutsLctopKpi.cxx:220
 AliRDHFCutsLctopKpi.cxx:221
 AliRDHFCutsLctopKpi.cxx:222
 AliRDHFCutsLctopKpi.cxx:223
 AliRDHFCutsLctopKpi.cxx:224
 AliRDHFCutsLctopKpi.cxx:225
 AliRDHFCutsLctopKpi.cxx:226
 AliRDHFCutsLctopKpi.cxx:227
 AliRDHFCutsLctopKpi.cxx:228
 AliRDHFCutsLctopKpi.cxx:229
 AliRDHFCutsLctopKpi.cxx:230
 AliRDHFCutsLctopKpi.cxx:231
 AliRDHFCutsLctopKpi.cxx:232
 AliRDHFCutsLctopKpi.cxx:233
 AliRDHFCutsLctopKpi.cxx:234
 AliRDHFCutsLctopKpi.cxx:235
 AliRDHFCutsLctopKpi.cxx:236
 AliRDHFCutsLctopKpi.cxx:237
 AliRDHFCutsLctopKpi.cxx:238
 AliRDHFCutsLctopKpi.cxx:239
 AliRDHFCutsLctopKpi.cxx:240
 AliRDHFCutsLctopKpi.cxx:241
 AliRDHFCutsLctopKpi.cxx:242
 AliRDHFCutsLctopKpi.cxx:243
 AliRDHFCutsLctopKpi.cxx:244
 AliRDHFCutsLctopKpi.cxx:245
 AliRDHFCutsLctopKpi.cxx:246
 AliRDHFCutsLctopKpi.cxx:247
 AliRDHFCutsLctopKpi.cxx:248
 AliRDHFCutsLctopKpi.cxx:249
 AliRDHFCutsLctopKpi.cxx:250
 AliRDHFCutsLctopKpi.cxx:251
 AliRDHFCutsLctopKpi.cxx:252
 AliRDHFCutsLctopKpi.cxx:253
 AliRDHFCutsLctopKpi.cxx:254
 AliRDHFCutsLctopKpi.cxx:255
 AliRDHFCutsLctopKpi.cxx:256
 AliRDHFCutsLctopKpi.cxx:257
 AliRDHFCutsLctopKpi.cxx:258
 AliRDHFCutsLctopKpi.cxx:259
 AliRDHFCutsLctopKpi.cxx:260
 AliRDHFCutsLctopKpi.cxx:261
 AliRDHFCutsLctopKpi.cxx:262
 AliRDHFCutsLctopKpi.cxx:263
 AliRDHFCutsLctopKpi.cxx:264
 AliRDHFCutsLctopKpi.cxx:265
 AliRDHFCutsLctopKpi.cxx:266
 AliRDHFCutsLctopKpi.cxx:267
 AliRDHFCutsLctopKpi.cxx:268
 AliRDHFCutsLctopKpi.cxx:269
 AliRDHFCutsLctopKpi.cxx:270
 AliRDHFCutsLctopKpi.cxx:271
 AliRDHFCutsLctopKpi.cxx:272
 AliRDHFCutsLctopKpi.cxx:273
 AliRDHFCutsLctopKpi.cxx:274
 AliRDHFCutsLctopKpi.cxx:275
 AliRDHFCutsLctopKpi.cxx:276
 AliRDHFCutsLctopKpi.cxx:277
 AliRDHFCutsLctopKpi.cxx:278
 AliRDHFCutsLctopKpi.cxx:279
 AliRDHFCutsLctopKpi.cxx:280
 AliRDHFCutsLctopKpi.cxx:281
 AliRDHFCutsLctopKpi.cxx:282
 AliRDHFCutsLctopKpi.cxx:283
 AliRDHFCutsLctopKpi.cxx:284
 AliRDHFCutsLctopKpi.cxx:285
 AliRDHFCutsLctopKpi.cxx:286
 AliRDHFCutsLctopKpi.cxx:287
 AliRDHFCutsLctopKpi.cxx:288
 AliRDHFCutsLctopKpi.cxx:289
 AliRDHFCutsLctopKpi.cxx:290
 AliRDHFCutsLctopKpi.cxx:291
 AliRDHFCutsLctopKpi.cxx:292
 AliRDHFCutsLctopKpi.cxx:293
 AliRDHFCutsLctopKpi.cxx:294
 AliRDHFCutsLctopKpi.cxx:295
 AliRDHFCutsLctopKpi.cxx:296
 AliRDHFCutsLctopKpi.cxx:297
 AliRDHFCutsLctopKpi.cxx:298
 AliRDHFCutsLctopKpi.cxx:299
 AliRDHFCutsLctopKpi.cxx:300
 AliRDHFCutsLctopKpi.cxx:301
 AliRDHFCutsLctopKpi.cxx:302
 AliRDHFCutsLctopKpi.cxx:303
 AliRDHFCutsLctopKpi.cxx:304
 AliRDHFCutsLctopKpi.cxx:305
 AliRDHFCutsLctopKpi.cxx:306
 AliRDHFCutsLctopKpi.cxx:307
 AliRDHFCutsLctopKpi.cxx:308
 AliRDHFCutsLctopKpi.cxx:309
 AliRDHFCutsLctopKpi.cxx:310
 AliRDHFCutsLctopKpi.cxx:311
 AliRDHFCutsLctopKpi.cxx:312
 AliRDHFCutsLctopKpi.cxx:313
 AliRDHFCutsLctopKpi.cxx:314
 AliRDHFCutsLctopKpi.cxx:315
 AliRDHFCutsLctopKpi.cxx:316
 AliRDHFCutsLctopKpi.cxx:317
 AliRDHFCutsLctopKpi.cxx:318
 AliRDHFCutsLctopKpi.cxx:319
 AliRDHFCutsLctopKpi.cxx:320
 AliRDHFCutsLctopKpi.cxx:321
 AliRDHFCutsLctopKpi.cxx:322
 AliRDHFCutsLctopKpi.cxx:323
 AliRDHFCutsLctopKpi.cxx:324
 AliRDHFCutsLctopKpi.cxx:325
 AliRDHFCutsLctopKpi.cxx:326
 AliRDHFCutsLctopKpi.cxx:327
 AliRDHFCutsLctopKpi.cxx:328
 AliRDHFCutsLctopKpi.cxx:329
 AliRDHFCutsLctopKpi.cxx:330
 AliRDHFCutsLctopKpi.cxx:331
 AliRDHFCutsLctopKpi.cxx:332
 AliRDHFCutsLctopKpi.cxx:333
 AliRDHFCutsLctopKpi.cxx:334
 AliRDHFCutsLctopKpi.cxx:335
 AliRDHFCutsLctopKpi.cxx:336
 AliRDHFCutsLctopKpi.cxx:337
 AliRDHFCutsLctopKpi.cxx:338
 AliRDHFCutsLctopKpi.cxx:339
 AliRDHFCutsLctopKpi.cxx:340
 AliRDHFCutsLctopKpi.cxx:341
 AliRDHFCutsLctopKpi.cxx:342
 AliRDHFCutsLctopKpi.cxx:343
 AliRDHFCutsLctopKpi.cxx:344
 AliRDHFCutsLctopKpi.cxx:345
 AliRDHFCutsLctopKpi.cxx:346
 AliRDHFCutsLctopKpi.cxx:347
 AliRDHFCutsLctopKpi.cxx:348
 AliRDHFCutsLctopKpi.cxx:349
 AliRDHFCutsLctopKpi.cxx:350
 AliRDHFCutsLctopKpi.cxx:351
 AliRDHFCutsLctopKpi.cxx:352
 AliRDHFCutsLctopKpi.cxx:353
 AliRDHFCutsLctopKpi.cxx:354
 AliRDHFCutsLctopKpi.cxx:355
 AliRDHFCutsLctopKpi.cxx:356
 AliRDHFCutsLctopKpi.cxx:357
 AliRDHFCutsLctopKpi.cxx:358
 AliRDHFCutsLctopKpi.cxx:359
 AliRDHFCutsLctopKpi.cxx:360
 AliRDHFCutsLctopKpi.cxx:361
 AliRDHFCutsLctopKpi.cxx:362
 AliRDHFCutsLctopKpi.cxx:363
 AliRDHFCutsLctopKpi.cxx:364
 AliRDHFCutsLctopKpi.cxx:365
 AliRDHFCutsLctopKpi.cxx:366
 AliRDHFCutsLctopKpi.cxx:367
 AliRDHFCutsLctopKpi.cxx:368
 AliRDHFCutsLctopKpi.cxx:369
 AliRDHFCutsLctopKpi.cxx:370
 AliRDHFCutsLctopKpi.cxx:371
 AliRDHFCutsLctopKpi.cxx:372
 AliRDHFCutsLctopKpi.cxx:373
 AliRDHFCutsLctopKpi.cxx:374
 AliRDHFCutsLctopKpi.cxx:375
 AliRDHFCutsLctopKpi.cxx:376
 AliRDHFCutsLctopKpi.cxx:377
 AliRDHFCutsLctopKpi.cxx:378
 AliRDHFCutsLctopKpi.cxx:379
 AliRDHFCutsLctopKpi.cxx:380
 AliRDHFCutsLctopKpi.cxx:381
 AliRDHFCutsLctopKpi.cxx:382
 AliRDHFCutsLctopKpi.cxx:383
 AliRDHFCutsLctopKpi.cxx:384
 AliRDHFCutsLctopKpi.cxx:385
 AliRDHFCutsLctopKpi.cxx:386
 AliRDHFCutsLctopKpi.cxx:387
 AliRDHFCutsLctopKpi.cxx:388
 AliRDHFCutsLctopKpi.cxx:389
 AliRDHFCutsLctopKpi.cxx:390
 AliRDHFCutsLctopKpi.cxx:391
 AliRDHFCutsLctopKpi.cxx:392
 AliRDHFCutsLctopKpi.cxx:393
 AliRDHFCutsLctopKpi.cxx:394
 AliRDHFCutsLctopKpi.cxx:395
 AliRDHFCutsLctopKpi.cxx:396
 AliRDHFCutsLctopKpi.cxx:397
 AliRDHFCutsLctopKpi.cxx:398
 AliRDHFCutsLctopKpi.cxx:399
 AliRDHFCutsLctopKpi.cxx:400
 AliRDHFCutsLctopKpi.cxx:401
 AliRDHFCutsLctopKpi.cxx:402
 AliRDHFCutsLctopKpi.cxx:403
 AliRDHFCutsLctopKpi.cxx:404
 AliRDHFCutsLctopKpi.cxx:405
 AliRDHFCutsLctopKpi.cxx:406
 AliRDHFCutsLctopKpi.cxx:407
 AliRDHFCutsLctopKpi.cxx:408
 AliRDHFCutsLctopKpi.cxx:409
 AliRDHFCutsLctopKpi.cxx:410
 AliRDHFCutsLctopKpi.cxx:411
 AliRDHFCutsLctopKpi.cxx:412
 AliRDHFCutsLctopKpi.cxx:413
 AliRDHFCutsLctopKpi.cxx:414
 AliRDHFCutsLctopKpi.cxx:415
 AliRDHFCutsLctopKpi.cxx:416
 AliRDHFCutsLctopKpi.cxx:417
 AliRDHFCutsLctopKpi.cxx:418
 AliRDHFCutsLctopKpi.cxx:419
 AliRDHFCutsLctopKpi.cxx:420
 AliRDHFCutsLctopKpi.cxx:421
 AliRDHFCutsLctopKpi.cxx:422
 AliRDHFCutsLctopKpi.cxx:423
 AliRDHFCutsLctopKpi.cxx:424
 AliRDHFCutsLctopKpi.cxx:425
 AliRDHFCutsLctopKpi.cxx:426
 AliRDHFCutsLctopKpi.cxx:427
 AliRDHFCutsLctopKpi.cxx:428
 AliRDHFCutsLctopKpi.cxx:429
 AliRDHFCutsLctopKpi.cxx:430
 AliRDHFCutsLctopKpi.cxx:431
 AliRDHFCutsLctopKpi.cxx:432
 AliRDHFCutsLctopKpi.cxx:433
 AliRDHFCutsLctopKpi.cxx:434
 AliRDHFCutsLctopKpi.cxx:435
 AliRDHFCutsLctopKpi.cxx:436
 AliRDHFCutsLctopKpi.cxx:437
 AliRDHFCutsLctopKpi.cxx:438
 AliRDHFCutsLctopKpi.cxx:439
 AliRDHFCutsLctopKpi.cxx:440
 AliRDHFCutsLctopKpi.cxx:441
 AliRDHFCutsLctopKpi.cxx:442
 AliRDHFCutsLctopKpi.cxx:443
 AliRDHFCutsLctopKpi.cxx:444
 AliRDHFCutsLctopKpi.cxx:445
 AliRDHFCutsLctopKpi.cxx:446
 AliRDHFCutsLctopKpi.cxx:447
 AliRDHFCutsLctopKpi.cxx:448
 AliRDHFCutsLctopKpi.cxx:449
 AliRDHFCutsLctopKpi.cxx:450
 AliRDHFCutsLctopKpi.cxx:451
 AliRDHFCutsLctopKpi.cxx:452
 AliRDHFCutsLctopKpi.cxx:453
 AliRDHFCutsLctopKpi.cxx:454
 AliRDHFCutsLctopKpi.cxx:455
 AliRDHFCutsLctopKpi.cxx:456
 AliRDHFCutsLctopKpi.cxx:457
 AliRDHFCutsLctopKpi.cxx:458
 AliRDHFCutsLctopKpi.cxx:459
 AliRDHFCutsLctopKpi.cxx:460
 AliRDHFCutsLctopKpi.cxx:461
 AliRDHFCutsLctopKpi.cxx:462
 AliRDHFCutsLctopKpi.cxx:463
 AliRDHFCutsLctopKpi.cxx:464
 AliRDHFCutsLctopKpi.cxx:465
 AliRDHFCutsLctopKpi.cxx:466
 AliRDHFCutsLctopKpi.cxx:467
 AliRDHFCutsLctopKpi.cxx:468
 AliRDHFCutsLctopKpi.cxx:469
 AliRDHFCutsLctopKpi.cxx:470
 AliRDHFCutsLctopKpi.cxx:471
 AliRDHFCutsLctopKpi.cxx:472
 AliRDHFCutsLctopKpi.cxx:473
 AliRDHFCutsLctopKpi.cxx:474
 AliRDHFCutsLctopKpi.cxx:475
 AliRDHFCutsLctopKpi.cxx:476
 AliRDHFCutsLctopKpi.cxx:477
 AliRDHFCutsLctopKpi.cxx:478
 AliRDHFCutsLctopKpi.cxx:479
 AliRDHFCutsLctopKpi.cxx:480
 AliRDHFCutsLctopKpi.cxx:481
 AliRDHFCutsLctopKpi.cxx:482
 AliRDHFCutsLctopKpi.cxx:483
 AliRDHFCutsLctopKpi.cxx:484
 AliRDHFCutsLctopKpi.cxx:485
 AliRDHFCutsLctopKpi.cxx:486
 AliRDHFCutsLctopKpi.cxx:487
 AliRDHFCutsLctopKpi.cxx:488
 AliRDHFCutsLctopKpi.cxx:489
 AliRDHFCutsLctopKpi.cxx:490
 AliRDHFCutsLctopKpi.cxx:491
 AliRDHFCutsLctopKpi.cxx:492
 AliRDHFCutsLctopKpi.cxx:493
 AliRDHFCutsLctopKpi.cxx:494
 AliRDHFCutsLctopKpi.cxx:495
 AliRDHFCutsLctopKpi.cxx:496
 AliRDHFCutsLctopKpi.cxx:497
 AliRDHFCutsLctopKpi.cxx:498
 AliRDHFCutsLctopKpi.cxx:499
 AliRDHFCutsLctopKpi.cxx:500
 AliRDHFCutsLctopKpi.cxx:501
 AliRDHFCutsLctopKpi.cxx:502
 AliRDHFCutsLctopKpi.cxx:503
 AliRDHFCutsLctopKpi.cxx:504
 AliRDHFCutsLctopKpi.cxx:505
 AliRDHFCutsLctopKpi.cxx:506
 AliRDHFCutsLctopKpi.cxx:507
 AliRDHFCutsLctopKpi.cxx:508
 AliRDHFCutsLctopKpi.cxx:509
 AliRDHFCutsLctopKpi.cxx:510
 AliRDHFCutsLctopKpi.cxx:511
 AliRDHFCutsLctopKpi.cxx:512
 AliRDHFCutsLctopKpi.cxx:513
 AliRDHFCutsLctopKpi.cxx:514
 AliRDHFCutsLctopKpi.cxx:515
 AliRDHFCutsLctopKpi.cxx:516
 AliRDHFCutsLctopKpi.cxx:517
 AliRDHFCutsLctopKpi.cxx:518
 AliRDHFCutsLctopKpi.cxx:519
 AliRDHFCutsLctopKpi.cxx:520
 AliRDHFCutsLctopKpi.cxx:521
 AliRDHFCutsLctopKpi.cxx:522
 AliRDHFCutsLctopKpi.cxx:523
 AliRDHFCutsLctopKpi.cxx:524
 AliRDHFCutsLctopKpi.cxx:525
 AliRDHFCutsLctopKpi.cxx:526
 AliRDHFCutsLctopKpi.cxx:527
 AliRDHFCutsLctopKpi.cxx:528
 AliRDHFCutsLctopKpi.cxx:529
 AliRDHFCutsLctopKpi.cxx:530
 AliRDHFCutsLctopKpi.cxx:531
 AliRDHFCutsLctopKpi.cxx:532
 AliRDHFCutsLctopKpi.cxx:533
 AliRDHFCutsLctopKpi.cxx:534
 AliRDHFCutsLctopKpi.cxx:535
 AliRDHFCutsLctopKpi.cxx:536
 AliRDHFCutsLctopKpi.cxx:537
 AliRDHFCutsLctopKpi.cxx:538
 AliRDHFCutsLctopKpi.cxx:539
 AliRDHFCutsLctopKpi.cxx:540
 AliRDHFCutsLctopKpi.cxx:541
 AliRDHFCutsLctopKpi.cxx:542
 AliRDHFCutsLctopKpi.cxx:543
 AliRDHFCutsLctopKpi.cxx:544
 AliRDHFCutsLctopKpi.cxx:545
 AliRDHFCutsLctopKpi.cxx:546
 AliRDHFCutsLctopKpi.cxx:547
 AliRDHFCutsLctopKpi.cxx:548
 AliRDHFCutsLctopKpi.cxx:549
 AliRDHFCutsLctopKpi.cxx:550
 AliRDHFCutsLctopKpi.cxx:551
 AliRDHFCutsLctopKpi.cxx:552
 AliRDHFCutsLctopKpi.cxx:553
 AliRDHFCutsLctopKpi.cxx:554
 AliRDHFCutsLctopKpi.cxx:555
 AliRDHFCutsLctopKpi.cxx:556
 AliRDHFCutsLctopKpi.cxx:557
 AliRDHFCutsLctopKpi.cxx:558
 AliRDHFCutsLctopKpi.cxx:559
 AliRDHFCutsLctopKpi.cxx:560
 AliRDHFCutsLctopKpi.cxx:561
 AliRDHFCutsLctopKpi.cxx:562
 AliRDHFCutsLctopKpi.cxx:563
 AliRDHFCutsLctopKpi.cxx:564
 AliRDHFCutsLctopKpi.cxx:565
 AliRDHFCutsLctopKpi.cxx:566
 AliRDHFCutsLctopKpi.cxx:567
 AliRDHFCutsLctopKpi.cxx:568
 AliRDHFCutsLctopKpi.cxx:569
 AliRDHFCutsLctopKpi.cxx:570
 AliRDHFCutsLctopKpi.cxx:571
 AliRDHFCutsLctopKpi.cxx:572
 AliRDHFCutsLctopKpi.cxx:573
 AliRDHFCutsLctopKpi.cxx:574
 AliRDHFCutsLctopKpi.cxx:575
 AliRDHFCutsLctopKpi.cxx:576
 AliRDHFCutsLctopKpi.cxx:577
 AliRDHFCutsLctopKpi.cxx:578
 AliRDHFCutsLctopKpi.cxx:579
 AliRDHFCutsLctopKpi.cxx:580
 AliRDHFCutsLctopKpi.cxx:581
 AliRDHFCutsLctopKpi.cxx:582
 AliRDHFCutsLctopKpi.cxx:583
 AliRDHFCutsLctopKpi.cxx:584
 AliRDHFCutsLctopKpi.cxx:585
 AliRDHFCutsLctopKpi.cxx:586
 AliRDHFCutsLctopKpi.cxx:587
 AliRDHFCutsLctopKpi.cxx:588
 AliRDHFCutsLctopKpi.cxx:589
 AliRDHFCutsLctopKpi.cxx:590
 AliRDHFCutsLctopKpi.cxx:591
 AliRDHFCutsLctopKpi.cxx:592
 AliRDHFCutsLctopKpi.cxx:593
 AliRDHFCutsLctopKpi.cxx:594
 AliRDHFCutsLctopKpi.cxx:595
 AliRDHFCutsLctopKpi.cxx:596
 AliRDHFCutsLctopKpi.cxx:597
 AliRDHFCutsLctopKpi.cxx:598
 AliRDHFCutsLctopKpi.cxx:599
 AliRDHFCutsLctopKpi.cxx:600
 AliRDHFCutsLctopKpi.cxx:601
 AliRDHFCutsLctopKpi.cxx:602
 AliRDHFCutsLctopKpi.cxx:603
 AliRDHFCutsLctopKpi.cxx:604
 AliRDHFCutsLctopKpi.cxx:605
 AliRDHFCutsLctopKpi.cxx:606
 AliRDHFCutsLctopKpi.cxx:607
 AliRDHFCutsLctopKpi.cxx:608
 AliRDHFCutsLctopKpi.cxx:609
 AliRDHFCutsLctopKpi.cxx:610
 AliRDHFCutsLctopKpi.cxx:611
 AliRDHFCutsLctopKpi.cxx:612
 AliRDHFCutsLctopKpi.cxx:613
 AliRDHFCutsLctopKpi.cxx:614
 AliRDHFCutsLctopKpi.cxx:615
 AliRDHFCutsLctopKpi.cxx:616
 AliRDHFCutsLctopKpi.cxx:617
 AliRDHFCutsLctopKpi.cxx:618
 AliRDHFCutsLctopKpi.cxx:619
 AliRDHFCutsLctopKpi.cxx:620
 AliRDHFCutsLctopKpi.cxx:621
 AliRDHFCutsLctopKpi.cxx:622
 AliRDHFCutsLctopKpi.cxx:623
 AliRDHFCutsLctopKpi.cxx:624
 AliRDHFCutsLctopKpi.cxx:625
 AliRDHFCutsLctopKpi.cxx:626
 AliRDHFCutsLctopKpi.cxx:627
 AliRDHFCutsLctopKpi.cxx:628
 AliRDHFCutsLctopKpi.cxx:629
 AliRDHFCutsLctopKpi.cxx:630
 AliRDHFCutsLctopKpi.cxx:631
 AliRDHFCutsLctopKpi.cxx:632
 AliRDHFCutsLctopKpi.cxx:633
 AliRDHFCutsLctopKpi.cxx:634
 AliRDHFCutsLctopKpi.cxx:635
 AliRDHFCutsLctopKpi.cxx:636
 AliRDHFCutsLctopKpi.cxx:637
 AliRDHFCutsLctopKpi.cxx:638
 AliRDHFCutsLctopKpi.cxx:639
 AliRDHFCutsLctopKpi.cxx:640
 AliRDHFCutsLctopKpi.cxx:641
 AliRDHFCutsLctopKpi.cxx:642
 AliRDHFCutsLctopKpi.cxx:643
 AliRDHFCutsLctopKpi.cxx:644
 AliRDHFCutsLctopKpi.cxx:645
 AliRDHFCutsLctopKpi.cxx:646
 AliRDHFCutsLctopKpi.cxx:647
 AliRDHFCutsLctopKpi.cxx:648
 AliRDHFCutsLctopKpi.cxx:649
 AliRDHFCutsLctopKpi.cxx:650
 AliRDHFCutsLctopKpi.cxx:651
 AliRDHFCutsLctopKpi.cxx:652
 AliRDHFCutsLctopKpi.cxx:653
 AliRDHFCutsLctopKpi.cxx:654
 AliRDHFCutsLctopKpi.cxx:655
 AliRDHFCutsLctopKpi.cxx:656
 AliRDHFCutsLctopKpi.cxx:657
 AliRDHFCutsLctopKpi.cxx:658
 AliRDHFCutsLctopKpi.cxx:659
 AliRDHFCutsLctopKpi.cxx:660
 AliRDHFCutsLctopKpi.cxx:661
 AliRDHFCutsLctopKpi.cxx:662
 AliRDHFCutsLctopKpi.cxx:663
 AliRDHFCutsLctopKpi.cxx:664
 AliRDHFCutsLctopKpi.cxx:665
 AliRDHFCutsLctopKpi.cxx:666
 AliRDHFCutsLctopKpi.cxx:667
 AliRDHFCutsLctopKpi.cxx:668
 AliRDHFCutsLctopKpi.cxx:669
 AliRDHFCutsLctopKpi.cxx:670
 AliRDHFCutsLctopKpi.cxx:671
 AliRDHFCutsLctopKpi.cxx:672
 AliRDHFCutsLctopKpi.cxx:673
 AliRDHFCutsLctopKpi.cxx:674
 AliRDHFCutsLctopKpi.cxx:675
 AliRDHFCutsLctopKpi.cxx:676
 AliRDHFCutsLctopKpi.cxx:677
 AliRDHFCutsLctopKpi.cxx:678
 AliRDHFCutsLctopKpi.cxx:679
 AliRDHFCutsLctopKpi.cxx:680
 AliRDHFCutsLctopKpi.cxx:681
 AliRDHFCutsLctopKpi.cxx:682
 AliRDHFCutsLctopKpi.cxx:683
 AliRDHFCutsLctopKpi.cxx:684
 AliRDHFCutsLctopKpi.cxx:685
 AliRDHFCutsLctopKpi.cxx:686
 AliRDHFCutsLctopKpi.cxx:687
 AliRDHFCutsLctopKpi.cxx:688
 AliRDHFCutsLctopKpi.cxx:689
 AliRDHFCutsLctopKpi.cxx:690
 AliRDHFCutsLctopKpi.cxx:691
 AliRDHFCutsLctopKpi.cxx:692
 AliRDHFCutsLctopKpi.cxx:693
 AliRDHFCutsLctopKpi.cxx:694
 AliRDHFCutsLctopKpi.cxx:695
 AliRDHFCutsLctopKpi.cxx:696
 AliRDHFCutsLctopKpi.cxx:697
 AliRDHFCutsLctopKpi.cxx:698
 AliRDHFCutsLctopKpi.cxx:699
 AliRDHFCutsLctopKpi.cxx:700
 AliRDHFCutsLctopKpi.cxx:701
 AliRDHFCutsLctopKpi.cxx:702
 AliRDHFCutsLctopKpi.cxx:703
 AliRDHFCutsLctopKpi.cxx:704
 AliRDHFCutsLctopKpi.cxx:705
 AliRDHFCutsLctopKpi.cxx:706
 AliRDHFCutsLctopKpi.cxx:707
 AliRDHFCutsLctopKpi.cxx:708
 AliRDHFCutsLctopKpi.cxx:709
 AliRDHFCutsLctopKpi.cxx:710
 AliRDHFCutsLctopKpi.cxx:711
 AliRDHFCutsLctopKpi.cxx:712
 AliRDHFCutsLctopKpi.cxx:713
 AliRDHFCutsLctopKpi.cxx:714
 AliRDHFCutsLctopKpi.cxx:715
 AliRDHFCutsLctopKpi.cxx:716
 AliRDHFCutsLctopKpi.cxx:717
 AliRDHFCutsLctopKpi.cxx:718
 AliRDHFCutsLctopKpi.cxx:719
 AliRDHFCutsLctopKpi.cxx:720
 AliRDHFCutsLctopKpi.cxx:721
 AliRDHFCutsLctopKpi.cxx:722
 AliRDHFCutsLctopKpi.cxx:723
 AliRDHFCutsLctopKpi.cxx:724
 AliRDHFCutsLctopKpi.cxx:725
 AliRDHFCutsLctopKpi.cxx:726
 AliRDHFCutsLctopKpi.cxx:727
 AliRDHFCutsLctopKpi.cxx:728
 AliRDHFCutsLctopKpi.cxx:729
 AliRDHFCutsLctopKpi.cxx:730
 AliRDHFCutsLctopKpi.cxx:731
 AliRDHFCutsLctopKpi.cxx:732
 AliRDHFCutsLctopKpi.cxx:733
 AliRDHFCutsLctopKpi.cxx:734
 AliRDHFCutsLctopKpi.cxx:735
 AliRDHFCutsLctopKpi.cxx:736
 AliRDHFCutsLctopKpi.cxx:737
 AliRDHFCutsLctopKpi.cxx:738
 AliRDHFCutsLctopKpi.cxx:739
 AliRDHFCutsLctopKpi.cxx:740
 AliRDHFCutsLctopKpi.cxx:741
 AliRDHFCutsLctopKpi.cxx:742
 AliRDHFCutsLctopKpi.cxx:743
 AliRDHFCutsLctopKpi.cxx:744
 AliRDHFCutsLctopKpi.cxx:745
 AliRDHFCutsLctopKpi.cxx:746
 AliRDHFCutsLctopKpi.cxx:747
 AliRDHFCutsLctopKpi.cxx:748
 AliRDHFCutsLctopKpi.cxx:749
 AliRDHFCutsLctopKpi.cxx:750
 AliRDHFCutsLctopKpi.cxx:751
 AliRDHFCutsLctopKpi.cxx:752
 AliRDHFCutsLctopKpi.cxx:753
 AliRDHFCutsLctopKpi.cxx:754
 AliRDHFCutsLctopKpi.cxx:755
 AliRDHFCutsLctopKpi.cxx:756
 AliRDHFCutsLctopKpi.cxx:757
 AliRDHFCutsLctopKpi.cxx:758
 AliRDHFCutsLctopKpi.cxx:759
 AliRDHFCutsLctopKpi.cxx:760
 AliRDHFCutsLctopKpi.cxx:761
 AliRDHFCutsLctopKpi.cxx:762
 AliRDHFCutsLctopKpi.cxx:763
 AliRDHFCutsLctopKpi.cxx:764
 AliRDHFCutsLctopKpi.cxx:765
 AliRDHFCutsLctopKpi.cxx:766
 AliRDHFCutsLctopKpi.cxx:767
 AliRDHFCutsLctopKpi.cxx:768
 AliRDHFCutsLctopKpi.cxx:769
 AliRDHFCutsLctopKpi.cxx:770
 AliRDHFCutsLctopKpi.cxx:771
 AliRDHFCutsLctopKpi.cxx:772
 AliRDHFCutsLctopKpi.cxx:773
 AliRDHFCutsLctopKpi.cxx:774
 AliRDHFCutsLctopKpi.cxx:775
 AliRDHFCutsLctopKpi.cxx:776
 AliRDHFCutsLctopKpi.cxx:777
 AliRDHFCutsLctopKpi.cxx:778
 AliRDHFCutsLctopKpi.cxx:779
 AliRDHFCutsLctopKpi.cxx:780
 AliRDHFCutsLctopKpi.cxx:781
 AliRDHFCutsLctopKpi.cxx:782
 AliRDHFCutsLctopKpi.cxx:783
 AliRDHFCutsLctopKpi.cxx:784
 AliRDHFCutsLctopKpi.cxx:785
 AliRDHFCutsLctopKpi.cxx:786
 AliRDHFCutsLctopKpi.cxx:787
 AliRDHFCutsLctopKpi.cxx:788
 AliRDHFCutsLctopKpi.cxx:789
 AliRDHFCutsLctopKpi.cxx:790
 AliRDHFCutsLctopKpi.cxx:791
 AliRDHFCutsLctopKpi.cxx:792
 AliRDHFCutsLctopKpi.cxx:793
 AliRDHFCutsLctopKpi.cxx:794
 AliRDHFCutsLctopKpi.cxx:795
 AliRDHFCutsLctopKpi.cxx:796
 AliRDHFCutsLctopKpi.cxx:797
 AliRDHFCutsLctopKpi.cxx:798
 AliRDHFCutsLctopKpi.cxx:799
 AliRDHFCutsLctopKpi.cxx:800
 AliRDHFCutsLctopKpi.cxx:801
 AliRDHFCutsLctopKpi.cxx:802
 AliRDHFCutsLctopKpi.cxx:803
 AliRDHFCutsLctopKpi.cxx:804
 AliRDHFCutsLctopKpi.cxx:805
 AliRDHFCutsLctopKpi.cxx:806
 AliRDHFCutsLctopKpi.cxx:807
 AliRDHFCutsLctopKpi.cxx:808
 AliRDHFCutsLctopKpi.cxx:809
 AliRDHFCutsLctopKpi.cxx:810
 AliRDHFCutsLctopKpi.cxx:811
 AliRDHFCutsLctopKpi.cxx:812
 AliRDHFCutsLctopKpi.cxx:813
 AliRDHFCutsLctopKpi.cxx:814
 AliRDHFCutsLctopKpi.cxx:815
 AliRDHFCutsLctopKpi.cxx:816
 AliRDHFCutsLctopKpi.cxx:817
 AliRDHFCutsLctopKpi.cxx:818
 AliRDHFCutsLctopKpi.cxx:819
 AliRDHFCutsLctopKpi.cxx:820
 AliRDHFCutsLctopKpi.cxx:821
 AliRDHFCutsLctopKpi.cxx:822
 AliRDHFCutsLctopKpi.cxx:823
 AliRDHFCutsLctopKpi.cxx:824
 AliRDHFCutsLctopKpi.cxx:825
 AliRDHFCutsLctopKpi.cxx:826
 AliRDHFCutsLctopKpi.cxx:827
 AliRDHFCutsLctopKpi.cxx:828
 AliRDHFCutsLctopKpi.cxx:829
 AliRDHFCutsLctopKpi.cxx:830
 AliRDHFCutsLctopKpi.cxx:831
 AliRDHFCutsLctopKpi.cxx:832
 AliRDHFCutsLctopKpi.cxx:833
 AliRDHFCutsLctopKpi.cxx:834
 AliRDHFCutsLctopKpi.cxx:835
 AliRDHFCutsLctopKpi.cxx:836
 AliRDHFCutsLctopKpi.cxx:837
 AliRDHFCutsLctopKpi.cxx:838
 AliRDHFCutsLctopKpi.cxx:839
 AliRDHFCutsLctopKpi.cxx:840
 AliRDHFCutsLctopKpi.cxx:841
 AliRDHFCutsLctopKpi.cxx:842
 AliRDHFCutsLctopKpi.cxx:843
 AliRDHFCutsLctopKpi.cxx:844
 AliRDHFCutsLctopKpi.cxx:845
 AliRDHFCutsLctopKpi.cxx:846
 AliRDHFCutsLctopKpi.cxx:847
 AliRDHFCutsLctopKpi.cxx:848
 AliRDHFCutsLctopKpi.cxx:849
 AliRDHFCutsLctopKpi.cxx:850
 AliRDHFCutsLctopKpi.cxx:851
 AliRDHFCutsLctopKpi.cxx:852
 AliRDHFCutsLctopKpi.cxx:853
 AliRDHFCutsLctopKpi.cxx:854
 AliRDHFCutsLctopKpi.cxx:855
 AliRDHFCutsLctopKpi.cxx:856
 AliRDHFCutsLctopKpi.cxx:857
 AliRDHFCutsLctopKpi.cxx:858
 AliRDHFCutsLctopKpi.cxx:859
 AliRDHFCutsLctopKpi.cxx:860
 AliRDHFCutsLctopKpi.cxx:861
 AliRDHFCutsLctopKpi.cxx:862
 AliRDHFCutsLctopKpi.cxx:863
 AliRDHFCutsLctopKpi.cxx:864
 AliRDHFCutsLctopKpi.cxx:865
 AliRDHFCutsLctopKpi.cxx:866
 AliRDHFCutsLctopKpi.cxx:867
 AliRDHFCutsLctopKpi.cxx:868
 AliRDHFCutsLctopKpi.cxx:869
 AliRDHFCutsLctopKpi.cxx:870
 AliRDHFCutsLctopKpi.cxx:871
 AliRDHFCutsLctopKpi.cxx:872
 AliRDHFCutsLctopKpi.cxx:873
 AliRDHFCutsLctopKpi.cxx:874
 AliRDHFCutsLctopKpi.cxx:875
 AliRDHFCutsLctopKpi.cxx:876
 AliRDHFCutsLctopKpi.cxx:877
 AliRDHFCutsLctopKpi.cxx:878
 AliRDHFCutsLctopKpi.cxx:879
 AliRDHFCutsLctopKpi.cxx:880
 AliRDHFCutsLctopKpi.cxx:881
 AliRDHFCutsLctopKpi.cxx:882
 AliRDHFCutsLctopKpi.cxx:883
 AliRDHFCutsLctopKpi.cxx:884
 AliRDHFCutsLctopKpi.cxx:885
 AliRDHFCutsLctopKpi.cxx:886
 AliRDHFCutsLctopKpi.cxx:887
 AliRDHFCutsLctopKpi.cxx:888
 AliRDHFCutsLctopKpi.cxx:889
 AliRDHFCutsLctopKpi.cxx:890
 AliRDHFCutsLctopKpi.cxx:891
 AliRDHFCutsLctopKpi.cxx:892
 AliRDHFCutsLctopKpi.cxx:893
 AliRDHFCutsLctopKpi.cxx:894
 AliRDHFCutsLctopKpi.cxx:895
 AliRDHFCutsLctopKpi.cxx:896
 AliRDHFCutsLctopKpi.cxx:897
 AliRDHFCutsLctopKpi.cxx:898
 AliRDHFCutsLctopKpi.cxx:899
 AliRDHFCutsLctopKpi.cxx:900
 AliRDHFCutsLctopKpi.cxx:901
 AliRDHFCutsLctopKpi.cxx:902
 AliRDHFCutsLctopKpi.cxx:903
 AliRDHFCutsLctopKpi.cxx:904
 AliRDHFCutsLctopKpi.cxx:905
 AliRDHFCutsLctopKpi.cxx:906
 AliRDHFCutsLctopKpi.cxx:907
 AliRDHFCutsLctopKpi.cxx:908
 AliRDHFCutsLctopKpi.cxx:909
 AliRDHFCutsLctopKpi.cxx:910
 AliRDHFCutsLctopKpi.cxx:911
 AliRDHFCutsLctopKpi.cxx:912
 AliRDHFCutsLctopKpi.cxx:913
 AliRDHFCutsLctopKpi.cxx:914
 AliRDHFCutsLctopKpi.cxx:915
 AliRDHFCutsLctopKpi.cxx:916
 AliRDHFCutsLctopKpi.cxx:917
 AliRDHFCutsLctopKpi.cxx:918
 AliRDHFCutsLctopKpi.cxx:919
 AliRDHFCutsLctopKpi.cxx:920
 AliRDHFCutsLctopKpi.cxx:921
 AliRDHFCutsLctopKpi.cxx:922
 AliRDHFCutsLctopKpi.cxx:923
 AliRDHFCutsLctopKpi.cxx:924
 AliRDHFCutsLctopKpi.cxx:925
 AliRDHFCutsLctopKpi.cxx:926
 AliRDHFCutsLctopKpi.cxx:927
 AliRDHFCutsLctopKpi.cxx:928
 AliRDHFCutsLctopKpi.cxx:929
 AliRDHFCutsLctopKpi.cxx:930
 AliRDHFCutsLctopKpi.cxx:931
 AliRDHFCutsLctopKpi.cxx:932
 AliRDHFCutsLctopKpi.cxx:933
 AliRDHFCutsLctopKpi.cxx:934
 AliRDHFCutsLctopKpi.cxx:935
 AliRDHFCutsLctopKpi.cxx:936
 AliRDHFCutsLctopKpi.cxx:937
 AliRDHFCutsLctopKpi.cxx:938
 AliRDHFCutsLctopKpi.cxx:939
 AliRDHFCutsLctopKpi.cxx:940
 AliRDHFCutsLctopKpi.cxx:941
 AliRDHFCutsLctopKpi.cxx:942
 AliRDHFCutsLctopKpi.cxx:943
 AliRDHFCutsLctopKpi.cxx:944
 AliRDHFCutsLctopKpi.cxx:945
 AliRDHFCutsLctopKpi.cxx:946
 AliRDHFCutsLctopKpi.cxx:947
 AliRDHFCutsLctopKpi.cxx:948
 AliRDHFCutsLctopKpi.cxx:949
 AliRDHFCutsLctopKpi.cxx:950
 AliRDHFCutsLctopKpi.cxx:951
 AliRDHFCutsLctopKpi.cxx:952
 AliRDHFCutsLctopKpi.cxx:953
 AliRDHFCutsLctopKpi.cxx:954
 AliRDHFCutsLctopKpi.cxx:955
 AliRDHFCutsLctopKpi.cxx:956
 AliRDHFCutsLctopKpi.cxx:957
 AliRDHFCutsLctopKpi.cxx:958
 AliRDHFCutsLctopKpi.cxx:959
 AliRDHFCutsLctopKpi.cxx:960
 AliRDHFCutsLctopKpi.cxx:961
 AliRDHFCutsLctopKpi.cxx:962
 AliRDHFCutsLctopKpi.cxx:963
 AliRDHFCutsLctopKpi.cxx:964
 AliRDHFCutsLctopKpi.cxx:965
 AliRDHFCutsLctopKpi.cxx:966
 AliRDHFCutsLctopKpi.cxx:967
 AliRDHFCutsLctopKpi.cxx:968
 AliRDHFCutsLctopKpi.cxx:969
 AliRDHFCutsLctopKpi.cxx:970
 AliRDHFCutsLctopKpi.cxx:971
 AliRDHFCutsLctopKpi.cxx:972
 AliRDHFCutsLctopKpi.cxx:973
 AliRDHFCutsLctopKpi.cxx:974
 AliRDHFCutsLctopKpi.cxx:975
 AliRDHFCutsLctopKpi.cxx:976
 AliRDHFCutsLctopKpi.cxx:977
 AliRDHFCutsLctopKpi.cxx:978
 AliRDHFCutsLctopKpi.cxx:979
 AliRDHFCutsLctopKpi.cxx:980
 AliRDHFCutsLctopKpi.cxx:981
 AliRDHFCutsLctopKpi.cxx:982
 AliRDHFCutsLctopKpi.cxx:983
 AliRDHFCutsLctopKpi.cxx:984
 AliRDHFCutsLctopKpi.cxx:985
 AliRDHFCutsLctopKpi.cxx:986
 AliRDHFCutsLctopKpi.cxx:987
 AliRDHFCutsLctopKpi.cxx:988
 AliRDHFCutsLctopKpi.cxx:989
 AliRDHFCutsLctopKpi.cxx:990
 AliRDHFCutsLctopKpi.cxx:991
 AliRDHFCutsLctopKpi.cxx:992
 AliRDHFCutsLctopKpi.cxx:993
 AliRDHFCutsLctopKpi.cxx:994
 AliRDHFCutsLctopKpi.cxx:995
 AliRDHFCutsLctopKpi.cxx:996
 AliRDHFCutsLctopKpi.cxx:997
 AliRDHFCutsLctopKpi.cxx:998
 AliRDHFCutsLctopKpi.cxx:999
 AliRDHFCutsLctopKpi.cxx:1000
 AliRDHFCutsLctopKpi.cxx:1001
 AliRDHFCutsLctopKpi.cxx:1002
 AliRDHFCutsLctopKpi.cxx:1003
 AliRDHFCutsLctopKpi.cxx:1004
 AliRDHFCutsLctopKpi.cxx:1005
 AliRDHFCutsLctopKpi.cxx:1006
 AliRDHFCutsLctopKpi.cxx:1007
 AliRDHFCutsLctopKpi.cxx:1008
 AliRDHFCutsLctopKpi.cxx:1009
 AliRDHFCutsLctopKpi.cxx:1010
 AliRDHFCutsLctopKpi.cxx:1011
 AliRDHFCutsLctopKpi.cxx:1012
 AliRDHFCutsLctopKpi.cxx:1013
 AliRDHFCutsLctopKpi.cxx:1014
 AliRDHFCutsLctopKpi.cxx:1015
 AliRDHFCutsLctopKpi.cxx:1016
 AliRDHFCutsLctopKpi.cxx:1017
 AliRDHFCutsLctopKpi.cxx:1018
 AliRDHFCutsLctopKpi.cxx:1019
 AliRDHFCutsLctopKpi.cxx:1020
 AliRDHFCutsLctopKpi.cxx:1021
 AliRDHFCutsLctopKpi.cxx:1022
 AliRDHFCutsLctopKpi.cxx:1023
 AliRDHFCutsLctopKpi.cxx:1024
 AliRDHFCutsLctopKpi.cxx:1025
 AliRDHFCutsLctopKpi.cxx:1026
 AliRDHFCutsLctopKpi.cxx:1027
 AliRDHFCutsLctopKpi.cxx:1028
 AliRDHFCutsLctopKpi.cxx:1029
 AliRDHFCutsLctopKpi.cxx:1030
 AliRDHFCutsLctopKpi.cxx:1031
 AliRDHFCutsLctopKpi.cxx:1032
 AliRDHFCutsLctopKpi.cxx:1033
 AliRDHFCutsLctopKpi.cxx:1034
 AliRDHFCutsLctopKpi.cxx:1035
 AliRDHFCutsLctopKpi.cxx:1036
 AliRDHFCutsLctopKpi.cxx:1037
 AliRDHFCutsLctopKpi.cxx:1038
 AliRDHFCutsLctopKpi.cxx:1039
 AliRDHFCutsLctopKpi.cxx:1040
 AliRDHFCutsLctopKpi.cxx:1041
 AliRDHFCutsLctopKpi.cxx:1042
 AliRDHFCutsLctopKpi.cxx:1043
 AliRDHFCutsLctopKpi.cxx:1044
 AliRDHFCutsLctopKpi.cxx:1045
 AliRDHFCutsLctopKpi.cxx:1046
 AliRDHFCutsLctopKpi.cxx:1047
 AliRDHFCutsLctopKpi.cxx:1048
 AliRDHFCutsLctopKpi.cxx:1049
 AliRDHFCutsLctopKpi.cxx:1050
 AliRDHFCutsLctopKpi.cxx:1051
 AliRDHFCutsLctopKpi.cxx:1052
 AliRDHFCutsLctopKpi.cxx:1053
 AliRDHFCutsLctopKpi.cxx:1054
 AliRDHFCutsLctopKpi.cxx:1055
 AliRDHFCutsLctopKpi.cxx:1056
 AliRDHFCutsLctopKpi.cxx:1057
 AliRDHFCutsLctopKpi.cxx:1058
 AliRDHFCutsLctopKpi.cxx:1059
 AliRDHFCutsLctopKpi.cxx:1060
 AliRDHFCutsLctopKpi.cxx:1061
 AliRDHFCutsLctopKpi.cxx:1062
 AliRDHFCutsLctopKpi.cxx:1063
 AliRDHFCutsLctopKpi.cxx:1064
 AliRDHFCutsLctopKpi.cxx:1065
 AliRDHFCutsLctopKpi.cxx:1066
 AliRDHFCutsLctopKpi.cxx:1067
 AliRDHFCutsLctopKpi.cxx:1068
 AliRDHFCutsLctopKpi.cxx:1069
 AliRDHFCutsLctopKpi.cxx:1070
 AliRDHFCutsLctopKpi.cxx:1071
 AliRDHFCutsLctopKpi.cxx:1072
 AliRDHFCutsLctopKpi.cxx:1073
 AliRDHFCutsLctopKpi.cxx:1074
 AliRDHFCutsLctopKpi.cxx:1075
 AliRDHFCutsLctopKpi.cxx:1076
 AliRDHFCutsLctopKpi.cxx:1077
 AliRDHFCutsLctopKpi.cxx:1078
 AliRDHFCutsLctopKpi.cxx:1079
 AliRDHFCutsLctopKpi.cxx:1080
 AliRDHFCutsLctopKpi.cxx:1081
 AliRDHFCutsLctopKpi.cxx:1082
 AliRDHFCutsLctopKpi.cxx:1083
 AliRDHFCutsLctopKpi.cxx:1084
 AliRDHFCutsLctopKpi.cxx:1085
 AliRDHFCutsLctopKpi.cxx:1086
 AliRDHFCutsLctopKpi.cxx:1087
 AliRDHFCutsLctopKpi.cxx:1088
 AliRDHFCutsLctopKpi.cxx:1089
 AliRDHFCutsLctopKpi.cxx:1090
 AliRDHFCutsLctopKpi.cxx:1091
 AliRDHFCutsLctopKpi.cxx:1092
 AliRDHFCutsLctopKpi.cxx:1093
 AliRDHFCutsLctopKpi.cxx:1094
 AliRDHFCutsLctopKpi.cxx:1095
 AliRDHFCutsLctopKpi.cxx:1096
 AliRDHFCutsLctopKpi.cxx:1097
 AliRDHFCutsLctopKpi.cxx:1098
 AliRDHFCutsLctopKpi.cxx:1099
 AliRDHFCutsLctopKpi.cxx:1100
 AliRDHFCutsLctopKpi.cxx:1101
 AliRDHFCutsLctopKpi.cxx:1102
 AliRDHFCutsLctopKpi.cxx:1103
 AliRDHFCutsLctopKpi.cxx:1104
 AliRDHFCutsLctopKpi.cxx:1105
 AliRDHFCutsLctopKpi.cxx:1106
 AliRDHFCutsLctopKpi.cxx:1107
 AliRDHFCutsLctopKpi.cxx:1108
 AliRDHFCutsLctopKpi.cxx:1109
 AliRDHFCutsLctopKpi.cxx:1110
 AliRDHFCutsLctopKpi.cxx:1111
 AliRDHFCutsLctopKpi.cxx:1112
 AliRDHFCutsLctopKpi.cxx:1113
 AliRDHFCutsLctopKpi.cxx:1114
 AliRDHFCutsLctopKpi.cxx:1115
 AliRDHFCutsLctopKpi.cxx:1116
 AliRDHFCutsLctopKpi.cxx:1117
 AliRDHFCutsLctopKpi.cxx:1118
 AliRDHFCutsLctopKpi.cxx:1119
 AliRDHFCutsLctopKpi.cxx:1120
 AliRDHFCutsLctopKpi.cxx:1121
 AliRDHFCutsLctopKpi.cxx:1122
 AliRDHFCutsLctopKpi.cxx:1123
 AliRDHFCutsLctopKpi.cxx:1124
 AliRDHFCutsLctopKpi.cxx:1125
 AliRDHFCutsLctopKpi.cxx:1126
 AliRDHFCutsLctopKpi.cxx:1127
 AliRDHFCutsLctopKpi.cxx:1128
 AliRDHFCutsLctopKpi.cxx:1129
 AliRDHFCutsLctopKpi.cxx:1130
 AliRDHFCutsLctopKpi.cxx:1131
 AliRDHFCutsLctopKpi.cxx:1132
 AliRDHFCutsLctopKpi.cxx:1133
 AliRDHFCutsLctopKpi.cxx:1134
 AliRDHFCutsLctopKpi.cxx:1135
 AliRDHFCutsLctopKpi.cxx:1136
 AliRDHFCutsLctopKpi.cxx:1137
 AliRDHFCutsLctopKpi.cxx:1138
 AliRDHFCutsLctopKpi.cxx:1139
 AliRDHFCutsLctopKpi.cxx:1140
 AliRDHFCutsLctopKpi.cxx:1141
 AliRDHFCutsLctopKpi.cxx:1142
 AliRDHFCutsLctopKpi.cxx:1143
 AliRDHFCutsLctopKpi.cxx:1144
 AliRDHFCutsLctopKpi.cxx:1145
 AliRDHFCutsLctopKpi.cxx:1146
 AliRDHFCutsLctopKpi.cxx:1147
 AliRDHFCutsLctopKpi.cxx:1148
 AliRDHFCutsLctopKpi.cxx:1149
 AliRDHFCutsLctopKpi.cxx:1150
 AliRDHFCutsLctopKpi.cxx:1151
 AliRDHFCutsLctopKpi.cxx:1152
 AliRDHFCutsLctopKpi.cxx:1153
 AliRDHFCutsLctopKpi.cxx:1154
 AliRDHFCutsLctopKpi.cxx:1155
 AliRDHFCutsLctopKpi.cxx:1156
 AliRDHFCutsLctopKpi.cxx:1157
 AliRDHFCutsLctopKpi.cxx:1158
 AliRDHFCutsLctopKpi.cxx:1159
 AliRDHFCutsLctopKpi.cxx:1160
 AliRDHFCutsLctopKpi.cxx:1161
 AliRDHFCutsLctopKpi.cxx:1162
 AliRDHFCutsLctopKpi.cxx:1163
 AliRDHFCutsLctopKpi.cxx:1164
 AliRDHFCutsLctopKpi.cxx:1165
 AliRDHFCutsLctopKpi.cxx:1166
 AliRDHFCutsLctopKpi.cxx:1167
 AliRDHFCutsLctopKpi.cxx:1168
 AliRDHFCutsLctopKpi.cxx:1169
 AliRDHFCutsLctopKpi.cxx:1170
 AliRDHFCutsLctopKpi.cxx:1171
 AliRDHFCutsLctopKpi.cxx:1172
 AliRDHFCutsLctopKpi.cxx:1173
 AliRDHFCutsLctopKpi.cxx:1174
 AliRDHFCutsLctopKpi.cxx:1175
 AliRDHFCutsLctopKpi.cxx:1176
 AliRDHFCutsLctopKpi.cxx:1177
 AliRDHFCutsLctopKpi.cxx:1178
 AliRDHFCutsLctopKpi.cxx:1179
 AliRDHFCutsLctopKpi.cxx:1180
 AliRDHFCutsLctopKpi.cxx:1181
 AliRDHFCutsLctopKpi.cxx:1182
 AliRDHFCutsLctopKpi.cxx:1183
 AliRDHFCutsLctopKpi.cxx:1184
 AliRDHFCutsLctopKpi.cxx:1185
 AliRDHFCutsLctopKpi.cxx:1186
 AliRDHFCutsLctopKpi.cxx:1187
 AliRDHFCutsLctopKpi.cxx:1188
 AliRDHFCutsLctopKpi.cxx:1189
 AliRDHFCutsLctopKpi.cxx:1190
 AliRDHFCutsLctopKpi.cxx:1191
 AliRDHFCutsLctopKpi.cxx:1192
 AliRDHFCutsLctopKpi.cxx:1193
 AliRDHFCutsLctopKpi.cxx:1194
 AliRDHFCutsLctopKpi.cxx:1195
 AliRDHFCutsLctopKpi.cxx:1196
 AliRDHFCutsLctopKpi.cxx:1197
 AliRDHFCutsLctopKpi.cxx:1198
 AliRDHFCutsLctopKpi.cxx:1199
 AliRDHFCutsLctopKpi.cxx:1200
 AliRDHFCutsLctopKpi.cxx:1201
 AliRDHFCutsLctopKpi.cxx:1202
 AliRDHFCutsLctopKpi.cxx:1203
 AliRDHFCutsLctopKpi.cxx:1204
 AliRDHFCutsLctopKpi.cxx:1205
 AliRDHFCutsLctopKpi.cxx:1206
 AliRDHFCutsLctopKpi.cxx:1207
 AliRDHFCutsLctopKpi.cxx:1208
 AliRDHFCutsLctopKpi.cxx:1209
 AliRDHFCutsLctopKpi.cxx:1210
 AliRDHFCutsLctopKpi.cxx:1211
 AliRDHFCutsLctopKpi.cxx:1212
 AliRDHFCutsLctopKpi.cxx:1213
 AliRDHFCutsLctopKpi.cxx:1214
 AliRDHFCutsLctopKpi.cxx:1215
 AliRDHFCutsLctopKpi.cxx:1216
 AliRDHFCutsLctopKpi.cxx:1217
 AliRDHFCutsLctopKpi.cxx:1218
 AliRDHFCutsLctopKpi.cxx:1219
 AliRDHFCutsLctopKpi.cxx:1220
 AliRDHFCutsLctopKpi.cxx:1221
 AliRDHFCutsLctopKpi.cxx:1222
 AliRDHFCutsLctopKpi.cxx:1223
 AliRDHFCutsLctopKpi.cxx:1224
 AliRDHFCutsLctopKpi.cxx:1225
 AliRDHFCutsLctopKpi.cxx:1226
 AliRDHFCutsLctopKpi.cxx:1227
 AliRDHFCutsLctopKpi.cxx:1228
 AliRDHFCutsLctopKpi.cxx:1229
 AliRDHFCutsLctopKpi.cxx:1230
 AliRDHFCutsLctopKpi.cxx:1231
 AliRDHFCutsLctopKpi.cxx:1232
 AliRDHFCutsLctopKpi.cxx:1233
 AliRDHFCutsLctopKpi.cxx:1234
 AliRDHFCutsLctopKpi.cxx:1235
 AliRDHFCutsLctopKpi.cxx:1236
 AliRDHFCutsLctopKpi.cxx:1237
 AliRDHFCutsLctopKpi.cxx:1238
 AliRDHFCutsLctopKpi.cxx:1239
 AliRDHFCutsLctopKpi.cxx:1240
 AliRDHFCutsLctopKpi.cxx:1241
 AliRDHFCutsLctopKpi.cxx:1242
 AliRDHFCutsLctopKpi.cxx:1243
 AliRDHFCutsLctopKpi.cxx:1244
 AliRDHFCutsLctopKpi.cxx:1245
 AliRDHFCutsLctopKpi.cxx:1246
 AliRDHFCutsLctopKpi.cxx:1247
 AliRDHFCutsLctopKpi.cxx:1248
 AliRDHFCutsLctopKpi.cxx:1249
 AliRDHFCutsLctopKpi.cxx:1250
 AliRDHFCutsLctopKpi.cxx:1251
 AliRDHFCutsLctopKpi.cxx:1252
 AliRDHFCutsLctopKpi.cxx:1253
 AliRDHFCutsLctopKpi.cxx:1254
 AliRDHFCutsLctopKpi.cxx:1255
 AliRDHFCutsLctopKpi.cxx:1256
 AliRDHFCutsLctopKpi.cxx:1257
 AliRDHFCutsLctopKpi.cxx:1258
 AliRDHFCutsLctopKpi.cxx:1259
 AliRDHFCutsLctopKpi.cxx:1260
 AliRDHFCutsLctopKpi.cxx:1261
 AliRDHFCutsLctopKpi.cxx:1262
 AliRDHFCutsLctopKpi.cxx:1263
 AliRDHFCutsLctopKpi.cxx:1264
 AliRDHFCutsLctopKpi.cxx:1265
 AliRDHFCutsLctopKpi.cxx:1266
 AliRDHFCutsLctopKpi.cxx:1267
 AliRDHFCutsLctopKpi.cxx:1268
 AliRDHFCutsLctopKpi.cxx:1269
 AliRDHFCutsLctopKpi.cxx:1270
 AliRDHFCutsLctopKpi.cxx:1271
 AliRDHFCutsLctopKpi.cxx:1272
 AliRDHFCutsLctopKpi.cxx:1273
 AliRDHFCutsLctopKpi.cxx:1274
 AliRDHFCutsLctopKpi.cxx:1275
 AliRDHFCutsLctopKpi.cxx:1276
 AliRDHFCutsLctopKpi.cxx:1277
 AliRDHFCutsLctopKpi.cxx:1278
 AliRDHFCutsLctopKpi.cxx:1279
 AliRDHFCutsLctopKpi.cxx:1280
 AliRDHFCutsLctopKpi.cxx:1281
 AliRDHFCutsLctopKpi.cxx:1282
 AliRDHFCutsLctopKpi.cxx:1283
 AliRDHFCutsLctopKpi.cxx:1284
 AliRDHFCutsLctopKpi.cxx:1285
 AliRDHFCutsLctopKpi.cxx:1286
 AliRDHFCutsLctopKpi.cxx:1287
 AliRDHFCutsLctopKpi.cxx:1288
 AliRDHFCutsLctopKpi.cxx:1289
 AliRDHFCutsLctopKpi.cxx:1290
 AliRDHFCutsLctopKpi.cxx:1291
 AliRDHFCutsLctopKpi.cxx:1292
 AliRDHFCutsLctopKpi.cxx:1293
 AliRDHFCutsLctopKpi.cxx:1294
 AliRDHFCutsLctopKpi.cxx:1295
 AliRDHFCutsLctopKpi.cxx:1296
 AliRDHFCutsLctopKpi.cxx:1297
 AliRDHFCutsLctopKpi.cxx:1298
 AliRDHFCutsLctopKpi.cxx:1299
 AliRDHFCutsLctopKpi.cxx:1300
 AliRDHFCutsLctopKpi.cxx:1301
 AliRDHFCutsLctopKpi.cxx:1302
 AliRDHFCutsLctopKpi.cxx:1303
 AliRDHFCutsLctopKpi.cxx:1304
 AliRDHFCutsLctopKpi.cxx:1305
 AliRDHFCutsLctopKpi.cxx:1306
 AliRDHFCutsLctopKpi.cxx:1307
 AliRDHFCutsLctopKpi.cxx:1308
 AliRDHFCutsLctopKpi.cxx:1309
 AliRDHFCutsLctopKpi.cxx:1310
 AliRDHFCutsLctopKpi.cxx:1311
 AliRDHFCutsLctopKpi.cxx:1312
 AliRDHFCutsLctopKpi.cxx:1313
 AliRDHFCutsLctopKpi.cxx:1314
 AliRDHFCutsLctopKpi.cxx:1315
 AliRDHFCutsLctopKpi.cxx:1316
 AliRDHFCutsLctopKpi.cxx:1317
 AliRDHFCutsLctopKpi.cxx:1318
 AliRDHFCutsLctopKpi.cxx:1319
 AliRDHFCutsLctopKpi.cxx:1320
 AliRDHFCutsLctopKpi.cxx:1321
 AliRDHFCutsLctopKpi.cxx:1322
 AliRDHFCutsLctopKpi.cxx:1323
 AliRDHFCutsLctopKpi.cxx:1324
 AliRDHFCutsLctopKpi.cxx:1325
 AliRDHFCutsLctopKpi.cxx:1326
 AliRDHFCutsLctopKpi.cxx:1327
 AliRDHFCutsLctopKpi.cxx:1328
 AliRDHFCutsLctopKpi.cxx:1329
 AliRDHFCutsLctopKpi.cxx:1330
 AliRDHFCutsLctopKpi.cxx:1331
 AliRDHFCutsLctopKpi.cxx:1332
 AliRDHFCutsLctopKpi.cxx:1333
 AliRDHFCutsLctopKpi.cxx:1334
 AliRDHFCutsLctopKpi.cxx:1335
 AliRDHFCutsLctopKpi.cxx:1336
 AliRDHFCutsLctopKpi.cxx:1337
 AliRDHFCutsLctopKpi.cxx:1338
 AliRDHFCutsLctopKpi.cxx:1339
 AliRDHFCutsLctopKpi.cxx:1340
 AliRDHFCutsLctopKpi.cxx:1341
 AliRDHFCutsLctopKpi.cxx:1342
 AliRDHFCutsLctopKpi.cxx:1343
 AliRDHFCutsLctopKpi.cxx:1344
 AliRDHFCutsLctopKpi.cxx:1345
 AliRDHFCutsLctopKpi.cxx:1346
 AliRDHFCutsLctopKpi.cxx:1347
 AliRDHFCutsLctopKpi.cxx:1348
 AliRDHFCutsLctopKpi.cxx:1349
 AliRDHFCutsLctopKpi.cxx:1350
 AliRDHFCutsLctopKpi.cxx:1351
 AliRDHFCutsLctopKpi.cxx:1352
 AliRDHFCutsLctopKpi.cxx:1353
 AliRDHFCutsLctopKpi.cxx:1354
 AliRDHFCutsLctopKpi.cxx:1355
 AliRDHFCutsLctopKpi.cxx:1356
 AliRDHFCutsLctopKpi.cxx:1357
 AliRDHFCutsLctopKpi.cxx:1358
 AliRDHFCutsLctopKpi.cxx:1359
 AliRDHFCutsLctopKpi.cxx:1360
 AliRDHFCutsLctopKpi.cxx:1361
 AliRDHFCutsLctopKpi.cxx:1362
 AliRDHFCutsLctopKpi.cxx:1363
 AliRDHFCutsLctopKpi.cxx:1364
 AliRDHFCutsLctopKpi.cxx:1365
 AliRDHFCutsLctopKpi.cxx:1366
 AliRDHFCutsLctopKpi.cxx:1367
 AliRDHFCutsLctopKpi.cxx:1368
 AliRDHFCutsLctopKpi.cxx:1369
 AliRDHFCutsLctopKpi.cxx:1370
 AliRDHFCutsLctopKpi.cxx:1371
 AliRDHFCutsLctopKpi.cxx:1372
 AliRDHFCutsLctopKpi.cxx:1373
 AliRDHFCutsLctopKpi.cxx:1374
 AliRDHFCutsLctopKpi.cxx:1375
 AliRDHFCutsLctopKpi.cxx:1376
 AliRDHFCutsLctopKpi.cxx:1377
 AliRDHFCutsLctopKpi.cxx:1378
 AliRDHFCutsLctopKpi.cxx:1379
 AliRDHFCutsLctopKpi.cxx:1380
 AliRDHFCutsLctopKpi.cxx:1381
 AliRDHFCutsLctopKpi.cxx:1382
 AliRDHFCutsLctopKpi.cxx:1383
 AliRDHFCutsLctopKpi.cxx:1384
 AliRDHFCutsLctopKpi.cxx:1385
 AliRDHFCutsLctopKpi.cxx:1386
 AliRDHFCutsLctopKpi.cxx:1387
 AliRDHFCutsLctopKpi.cxx:1388
 AliRDHFCutsLctopKpi.cxx:1389
 AliRDHFCutsLctopKpi.cxx:1390
 AliRDHFCutsLctopKpi.cxx:1391
 AliRDHFCutsLctopKpi.cxx:1392
 AliRDHFCutsLctopKpi.cxx:1393
 AliRDHFCutsLctopKpi.cxx:1394
 AliRDHFCutsLctopKpi.cxx:1395
 AliRDHFCutsLctopKpi.cxx:1396
 AliRDHFCutsLctopKpi.cxx:1397
 AliRDHFCutsLctopKpi.cxx:1398
 AliRDHFCutsLctopKpi.cxx:1399
 AliRDHFCutsLctopKpi.cxx:1400
 AliRDHFCutsLctopKpi.cxx:1401
 AliRDHFCutsLctopKpi.cxx:1402
 AliRDHFCutsLctopKpi.cxx:1403