void myPISAtoDST(const Int_t firstEvent=6000, 
		 const Int_t maxEvents=-1, 
		 Int_t simFlag=2,
		 const char *pisaIFile="PISAEvent.root", 
		 const char *pisaOFile="simDST.root", 
		 const char *pisaRFile="rawrel.root", 
		 const Int_t forced_acc = 1,
		 const Int_t runNumber=28163)
{

  Int_t dchReqFlag=0, pc1ReqFlag=0, pc2ReqFlag=0, pc3ReqFlag=0, tofReqFlag=0, relOutFlag=0;

  Int_t writeNullFlag=0; // don't write out empty (skipped) events

  timecontrol = new TStopwatch(); timecontrol.Start();
  //getchar();

  gSystem->Load("libpreco.so"); 

  Int_t gdb_pause ;
  cin >>  gdb_pause ;
  getchar();

  recoConsts *rc = recoConsts::instance(); // instantiate flags
  setSimulationFlag(simFlag);  // simFlag = 2 sets PISA-to-DST mode in PRECO
  setMapFileFlag(2);  // using November 2001 3D magnetic field map (default was Feb. '97 2D map)
  setMapFileScale(1.0); // set map file scale factor as 1.0 

  setEvaluationFlag();
  setEMCalForcedAcceptance(forced_acc); // DO NOT use forced acceptance with HIJING simulations

// // simVertexFlag = 0 (default) means that the BBC Z0 value will be used
// // simVertexFlag = 1 means that the same simZ0Vertex value is used for all events
// // simVertexFlag = 2 means that the Z0 is taken from the PISA event header for each event
//   Int_t simVertexFlag=0;
//   Float_t simZ0Vertex=0.0, simT0Vertex=0.0;
//   Float_t simZ0VertexWidth=-9999.0, simT0VertexWidth=-9999.0;
//   if(simVertexFlag>0) {
//     setSimVertexFlag(simVertexFlag, simZ0Vertex, simT0Vertex, simZ0VertexWidth, simT0VertexWidth);
//   } // check if user wants to set the Z0 and the T0 in the BBC output

  rc->set_SysFlag("dchReqFlag", dchReqFlag);
  rc->set_SysFlag("pc1ReqFlag", pc1ReqFlag);
  rc->set_SysFlag("pc2ReqFlag", pc2ReqFlag);
  rc->set_SysFlag("pc3ReqFlag", pc3ReqFlag);
  rc->set_SysFlag("tofReqFlag", tofReqFlag);

  rc->set_SysFlag("relOutFlag", relOutFlag);
  rc->set_SysFlag("writeNullFlag", writeNullFlag);
  
  rc->set_SysFlag("MVD",0);  // 0 means the system is OFF, 1 means the system is ON
  rc->set_SysFlag("ZDC",0);
  rc->set_SysFlag("BBC",0);   
  rc->set_SysFlag("DCH",0);
  rc->set_SysFlag("PAD",1);
  rc->set_SysFlag("CRK",0);
  rc->set_SysFlag("TEC",0);
  rc->set_SysFlag("TOF",0);
  rc->set_SysFlag("EMC",1);
  rc->set_SysFlag("VTX",0);
  rc->set_SysFlag("CGL",0);
  rc->set_SysFlag("MOM",0);
  rc->set_SysFlag("MUTR",0);
  rc->set_SysFlag("MUID",0);
  rc->set_SysFlag("NTC",0);
  rc->set_SysFlag("LVL2",0);
  rc->set_SysFlag("SPIN",0);
  rc->set_SysFlag("TZR",0);
  rc->set_SysFlag("ERT",0);

  rc->set_SysFlag("HEAD",0);
  rc->set_SysFlag("TRIG",0);
  rc->set_SysFlag("T0",0);

  pisaFileOpen(pisaIFile);
  //  simSetRunNumber(-1);  // Full field, normal geometry
  simSetRunNumber(runNumber); // For run-dependent dead-map propagation 
  simDstFileOpen(pisaOFile);
  relFileOpen(pisaRFile);

  gBenchmark->Start("eventLoop");
  pisaFileRun(firstEvent,maxEvents);
  gBenchmark->Show("eventLoop");
  pisaFileEnd();

  timecontrol.Stop(); timecontrol.Print();

  return;
}