EcalDarkBremFilter.h

Go to the documentation of this file.

 
#ifndef BIASING_ECALDARKBREMFILTER_H_
#define BIASING_ECALDARKBREMFILTER_H_
 
//----------------//
//   C++ StdLib   //
//----------------//
#include <algorithm>
 
//------------//
//   Geant4   //
//------------//
#include "G4RunManager.hh"
 
/*~~~~~~~~~~~~*/
/*   SimCore  */
/*~~~~~~~~~~~~*/
 
#include "G4DarkBreM/G4APrime.h"              //checking if particles match A'
#include "G4DarkBreM/G4DarkBremsstrahlung.h"  //checking for dark brem secondaries
#include "G4LogicalVolumeStore.hh"            //for the store
#include "SimCore/G4User/UserAction.h"
#include "SimCore/G4User/UserTrackInformation.h"  //make sure A' is saved
 
namespace biasing {
 
class EcalDarkBremFilter : public simcore::UserAction {
 public:
  EcalDarkBremFilter(const std::string& name,
                     framework::config::Parameters& parameters);
 
  ~EcalDarkBremFilter() {}
 
  std::vector<simcore::TYPE> getTypes() override {
    return {simcore::TYPE::STACKING, simcore::TYPE::EVENT,
            simcore::TYPE::TRACKING};
  }
 
  void BeginOfEventAction(const G4Event* event) override;
 
  G4ClassificationOfNewTrack ClassifyNewTrack(
      const G4Track* aTrack,
      const G4ClassificationOfNewTrack& currentTrackClass) override;
 
  void NewStage() override;
 
  void PostUserTrackingAction(const G4Track* track) override;
 
 private:
  bool inDesiredVolume(const G4Track*) const;
 
  void AbortEvent(const std::string& reason) const;
 
 private:
  double threshold_;
 
  std::vector<G4LogicalVolume*> volumes_;
 
  bool found_ap_;
 
};  // EcalDarkBremFilter
}  // namespace biasing
 
#endif  // BIASING_ECALDARKBREMFILTER_H__