llbar/run_reco_tpc_evt.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)
fRun	AddTask (Geane)
parInput1	open (parFile.Data())
parIo1	open (allDigiFile.Data(),"in")
rtdb	setFirstInput (parInput1)
rtdb	setSecondInput (parIo1)
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 (0)
corr	SetOutTrackBranchName ("TrackPreFitMVD")
corr	SetPersistence (true)
fRun	AddTask (corr)
kalman2	SetTrackBranchName ("TrackPreFitMVD")
kalman2	SetOutBranchName ("TrackPostFitMVD")
fRun	AddTask (kalman2)
corrG	SetMinGEMHits (0)
corrG	SetTrackBranchName ("TrackPostFitMVD")
corrG	SetOutTrackBranchName ("TrackPreFitGEM")
fRun	AddTask (corrG)
kalman3	SetTrackBranchName ("TrackPreFitGEM")
kalman3	SetOutBranchName ("TrackPostFitComplete")
fRun	AddTask (kalman3)
converter	SetTrackInBranchName ("TrackPostFitComplete")
converter	SetTrackOutBranchName ("PndTrackPostFitComplete")
fRun	AddTask (converter)
trackMC	SetTrackInBranchName ("PndTrackPostFitComplete")
trackMC	SetTrackOutBranchName ("TrackPostFitCompleteID")
fRun	AddTask (trackMC)
fRun	Init ()
fRun	Run (0, nEvents)
rtdb	saveOutput ()
rtdb	print ()
corr	WriteHistograms ("MVDRes.root")
corrG	WriteHistograms ("GEMRes.root")
timer	Stop ()

Variables
TString	inDigiFile = "evt_digi_tpc.root"
TString	inSimFile = "evt_points_tpc.root"
TString	parFile = "evt_params_tpc.root"
TString	outFile = "evt_reco_tpc.root"
Int_t	nEvents = 0
TString	sysFile = gSystem->Getenv("VMCWORKDIR")
TStopwatch	timer
FairRunAna *	fRun = new FairRunAna()
FairGeane *	Geane = new FairGeane()
TString	allDigiFile = sysFile+"/macro/params/all.par"
FairRuntimeDb *	rtdb = fRun->GetRuntimeDb()
FairParRootFileIo *	parInput1 = new FairParRootFileIo()
FairParAsciiFileIo *	parIo1 = new FairParAsciiFileIo()
PndGeoHandling *	geoH = PndGeoHandling::Instance()
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()
PndGFTrackToPndTrackConvertorTask *	converter =new PndGFTrackToPndTrackConvertorTask()
PndMCTrackAssociator *	trackMC = new PndMCTrackAssociator()
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 AddTask

(

converter

)

fRun AddTask

(

trackMC

)

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

(

0

)

corr SetMinMVDHits

(

0

)

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"

)

converter SetTrackInBranchName

(

"TrackPostFitComplete"

)

trackMC SetTrackInBranchName

(

"PndTrackPostFitComplete"

)

converter SetTrackOutBranchName

(

"PndTrackPostFitComplete"

)

trackMC SetTrackOutBranchName

(

"TrackPostFitCompleteID"

)

timer Start

(

)

timer Stop

(

)

corr WriteHistograms

(

"MVDRes.root"

)

corrG WriteHistograms

(

"GEMRes.root"

)

Variable Documentation

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

Definition at line 42 of file llbar/run_reco_tpc_evt.C.

PndGFTrackToPndTrackConvertorTask* converter =new PndGFTrackToPndTrackConvertorTask()

Definition at line 119 of file llbar/run_reco_tpc_evt.C.

PndTpcMVDCorrelatorTask* corr = new PndTpcMVDCorrelatorTask()

Definition at line 87 of file llbar/run_reco_tpc_evt.C.

PndTpcGEMCorrelatorTask* corrG = new PndTpcGEMCorrelatorTask()

Definition at line 103 of file llbar/run_reco_tpc_evt.C.

Double_t ctime = timer.CpuTime()

Definition at line 146 of file llbar/run_reco_tpc_evt.C.

FairRunAna* fRun = new FairRunAna()

Definition at line 33 of file llbar/run_reco_tpc_evt.C.

FairGeane* Geane = new FairGeane()

Definition at line 37 of file llbar/run_reco_tpc_evt.C.

PndGeoHandling* geoH = PndGeoHandling::Instance()

Definition at line 53 of file llbar/run_reco_tpc_evt.C.

TString inDigiFile = "evt_digi_tpc.root"

Definition at line 7 of file llbar/run_reco_tpc_evt.C.

TString inSimFile = "evt_points_tpc.root"

Definition at line 8 of file llbar/run_reco_tpc_evt.C.

KalmanTask* kalman =new KalmanTask()

Definition at line 81 of file llbar/run_reco_tpc_evt.C.

KalmanTask* kalman2 =new KalmanTask()

Definition at line 96 of file llbar/run_reco_tpc_evt.C.

KalmanTask* kalman3 =new KalmanTask()

Definition at line 112 of file llbar/run_reco_tpc_evt.C.

Int_t nEvents = 0

Definition at line 17 of file llbar/run_reco_tpc_evt.C.

TString outFile = "evt_reco_tpc.root"

Definition at line 14 of file llbar/run_reco_tpc_evt.C.

TString parFile = "evt_params_tpc.root"

Definition at line 11 of file llbar/run_reco_tpc_evt.C.

FairParRootFileIo* parInput1 = new FairParRootFileIo()

Definition at line 45 of file llbar/run_reco_tpc_evt.C.

FairParAsciiFileIo* parIo1 = new FairParAsciiFileIo()

Definition at line 48 of file llbar/run_reco_tpc_evt.C.

FairRuntimeDb* rtdb = fRun->GetRuntimeDb()

Definition at line 44 of file llbar/run_reco_tpc_evt.C.

Double_t rtime = timer.RealTime()

Definition at line 145 of file llbar/run_reco_tpc_evt.C.

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

Definition at line 22 of file llbar/run_reco_tpc_evt.C.

TStopwatch timer

Definition at line 28 of file llbar/run_reco_tpc_evt.C.

PndTpcRiemannTrackingTask* tpcSPR = new PndTpcRiemannTrackingTask()

Definition at line 66 of file llbar/run_reco_tpc_evt.C.

PndTpcTrackInitTask* trackInit =new PndTpcTrackInitTask()

Definition at line 72 of file llbar/run_reco_tpc_evt.C.

PndMCTrackAssociator* trackMC = new PndMCTrackAssociator()

Definition at line 124 of file llbar/run_reco_tpc_evt.C.