AliJetFinder* ConfigJetAnalysis() { // // Configuration goes here // printf("ConfigJetAnalysis() \n"); // Define the grids AliJetGrid *grid = new AliJetGrid(419,119,0.,2*TMath::Pi(),-0.9,0.9); grid->SetGridType(1); grid->InitParams(80.*TMath::Pi()/180,190.*TMath::Pi()/180,-0.7,0.7); grid->SetMatrixIndexes(); grid->SetIndexIJ(); AliJetGrid *grid2 = new AliJetGrid(131,95,80.*TMath::Pi()/180.,190.*TMath::Pi()/180.,-0.7,0.7); grid2->SetGridType(0); grid2->SetMatrixIndexes(); grid2->SetIndexIJ(); // Define ESD reader header AliJetESDReaderHeader *jrh = new AliJetESDReaderHeader(); jrh->SetComment("Testing"); jrh->SetReadSignalOnly(kFALSE); // Detector options: 0 = Charged particles only (MomentumArray) // 1 = Charged particles only (UnitArray) // 2 = Neutral cells only (UnitArray) // 3 = Charged particles + neutral cells (UnitArray) jrh->SetDetector(0); jrh->SetDebug(0); jrh->SetFiducialEta(-0.9,0.9); jrh->SetFiducialPhi(0,2*TMath::Pi()); // Define reader and set its header AliJetESDReader *er = new AliJetESDReader(); er->SetReaderHeader(jrh); er->SetTPCGrid(grid); er->SetEMCalGrid(grid2); // Define jet header AliFastJetHeaderV1 *jh=new AliFastJetHeaderV1(); jh->SetComment("Fast jet code with default parameters"); jh->SetDebug(0); jh->SetBGMode(0); jh->SetRparam(1.0); // setup parameters jh->SetPtMin(20.); // Define jet finder. Set its header and reader jetFinder = new AliFastJetFinder(); jetFinder->SetJetHeader(jh); jetFinder->SetJetReader(er); jetFinder->SetPlotMode(kTRUE); return jetFinder; }