21d4a8d1 |
1 | void rec(const char *filename="raw.root")\r |
40d0f38f |
2 | {\r |
40d0f38f |
3 | /////////////////////////////////////////////////////////////////////////////////////////\r |
4 | //\r |
5 | // First version of the reconstruction\r |
6 | // script for the FDR'08\r |
7 | \r |
8 | // Set the CDB storage location\r |
9 | // AliLog::SetModuleDebugLevel("STEER",2);\r |
10 | AliCDBManager * man = AliCDBManager::Instance();\r |
11 | // man->SetDefaultStorage("local://LocalCDB");\r |
12 | man->SetDefaultStorage("alien://folder=/alice/data/2008/LHC08a/OCDB/");\r |
13 | \r |
462112b2 |
14 | // Files that we can not read from alien...solved\r |
15 | // man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local://$ALICE_ROOT");\r |
16 | // man->SetSpecificStorage("ITS/Calib/MapsTimeSDD","local://$ALICE_ROOT");\r |
17 | // man->SetSpecificStorage("TPC/Calib/ExB","local://$ALICE_ROOT");\r |
40d0f38f |
18 | \r |
462112b2 |
19 | // Objects not found if using LHC07w database...solved\r |
40d0f38f |
20 | // man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local:///afs/cern.ch/user/c/cheshkov/public/OCDB");\r |
21 | // man->SetSpecificStorage("GRP/GRP/Data","local://$ALICE_ROOT");\r |
22 | // man->SetSpecificStorage("ITS/Calib/DDLMapSDD","local://$ALICE_ROOT");\r |
23 | // man->SetSpecificStorage("MUON/Calib/Mapping","local://$ALICE_ROOT");\r |
24 | // man->SetSpecificStorage("MUON/Calib/DDLStore","local://$ALICE_ROOT");\r |
25 | \r |
26 | // ITS settings\r |
27 | AliITSRecoParam * itsRecoParam = AliITSRecoParam::GetCosmicTestParam();\r |
28 | itsRecoParam->SetClusterErrorsParam(2);\r |
29 | itsRecoParam->SetFindV0s(kFALSE);\r |
30 | itsRecoParam->SetAddVirtualClustersInDeadZone(kFALSE);\r |
31 | itsRecoParam->SetUseAmplitudeInfo(kFALSE);\r |
32 | // In case we want to switch off a layer\r |
33 | // itsRecoParam->SetLayerToSkip(<N>);\r |
34 | itsRecoParam->SetLayerToSkip(4);\r |
35 | itsRecoParam->SetLayerToSkip(5);\r |
36 | itsRecoParam->SetLayerToSkip(2);\r |
37 | itsRecoParam->SetLayerToSkip(3);\r |
38 | AliITSReconstructor::SetRecoParam(itsRecoParam);\r |
39 | \r |
40 | // TPC settings\r |
41 | AliLog::SetClassDebugLevel("AliTPCclustererMI",2);\r |
42 | AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kTRUE);\r |
43 | tpcRecoParam->SetTimeInterval(60,940);\r |
44 | tpcRecoParam->Dump();\r |
45 | AliTPCReconstructor::SetRecoParam(tpcRecoParam);\r |
46 | AliTPCReconstructor::SetStreamLevel(1);\r |
47 | \r |
462112b2 |
48 | // TRD setting\r |
49 | AliTRDrawStreamBase::SetRawStreamVersion("TB");\r |
50 | \r |
40d0f38f |
51 | // PHOS settings\r |
52 | AliPHOSRecoParam* recEmc = new AliPHOSRecoParamEmc();\r |
53 | recEmc->SetSubtractPedestals(kTRUE);\r |
54 | recEmc->SetMinE(0.05);\r |
55 | recEmc->SetClusteringThreshold(0.10);\r |
56 | AliPHOSReconstructor::SetRecoParamEmc(recEmc);\r |
57 | \r |
58 | // T0 settings\r |
59 | AliLog::SetModuleDebugLevel("T0", 10);\r |
60 | \r |
61 | // MUON settings\r |
62 | AliLog::SetClassDebugLevel("AliMUONRawStreamTracker",3);\r |
63 | AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();\r |
64 | muonRecoParam->CombineClusterTrackReco(kTRUE);\r |
65 | muonRecoParam->SetCalibrationMode("NOGAIN");\r |
66 | //muonRecoParam->SetClusteringMode("PEAKFIT");\r |
67 | //muonRecoParam->SetClusteringMode("PEAKCOG");\r |
68 | muonRecoParam->Print("FULL");\r |
69 | AliRecoParam::Instance()->RegisterRecoParam(muonRecoParam);\r |
70 | \r |
71 | // Tracking settings\r |
72 | // AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);\r |
73 | AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 0., 10., 2);\r |
74 | AliTracker::SetFieldMap(field,1);\r |
75 | Double_t mostProbPt=0.35;\r |
76 | AliExternalTrackParam::SetMostProbablePt(mostProbPt);\r |
77 | \r |
78 | // AliReconstruction settings\r |
79 | AliReconstruction rec;\r |
80 | rec.SetUniformFieldTracking(kFALSE);\r |
81 | rec.SetWriteESDfriend(kTRUE);\r |
82 | rec.SetWriteAlignmentData();\r |
21d4a8d1 |
83 | rec.SetInput(filename);\r |
40d0f38f |
84 | rec.SetRunReconstruction("ALL");\r |
85 | rec.SetUseTrackingErrorsForAlignment("ITS");\r |
86 | \r |
87 | // In case some detectors have to be switched off...\r |
88 | // rec.SetRunLocalReconstruction("ALL");\r |
89 | // rec.SetRunTracking("ALL");\r |
90 | // rec.SetFillESD("ALL");\r |
62b2441a |
91 | // Enable vertex finder - it is needed for cosmic track reco\r |
92 | rec.SetRunVertexFinder(kTRUE);\r |
40d0f38f |
93 | \r |
94 | // To be enabled if some equipment IDs are not set correctly by DAQ\r |
95 | // rec.SetEquipmentIdMap("EquipmentIdMap.data");\r |
96 | \r |
97 | // Detector options if any\r |
98 | rec.SetOption("ITS","cosmics,onlyITS");\r |
99 | rec.SetOption("MUON","SAVEDIGITS");\r |
100 | rec.SetOption("TPC","OldRCUFormat");\r |
101 | rec.SetOption("PHOS","OldRCUFormat");\r |
cbc0d7f8 |
102 | rec.SetOption("T0","cosmic");\r |
40d0f38f |
103 | \r |
104 | // To be enabled when CTP readout starts\r |
105 | rec.SetFillTriggerESD(kFALSE);\r |
106 | \r |
107 | // all events in one single file\r |
108 | rec.SetNumberOfEventsPerFile(-1);\r |
109 | \r |
110 | // switch off cleanESD\r |
111 | rec.SetCleanESD(kFALSE);\r |
112 | \r |
113 | // rec.SetEventRange(0,15);\r |
114 | // AliLog::SetGlobalDebugLevel(2);\r |
115 | \r |
116 | rec.SetRunQA(kFALSE);\r |
117 | AliLog::Flush();\r |
118 | rec.Run();\r |
119 | \r |
120 | //cout << "-----------------------------------------------------------------" << endl;\r |
121 | //cout << "-----------------------------------------------------------------" << endl;\r |
122 | //cout << "--------- Reconstruction Completed. Start merging QAs -----------" << endl;\r |
123 | //cout << "-----------------------------------------------------------------" << endl;\r |
124 | //cout << "-----------------------------------------------------------------" << endl;\r |
125 | //AliQADataMakerSteer qas;\r |
126 | //qas.Merge();\r |
127 | }\r |