reco_template.C File Reference

Go to the source code of this file.

Functions
gROOT	LoadMacro ("$VMCWORKDIR/gconfig/rootlogon.C")
	rootlogon ()
timer	Start ()
fRun	SetInputFile (inDigiFile)
fRun	AddFriend (inSimFile)
fRun	SetOutputFile (outFile)
parInput1	open (parFile.Data())
parIo1	open (allDigiFile.Data(),"in")
rtdb	setFirstInput (parInput1)
rtdb	setSecondInput (parIo1)
fRun	AddTask (Geane)
fRun	AddTask (tpcSPR)
trackInit	SetPersistence ()
trackInit	SetMCPid ()
trackInit	SetSmoothing (true)
fRun	AddTask (trackInit)
kalman	SetNumIterations (3)
fRun	AddTask (kalman)
corr	SetMatchDistance (0.18)
corr	SetMinMVDHits (3)
corr	SetOutTrackBranchName ("TrackPreFitMVD")
corr	SetPersistence (true)
fRun	AddTask (corr)
kalman2	SetTrackBranchName ("TrackPreFitMVD")
kalman2	SetOutBranchName ("TrackPostFitMVD")
fRun	AddTask (kalman2)
corrG	SetMinGEMHits (2)
corrG	SetTrackBranchName ("TrackPostFitMVD")
corrG	SetOutTrackBranchName ("TrackPreFitGEM")
fRun	AddTask (corrG)
kalman3	SetTrackBranchName ("TrackPreFitGEM")
kalman3	SetOutBranchName ("TrackPostFitComplete")
fRun	AddTask (kalman3)
fRun	Init ()
fRun	Run (0, nEvents)
rtdb	saveOutput ()
rtdb	print ()
corr	WriteHistograms ("MVDRes.root")
corrG	WriteHistograms ("GEMRes.root")
timer	Stop ()

Variables
TString	inDigiFile = "dpm_digi_tpc.root"
TString	inSimFile = "dpm_points_tpc.root"
TString	parFile = "dpm_params_tpc.root"
TString	outFile = "dpm_reco_tpc.root"
Int_t	nEvents = 0
TString	sysFile = gSystem->Getenv("VMCWORKDIR")
TStopwatch	timer
FairRunAna *	fRun = new FairRunAna()
TString	allDigiFile = sysFile+"/macro/params/all.par"
FairRuntimeDb *	rtdb = fRun->GetRuntimeDb()
FairParRootFileIo *	parInput1 = new FairParRootFileIo()
FairParAsciiFileIo *	parIo1 = new FairParAsciiFileIo()
PndGeoHandling *	geoH = PndGeoHandling::Instance()
bool	SimpleClustering =true
FairGeane *	Geane = new FairGeane()
PndTpcRiemannTrackingTask *	tpcSPR = new PndTpcRiemannTrackingTask()
PndTpcTrackInitTask *	trackInit =new PndTpcTrackInitTask()
KalmanTask *	kalman =new KalmanTask()
PndTpcMVDCorrelatorTask *	corr = new PndTpcMVDCorrelatorTask()
KalmanTask *	kalman2 =new KalmanTask()
PndTpcGEMCorrelatorTask *	corrG = new PndTpcGEMCorrelatorTask()
KalmanTask *	kalman3 =new KalmanTask()
Double_t	rtime = timer.RealTime()
Double_t	ctime = timer.CpuTime()

Function Documentation

fRun AddFriend

(

inSimFile

)

fRun AddTask

(

Geane

)

fRun AddTask

(

tpcSPR

)

fRun AddTask

(

trackInit

)

fRun AddTask

(

kalman

)

fRun AddTask

(

corr

)

fRun AddTask

(

kalman2

)

fRun AddTask

(

corrG

)

fRun AddTask

(

kalman3

)

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);
}

gROOT LoadMacro

(

"$VMCWORKDIR/gconfig/rootlogon.C"

)

parInput1 open

(

parFile.

Data()

)

parIo1 open	(	allDigiFile.	Data(),
		"in"
	)

rtdb print

(

)

rootlogon

(

)

Definition at line 1 of file outdated/mpiTools/macros/emc/rootlogon.C.

