region_cdet.h

Go to the documentation of this file.

/*
 * To change this license header, choose License Headers in Project Properties.
 * To change this template file, choose Tools | Templates
 * and open the template in the editor.
 */
 
/*
 * File:   region_cdet.h
 * Author: dglazier
 *
 * Created on April 27, 2017, 10:01 AM
 */
 
#ifndef REGION_CDET_H
#define REGION_CDET_H
 
#include "clas12defs.h"
#include "region_particle.h"
 
 
namespace clas12 {
 
 
  class region_cdet  : public region_particle {
 
 
  public:
 
    // region_cdet(par_ptr pars,covmat_ptr cm, scint_ptr scp,
    //      trck_ptr trp,traj_ptr trj);
    region_cdet(par_ptr pars,ftbpar_ptr ftbpars,covmat_ptr cm, cal_ptr calp, scint_ptr scp,
        trck_ptr trp,utrck_ptr utrp, traj_ptr trj, cher_ptr chp, ft_ptr ftp,event_ptr event,
        rich_ptr rich,mcpar_ptr mcp=nullptr);
    ~region_cdet() override=default;
 
    
    bool sort() final;
 
    scint_ptr sci(ushort lay) const final;
    traj_ptr traj(ushort det,ushort layer=0) const final;
    trck_ptr trk(ushort lay) const final{_trck->setIndex(_ptrck);return _trck;};
    utrck_ptr utrk(ushort lay) const final{_utrck->setIndex(_ptrck);return _utrck;};
    
     
    
    double getTime() const final{
      if(_ptof>-1){
    _scint->setIndex(_ptof);
    return _scint->getTime();
      }
      _scint->setIndex(_pcnd);
      return _scint->getTime();
      
    }
    double getPath() const final{
      if(_ptof>-1){
    _scint->setIndex(_ptof);
    return _scint->getPath();
      }
      _scint->setIndex(_pcnd);
      return _scint->getPath();
    }
 
    double getDetEnergy() const final{
      //sum of cnd layers
      _scint->setIndex(_pcnd1);
      double energy=_scint->getEnergy();
      _scint->setIndex(_pcnd2);
      energy+=_scint->getEnergy();
      _scint->setIndex(_pcnd3);
      energy+=_scint->getEnergy();
      return energy;
    }
 
    double getDeltaEnergy() const final{
      _scint->setIndex(_ptof);
      return _scint->getEnergy();
    }
 
    short getSector() const  final{
      return 0;
    }
 
    
    private:
 
    
  
    //scintillator indices
    short _pcnd=-1;
    short _pcnd1=-1;
    short _pcnd2=-1;
    short _pcnd3=-1;
    short _ptof=-1;
    //track index
    short _ptrck=-1;
    
   
    
  };
  
  using region_cdet_ptr=clas12::region_cdet*;
  using region_cdet_uptr=std::unique_ptr<clas12::region_cdet>;
 
}
 
#endif /* REGION_CDET_H */