CbmL1MCTrack.h

Go to the documentation of this file.

/*
 *====================================================================
 *
 *  CBM Level 1 Reconstruction 
 *  
 *  Authors: I.Kisel,  S.Gorbunov
 *
 *  e-mail : ikisel@kip.uni-heidelberg.de 
 *
 *====================================================================
 *
 *  L1 Monte Carlo information
 *
 *====================================================================
 */
 
#ifndef CbmL1MCTrack_H
#define CbmL1MCTrack_H
 
#include "CbmL1MCPoint.h"
#include "TLorentzVector.h"
#include "TVector3.h"
#include <iostream>
#include <vector>
using std::vector;
 
class CbmL1Track;
 
class CbmL1MCTrack {
public:
  double mass, q, p, x, y, z, px, py, pz, time;
  int ID, iFile, iEvent, mother_ID, pdg;
  vector<int> Points;   // indices of pints in L1::vMCPoints
  vector<int> StsHits;  // indices of hits in algo->vStsHits or L1::vStsHits
 
  CbmL1MCTrack()
    : mass(0.)
    , q(0.)
    , p(0.)
    , x(0.)
    , y(0.)
    , z(0.)
    , px(0.)
    , py(0.)
    , pz(0.)
    , time(0.)
    , ID(-1)
    , iFile(-1)
    , iEvent(-1)
    , mother_ID(-1)
    , pdg(-1)
    , Points()
    , StsHits()
    , nMCContStations(0)
    , nHitContStations(0)
    , maxNStaMC(0)
    , maxNSensorMC(0)
    , maxNStaHits(0)
    , nStations(0)
    , nMCStations(0)
    , isReconstructable(0)
    , isAdditional(0)
    , rTracks()
    , tTracks() {};
 
  CbmL1MCTrack(int _ID)
    : mass(0.)
    , q(0.)
    , p(0.)
    , x(0.)
    , y(0.)
    , z(0.)
    , px(0.)
    , py(0.)
    , pz(0.)
    , time(0.)
    , ID(_ID)
    , iFile(-1)
    , iEvent(-1)
    , mother_ID(-1)
    , pdg(-1)
    , Points()
    , StsHits()
    , nMCContStations(0)
    , nHitContStations(0)
    , maxNStaMC(0)
    , maxNSensorMC(0)
    , maxNStaHits(0)
    , nStations(0)
    , nMCStations()
    , isReconstructable(0)
    , isAdditional(0)
    , rTracks()
    , tTracks() {};
  CbmL1MCTrack(double mass,
               double q,
               TVector3 vr,
               TLorentzVector vp,
               int ID,
               int mother_ID,
               int pdg);
  //   CbmL1MCTrack(TmpMCPoints &mcPoint, TVector3 vr, TLorentzVector vp, int ID, int mother_ID);
 
  bool IsPrimary() const { return mother_ID < 0; };
  bool IsReconstructable() const { return isReconstructable; };
  bool IsAdditional() const { return isAdditional; }
  int NStations() const { return nStations; };
  int NMCStations() const { return nMCStations; };
  int NMCContStations() const { return nMCContStations; };
  float Fraction_MC();
 
  void Init();
 
  void AddRecoTrack(CbmL1Track* rTr) { rTracks.push_back(rTr); }
  vector<CbmL1Track*>& GetRecoTracks() { return rTracks; }
  int GetNClones() const { return rTracks.size() - 1; }
  bool IsReconstructed() const { return rTracks.size(); }
 
  void AddTouchTrack(CbmL1Track* tTr) { tTracks.push_back(tTr); }
  bool IsDisturbed() const { return tTracks.size(); }
 
  friend class CbmL1;
 
private:
  int nMCContStations;   // number of consecutive stations with mcPoints
  int nHitContStations;  // number of consecutive stations with hits
  int maxNStaMC;         // max number of mcPoints on station
  int maxNSensorMC;      // max number of mcPoints with same z
  int maxNStaHits;       // max number of hits on station
 
  int nStations;    // number of stations with hits
  int nMCStations;  // number of stations with MCPoints
 
  bool isReconstructable;
  bool isAdditional;  // is not reconstructable, but stil interesting
 
  void CalculateMCCont();
  void CalculateHitCont();
  void CalculateMaxNStaHits();
  void CalculateMaxNStaMC();
  void CalculateIsReconstructable();
 
  // next members filled and used in Performance
  vector<CbmL1Track*> rTracks;  // array of assosiated recoTracks
  vector<CbmL1Track*>
    tTracks;  // array of recoTracks wich aren't assosiated with this mcTrack, but use some hits from it.
};
 
 
#endif