{
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcDigi");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcMC");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcTools");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcReco");
  gSystem->AddIncludePath("-I$VMCWORKDIR/emc/EmcData");
  gSystem->AddIncludePath("-I$VMCWORKDIR/pnddata");
  gSystem->AddIncludePath("-I$VMCWORKDIR/base");
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
  basiclibs();

  gSystem->Load("libDpmEvtGen.so");
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  gSystem->Load("libTrkBase");
  gSystem->Load("libPndData");
  gSystem->Load("libField");
  gSystem->Load("libPassive");
  gSystem->Load("libGen");
  gSystem->Load("libPGen");
  gSystem->Load("libEmc"); 

  delete gRandom;
  gRandom=new TRandom3(0);
}

fRun Run	(	0	,
		nEvents
	)

rtdb saveOutput

(

)

rtdb setFirstInput

(

parInput1

)

fRun SetInputFile

(

inDigiFile

)

corr SetMatchDistance

(

0.

18

)

trackInit SetMCPid

(

)

corrG SetMinGEMHits

(

2

)

corr SetMinMVDHits

(

3

)

kalman SetNumIterations

(

3

)

kalman2 SetOutBranchName

(

"TrackPostFitMVD"

)

kalman3 SetOutBranchName

(

"TrackPostFitComplete"

)

fRun SetOutputFile

(

outFile

)

corr SetOutTrackBranchName

(

"TrackPreFitMVD"

)

corrG SetOutTrackBranchName

(

"TrackPreFitGEM"

)

trackInit SetPersistence

(

)

corr SetPersistence

(

true

)

rtdb setSecondInput

(

parIo1

)

trackInit SetSmoothing

(

true

)

kalman2 SetTrackBranchName

(

"TrackPreFitMVD"

)

corrG SetTrackBranchName

(

"TrackPostFitMVD"

)

kalman3 SetTrackBranchName

(

"TrackPreFitGEM"

)

timer Start

(

)

timer Stop

(

)

corr WriteHistograms

(

"MVDRes.root"

)

corrG WriteHistograms

(

"GEMRes.root"

)

Variable Documentation

TString allDigiFile = sysFile+"/macro/params/all.par"

Definition at line 41 of file reco_template.C.

PndTpcMVDCorrelatorTask* corr = new PndTpcMVDCorrelatorTask()

Definition at line 94 of file reco_template.C.

PndTpcGEMCorrelatorTask* corrG = new PndTpcGEMCorrelatorTask()

Definition at line 110 of file reco_template.C.

Double_t ctime = timer.CpuTime()

Definition at line 145 of file reco_template.C.

FairRunAna* fRun = new FairRunAna()

Definition at line 33 of file reco_template.C.

FairGeane* Geane = new FairGeane()

Definition at line 60 of file reco_template.C.

PndGeoHandling* geoH = PndGeoHandling::Instance()

Definition at line 52 of file reco_template.C.

TString inDigiFile = "dpm_digi_tpc.root"

Definition at line 7 of file reco_template.C.

TString inSimFile = "dpm_points_tpc.root"

Definition at line 8 of file reco_template.C.

KalmanTask* kalman =new KalmanTask()

Definition at line 88 of file reco_template.C.

KalmanTask* kalman2 =new KalmanTask()

Definition at line 103 of file reco_template.C.

KalmanTask* kalman3 =new KalmanTask()

Definition at line 119 of file reco_template.C.

Int_t nEvents = 0

Definition at line 17 of file reco_template.C.

TString outFile = "dpm_reco_tpc.root"

Definition at line 14 of file reco_template.C.

TString parFile = "dpm_params_tpc.root"

Definition at line 11 of file reco_template.C.

FairParRootFileIo* parInput1 = new FairParRootFileIo()

Definition at line 44 of file reco_template.C.

FairParAsciiFileIo* parIo1 = new FairParAsciiFileIo()

Definition at line 47 of file reco_template.C.

FairRuntimeDb* rtdb = fRun->GetRuntimeDb()

Definition at line 43 of file reco_template.C.

Double_t rtime = timer.RealTime()

Definition at line 144 of file reco_template.C.

bool SimpleClustering =true

Definition at line 54 of file reco_template.C.

TString sysFile = gSystem->Getenv("VMCWORKDIR")

Definition at line 22 of file reco_template.C.

TStopwatch timer

Definition at line 28 of file reco_template.C.

PndTpcRiemannTrackingTask* tpcSPR = new PndTpcRiemannTrackingTask()

Definition at line 73 of file reco_template.C.

PndTpcTrackInitTask* trackInit =new PndTpcTrackInitTask()

Definition at line 79 of file reco_template.C.