- //
- // Initialise the sub objects and stuff. Called on first event
- //
- //
- DGUARD(fDebug,1,"Initialize sub-algorithms");
- const TAxis* pe = 0;
- const TAxis* pv = 0;
-
- if (!ReadCorrections(pe,pv)) return false;
-
- fHistos.Init(*pe);
- fAODFMD.Init(*pe);
- fAODEP.Init(*pe);
- fRingSums.Init(*pe);
-
- fHData = static_cast<TH2D*>(fAODFMD.GetHistogram().Clone("d2Ndetadphi"));
- fHData->SetStats(0);
- fHData->SetDirectory(0);
- fList->Add(fHData);
-
- TList* rings = new TList;
- rings->SetName("ringSums");
- rings->SetOwner();
- fList->Add(rings);
-
- rings->Add(fRingSums.Get(1, 'I'));
- rings->Add(fRingSums.Get(2, 'I'));
- rings->Add(fRingSums.Get(2, 'O'));
- rings->Add(fRingSums.Get(3, 'I'));
- rings->Add(fRingSums.Get(3, 'O'));
- fRingSums.Get(1, 'I')->SetMarkerColor(AliForwardUtil::RingColor(1, 'I'));
- fRingSums.Get(2, 'I')->SetMarkerColor(AliForwardUtil::RingColor(2, 'I'));
- fRingSums.Get(2, 'O')->SetMarkerColor(AliForwardUtil::RingColor(2, 'O'));
- fRingSums.Get(3, 'I')->SetMarkerColor(AliForwardUtil::RingColor(3, 'I'));
- fRingSums.Get(3, 'O')->SetMarkerColor(AliForwardUtil::RingColor(3, 'O'));
-
- for(int i=1;i<=pv->GetNbins();i++)
- {
- TString nametmp=Form("vtxbin%03d",i);
- //TList* lbin= new TList();
- //lbin->SetName(nametmp.Data());
- //lbin->SetOwner();
- //fListVertexBins->Add(lbin);
- AliForwardUtil::Histos* bin=new AliForwardUtil::Histos();
- bin->Init(*pe);
- bin->Get(1, 'I')->SetName(Form("%s%s",bin->Get(1, 'I')->GetName(),nametmp.Data()));
- bin->Get(2, 'I')->SetName(Form("%s%s",bin->Get(2, 'I')->GetName(),nametmp.Data()));
- bin->Get(2, 'O')->SetName(Form("%s%s",bin->Get(2, 'O')->GetName(),nametmp.Data()));
- bin->Get(3, 'I')->SetName(Form("%s%s",bin->Get(3, 'I')->GetName(),nametmp.Data()));
- bin->Get(3, 'O')->SetName(Form("%s%s",bin->Get(3, 'O')->GetName(),nametmp.Data()));
- fList->Add(bin->Get(1, 'I'));
- fList->Add(bin->Get(2, 'I'));
- fList->Add(bin->Get(2, 'O'));
- fList->Add(bin->Get(3, 'I'));
- fList->Add(bin->Get(3, 'O'));
- fListVertexBins->Add(bin);
-
-}
-
-
- fEventInspector.SetupForData(*pv);
- fSharingFilter.SetupForData(*pe);
- fDensityCalculator.SetupForData(*pe);
- fCorrections.SetupForData(*pe);
- fHistCollector.SetupForData(*pv,*pe);
- fEventPlaneFinder.SetupForData(*pe);
-
- fFMD1icent=new TH3D("FMD1Ietavcent","FMD1ietavcent;#eta;cent",pe->GetNbins(),pe->GetXmin(),pe->GetXmax(),101,-0.5,100.5,1,0,1);
- fFMD2icent=new TH3D("FMD2Ietavcent","FMD2ietavcent;#eta;cent",pe->GetNbins(),pe->GetXmin(),pe->GetXmax(),101,-0.5,100.5,1,0,1);
- fFMD2ocent=new TH3D("FMD2Oetavcent","FMD2oetavcent;#eta;cent",pe->GetNbins(),pe->GetXmin(),pe->GetXmax(),101,-0.5,100.5,1,0,1);
- fFMD3icent=new TH3D("FMD3Ietavcent","FMD3ietavcent;#eta;cent",pe->GetNbins(),pe->GetXmin(),pe->GetXmax(),101,-0.5,100.5,1,0,1);
- fFMD3ocent=new TH3D("FMD3Oetavcent","FMD3oetavcent;#eta;cent",pe->GetNbins(),pe->GetXmin(),pe->GetXmax(),101,-0.5,100.5,1,0,1);
- fList->Add(fFMD1icent);
- fList->Add(fFMD2icent);
- fList->Add(fFMD2ocent);
- fList->Add(fFMD3icent);
- fList->Add(fFMD3ocent);
-
-
-
- this->Print();