PWG1/RunAlignmentDataFilterITS.C

   1 void RunAlignmentDataFilterITS() {
   2   //
   3   // Macro to extract AliTrackPoints for ITS
   4   // A.Dainese, andrea.dainese@pd.infn.it
   5   //
   6
   7   // Input
   8   Bool_t singlefile=kFALSE;
   9   //TString esdpath="/home/dainesea/alignData/RAWdata_CosmicsSum09/RecoSPD/chunk.";
  10   TString esdpath="/home/dainesea/alignData/RAWdata_CosmicsSum09/RecoITS_B_mille_SPD_SDDSSDsurvey_SSDHLayer_th50_130709/chunk.";
  11   Int_t ifirst=1, ilast=11;
  12   //
  13   Int_t nentries=1234567890;
  14   Int_t firstentry=0;
  15
  16
  17   // Load PWG1 library
  18   gSystem->Load("libANALYSIS.so");
  19   gSystem->Load("libANALYSISalice.so");
  20   gSystem->Load("libPWG1.so");
  21
  22
  23   // Create the analysis manager
  24   AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis");
  25
  26   // Create the task
  27   AliAlignmentDataFilterITS *taskFilter = new AliAlignmentDataFilterITS("filterITS");
  28   AliLog::SetClassDebugLevel("AliAlignmentDataFilterITS",10);
  29   // configuration via AliITSRecoParam (should be taken from OCDB)
  30   AliITSRecoParam *itsRecoParam = AliITSRecoParam::GetCosmicTestParam();
  31   itsRecoParam->SetAlignFilterUseLayer(0,kTRUE);
  32   itsRecoParam->SetAlignFilterUseLayer(1,kTRUE);
  33   itsRecoParam->SetAlignFilterUseLayer(2,kFALSE);
  34   itsRecoParam->SetAlignFilterUseLayer(3,kFALSE);
  35   itsRecoParam->SetAlignFilterUseLayer(4,kFALSE);
  36   itsRecoParam->SetAlignFilterUseLayer(5,kFALSE);
  37   taskFilter->SetITSRecoParam(itsRecoParam);
  38   taskFilter->SetOnlySPDFO();
  39
  40   // Add ESD handler
  41   AliESDInputHandler *esdH = new AliESDInputHandler();
  42   mgr->SetInputEventHandler(esdH);
  43
  44   TChain *chainESD = new TChain("esdTree");
  45   if(singlefile) {
  46     chainESD->Add("AliESDs.root");
  47   } else {
  48     for(Int_t i=ifirst; i<=ilast; i++) {
  49       TString esdfile=esdpath; esdfile+=i; esdfile.Append("/AliESDs.root");
  50       chainESD->Add(esdfile.Data());
  51     }
  52   }
  53
  54   // Attach input
  55   cInput = mgr->CreateContainer("cInput",TChain::Class(),AliAnalysisManager::kInputContainer);
  56   //mgr->ConnectInput(taskFilter, 0, cInput); // v4-16-Release
  57   mgr->ConnectInput(taskFilter,0,mgr->GetCommonInputContainer());
  58
  59   // Attach output
  60   cOutput0= mgr->CreateContainer("cOutput0",TTree::Class(),AliAnalysisManager::kOutputContainer,"AliTrackPoints.root");
  61   mgr->ConnectOutput(taskFilter,0,cOutput0);
  62   cOutput1= mgr->CreateContainer("cOutput1",TList::Class(),AliAnalysisManager::kOutputContainer,"AliTrackPoints.root");
  63   mgr->ConnectOutput(taskFilter,1,cOutput1);
  64
  65   // Enable debug printouts
  66   mgr->SetDebugLevel(10);
  67
  68   // Run analysis
  69   mgr->InitAnalysis();
  70   mgr->PrintStatus();
  71   mgr->StartAnalysis("local",chainESD,nentries,firstentry);
  72
  73  return;
  74
  75 }