]>
Commit | Line | Data |
---|---|---|
5225c9db | 1 | void 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 |