3 # This script processes the following steps for n events (AliRoot v3-09-08,
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
22 # See ITSPIDplot.ps after run of this script.
24 if [ $# = 0 ]; then nev=1; else nev=$1; fi
25 if [ $nev = 0 ]; then nev=1; fi
27 # delete eventual old files from the last run
28 echo "Start simulation for " $nev " event(s)"
29 $ALICE_ROOT/ITS/AliITSDeleteOldFiles.sh
31 # run the hit generation
32 aliroot -q -b "$ALICE_ROOT/macros/grun.C($nev)"
34 aliroot -q -b "$ALICE_ROOT/TPC/AliTPCHits2Digits.C($nev)"
36 ln -s galice.root digits.root
37 aliroot -q -b "$ALICE_ROOT/TPC/AliTPCFindClustersMI.C($nev)"
39 .L $ALICE_ROOT/TPC/AliTPCFindTracksMI.C
40 AliTPCFindTracks($nev);
45 aliroot -q -b "$ALICE_ROOT/ITS/AliITSHits2SDigits.C"
46 aliroot -q -b "$ALICE_ROOT/ITS/AliITSSDigits2Digits.C"
48 aliroot -q -b "$ALICE_ROOT/ITS/AliITSFindClustersV2.C('s',$nev)"
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;}
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());
59 .x $ALICE_ROOT/ITS/AliITSFindTracksV2.C($nev);
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])"