MuonSensitiveDetector.h

Go to the documentation of this file.

#ifndef SimG4CMS_Muon_MuonSensitiveDetector_h
#define SimG4CMS_Muon_MuonSensitiveDetector_h

#include "SimG4Core/Notification/interface/Observer.h"
#include "SimG4Core/Notification/interface/BeginOfEvent.h"
#include "SimG4Core/Notification/interface/EndOfEvent.h"
#include "SimG4Core/SensitiveDetector/interface/SensitiveTkDetector.h"
#include "DataFormats/GeometryVector/interface/GlobalPoint.h"
#include "DataFormats/GeometryVector/interface/LocalPoint.h"

#include "G4Step.hh"
#include "G4StepPoint.hh"
#include "G4Track.hh"

#include <string>

class MuonSlaveSD;
class MuonSimHitNumberingScheme;
class MuonFrameRotation;
class UpdatablePSimHit;
class MuonSubDetector;
class MuonG4Numbering;
class SimHitPrinter;
class TrackInformation;
class G4Track;
class G4ProcessTypeEnumerator;
class G4TrackToParticleID;
class SimTrackManager;

class MuonSensitiveDetector : 
public SensitiveTkDetector,
public Observer<const BeginOfEvent*>,
public Observer<const EndOfEvent*>
 {

 public:    
  MuonSensitiveDetector(std::string, const DDCompactView &,
            const SensitiveDetectorCatalog &, edm::ParameterSet const &,
            const SimTrackManager*);
  virtual ~MuonSensitiveDetector();
  virtual G4bool ProcessHits(G4Step *,G4TouchableHistory *);
  virtual uint32_t setDetUnitId(G4Step *);
  virtual void EndOfEvent(G4HCofThisEvent*);

  void fillHits(edm::PSimHitContainer&, std::string use);
  std::vector<std::string> getNames();
  std::string type();

  const MuonSlaveSD* GetSlaveMuon() const {
    return slaveMuon; }
  
 private:
  void update(const BeginOfEvent *);
  void update(const ::EndOfEvent *);
  virtual void clearHits();

  Local3DPoint toOrcaRef(Local3DPoint in ,G4Step * s);
  Local3DPoint toOrcaUnits(Local3DPoint);
  Global3DPoint toOrcaUnits(Global3DPoint);

  TrackInformation* getOrCreateTrackInformation( const G4Track* theTrack );

 private:
  MuonSlaveSD* slaveMuon;
  MuonSimHitNumberingScheme* numbering;
  MuonSubDetector* detector;
  MuonFrameRotation* theRotation;
  MuonG4Numbering* g4numbering;

  void storeVolumeAndTrack(G4Step *);
  bool newHit(G4Step *);
  void createHit(G4Step *);
  void updateHit(G4Step *);
  void saveHit();
  
  Local3DPoint InitialStepPositionVsParent(G4Step * currentStep, G4int levelsUp);
  Local3DPoint FinalStepPositionVsParent(G4Step * currentStep, G4int levelsUp);

  G4VPhysicalVolume * thePV;
  UpdatablePSimHit* theHit;
  uint32_t theDetUnitId; 
  unsigned int theTrackID;
 
  bool printHits;
  SimHitPrinter* thePrinter;
  Global3DPoint theGlobalEntry;

  //--- SimTracks cuts
  double STenergyPersistentCut;
  bool STallMuonsPersistent;

  G4ProcessTypeEnumerator* theG4ProcessTypeEnumerator;

  G4TrackToParticleID* myG4TrackToParticleID;
  const SimTrackManager* theManager;
};

#endif // MuonSensitiveDetector_h