first commit

Author: xiangliu535

Makefile

@@ -0,0 +1,27 @@
+# --------------------------------------------------------------
+# GNUmakefile
+# --------------------------------------------------------------
+
+name := minidex_sim
+G4TARGET := $(name)
+G4EXLIB := true
+
+G4DEBUG := 0
+
+ROOTCFLAGS      = $(shell root-config --cflags)
+ROOTLIBS        = $(shell root-config --libs)
+ROOTGLIBS       = $(shell root-config --glibs)
+
+EXTRALIBS +=$(ROOTLIBS)
+EXTRALIBS +=$(ROOTGLIBS)
+EXTRALIBS +=-L/usr/lib
+
+
+CPPFLAGS += $(ROOTCFLAGS)
+
+.PHONY: all
+all: lib bin
+print-%  : ; @echo $* = $($*)
+
+include $(G4INSTALL)/config/binmake.gmk
+

include/HPGeAnalysisManager.hh

@@ -0,0 +1,68 @@
+#ifndef __XENON10PANALYSISMANAGER_H__
+#define __XENON10PANALYSISMANAGER_H__
+
+#include <globals.hh>
+ #include "G4PhysicalConstants.hh"
+ #include "G4SystemOfUnits.hh"
+#include <G4Navigator.hh>
+ 
+#include <TParameter.h>
+
+class G4Run;
+class G4Event;
+class G4Step;
+
+class TFile;
+class TTree;
+
+class HPGeEventData;
+class HPGePrimaryGeneratorAction;
+
+class HPGeAnalysisManager
+{
+public:
+	HPGeAnalysisManager(HPGePrimaryGeneratorAction *pPrimaryGeneratorAction);
+	virtual ~HPGeAnalysisManager();
+
+public:
+	virtual void BeginOfRun(const G4Run *pRun); 
+	virtual void EndOfRun(const G4Run *pRun); 
+	virtual void BeginOfEvent(const G4Event *pEvent); 
+	virtual void EndOfEvent(const G4Event *pEvent); 
+	virtual void Step(const G4Step *pStep);	
+
+	void SetDataFilename(const G4String &hFilename) { m_hDataFilename = hFilename; }
+	void SetNbEventsToSimulate(G4int iNbEventsToSimulate) { m_iNbEventsToSimulate = iNbEventsToSimulate; }
+
+    void FillHitsToNtuple(const G4Event *pEvent);
+    void WritePartialEvent(const G4Event *pEvent);
+    void ResetPartialEvent(const G4Event *pEvent); // remove all hits collections
+    //void SaveBremsstrahlungGamma(G4float px, G4float py, G4float pz);
+
+        bool GetNeedTrajectoryOrNot() {return NeedTrajectory;}
+        bool GetNeedTrajectoryPointOrNot() {return NeedTrajectoryPoint;}
+
+private:
+	//G4bool FilterEvent(HPGeEventData *pEventData);
+
+private:
+	G4int m_iHPGeHitsCollectionID;
+
+	G4String m_hDataFilename;
+	G4int m_iNbEventsToSimulate;
+
+	TFile *m_pTreeFile;
+	TTree *m_pTree;
+	TParameter<int> *m_pNbEventsToSimulateParameter;
+
+	HPGePrimaryGeneratorAction *m_pPrimaryGeneratorAction;
+
+	HPGeEventData *m_pEventData;
+
+    bool NeedTrajectory;
+    bool NeedTrajectoryPoint;
+
+};
+
+#endif // __XENON10PANALYSISMANAGER_H__
+

include/HPGeDetectorConstruction.hh

