Recent first physics productions.
[u/mrichter/AliRoot.git] / prod / LHC09a14 / rec.C
index a7b30f16d9f177bf3377685b779cc75c8ed9f616..a632e08918f1edf36c63ff416480e02b085b4897 100644 (file)
@@ -1,69 +1,48 @@
 void rec() {\r
-  //AliLog::SetGlobalLogLevel(AliLog::kError);\r
 \r
-  AliReconstruction reco;\r
-\r
-  reco.SetWriteESDfriend();\r
-  reco.SetWriteAlignmentData();\r
-\r
-\r
-  AliITSRecoParam * itspar = AliITSRecoParam::GetLowFluxParam();\r
-  itspar->SetStoreLikeSignV0s(kTRUE);\r
-  reco.SetRecoParam("ITS",itspar);\r
-  reco.SetRecoParam("TPC",AliTPCRecoParam::GetLowFluxParam());\r
-  reco.SetRecoParam("TRD",AliTRDrecoParam::GetLowFluxParam());\r
-  reco.SetRecoParam("PHOS",AliPHOSRecoParam::GetDefaultParameters());\r
-  reco.SetRecoParam("MUON",AliMUONRecoParam::GetLowFluxParam());\r
-  reco.SetRecoParam("EMCAL",AliEMCALRecParam::GetLowFluxParam());\r
-  reco.SetRecoParam("GRP",AliGRPRecoParam::GetLowFluxParam());\r
-\r
-  reco.SetOption("TRD","sl_tr_1");            // Stream Level for the tracker equal to 1\r
-\r
-  // Only in case of Full misalignment\r
-//   AliGRPRecoParam *grpRecoParam = AliGRPRecoParam::GetLowFluxParam();\r
-//   grpRecoParam->SetVertexerTracksConstraintITS(kFALSE);\r
-//   grpRecoParam->SetVertexerTracksConstraintTPC(kFALSE);\r
-//   reco.SetRecoParam("GRP",grpRecoParam);\r
-\r
-  reco.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");\r
-  //reco.SetSpecificStorage("GRP/GRP/Data/","alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
-  reco.SetSpecificStorage("GRP/GRP/Data",Form("local://%s",gSystem->pwd()));
-  // No write access to the OCDB => local specific storage\r
-//   reco.SetSpecificStorage("GRP/GRP/Data",\r
-//                        Form("local://%s",gSystem->pwd()));\r
-\r
-  //-------------------------------------------------------------------------\r
-  // Setting the cuts for the V0 and cascade finding\r
-  // The values of the cuts below are "reasonable" for pp events\r
-  //-------------------------------------------------------------------------\r
-\r
-  Double_t v0sels[]={33,    // max allowed chi2\r
-                     0.05,  // min allowed impact parameter for the 1st daughter\r
-                     0.05,  // min allowed impact parameter for the 2nd daughter\r
-                     0.5,   // max allowed DCA between the daughter tracks\r
-                     0.99,  // max allowed cosine of V0's pointing angle\r
-                     0.2,   // min radius of the fiducial volume\r
-                     100    // max radius of the fiducial volume\r
-  };\r
-  AliV0vertexer::SetDefaultCuts(v0sels);\r
-\r
-  Double_t xisels[]={33.,   // max allowed chi2 (same as PDC07)\r
-                     0.025, // min allowed V0 impact parameter (PDC07 was 0.05)\r
-                     0.010, // "window" around the Lambda mass (PDC07 was 0.008)\r
-                     0.025, // min allowed bachelor's impact parameter (PDC07 was 0.035)\r
-                     0.2,   // max allowed DCA between the V0 and the bachelor (PDC07 was 0.1)\r
-                     0.998, // max allowed cosine of the cascade pointing angle (PDC07 was 0.9985)\r
-                     0.2,   // min radius of the fiducial volume (PDC07 was 0.9)\r
-                     100    // max radius of the fiducial volume (same as PDC07)\r
-  };\r
-  AliCascadeVertexer::SetDefaultCuts(xisels);\r
-\r
-  reco.SetRunQA("ALL:ALL");\r
+  AliReconstruction MuonRec;\r
+\r
+  MuonRec.SetInput("raw.root");\r
+  MuonRec.SetRunLocalReconstruction("MUON ITS VZERO");\r
+  MuonRec.SetRunTracking("MUON ITS VZERO");\r
+  MuonRec.SetRunVertexFinder(kTRUE);\r
+  MuonRec.SetFillESD("MUON VZERO HLT");\r
+  MuonRec.SetRunQA("MUON:ALL");\r
+  \r
+  MuonRec.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");\r
+\r
+  // QA reference\r
+  MuonRec.SetQARefDefaultStorage("local://$ALICE_ROOT/QAref") ;\r
+\r
+  // HLT\r
+  MuonRec.SetSpecificStorage("HLT/ConfigMUON/DecisionComponent","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");\r
+  MuonRec.SetSpecificStorage("HLT/ConfigMUON/HitReconstructor","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");\r
+  MuonRec.SetSpecificStorage("HLT/ConfigMUON/MansoTrackerFSM","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");\r
+  MuonRec.SetSpecificStorage("HLT/ConfigMUON/TriggerReconstructor","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");\r
+  \r
+  // CTP\r
+  MuonRec.SetSpecificStorage("GRP/CTP/Config","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb");\r
+  \r
+  // tracker masks\r
+  MuonRec.SetSpecificStorage("MUON/Calib/Gains","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb/rec");\r
+  MuonRec.SetSpecificStorage("MUON/Calib/Pedestals","alien://Folder=/alice/cern.ch/user/b/bogdan/prod2009/cdb/rec");\r
+\r
+  // reconstruction parameters\r
+  AliMUONRecoParam *muonRecoParam = AliMUONRecoParam::GetLowFluxParam();\r
+  muonRecoParam->SaveFullClusterInESD(kTRUE,100.);\r
+  for (Int_t i = 0; i < 10; i++) {\r
+    muonRecoParam->SetDefaultNonBendingReso(i, 0.1524);\r
+    muonRecoParam->SetDefaultBendingReso(i, 0.05099);\r
+  }\r
+  //muonRecoParam->Print("FULL");\r
+  MuonRec.SetRecoParam("MUON",muonRecoParam);\r
+\r
+  MuonRec.SetNumberOfEventsPerFile(500);\r
 \r
   TStopwatch timer;\r
   timer.Start();\r
-  reco.Run();\r
+  MuonRec.Run();\r
   timer.Stop();\r
   timer.Print();\r
+  \r
 }\r
-\r