update of package

[u/mrichter/AliRoot.git] / test / embedding / simrun.C

// #define VERBOSEARGS\r
// simrun.C\r
{\r
// set job and simulation variables as :\r
// root.exe -b -q simrun.C  --run <x> --event <y> --process <kPythia6/kPhojet> --field <kNoField/k5kG> --energy <900/10000>\r
\r
  int nrun = 0;\r
  int nevent = 0;\r
  int seed = 0;\r
\r
  char sseed[1024];\r
  char srun[1024];\r
  char sevent[1024];\r
  char sprocess[1024];\r
  char sfield[1024];\r
  char senergy[1024];\r
\r
  sprintf(srun,"");\r
  sprintf(sevent,"");\r
  sprintf(sprocess,"");\r
  sprintf(sfield,"");\r
  sprintf(senergy,"");\r
\r
  for (int i=0; i< gApplication->Argc();i++){\r
#ifdef VERBOSEARGS\r
    printf("Arg  %d:  %s\n",i,gApplication->Argv(i));\r
#endif\r
    if (!(strcmp(gApplication->Argv(i),"--run")))\r
      nrun = atoi(gApplication->Argv(i+1));\r
    sprintf(srun,"%d",nrun);\r
\r
    if (!(strcmp(gApplication->Argv(i),"--event")))\r
      nevent = atoi(gApplication->Argv(i+1));\r
    sprintf(sevent,"%d",nevent);\r
\r
    if (!(strcmp(gApplication->Argv(i),"--process")))\r
      sprintf(sprocess, gApplication->Argv(i+1));\r
\r
    if (!(strcmp(gApplication->Argv(i),"--field")))\r
      sprintf(sfield,gApplication->Argv(i+1));\r
\r
    if (!(strcmp(gApplication->Argv(i),"--energy")))\r
      sprintf(senergy,gApplication->Argv(i+1));\r
\r
  }\r
\r
  seed = nrun * 100000 + nevent;\r
  sprintf(sseed,"%d",seed);\r
\r
  if (seed==0) {\r
    fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
    fprintf(stderr,"!!!!  WARNING! Seeding variable for MC is 0          !!!!\n");\r
    fprintf(stderr,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
  } else {\r
    fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
    fprintf(stdout,"!!!  MC Seed is %d \n",seed);\r
    fprintf(stdout,"!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!\n");\r
  }\r
  \r
// set the seed environment variable\r
  gSystem->Setenv("CONFIG_SEED",sseed);\r
  gSystem->Setenv("CONFIG_RUN_TYPE",sprocess); // kPythia6 or kPhojet\r
  gSystem->Setenv("CONFIG_FIELD",sfield);      // kNoField or k5kG\r
  gSystem->Setenv("CONFIG_ENERGY",senergy);    // 900 or 10000 (GeV)\r
  gSystem->Setenv("DC_RUN",srun); // Not used in Config.C\r
  gSystem->Setenv("DC_EVENT",sevent); // Not used in Config.C\r
  \r
// Needed to produce simulated RAW data\r
  gSystem->Setenv("ALIMDC_RAWDB1","./mdc1");\r
  gSystem->Setenv("ALIMDC_RAWDB2","./mdc2");\r
  gSystem->Setenv("ALIMDC_TAGDB","./mdc1/tag");\r
  gSystem->Setenv("ALIMDC_RUNDB","./mdc1/meta");\r
  cout    << "SIMRUN:: Run " << gSystem->Getenv("DC_RUN") << " Event " << gSystem->Getenv("DC_EVENT")\r
	  << endl;\r
\r
\r
  // Background simulation\r
  gSystem->Setenv("CONFIG_EMBEDDING","kBackground");\r
\r
  cout<<">>>>> BACKGROUND SIMULATION <<<<<"<<endl;\r
  gSystem->Exec("mkdir BackgroundFull");\r
  gSystem->Exec("cp Config.C BackgroundFull/");\r
  gSystem->Exec("cp sim.C BackgroundFull/");\r
  gSystem->Exec("cp rec.C BackgroundFull/");\r
  gSystem->ChangeDirectory("BackgroundFull/");\r
  gSystem->Exec("aliroot -b -q 'sim.C(0)' > sim.log 2>&1");\r
  cout<<">>>>> BACKGROUND RECONSTRUCTION <<<<<"<<endl;\r
  gSystem->Exec("aliroot -b -q 'rec.C(0)' > rec.log 2>&1");\r
  gSystem->ChangeDirectory("../");\r
\r
  // Convert Raw to SDigits\r
  cout << ">>>>> CONVERTING RAW 2 SDIGITS <<<<<" << endl;\r
  gSystem->Exec("mkdir Background");\r
  gSystem->Exec("cp BackgroundFull/raw.root Background/");\r
  gSystem->Exec("cp BackgroundFull/AliESDs.root Background/");\r
  gSystem->Exec("cp -a BackgroundFull/GRP Background/");\r
  gSystem->Exec("cp sim.C Background/");\r
  gSystem->Exec("cp Config.C Background/");\r
  gSystem->ChangeDirectory("Background/");\r
  gSystem->Exec("aliroot -b -q 'sim.C(4)' > sim.log 2>&1");\r
  gSystem->ChangeDirectory("../");\r
  gSystem->Exec("mkdir BackgroundSDigits");\r
  gSystem->Exec("cp Background/*SDigits.root BackgroundSDigits");\r
  gSystem->Exec("cp BackgroundFull/galice.root BackgroundSDigits/");\r
  gSystem->Exec("cp BackgroundFull/AliESDs.root BackgroundSDigits/");\r
\r
  // Merged simulation\r
  gSystem->Setenv("CONFIG_EMBEDDING","kMerged");\r
\r
  cout<<">>>>> MERGED SIMULATION <<<<<<"<< endl;\r
  gSystem->Exec("mkdir Merged");\r
  gSystem->Exec("cp Config.C Merged/");\r
  gSystem->Exec("cp sim.C Merged/");\r
  gSystem->Exec("cp rec.C Merged/");\r
  gSystem->ChangeDirectory("Merged/");\r
  gSystem->Exec("aliroot -b -q 'sim.C(1)' > sim.log 2>&1");\r
  cout<<">>>>> MERGED RECONSTRUCTION <<<<<"<<endl;\r
  gSystem->Exec("aliroot -b -q 'rec.C(1)' > rec.log 2>&1");\r
  gSystem->ChangeDirectory("../");\r
\r
  // Pure signal re-reconstruction\r
  gSystem->Setenv("CONFIG_EMBEDDING","kSignal");\r
  \r
  cout<<">>>>> SIGNAL SIMULATION <<<<<<"<< endl;\r
  gSystem->Exec("mkdir Signal");\r
  gSystem->Exec("cp Config.C Signal/");\r
  gSystem->Exec("cp sim.C Signal/");\r
  gSystem->Exec("cp rec.C Signal/");\r
  gSystem->Exec("cp Merged/*SDigits*.root Signal/");\r
  gSystem->Exec("cp Merged/galice.root Signal/");\r
  gSystem->Exec("cp Merged/Kinematics.root Signal/");\r
  gSystem->Exec("cp -a Merged/GRP Signal/");\r
  gSystem->ChangeDirectory("Signal/");\r
  gSystem->Exec("aliroot -b -q 'sim.C(2)' > sim.log 2>&1");\r
  cout<<">>>>> SIGNAL RECONSTRUCTION <<<<<"<<endl;\r
  gSystem->Exec("aliroot -b -q 'rec.C(2)' > rec.log 2>&1");\r
  gSystem->ChangeDirectory("../");\r
  \r
  //  cout<<">>>>> TAG <<<<<"<<endl;\r
  //  gSystem->Exec("aliroot -b -q tag.C > tag.log 2>&1");\r
  //  cout<<">>>>> CHECK ESD <<<<<"<<endl;\r
  //  gSystem->Exec("aliroot -b -q CheckESD.C > check.log 2>&1");\r
  //  cout<<">>>>> AOD <<<<<"<<endl;\r
  //  gSystem->Exec("aliroot -b -q CreateAODfromESD.C > aod.log 2>&1");\r
}\r