Remove unnecessary classes
authorsma <sma@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 21 Apr 2011 12:40:11 +0000 (12:40 +0000)
committersma <sma@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 21 Apr 2011 12:40:11 +0000 (12:40 +0000)
PWG3/hfe/AliAnalysisTaskCheckV0tender.cxx [deleted file]
PWG3/hfe/AliAnalysisTaskCheckV0tender.h [deleted file]
PWG3/hfe/AliAnalysisTaskCheckV0tenderII.cxx [deleted file]
PWG3/hfe/AliAnalysisTaskCheckV0tenderII.h [deleted file]
PWG3/hfe/AliESDv0KineCuts.cxx [deleted file]
PWG3/hfe/AliESDv0KineCuts.h [deleted file]

diff --git a/PWG3/hfe/AliAnalysisTaskCheckV0tender.cxx b/PWG3/hfe/AliAnalysisTaskCheckV0tender.cxx
deleted file mode 100644 (file)
index 508213a..0000000
+++ /dev/null
@@ -1,478 +0,0 @@
-/**************************************************************************
-* Copyright(c) 1998-1999, 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$ */
-
-//
-// Task fir checking the performance of the V0 tender
-// 
-// 
-// Authors
-//   Matus Kalisky <matus.kalisky@cern.ch>
-//
-
-#include <sstream>
-
-#include <TH1F.h>
-#include <TList.h>
-
-#include "AliAnalysisManager.h"
-#include "AliMCEventHandler.h"
-#include "AliESDInputHandler.h"
-#include "AliESDv0.h"
-#include "AliESDtrack.h"
-#include "AliMCParticle.h"
-#include "AliMCEvent.h"
-
-#include "AliHFEtools.h"
-#include "AliHFEcollection.h"
-
-#include "AliAnalysisTaskCheckV0tender.h"
-
-ClassImp(AliAnalysisTaskCheckV0tender)
-
-//__________________________________________________________
-AliAnalysisTaskCheckV0tender::AliAnalysisTaskCheckV0tender():
-  AliAnalysisTaskSE("CheckV0tenderTask")
-  , fOutput(0x0)
-  , fColl(0x0)
-  , fCollMC(0x0)
-  , fEvents(0x0)
-{
-  //
-  // Default Constructor
-  //
-}
-//__________________________________________________________
-AliAnalysisTaskCheckV0tender::AliAnalysisTaskCheckV0tender(const Char_t *name):
-  AliAnalysisTaskSE(name)
-  , fOutput(0x0)
-  , fColl(0x0)
-  , fCollMC(0x0)
-  , fEvents(0x0)
-{
-  //
-  // Default Constructor
-  //
-  DefineOutput(1, TH1F::Class());
-  DefineOutput(2, TList::Class());
-
-}
-//__________________________________________________________
-AliAnalysisTaskCheckV0tender::~AliAnalysisTaskCheckV0tender(){
-  //
-  // Destructor
-  //
-
-  if (fOutput) delete fOutput;
-  //if (fColl) delete fColl;
-  //if (fCollMC) delete fCollMC;
-  if (fEvents) delete fEvents;
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tender::UserCreateOutputObjects(){
-  //
-  // prepare output objects
-  //
-
-  fOutput = new TList();
-  fOutput->SetOwner();
-  // Counter for number of events
-  fEvents = new TH1I("nEvents", "NumberOfEvents", 1, 1, 2);
-
-  fColl = new AliHFEcollection("V0_QA", "tender V0s for data");
-  fCollMC = new AliHFEcollection("V0_MC_QA", "tender V0s for MC");
-
-  // 
-  // Data histos
-  //
-  
-  // total number of tagged V0s
-  fColl->CreateTH1F("h_NumberOf_V0s", "Number of tagged V0s; type; counts", 4, -0.5, 3.5);
-  // pT spectra of the tagged V0s
-  fColl->CreateTH1F("h_Gamma_pt", "p_{T} spectrum of tagged gammas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_K0_pt", "p_{T} spectrum of tagged K0s; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_Lambda_pt", "p_{T} spectrum of tagged Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_ALambda_pt", "p_{T} spectrum of tagged A-Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  // invariant mass of the V0s
-  fColl->CreateTH1F("h_Gamma_mass", "Inv. Mass of the ", 100, 0., 0.1);
-  fColl->CreateTH1F("h_K0_mass", "Inv. Mass of the ", 100, 0.45, 0.55);
-  fColl->CreateTH1F("h_Lambda_mass", "Inv. Mass of the ", 100, 1.08, 1.14);
-  fColl->CreateTH1F("h_ALambda_mass", "Inv. Mass of the ", 100, 1.08, 1.14);
-
-  // total number of tagged daughter particles (should correlate with number of V0s !
-  fColl->CreateTH1F("h_NumberOfDaughters", "Number of tagged daughters; type; counts", 3, -0.5, 2.5);
-  // pT spectra of tagged daughter particles
-  fColl->CreateTH1F("h_Electron_pt", "p_{T} spectrum of tagged; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_Pion_pt", "p_{T} spectrum of tagged; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_Proton_pt", "p_{T} spectrum of tagged; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-
-  //
-  // MC histos
-  //
-
-  // pT spectra of the tagged V0s
-  fCollMC->CreateTH1F("h_Gamma_pt_S", "MC-S: p_{T} spectrum of tagged gammas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_K0_pt_S", "MC-S: p_{T} spectrum of tagged K0s; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Lambda_pt_S", "MC-S: p_{T} spectrum of tagged Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_ALambda_pt_S", "MC-S: p_{T} spectrum of tagged A-Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0); 
-  fCollMC->CreateTH1F("h_Gamma_pt_B", "MC-B: p_{T} spectrum of tagged gammas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_K0_pt_B", "MC-B: p_{T} spectrum of tagged K0s; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Lambda_pt_B", "MC-B: p_{T} spectrum of tagged Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_ALambda_pt_B", "MC-B: p_{T} spectrum of tagged A-Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0); 
-  // invariant mass of the V0s
-  fCollMC->CreateTH1F("h_Gamma_mass_S", "MC-S: Inv. Mass of the gamma; m (GeV/c^{2}); counts", 100, 0., 0.1);
-  fCollMC->CreateTH1F("h_K0_mass_S", "MC-S: Inv. Mass of the K0; m (GeV/c^{2}); counts", 100, 0.45, 0.55);
-  fCollMC->CreateTH1F("h_Lambda_mass_S", "MC-S: Inv. Mass of the Lambda; m (GeV/c^{2}); counts", 100, 1.08, 1.14);
-  fCollMC->CreateTH1F("h_ALambda_mass_S", "MC-S: Inv. Mass of the A-Lambda; m (GeV/c^{2}); counts", 100, 1.08, 1.14);
-  fCollMC->CreateTH1F("h_Gamma_mass_B", "MC-B: Inv. Mass of the gamma; m (GeV/c^{2}); counts", 100, 0., 0.1);
-  fCollMC->CreateTH1F("h_K0_mass_B", "MC-B: Inv. Mass of the K0; m (GeV/c^{2}); counts", 100, 0.45, 0.55);
-  fCollMC->CreateTH1F("h_Lambda_mass_B", "MC-B: Inv. Mass of the Lambda; m (GeV/c^{2}); counts", 100, 1.08, 1.14);
-  fCollMC->CreateTH1F("h_ALambda_mass_B", "MC-B: Inv. Mass of the A-Lambda; m (GeV/c^{2}); counts", 100, 1.08, 1.14);
-  // pT spectra of tagged daughter particles
-  fCollMC->CreateTH1F("h_Electron_pt_S", "MC-S: p_{T} spectrum of tagged electrons; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Pion_pt_S", "MC-S: p_{T} spectrum of tagged pions; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Proton_pt_S", "MC-S: p_{T} spectrum of tagged protons; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Electron_pt_B", "MC-B: p_{T} spectrum of tagged electrons; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Pion_pt_B", "MC-B: p_{T} spectrum of tagged pions; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Proton_pt_B", "MC-B: p_{T} spectrum of tagged protons; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-
-
-  TList *tmp = fColl->GetList();
-  tmp->SetName(fColl->GetName());
-  fOutput->Add(tmp);
-  tmp = 0x0;
-  tmp = fCollMC->GetList();
-  tmp->SetName(fCollMC->GetName());
-  fOutput->Add(tmp);
-  
-  
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tender::UserExec(Option_t *){
-  //
-  // Event Loop
-  // 
-
-  AliMCEventHandler* mcHandler = (dynamic_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()));
-  //AliESDInputHandler *inh = dynamic_cast<AliESDInputHandler *>(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
-  //AliESDpid *workingPID = NULL;
-  //if(inh && (workingPID = inh->GetESDpid())) workingPID = inh->GetESDpid();
-  //else workingPID = AliHFEtools::GetDefaultPID(mcHandler ? kTRUE : kFALSE);
-     
-  // check the MC data
-  if(fMCEvent && !mcHandler ) return;
-  if(fMCEvent &&  !mcHandler->InitOk() ) return;
-  if(fMCEvent &&  !mcHandler->TreeK() ) return;
-  if(fMCEvent &&  !mcHandler->TreeTR() ) return;
-
-  ProcessV0s();
-  ProcessDaughters();
-
-  if(fMCEvent){
-    ProcessV0sMC();
-    ProcessDaughtersMC();
-  }
-
-  fEvents->Fill(1.1);
-  PostData(1, fEvents);
-  PostData(2, fOutput);
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tender::Terminate(Option_t *){
-  //
-  // Do Post Processing
-  //
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tender::ProcessV0s(){
-  //
-  // loop over the V0s and extract the information about
-  // the V0s tagged by the V0 tender
-  //
-
-  const TString type[4] = {"Gamma", "K0", "Lambda", "ALambda"};
-  TString name;
-
-  Int_t nV0s = fInputEvent->GetNumberOfV0s();
-  for(Int_t i=0; i<nV0s; ++i){
-    AliESDv0 *esdV0 = (static_cast<AliESDEvent *>(fInputEvent))->GetV0(i);
-    if(!esdV0) continue;
-    if(!esdV0->GetOnFlyStatus()) continue; // Take only V0s from the On-the-fly v0 finder
-    Int_t pid = GetTenderPidV0(esdV0);
-    if(pid < 0) continue;
-    fColl->Fill("h_NumberOf_V0s", pid);
-    Float_t pT = esdV0->Pt();
-    name = "h_" + type[pid] + "_pt";
-    fColl->Fill(name, pT);
-    Float_t mass = MassV0(esdV0, pid);    
-    name = "h_" + type[pid] + "_mass";
-    //printf(" -D: name: %s \n", name.Data());
-    fColl->Fill(name, mass);
-  }
-
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tender::ProcessDaughters(){
-  //
-  // produce some check plots for V0 tender tagged single tracks
-  //
-
-  const TString type[3] = {"Electron", "Pion", "Proton"};
-  TString name;
-
-  Int_t nTracks = fInputEvent->GetNumberOfTracks();
-  for(Int_t i=0; i<nTracks; ++i){
-    AliESDtrack *track = dynamic_cast<AliESDtrack*>(fInputEvent->GetTrack(i));
-    if(!track) continue;
-    Int_t  pid = GetTenderPidDaughter(track);
-    if(pid < 0) continue;
-    fColl->Fill("h_NumberOfDaughters", pid*1.0);
-    Float_t pT = track->Pt();
-    name = "h_" + type[pid] + "_pt";
-    fColl->Fill(name, pT);
-    
-  }
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tender::ProcessV0sMC(){
-  //
-  // check all V0tender selected V0 on their true identity
-  // 
-
-  const Int_t pid2pdg[4] = {22, 310, 3122, -3122};
-  const TString type[4] = {"Gamma", "K0", "Lambda", "ALambda"};
-  TString name;
-  Int_t nV0s = fInputEvent->GetNumberOfV0s();
-
-  Int_t nTracks = fInputEvent->GetNumberOfTracks();
-
-  // V0 loop
-  for(Int_t i=0; i<nV0s; ++i){
-    Bool_t id = kFALSE;
-    AliESDv0 *esdV0 = (static_cast<AliESDEvent *>(fInputEvent))->GetV0(i);
-    if(!esdV0) continue;
-    if(!esdV0->GetOnFlyStatus()) continue; // Take only V0s from the On-the-fly v0 finder
-    Int_t pid = GetTenderPidV0(esdV0);
-    if(pid < 0) continue;
-
-    // both ESD daughtr tracks
-    Int_t iN, iP;
-    iN = iP = -1;
-    iP = esdV0->GetPindex();
-    iN = esdV0->GetNindex();    
-    if(iN < 0 || iP < 0) continue;
-    if(iN >= nTracks || iP >= nTracks) continue;    
-    AliESDtrack *dP, *dN;
-    dP = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(iP));
-    dN = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(iN));  
-    if(!dN || !dP) continue;
-
-    // MC labels of the daughter tracks
-    Int_t lN, lP;
-    lN = dN->GetLabel();
-    lP = dP->GetLabel();
-    if(lN < 0 || lP < 0) continue;
-
-    // MC daughter particles
-    AliMCParticle *mcP, *mcN;
-    mcP = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(lP));
-    mcN = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(lN));
-    if(!mcP || !mcN) continue;
-
-    // labels of the mother particles
-    Int_t lPm, lNm;
-    lPm = mcP->GetMother();
-    lNm = mcN->GetMother();
-    if(lPm < 0) continue;
-    AliMCParticle *m = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(lPm)); 
-    if(!m) continue;
-    Int_t pdg = m->PdgCode();
-    if((lPm == lNm) && (pdg == pid2pdg[pid])) id = kTRUE;
-
-    if(id) name = "h_" + type[pid] + "_pt_S"; 
-    else name = "h_" + type[pid] + "_pt_B"; 
-    Float_t pT = esdV0->Pt();
-    fCollMC->Fill(name, pT);
-
-    if(id) name = "h_" + type[pid] + "_mass_S";
-    else name = "h_" + type[pid] + "_mass_B";
-    Float_t mass = MassV0(esdV0, pid);
-    fCollMC->Fill(name, mass);
-
-   }
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tender::ProcessDaughtersMC(){
-  //
-  // check the identity of the V0tender selected V0 daughters
-  // !!! for positive check only the true identity plays a role here, 
-  // not the true V0 mother identity (e.g. selected electron could come
-  // from primary vertex or pion dalitz deca or true gamma conversion) !!!
-  //
-
-  const Int_t pid2pdg [3] = {11, 211, 2212};
-  const TString type[3] = {"Electron", "Pion", "Proton"};
-  TString name;
-
-
-  Int_t nTracks = fInputEvent->GetNumberOfTracks();
-  for(Int_t i=0; i<nTracks; ++i){
-    Bool_t id = kFALSE;
-    AliESDtrack *track = dynamic_cast<AliESDtrack*>(fInputEvent->GetTrack(i));
-    if(!track) continue;
-    Int_t  pid = GetTenderPidDaughter(track);
-    if(pid < 0) continue;
-    Float_t pT = track->Pt();
-    Int_t label = track->GetLabel();
-    if(label < 0) continue;
-    AliMCParticle *mcp = 0x0;
-    mcp = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(label));
-    if(!mcp) continue;
-    Int_t pdg = TMath::Abs(mcp->PdgCode());
-    if(pdg == pid2pdg[pid]) id = kTRUE;
-    if(id) name = "h_" + type[pid] + "_pt_S";
-    else name = "h_" + type[pid] + "_pt_B";
-    fCollMC->Fill(name, pT);
-  }
-}
-//__________________________________________________________
-Int_t AliAnalysisTaskCheckV0tender::GetTenderPidV0(AliESDv0 * const v0){
-  //
-  // retrieve the PID nformation stored in the status flags by the train tender
-  // 
-  Int_t pid = -1;
-  if(!v0){
-    return pid;
-  }
-  Int_t nTimes = 0;
-  if(v0->TestBit(BIT(14))){
-    pid = 0;
-    nTimes++;
-  }
-  if(v0->TestBit(BIT(15))){
-    pid = 1;
-    nTimes++;
-  }
-  if(v0->TestBit(BIT(16))){
-    pid = 2;
-    nTimes++;
-  }
-  if(v0->TestBit(BIT(17))){
-    pid = 3;
-    nTimes++;
-  }
-  if(nTimes > 1){
-    AliWarning("V0 track labeled multiple times by the V0 tender");
-    pid = -1;
-  }    
-
-  //printf(" -D: pid: %i \n", pid);
-
-  return pid;
-}
-//__________________________________________________________
-Int_t AliAnalysisTaskCheckV0tender::GetTenderPidDaughter(AliESDtrack * const track){
-  //
-  // retrieve the PID nformation stored in the status flags by the train tender
-  // 
-
-  Int_t pid = -1;
-  if(!track){
-    return pid;
-  }
-  Int_t nTimes = 0;
-  if(track->TestBit(BIT(14))){
-    pid = 0;
-    nTimes++;
-  }
-  if(track->TestBit(BIT(15))){
-    pid = 1;
-    nTimes++;
-  }
-  if(track->TestBit(BIT(16))){
-    pid = 2;
-    nTimes++;
-  }
-  if(nTimes > 1){
-    AliWarning("V0 track labeled multiple times by the V0 tender");
-    pid = -1;
-  }    
-  return pid;
-}
-//__________________________________________________________
-Float_t AliAnalysisTaskCheckV0tender::MassV0(AliESDv0 * const v0, Int_t id){
-  //
-  // Get the V0 effective mass
-  //
-
-  Float_t mass = -0.1;
-  Bool_t sign = CheckSigns(v0);
-  if(0 == id){
-    mass = v0->GetEffMass(0, 0);
-  }
-  else if(1 == id){
-    mass = v0->GetEffMass(2, 2);
-  }
-  else if(2 == id){
-    mass = (sign) ? v0->GetEffMass(4, 2) : v0->GetEffMass(2, 4);
-  }
-  else if(3 == id){
-    mass = (sign) ? v0->GetEffMass(2, 4) : v0->GetEffMass(4, 2);
-  }
-  else{
-    AliWarning(Form("Unrecognized V0 id: %i", id));
-  }
-
-  return mass;
-
-}
-//__________________________________________________________
-Bool_t AliAnalysisTaskCheckV0tender::CheckSigns(AliESDv0 * const v0){
-  //
-  // check wheter the sign was correctly applied to 
-  // V0 daughter tracks
-  // This function should become obsolete once the V0 finder will be updated (fixed)
-  //
-  
-  Bool_t correct = kFALSE;
-
-  Int_t pIndex = 0, nIndex = 0;
-  pIndex = v0->GetPindex();
-  nIndex = v0->GetNindex();
-  
-  AliESDtrack* d[2];
-  d[0] = dynamic_cast<AliESDtrack*>(fInputEvent->GetTrack(pIndex));
-  d[1] = dynamic_cast<AliESDtrack*>(fInputEvent->GetTrack(nIndex));
-
-  Int_t sign[2];
-  sign[0] = (int)d[0]->GetSign();
-  sign[1] = (int)d[1]->GetSign();
-  
-  if(-1 == sign[0] && 1 == sign[1]){
-    correct = kFALSE;
-    //v0->SetIndex(0, pIndex);  // set the index of the negative v0 track
-    //v0->SetIndex(1, nIndex);  // set the index of the positive v0 track
-  }
-  else{
-    correct = kTRUE;
-  }
-  
-  //pIndex = v0->GetPindex();
-  //nIndex = v0->GetNindex();
-  //printf("-D2: P: %i, N: %i\n", pIndex, nIndex);
-
-  return correct;
-}
diff --git a/PWG3/hfe/AliAnalysisTaskCheckV0tender.h b/PWG3/hfe/AliAnalysisTaskCheckV0tender.h
deleted file mode 100644 (file)
index dd015f7..0000000
+++ /dev/null
@@ -1,81 +0,0 @@
-#ifndef ALIANALYSISTASKCHECKV0TENDER_H
-#define ALIANALYSISTASKCHECKV0TENDER_H
-
-/**************************************************************************
-* Copyright(c) 1998-1999, 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$ */ 
-
-//
-// Task for PID QA
-// Using AliHFEpidQA and AliHFEMCpidQA
-// More information can be found in the source file
-//
-#ifndef ALIANALYSISTASKSE_H
-#include "AliAnalysisTaskSE.h"
-#endif
-
-class TH1F;
-class TList;
-class AliHFEcollection;
-
-class AliAnalysisTaskCheckV0tender : public AliAnalysisTaskSE{
- public:
-   enum{ // Reconstructed V0
-      kRecoGamma = 0,
-      kRecoK0 = 1,
-      kRecoLambda = 2,
-      kRecoALambda = 3
-      
-      };
-  enum{ // Identified Daughter particles
-    kRecoElectron = 0,
-      kRecoPionK0 = 1,
-      kRecoPionL = 2,
-      kRecoProton = 3
-      };
-
-  AliAnalysisTaskCheckV0tender();
-  AliAnalysisTaskCheckV0tender(const Char_t *name);
-  ~AliAnalysisTaskCheckV0tender();
-  
-  virtual void UserCreateOutputObjects();
-  virtual void UserExec(Option_t *);
-  virtual void Terminate(Option_t *);
-  
-  void ProcessV0s();
-  void ProcessDaughters();
-  void ProcessV0sMC();
-  void ProcessDaughtersMC();
-
-  Int_t GetTenderPidV0(AliESDv0 * const v0);
-  Int_t GetTenderPidDaughter(AliESDtrack * const track);
-
- private:
-  AliAnalysisTaskCheckV0tender(const AliAnalysisTaskCheckV0tender &ref);
-  AliAnalysisTaskCheckV0tender &operator=(const AliAnalysisTaskCheckV0tender &ref);
-
-  Float_t MassV0(AliESDv0 * const v0, Int_t id);
-  Bool_t CheckSigns(AliESDv0 * const v0);
-  TList *fOutput;            //! Container for output histos
-  AliHFEcollection *fColl;   //! collection of Data output
-  AliHFEcollection *fCollMC; //! collection of MC output
-  TH1 *fEvents;              //! Number of Events
-  ClassDef(AliAnalysisTaskCheckV0tender, 1)
-
-};
-
-#endif
diff --git a/PWG3/hfe/AliAnalysisTaskCheckV0tenderII.cxx b/PWG3/hfe/AliAnalysisTaskCheckV0tenderII.cxx
deleted file mode 100644 (file)
index 9ce543d..0000000
+++ /dev/null
@@ -1,603 +0,0 @@
-/**************************************************************************
-* Copyright(c) 1998-1999, 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$ */
-
-//
-// Task for checking the performance of the V0 tender
-// 
-// 
-// Authors:
-//   Matus Kalisky <matus.kalisky@cern.ch>
-//
-
-#include <TH1F.h>
-#include <TList.h>
-
-#include "AliAnalysisManager.h"
-#include "AliMCEventHandler.h"
-#include "AliESDInputHandler.h"
-#include "AliESDv0.h"
-#include "AliESDtrack.h"
-#include "AliMCParticle.h"
-#include "AliMCEvent.h"
-#include "AliESDv0KineCuts.h"
-#include "AliKFVertex.h"
-
-#include "AliHFEtools.h"
-#include "AliHFEcollection.h"
-
-#include "AliAnalysisTaskCheckV0tenderII.h"
-
-ClassImp(AliAnalysisTaskCheckV0tenderII)
-
-//__________________________________________________________
-AliAnalysisTaskCheckV0tenderII::AliAnalysisTaskCheckV0tenderII():
-  AliAnalysisTaskSE("CheckV0tenderTask")
-  , fOutput(0x0)
-  , fColl(0x0)
-  , fCollMC(0x0)
-  , fV0cuts(0x0)
-  , fPrimaryVertex(0x0)
-  , fpdgV0(-1)
-  , fpdgP(-1)
-  , fpdgN(-1)
-  , fEvents(0x0)
-{
-  //
-  // Default Constructor
-  //
-
-
-  DefineOutput(1, TH1F::Class());
-  DefineOutput(2, TList::Class());
-
-
-}
-//__________________________________________________________
-AliAnalysisTaskCheckV0tenderII::AliAnalysisTaskCheckV0tenderII(const Char_t *name):
-  AliAnalysisTaskSE(name)
-  , fOutput(0x0)
-  , fColl(0x0)
-  , fCollMC(0x0)
-  , fV0cuts(0x0)
-  , fPrimaryVertex(0x0)
-  , fpdgV0(0)
-  , fpdgP(0)
-  , fpdgN(0)
-  , fEvents(0x0)
-{
-  //
-  // Constructor
-  //
-
-  DefineOutput(1, TH1F::Class());
-  DefineOutput(2, TList::Class());
-
-}
-//__________________________________________________________
-AliAnalysisTaskCheckV0tenderII::~AliAnalysisTaskCheckV0tenderII(){
-  //
-  // Destructor
-  //
-
-  if (fOutput) delete fOutput;
-  //if (fColl) delete fColl;
-  //if (fCollMC) delete fCollMC;
-  if (fV0cuts) delete fV0cuts;
-  if (fPrimaryVertex) delete fPrimaryVertex;
-  if (fEvents) delete fEvents;
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tenderII::UserCreateOutputObjects(){
-  //
-  // prepare output objects
-  //
-
-  fOutput = new TList();
-  //fOutput->SetOwner();
-  // Counter for number of events
-  fEvents = new TH1I("nEvents", "NumberOfEvents", 1, 1, 2);
-
-  fColl = new AliHFEcollection("V0_QA", "tender V0s for data");
-  fCollMC = new AliHFEcollection("V0_MC_QA", "tender V0s for MC");
-
-  fV0cuts = new AliESDv0KineCuts();
-  if(!fV0cuts){
-    AliError("Failed to initialize AliESDv0KineCuts instance");
-    return;
-  }
-
-  // 
-  // Data histos
-  //
-  
-  // total number of tagged V0s
-  fColl->CreateTH1F("h_NumberOf_V0s", "Number of tagged V0s; type; counts", 4, -0.5, 3.5);
-  // pT spectra of the tagged V0s
-  fColl->CreateTH1F("h_Gamma_pt", "p_{T} spectrum of tagged gammas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_K0_pt", "p_{T} spectrum of tagged K0s; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_Lambda_pt", "p_{T} spectrum of tagged Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_ALambda_pt", "p_{T} spectrum of tagged A-Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  // invariant mass of the V0s
-  fColl->CreateTH1F("h_Gamma_mass", "Inv. Mass of the ", 100, 0., 0.1);
-  fColl->CreateTH1F("h_K0_mass", "Inv. Mass of the ", 100, 0.45, 0.55);
-  fColl->CreateTH1F("h_Lambda_mass", "Inv. Mass of the ", 100, 1.08, 1.14);
-  fColl->CreateTH1F("h_ALambda_mass", "Inv. Mass of the ", 100, 1.08, 1.14);
-
-  // total number of tagged daughter particles (should correlate with number of V0s !
-  fColl->CreateTH1F("h_NumberOfDaughters", "Number of tagged daughters; type; counts", 3, -0.5, 2.5);
-  // pT spectra of tagged daughter particles
-  fColl->CreateTH1F("h_Electron_pt", "p_{T} spectrum of tagged; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_Pion_pt", "p_{T} spectrum of tagged; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fColl->CreateTH1F("h_Proton_pt", "p_{T} spectrum of tagged; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-
-  //
-  // MC histos
-  //
-
-  // pT spectra of the tagged V0s
-  fCollMC->CreateTH1F("h_Gamma_pt_S", "MC-S: p_{T} spectrum of tagged gammas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_K0_pt_S", "MC-S: p_{T} spectrum of tagged K0s; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Lambda_pt_S", "MC-S: p_{T} spectrum of tagged Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_ALambda_pt_S", "MC-S: p_{T} spectrum of tagged A-Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0); 
-  fCollMC->CreateTH1F("h_Gamma_pt_B", "MC-B: p_{T} spectrum of tagged gammas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_K0_pt_B", "MC-B: p_{T} spectrum of tagged K0s; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Lambda_pt_B", "MC-B: p_{T} spectrum of tagged Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_ALambda_pt_B", "MC-B: p_{T} spectrum of tagged A-Lambdas; p_{T} (GeV/c); counts", 20, 0.1, 20, 0); 
-  // invariant mass of the V0s
-  fCollMC->CreateTH1F("h_Gamma_mass_S", "MC-S: Inv. Mass of the gamma; m (GeV/c^{2}); counts", 100, 0., 0.1);
-  fCollMC->CreateTH1F("h_K0_mass_S", "MC-S: Inv. Mass of the K0; m (GeV/c^{2}); counts", 100, 0.45, 0.55);
-  fCollMC->CreateTH1F("h_Lambda_mass_S", "MC-S: Inv. Mass of the Lambda; m (GeV/c^{2}); counts", 100, 1.08, 1.14);
-  fCollMC->CreateTH1F("h_ALambda_mass_S", "MC-S: Inv. Mass of the A-Lambda; m (GeV/c^{2}); counts", 100, 1.08, 1.14);
-  fCollMC->CreateTH1F("h_Gamma_mass_B", "MC-B: Inv. Mass of the gamma; m (GeV/c^{2}); counts", 100, 0., 0.1);
-  fCollMC->CreateTH1F("h_K0_mass_B", "MC-B: Inv. Mass of the K0; m (GeV/c^{2}); counts", 100, 0.45, 0.55);
-  fCollMC->CreateTH1F("h_Lambda_mass_B", "MC-B: Inv. Mass of the Lambda; m (GeV/c^{2}); counts", 100, 1.08, 1.14);
-  fCollMC->CreateTH1F("h_ALambda_mass_B", "MC-B: Inv. Mass of the A-Lambda; m (GeV/c^{2}); counts", 100, 1.08, 1.14);
-  // pT spectra of tagged daughter particles
-  fCollMC->CreateTH1F("h_Electron_pt_S", "MC-S: p_{T} spectrum of tagged electrons; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Pion_pt_S", "MC-S: p_{T} spectrum of tagged pions; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Proton_pt_S", "MC-S: p_{T} spectrum of tagged protons; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Electron_pt_B", "MC-B: p_{T} spectrum of tagged electrons; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Pion_pt_B", "MC-B: p_{T} spectrum of tagged pions; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1F("h_Proton_pt_B", "MC-B: p_{T} spectrum of tagged protons; p_{T} (GeV/c); counts", 20, 0.1, 20, 0);
-  // background study - MC only
-  fCollMC->CreateTH1Fvector1(6, "h_Gamma_Bg", "MC gamma bg. - different sources; p_{T} (GeV/c), counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1Fvector1(6, "h_K0_Bg", "MC K0 bg. - different sources; p_{T} (GeV/c), counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1Fvector1(6, "h_Lambda_Bg", "MC gamma bg. - different sources; p_{T} (GeV/c), counts", 20, 0.1, 20, 0);
-  fCollMC->CreateTH1Fvector1(6, "h_ALambda_Bg", "MC gamma bg. - different sources; p_{T} (GeV/c), counts", 20, 0.1, 20, 0);
-  
-
-  TList *tmp = fColl->GetList();
-  tmp->SetName(fColl->GetName());
-  fOutput->Add(tmp);
-  tmp = 0x0;
-  tmp = fCollMC->GetList();
-  tmp->SetName(fCollMC->GetName());
-  fOutput->Add(tmp);
-  
-  
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tenderII::UserExec(Option_t *){
-  //
-  // Event Loop
-  // 
-  AliMCEventHandler* mcHandler = (dynamic_cast<AliMCEventHandler*>(AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler()));
-
-  //AliESDInputHandler *inh = dynamic_cast<AliESDInputHandler *>(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
-  //AliESDpid *workingPID = NULL;
-  //if(inh && (workingPID = inh->GetESDpid())) workingPID = inh->GetESDpid();
-  //else workingPID = AliHFEtools::GetDefaultPID(mcHandler ? kTRUE : kFALSE);
-     
-  // check the MC data
-  if(fMCEvent && !mcHandler ) return;
-  if(fMCEvent &&  !mcHandler->InitOk() ) return;
-  if(fMCEvent &&  !mcHandler->TreeK() ) return;
-  if(fMCEvent &&  !mcHandler->TreeTR() ) return;
-
-  AliKFVertex *primaryVertex = new AliKFVertex(*(fInputEvent->GetPrimaryVertex()));
-  if(!primaryVertex) return;
-  
-  fV0cuts->SetPrimaryVertex(primaryVertex);
-  fV0cuts->SetEvent(fInputEvent);
-  
-  ProcessV0s();
-  
-
-  if (primaryVertex) delete primaryVertex;
-  fEvents->Fill(1.1);
-  PostData(1, fEvents);
-  PostData(2, fOutput);
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tenderII::Terminate(Option_t *){
-  //
-  // Do Post Processing
-  //
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tenderII::ProcessV0s(){
-  //
-  // loop over the V0s and extract the information about
-  // the V0s tagged by the V0 tender
-  //
-
-  const TString type[4] = {"Gamma", "K0", "Lambda", "ALambda"};
-  TString name;
-
-  Int_t nV0s = fInputEvent->GetNumberOfV0s();
-  for(Int_t i=0; i<nV0s; ++i){
-    AliESDv0 *esdV0 = (static_cast<AliESDEvent *>(fInputEvent))->GetV0(i);
-    if(!esdV0) continue;
-    if(!esdV0->GetOnFlyStatus()) continue; // Take only V0s from the On-the-fly v0 finder
-
-    ResetPDGcodes();
-
-    // New standalone V0 selection software
-    if( ! (fV0cuts->ProcessV0(esdV0, fpdgV0, fpdgP, fpdgN)) ) continue;
-    
-    Int_t pid = PDGtoPIDv0(fpdgV0);
-    //printf(" -D: pdg: %i, pid: %i\n", fpdgV0, pid);
-    if(pid <= -1) continue;
-    
-    fColl->Fill("h_NumberOf_V0s", pid);
-    name = "h_" + type[pid] + "_pt";
-    Float_t pT = esdV0->Pt();
-    fColl->Fill(name, pT);
-    Float_t mass = MassV0(esdV0, pid);
-    name = "h_" + type[pid] + "_mass";
-    fColl->Fill(name, mass);
-
-    ProcessDaughters(esdV0);
-    if(fMCEvent){
-      ProcessV0sMC(esdV0);
-      ProcessDaughtersMC(esdV0);
-      ProcessBackground(esdV0);
-    }
-
-  }
-
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tenderII::ProcessDaughters(AliESDv0 * const v0){
-  //
-  // produce some check plots for V0 tender tagged single tracks
-  //
-
-  const TString type[3] = {"Electron", "Pion", "Proton"};
-  TString name;
-
-  if(!v0) return;
-
-  // daughter tracks
-  const Int_t nTracks = fInputEvent->GetNumberOfTracks();
-  AliESDtrack *d[2];
-  Int_t iN, iP;
-  iN = iP = -1;
-  iP = v0->GetPindex();
-  iN = v0->GetNindex();
-  if(iN < 0 || iP < 0) return;
-  if(iN >= nTracks || iP >= nTracks) return;
-  d[0] = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(iP));
-  d[1] = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(iN));  
-  if(!d[0] || !d[1]) return;
-
-  for(Int_t i=0; i<2; ++i){
-    Int_t  pid = (i == 0) ? PDGtoPID(fpdgP) : PDGtoPID(fpdgN);
-    if(pid < 0) continue;
-    fColl->Fill("h_NumberOfDaughters", pid*1.0);
-    Float_t pT = d[i]->Pt();
-    name = "h_" + type[pid] + "_pt";
-    fColl->Fill(name, pT);
-  }
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tenderII::ProcessV0sMC(AliESDv0 * const v0){
-  //
-  // check all V0tender selected V0 on their true identity
-  // 
-
-  const TString type[4] = {"Gamma", "K0", "Lambda", "ALambda"};
-  TString name;
-  
-  Int_t pid = PDGtoPIDv0(fpdgV0);
-  if(pid < 0) return;
-
-  // true if fpdgV0 agrees with MC pdg of the mother
-  Bool_t id = kFALSE;
-  const Int_t nTracks = fInputEvent->GetNumberOfTracks();
-  // both ESD daughtr tracks
-  Int_t iN, iP;
-  iN = iP = -1;
-  iP = v0->GetPindex();
-  iN = v0->GetNindex();    
-  if(iN < 0 || iP < 0) return;
-  if(iN >= nTracks || iP >= nTracks) return;    
-  AliESDtrack *dP, *dN;
-  dP = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(iP));
-  dN = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(iN));  
-  if(!dN || !dP) return;
-
-  // MC labels of the daughter tracks
-  Int_t lN, lP;
-  lN = dN->GetLabel();
-  lP = dP->GetLabel();
-  if(lN < 0 || lP < 0) return;
-
-  // MC daughter particles
-  AliMCParticle *mcP, *mcN;
-  mcP = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(lP));
-  mcN = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(lN));
-  if(!mcP || !mcN) return;
-
-  // labels of the mother particles
-  Int_t lPm, lNm;
-  lPm = mcP->GetMother();
-  lNm = mcN->GetMother();
-  AliMCParticle *m = 0x0;
-  // particles without mother particle are most probably
-  // primary particles
-  if(lPm >= 0){
-    m = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(lPm)); 
-    if(!m) return;
-  }
-  if(m){
-    Int_t pdg = m->PdgCode();
-    //if(lPm == lNm) printf(" -D: pdg: %i, fpdgV0: %i \n", pdg, fpdgV0);
-    if((lPm == lNm) && (pdg == fpdgV0)) id = kTRUE;
-  }
-  
-  if(id) name = "h_" + type[pid] + "_pt_S";
-  else  name = "h_" + type[pid] + "_pt_B";
-  Float_t pT = v0->Pt();
-  fCollMC->Fill(name, pT);
-
-  if(id) name = "h_" + type[pid] + "_mass_S";
-  else name = "h_" + type[pid] + "_mass_B";
-  Float_t mass = MassV0(v0, pid);
-  fCollMC->Fill(name, mass);
-
-  
-}
-//__________________________________________________________
-void AliAnalysisTaskCheckV0tenderII::ProcessDaughtersMC(AliESDv0 * const v0){
-  //
-  // check the identity of the V0tender selected V0 daughters
-  // !!! for positive check only the true identity plays a role here, 
-  // not the true V0 mother identity (e.g. selected electron could come
-  // from primary vertex or pion dalitz deca or true gamma conversion) !!!!
-  //
-
-  const TString type[3] = {"Electron", "Pion", "Proton"};
-  TString name;
-
-  if(!v0) return;
-
-  // daughter tracks
-  const Int_t nTracks = fInputEvent->GetNumberOfTracks();
-  AliESDtrack *d[2];
-  Int_t iN, iP;
-  iN = iP = -1;
-  iP = v0->GetPindex();
-  iN = v0->GetNindex();
-  if(iN < 0 || iP < 0) return;
-  if(iN >= nTracks || iP >= nTracks) return;
-  d[0] = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(iP));
-  d[1] = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(iN));  
-  if(!d[0] || !d[1]) return;
-
-  //printf(" *** fpdgV0: %i, fpdgP: %i, fpdgN: %i \n", fpdgV0, fpdgP, fpdgN);
-
-  for(Int_t i=0; i<2; ++i){
-    Bool_t id = kFALSE;
-    Int_t  pid = (i == 0) ? PDGtoPID(fpdgP) : PDGtoPID(fpdgN);
-    Int_t  pdg  = (i == 0) ? fpdgP : fpdgN;
-    if(pid < 0) continue;
-    Float_t pT = d[i]->Pt();
-    Int_t label = d[i]->GetLabel();
-    if(label < 0) continue;
-    AliMCParticle *mcp = 0x0;
-    mcp = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(label));
-    if(!mcp) continue;
-    Int_t pdgMC = mcp->PdgCode();
-    //printf("   * pid: %i, pdg: %i, pdgMC: %i \n", pid, pdg, pdgMC);
-    if(pdgMC == pdg) id = kTRUE;
-    if(id) name =  "h_" + type[pid] + "_pt_S";
-    else name =  "h_" + type[pid] + "_pt_B";
-    fCollMC->Fill(name, pT);
-  }
-}
-//__________________________________________________________
-void  AliAnalysisTaskCheckV0tenderII::ProcessBackground(AliESDv0 * const v0){
-  //
-  // look at the compostition and the properties of the
-  // miss-identified V0
-  //
-
-  if(!v0) return;
-  Float_t pt = v0->Pt();
-  const TString type[4] = {"Gamma", "K0", "Lambda", "ALambda"};
-  TString name;
-  // daughter tracks
-  const Int_t nTracks = fInputEvent->GetNumberOfTracks();
-  AliESDtrack *d[2];
-  Int_t index[2] = {-1, -1};
-  index[0] = v0->GetPindex();
-  index[1] = v0->GetNindex();
-  if(index[0] < 0 || index[1] < 0) return;
-  if(index[0] >= nTracks || index[1] >= nTracks) return;
-  d[0] = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(index[0]));
-  d[1] = dynamic_cast<AliESDtrack *>(fInputEvent->GetTrack(index[1]));  
-  if(!d[0] || !d[1]) return;
-
-  // daughter MC particles
-  Int_t label[2] = {d[0]->GetLabel(), d[1]->GetLabel()};
-  if(label[0] < 0 || label[1] < 0) return;
-  AliMCParticle *dmc[2] = {0x0, 0x0};
-  dmc[0] = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(label[0]));
-  dmc[1] = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(label[1]));
-  if(!dmc[0]  || !dmc[1]) return;
-
-  Bool_t primary[2] = {dmc[0]->Particle()->IsPrimary(), dmc[1]->Particle()->IsPrimary()};
-
-  // mother MC particles
-  Int_t labelM[2] = {dmc[0]->GetMother(), dmc[1]->GetMother()};
-  AliMCParticle *mmc[2] = {0x0, 0x0};
-  if(!primary[0])  mmc[0] = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(labelM[0]));
-  if(!primary[1])  mmc[1] = dynamic_cast<AliMCParticle*>(fMCEvent->GetTrack(labelM[1]));
-
-  // if particle is not primary it MUST have a mother particle
-  if(!primary[0] && !mmc[0]) return;
-  if(!primary[1] && !mmc[1]) return;
-
-  Int_t pdgM[2] = {0, 0};
-  if(mmc[0]) pdgM[0] = mmc[0]->PdgCode();
-  if(mmc[1]) pdgM[1] = mmc[1]->PdgCode();
-
-  // select only miss-identified V0s
-  if(labelM[0] == labelM[1] && pdgM[0] == fpdgV0) return;
-
-  // indices for V0 background histograms
-  // 0 - random match
-  // 1 - gamma
-  // 2 - K0
-  // 3 - lambda 
-  // 4 - a-lambda
-  // 5 - other dacay or interaction
-
-  Int_t ixM = -1;  
-
-  Int_t typeV0  = (labelM[0] != labelM[1]) ? 0 : 1;
-  if(0 == typeV0) ixM = 0;
-  else{
-    ixM = PDGtoPIDv0(pdgM[0]) + 1;  
-    if(ixM < 0) ixM = 5; 
-  }
-  Int_t ix = PDGtoPIDv0(fpdgV0);
-  if(0 <= ix ){
-    name = "h_" + type[PDGtoPIDv0(fpdgV0)] + "_Bg"; 
-    fCollMC->Fill(name, ixM, pt);
-  }
-  
-  // now look at the daughter tracks
-  
-
-}
-//__________________________________________________________
-Float_t AliAnalysisTaskCheckV0tenderII::MassV0(AliESDv0 * const v0, Int_t id){
-  //
-  // Get the V0 effective mass
-  //
-
-  Float_t mass = -0.1;
-  Bool_t sign = CheckSigns(v0);
-  if(0 == id){
-    mass = v0->GetEffMass(0, 0);
-  }
-  else if(1 == id){
-    mass = v0->GetEffMass(2, 2);
-  }
-  else if(2 == id){
-    mass = (sign) ? v0->GetEffMass(4, 2) : v0->GetEffMass(2, 4);
-  }
-  else if(3 == id){
-    mass = (sign) ? v0->GetEffMass(2, 4) : v0->GetEffMass(4, 2);
-  }
-  else{
-    AliWarning(Form("Unrecognized V0 id: %i", id));
-  }
-
-  return mass;
-
-}
-//__________________________________________________________
-Bool_t AliAnalysisTaskCheckV0tenderII::CheckSigns(AliESDv0 * const v0){
-  //
-  // check wheter the sign was correctly applied to 
-  // V0 daughter tracks
-  // This function should become obsolete once the V0 finder will be updated (fixed)
-  //
-  
-  Bool_t correct = kFALSE;
-
-  Int_t pIndex = 0, nIndex = 0;
-  pIndex = v0->GetPindex();
-  nIndex = v0->GetNindex();
-  
-  AliESDtrack* d[2];
-  d[0] = dynamic_cast<AliESDtrack*>(fInputEvent->GetTrack(pIndex));
-  d[1] = dynamic_cast<AliESDtrack*>(fInputEvent->GetTrack(nIndex));
-
-  Int_t sign[2];
-  sign[0] = (int)d[0]->GetSign();
-  sign[1] = (int)d[1]->GetSign();
-  
-  if(-1 == sign[0] && 1 == sign[1]){
-    correct = kFALSE;
-    //v0->SetIndex(0, pIndex);  // set the index of the negative v0 track
-    //v0->SetIndex(1, nIndex);  // set the index of the positive v0 track
-  }
-  else{
-    correct = kTRUE;
-  }
-  
-  //pIndex = v0->GetPindex();
-  //nIndex = v0->GetNindex();
-  //printf("-D2: P: %i, N: %i\n", pIndex, nIndex);
-
-  return correct;
-}
-//__________________________________________________________
-Int_t  AliAnalysisTaskCheckV0tenderII::PDGtoPIDv0(Int_t pdgV0) const {
-  //
-  // convert thereconstructed V0 pdg to local pid
-  //
-
-  switch(pdgV0){
-  case 22: return 0; break;
-  case 310: return 1; break;
-  case 3122: return 2; break;
-  case -3122: return 3; break;
-  }
-  
-  return -1;
-
-}
-//__________________________________________________________
-Int_t  AliAnalysisTaskCheckV0tenderII::PDGtoPID(Int_t pdg) const {
-  //
-  // convert daughter pdg code to local pid
-  //
-  switch(TMath::Abs(pdg)){
-  case 11: return 0; break;
-  case 211: return 1; break;
-  case 2212: return 2; break;
-  }
-  return -1;
-
-}
-//__________________________________________________________
-void  AliAnalysisTaskCheckV0tenderII::ResetPDGcodes(){
-  //
-  // reset the PDG codes of the V0 and daughter
-  // particles to default values
-  //
-
-  fpdgV0 = -1;
-  fpdgP = -1;
-  fpdgN = -1;
-
-}
diff --git a/PWG3/hfe/AliAnalysisTaskCheckV0tenderII.h b/PWG3/hfe/AliAnalysisTaskCheckV0tenderII.h
deleted file mode 100644 (file)
index ee94bda..0000000
+++ /dev/null
@@ -1,95 +0,0 @@
-#ifndef ALIANALYSISTASKCHECKV0TENDERII_H
-#define ALIANALYSISTASKCHECKV0TENDERII_H
-
-/**************************************************************************
-* Copyright(c) 1998-1999, 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$ */ 
-
-//
-// Task for perfomance studies of V0 selection code
-// More information can be found in the source file
-//
-#ifndef ALIANALYSISTASKSE_H
-#include "AliAnalysisTaskSE.h"
-#endif
-
-class TH1F;
-class TList;
-class AliHFEcollection;
-class AliESDv0KineCuts;
-class AliKFVertex;
-
-class AliAnalysisTaskCheckV0tenderII : public AliAnalysisTaskSE{
- public:
-   enum{ // Reconstructed V0
-      kRecoGamma = 0,
-      kRecoK0 = 1,
-      kRecoLambda = 2,
-      kRecoALambda = 3
-      
-      };
-  enum{ // Identified Daughter particles
-    kRecoElectron = 0,
-      kRecoPionK0 = 1,
-      kRecoPionL = 2,
-      kRecoProton = 3
-      };
-
-  AliAnalysisTaskCheckV0tenderII();
-  AliAnalysisTaskCheckV0tenderII(const Char_t *name);
-  ~AliAnalysisTaskCheckV0tenderII();
-  
-  virtual void UserCreateOutputObjects();
-  virtual void UserExec(Option_t *);
-  virtual void Terminate(Option_t *);
-  
-  void ProcessV0s();
-  void ProcessDaughters(AliESDv0 * const v0);
-  void ProcessV0sMC(AliESDv0 * const v0);
-  void ProcessDaughtersMC(AliESDv0 * const v0);
-  void ProcessBackground(AliESDv0 * const v0);
-
- private:
-  AliAnalysisTaskCheckV0tenderII(const AliAnalysisTaskCheckV0tenderII &ref);
-  AliAnalysisTaskCheckV0tenderII &operator=(const AliAnalysisTaskCheckV0tenderII &ref);
-
-  Float_t MassV0(AliESDv0 * const v0, Int_t id);
-  Bool_t  CheckSigns(AliESDv0 * const v0);
-
-  Int_t   PDGtoPIDv0(Int_t pdgV0) const;
-  Int_t   PDGtoPID(Int_t pdg) const;  
-  
-  void    ResetPDGcodes();
-  
-  TList              *fOutput;        //! Container for output histos
-  AliHFEcollection   *fColl;          //! collection of Data output
-  AliHFEcollection   *fCollMC;        //! collection of MC output
-  AliESDv0KineCuts   *fV0cuts;        //! standalone V0 selection class
-  AliKFVertex        *fPrimaryVertex; //! primary vertex of the current event
-
-  Int_t               fpdgV0;  // PDG code of the reconstructed V0
-  Int_t               fpdgP;   // PDG code of the positive daughter (sign corrected)
-  Int_t               fpdgN;   // PDG code of the negative daughter (sign coreccted)
-  
-  TH1 *fEvents;              //! Number of Events
-
-  ClassDef(AliAnalysisTaskCheckV0tenderII, 1)
-
-};
-
-
-#endif
diff --git a/PWG3/hfe/AliESDv0KineCuts.cxx b/PWG3/hfe/AliESDv0KineCuts.cxx
deleted file mode 100644 (file)
index 977ae7e..0000000
+++ /dev/null
@@ -1,1036 +0,0 @@
-/**************************************************************************
- * Copyright(c) 1998-1999, 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$ */
-
-/*
- * author: M.Kalisky@gsi.de
- * 08/Dec/2010
- *
- * Description: This class allows with purely kinematical cuts
- * to select clean samples of electrons, pions and protons from the
- * V0 online finder ESD V0 candidates for PID and dectector resonse
- * studies.
- */
-
-#include <TVector3.h>
-#include <TDatabasePDG.h>
-
-#include "AliESDv0.h"
-#include "AliESDtrack.h"
-#include "AliESDEvent.h"
-#include "AliVEvent.h"
-#include "AliLog.h"
-#include "AliKFParticle.h"
-#include "AliVTrack.h"
-#include "AliKFVertex.h"
-
-#include "AliESDv0KineCuts.h"
-
-ClassImp(AliESDv0KineCuts)
-
-//____________________________________________________________________
-AliESDv0KineCuts::AliESDv0KineCuts() :
-  fEvent(0x0)
-  , fPrimaryVertex(0x0)
-  , fType(0)
-  , fMode(0)
-  , fTPCNcls(1)
-  , fTPCrefit(kTRUE)
-  , fTPCchi2perCls(4.0)
-  , fTPCclsRatio(0.6)
-  , fNoKinks(kTRUE)
-  , fGcutChi2NDF(10)
-  , fGcutInvMass(0.05)
-  , fK0cutChi2NDF(10)
-  , fLcutChi2NDF(10)
-{
-  //
-  // Default constructor
-  //
-
-  // default single track cuts
-  fTPCNcls = 1;                // minimal number of the TPC clusters
-  fTPCrefit = kTRUE;           // TPC refit
-  fTPCchi2perCls = 4.0;        // chi2 per TPC cluster
-  fTPCclsRatio = 0.6;          // minimal foun/findable TPC cluster ratio
-  fNoKinks = kTRUE;            // kinks - no [kTRUE] or do not care [kFalse]
-
-
-  // default gamma cuts values
-  fGcutChi2NDF = 10;           // Chi2NF cut value for the AliKFparticle gamma
-  fGcutCosPoint[0] = 0;        // cos of the pointing angle [min, max]
-  fGcutCosPoint[1] = 0.02;     // cos of the pointing angle [min, max]
-  fGcutDCA[0] = 0.;            // DCA between the daughter tracks [min, max]
-  fGcutDCA[1] = 0.25;          // DCA between the daughter tracks [min, max]
-  fGcutVertexR[0] = 3.;        // radius of the conversion point [min, max]
-  fGcutVertexR[1] = 90.;       // radius of the conversion point [min, max]
-  fGcutPsiPair[0] = 0.;        // value of the psi pair cut [min, max]
-  fGcutPsiPair[1] = 0.05;      // value of the psi pair cut [min, max]
-  fGcutInvMass = 0.05;         // upper value on the gamma invariant mass
-  // default K0 cuts
-  fK0cutChi2NDF = 10;          // Chi2NF cut value for the AliKFparticle K0
-  fK0cutCosPoint[0] = 0.;      // cos of the pointing angle [min, max]
-  fK0cutCosPoint[1] = 0.02;    // cos of the pointing angle [min, max]
-  fK0cutDCA[0] = 0.;           // DCA between the daughter tracks [min, max]
-  fK0cutDCA[1] = 0.2;          // DCA between the daughter tracks [min, max]
-  fK0cutVertexR[0] = 2.0;      // radius of the decay point [min, max]
-  fK0cutVertexR[1] = 30.0;     // radius of the decay point [min, max]
-  fK0cutInvMass[0] = 0.486;    // invariant mass window
-  fK0cutInvMass[1] = 0.508;    // invariant mass window
-  // Lambda & anti-Lambda cut values
-  fLcutChi2NDF = 10;           // Chi2NF cut value for the AliKFparticle K0
-  fLcutCosPoint[0] = 0.;       // cos of the pointing angle [min, max]
-  fLcutCosPoint[1] = 0.02;     // cos of the pointing angle [min, max]
-  fLcutDCA[0] = 0.;            // DCA between the daughter tracks [min, max]
-  fLcutDCA[1] = 0.2;           // DCA between the daughter tracks [min, max]
-  fLcutVertexR[0] = 2.0;       // radius of the decay point [min, max]
-  fLcutVertexR[1] = 40.0;      // radius of the decay point [min, max]
-  fLcutInvMass[0] = 1.11;      // invariant mass window
-  fLcutInvMass[1] = 1.12;      // invariant mass window
-    
-}
-//____________________________________________________________________
-AliESDv0KineCuts::~AliESDv0KineCuts(){
-  //
-  // Destructor
-  //
-
-
-}
-//____________________________________________________________________
-AliESDv0KineCuts::AliESDv0KineCuts(const AliESDv0KineCuts &ref):
-  TObject(ref)
-  , fEvent(0x0)
-  , fPrimaryVertex(0x0)
-  , fType(0)
-  , fMode(0)
-  , fTPCNcls(1)
-  , fTPCrefit(kTRUE)
-  , fTPCchi2perCls(4.0)
-  , fTPCclsRatio(0.6)
-  , fNoKinks(kTRUE)
-  , fGcutChi2NDF(10)
-  , fGcutInvMass(0.05)
-  , fK0cutChi2NDF(10)
-  , fLcutChi2NDF(10)
-{
-  //
-  // Copy operator
-  //
-
-  ref.Copy(*this);
-}
-//____________________________________________________________________
-AliESDv0KineCuts &AliESDv0KineCuts::operator=(const AliESDv0KineCuts &ref){
-  //
-  // assignment operator
-  //
-  if(this != &ref)
-    ref.Copy(*this);
-  return *this; 
-}
-//____________________________________________________________________
-void AliESDv0KineCuts::Copy(TObject &ref) const {
-  //
-  // Performs the copying of the object
-  //
-
-  TObject::Copy(ref);
-
-  AliESDv0KineCuts &target = dynamic_cast<AliESDv0KineCuts &>(ref);
-
-  // default single track cuts
-  target.fTPCNcls = fTPCNcls;
-  target.fTPCrefit = fTPCrefit;
-  target.fTPCchi2perCls = fTPCchi2perCls;
-  target.fTPCclsRatio = fTPCclsRatio;
-  target.fNoKinks = fNoKinks;
-
-
-  // default gamma cuts values
-  target.fGcutChi2NDF = fGcutChi2NDF;
-  memcpy(target.fGcutCosPoint, fGcutCosPoint, sizeof(Float_t) * 2);
-  memcpy(target.fGcutDCA, fGcutDCA, sizeof(Float_t) * 2); 
-  memcpy(target.fGcutVertexR, fGcutVertexR, sizeof(Float_t) * 2);
-  memcpy(target.fGcutPsiPair, fGcutPsiPair, sizeof(Float_t) * 2);
-  target.fGcutInvMass = fGcutInvMass;
-  // default K0 cuts
-  target.fK0cutChi2NDF = fK0cutChi2NDF;
-  memcpy(target.fK0cutCosPoint, fK0cutCosPoint, sizeof(Float_t) * 2);
-  memcpy(target.fK0cutDCA, fK0cutDCA, sizeof(Float_t) * 2);
-  memcpy(target.fK0cutVertexR, fK0cutVertexR, sizeof(Float_t) * 2);
-  memcpy(target.fK0cutInvMass, fK0cutInvMass, sizeof(Float_t) * 2);
-  // Lambda & anti-Lambda cut values
-  target.fLcutChi2NDF = fLcutChi2NDF;
-  memcpy(target.fLcutCosPoint, fLcutCosPoint, sizeof(Float_t) * 2);
-  memcpy(target.fLcutDCA, fLcutDCA, sizeof(Float_t) * 2);
-  memcpy(target.fLcutVertexR, fLcutVertexR, sizeof(Float_t) * 2);
-  memcpy(target.fLcutInvMass, fLcutInvMass, sizeof(Float_t) * 2);
-  
-}
-//____________________________________________________________________
-Bool_t AliESDv0KineCuts::ProcessV0(AliESDv0* const v0, Int_t &pdgV0, Int_t &pdgP, Int_t &pdgN){
-  //
-  // main user function
-  //
-
-  if(!v0) return kFALSE;
-  if(!fEvent){
-    AliErrorClass("No valid Event pointer available, provide it first");
-    return kFALSE;
-  }
-
-  if(!V0CutsCommon(v0)) return kFALSE;
-
-  const Int_t id = PreselectV0(v0);
-
-  if(!SingleTrackCuts(v0)) return kFALSE;
-
-  switch(id){
-  case kUndef:
-    return kFALSE;
-  case kGamma:
-    return CaseGamma(v0, pdgV0, pdgP, pdgN);
-  case kK0:
-    return CaseK0(v0, pdgV0, pdgP, pdgN);
-  case kLambda:
-    return CaseLambda(v0, pdgV0, pdgP, pdgN, 0);
-  case kALambda:
-    return CaseLambda(v0, pdgV0, pdgP, pdgN, 1);
-  default:
-    return kFALSE; 
-  }
-
-  return kFALSE;
-}
-//____________________________________________________________________
-Bool_t AliESDv0KineCuts::ProcessV0(AliESDv0* const v0, Int_t &pdgP, Int_t &pdgN){
-  //
-  // main user function, simplified if the V0 identity is not necessary
-  //
-
-  if(!v0) return kFALSE;
-  if(!fEvent){
-    AliErrorClass("No valid Event pointer available, provide it first");
-    return kFALSE;
-  }
-
-  Int_t idV0 = -1;
-  return ProcessV0(v0, idV0, pdgP, pdgN);
-
-}
-//____________________________________________________________________
-Int_t  AliESDv0KineCuts::PreselectV0(AliESDv0* const v0){
-  //
-  // Make a preselection (exclusive) of the V0 cadidates based on
-  // Armenteros plot
-  // the armenteros cut values are currently fixed and user is not able to set them via
-  // set funcions. The reason is that these cuts are optimized and furneter changes should 
-  // not be necessary. To prove otherwise please study in detail before changing the values
-  //
-  Float_t ap[2] = {-1., -1.};
-  Armenteros(v0, ap);
-  // for clarity
-  const Float_t alpha = ap[0];
-  const Float_t qt = ap[1];
-
-  // selection cuts 
-  // - the reagions for different candidates must not overlap 
-
-  // Gamma cuts
-  const Double_t cutAlphaG = 0.35; 
-  const Double_t cutQTG = 0.05;
-  const Double_t cutAlphaG2[2] = {0.6, 0.8};
-  const Double_t cutQTG2 = 0.04;
-
-  // K0 cuts
-  const Float_t cutQTK0[2] = {0.1075, 0.215};
-  const Float_t cutAPK0[2] = {0.199, 0.8};   // parameters for curved QT cut
-  
-  // Lambda & A-Lambda cuts
-  const Float_t cutQTL = 0.03;
-  const Float_t cutAlphaL[2] = {0.35, 0.7};
-  const Float_t cutAlphaAL[2] = {-0.7,  -0.35};
-  const Float_t cutAPL[3] = {0.107, -0.69, 0.5};  // parameters fir curved QT cut
-
-
-  if(kPurity == fMode){
-  // Check for Gamma candidates
-    if(qt < cutQTG){
-      if( (TMath::Abs(alpha) < cutAlphaG) ) return kGamma;
-    }
-    // additional region - should help high pT gammas
-    if(qt < cutQTG2){
-      if( (TMath::Abs(alpha) > cutAlphaG2[0]) &&  (TMath::Abs(alpha) < cutAlphaG2[1]) ) return kGamma;
-    }
-  }
-  if(kEffGamma == fMode){
-    if(qt < cutQTG) return kGamma;
-  }
-
-  
-  // Check for K0 candidates
-  Float_t q = cutAPK0[0] * TMath::Sqrt(TMath::Abs(1 - alpha*alpha/(cutAPK0[1]*cutAPK0[1])));
-  if( (qt > cutQTK0[0]) && (qt < cutQTK0[1]) && (qt > q) ){
-    return kK0;
-  }
-
-  // Check for Lambda candidates
-  q = cutAPL[0] * TMath::Sqrt(TMath::Abs(1 - ( (alpha + cutAPL[1]) * (alpha + cutAPL[1]) ) / (cutAPL[2]*cutAPL[2]) ));
-  if( (alpha > cutAlphaL[0]) && (alpha < cutAlphaL[1]) && (qt > cutQTL) && (qt < q)  ){
-    return kLambda;
-  }
-
-  // Check for A-Lambda candidates
-  q = cutAPL[0] * TMath::Sqrt(TMath::Abs(1 - ( (alpha - cutAPL[1]) * (alpha - cutAPL[1]) ) / (cutAPL[2]*cutAPL[2]) ));
-  if( (alpha > cutAlphaAL[0]) && (alpha < cutAlphaAL[1]) && (qt > cutQTL) && (qt < q)  ){
-    return kALambda;
-  }
-  
-  return kUndef;
-}
-//____________________________________________________________________
-Bool_t AliESDv0KineCuts::SingleTrackCuts(AliESDv0 * const v0){
-  //
-  // apply single track cuts
-  // correct sign not relevat here
-  //
-
-  if(!v0) return kFALSE;
-  
-  Int_t pIndex = 0, nIndex = 0;
-  pIndex = v0->GetPindex();
-  nIndex = v0->GetNindex();
-  AliESDtrack* d[2];
-  d[0] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(pIndex));
-  d[1] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(nIndex));
-  
-  for(Int_t i=0; i<2; ++i){
-    if(!d[i]) return kFALSE;
-    
-    // status word
-    ULong_t status = d[i]->GetStatus();
-
-    // No. of TPC clusters leave to the users
-    if(d[i]->GetTPCNcls() < 1) return kFALSE;
-
-    // TPC refit
-    if(!(status & AliESDtrack::kTPCrefit)) return kFALSE;
-  
-    // Chi2 per TPC cluster
-    Int_t nTPCclusters = d[i]->GetTPCNcls();
-    Float_t chi2perTPCcluster = d[i]->GetTPCchi2()/Float_t(nTPCclusters);
-    if(chi2perTPCcluster > 4) return kFALSE;
-
-    // TPC cluster ratio
-    Float_t cRatioTPC = d[i]->GetTPCNclsF() > 0. ? static_cast<Float_t>(d[i]->GetTPCNcls())/static_cast<Float_t> (d[i]->GetTPCNclsF()) : 1.;
-    if(cRatioTPC < 0.6) return kFALSE;
-    
-    // kinks
-    if(d[i]->GetKinkIndex(0) != 0) return kFALSE;
-    
-  }
-
-  return kTRUE;
-}
-//____________________________________________________________________
-Bool_t  AliESDv0KineCuts::CaseGamma(AliESDv0* const v0, Int_t &pdgV0, Int_t &pdgP, Int_t &pdgN){
-  //
-  // process the gamma conversion candidate
-  //
-
-  if(!v0) return kFALSE;
-
-  AliVTrack* daughter[2];
-  Int_t pIndex = 0, nIndex = 0;
-
-  Bool_t sign = CheckSigns(v0);
-  if(sign){
-    pIndex = v0->GetPindex();
-    nIndex = v0->GetNindex();
-  }
-  else{
-    pIndex = v0->GetNindex();
-    nIndex = v0->GetPindex();    
-  }
-  daughter[0] = dynamic_cast<AliVTrack *>(fEvent->GetTrack(pIndex));
-  daughter[1] = dynamic_cast<AliVTrack *>(fEvent->GetTrack(nIndex));
-  if(!daughter[0] || !daughter[1]) return kFALSE;
-
-  AliKFParticle *kfMother = CreateMotherParticle(daughter[0], daughter[1], TMath::Abs(kElectron), TMath::Abs(kElectron));
-  if(!kfMother) return kFALSE;
-
-  AliESDtrack* d[2];
-  d[0] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(pIndex));
-  d[1] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(nIndex));
-
-  Float_t iMass = v0->GetEffMass(0, 0);
-
-  // cos pointing angle
-  Double_t cosPoint = v0->GetV0CosineOfPointingAngle();
-  cosPoint = TMath::ACos(cosPoint);
-
-  // DCA between daughters
-  Double_t dca = v0->GetDcaV0Daughters();
-
-  // Production vertex
-  Double_t x, y, z; 
-  v0->GetXYZ(x,y,z);
-  Double_t r = TMath::Sqrt(x*x + y*y);
-
-  Double_t xy[2];
-  Double_t r2 = -1.;
-  if ( GetConvPosXY(d[0], d[1], xy) ){
-    r2 = TMath::Sqrt(xy[0]*xy[0] + xy[1]*xy[1]);
-  }
-
-  // psi pair 
-  Double_t psiPair = PsiPair(v0);
-  
-  // V0 chi2/ndf
-  Double_t chi2ndf = kfMother->GetChi2()/kfMother->GetNDF();
-
-  if(kfMother) delete kfMother; 
-  
-  // apply the cuts
-
-  if(iMass > fGcutInvMass) return kFALSE;
-
-  if(chi2ndf > fGcutChi2NDF) return kFALSE;
-
-  if(cosPoint < fGcutCosPoint[0] || cosPoint > fGcutCosPoint[1]) return kFALSE;
-
-  if(dca < fGcutDCA[0] || dca > fGcutDCA[1]) return kFALSE;
-
-  if(r < fGcutVertexR[0] || r > fGcutVertexR[1]) return kFALSE;
-
-  if(psiPair < fGcutPsiPair[0] || psiPair > fGcutPsiPair[1]) return kFALSE;
-  
-  // all cuts passed
-
-  pdgV0 = 22;
-  if(sign){
-    pdgP = -11;
-    pdgN = 11;
-  }
-  else{
-    pdgP = 11;
-    pdgN = -11;
-  }
-
-  return kTRUE;
-}
-//____________________________________________________________________
-Bool_t  AliESDv0KineCuts::CaseK0(AliESDv0* const v0, Int_t &pdgV0, Int_t &pdgP, Int_t &pdgN){
-  //
-  // process the K0 candidate
-  //
-
-  if(!v0) return kFALSE;
-  
-  AliVTrack* daughter[2];
-  Int_t pIndex = 0, nIndex = 0;
-  Bool_t sign = CheckSigns(v0);
-  if(sign){
-    pIndex = v0->GetPindex();
-    nIndex = v0->GetNindex();
-  }
-  else{
-    pIndex = v0->GetNindex();
-    nIndex = v0->GetPindex();    
-  }
-  daughter[0] = dynamic_cast<AliVTrack *>(fEvent->GetTrack(pIndex));
-  daughter[1] = dynamic_cast<AliVTrack *>(fEvent->GetTrack(nIndex));
-  if(!daughter[0] || !daughter[1]) return kFALSE;
-
-  AliKFParticle *kfMother = CreateMotherParticle(daughter[0], daughter[1], TMath::Abs(kPiPlus), TMath::Abs(kPiPlus));
-  if(!kfMother) return kFALSE;
-
-  AliESDtrack* d[2];
-  d[0] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(pIndex));
-  d[1] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(nIndex));
-
-  Float_t iMass = v0->GetEffMass(2, 2);
-
-  // cos pointing angle
-  Double_t cosPoint = v0->GetV0CosineOfPointingAngle();
-  cosPoint = TMath::ACos(cosPoint);
-
-  // DCA between daughters
-  Double_t dca = v0->GetDcaV0Daughters();
-
-  // Production vertex
-  Double_t x, y, z; 
-  v0->GetXYZ(x,y,z);
-
-  Double_t r = TMath::Sqrt(x*x + y*y);  
-
-  // V0 chi2/ndf
-  Double_t chi2ndf = kfMother->GetChi2()/kfMother->GetNDF();
-  
-  if(kfMother) delete kfMother; 
-
-  //
-  // apply the cuts
-  //
-  if(iMass < fK0cutInvMass[0] || iMass > fK0cutInvMass[1]) return kFALSE;
-
-  if(chi2ndf > fK0cutChi2NDF) return kFALSE;
-
-  if(cosPoint < fK0cutCosPoint[0] || cosPoint > fK0cutCosPoint[1]) return kFALSE;
-
-  if(dca < fK0cutDCA[0] || dca > fK0cutDCA[1]) return kFALSE;
-
-  if(r < fK0cutVertexR[0] || r > fK0cutVertexR[1]) return kFALSE;
-
-  // all cuts passed
-  pdgV0 = 310;
-  if(sign){
-    pdgP = 211;
-    pdgN = -211;
-  }
-  else{
-    pdgP = -211;
-    pdgN = 211;
-  }
-
-  return kTRUE;
-}
-//____________________________________________________________________
-Bool_t  AliESDv0KineCuts::CaseLambda(AliESDv0* const v0, Int_t &pdgV0, Int_t &pdgP, Int_t &pdgN, Int_t id){
-  //
-  // process teh Lambda and Anti-Lambda candidate
-  //
-  
-  if(!v0) return kFALSE;
-
-    const Double_t cL0mass=TDatabasePDG::Instance()->GetParticle(kLambda0)->Mass();  // PDG lambda mass
-
-  AliVTrack* daughter[2];
-  Int_t pIndex = 0, nIndex = 0;
-  Float_t mMass[2] = {-1., -1.};
-  Bool_t sign = CheckSigns(v0);
-  if(sign){
-    pIndex = v0->GetPindex();
-    nIndex = v0->GetNindex();
-    mMass[0] = v0->GetEffMass(4, 2);
-    mMass[1] = v0->GetEffMass(2, 4);
-  }
-  else{
-    pIndex = v0->GetNindex();
-    nIndex = v0->GetPindex();    
-    mMass[0] = v0->GetEffMass(2, 4);
-    mMass[1] = v0->GetEffMass(4, 2);
-  }
-  daughter[0] = dynamic_cast<AliVTrack *>(fEvent->GetTrack(pIndex));
-  daughter[1] = dynamic_cast<AliVTrack *>(fEvent->GetTrack(nIndex));
-  if(!daughter[0] || !daughter[1]) return kFALSE;
-
-  AliKFParticle *kfMother[2] = {0x0, 0x0};
-  // Lambda
-  kfMother[0] = CreateMotherParticle(daughter[0], daughter[1], TMath::Abs(kProton), TMath::Abs(kPiPlus));
-  if(!kfMother[0]) return kFALSE;
-  
-  // Anti-Lambda
-  kfMother[1] = CreateMotherParticle(daughter[0], daughter[1], TMath::Abs(kPiPlus), TMath::Abs(kProton));
-  if(!kfMother[1]) return kFALSE;
-
-  Float_t dMass[2] = {TMath::Abs(mMass[0] - cL0mass), TMath::Abs(mMass[1] - cL0mass)};
-  
-  AliESDtrack* d[2];
-  d[0] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(pIndex));
-  d[1] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(nIndex));
-  if(!d[0] || !d[1])    return kFALSE;
-  
-  Float_t p[2] = {d[0]->GetP(), d[1]->GetP()}; 
-
-  // check the 3 lambda - antilambda variables
-  Int_t check[2] = {-1, -1};   // 0 : lambda, 1 : antilambda
-  // 1) momentum of the daughter particles - proton is expected to have higher momentum than pion
-  check[0] = (p[0] > p[1]) ? 0 : 1;
-  // 2) mass of the mother particle
-  check[1] = (dMass[0] < dMass[1]) ? 0 : 1;
-  // require positive correlation of (1) and (2)
-  if(check[0] != check[1]){
-    if(kfMother[0]) delete kfMother[0]; 
-    if(kfMother[1]) delete kfMother[1]; 
-    return kFALSE;
-  }
-
-  // now that the check[0] == check[1]
-  const Int_t type = check[0];
-
-  // require that the input armenteros preselection agree:
-  if(type != id) return kFALSE;
-
-  Float_t iMass =0.;
-  if(sign){
-    iMass = (type == 0) ? v0->GetEffMass(4, 2) : v0->GetEffMass(2, 4);
-  }
-  else{
-    iMass = (type == 0) ? v0->GetEffMass(2, 4) : v0->GetEffMass(4, 2);
-  }
-
-  // cos pointing angle
-  Double_t cosPoint = v0->GetV0CosineOfPointingAngle();
-  cosPoint = TMath::ACos(cosPoint);
-
-  // DCA between daughters
-  Double_t dca = v0->GetDcaV0Daughters();
-  
-  // Production vertex
-  Double_t x, y, z; 
-  v0->GetXYZ(x,y,z);
-  Double_t r = TMath::Sqrt(x*x + y*y);
-
-  // proton - pion indices
-  Int_t ix[2] = {0, 1};
-  if(1 == type){
-    ix[0] = 1;
-    ix[1] = 0;
-  }
-
-  // V0 chi2/ndf
-  Double_t chi2ndf = kfMother[type]->GetChi2()/kfMother[type]->GetNDF();
-
-  if(kfMother[0]) delete kfMother[0]; 
-  if(kfMother[1]) delete kfMother[1]; 
-
-  //
-  // apply the cuts
-  //
-
-  if(iMass < fLcutInvMass[0] || iMass > fLcutInvMass[1]) return kFALSE;
-
-  if(chi2ndf > fLcutChi2NDF) return kFALSE;
-
-  if(cosPoint < fLcutCosPoint[0] || cosPoint > fLcutCosPoint[1]) return kFALSE;
-
-  if(dca < fLcutDCA[0] || dca > fLcutDCA[1]) return kFALSE;
-
-  if(r < fLcutVertexR[0] || r > fLcutVertexR[1]) return kFALSE;
-
-  // all cuts passed
-
-  if(0 == type){
-    pdgV0 = 3122;
-    if(sign){
-      pdgP = 2212;
-      pdgN = -211;
-    }
-    else{
-      pdgP = -211;
-      pdgN = 2212;
-    }
-  }
-  else{
-    pdgV0 = -3122;
-    if(sign){
-      pdgP = 211;
-      pdgN = -2212;
-    }
-    else{
-      pdgP = -2212;
-      pdgN = 211;
-    }
-  }
-
-  return kTRUE;
-}
-//____________________________________________________________________
-Bool_t AliESDv0KineCuts::V0CutsCommon(const AliESDv0 * const v0){
-  //
-  // V0 cuts common to all V0s
-  //
-
-  AliESDtrack* dN, *dP; 
-  dP = dynamic_cast<AliESDtrack *>(fEvent->GetTrack(v0->GetPindex()));
-  dN = dynamic_cast<AliESDtrack *>(fEvent->GetTrack(v0->GetNindex())); 
-  
-  if(!dN || !dP) return kFALSE;
-
-  Int_t qP = dP->Charge();
-  Int_t qN = dN->Charge();
-
-  if((qP*qN) != -1) return kFALSE;
-
-  return kTRUE;
-}
-//____________________________________________________________________
-void   AliESDv0KineCuts::Armenteros(AliESDv0* const v0, Float_t val[2]){
-  //
-  // computes the Armenteros variables for given V0
-  // fills the histogram
-  // returns the values via "val"
-  //
-  
-  Double_t mn[3] = {0,0,0};
-  Double_t mp[3] = {0,0,0};  
-  Double_t mm[3] = {0,0,0};  
-
-  if(CheckSigns(v0)){
-    v0->GetNPxPyPz(mn[0],mn[1],mn[2]); //reconstructed cartesian momentum components of negative daughter
-    v0->GetPPxPyPz(mp[0],mp[1],mp[2]); //reconstructed cartesian momentum components of positive daughter
-  }
-  else{
-    v0->GetPPxPyPz(mn[0],mn[1],mn[2]); //reconstructed cartesian momentum components of negative daughter
-    v0->GetNPxPyPz(mp[0],mp[1],mp[2]); //reconstructed cartesian momentum components of positive daughter
-  }
-  v0->GetPxPyPz(mm[0],mm[1],mm[2]); //reconstructed cartesian momentum components of mother
-
-  TVector3 vecN(mn[0],mn[1],mn[2]);
-  TVector3 vecP(mp[0],mp[1],mp[2]);
-  TVector3 vecM(mm[0],mm[1],mm[2]);
-  
-  Double_t thetaP = acos((vecP * vecM)/(vecP.Mag() * vecM.Mag()));
-  Double_t thetaN = acos((vecN * vecM)/(vecN.Mag() * vecM.Mag()));
-  
-  Double_t alfa = ((vecP.Mag())*cos(thetaP)-(vecN.Mag())*cos(thetaN))/
-    ((vecP.Mag())*cos(thetaP)+(vecN.Mag())*cos(thetaN)) ;
-  Double_t qt = vecP.Mag()*sin(thetaP);
-
-  val[0] = alfa;
-  val[1] = qt;
-}
-//____________________________________________________________________
-Bool_t AliESDv0KineCuts::CheckSigns(AliESDv0* const v0){
-  //
-  // check wheter the sign was correctly applied to 
-  // V0 daughter tracks
-  //
-  
-  Bool_t correct = kFALSE;
-
-  Int_t pIndex = 0, nIndex = 0;
-  pIndex = v0->GetPindex();
-  nIndex = v0->GetNindex();
-  
-  AliESDtrack* d[2];
-  d[0] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(pIndex));
-  d[1] = dynamic_cast<AliESDtrack*>(fEvent->GetTrack(nIndex));
-
-  Int_t sign[2];
-  sign[0] = (int)d[0]->GetSign();
-  sign[1] = (int)d[1]->GetSign();
-  
-  if(-1 == sign[0] && 1 == sign[1]){
-    correct = kFALSE;
-  }
-  else{
-    correct = kTRUE;
-  }
-  
-  return correct;
-}
-//________________________________________________________________
-Double_t AliESDv0KineCuts::PsiPair(AliESDv0* const v0) {
-  //
-  // Angle between daughter momentum plane and plane 
-  // 
-
-  if(!fEvent) return -1.;
-
-  Float_t magField = fEvent->GetMagneticField();
-
-  Int_t pIndex = -1;
-  Int_t nIndex = -1;
-  if(CheckSigns(v0)){
-    pIndex = v0->GetPindex();
-    nIndex = v0->GetNindex();
-  }
-  else{
-    pIndex = v0->GetNindex();
-    nIndex = v0->GetPindex();    
-  }
-
-  AliESDtrack* daughter[2];
-
-  daughter[0] = dynamic_cast<AliESDtrack *>(fEvent->GetTrack(pIndex));
-  daughter[1] = dynamic_cast<AliESDtrack *>(fEvent->GetTrack(nIndex));
-
-  Double_t x, y, z;
-  v0->GetXYZ(x,y,z);//Reconstructed coordinates of V0; to be replaced by Markus Rammler's method in case of conversions!
-  
-  Double_t mn[3] = {0,0,0};
-  Double_t mp[3] = {0,0,0};
-  
-
-  v0->GetNPxPyPz(mn[0],mn[1],mn[2]);//reconstructed cartesian momentum components of negative daughter;
-  v0->GetPPxPyPz(mp[0],mp[1],mp[2]);//reconstructed cartesian momentum components of positive daughter; 
-
-
-  Double_t deltat = 1.;
-  deltat = TMath::ATan(mp[2]/(TMath::Sqrt(mp[0]*mp[0] + mp[1]*mp[1])+1.e-13)) -  TMath::ATan(mn[2]/(TMath::Sqrt(mn[0]*mn[0] + mn[1]*mn[1])+1.e-13));//difference of angles of the two daughter tracks with z-axis
-
-  Double_t radiussum = TMath::Sqrt(x*x + y*y) + 50;//radius to which tracks shall be propagated
-
-  Double_t momPosProp[3];
-  Double_t momNegProp[3];
-    
-  AliExternalTrackParam pt(*daughter[0]), nt(*daughter[1]);
-    
-  Double_t psiPair = 4.;
-
-  if(nt.PropagateTo(radiussum,magField) == 0)//propagate tracks to the outside
-    psiPair =  -5.;
-  if(pt.PropagateTo(radiussum,magField) == 0)
-    psiPair = -5.;
-  pt.GetPxPyPz(momPosProp);//Get momentum vectors of tracks after propagation
-  nt.GetPxPyPz(momNegProp);
-  
-  Double_t pEle =
-    TMath::Sqrt(momNegProp[0]*momNegProp[0]+momNegProp[1]*momNegProp[1]+momNegProp[2]*momNegProp[2]);//absolute momentum value of negative daughter
-  Double_t pPos =
-    TMath::Sqrt(momPosProp[0]*momPosProp[0]+momPosProp[1]*momPosProp[1]+momPosProp[2]*momPosProp[2]);//absolute momentum value of positive daughter
-    
-  Double_t scalarproduct =
-    momPosProp[0]*momNegProp[0]+momPosProp[1]*momNegProp[1]+momPosProp[2]*momNegProp[2];//scalar product of propagated positive and negative daughters' momenta
-    
-  Double_t chipair = TMath::ACos(scalarproduct/(pEle*pPos));//Angle between propagated daughter tracks
-
-  psiPair =  TMath::Abs(TMath::ASin(deltat/chipair));  
-
-  return psiPair; 
-}
-//___________________________________________________________________
-Bool_t AliESDv0KineCuts::GetConvPosXY(AliESDtrack * const ptrack, AliESDtrack * const ntrack, Double_t convpos[2]){
-  //
-  // recalculate the gamma conversion XY postition
-  //
-
-  const Double_t b = fEvent->GetMagneticField();
-
-  Double_t helixcenterpos[2];
-  GetHelixCenter(ptrack,b,ptrack->Charge(),helixcenterpos);
-
-  Double_t helixcenterneg[2];
-  GetHelixCenter(ntrack,b,ntrack->Charge(),helixcenterneg);
-
-  Double_t  poshelix[6];
-  ptrack->GetHelixParameters(poshelix,b);
-  Double_t posradius = TMath::Abs(1./poshelix[4]);
-
-  Double_t  neghelix[6];
-  ntrack->GetHelixParameters(neghelix,b);
-  Double_t negradius = TMath::Abs(1./neghelix[4]);
-
-  Double_t xpos = helixcenterpos[0];
-  Double_t ypos = helixcenterpos[1];
-  Double_t xneg = helixcenterneg[0];
-  Double_t yneg = helixcenterneg[1];
-
-  convpos[0] = (xpos*negradius + xneg*posradius)/(negradius+posradius);
-  convpos[1] = (ypos*negradius+  yneg*posradius)/(negradius+posradius);
-
-  return 1;
-}
-//___________________________________________________________________
-Bool_t AliESDv0KineCuts::GetHelixCenter(AliESDtrack * const track, Double_t b,Int_t charge, Double_t center[2]){
-  //
-  // computes the center of the track helix
-  //
-  
-  Double_t pi = TMath::Pi();
-  
-  Double_t  helix[6];
-  track->GetHelixParameters(helix,b);
-  
-  Double_t xpos =  helix[5];
-  Double_t ypos =  helix[0];
-  Double_t radius = TMath::Abs(1./helix[4]);
-  Double_t phi = helix[2];
-
-  if(phi < 0){
-    phi = phi + 2*pi;
-  }
-
-  phi -= pi/2.;
-  Double_t xpoint =  radius * TMath::Cos(phi);
-  Double_t ypoint =  radius * TMath::Sin(phi);
-
-  if(b<0){
-    if(charge > 0){
-      xpoint = - xpoint;
-      ypoint = - ypoint;
-    }
-
-    if(charge < 0){
-      xpoint =  xpoint;
-      ypoint =  ypoint;
-    }
-  }
-  if(b>0){
-    if(charge > 0){
-      xpoint =  xpoint;
-      ypoint =  ypoint;
-    }
-
-    if(charge < 0){
-      xpoint = - xpoint;
-      ypoint = - ypoint;
-    }
-  }
-  center[0] =  xpos + xpoint;
-  center[1] =  ypos + ypoint;
-
-  return 1;
-}
-//___________________________________________________________________
-AliKFParticle *AliESDv0KineCuts::CreateMotherParticle(const AliVTrack* const pdaughter, const AliVTrack* const ndaughter, Int_t pspec, Int_t nspec){
-  //
-  // Creates a mother particle
-  //
-  AliKFParticle pkfdaughter(*pdaughter, pspec);
-  AliKFParticle nkfdaughter(*ndaughter, nspec);
-  
-  
-  // Create the mother particle 
-  AliKFParticle *m = new AliKFParticle(pkfdaughter, nkfdaughter);
-  // DEBUG - testing
-  if(TMath::Abs(kElectron) == pspec && TMath::Abs(kElectron) == nspec) m->SetMassConstraint(0, 0.001);
-  else if(TMath::Abs(kPiPlus) == pspec && TMath::Abs(kPiPlus) == nspec) m->SetMassConstraint(TDatabasePDG::Instance()->GetParticle(kK0Short)->Mass(), 0.);
-  else if(TMath::Abs(kProton) == pspec && TMath::Abs(kPiPlus) == nspec) m->SetMassConstraint(TDatabasePDG::Instance()->GetParticle(kLambda0)->Mass(), 0.);
-  else if(TMath::Abs(kPiPlus) == pspec && TMath::Abs(kProton) == nspec) m->SetMassConstraint(TDatabasePDG::Instance()->GetParticle(kLambda0)->Mass(), 0.);
-  else{
-    AliErrorClass("Wrong daughter ID - mass constraint can not be set");
-  }
-
-  AliKFVertex improvedVertex = *fPrimaryVertex;
-  improvedVertex += *m;
-  m->SetProductionVertex(improvedVertex);
-  
-  // update 15/06/2010
-  // mother particle will not be added to primary vertex but only to its copy 
-  // as this confilcts with calling
-  // m->SetPrimaryVertex() function and
-  // subsequently removing the mother particle afterwards
-  // Source: Sergey Gorbunov
-
-  return m;
-}
-//____________________________________________________________________
-void  AliESDv0KineCuts::SetEvent(AliESDEvent* const event){
-  //
-  // direct setter of ESD event
-  //
-  fEvent = event;
-  if(!fEvent){
-    AliErrorClass("Invalid input event pointer");
-    return;
-  }
-
-}
-//____________________________________________________________________
-void  AliESDv0KineCuts::SetEvent(AliVEvent* const event){
-  //
-  // direct setter of ESD event
-  //
-  if(event)
-    fEvent = static_cast<AliESDEvent*>(event);
-  if(!fEvent){
-    AliErrorClass("Invalid input event pointer");
-    return;
-  }
-
-}
-//________________________________________________________________
-void AliESDv0KineCuts::SetPrimaryVertex(AliKFVertex* const v){
-  //
-  // set the primary vertex of the event
-  //
-  fPrimaryVertex = v;
-  if(!fPrimaryVertex){
-    AliErrorClass("Failed to initialize the primary vertex");
-    return;
-  }
-}
-//___________________________________________________________________
-void AliESDv0KineCuts::SetMode(Int_t mode, Int_t type){
-  //
-  // this function allows the user to select (prior running the 'ProcessV0' function)
-  // to select different approaches to V0 selection - the 'mode'
-  // - and -
-  // different systems (pp, PbPb) - 'type' 
-  //
-  // To see the cut values for different modes please refer to the
-  // function SetCuts()
-  //
-  // Important notice: based on the parameters particular sets of cuts will
-  // be activated for teh V0 selection. If some additional changes to single
-  // cuts are needed please us the SetXXXcut function (see the header file)
-  // 
-
-  switch(mode){
-  case kPurity:
-    fMode = kPurity;  // used to obtain highest purity possible - the efficiency may be low
-  case kEffGamma:
-    fMode = kEffGamma; // used to obtain highes efficiency possible - the purity may be worse
-  default:
-    AliError("V0 selection mode not recognozed, setting 'kPurity'");
-    fMode = kPurity;
-  }
-
-  switch(type){
-  case kPP:
-    fType = kPP;  // cuts optimized for low multiplicity 
-  case kPbPb:
-    fType = kPbPb;  // cuts optimized for high multiplicity
-  }
-  
-  // setup the cut values for selected mode & type
-  SetCuts();
-
-}
-//___________________________________________________________________
-void AliESDv0KineCuts::SetMode(Int_t mode, const char* type){
-  //
-  // overloaded function - please see above
-  // 
-  
-  Int_t t = -1;
-
-  if(!strcmp("pp", type)) t = kPP;
-  else if(!(strcmp("PbPb", type))) t = kPbPb;
-  else{
-    AliError("data type not recognized, setting 'pp'");
-    t = kPP;    
-  }
-
-  SetMode(mode, t);
-
-}
-//___________________________________________________________________
-void AliESDv0KineCuts::SetCuts(){
-  //
-  // this funciton sets the default cut values based on the selected
-  // fMode and fType.
-  // please note that only the cuts that have different values than the default
-  // cuts are updated here
-  //
-  
-  // last update: 14/02/2011
-  // as a very preliminary  - the only change to default cuts is to apply
-  // less restricting gamma conversion selection in PreselectV0() function
-  
-
-  
-}
diff --git a/PWG3/hfe/AliESDv0KineCuts.h b/PWG3/hfe/AliESDv0KineCuts.h
deleted file mode 100644 (file)
index d59233f..0000000
+++ /dev/null
@@ -1,195 +0,0 @@
-#ifndef ALIESDV0KINECUTS_H
-#define ALIESDV0KINECUTS_H
-
-/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
- * See cxx source for full Copyright notice                               */
-
-/* $Id$ */ 
-
-/*
- * plesae see source file for more details
- */
-
-#include <TObject.h>
-
-class AliESDv0;
-class AliESDEvent;
-class AliVEvent;
-class AliESDtrack;
-class AliVTrack;
-class AliKFParticle;
-class AliKFVertex;
-
-class AliESDv0KineCuts : public TObject{
- public:
-  enum{ // Reconstructed V0
-    kUndef = -1,
-      kGamma = 0,
-      kK0 = 1,
-      kLambda = 2,
-      kALambda = 3
-      };
-  enum{ // data types
-    kPP = 0,
-      kPbPb = 1,  // not yet implemented
-      };
-  enum{ // operation modes
-    kPurity = 0, // purely kinematical selection
-      kEffGamma = 1  // !!! involves TPC dEdx or nSimga cuts !!!
-      };
-  
-  AliESDv0KineCuts();
-  virtual ~AliESDv0KineCuts();
-
-  AliESDv0KineCuts(const AliESDv0KineCuts &ref);
-  AliESDv0KineCuts &operator=(const AliESDv0KineCuts &ref);
-
-  // main selection function - called once per V0 candidate
-  Bool_t ProcessV0(AliESDv0* const v0, Int_t &pdgV0, Int_t &pdgP, Int_t &pdgN);
-  Bool_t ProcessV0(AliESDv0* const v0, Int_t &pdgP, Int_t &pdgN);
-
-  // must be called by the user
-  void SetEvent(AliESDEvent* const event);
-  void SetEvent(AliVEvent* const event);
-  void SetPrimaryVertex(AliKFVertex* const v);
-
-  // user can select an operation modes [see .cxx for details]
-  void   SetMode(Int_t mode, Int_t type);
-  void   SetMode(Int_t mode, const char* type);
-
-  //
-  // setter functions for V0 cut values
-  // for default values see the constructor
-  // see the default contructor for comments
-  //
-
-  // single track cuts
-  void   SetNTPCclusters(Int_t n) { fTPCNcls = n; };
-  void   SetTPCrefit(Bool_t r = kTRUE) { fTPCrefit = r; };
-  void   SetTPCchi2perCls(Float_t chi2) { fTPCchi2perCls = chi2; };
-  void   SetTPCclusterratio(Float_t r) { fTPCclsRatio = r; };
-  void   SetNoKinks(Bool_t k = kTRUE) { fNoKinks = k; };
-
-  // gamma cuts
-  void   SetGammaCutChi2NDF(Float_t val)  { fGcutChi2NDF = val; };
-  void   SetGammaCutCosPoint(Float_t *val) { 
-    fGcutCosPoint[0] = val[0];
-    fGcutCosPoint[1] = val[1];
-  };
-  void   SetGammaCutDCA(Float_t *val){
-    fGcutDCA[0] = val[0];
-    fGcutDCA[1] = val[1];
-  };
-  void   SetGammaCutVertexR(Float_t *val){
-    fGcutVertexR[0] = val[0];
-    fGcutVertexR[1] = val[1];
-  };
-  void   SetGammaCutPsiPair(Float_t *val){
-    fGcutPsiPair[0] = val[0];
-    fGcutPsiPair[1] = val[1];
-  };
-  void   SetGammaCutInvMass(Float_t val){
-    fGcutInvMass = val;
-  };
-  // K0 cuts
-  void   SetK0CutChi2NDF(Float_t val)  { fK0cutChi2NDF = val; };
-  void   SetK0CutCosPoint(Float_t *val) { 
-    fK0cutCosPoint[0] = val[0];
-    fK0cutCosPoint[1] = val[1];
-  };
-  void   SetK0CutDCA(Float_t *val){
-    fK0cutDCA[0] = val[0];
-    fK0cutDCA[1] = val[1];
-  };
-  void   SetK0CutVertexR(Float_t *val){
-    fK0cutVertexR[0] = val[0];
-    fK0cutVertexR[1] = val[1];
-  };
-  void   SetK0CutInvMass(Float_t *val){
-    fK0cutInvMass[0] = val[0];
-    fK0cutInvMass[1] = val[1];
-  };
-  // lambda & anti-lambda cuts
-  void   SetLambdaCutChi2NDF(Float_t val)  { fLcutChi2NDF = val; };
-  void   SetLambdaCutCosPoint(Float_t *val) { 
-    fLcutCosPoint[0] = val[0];
-    fLcutCosPoint[1] = val[1];
-  };
-  void   SetLambdaCutDCA(Float_t *val){
-    fLcutDCA[0] = val[0];
-    fLcutDCA[1] = val[1];
-  };
-  void   SetLambdaCutVertexR(Float_t *val){
-    fLcutVertexR[0] = val[0];
-    fLcutVertexR[1] = val[1];
-  };
-  void   SetLambdaCutInvMass(Float_t *val){
-    fLcutInvMass[0] = val[0];
-    fLcutInvMass[1] = val[1];
-  };
-  
-
-  Int_t  PreselectV0(AliESDv0* const v0);
-
-  Bool_t CaseGamma(AliESDv0* const v0, Int_t &pdgV0, Int_t &pdgP, Int_t &pdgN);
-  Bool_t CaseK0(AliESDv0* const v0, Int_t &pdgV0, Int_t &pdgP, Int_t &pdgN);
-  Bool_t CaseLambda(AliESDv0* const v0, Int_t &pdgV0, Int_t &pdgP, Int_t &pdgN, Int_t id);
-
-  Bool_t V0CutsCommon(const AliESDv0 * const v0);
-  Bool_t SingleTrackCuts(AliESDv0 * const v0);
-  void   Armenteros(AliESDv0* const v0, Float_t val[2]);
-  Bool_t CheckSigns(AliESDv0* const v0);
-
-  Double_t PsiPair(AliESDv0* const v0);
-  Bool_t   GetConvPosXY(AliESDtrack * const ptrack, AliESDtrack * const ntrack, Double_t convpos[2]);
-  Bool_t   GetHelixCenter(AliESDtrack * const track, Double_t b, Int_t charge, Double_t center[2]);
-
- protected:
-  void Copy(TObject &ref) const;
-
- private:
-
-  AliKFParticle *CreateMotherParticle(const AliVTrack* const pdaughter, const AliVTrack* const ndaughter, Int_t pspec, Int_t nspec);
-  void      SetCuts();                          // setup cuts for selected fMode and fType, see source file for details
-  Bool_t    GammaEffCuts(AliESDv0 * const v0);  // set of cuts optimized for high gamma efficiency
-
- private:
-  AliESDEvent           *fEvent;          // current event
-  AliKFVertex           *fPrimaryVertex;  // primary vertex
-
-  Int_t                 fType;            // data type: p-p or Pb-Pb
-  Int_t                 fMode;            // current operation mode
-
-  // single track cuts
-  Int_t                 fTPCNcls;          // number of TPC clusters
-  Bool_t                fTPCrefit;         // TPC refit - yes [kTRUE] or do not care [kFALSE]
-  Float_t               fTPCchi2perCls;    // max. chi2 per TPC cluster
-  Float_t               fTPCclsRatio;      // min. TPC cluster ratio
-  Bool_t                fNoKinks;          // kinks - no [kTRUE] or do not care [kFalse]
-
-  // gamma cut values
-  Float_t               fGcutChi2NDF;      // Chi2NF cut value for the AliKFparticle gamma
-  Float_t               fGcutCosPoint[2];  // cos of the pointing angle [min, max]
-  Float_t               fGcutDCA[2];       // DCA between the daughter tracks [min, max]
-  Float_t               fGcutVertexR[2];   // radius of the conversion point [min, max]
-  Float_t               fGcutPsiPair[2];   // value of the psi pair cut [min, max]
-  Float_t               fGcutInvMass;      // upper value on the gamma invariant mass
-  // K0 cut values
-  Float_t               fK0cutChi2NDF;     // Chi2NF cut value for the AliKFparticle K0
-  Float_t               fK0cutCosPoint[2]; // cos of the pointing angle [min, max]
-  Float_t               fK0cutDCA[2];      // DCA between the daughter tracks [min, max]
-  Float_t               fK0cutVertexR[2];  // radius of the decay point [min, max]
-  Float_t               fK0cutInvMass[2];  // invariant mass window
-  // Lambda & anti-Lambda cut values
-  Float_t               fLcutChi2NDF;      // Chi2NF cut value for the AliKFparticle K0
-  Float_t               fLcutCosPoint[2];  // cos of the pointing angle [min, max]
-  Float_t               fLcutDCA[2];       // DCA between the daughter tracks [min, max]
-  Float_t               fLcutVertexR[2];   // radius of the decay point [min, max]
-  Float_t               fLcutInvMass[2];   // invariant mass window
-  
-  
-  ClassDef(AliESDv0KineCuts, 0);
-
-};
-
-#endif