MUONRecoCheck.C File Reference

Utility macro to check the muon reconstruction. More...

#include <Riostream.h>
#include "TMath.h"
#include "TObjArray.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TGraphErrors.h"
#include "TGraphAsymmErrors.h"
#include "TF1.h"
#include "TFile.h"
#include "TCanvas.h"
#include "TLegend.h"
#include "TGeoManager.h"
#include "AliCDBManager.h"
#include "AliGeomManager.h"
#include "AliLog.h"
#include "AliMUONCDB.h"
#include "AliMUONConstants.h"
#include "AliMUONTrack.h"
#include "AliMUONRecoCheck.h"
#include "AliMUONTrackParam.h"
#include "AliMUONRecoParam.h"
#include "AliMUONVTrackStore.h"
#include "AliMUONVCluster.h"
#include "AliMUONTrackExtrap.h"
#include "AliMUONESDInterface.h"
#include "AliMUONVTriggerTrackStore.h"
#include "AliMUONTriggerTrack.h"
#include "AliMpDEIterator.h"

Include dependency graph for MUONRecoCheck.C:

Go to the source code of this file.

Functions
Double_t	langaufun (Double_t *x, Double_t *par)
void	FitGausResVsMom (TH2 *h, Int_t nBins, const Double_t mean0, const Double_t sigma0, const char *fitting, TGraphAsymmErrors *gMean, TGraphAsymmErrors *gSigma)
void	FitPDCAVsMom (TH2 *h, Int_t nBins, const char *fitting, TGraphAsymmErrors *gMean, TGraphAsymmErrors *gSigma)
void	FillResidual (TH1 *h, Int_t i, Double_t &sigma, TGraphErrors *gMean, TGraphErrors *gSigma, Bool_t correctForSystematics, Bool_t fitResiduals)
TCanvas *	DrawVsAng (const char *name, const char *title, TH1 *h1, TH2 *h2)
TCanvas *	DrawVsPos (const char *name, const char *title, TH2 *h1, TH2 *h2, TH2 *h3)
TCanvas *	DrawResMomVsMom (const char *name, const char *title, TH2 *h, Int_t nBins, TF1 *f2=0x0, const char *fitting="")
void	Zoom (TH1 *h, Double_t fractionCut=0.01)
void	MUONRecoCheck (Int_t nEvent=-1, const char *pathSim="./generated/", const char *esdFileName="AliESDs.root", const char *ocdbPath="local://$ALICE_ROOT/OCDB", const Double_t pMin=0., const Double_t pMax=300., const Int_t pNBins=30, Int_t absorberRegion=-1, Bool_t fitResiduals=kTRUE)

Detailed Description

Utility macro to check the muon reconstruction.

Reconstructed tracks are compared to reference tracks. The reference tracks are built from AliTrackReference for the hit in chamber (0..9) and from kinematics (TreeK) for the vertex parameters.

Author

Jean-Pierre Cussonneau, Philippe Pillot, Subatech

Definition in file MUONRecoCheck.C.

Function Documentation

TCanvas * DrawResMomVsMom	(	const char *	name,
		const char *	title,
		TH2 *	h,
		Int_t	nBins,
		TF1 *	f2 = 0x0,
		const char *	fitting = ""
	)

generic function to draw and eventually fit momentum residuals versus momentum

Definition at line 1189 of file MUONRecoCheck.C.

Referenced by MUONRecoCheck().

TCanvas * DrawVsAng	(	const char *	name,
		const char *	title,
		TH1 *	h1,
		TH2 *	h2
	)

generic function to draw histograms versus absorber angular region

Definition at line 1155 of file MUONRecoCheck.C.

Referenced by MUONRecoCheck().

TCanvas * DrawVsPos	(	const char *	name,
		const char *	title,
		TH2 *	h1,
		TH2 *	h2,
		TH2 *	h3
	)

generic function to draw histograms versus position at absorber end

Definition at line 1174 of file MUONRecoCheck.C.

Referenced by MUONRecoCheck().

void FillResidual	(	TH1 *	h,
		Int_t	i,
		Double_t &	sigma,
		TGraphErrors *	gMean,
		TGraphErrors *	gSigma,
		Bool_t	correctForSystematics,
		Bool_t	fitResiduals
	)

fill graphs with residual mean and sigma

Definition at line 1116 of file MUONRecoCheck.C.

Referenced by MUONRecoCheck().

void FitGausResVsMom	(	TH2 *	h,
		Int_t	nBins,
		const Double_t	mean0,
		const Double_t	sigma0,
		const char *	fitting,
		TGraphAsymmErrors *	gMean,
		TGraphAsymmErrors *	gSigma
	)

generic function to fit residuals versus momentum with a gaussian

Definition at line 1056 of file MUONRecoCheck.C.

Referenced by MUONRecoCheck().

void FitPDCAVsMom	(	TH2 *	h,
		Int_t	nBins,
		const char *	fitting,
		TGraphAsymmErrors *	gMean,
		TGraphAsymmErrors *	gSigma
	)

generic function to fit p*DCA distributions

Definition at line 1087 of file MUONRecoCheck.C.

Referenced by MUONRecoCheck().

Double_t langaufun	(	Double_t *	x,
		Double_t *	par
	)

Definition at line 999 of file MUONRecoCheck.C.

Referenced by MUONRecoCheck().

void MUONRecoCheck	(	Int_t	nEvent = -1,
		const char *	pathSim = "./generated/",
		const char *	esdFileName = "AliESDs.root",
		const char *	ocdbPath = "local://$ALICE_ROOT/OCDB",
		const Double_t	pMin = 0.,
		const Double_t	pMax = 300.,
		const Int_t	pNBins = 30,
		Int_t	absorberRegion = -1,
		Bool_t	fitResiduals = kTRUE
	)

Associate the reconstructed tracks with the simulated ones and check the quality of the reconstruction (tracking/trigger efficiency; momentum, slope,... resolutions at first cluster and at vertex; cluster resolution).

• You can choose the momentum range and number of bins used to study the track resolution versus momentum.
• You can limit the calculation of track resolution at vertex to the tracks crossing the absorber in a given region with the flag "absorberRegion": -1=all, 1=[2,3]deg, 2=[3,10]deg.

Definition at line 73 of file MUONRecoCheck.C.

void Zoom	(	TH1 *	h,
		Double_t	fractionCut = 0.01
	)

Reduce the range of the histogram by removing a given fration of the statistic at each edge

Definition at line 1226 of file MUONRecoCheck.C.

Referenced by FillResidual().