[u/mrichter/AliRoot.git] / ITS / oldmacros / AliITStestPIDV2.sh

#!/bin/sh
# 12 May 2003,Dubna
# This script processes the following steps for n events (AliRoot v3-09-08,
# root v3.05/04)
# (use the parameter n):
# - the TPC and ITS slow simulation (hits+digits+clusters),
# - the TPC+ITS tracking,
# - the TPC and ITS PID
# - the AliITStracksPidV2.root is created with the following information:
#    fGlab - track number
#    fPcode - particle code after the TPC PID
#    fMom   - particle momentum      (from the AliTPCtracks.root)
#    fLam   - particle lambder angle (from the AliTPCtracks.root)
#    fPhi   - particle phi angle     (from the AliTPCtracks.root)
#    fSignal- TPC signal (mips)
#    fWpi   - the PID weight for the identified pion
#    fWk    - the PID weight for the identified kaon
#    fWp    - the PID weight for the identified proton
#       (the weghts are tuned for the HIJING events)
# and the AliITSScanPIDV2.C is the example of macros to read the PID
# information.
# See ITSPIDplot.ps after run of this script.

if [ $# = 0   ]; then nev=1; else nev=$1; fi
if [ $nev = 0 ]; then nev=1; fi
#
# delete eventual old files from the last run
echo "Start simulation for " $nev " event(s)"
$ALICE_ROOT/ITS/AliITSDeleteOldFiles.sh
#
# run the hit generation
aliroot -q -b "$ALICE_ROOT/macros/grun.C($nev)"  
# digitize TPC
aliroot -q -b "$ALICE_ROOT/TPC/AliTPCHits2Digits.C($nev)" 
# TPC tracking
ln -s galice.root digits.root
aliroot -q -b "$ALICE_ROOT/TPC/AliTPCFindClustersMI.C($nev)" 
aliroot -b <<EOI
.L $ALICE_ROOT/TPC/AliTPCFindTracksMI.C
AliTPCFindTracks($nev);
.q
EOI
#
# digitize ITS
aliroot -q -b "$ALICE_ROOT/ITS/AliITSHits2SDigits.C"
aliroot -q -b "$ALICE_ROOT/ITS/AliITSSDigits2Digits.C"
# ITS-TPC tracking
aliroot -q -b "$ALICE_ROOT/ITS/AliITSFindClustersV2.C('s',$nev)"
aliroot -b <<EOI
TFile *inkin = TFile::Open("galice.root");
if (!inkin->IsOpen()) cerr<<"Can't open galice.root !\n";               
if (gAlice) {delete gAlice; gAlice=0;}
                                  
gAlice = (AliRun*)inkin->Get("gAlice");
cout<<"AliRun object found on file "<<gAlice<<endl;
cout<<"!!!! field ="<<gAlice->Field()->SolenoidField()<<endl;
AliKalmanTrack::SetConvConst(1000/0.299792458/gAlice->Field()->SolenoidField());
inkin->Close();
.x $ALICE_ROOT/ITS/AliITSFindTracksV2.C($nev);
.q
EOI
#
# Do the PID procedure for the ITS. The output file AliITStrackV2Pid.root is
# created with the information for each track:
aliroot -q -b "$ALICE_ROOT/ITS/AliITSSavePIDV2.C(0,$[$nev-1])"
aliroot -q -b "$ALICE_ROOT/ITS/AliITSScanPIDV2.C(0,$[$nev-1])"
Commit	Line	Data
b3d9d240	1	#!/bin/sh
	2	# 12 May 2003,Dubna
	3	# This script processes the following steps for n events (AliRoot v3-09-08,
	4	# root v3.05/04)
	5	# (use the parameter n):
	6	# - the TPC and ITS slow simulation (hits+digits+clusters),
	7	# - the TPC+ITS tracking,
	8	# - the TPC and ITS PID
	9	# - the AliITStracksPidV2.root is created with the following information:
	10	# fGlab - track number
	11	# fPcode - particle code after the TPC PID
	12	# fMom - particle momentum (from the AliTPCtracks.root)
	13	# fLam - particle lambder angle (from the AliTPCtracks.root)
	14	# fPhi - particle phi angle (from the AliTPCtracks.root)
	15	# fSignal- TPC signal (mips)
	16	# fWpi - the PID weight for the identified pion
	17	# fWk - the PID weight for the identified kaon
	18	# fWp - the PID weight for the identified proton
	19	# (the weghts are tuned for the HIJING events)
	20	# and the AliITSScanPIDV2.C is the example of macros to read the PID
	21	# information.
	22	# See ITSPIDplot.ps after run of this script.
	23
	24	if [ $# = 0 ]; then nev=1; else nev=$1; fi
	25	if [ $nev = 0 ]; then nev=1; fi
	26	#
	27	# delete eventual old files from the last run
	28	echo "Start simulation for " $nev " event(s)"
	29	$ALICE_ROOT/ITS/AliITSDeleteOldFiles.sh
	30	#
	31	# run the hit generation
	32	aliroot -q -b "$ALICE_ROOT/macros/grun.C($nev)"
	33	# digitize TPC
	34	aliroot -q -b "$ALICE_ROOT/TPC/AliTPCHits2Digits.C($nev)"
	35	# TPC tracking
	36	ln -s galice.root digits.root
	37	aliroot -q -b "$ALICE_ROOT/TPC/AliTPCFindClustersMI.C($nev)"
	38	aliroot -b <<EOI
	39	.L $ALICE_ROOT/TPC/AliTPCFindTracksMI.C
	40	AliTPCFindTracks($nev);
	41	.q
	42	EOI
	43	#
	44	# digitize ITS
	45	aliroot -q -b "$ALICE_ROOT/ITS/AliITSHits2SDigits.C"
	46	aliroot -q -b "$ALICE_ROOT/ITS/AliITSSDigits2Digits.C"
	47	# ITS-TPC tracking
	48	aliroot -q -b "$ALICE_ROOT/ITS/AliITSFindClustersV2.C('s',$nev)"
	49	aliroot -b <<EOI
	50	TFile *inkin = TFile::Open("galice.root");
	51	if (!inkin->IsOpen()) cerr<<"Can't open galice.root !\n";
	52	if (gAlice) {delete gAlice; gAlice=0;}
	53
	54	gAlice = (AliRun*)inkin->Get("gAlice");
	55	cout<<"AliRun object found on file "<<gAlice<<endl;
	56	cout<<"!!!! field ="<<gAlice->Field()->SolenoidField()<<endl;
	57	AliKalmanTrack::SetConvConst(1000/0.299792458/gAlice->Field()->SolenoidField());
	58	inkin->Close();
	59	.x $ALICE_ROOT/ITS/AliITSFindTracksV2.C($nev);
	60	.q
	61	EOI
	62	#
	63	# Do the PID procedure for the ITS. The output file AliITStrackV2Pid.root is
	64	# created with the information for each track:
65	aliroot -q -b "$ALICE_ROOT/ITS/AliITSSavePIDV2.C(0,$[$nev-1])"
66	aliroot -q -b "$ALICE_ROOT/ITS/AliITSScanPIDV2.C(0,$[$nev-1])"
67
68
69