]>
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 |