Updated TRD reconstruction, modifications for LHC08d period
[u/mrichter/AliRoot.git] / test / cosmic / rec.C
CommitLineData
5225c9db 1void rec(const char *filename="raw.root", const Int_t mfield=1)\r
40d0f38f 2{\r
40d0f38f 3 /////////////////////////////////////////////////////////////////////////////////////////\r
4 //\r
09244196 5 // Second version of the reconstruction\r
6 // script for the 2008 cosmic data (LHC08b) \r
7 //\r
8 /////////////////////////////////////////////////////////////////////////////////////////\r
9 //AliLog::SetGlobalLogLevel(AliLog::kWarning);\r
5225c9db 10 AliLog::SetGlobalLogLevel(AliLog::kError);\r
09244196 11\r
12 gSystem->Load("libRAliEn.so");\r
13 gSystem->Load("libNet.so");\r
5225c9db 14 gSystem->Load("libMonaLisa.so");\r
15 new TMonaLisaWriter(0, "GridAliRoot-rec.C", 0, 0, "global");\r
16 gSystem->Setenv("APMON_INTERVAL", "120");\r
40d0f38f 17\r
18 // Set the CDB storage location\r
40d0f38f 19 AliCDBManager * man = AliCDBManager::Instance();\r
b8df7b30 20// man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");\r
21 man->SetDefaultStorage("alien://folder=/alice/data/2008/LHC08d/OCDB/");\r
22 man->SetSpecificStorage("ITS/Calib/*","local://$ALICE_ROOT/OCDB");\r
40d0f38f 23 \r
09244196 24 // Example in case a specific CDB storage is needed\r
162637e4 25 // man->SetSpecificStorage("ITS/Calib/MapsAnodeSDD","local://$ALICE_ROOT/OCDB");\r
40d0f38f 26\r
b8df7b30 27 AliQA::SetQARefStorage("local://$ALICE_ROOT/QAref") ;\r
28\r
09244196 29 // Reconstruction settings\r
30 AliReconstruction rec;\r
40d0f38f 31\r
32 // ITS settings\r
33 AliITSRecoParam * itsRecoParam = AliITSRecoParam::GetCosmicTestParam();\r
09244196 34 itsRecoParam->SetFactorSAWindowSizes(20);\r
40d0f38f 35 itsRecoParam->SetClusterErrorsParam(2);\r
36 itsRecoParam->SetFindV0s(kFALSE);\r
37 itsRecoParam->SetAddVirtualClustersInDeadZone(kFALSE);\r
38 itsRecoParam->SetUseAmplitudeInfo(kFALSE);\r
39 // In case we want to switch off a layer\r
40 // itsRecoParam->SetLayerToSkip(<N>);\r
09244196 41 // itsRecoParam->SetLayerToSkip(4);\r
42 // itsRecoParam->SetLayerToSkip(5);\r
40d0f38f 43 itsRecoParam->SetLayerToSkip(2);\r
44 itsRecoParam->SetLayerToSkip(3);\r
09244196 45 //itsRecoParam->SetSAOnePointTracks();\r
46 itsRecoParam->SetClusterMisalError(0.1); // [cm]\r
47 itsRecoParam->SetSAUseAllClusters();\r
48 rec.SetRecoParam("ITS",itsRecoParam);\r
40d0f38f 49\r
50 // TPC settings\r
09244196 51 //AliLog::SetClassDebugLevel("AliTPCclustererMI",2);\r
52 AliTPCRecoParam * tpcRecoParam = AliTPCRecoParam::GetCosmicTestParam(kFALSE);\r
40d0f38f 53 tpcRecoParam->SetTimeInterval(60,940);\r
09244196 54 Double_t sysError[5]={0.3,1, 0.3/150., 1./150.,0.3/(150*150.)};\r
55 tpcRecoParam->SetSystematicError(sysError);\r
56 tpcRecoParam->SetMinMaxCutAbs(4.);\r
57 tpcRecoParam->SetMinLeftRightCutAbs(6.);\r
58 tpcRecoParam->SetMinUpDownCutAbs(6.);\r
59 // tpcRecoParam->Dump();\r
60 rec.SetRecoParam("TPC",tpcRecoParam);\r
40d0f38f 61 AliTPCReconstructor::SetStreamLevel(1);\r
62\r
462112b2 63 // TRD setting\r
09244196 64 // Settings for the TRD Raw Reader\r
462112b2 65 AliTRDrawStreamBase::SetRawStreamVersion("TB");\r
09244196 66 \r
67 // TRD reconstruction params\r
68 AliTRDrecoParam *fTRDrecoParam = AliTRDrecoParam::GetCosmicTestParam();\r
69 rec.SetRecoParam("TRD",fTRDrecoParam);\r
70 AliTRDtrackerV1::SetNTimeBins(30);\r
462112b2 71\r
40d0f38f 72 // PHOS settings\r
09244196 73 AliPHOSRecoParam* recPHOS = new AliPHOSRecoParam();\r
74 recPHOS->SetEMCSubtractPedestals(kTRUE);\r
75 recPHOS->SetEMCMinE(0.05);\r
76 recPHOS->SetEMCClusteringThreshold(0.10);\r
77 rec.SetRecoParam("PHOS",recPHOS);\r
40d0f38f 78\r
79 // T0 settings\r
09244196 80 //AliLog::SetModuleDebugLevel("T0", 10);\r
40d0f38f 81\r
82 // MUON settings\r
09244196 83 //AliLog::SetClassDebugLevel("AliMUONRawStreamTracker",3);\r
40d0f38f 84 AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();\r
85 muonRecoParam->CombineClusterTrackReco(kTRUE);\r
86 muonRecoParam->SetCalibrationMode("NOGAIN");\r
87 //muonRecoParam->SetClusteringMode("PEAKFIT");\r
88 //muonRecoParam->SetClusteringMode("PEAKCOG");\r
89 muonRecoParam->Print("FULL");\r
09244196 90 rec.SetRecoParam("MUON",muonRecoParam);\r
40d0f38f 91 \r
92 // Tracking settings\r
f7a1cc68 93 AliMagF* field;\r
5225c9db 94 if (mfield)\r
f7a1cc68 95 field = new AliMagF("Maps","Maps", 2, 1., 1., 10., AliMagF::k5kG);\r
5225c9db 96 else\r
f7a1cc68 97 field = new AliMagF("Maps","Maps", 2, 0., 0., 10., AliMagF::k2kG);\r
98\r
40d0f38f 99 Double_t mostProbPt=0.35;\r
100 AliExternalTrackParam::SetMostProbablePt(mostProbPt);\r
101\r
102 // AliReconstruction settings\r
40d0f38f 103 rec.SetUniformFieldTracking(kFALSE);\r
104 rec.SetWriteESDfriend(kTRUE);\r
105 rec.SetWriteAlignmentData();\r
21d4a8d1 106 rec.SetInput(filename);\r
09244196 107 // rec.SetRunReconstruction("ALL");\r
40d0f38f 108 rec.SetUseTrackingErrorsForAlignment("ITS");\r
109\r
110 // In case some detectors have to be switched off...\r
5225c9db 111 rec.SetRunReconstruction("ITS TPC TRD TOF HMPID PHOS MUON FMD PMD T0 VZERO ZDC ACORDE");\r
09244196 112\r
62b2441a 113 // Enable vertex finder - it is needed for cosmic track reco\r
114 rec.SetRunVertexFinder(kTRUE);\r
40d0f38f 115\r
116 // To be enabled if some equipment IDs are not set correctly by DAQ\r
117 // rec.SetEquipmentIdMap("EquipmentIdMap.data");\r
118\r
119 // Detector options if any\r
09244196 120 rec.SetOption("ITS","cosmics");\r
40d0f38f 121 rec.SetOption("MUON","SAVEDIGITS");\r
cbc0d7f8 122 rec.SetOption("T0","cosmic");\r
40d0f38f 123\r
09244196 124 // Enabled when CTP readout starts\r
125 rec.SetFillTriggerESD(kTRUE);\r
40d0f38f 126\r
127 // all events in one single file\r
128 rec.SetNumberOfEventsPerFile(-1);\r
129\r
130 // switch off cleanESD\r
131 rec.SetCleanESD(kFALSE);\r
132\r
09244196 133 //rec.SetEventRange(0,15);\r
40d0f38f 134\r
b8df7b30 135 rec.SetRunQA(":");\r
09244196 136 rec.SetRunGlobalQA(kTRUE);\r
40d0f38f 137 AliLog::Flush();\r
138 rec.Run();\r
139\r
140 //cout << "-----------------------------------------------------------------" << endl;\r
141 //cout << "-----------------------------------------------------------------" << endl;\r
142 //cout << "--------- Reconstruction Completed. Start merging QAs -----------" << endl;\r
143 //cout << "-----------------------------------------------------------------" << endl;\r
144 //cout << "-----------------------------------------------------------------" << endl;\r
145 //AliQADataMakerSteer qas;\r
146 //qas.Merge();\r
147}\r