HBTAN/AliHBTWriteInternFormat.C

   1 #if 0
   2   #include "/afs/cern.ch/user/s/skowron/aliroot/my/TPC/alles.h"
   3   #include "AliHBTReader.h"
   4   #include "AliHBTReaderKineTree.h"
   5   #include "AliHBTReaderITSv2.h"
   6   #include "AliHBTReaderITSv1.h"
   7   #include "AliHBTReaderTPC.h"
   8   #include "AliHBTParticleCut.h"
   9   #include "AliHBTEvent.h"
  10   #include "AliHBTPairCut.h"
  11   #include "AliHBTQResolutionFctns.h"
  12   #include "AliHBTTwoTrackEffFctn.h"
  13   #include "AliHBTCorrelFctn.h"
  14   #include "TSystem.h"
  15   #include "TObjString.h"
  16   #include "TString.h"
  17   #include "AliPDG.h"
  18 #endif
  19
  20
  21 void AliHBTWriteInternFormat(Option_t* datatype, Option_t* processopt="TracksAndParticles",
  22                 Int_t first = -1,Int_t last = -1,
  23                 char *outfile = "data.root")
  24  {
  25 //HBT Anlysis Macro
  26 //Anlyzes TPC recontructed tracks and simulated particles that corresponds to them
  27
  28 //datatype defines type of data to be read
  29 //  Kine  - analyzes Kine Tree: simulated particles
  30 //  TPC   - analyzes TPC   tracking + particles corresponding to these tracks
  31 //  ITSv1 - analyzes ITSv1 ----------------------//--------------------------
  32 //  ITSv2 - analyzes ITSv2 ----------------------//--------------------------
  33
  34 //processopt defines option passed to AliHBTAnlysis::Process method
  35 // default: TracksAndParticles - process both recontructed tracks and sim. particles corresponding to them
  36 //          Tracks - process only recontructed tracks
  37 //          Particles - process only simulated particles
  38
  39 //Reads data from diroctories from first to last(including)
  40 // For examples if first=3 and last=5 it reads from
  41 //  ./3/
  42 //  ./4/
  43 //  ./5/
  44 //if first or last is negative (or both), it reads from current directory
  45 //
  46 //these names I use when analysis is done directly from CASTOR files via RFIO
  47   const char* basedir="rfio:/castor/cern.ch/user/s/skowron";
  48   const char* serie="10k.12fm";
  49   const char* field = "0.2";
  50
  51   // Dynamically link some shared libs
  52
  53 //  gROOT->LoadMacro("loadlibs.C");
  54 //  loadlibs();
  55
  56   gSystem->Load("$(ALICE_ROOT)/lib/tgt_$(ALICE_TARGET)/libHBTAN");
  57
  58
  59   /***********************************************************/
  60
  61   AliHBTReader* reader;
  62   Int_t kine = strcmp(datatype,"Kine");
  63   Int_t TPC = strcmp(datatype,"TPC");
  64   Int_t ITSv1 = strcmp(datatype,"ITSv1");
  65   Int_t ITSv2 = strcmp(datatype,"ITSv2");
  66
  67   if(!kine)
  68    {
  69     reader = new AliHBTReaderKineTree();
  70     processopt="Particles"; //this reader by definition reads only simulated particles
  71    }
  72   else if(!TPC)
  73    {
  74     reader = new AliHBTReaderTPC();
  75    }
  76   else if(!ITSv1)
  77    {
  78     reader = new AliHBTReaderITSv1();
  79    }
  80   else if(!ITSv2)
  81    {
  82     reader = new AliHBTReaderITSv2();
  83    }
  84   else
  85    {
  86     cerr<<"Option "<<datatype<<"  not recognized. Exiting"<<endl;
  87     return;
  88    }
  89
  90   TObjArray* dirs=0;
  91   if ( (first >= 0) && (last>=0) && ( (last-first)>=0 ) )
  92    {//read from many dirs dirs
  93      char buff[50];
  94      dirs = new TObjArray(last-first+1);
  95      for (Int_t i = first; i<=last; i++)
  96       {
  97         sprintf(buff,"%s/%s/%s/%d",basedir,field,serie,i);
  98         TObjString *odir= new TObjString(buff);
  99         dirs->Add(odir);
 100       }
 101     }
 102    reader->SetDirs(dirs);
 103    AliHBTReaderInternal::Write(reader,outfile);
 104
 105  }
 106