#include <fstream>
#include "TSystem.h"
#include "TROOT.h"
#include "TRandom.h"
#include "TMath.h"
#include "TBits.h"
#include "TFormula.h"
#include "TF1.h"
#include "TLinearFitter.h"
#include "TFile.h"
#include "TChain.h"
#include "TCut.h"
#include "TEntryList.h"
#include "TH1F.h"
#include "THnSparse.h"
#include "AliSysInfo.h"
#include "AliExternalTrackParam.h"
#include "TTreeStream.h"
#include "AliTrackPointArray.h"
#include "AliLog.h"
#include "AliTPCTransformation.h"
#include "AliTPCkalmanFit.h"
#include "AliMathBase.h"
#include "AliXRDPROOFtoolkit.h"

Include dependency graph for CalibAlignKalman.C:

Functions
TMatrixD	cutMatrix (4 *7, 2)

AliTPCkalmanFit *	CalibAlignKalmanFit (Int_t maxTracks, Int_t trackDump)

void	FilterTracks ()

AliTPCkalmanFit *	SetupFit ()

void	AddFitFieldCage (AliTPCkalmanFit *kalmanFit)

void	AddPhiScaling (AliTPCkalmanFit *kalmanFit)

void	AddZShift (AliTPCkalmanFit *kalmanFit, Int_t ncos, Int_t nsin)

void	AddZTilting (AliTPCkalmanFit *kalmanFit, Int_t ncos, Int_t nsin)

void	AddLocalXYMisalignment (AliTPCkalmanFit *kalmanFit)

void	AddLocalXYMisalignmentSector (AliTPCkalmanFit *kalmanFit)

void	AddAlignOROCIROCFourier (AliTPCkalmanFit *kalmanFit, Int_t ncos, Int_t nsin)

void	AddAlignSectorFourier (AliTPCkalmanFit *kalmanFit, Int_t ncos, Int_t nsin)

void	AddDrift (AliTPCkalmanFit *kalmanFit)

AliTrackPointArray *	FilterPoints (AliTrackPointArray &points, Int_t dir, TTreeSRedirector *pcstream)

TVectorD *	EstimateScatteringKalmanLinear (AliTrackPointArray &points, AliExternalTrackParam &p0, AliExternalTrackParam &p1, TTreeSRedirector *pcstream)

AliTrackPointArray *	SkipPoints (AliTrackPointArray &points, Int_t nskip, Int_t nskipOffset)

AliTPCkalmanFit *	FitPointsLinear (Int_t maxTracks, Int_t trackDump)

void	CalibAlignKalman (Int_t npoints, Int_t maxFiles, Int_t startFile, Int_t trackDump, Int_t nSkipTrack, Int_t nSkipTrackOffset, Int_t nSkip, Int_t nSkipOffset, Int_t bfilterTest)

void	QAPointsLinear (Int_t maxTracks, Int_t trackDump)

void	TestScattering (Int_t maxTracks, Int_t trackDump)

AliTrackPointArray *	SortPoints (AliTrackPointArray &points)

void	SelectPixel ()

void	SelectNonPixelA ()

void	SelectNonPixelC ()

void	DumpQA1D (TObjArray &arrayOut)

void	MakeFolders (TObjArray *arrayOut)

Variables
const Float_t	krmsYcut = 0.2

const Float_t	krmsZcut = 0.2

const Float_t	krmsYcutGlobal = 0.2

const Float_t	krmsZcutGlobal = 2.0

const Float_t	kSigmaCut = 5.

const Int_t	knclCut = 80

const Double_t	kArmCut = 50.

const Double_t	kNsigma = 5.

const Double_t	rmsCut09 [6] ={0.15,0.9,0.015, 0.2, 0.9, 0.015}

Int_t	toSkip = 2

Int_t	toSkipOffset = 0

Int_t	toSkipTrack = 2

Int_t	toSkipTrackOffset = 0

Int_t	isFilterTest = 0

TCut *	cSide [4] ={0,0,0,0}

TCut *	cP3 [4] ={0,0,0,0}

TCut *	cSP3 [4] ={0,0,0,0}

TCut *	cP4 [4] ={0,0,0,0}

TCut *	cM4 [4] ={0,0,0,0}

TCut *	cA [4] ={0,0,0,0}

TCut	cutAll =""

