QAmacro_sim_G4.C File Reference

Go to the source code of this file.

Functions
timer	Start ()
fRun	SetName ("TGeant4")
fRun	SetOutputFile ("simfullG4.root")
fRun	SetMaterials ("media_pnd.geo")
Cave	SetGeometryFileName ("pndcave.geo")
fRun	AddModule (Cave)
Magnet	SetGeometryFileName ("FullSolenoid.root")
fRun	AddModule (Magnet)
Dipole	SetGeometryFileName ("dipole.geo")
fRun	AddModule (Dipole)
Pipe	SetGeometryFileName ("pipebeamtarget.geo")
fRun	AddModule (Pipe)
Tpc	SetGeometryFileName ("tpc.geo")
fRun	AddModule (Tpc)
Mvd	SetGeometryFileName ("MVD_v1.0_woPassiveTraps.root")
fRun	AddModule (Mvd)
Emc	SetGeometryFileNameDouble ("emc_module1245.dat","emc_module3new.root")
fRun	AddModule (Emc)
Tof	SetGeometryFileName ("tofbarrel.geo")
fRun	AddModule (Tof)
Drc	SetGeometryFileName ("dirc.geo")
fRun	AddModule (Drc)
Muo	SetBarrel ("fast")
Muo	SetEndcap ("fast")
Muo	SetMuonFilter ("fast")
Muo	SetMdtMagnet (kTRUE)
Muo	SetMdtMFIron (kTRUE)
fRun	AddModule (Muo)
Dch	SetGeometryFileName ("dch.root")
fRun	AddModule (Dch)
fRun	SetGenerator (primGen)
gRandom	SetSeed (3523)
boxGen	SetPRange (1., 10.)
boxGen	SetPhiRange (-180, 180)
boxGen	SetThetaRange (1, 5)
boxGen	SetXYZ (0., 0., 0.)
primGen	AddGenerator (boxGen)
boxGen1	SetPRange (1., 1.)
boxGen1	SetPhiRange (0., 360.)
boxGen1	SetThetaRange (0., 180.)
primGen	AddGenerator (boxGen1)
boxGen2	SetPRange (2.192840326, 2.192840326)
primGen	AddGenerator (boxGen2)
fBox	SetPRange (0.1, 2.5)
fBox	SetPhiRange (0., 360)
fBox	SetThetaRange (5, 150)
primGen	AddGenerator (fBox)
fRun	SetBeamMom (15)
fRun	SetField (fField)
fRun	Init ()
	if (fField)
Par	setInputVersion (fRun->GetRunId(), 1)
Par	setChanged ()
output	open ("g4params.root")
rtdb	setOutput (output)
rtdb	saveOutput ()
rtdb	print ()
fRun	Run (nEvents)
timer	Stop ()
	printf ("RealTime=%f seconds, CpuTime=%f seconds\n", rtime, ctime)

Variables
	gDebug =0
FairRunSim *	fRun = new FairRunSim()
FairModule *	Cave = new PndCave("CAVE")
FairModule *	Magnet = new PndMagnet("MAGNET")
FairModule *	Dipole = new PndMagnet("MAGNET")
FairModule *	Pipe = new PndPipe("PIPE")
FairDetector *	Tpc = new PndTpcDetector("TPC", kTRUE)
FairDetector *	Mvd = new PndMvdDetector("MVD", kTRUE)
PndEmc *	Emc = new PndEmc("EMC",kTRUE)
FairDetector *	Tof = new PndTof("TOF",kTRUE)
FairDetector *	Drc = new PndDrc("DIRC", kTRUE)
PndMdt *	Muo = new PndMdt("MDT",kTRUE)
FairDetector *	Dch = new PndDchDetector("DCH", kTRUE)
FairPrimaryGenerator *	primGen = new FairPrimaryGenerator()
FairBoxGenerator *	boxGen = new FairBoxGenerator(2212, 3)
FairBoxGenerator *	boxGen1 = new FairBoxGenerator(22, 10)
FairBoxGenerator *	boxGen2 = new FairBoxGenerator(-211, 5)
FairBoxGenerator *	fBox = new FairBoxGenerator(211, 5)
PndMultiField *	fField = new PndMultiField("FULL")
FairRuntimeDb *	rtdb =fRun->GetRuntimeDb()
Bool_t	kParameterMerged =kTRUE
PndMultiFieldPar *	Par = (PndMultiFieldPar*) rtdb->getContainer("PndMultiFieldPar")
FairParRootFileIo *	output =new FairParRootFileIo(kParameterMerged)
Int_t	nEvents = 1
Double_t	rtime = timer.RealTime()
Double_t	ctime = timer.CpuTime()

