- bug fixes in the CALO clusterizer - digits pushed from the clusterizer

[u/mrichter/AliRoot.git] / HLT / PHOS / macros / rec-hlt-calo-phos.C

//#include "AliHLTPHOSConstants.h"

//void rec_hlt_phos(const char* input="./")//, char* opt="decoder ESD")
void rec_hlt_phos()//, char* opt="decoder ESD")
{
  //  AliCDBManager::Instance()->SetDefaultStorage("raw://");
  //  AliCDBManager::Instance()->SetDefaultStorage("raw://");
  AliCDBManager::Instance()->SetRun(10);
  
  if(!gSystem->AccessPathName("galice.root")){
    cerr << "please delete the galice.root or run at different place." << endl;
    return;
  }
 
//   if (!input) {
//     cerr << "please specify input or run without arguments" << endl;
//     return;
//   }

  ///////////////////////////////////////////////////////////////////////////////////////////////////
  //
  // init the HLT system in order to define the analysis chain below
  //
  AliHLTSystem* gHLT=AliHLTPluginBase::GetInstance();
 
  ///////////////////////////////////////////////////////////////////////////////////////////////////
  //
  // define the analysis chain to be run
  //
  int moduleStart = 2;
  int moduleEnd = 2;
  int rcuStart = 0;
  int rcuEnd = 3;
  //  TString option="libAliHLTUtil.so libAliHLTRCU.so libAliHLTPHOS.so libAliHLTGlobal.so loglevel=0x7f chains=ESD-FILE";
  //TString option="libAliHLTUtil.so libAliHLTRCU.so libAliHLTPHOS.so libAliHLTGlobal.so loglevel=0x7f chains=PHS-CA_02";
  TString option="libAliHLTUtil.so libAliHLTRCU.so libAliHLTCalo.so libAliHLTPHOS.so libAliHLTGlobal.so loglevel=0x7f chains=PHS-CL_02";
  TString ecInput;
  TString emInput;
  
  for (int module = moduleStart; module <= moduleEnd; module++) 
    {
      TString clInput;
      
      for(int rcu = rcuStart; rcu <= rcuEnd; rcu++) 
	{
	  TString arg, publisher, ra, dm;
	  // raw data publisher components
	  publisher.Form("PHS-RP_%02d_%d", module, rcu);
	  arg.Form("-minid %d -datatype 'DDL_RAW ' 'PHOS'  -dataspec 0x%x ", 1792 + module*(4) + rcu, 0x1 << (module*4 + rcu));
	  AliHLTConfiguration pubConf(publisher.Data(), "AliRawReaderPublisher", NULL , arg.Data());
	  
	  // Raw analyzer
	  arg = "";
	  ra.Form("PHS-RA_%02d_%d", module, rcu);
	  //	  AliHLTConfiguration rawConf(ra.Data(), "PhosRawCrudev2", publisher.Data(), arg.Data());
	  AliHLTConfiguration rawConf(ra.Data(), "PhosRawCrude", publisher.Data(), arg.Data());
	  
	  // digit maker components
	  dm.Form("PHS-DM_%02d_%d", module, rcu);
	  arg="";
	  arg.Form("-sethighgainfactor 0.005 -setlowgainfactor 0.08 -setdigitthresholds 0.005 0.002");
	  AliHLTConfiguration dmConf(dm.Data(), "PhosDigitMaker", ra.Data(), arg.Data());

	  if(clInput.Length() > 0) clInput += " ";
	  clInput+=dm;
	}
        TString arg, cl, ca;

      cl.Form("PHS-CL_%02d", module);
      arg = "";
      arg.Form("-digitthreshold 0.005 -recpointthreshold 0.1 -modulemode");
      AliHLTConfiguration clConf(cl.Data(), "PhosClusterizer", clInput.Data(), arg.Data());
// 	
//       ca.Form("PHS-CA_%02d", module);
//       arg = "";
//       AliHLTConfiguration caConf(ca.Data(), "PhosClusterAnalyser", cl.Data(), arg.Data());

      if(ecInput.Length() > 0) ecInput += " ";
      ecInput += ca;
    }
      
  emInput = ecInput;

  TString arg, ec, em, hp, ef;
  

//   ec.Form("ESD-CONVERTER");
//   arg = "";
// 
//   //  AliHLTConfiguration emConf(em.Data(), "PhosEsdEntriesMaker", ecInput.Data(), "");
// 
//   AliHLTConfiguration esdcconf(ec.Data(), "GlobalEsdConverter"   , ecInput.Data(), "");
//   
//   // the root file writer configuration
//   ef.Form("ESD-FILE");
//   AliHLTConfiguration sink(ef.Data(), "EsdCollector", ec.Data(), "-directory hlt-phos-esd");
// 

  ///////////////////////////////////////////////////////////////////////////////////////////////////
  //
  // Init and run the reconstruction
  // All but HLT reconstructio is switched off 
  //
  AliReconstruction rec;
  //  rec.SetInput(input);
  rec.SetInput("./");
  rec.SetRunVertexFinder(kFALSE);
  rec.SetRunReconstruction("HLT");
  rec.SetRunTracking(":");
  rec.SetLoadAlignFromCDB(0);
  rec.SetRunQA(":");

  // NOTE: FillESD is a step in the AliReconstruction sequence and has
  // nothing to do with the fact that this macro writes ESD output
  // HLT processes the HLTOUT during FillESD and extracts data which
  // has already been prepared. This step is currently not necessary for
  // this macro
  //  rec.SetRunLocalReconstruction("PHOS") ;
    rec.SetOption("HLT", option);
  //  rec.SetOption("HLT", "libAliHLTUtil.so libAliHLTRCU.so libAliHLTPHOS.so loglevel=0x7f chains=ESD_WRITER" )
    rec.SetSpecificStorage("GRP/GRP/Data", Form("local://%s", gSystem->pwd()));
  rec.Run();
}