RhoGoodPhotonSelector.h

Go to the documentation of this file.

#ifndef RHOGOODPHOTONSELECTOR_H
#define RHOGOODPHOTONSELECTOR_H
//                                                                      //
// RhoGoodPhotonSelector                                                //
//                                                                      //
// Selector classes for photon quality selections                       //
//                                                                      //
// Author List:                                                         //
// Marcel Kunze,  RUB, Feb. 99                                          //
// Copyright (C) 1999-2001, Ruhr-University Bochum.                     //
// Ralf Kliemt, HIM/GSI Feb.2013 (Cleanup & Restructuring)              //
// Klaus Goetzen, GSI, May. 2020 (Updated/modified criteria)            //
//                                                                      //

#include "RhoParticleSelectorBase.h"
#include "PndPidCandidate.h"

class RhoGoodPhotonSelector : public RhoParticleSelectorBase {

 public:
  // Default constructor installs the currently accepted general setting
  RhoGoodPhotonSelector(const char *name = "RhoGoodPhotonSelector", Double_t e=0.0, Int_t nCrys=0, Double_t epx=0., Double_t e25e=0., Double_t z20=0., Double_t z53=1.);
  RhoGoodPhotonSelector(const char *name, criterion crit);
  RhoGoodPhotonSelector(const char *name, const char *crit);

  // Destructor
  virtual ~RhoGoodPhotonSelector() {}

  // operations
  virtual void SetCriterion(const char *crit);
  virtual void SetCriterion(criterion crit);
  virtual Bool_t Accept(RhoCandidate *b);
  virtual Bool_t Accept(PndPidCandidate *b);
  virtual void PrintOn(std::ostream &o = std::cout) const;
  
  void SetCrystals(Int_t min, Int_t max)
  {
    fNcMin = min;
    fNcMax = max;
  }
  
  void SetEnergyPerXtal(Double_t min, Double_t max)
  {
    fEperXtalMin = min;
    fEperXtalMax = max;
  }
  
  void SetEnergy25toEtotal(Double_t min, Double_t max)
  {
    fE25toEMin = min;
    fE25toEMax = max;
  }
    
  void SetZ20(Double_t min, Double_t max)
  {
    fZ20Min = min;
    fZ20Max = max;
  }
  
  void SetZ53(Double_t min, Double_t max)
  {
    fZ53Min = min;
    fZ53Max = max;
  }
  
  void SetTheta(Double_t min, Double_t max)
  {
    fThetaMin = min;
    fThetaMax = max;
  }
 
  void SetShowerEnergy(Double_t min, Double_t max = 1.E8)
  {
    fEmcMin = min;
    fEmcMax = max;
  }
  //    void SetLAT ( Double_t min, Double_t max ) {
  //        fLatMin = min;
  //        fLatMax = max;
  //    }
  //    void SetPhi ( Double_t min, Double_t max ) {
  //        fPhiMin = min;
  //        fPhiMax = max;
  //    }

 private:
  Int_t fNcMin, fNcMax; // Number of crystals
  Double_t fEperXtalMin, fEperXtalMax;  // energy per crystal
  Double_t fE25toEMin, fE25toEMax; // energy in 5x5 array to total energy
  Double_t fZ20Min, fZ20Max; // Zernike moment Z20
  Double_t fZ53Min, fZ53Max; // Zernike moment Z53
  Double_t fThetaMin, fThetaMax; // Angular cuts at poca
  Double_t fEmcMin, fEmcMax;     // In addition: associated bump energy
  // Double_t fLatMin, fLatMax;      // Cut // Angular cuts at poca on LAT
  // Double_t fPhiMin, fPhiMax;          // Angular cuts at poca

 public:
  ClassDef(RhoGoodPhotonSelector, 1) // Track selector
};

// standalone print
std::ostream &operator<<(std::ostream &o, const RhoGoodPhotonSelector &);

#endif