PWGGA/PHOSTasks/PHOS_embedding/sim.C

   1 void sim(Int_t nev=50) {
   2
   3   if (gSystem->Getenv("SIM_EVENTS"))
   4     nev = atoi(gSystem->Getenv("SIM_EVENTS"));
   5
   6   printf("GENERATE << %d >> events \n",nev);
   7
   8   gROOT->LoadMacro("IpPion.C++") ;
   9
  10
  11   AliSimulation simulator;
  12   simulator.SetMakeSDigits("PHOS");
  13   simulator.SetMakeDigits("PHOS");
  14 //
  15 // Ideal OCDB
  16 //  simulator.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  17   simulator.SetDefaultStorage("local://./OCDB");
  18 //  simulator.SetSpecificStorage("GRP/GRP/Data",
  19 //                               Form("local://%s",gSystem->pwd()));
  20
  21 //  simulator.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
  22
  23  //simulator.SetSpecificStorage("GRP/Calib/MeanVertexSPD", "alien://folder=/alice/data/2010/OCDB");
  24
  25   //PHOS bad map from RAW OCDB
  26   simulator.SetSpecificStorage("PHOS/*/*/","local://./OCDB");
  27 //  simulator.SetSpecificStorage("PHOS/Calib/EmcBadChannels/","local://./OCDB");
  28 //  simulator.SetSpecificStorage("PHOS/Calib/EmcGainPedestals/","local://./OCDB");
  29
  30   simulator.SetRunHLT("");
  31 //
  32
  33   simulator.SetSpecificStorage("GRP/GRP/Data", "alien://Folder=/alice/data/2010/OCDB");
  34
  35   AliPHOSSimParam *simParam =  AliPHOSSimParam::GetInstance() ;
  36   simParam->SetAPDNoise(0.000001) ;
  37   simParam->SetCellNonLineairyA(0.001) ;
  38 //  simParam->SetCellNonLineairyA(0.1) ; //Default
  39   simParam->SetCellNonLineairyB(0.2) ;
  40 //  simParam->SetCellNonLineairyC(0.989) ; //Jan4
  41 //  simParam->SetCellNonLineairyC(0.995) ; //Jan5 - 2GeV
  42   simParam->SetCellNonLineairyC(1.031) ; //no NL
  43
  44 // Vertex and Mag.field from OCDB
  45
  46 //  simulator.UseVertexFromCDB();
  47   simulator.UseMagFieldFromGRP();
  48   simulator.SetRunQA(":") ;
  49
  50 //
  51 // The rest
  52
  53   TStopwatch timer;
  54   timer.Start();
  55   simulator.Run(nev);
  56   timer.Stop();
  57   timer.Print();
  58 }