TChain *	chainPoints =0

TEntryList *	elist =0

AliTPCkalmanFit *	kalmanFitNew =0

AliTPCkalmanFit *	kalmanFitOrig =0

AliTPCkalmanFit *	kalmanFitApply =0

Detailed Description

gSystem->AddIncludePath("-I$ALICE_ROOT/TPC");
gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/macros");
gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+");
gROOT->LoadMacro("$ALICE_ROOT/TPC/CalibMacros/CalibAlignKalman.C+");
AliTPCTransformation::BuildBasicFormulas();
AliXRDPROOFtoolkit tool;
chainPoints = tool.MakeChainRandom("align.txt","trackPoints",0,50000);
chainPoints->Lookup();
chainMS = tool.MakeChainRandom("kalmanFit.list","kf",0,50000);
chainMS->Lookup();
chainFP = tool.MakeChainRandom("kalmanFit.list","filter",0,50000);
chainFP->Lookup();
CalibAlignKalmanFit(40000,1);
kalmanFit0->DumpCorelation(0.8);
TFile f("kalmanfitTPC.root");

Definition in file CalibAlignKalman.C.

Function Documentation

void AddAlignOROCIROCFourier	(	AliTPCkalmanFit *	kalmanFit,
		Int_t	ncos,
		Int_t	nsin
	)

Definition at line 1182 of file CalibAlignKalman.C.

Referenced by SetupFit().

void AddAlignSectorFourier	(	AliTPCkalmanFit *	kalmanFit,
		Int_t	ncos,
		Int_t	nsin
	)

Definition at line 1275 of file CalibAlignKalman.C.

Referenced by SetupFit().

void AddDrift ( AliTPCkalmanFit * kalmanFit )

Add drift velocity transformation

Definition at line 934 of file CalibAlignKalman.C.

Referenced by SetupFit().

void AddFitFieldCage ( AliTPCkalmanFit * kalmanFit )

Add radial scaling due field cage

Definition at line 858 of file CalibAlignKalman.C.

Referenced by SetupFit().

void AddLocalXYMisalignment ( AliTPCkalmanFit * kalmanFit )

Definition at line 1096 of file CalibAlignKalman.C.

void AddLocalXYMisalignmentSector ( AliTPCkalmanFit * kalmanFit )

Definition at line 1126 of file CalibAlignKalman.C.

void AddPhiScaling ( AliTPCkalmanFit * kalmanFit )

Add linear local phi scaling - separate IROC/OROC - A side/C side "tscalingLocalPhiIROC" "tscalingLocalPhiOROC"

Definition at line 906 of file CalibAlignKalman.C.

Referenced by SetupFit().

void AddZShift	(	AliTPCkalmanFit *	kalmanFit,
		Int_t	ncos,
		Int_t	nsin
	)

Definition at line 952 of file CalibAlignKalman.C.

Referenced by SetupFit().

void AddZTilting	(	AliTPCkalmanFit *	kalmanFit,
		Int_t	ncos,
		Int_t	nsin
	)

z tilting absolute (sector) and relative (IROC-OROC)

Definition at line 1025 of file CalibAlignKalman.C.

Referenced by SetupFit().

void CalibAlignKalman	(	Int_t	npoints,
		Int_t	maxFiles,
		Int_t	startFile,
		Int_t	trackDump,
		Int_t	nSkipTrack,
		Int_t	nSkipTrackOffset,
		Int_t	nSkip,
		Int_t	nSkipOffset,
		Int_t	bfilterTest
	)

Definition at line 138 of file CalibAlignKalman.C.

AliTPCkalmanFit * CalibAlignKalmanFit	(	Int_t	maxTracks,
		Int_t	trackDump
	)

Fitting procedure

Definition at line 172 of file CalibAlignKalman.C.

Referenced by CalibAlignKalman().

TMatrixD cutMatrix	(	4 *	7,
		2
	)

Referenced by FilterTracks(), and FitPointsLinear().

void DumpQA1D ( TObjArray & arrayOut )

Definition at line 1412 of file CalibAlignKalman.C.

TVectorD * EstimateScatteringKalmanLinear	(	AliTrackPointArray &	points,
		AliExternalTrackParam &	p0,
		AliExternalTrackParam &	p1,
		TTreeSRedirector *	pcstream
	)