PairAnalysisEvent.h

Go to the documentation of this file.

/* Copyright (C) 2015-2017 Justus-Liebig-Universitaet Giessen, Giessen
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Julian Book [committer] */
 
#ifndef PAIRANALYSISEVENT_H
#define PAIRANALYSISEVENT_H
 
//#############################################################
//#                                                           #
//#      Class PairAnalysisEvent                                #
//#      Class for event structure and references to tracks   #
//#                                                           #
//#  Authors:                                                 #
//#   Julian    Book,     Uni Ffm / Julian.Book@cern.ch       #
//#                                                           #
//#############################################################
#include "CbmDefs.h"
#include "CbmVertex.h"
 
#include "FairEventHeader.h"
 
#include <TClonesArray.h>
#include <TNamed.h>
#include <TObjArray.h>
 
class FairRootManager;
class FairMCEventHeader;
class CbmStsTrack;
class PairAnalysisTrack;
 
class PairAnalysisEvent : public TNamed {
public:
  PairAnalysisEvent();
  PairAnalysisEvent(const char* name, const char* title);
 
  virtual ~PairAnalysisEvent();
 
  void SetInput(FairRootManager* man);
  void SetANNId(Int_t Id) { fANNId = Id; }
  void Init();
 
 
  // getters
  FairMCEventHeader* GetMCHeader() const { return fMCHeader; }
  CbmVertex* GetPrimaryVertex() const { return fPrimVertex; }
  Int_t GetNumberOfTracks() const { return (fTracks ? fTracks->GetEntriesFast() : -1); }
  Int_t GetNumberOfMCTracks() const { return (fMCTracks ? fMCTracks->GetEntriesFast() : -1); }
  Int_t GetNumberOfVageMatches() const { return fMultiMatch; }
  Int_t GetNumberOfMatches(ECbmModuleId det) const;
  Int_t GetNumberOfHits(ECbmModuleId det) const;
  Int_t GetNumberOfHitMatches(ECbmModuleId det) const;
  Int_t GetNumberOfPoints(ECbmModuleId det) const;
  Double_t GetEvStartTime() const { return fEvHeader->GetEventTime(); }
 
  TClonesArray* GetHits(ECbmModuleId det) const;
  TClonesArray* GetHitMatches(ECbmModuleId det) const;
  TClonesArray* GetPoints(ECbmModuleId det) const;
 
  TClonesArray* GetCluster(ECbmModuleId det) const;
 
  TObjArray* GetTracks() const { return fTracks; }
  PairAnalysisTrack* GetTrack(UInt_t pos);
 
  virtual void Clear(Option_t* opt = "C");
 
 
protected:
  Int_t GetMvdMatchingIndex(CbmStsTrack* track) const;
 
 
private:
  FairMCEventHeader* fMCHeader = NULL;  // mc header
  FairEventHeader* fEvHeader   = NULL;  // rec header
  TClonesArray* fMCTracks      = NULL;  // mc tracks
  TClonesArray* fStsMatches    = NULL;  // STS matches
  TClonesArray* fMuchMatches   = NULL;  // MUCH matches
  TClonesArray* fTrdMatches    = NULL;  // TRD matches
  TClonesArray* fRichMatches   = NULL;  // RICH matches
 
  TClonesArray* fMvdPoints  = NULL;  // MVD points
  TClonesArray* fStsPoints  = NULL;  // STS points
  TClonesArray* fMuchPoints = NULL;  // MUCH points
  TClonesArray* fRichPoints = NULL;  // RICH points
  TClonesArray* fTrdPoints  = NULL;  // TRD points
  TClonesArray* fTofPoints  = NULL;  // TOF points
 
  TClonesArray* fGlobalTracks = NULL;  // global tracks
  TClonesArray* fTrdTracks    = NULL;  // TRD tracks
  TClonesArray* fStsTracks    = NULL;  // STS tracks
  TClonesArray* fMuchTracks   = NULL;  // MUCH tracks
  TClonesArray* fRichRings    = NULL;  // RICH rings
 
  TClonesArray* fMvdHits       = NULL;  // MVD hits
  TClonesArray* fStsHits       = NULL;  // STS hits
  TClonesArray* fMuchHits      = NULL;  // MUCH pixel hits
  TClonesArray* fMuchHitsStraw = NULL;  // MUCH straw hits
  TClonesArray* fTrdHits       = NULL;  // TRD hits
  TClonesArray* fRichHits      = NULL;  // RICH hits
  TClonesArray* fTofHits       = NULL;  // TOF hits
 
  TClonesArray* fRichProjection = NULL;  // RICH projection
 
  TClonesArray* fMvdHitMatches  = NULL;  // hit matches
  TClonesArray* fStsHitMatches  = NULL;  // hit matches
  TClonesArray* fRichHitMatches = NULL;  // hit matches
  TClonesArray* fMuchHitMatches = NULL;  // hit matches
  TClonesArray* fTrdHitMatches  = NULL;  // hit matches
  TClonesArray* fTofHitMatches  = NULL;  // hit matches
 
  TClonesArray* fTrdCluster = NULL;  // TRD cluster
 
  CbmVertex* fPrimVertex    = NULL;  // primary vertex
  TClonesArray* fFastTracks = NULL;  // fast(sim) tracks
  TObjArray* fTracks;                // papa tracks
 
  Int_t fMultiMatch = 0;  // number of mutiple matched tracks
  Bool_t fANN       = true;
  Int_t fANNId      = 1;
 
  PairAnalysisEvent(const PairAnalysisEvent& event);
  PairAnalysisEvent& operator=(const PairAnalysisEvent& c);
 
  ClassDef(PairAnalysisEvent, 7)  // Event structure
};
 
 
#endif