Config for charged + neutral

JETAN/ConfigJetAnalysisFastJet.C

@@ -0,0 +1,53 @@
+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;
+}