QAmacro_sim_G3.C File Reference

Go to the source code of this file.

Functions
timer	Start ()
fRun	SetName ("TGeant3")
fRun	SetOutputFile ("simfullG3.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)
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)
boxGen2	SetPtRange (0.798217349, 0.798217349)
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 ("g3params.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)
PndMdt *	Muo = new PndMdt("MDT",kTRUE)
FairDetector *	Dch = new PndDchDetector("DCH", kFALSE)
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

(

Muo

)

fRun AddModule

(

Dch

)

if

(

fField

)

Definition at line 116 of file QAmacro_sim_G3.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

(

"g3params.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"

)

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

(

"TGeant3"

)

rtdb setOutput

(

output

)

fRun SetOutputFile

(

"simfullG3.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
	)

boxGen2 SetPtRange	(	0.	798217349,
		0.	798217349
	)

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 73 of file QAmacro_sim_G3.C.

FairBoxGenerator* boxGen1 = new FairBoxGenerator(22, 10)

Definition at line 80 of file QAmacro_sim_G3.C.

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

Definition at line 87 of file QAmacro_sim_G3.C.

FairModule* Cave = new PndCave("CAVE")

Definition at line 19 of file QAmacro_sim_G3.C.

Double_t ctime = timer.CpuTime()

Definition at line 135 of file QAmacro_sim_G3.C.

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

Definition at line 62 of file QAmacro_sim_G3.C.

FairModule* Dipole = new PndMagnet("MAGNET")

Definition at line 27 of file QAmacro_sim_G3.C.

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

Definition at line 43 of file QAmacro_sim_G3.C.

FairBoxGenerator* fBox = new FairBoxGenerator(211, 5)

Definition at line 93 of file QAmacro_sim_G3.C.

PndMultiField* fField = new PndMultiField("FULL")

Definition at line 105 of file QAmacro_sim_G3.C.

FairRunSim* fRun = new FairRunSim()

Definition at line 7 of file QAmacro_sim_G3.C.

gDebug =0

Definition at line 6 of file QAmacro_sim_G3.C.

Bool_t kParameterMerged =kTRUE

Definition at line 114 of file QAmacro_sim_G3.C.

FairModule* Magnet = new PndMagnet("MAGNET")

Definition at line 23 of file QAmacro_sim_G3.C.

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

Definition at line 54 of file QAmacro_sim_G3.C.

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

Definition at line 39 of file QAmacro_sim_G3.C.

Int_t nEvents = 1

Definition at line 129 of file QAmacro_sim_G3.C.

FairParRootFileIo* output =new FairParRootFileIo(kParameterMerged)

Definition at line 120 of file QAmacro_sim_G3.C.

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

Definition at line 115 of file QAmacro_sim_G3.C.

FairModule* Pipe = new PndPipe("PIPE")

Definition at line 31 of file QAmacro_sim_G3.C.

FairPrimaryGenerator* primGen = new FairPrimaryGenerator()

Create and Set Event Generator

Definition at line 69 of file QAmacro_sim_G3.C.

FairRuntimeDb* rtdb =fRun->GetRuntimeDb()

Definition at line 113 of file QAmacro_sim_G3.C.

Double_t rtime = timer.RealTime()

Definition at line 134 of file QAmacro_sim_G3.C.

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

Definition at line 47 of file QAmacro_sim_G3.C.

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

Definition at line 35 of file QAmacro_sim_G3.C.