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

/*
 * PairFinderTask.h
 *
 *  Created on: Jul 22, 2014
 *      Author: R. Klasen, roklasen@uni-mainz.de or r.klasen@gsi.de or r.klasen@ep1.rub.de
 */

#ifndef PNDLMDPAIRFINDERTASK_H_
#define PNDLMDPAIRFINDERTASK_H_

#include "PndLmdAlignStructs.h"
#include <PndLmdHitPair.h>
#include <PndSdsHit.h>
#include "PndSdsTask.h"

#include "TClonesArray.h"

#include <iostream>
#include <string>
#include <vector>

class PndLmdGeometryHelper;

class PndLmdPairFinderTask : public PndSdsTask {

 private:
  // for statistics
  Int_t noOfGoodPairs, noOfEvents;
  Int_t plane0, plane1, plane2, plane3;
  Int_t eventMissedAllPlanes, noOverlap;
  Int_t distanceTooHigh;
  Int_t sumOfPixelHits;
  Int_t unsuitable;
  Int_t hitsClustered, hitsSinglePixel;
  Int_t noOfCombos;

  Double_t _maxDistance;

  Bool_t _ignoreClusters;

  pixelHit getPixelHitFromSdsHit(PndSdsHit *sdsHit);

 public:
  void SetParContainers();

  PndLmdPairFinderTask();
  PndLmdPairFinderTask(const char *name);
  virtual ~PndLmdPairFinderTask();

  virtual void SetBranchNames();

  virtual InitStatus Init();
  virtual InitStatus ReInit();

  virtual void Exec(Option_t *opt);

  virtual void FinishEvent();
  virtual void FinishTask();

  // apply distance cut, will be ignored when using dynamic cut
  void setMaxDistance(Double_t value) { _maxDistance = value; }

  void ignoreClusters(Bool_t value) { _ignoreClusters = value; }

 protected:
  // PndLmdDim* dimension;

  PndLmdGeometryHelper *helper;

  TClonesArray *digiArray;
  TClonesArray *recoArray;
  TClonesArray *hitPairArray;
  TClonesArray *clusterCandidateArray;

  TString fInRecoBranchName;
  TString fInClusterCandidates;

  std::map<int, TClonesArray *> hitPairMap;

  void Register();
  void Reset();

  // function to return result of all checks, distance cut etc.
  void getStatistics(PndLmdHitPair &candidate);

  bool pairDistanceValid(PndLmdHitPair &candidate);

  // self explanatory
  bool candHitsOverlappingArea(const PndLmdHitPair &candidate);

  ClassDef(PndLmdPairFinderTask, 20);
};

#endif /* PAIRFINDERTASK_H_ */