eta_c/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 (1)
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)
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

)

corrG SetMatchDistance

(

0.

18

)

trackInit SetMCPid

(

)

corrG SetMinGEMHits

(

2

)

corr SetMinMVDHits

(

1

)

kalman SetNumIterations

(

3

)

kalman2 SetOutBranchName

(

"TrackPostFitMVD"

)

kalman3 SetOutBranchName

(

"TrackPostFitComplete"

)

fRun SetOutputFile

(

outFile

)

corr SetOutTrackBranchName

(

"TrackPreFitMVD"

)

corrG SetOutTrackBranchName

(

"TrackPreFitGEM"

)

trackInit SetPersistence

(

)

corrG 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 eta_c/run_reco_tpc_evt.C.

PndGFTrackToPndTrackConvertorTask* converter =new PndGFTrackToPndTrackConvertorTask()

Definition at line 118 of file eta_c/run_reco_tpc_evt.C.

PndTpcMVDCorrelatorTask* corr = new PndTpcMVDCorrelatorTask()

Definition at line 86 of file eta_c/run_reco_tpc_evt.C.

PndTpcGEMCorrelatorTask* corrG = new PndTpcGEMCorrelatorTask()

Definition at line 102 of file eta_c/run_reco_tpc_evt.C.

Double_t ctime = timer.CpuTime()

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

FairRunAna* fRun = new FairRunAna()

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

FairGeane* Geane = new FairGeane()

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

PndGeoHandling* geoH = PndGeoHandling::Instance()

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

TString inDigiFile = "evt_digi_tpc.root"

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

TString inSimFile = "evt_points_tpc.root"

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

KalmanTask* kalman =new KalmanTask()

Definition at line 80 of file eta_c/run_reco_tpc_evt.C.

Referenced by DecayTreeFitter::Constraint::filter(), FairGeaneTrKalStt::Kalman(), RecoComplete(), and PndSttMvdGemTracking::Retrack().

KalmanTask* kalman2 =new KalmanTask()

Definition at line 95 of file eta_c/run_reco_tpc_evt.C.

KalmanTask* kalman3 =new KalmanTask()

Definition at line 111 of file eta_c/run_reco_tpc_evt.C.

Int_t nEvents = 0

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

TString outFile = "evt_reco_tpc.root"

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

TString parFile = "evt_params_tpc.root"

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

FairParRootFileIo* parInput1 = new FairParRootFileIo()

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

FairParAsciiFileIo* parIo1 = new FairParAsciiFileIo()

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

FairRuntimeDb* rtdb = fRun->GetRuntimeDb()

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

Double_t rtime = timer.RealTime()

Definition at line 144 of file eta_c/run_reco_tpc_evt.C.

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

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

TStopwatch timer

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

PndTpcRiemannTrackingTask* tpcSPR = new PndTpcRiemannTrackingTask()

Definition at line 65 of file eta_c/run_reco_tpc_evt.C.

PndTpcTrackInitTask* trackInit =new PndTpcTrackInitTask()

Definition at line 71 of file eta_c/run_reco_tpc_evt.C.

PndMCTrackAssociator* trackMC = new PndMCTrackAssociator()

Definition at line 123 of file eta_c/run_reco_tpc_evt.C.