New production macros (Yves)

[u/mrichter/AliRoot.git] / prod / LHC08d11 / rec.C

diff --git a/prod/LHC08d11/rec.C b/prod/LHC08d11/rec.C
new file mode 100644 (file)
index 0000000..a7b30f1
--- /dev/null
+++ b/prod/LHC08d11/rec.C
@@ -0,0 +1,69 @@
+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
+\r
+  TStopwatch timer;\r
+  timer.Start();\r
+  reco.Run();\r
+  timer.Stop();\r
+  timer.Print();\r
+}\r
+\r