Working on the electron cut
[u/mrichter/AliRoot.git] / prod / LHC09d6 / rec.C
CommitLineData
b0e91ec7 1void rec() {\r
2\r
3 AliReconstruction reco;\r
4 \r
5 reco.SetWriteESDfriend();\r
6 reco.SetWriteAlignmentData();\r
7 reco.SetRunReconstruction("ITS VZERO");\r
8\r
9 // ITS Efficiency\r
10 reco.SetRunPlaneEff(kTRUE);\r
11 AliITSRecoParam *itsRecoParam = AliITSRecoParam::GetPlaneEffParam(-1);\r
12 itsRecoParam->SetOptTrackletsPlaneEff(kTRUE);\r
13\r
14\r
15 //****** FIRST PHYSICS 2009 (same as COSMICS 2009) *********************\r
16 //AliITSRecoParam * itsRecoParam = AliITSRecoParam::GetLowFluxParam();\r
17 itsRecoParam->SetClusterErrorsParam(2);\r
18 // find independently ITS SA tracks\r
19 itsRecoParam->SetSAUseAllClusters();\r
20 itsRecoParam->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);\r
21 // to maximize efficiency\r
22 itsRecoParam->SetAllowProlongationWithEmptyRoad(); // larger seach windows for SA (in case of large misalignments)\r
23 itsRecoParam->SetNLoopsSA(33);\r
24 itsRecoParam->SetFactorSAWindowSizes(20);\r
25 // additional error due to misal (B off)\r
26 itsRecoParam->SetClusterMisalErrorY(1.0,1.0,1.0,1.0,1.0,1.0); // [cm]\r
27 itsRecoParam->SetClusterMisalErrorZ(1.0,1.0,1.0,1.0,1.0,1.0); // [cm]\r
28 // additional error due to misal (B on)\r
29 itsRecoParam->SetClusterMisalErrorYBOn(0.0,0.0,0.1,0.1,0.1,0.1); // [cm]\r
30 itsRecoParam->SetClusterMisalErrorZBOn(0.1,0.1,0.1,0.1,0.1,0.1); // [cm]\r
31 //***********************************************************************\r
32 itsRecoParam->SetEventSpecie(AliRecoParam::kLowMult);\r
33\r
34 reco.SetRecoParam("ITS",itsRecoParam);\r
35\r
36\r
37\r
38 reco.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");\r
39// reco.SetSpecificStorage("GRP/GRP/Data",\r
40// Form("local://%s",gSystem->pwd()));\r
41 // We store the object in AliEn during the simulation\r
42 reco.SetSpecificStorage("GRP/GRP/Data",\r
43 "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");\r
44\r
45 \r
46 reco.SetSpecificStorage("GRP/Calib/MeanVertexSPD","alien://folder=/alice/cern.ch/user/r/rgrosso/ShiftedVertex");\r
47\r
48 TStopwatch timer;\r
49 timer.Start();\r
50 reco.Run();\r
51 timer.Stop();\r
52 timer.Print();\r
53}\r