LmvmCand.h

Go to the documentation of this file.

/* Copyright (C) 2015-2021 Justus-Liebig-Universitaet Giessen, Giessen
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Elena Lebedeva [committer], Semen Lebedev */
 
#ifndef LMVM_CAND_H
#define LMVM_CAND_H
 
#include "LmvmDef.h"
#include "TVector3.h"
 
class LmvmCand {
 public:
  LmvmCand() {}
 
  void ResetMcParams()
  {
    fMcSrc         = ELmvmSrc::Undefined;
    fMcMotherId    = -1;
    fStsMcTrackId  = -1;
    fRichMcTrackId = -1;
    fTrdMcTrackId  = -1;
    fTofMcTrackId  = -1;
  }
 
  void SetIsTopologyCutElectron(ELmvmTopologyCut cut, bool value)
  {
    if (cut == ELmvmTopologyCut::TT) {
      fIsTtCut = value;
    }
    if (cut == ELmvmTopologyCut::ST) {
      fIsStCut = value;
    }
    if (cut == ELmvmTopologyCut::RT) {
      fIsRtCut = value;
    }
  }
 
  bool IsCutTill(ELmvmAnaStep step) const
  {
    if (step == ELmvmAnaStep::Mc) return true;
    if (step == ELmvmAnaStep::Reco && fIsRec) return true;
    if (step == ELmvmAnaStep::Acc && IsCutTill(ELmvmAnaStep::Reco) && fIsAcc) return true;
    if (step == ELmvmAnaStep::Chi2Prim && IsCutTill(ELmvmAnaStep::Acc) && fIsChi2Prim) return true;
    if (step == ELmvmAnaStep::ElId && IsCutTill(ELmvmAnaStep::Chi2Prim) && fIsElectron) return true;
    if (step == ELmvmAnaStep::GammaCut && IsCutTill(ELmvmAnaStep::ElId) && fIsGammaCut) return true;
    if (step == ELmvmAnaStep::Mvd1Cut && IsCutTill(ELmvmAnaStep::GammaCut) && fIsMvd1Cut) return true;
    if (step == ELmvmAnaStep::Mvd2Cut && IsCutTill(ELmvmAnaStep::Mvd1Cut) && fIsMvd2Cut) return true;
    if (step == ELmvmAnaStep::StCut && IsCutTill(ELmvmAnaStep::Mvd2Cut) && fIsStCut) return true;
    if (step == ELmvmAnaStep::RtCut && IsCutTill(ELmvmAnaStep::StCut) && fIsRtCut) return true;
    if (step == ELmvmAnaStep::TtCut && IsCutTill(ELmvmAnaStep::RtCut) && fIsTtCut) return true;
    if (step == ELmvmAnaStep::PtCut && IsCutTill(ELmvmAnaStep::TtCut) && fIsPtCut) return true;
    return false;
  }
 
  // track parameters
  TVector3 fPosition;
  TVector3 fMomentum;
  double fMass     = -1.;  // mass from MC
  double fMassSig  = -1.;  // mass of signal mother (remains '-1' if candidate is not signal electron)
  double fWeight   = -1.;  // fW * MinvScale(double minv) for PLUTOs. For UrQMDs: fWeight = 1;
  double fEnergy   = 0.;
  double fRapidity = 0.;
  int fCharge      = 0;
  double fChi2Prim = 0.;
  double fChi2Sts  = 0.;
  double fChi2Rich = 0.;
  double fChi2Trd  = 0.;
  double fChi2Tof  = -1.;
  double fLength   = 0.;  // length of according global track
  double fTime     = 0.;
  double fTofDist  = -1.;  // Hit-Track-Distance in ToF
 
  // To investigate misidentifications in single sub detectors
  bool fIsRichElectron = false;
  bool fIsTrdElectron  = false;
  bool fIsTofElectron  = false;
 
  int fMcPdg         = -1;
  int fMcMotherId    = -1;
  int fEventNumber   = 0;
  int fStsMcTrackId  = -1;
  int fRichMcTrackId = -1;
  int fTrdMcTrackId  = -1;
  int fTofMcTrackId  = -1;
  int fStsInd        = -1;
  int fRichInd       = -1;
  int fTrdInd        = -1;
  int fTofHitInd     = -1;
  int fTofTrackInd   = -1;
  int fGTrackInd     = -1;
  int fNofHitsSts    = -1;  // TODO: delete nofHits (necessary)?
  int fNofHitsMvd    = -1;
  int fNofHitsRich   = -1;
  int fNofHitsTrd    = -1;
  int fNofHitsTof    = -1;
  double fRichAnn    = 0.;
  double fTrdAnn     = 0.;
  double fMass2      = -1.;  // square mass calculated by TOF
  double fTrdLikeEl  = -1.;
  double fTrdLikePi  = -1.;
  double fELossSts   = -1.;
 
  // Cuts. If true then cut is passed
  bool fIsRec      = false;
  bool fIsAcc      = false;
  bool fIsChi2Prim = false;
  bool fIsElectron = false;
  bool fIsGammaCut = true;  // Will be set to 'false' as soon as a partner with minv < 25 MeV is found
  bool fIsMvd1Cut  = false;
  bool fIsMvd2Cut  = false;
  bool fIsTtCut    = false;
  bool fIsStCut    = false;
  bool fIsRtCut    = false;
  bool fIsPtCut    = false;
 
  std::string fName = "";  // to distinguish between data from different folders in LmvmEventMix
  std::string fTaskId =
    "";  // to distinguish between events with same event number but different tasksId in LmvmEventMix
 
  // MC
  ELmvmSrc fMcSrc = ELmvmSrc::Undefined;
  bool IsMcSignal() const { return fMcSrc == ELmvmSrc::Signal; }
  bool IsMcPi0() const { return fMcSrc == ELmvmSrc::Pi0; }
  bool IsMcGamma() const { return fMcSrc == ELmvmSrc::Gamma; }
  bool IsMcEta() const { return fMcSrc == ELmvmSrc::Eta; }
  // for candidates BG is all candidates which are not signal
  bool fIsMcBg() const { return fMcSrc != ELmvmSrc::Signal; }
};
 
#endif