Function Documentation

primGen AddGenerator

(

boxGen

)

primGen AddGenerator

(

boxGen1

)

primGen AddGenerator

(

boxGen2

)

primGen AddGenerator

(

fBox

)

fRun AddModule

(

Cave

)

fRun AddModule

(

Magnet

)

fRun AddModule

(

Dipole

)

fRun AddModule

(

Pipe

)

fRun AddModule

(

Tpc

)

fRun AddModule

(

Mvd

)

fRun AddModule

(

Emc

)

fRun AddModule

(

Tof

)

fRun AddModule

(

Drc

)

fRun AddModule

(

Muo

)

fRun AddModule

(

Dch

)

if

(

fField

)

Definition at line 118 of file QAmacro_sim_G4.C.

References PndMultiFieldPar::SetParameters().

{  Par->SetParameters(fField); }

fRun Init

(

)

Initialize the session

Definition at line 20 of file NHitsPerEvent.C.

{
        TFile* f = new TFile("2014-10-27-00-06-02_50_MHz_Filtered_PosCorrected_Event.root", "READ");
        t = (TTree*)f->Get("pndsim");
//      t->AddFriend("pndsim", "2014-10-27-00-06-02_50_MHz_Filtered_PosCorrected_Residuals.root");
//      t->AddFriend("pndsim", "2014-10-27-00-06-02_50_MHz_Filtered_PosCorrected_Track.root");
//      t->AddFriend("pndsim", "2014-10-27-00-06-02_50_MHz_Filtered_TimeCorrectedSorted_Event.root");
        t->StartViewer();

        t->SetBranchAddress("MVDHitsPixelPosCorrected", &mvdHitArray);
//      t->SetBranchAddress("MvdProjectedHit", &mvdProjectedArray);
//      t->SetBranchAddress("MvdTrack", &trackArray);

        gHitPoints->SetLineColor(kRed);
        gHitPoints->SetMarkerStyle(2);
        gHitPoints->SetMarkerColor(kRed);

        gProjections->SetLineColor(kBlue);
        gProjections->SetMarkerStyle(5);
        gProjections->SetMarkerColor(kBlue);


//      mg->Add(gProjections, "Projected");
        mg->Add(gHitPoints, "Hits");
        t->GetEntry(0);
}

output open

(

"g4params.root"

)

rtdb print

(

)

printf

(

"

RealTime = %f seconds

)

fRun Run

(

nEvents

)

rtdb saveOutput

(

)

Muo SetBarrel

(

"fast"

)

fRun SetBeamMom

(

15

)

Set the field

Par setChanged

(

)

Muo SetEndcap

(

"fast"

)

fRun SetField

(

fField

)

fRun SetGenerator

(

primGen

)

Cave SetGeometryFileName

(

"pndcave.geo"

)

Magnet SetGeometryFileName

(

"FullSolenoid.root"

)

Dipole SetGeometryFileName

(

"dipole.geo"

)

Pipe SetGeometryFileName

(

"pipebeamtarget.geo"

)

Tpc SetGeometryFileName

(

"tpc.geo"

)

Mvd SetGeometryFileName

(

"MVD_v1.0_woPassiveTraps.root"

)

Tof SetGeometryFileName

(

"tofbarrel.geo"

)

Drc SetGeometryFileName

(

"dirc.geo"

)

Dch SetGeometryFileName

(

"dch.root"

)

Emc SetGeometryFileNameDouble	(	"emc_module1245.dat"	,
		"emc_module3new.root"
	)

Par setInputVersion	(	fRun->	GetRunId(),
		1
	)

fRun SetMaterials

(

"media_pnd.geo"

)

Muo SetMdtMagnet

(

kTRUE

)

Muo SetMdtMFIron

(

kTRUE

)

Muo SetMuonFilter

(

"fast"

)

fRun SetName

(

"TGeant4"

)

rtdb setOutput

(

output

)

fRun SetOutputFile

