1 void RunAlignmentDataFilterITS() {
3 // Macro to extract AliTrackPoints for ITS
4 // A.Dainese, andrea.dainese@pd.infn.it
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;
13 Int_t nentries=1234567890;
18 gSystem->Load("libANALYSIS.so");
19 gSystem->Load("libANALYSISalice.so");
20 gSystem->Load("libPWG1.so");
23 // Create the analysis manager
24 AliAnalysisManager *mgr = new AliAnalysisManager("testAnalysis");
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();
41 AliESDInputHandler *esdH = new AliESDInputHandler();
42 mgr->SetInputEventHandler(esdH);
44 TChain *chainESD = new TChain("esdTree");
46 chainESD->Add("AliESDs.root");
48 for(Int_t i=ifirst; i<=ilast; i++) {
49 TString esdfile=esdpath; esdfile+=i; esdfile.Append("/AliESDs.root");
50 chainESD->Add(esdfile.Data());
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());
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);
65 // Enable debug printouts
66 mgr->SetDebugLevel(10);
71 mgr->StartAnalysis("local",chainESD,nentries,firstentry);