run_GainEvolution.C

Go to the documentation of this file.

// bool run_GainEvolution(const char* runlist = "runlist_63GeV_good_and_or_ppg.list",
//                     const char* outputfile = "gainEvolution.runlist_63GeV_good_and_or_ppg.root")

bool run_GainEvolution(const char* runlist = "test.list",
                       const char* outputfile = "gainEvolution.test.list.root")
{
  gSystem->Load("libPgCalInstance.so");
  gSystem->Load("libemcCalibGFX.so");
  gSystem->Load("libemcOMpg.so");

 //  emcDataManager* dm = emcDataManager::GetInstance();
//   dm->SetVerboseLevel(1);

//   int test;
//   cout << "PID = " << gSystem->GetPid() << endl;
//   cin >> test;

  emcGainEvolution ge(outputfile);

  //  ge.setHistogrammingType(emcGainEvolution::kOneFilePerRun);
  ge.setHistogrammingType(emcGainEvolution::kOneFilePer1KRange);
  //  ge.setHistogrammingType(emcGainEvolution::kOneFileForAll);

  ifstream in(runlist);
  if (!in)
    {
      return false;
    }

  int run;

  while ( in >> run )
    {
      ge.addRun(run);
    }

  in.close();

  ge.addTrimmingPercentage(0.00);
  ge.addTrimmingPercentage(0.05);

  ge.addSector("W0");
  ge.addSector("W1");
  ge.addSector("W2");
  ge.addSector("W3");
  ge.addSector("E2");
  ge.addSector("E3");

  ge.addDetails("-3:-9:-15:92446");
  ge.addDetails("ZS:92446");
  ge.addDetails("ZS:AVOFRATIO:92446");
  ge.addDetails("ZS:AVOFASYM:92446");

  ge.addDetails("-3:-9:-15:0");
  ge.addDetails("ZS:0");

  ge.run();

  return true;
}