effc++ warnings

[u/mrichter/AliRoot.git] / CORRFW / AliCFV0TopoCuts.cxx

/**************************************************************************
 * 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.                  *
 **************************************************************************/


///////////////////////////////////////////////////////////////////////////
//          ----   CORRECTION FRAMEWORK   ----
// Class to cut on V0 topology
//   -> support for :
//                    DCA between V0 daughters
//                    V0 daughter impact parameters wrt primary vertex
//                    cosine of V0 pointing angle
//
///////////////////////////////////////////////////////////////////////////
// author : R. Vernet (renaud.vernet@cern.ch)
///////////////////////////////////////////////////////////////////////////

#include "AliCFV0TopoCuts.h"
#include "AliESDv0.h"
#include "AliESDtrack.h"
#include "AliCFPair.h"

ClassImp(AliCFV0TopoCuts)

//______________________________________________________________
AliCFV0TopoCuts::AliCFV0TopoCuts() : 
  AliCFCutBase(),
  fMaxDcaDaughters(1.e99),
  fMinDcaNeg(0),
  fMinDcaPos(0),
  fMinCosP(0)
{
  //
  //default constructor
  //
}

//______________________________________________________________
AliCFV0TopoCuts::AliCFV0TopoCuts(const Char_t* name, const Char_t* title) : 
  AliCFCutBase(name,title),
  fMaxDcaDaughters(1.e99),
  fMinDcaNeg(0),
  fMinDcaPos(0),
  fMinCosP(0)
{
  //
}

//______________________________________________________________
AliCFV0TopoCuts::AliCFV0TopoCuts(const AliCFV0TopoCuts& c) : 
  AliCFCutBase(c),
  fMaxDcaDaughters(c.fMaxDcaDaughters),
  fMinDcaNeg(c.fMinDcaNeg),
  fMinDcaPos(c.fMinDcaPos),
  fMinCosP(c.fMinCosP)
{
  //
  // copy constructor
  //
}

//______________________________________________________________
AliCFV0TopoCuts& AliCFV0TopoCuts::operator=(const AliCFV0TopoCuts& c) 
{
  //
  // assignment operator
  //

  if (this != &c) {
    AliCFCutBase::operator=(c) ;
    fMaxDcaDaughters = c.fMaxDcaDaughters ;
    fMinDcaNeg       = c.fMinDcaNeg ;
    fMinDcaPos       = c.fMinDcaPos ;
    fMinCosP         = c.fMinCosP ;
  }
  return *this ;
}

//______________________________________________________________
Bool_t AliCFV0TopoCuts::IsSelected(TObject *obj) {
  //
  // computes V0 topological variables to cut on and return true 
  // in case the V0 is accepted
  //

  if (!obj) return kFALSE ;

  TString className(obj->ClassName());
  if (className.CompareTo("AliCFPair") != 0) {
    Error("IsSelected","obj must point to an AliCFPair !");
    return kFALSE ;
  }

  AliCFPair * pair = dynamic_cast<AliCFPair*>(obj);
  AliESDv0    * v0          = pair->GetV0();
  AliESDtrack * negDaughter = pair->GetNeg();
  AliESDtrack * posDaughter = pair->GetPos();

  if (v0->GetDcaV0Daughters() > fMaxDcaDaughters)  return kFALSE ;
  if (v0->GetV0CosineOfPointingAngle() < fMinCosP) return kFALSE ;


  Float_t tDca[2];
  if (negDaughter) negDaughter->GetImpactParameters(tDca[0],tDca[1]);
  else { tDca[0]=1.e99;  tDca[1]=1.e99;}
  Double32_t negDca = TMath::Sqrt(tDca[0]*tDca[0]+tDca[1]*tDca[1]);
  if (posDaughter) posDaughter->GetImpactParameters(tDca[0],tDca[1]);
  else { tDca[0]=1.e99;  tDca[1]=1.e99;}
  Double32_t posDca = TMath::Sqrt(tDca[0]*tDca[0]+tDca[1]*tDca[1]);

  if (negDca < fMinDcaNeg) return kFALSE ;
  if (posDca < fMinDcaPos) return kFALSE ; 

  return kTRUE ;
}