(

"simfullG4.root"

)

boxGen SetPhiRange	(	-	180,
		180
	)

boxGen1 SetPhiRange	(	0.	,
		360.
	)

fBox SetPhiRange	(	0.	,
		360
	)

boxGen SetPRange	(	1.	,
		10.
	)

boxGen1 SetPRange	(	1.	,
		1.
	)

boxGen2 SetPRange	(	2.	192840326,
		2.	192840326
	)

fBox SetPRange	(	0.	1,
		2.	5
	)

gRandom SetSeed

(

3523

)

boxGen SetThetaRange	(	1	,
		5
	)

boxGen1 SetThetaRange	(	0.	,
		180.
	)

fBox SetThetaRange	(	5	,
		150
	)

boxGen SetXYZ	(	0.	,
		0.	,
		0.
	)

timer Start

(

)

timer Stop

(

)

Variable Documentation

FairBoxGenerator* boxGen = new FairBoxGenerator(2212, 3)

Definition at line 75 of file QAmacro_sim_G4.C.

FairBoxGenerator* boxGen1 = new FairBoxGenerator(22, 10)

Definition at line 82 of file QAmacro_sim_G4.C.

FairBoxGenerator* boxGen2 = new FairBoxGenerator(-211, 5)

Definition at line 89 of file QAmacro_sim_G4.C.

FairModule* Cave = new PndCave("CAVE")

Definition at line 19 of file QAmacro_sim_G4.C.

Double_t ctime = timer.CpuTime()

Definition at line 137 of file QAmacro_sim_G4.C.

FairDetector* Dch = new PndDchDetector("DCH", kTRUE)

Definition at line 64 of file QAmacro_sim_G4.C.

FairModule* Dipole = new PndMagnet("MAGNET")

Definition at line 27 of file QAmacro_sim_G4.C.

FairDetector* Drc = new PndDrc("DIRC", kTRUE)

Definition at line 51 of file QAmacro_sim_G4.C.

PndEmc* Emc = new PndEmc("EMC",kTRUE)

Definition at line 43 of file QAmacro_sim_G4.C.

FairBoxGenerator* fBox = new FairBoxGenerator(211, 5)

Definition at line 94 of file QAmacro_sim_G4.C.

PndMultiField* fField = new PndMultiField("FULL")

Definition at line 106 of file QAmacro_sim_G4.C.

FairRunSim* fRun = new FairRunSim()

Definition at line 7 of file QAmacro_sim_G4.C.

gDebug =0

Definition at line 6 of file QAmacro_sim_G4.C.

Bool_t kParameterMerged =kTRUE

Definition at line 116 of file QAmacro_sim_G4.C.

FairModule* Magnet = new PndMagnet("MAGNET")

Definition at line 23 of file QAmacro_sim_G4.C.

PndMdt* Muo = new PndMdt("MDT",kTRUE)

Definition at line 55 of file QAmacro_sim_G4.C.

FairDetector* Mvd = new PndMvdDetector("MVD", kTRUE)

Definition at line 39 of file QAmacro_sim_G4.C.

Int_t nEvents = 1

Definition at line 131 of file QAmacro_sim_G4.C.

FairParRootFileIo* output =new FairParRootFileIo(kParameterMerged)

Definition at line 122 of file QAmacro_sim_G4.C.

PndMultiFieldPar* Par = (PndMultiFieldPar*) rtdb->getContainer("PndMultiFieldPar")

Definition at line 117 of file QAmacro_sim_G4.C.

FairModule* Pipe = new PndPipe("PIPE")

Definition at line 31 of file QAmacro_sim_G4.C.

FairPrimaryGenerator* primGen = new FairPrimaryGenerator()

Create and Set Event Generator

Definition at line 71 of file QAmacro_sim_G4.C.

FairRuntimeDb* rtdb =fRun->GetRuntimeDb()

Definition at line 115 of file QAmacro_sim_G4.C.

Double_t rtime = timer.RealTime()

Definition at line 136 of file QAmacro_sim_G4.C.

FairDetector* Tof = new PndTof("TOF",kTRUE)

Definition at line 47 of file QAmacro_sim_G4.C.

FairDetector* Tpc = new PndTpcDetector("TPC", kTRUE)

Definition at line 35 of file QAmacro_sim_G4.C.