Recent first physics productions.

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

void rec() {\r
\r
  AliReconstruction reco;\r
  \r
  reco.SetWriteESDfriend();\r
  reco.SetWriteAlignmentData();\r
  reco.SetRunReconstruction("ITS VZERO");\r
\r
  // ITS Efficiency\r
  reco.SetRunPlaneEff(kTRUE);\r
  AliITSRecoParam *itsRecoParam = AliITSRecoParam::GetPlaneEffParam(-1);\r
  itsRecoParam->SetOptTrackletsPlaneEff(kTRUE);\r
\r
\r
 //****** FIRST PHYSICS 2009 (same as COSMICS 2009) *********************\r
  //AliITSRecoParam * itsRecoParam = AliITSRecoParam::GetLowFluxParam();\r
   itsRecoParam->SetClusterErrorsParam(2);\r
   // find independently ITS SA tracks\r
   itsRecoParam->SetSAUseAllClusters();\r
   itsRecoParam->SetOuterStartLayerSA(AliITSgeomTGeo::GetNLayers()-2);\r
   // to maximize efficiency\r
   itsRecoParam->SetAllowProlongationWithEmptyRoad();     // larger seach windows for SA (in case of large misalignments)\r
   itsRecoParam->SetNLoopsSA(33);\r
   itsRecoParam->SetFactorSAWindowSizes(20);\r
   // additional error due to misal (B off)\r
   itsRecoParam->SetClusterMisalErrorY(1.0,1.0,1.0,1.0,1.0,1.0); // [cm]\r
   itsRecoParam->SetClusterMisalErrorZ(1.0,1.0,1.0,1.0,1.0,1.0); // [cm]\r
   // additional error due to misal (B on)\r
   itsRecoParam->SetClusterMisalErrorYBOn(0.0,0.0,0.1,0.1,0.1,0.1); // [cm]\r
   itsRecoParam->SetClusterMisalErrorZBOn(0.1,0.1,0.1,0.1,0.1,0.1); // [cm]\r
   //***********************************************************************\r
   itsRecoParam->SetEventSpecie(AliRecoParam::kLowMult);\r
\r
   reco.SetRecoParam("ITS",itsRecoParam);\r
\r
\r
\r
  reco.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");\r
//   reco.SetSpecificStorage("GRP/GRP/Data",\r
// 			  Form("local://%s",gSystem->pwd()));\r
  // We store the object in AliEn during the simulation\r
  reco.SetSpecificStorage("GRP/GRP/Data",\r
			  "alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");\r
\r
			  \r
  reco.SetSpecificStorage("GRP/Calib/MeanVertexSPD","alien://folder=/alice/cern.ch/user/r/rgrosso/ShiftedVertex");\r
\r
  TStopwatch timer;\r
  timer.Start();\r
  reco.Run();\r
  timer.Stop();\r
  timer.Print();\r
}\r