PndGemFindHits.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

//* $Id: */

// -------------------------------------------------------------------------
// -----                     PndGemFindHits header file                -----
// -----                  Created 15/02/2009 by R. Karabowicz          -----
// -------------------------------------------------------------------------

#ifndef PNDGEMFINDHITS_H
#define PNDGEMFINDHITS_H 1

#include <PndPersistencyTask.h>
#include "TStopwatch.h"
#include "TH2F.h"
#include "TVector3.h"

#include <map>
#include <set>
#include <vector>

class TH1F;
class TClonesArray;
class PndGemMonitor;
class PndGemDigiPar;
class PndGemSensor;
class PndGemStation;
class PndGemCluster;
class PndGemDigitize;

class PndGemFindHits : public PndPersistencyTask {

 public:
  PndGemFindHits();

  PndGemFindHits(Int_t iVerbose);

  PndGemFindHits(const char *name, Int_t iVerbose);

  virtual ~PndGemFindHits();
  void ConfirmHits();   // why in public??
  void ActivateDigis(); // why in public??

  virtual void Exec(Option_t *opt);

  void SetUseClusters(Bool_t bt = kTRUE) { fUseClusters = bt; }

  void RunTimeBased(Bool_t bt = kTRUE) { fTimeOrderedDigi = bt; }

  Bool_t GetTimeBased() const { return fTimeOrderedDigi; }

  void SetHitWindow(Double_t dd = 1.5) { fHitWindow = dd; }

 private:
  PndGemMonitor *fMonitor;                             
  PndGemDigiPar *fDigiPar;                             
  TClonesArray *fDigis;                                
  TClonesArray *fHits;                                 
  TClonesArray *fHitsTemp;                             
  std::map<PndGemSensor *, std::set<Int_t>> fDigiMapF; 
  std::map<PndGemSensor *, std::set<Int_t>> fDigiMapB; 
  Int_t fMCPointBranchId;
  Bool_t fUseClusters;
  Bool_t fTimeOrderedDigi;

  Double_t fPrepTime;
  Double_t fSortTime;
  Double_t fCreateTime;
  Double_t fConfirmTime;
  Double_t fActivateTime;
  Double_t fAllTime;
  Double_t fHitWindow;

  Int_t fTNofEvents;
  Int_t fTNofDigis;
  Int_t fTNofHits;
  Int_t fTNofHitsTemp;

  TStopwatch fTimer;

  virtual void SetParContainers();

  virtual InitStatus Init();

  virtual InitStatus ReInit();

  void MakeSets();

  void SortDigis();

  Int_t FindHits(PndGemSensor *sensor, std::set<Int_t> &fSet, std::set<Int_t> &bSet);

  Int_t FindHits2(PndGemSensor *sensor, std::set<Int_t> &fSet, std::set<Int_t> &bSet);

  void ConfirmHits2();

  virtual void Finish();

  ClassDef(PndGemFindHits, 1);
};

#endif