@@ -0,0 +1,66 @@
+#ifndef __XENON10PDETECTORCONSTRUCTION_H__
+#define __XENON10PDETECTORCONSTRUCTION_H__
+
+//----------------------------------------------------------
+// from outside inwards
+//   20cm Pb, 10cmCu, inner chamber 20cmx20cm space, 40cm high
+//   HPGe at the very origin
+//   all these logical volumes should be put directly into Mother volume
+//   therefore, all shield volumes are subtracted volumes
+//
+//  2011.8.7 move the detector 9mm away from the copper wall and add PE hat
+//           following 3 components are moved towards +x direction by 9 mm within the Chamber
+//           Cryostat
+//           ChamberCoolingFinger
+//           CryostatinnerVolume 
+//
+// 2013.5.18 construct GEMMX detector with modified shielding
+//           the counting chamber volume is 20cm x 20cm x 35cm
+//           the copper shield is 40cm x 40cm, total height is 70cm, top 10cm thick, bottom 25cm thick,
+//           the lead shield is 80cm x 80cm, total height is 110cm,
+//
+// 2013.6.20 put eight blocks of stainlesssteel samples ,3.844kg
+//           100x100x50 mm3
+//----------------------------------------------------------
+
+#include <globals.hh>
+
+#include <vector>
+#include <map>
+
+using std::vector;
+using std::map;
+
+class G4LogicalVolume;
+class G4VPhysicalVolume;
+class HPGeDetectorMessenger;
+
+#include <G4VUserDetectorConstruction.hh>
+
+class HPGeDetectorConstruction: public G4VUserDetectorConstruction
+{
+public:
+	HPGeDetectorConstruction();
+	~HPGeDetectorConstruction();
+
+	G4VPhysicalVolume* Construct();
+
+        void setGDMLsourceFile( G4String );
+
+private:
+    G4LogicalVolume *m_pMotherLogicalVolume;
+    G4LogicalVolume *m_pLabLogicalVolume;
+    G4VPhysicalVolume *m_pMotherPhysicalVolume;
+    G4VPhysicalVolume *m_pLabPhysicalVolume;
+
+	void CheckOverlapping();
+	void PrintGeometryInformation();
+
+private:
+	
+	HPGeDetectorMessenger *m_pDetectorMessenger;
+        G4String fGDMLsourceFile;
+};
+
+#endif // __XENON10PDETECTORCONSTRUCTION_H__
+

include/HPGeDetectorMessenger.hh

@@ -0,0 +1,38 @@
+#ifndef HPGeDetectorMessenger_h
+#define HPGeDetectorMessenger_h 1
+
+#include <G4UImessenger.hh>
+#include <globals.hh>
+
+class HPGeDetectorConstruction;
+
+class G4UIcommand;
+class G4UIdirectory;
+class G4UIcmdWithoutParameter;
+class G4UIcmdWithAString;
+class G4UIcmdWithADoubleAndUnit;
+class G4UIcmdWith3Vector;
+class G4UIcmdWith3VectorAndUnit;
+class G4UIcmdWithAnInteger;
+class G4UIcmdWithADouble;
+class G4UIcmdWithABool;
+class G4UIcmdWithoutParameter;
+
+class HPGeDetectorMessenger: public G4UImessenger
+{
+public:
+	HPGeDetectorMessenger(HPGeDetectorConstruction *pHPGeDetector);
+	~HPGeDetectorMessenger();
+
+	void SetNewValue(G4UIcommand *pUIcommand, G4String hString);
+
+private:
+	HPGeDetectorConstruction* m_pHPGeDetector;
+
+	G4UIdirectory *m_pDetectorDir;
+        G4UIcmdWithAString *fSetGDMLsourceFileCmd;
+
+};
+
+#endif
+

include/HPGeEventAction.hh

@@ -0,0 +1,25 @@
+#ifndef __XENON10PEVENTACTION_H__
+#define __XENON10PEVENTACTION_H__
+
+#include <G4UserEventAction.hh>
+
+#include "HPGeAnalysisManager.hh"
+
+class G4Event;
+
+class HPGeEventAction : public G4UserEventAction
+{
+public:
+	HPGeEventAction(HPGeAnalysisManager *pAnalysisManager = 0);
+	~HPGeEventAction();
+
+public:
+	void BeginOfEventAction(const G4Event *pEvent);
+	void EndOfEventAction(const G4Event *pEvent);
+
+private:
+	HPGeAnalysisManager *m_pAnalysisManager;
+};
+
+#endif // __XENON10PEVENTACTION_H__
+