AliRoot Core  3dc7879 (3dc7879)
makeBBFit.C File Reference

Macro to make the MC based dEdx fits, and study the influence of thre track selection to the PID. More...

#include "TCut.h"
#include "TFile.h"
#include "TTree.h"
#include "TH1.h"
#include "TH2.h"
#include "TH3.h"
#include "TF1.h"
#include "TCanvas.h"
#include "TDatabasePDG.h"
#include "TStatToolkit.h"
#include "TGraphErrors.h"
#include "TStopwatch.h"
#include "TLegend.h"
#include "TLatex.h"
#include "TTreeStream.h"
#include "AliTPCParam.h"
#include "AliTPCcalibBase.h"
Include dependency graph for makeBBFit.C:

Go to the source code of this file.

Functions

void Init ()
 
void FitTransferFunctionTheta (Bool_t isMax, Int_t dType, Int_t tgl, Bool_t skipKaon, Double_t maxP, TTreeSRedirector *pcstream, Bool_t doDraw)
 
void FitTransferFunctionAll (Bool_t isMax, Int_t dType, Bool_t skipKaon, Double_t maxP, TTreeSRedirector *pcstream, Bool_t doDraw)
 
void MakeESDCutsPID (Double_t fgeom=1, Double_t fcr=0.85, Double_t fcl=0.7)
 
void makeBBfit (Int_t ntracks=100000)
 
void FitTransferFunctionScanAll ()
 
void FitTransferFunctionScanTheta ()
 
void DrawFit ()
 

Variables

TCut cutGeom
 
TCut cutNcr
 
TCut cutNcl
 
TCut cutFiducial
 
Int_t pdgs [4] ={11,211,321,2212}
 
Double_t massPDG [4] ={0}
 
const char * partName [4] ={"Electron","Pion","Kaon","Proton"}
 
const Int_t colors [5] ={kGreen,kBlack,kRed,kBlue,5}
 
const Int_t markers [5] ={24,20,21,25,25}
 
const Float_t markerSize [5] ={1.2,1.2,1.0,1.2,1}
 
TTree * treeFit =0
 

Detailed Description

Macro to make the MC based dEdx fits, and study the influence of thre track selection to the PID.

Motivation

In the ALICE Geant3 MC the input Bethe-Bloch parameterization of the primary ionization can be parameterized by user defined formula.

In detector Input \( dE/dx(BG)_{in} \) (more exact \( dN_{prim}/dx \)) is transformed to the output reconstructed \( dE/dx_{rec} \). While original input function is just function of particle , random variable, reconstructed dEdx estimate, is influenced by detection processes and is sensitive to other aspects namily diffusion, track inclination angle ( \(\phi\), \(\theta\)), gain ...

In the following we will calibrate transform function:

\[ \frac{dE}{dx}_{BB}= f_{tr}\left(\frac{dE}{dx}_{rec}, \phi, \eta\right) \]

Example code for author:

.x $HOME/rootlogon.C
.x $HOME/NimStyle.C
.L $ALICE_ROOT/TPC/macros/data2011/makeBBFit.C+
Init();
makeBBfit(10000000);

Definition in file makeBBFit.C.

Function Documentation

void DrawFit ( )

Definition at line 609 of file makeBBFit.C.

void FitTransferFunctionAll ( Bool_t  isMax,
Int_t  dType,
Bool_t  skipKaon,
Double_t  maxP,
TTreeSRedirector pcstream,
Bool_t  doDraw 
)

dEdx_{BB}= f_{tr}(dEdx_{rec}, #phi,#eta)

Example usage: FitTransferFunctionAll(1,1,1,15,0,1) isMax=1; dType=1; skipKaon=0; Double_t maxP=10

Definition at line 440 of file makeBBFit.C.

Referenced by FitTransferFunctionScanAll().

void FitTransferFunctionScanAll ( )

Make fit of transfer functions - parabolic in theta

Definition at line 250 of file makeBBFit.C.

void FitTransferFunctionScanTheta ( )

Make fit of transfer functions - bin by bin in Theta

Definition at line 272 of file makeBBFit.C.

void FitTransferFunctionTheta ( Bool_t  isMax,
Int_t  dType,
Int_t  tgl,
Bool_t  skipKaon,
Double_t  maxP,
TTreeSRedirector pcstream,
Bool_t  doDraw 
)

dEdx_{BB}= f_{tr}(dEdx_{rec}, #phi,#eta)

Fit the parematers of transfer function 4 models of transfer function used:

Example usage: FitTransferFunctionTheta(1,3,8,1,5,0) isMax=1; dType=1; tgl=5; skipKaon=0; Double_t maxP=10

Definition at line 297 of file makeBBFit.C.

Referenced by FitTransferFunctionScanTheta().

void Init ( )

1.) Register default ALICE parametrization

Definition at line 68 of file makeBBFit.C.

void makeBBfit ( Int_t  ntracks = 100000)

Make dEdx fits of indiviadual particle species in bins: a.) momenta b.) tan(theta) -tgl c.) per detector segmen d.) Qmax/Qtot

Definition at line 119 of file makeBBFit.C.

void MakeESDCutsPID ( Double_t  fgeom = 1,
Double_t  fcr = 0.85,
Double_t  fcl = 0.7 
)

Cuts to be used in the sequenece 1.) cutGeom - stronger cut on geometry - perfect agreement data-MC

  • Default length factor 1 2.) cutNcr - relativally stong cut on the number of crossed rows to gaurantee tracking performance relativally good agrement in the MC except of the week decay propability
  • Default length factor 0.85 3.) cutNcl - very week cut on the numebr of clusters
  • Default length factor 0.7

Combined cuts should be combination: cutGeom+cutNcr+cutNcl Default factors 1, 0.85 and 0.7 should be suited for most of analysis Modification can be expected for the Jet analysis and further tuning of the parameters for the relativistic PID analysis

Definition at line 77 of file makeBBFit.C.

Variable Documentation

const Int_t colors[5] ={kGreen,kBlack,kRed,kBlue,5}

Definition at line 57 of file makeBBFit.C.

Referenced by FitTransferFunctionAll(), and FitTransferFunctionTheta().

TCut cutFiducial

Definition at line 53 of file makeBBFit.C.

Referenced by makeBBfit(), and MakeESDCutsPID().

TCut cutGeom

Definition at line 52 of file makeBBFit.C.

Referenced by makeBBfit(), and MakeESDCutsPID().

TCut cutNcl

Definition at line 52 of file makeBBFit.C.

Referenced by makeBBfit(), and MakeESDCutsPID().

TCut cutNcr

Definition at line 52 of file makeBBFit.C.

Referenced by makeBBfit(), and MakeESDCutsPID().

const Int_t markers[5] ={24,20,21,25,25}
const Float_t markerSize[5] ={1.2,1.2,1.0,1.2,1}
Double_t massPDG[4] ={0}

Definition at line 55 of file makeBBFit.C.

Referenced by Init(), and makeBBfit().

const char* partName[4] ={"Electron","Pion","Kaon","Proton"}

Definition at line 56 of file makeBBFit.C.

Referenced by FitTransferFunctionAll(), FitTransferFunctionTheta(), and makeBBfit().

Int_t pdgs[4] ={11,211,321,2212}

Definition at line 54 of file makeBBFit.C.

Referenced by Init(), and makeBBfit().