// ------------------------------------------
TString publisherId("ESDMCEventPublisher");
// ------------------------------------------
-
- TString publisherArg( Form("-entrytype MCFAST -dataspec 0x0000001F -datapath /home/jthaeder/jet/data/HEAD_2009-10-26/FastGen/kPythia6Jets125_150_14TeV/JET-ETA=-0.2,0.2_JET-ET=50,1000_R=0.4_200ev") );
+ // MCFAST - MC - ESD
+ TString publisherArg("-entrytype ESD -dataspec 0x0000001F -datapath /lustre/alice/jthaeder/data/HEAD_2010-07-09/7TeV/pp_Perugia0/014000");
+
AliHLTConfiguration ESDMCEventPublisher(publisherId.Data(), "ESDMCEventPublisher", NULL, publisherArg.Data() );
if (!analysisInput.IsNull()) analysisInput+=" ";
// -- Processing Components -- //
// - - //
// ----------------------------//
+
#if 1
+
// ------------------------------------------
// -- ConeJetFinder
// ------------------------------------------
TString jetId("JETConeJet");
// ------------------------------------------
- TString jetArg( Form("-algorithm FSCSquareCell -leading 1 -coneRadius %.1f -trackCutMinPt 0.0 -seedCutMinPt %.1f -jetCutMinEt %.1f",
- coneRadius, cutPtSeed, cutEtJet) );
+ TString jetArg(
+ Form("-algorithm FSCSquareCell -leading 1 -coneRadius %.1f -trackCutMinPt 0.0 -seedCutMinPt %.1f -jetCutMinEt %.1f", coneRadius, cutPtSeed, cutEtJet)
+ );
AliHLTConfiguration jetCone(jetId.Data(), "JETConeJetFinder", jetInput.Data(), jetArg.Data());
analysisInput += jetId;
#else
+
// ------------------------------------------
// -- FastJetFinder
// ------------------------------------------
- AliHLTConfiguration jetFinder("JETFastJet_Kt", "JETFastJetFinder",
- jetInput.Data(),"-finderType kt");
+ AliHLTConfiguration jetFinder("JETFastJet_Kt", "JETFastJetFinder", jetInput.Data(),"-finderType kt");
if (!writerInput.IsNull()) writerInput+=" ";
writerInput+="JETFastJet_Kt";
-
- AliHLTConfiguration jetFinder("JETFastJet_AntiKt", "JETFastJetFinder",
- jetInput.Data(),"-finderType antikt");
+ AliHLTConfiguration jetFinder("JETFastJet_AntiKt", "JETFastJetFinder", jetInput.Data(),"-finderType antikt");
if (!writerInput.IsNull()) writerInput+=" ";
writerInput+="JETFastJet_AntiKt";
#endif
-
+
// ------------------------------------------
// -- Jet Analysis
// ------------------------------------------
TString writerArg( Form("-directory analysis -datafile analyze_%d_%s -write-all-events", nEvents, pprRunName[runType] ));
// -- The RootFileWriter
- AliHLTConfiguration rootWriter("RootWriter", "ROOTFileWriter",
- writerInput.Data(), writerArg.Data() );
+ AliHLTConfiguration rootWriter("RootWriter", "ROOTFileWriter", writerInput.Data(), writerArg.Data() );
// --------------------------- //
// - - //
AliHLTJETAnalysisMerge* merge = new AliHLTJETAnalysisMerge();
merge->Initialize();
-
+
merge->AddJets( jets );
merge->CreateCanvas();
#/bin/bash
+NEVENTS=200
+Pt_BIN=kPythia6Jets104_125
+
if [ ! -d ./analysis ] ; then
mkdir analysis
fi
if [ -n "${ROOTFILES}" ] ; then
rm ${ROOTFILES}
fi
+
-aliroot -b -l -q 'HLTJetReconstruction.C(20,kTRUE)' 2>&1 | tee ChainLog.log
+aliroot -b -l -q HLTJetReconstruction.C'('${NEVENTS}', kFALSE, '${Pt_BIN}')' 2>&1 | tee ChainLog.log
ROOTFILES=`find . -maxdepth 1 -name "*.root" `
rm ${ROOTFILES}
fi
-aliroot -l 'readJets.C("./analysis/EOR_analyze_20_kPythia6Jets104_125.root")'
+aliroot -l -b -q readJets.C'("./analysis/EOR_analyze_'${NEVENTS}'_'${Pt_BIN}'.root")'
-#valgrind --error-limit=no --leak-check=full --show-reachable=yes aliroot -b -l -q 'HLTJetReconstruction.C(10)'
+# valgrind --error-limit=no --leak-check=full --show-reachable=yes aliroot -b -l -q 'HLTJetReconstruction.C(10)'
#/bin/bash
-if [ ! -d ./analysis ] ; then
- mkdir analysis
+if [ "$1" == "MC" ] ; then
+ MODE=MC
+elif [ "$1" == "ESD" ] ; then
+ MODE=ESD
+else
+ MODE=MC
fi
ROOTFILES=`find . -maxdepth 1 -name "*.root" `
rm ${ROOTFILES}
fi
-aliroot -b -l -q './tasks/JetAnalysisManagerHLT.C' 2>&1 | tee TaskLog.log
+if [ "$MODE" == "ESD" ] ; then
+ echo " -= Load : JetAnalysisManagerHLT.C =- "
+ aliroot -b -l -q './tasks/JetAnalysisManagerHLT.C' 2>&1 | tee TaskLog.log
+else
+ echo " -= Load : JetAnalysisManagerHLTMC.C =- "
+ aliroot -b -l -q './tasks/JetAnalysisManagerHLTMC.C' 2>&1 | tee TaskLog.log
+fi
ROOTFILES=`find . -maxdepth 1 -name "*.root" `
/** ---------------------------------------------------------------------
* @file ConfigJetAnalysisHLT.C
- * @author Jochen Thaeder <thaeder@kip.uni-heidelberg.de>
+ * @author Jochen Thaeder <jochen@thaeder.de>
* @brief Run HLT cone finder in analysis framework, on ESDs
*
* ---------------------------------------------------------------------
jetReaderHeader->SetJetAlgorithm(algorithm);
// Set prt to track cuts
- jetReaderHeader->SetTrackCuts( trackCuts );
- jetReaderHeader->SetSeedCuts( seedCuts );
+ jetReaderHeader->SetTrackCuts(trackCuts);
+ jetReaderHeader->SetSeedCuts(seedCuts);
// Set Eta min/max and Phi min/max
jetReaderHeader->SetFiducialEta( -0.9, 0.9) ;
/** ---------------------------------------------------------------------
* @file ConfigJetAnalysisHLTMC.C
- * @author Jochen Thaeder <thaeder@kip.uni-heidelberg.de>
+ * @author Jochen Thaeder <jochen@thaeder.de>
* @brief Run HLT cone finder in analysis framework, on Kinematics
*
* ---------------------------------------------------------------------
jetReaderHeader->SetJetAlgorithm(algorithm);
// Set prt to track cuts
- jetReaderHeader->SetTrackCuts( trackCuts );
- jetReaderHeader->SetSeedCuts( seedCuts );
+ jetReaderHeader->SetTrackCuts(trackCuts);
+ jetReaderHeader->SetSeedCuts(seedCuts);
// Set Eta min/max and Phi min/max
jetReaderHeader->SetFiducialEta( -0.9, 0.9) ;
//
// --------------------------------------------------------------------------------
- gROOT->LoadMacro("$ALICE_ROOT/PWG0/CreateESDChain.C");
+ gROOT->LoadMacro("${ALICE_ROOT}/PWG0/CreateESDChain.C");
TChain* chain = new TChain("esdTree");
- chain->Add("~/jet/data/HEAD_2009-06-04/kPythia6Jets104_125_14TeV/JET-ETA=-0.2,0.2_JET-ET=10,1000_R=0.7_100ev/AliESDs.root");
+ chain->Add("~/jet/data/HEAD_2010-01-08/Gen/kPythia6Jets86_104_14TeV/JET-ETA=-0.2,0.2_JET-ET=10,1000_R=0.4_20ev/AliESDs.root");
+
+
+ // chain->Add("~/jet/data/HEAD_2009-06-04/kPythia6Jets104_125_14TeV/JET-ETA=-0.2,0.2_JET-ET=10,1000_R=0.7_100ev/AliESDs.root");
//chain->Add("~/jet/data/HEAD_2009-06-04/kPythia6Jets104_125_14TeV/JET-ETA=-0.2,0.2_JET-ET=10,1000_R=0.7_1ev/AliESDs.root");
// --------------------------------------------------------------------------------
// MC Truth
AliMCEventHandler* mcHandler = new AliMCEventHandler();
-
+
AliAnalysisManager *mgr = new AliAnalysisManager("Jet Manager", "Jet Manager");
mgr->SetInputEventHandler (inpHandler);
mgr->SetOutputEventHandler (aodHandler);
mgr->SetMCtruthEventHandler(mcHandler);
- mgr->SetDebugLevel(0);
+ mgr->SetDebugLevel(10);
// --------------------------------------------------------------------------------
//
jetana->SetConfigFile("./tasks/ConfigJetAnalysisHLT.C");
jetana->SetNonStdBranch("jetsHLT");
jetana->SetDebugLevel(10);
- mgr->AddTask(jetana);
-
- // -- HLT FFSC - MC
- AliAnalysisTaskJets *jetanaMC = new AliAnalysisTaskJets("JetAnalysisHLTMC");
- jetanaMC->SetConfigFile("./tasks/ConfigJetAnalysisHLTMC.C");
- jetanaMC->SetNonStdBranch("jetsHLTMC");
- jetanaMC->SetDebugLevel(0);
- mgr->AddTask(jetanaMC);
+ mgr->AddTask(jetana);
// --------------------------------------------------------------------------------
//
AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(), AliAnalysisManager::kOutputContainer, "histos.root");
- AliAnalysisDataContainer *coutputMC2 = mgr->CreateContainer("histosMC", TList::Class(), AliAnalysisManager::kOutputContainer, "histosMC.root");
-
+
mgr->ConnectInput (jetana, 0, cinput1 );
mgr->ConnectOutput (jetana, 0, coutput1 );
mgr->ConnectOutput (jetana, 1, coutput2 );
-
- mgr->ConnectInput (jetanaMC, 0, cinput1 );
- mgr->ConnectOutput (jetanaMC, 0, coutput1 );
- mgr->ConnectOutput (jetanaMC, 1, coutputMC2 );
-
+
+
// --------------------------------------------------------------------------------
//
// Run the analysis
//
// --------------------------------------------------------------------------------
-
+
mgr->InitAnalysis();
mgr->PrintStatus();
mgr->StartAnalysis("local",chain);
+++ /dev/null
-void JetAnalysisManagerHLTKine()
-{
- gSystem->Load("libCGAL.so");
- gSystem->Load("${FASTJET}/lib/libfastjet.so");
-
- gSystem->Load("libTree.so");
- gSystem->Load("libPhysics.so");
- gSystem->Load("libGeom.so");
- gSystem->Load("libVMC.so");
- gSystem->Load("libANALYSIS.so");
- gSystem->Load("libSTEERBase.so");
- gSystem->Load("libAOD.so");
- gSystem->Load("libESD.so");
- gSystem->Load("libANALYSISalice.so");
- gSystem->Load("libJETAN.so");
-
- //
- if (gApplication) gApplication->InitializeGraphics();
- // Create the chain
- //
-
- gROOT->LoadMacro("${ALICE_ROOT}/PWG0/CreateESDChain.C");
- TChain* chain = new TChain("TE");
- chain->Add("/home/jthaeder/jet/data/HEAD_2009-03-17/FastGen/kPythia6Jets104_125_14TeV/JET-ETA=-0.2,0.2_JET-ET=50,1000_R=0.4_10ev/galice.root");
-
- ///////////////////////////////////////////////////////////////////////////////////
- // Create the analysis manager
- //
- // Output
- AliAODHandler* aodHandler = new AliAODHandler();
- aodHandler->SetOutputFileName("aod.root");
- // MC Truth
- AliMCEventHandler* mcHandler = new AliMCEventHandler();
- mcHandler->SetReadTR(kFALSE);
-
- AliAnalysisManager *mgr = new AliAnalysisManager("Jet Manager", "Jet Manager");
- mgr->SetOutputEventHandler (aodHandler);
- mgr->SetMCtruthEventHandler(mcHandler);
- mgr->SetDebugLevel(10);
-
- AliAnalysisTaskKineFilter *kinefilter = new AliAnalysisTaskKineFilter("Kine Filter");
- mgr->AddTask(kinefilter);
-
-
- AliAnalysisTaskJets *jetana = new AliAnalysisTaskJets("JetAnalysis");
- //jetana->SetConfigFile("ConfigJetAnalysisFastjetMC.C");
- jetana->SetConfigFile("${ALICE_ROOT}/HLT/JET/macros/tasks/ConfigJetAnalysisHLTMC.C");
- jetana->SetDebugLevel(10);
- mgr->AddTask(jetana);
-
- //
- // Create containers for input/output
- AliAnalysisDataContainer *cinput1 = mgr->CreateContainer("cchain",TChain::Class(),
- AliAnalysisManager::kInputContainer);
-
- AliAnalysisDataContainer *coutput1 = mgr->CreateContainer("tree", TTree::Class(),
- AliAnalysisManager::kOutputContainer, "default");
-
- AliAnalysisDataContainer *coutput2 = mgr->CreateContainer("histos", TList::Class(),
- AliAnalysisManager::kOutputContainer, "histos.root");
-
- mgr->ConnectInput (kinefilter, 0, cinput1 );
- mgr->ConnectOutput (kinefilter, 0, coutput1 );
-
- mgr->ConnectInput (jetana, 0, cinput1 );
- mgr->ConnectOutput (jetana, 0, coutput1 );
- mgr->ConnectOutput (jetana, 1, coutput2 );
-
- //
- // Run the analysis
- //
-
- mgr->InitAnalysis();
- mgr->PrintStatus();
- mgr->StartAnalysis("local",chain);
-}
--- /dev/null
+void JetAnalysisManagerHLTMC() {
+
+ if ( getenv("FASTJET") ) {
+ gSystem->Load("libCGAL.so");
+ gSystem->Load("libfastjet.so");
+ gSystem->Load("libSISConePlugin.so");
+ }
+
+ gSystem->Load("libTree.so");
+ gSystem->Load("libPhysics.so");
+ gSystem->Load("libGeom.so");
+ gSystem->Load("libVMC.so");
+ gSystem->Load("libANALYSIS.so");
+ gSystem->Load("libSTEERBase.so");
+ gSystem->Load("libAOD.so");
+ gSystem->Load("libESD.so");
+ gSystem->Load("libANALYSISalice.so");
+ gSystem->Load("libJETAN.so");
+
+ gSystem->Load("libHLTbase.so");
+ gSystem->Load("libAliHLTUtil.so");
+ gSystem->Load("libAliHLTJET.so");
+
+ Int_t debugLevel = 10;
+
+ // --------------------------------------------------------------------------------
+
+ if (gApplication) gApplication->InitializeGraphics();
+
+ // --------------------------------------------------------------------------------
+ //
+ // Create the chain
+ //
+ // --------------------------------------------------------------------------------
+
+ gROOT->LoadMacro("${ALICE_ROOT}/PWG0/CreateESDChain.C");
+ TChain* chain = new TChain("TE");
+
+ chain->Add("/home/jthaeder/jet/data/HEAD_2009-10-26/FastGen/kPythia6Jets125_150_14TeV/JET-ETA=-0.2,0.2_JET-ET=50,1000_R=0.4_200ev/galice.root");
+
+ // --------------------------------------------------------------------------------
+ //
+ // Create the analysis manager
+ //
+ // --------------------------------------------------------------------------------
+
+ // Output
+ AliAODHandler* aodHandler = new AliAODHandler();
+ aodHandler->SetOutputFileName("aod.root");
+
+ // MC Truth
+ AliMCEventHandler* mcHandler = new AliMCEventHandler();
+ mcHandler->SetReadTR(kFALSE); // for fastgen
+
+ AliAnalysisManager *mgr = new AliAnalysisManager("Jet Manager", "Jet Manager");
+ mgr->SetOutputEventHandler (aodHandler);
+ mgr->SetMCtruthEventHandler(mcHandler);
+ mgr->SetDebugLevel(debugLevel);
+
+ // --------------------------------------------------------------------------------
+ //
+ // Set Configfiles
+ //
+ // --------------------------------------------------------------------------------
+
+ // -- HLT FFSC - MC
+#if FFSC
+ AliAnalysisTaskJets *taskFFSC = new AliAnalysisTaskJets("JetAnalysisHLTMC_FFSC");
+ taskFFSC->SetConfigFile("./tasks/ConfigJetAnalysisHLTMC.C");
+ taskFFSC->SetNonStdBranch("jetsHLTMC_FFSC");
+ taskFFSC->SetDebugLevel(debugLevel);
+ mgr->AddTask(taskFFSC);
+#endif
+
+ // -- HLT FastJet - Kt
+#if FASTJET_KT
+ if ( getenv("FASTJET") ) {
+ AliAnalysisTaskJets *taskKt = new AliAnalysisTaskJets("JetAnalysisHLTMC_Kt");
+ taskKt->SetConfigFile("./tasks/ConfigJetAnalysisHLTMCKt.C");
+ taskKt->SetNonStdBranch("jetsHLTMC_Kt");
+ taskKt->SetDebugLevel(debugLevel);
+ mgr->AddTask(taskKt);
+ }
+#endif
+
+#if FASTJET_ANTIKT
+ // -- HLT FastJet - AntiKt
+ if ( getenv("FASTJET") ) {
+ AliAnalysisTaskJets *taskAntiKt = new AliAnalysisTaskJets("JetAnalysisHLTMC_AntiKt");
+ taskAntiKt->SetConfigFile("./tasks/ConfigJetAnalysisHLTMCAntiKt.C");
+ taskAntiKt->SetNonStdBranch("jetsHLTMC_AntiKt");
+ taskAntiKt->SetDebugLevel(debugLevel);
+ mgr->AddTask(taskAntiKt);
+ }
+#endif
+
+ // --------------------------------------------------------------------------------
+ //
+ // Create containers for input/output
+ //
+ // --------------------------------------------------------------------------------
+
+ AliAnalysisDataContainer *cinput1 = mgr->GetCommonInputContainer();
+ AliAnalysisDataContainer *coutput1 = mgr->GetCommonOutputContainer();
+
+#if FFSC
+ AliAnalysisDataContainer *coutputFFSC = mgr->CreateContainer("histos", TList::Class(),
+ AliAnalysisManager::kOutputContainer, "histos_MC_FFSC.root");
+ mgr->ConnectInput (taskFFSC, 0, cinput1 );
+ mgr->ConnectOutput (taskFFSC, 0, coutput1 );
+ mgr->ConnectOutput (taskFFSC, 1, coutputFFSC );
+#endif
+
+#if FASTJET_KT
+ if ( getenv("FASTJET") ) {
+ AliAnalysisDataContainer *coutputKt = mgr->CreateContainer("histos", TList::Class(),
+ AliAnalysisManager::kOutputContainer, "histos_MC_KT.root");
+ mgr->ConnectInput (taskKt, 0, cinput1 );
+ mgr->ConnectOutput (taskKt, 0, coutput1 );
+ mgr->ConnectOutput (taskKt, 1, coutputKt );
+ }
+#endif
+
+#if FASTJET_ANTIKT
+ if ( getenv("FASTJET") ) {
+ AliAnalysisDataContainer *coutputAntiKt = mgr->CreateContainer("histos", TList::Class(),
+ AliAnalysisManager::kOutputContainer, "histos_MC_ANTIKT.root");
+ mgr->ConnectInput (taskAntiKt, 0, cinput1 );
+ mgr->ConnectOutput (taskAntiKt, 0, coutput1 );
+ mgr->ConnectOutput (taskAntiKt, 1, coutputAntiKt );
+ }
+#endif
+
+ // AliAnalysisTaskKineFilter *kinefilter = new AliAnalysisTaskKineFilter("Kine Filter");
+ // mgr->AddTask(kinefilter);
+ // mgr->ConnectInput (kinefilter, 0, cinput1 );
+ // mgr->ConnectOutput (kinefilter, 0, coutput1 );
+
+ // --------------------------------------------------------------------------------
+ //
+ // Run the analysis
+ //
+ // --------------------------------------------------------------------------------
+
+ mgr->InitAnalysis();
+ mgr->PrintStatus();
+ mgr->StartAnalysis("local",chain);
+}