#include "AliMCEvent.h"
#include "AliHeader.h"
+#include "AliStack.h"
#include "TString.h"
#include "TObjString.h"
// see header file for class documentation
if ( fJetFinder || fJetReader || fJetHeader || fJetReader )
- return EINPROGRESS;
+ return -EINPROGRESS;
-#ifdef HAVE_FASTJET
- printf("HAVE FASTJET");
-#else
- printf("NO FASTJET");
-#endif
-#if 1
- fJetReaderHeader = new AliJetKineReaderHeader();
+ // -- Jet Reader Header
+ // -------------------------------------------
+ if ( ! (fJetReaderHeader = new AliJetKineReaderHeader()) ) {
+ HLTError("Error initializing Jet Reader Header");
+ return -EINPROGRESS;
+ }
+
fJetReaderHeader->SetComment("MC full Kinematics");
fJetReaderHeader->SetFastSimTPC(kFALSE);
fJetReaderHeader->SetFastSimEMCAL(kFALSE);
fJetReaderHeader->SetPtCut(0.);
-
- // Define reader and set its header
- fJetReader = new AliJetKineReader();
+
+ // -- Jet Reader
+ // -------------------------------------------
+ if ( ! (fJetReader = new AliJetKineReader()) ) {
+ HLTError("Error initializing Jet Reader");
+ return -EINPROGRESS;
+ }
+
fJetReader->SetReaderHeader(fJetReaderHeader);
-#else
+ // -- Jet Header
+ // -------------------------------------------
+ if ( ! (fJetHeader = new AliFastJetHeader()) ) {
+ HLTError("Error initializing Jet Header");
+ return -EINPROGRESS;
+ }
+
+ fJetHeader->SetRparam(0.7);
+
+ // -- Jet Finder
+ // -------------------------------------------
+ if ( ! (fJetFinder = new AliFastJetFinder()) ) {
+ HLTError("Error initializing Jet Finder");
+ return -EINPROGRESS;
+ }
+
+ fJetFinder->SetJetHeader(fJetHeader);
+ fJetFinder->SetJetReader(fJetReader);
+ fJetFinder->SetOutputFile("jets.root");
+
+ // -- Initialize Jet Finder
+ // -------------------------------------------
+ fJetFinder->Init();
+
+
+#if 0
fJetReaderHeader = new AliJetESDReaderHeader();
fJetReaderHeader->SetComment("Testing");
fJetReaderHeader->SetFirstEvent(0);
fJetReaderHeader->SetLastEvent(4);
-
fJetReader = new AliJetESDReader();
fJetReader->SetReaderHeader(fJetReaderHeader);
#endif
-
- fJetHeader = new AliFastJetHeader();
- fJetHeader->SetRparam(0.7);
-
- fJetFinder = new AliFastJetFinder();
- fJetFinder->SetJetHeader(fJetHeader);
- fJetFinder->SetJetReader(fJetReader);
- fJetFinder->SetOutputFile("jets.root");
return 0;
}
Int_t AliHLTJETFastJetComponent::DoDeinit() {
// see header file for class documentation
- /*
if ( fJetFinder )
delete fJetFinder;
fJetFinder = NULL;
if ( fJetHeader )
delete fJetHeader;
fJetHeader = NULL;
- */
- /*
+
if ( fJetReader )
delete fJetReader;
fJetReader = NULL;
-
+
if ( fJetReaderHeader )
delete fJetReaderHeader;
fJetReaderHeader = NULL;
- */
+
return 0;
}
// ADD MC Object
+ // -- Set Input Event
+ fJetFinder->GetReader()->SetInputEvent( NULL, NULL, const_cast<TObject*>(iter) );
+
AliMCEvent* foo = ( AliMCEvent* ) iter;
cout << foo->GetNumberOfTracks() << " -- "
- << foo->Stack() << " -- "
- << foo->Header()->Stack() << " -- "
- << foo->Header() << endl;
-
- printf (" --- MC --- \n");
+ << foo->Stack() << " -- "
+ << foo->Header()->Stack() << " -- "
+ << foo->Header() << endl;
+
+ AliStack* stack = foo->Stack();
+ cout << "N tracks" << stack->GetNtrack() << endl;
+
+ for (Int_t iterStack = 0; iterStack < stack->GetNtrack(); iterStack++) {
+ cout << iterStack << " -- " << stack->Particle(iterStack) << endl;
+
+ }
+ // -- Process Event
+ // fJetFinder->ProcessEvent();
+ }
+
+
+ if ( GetFirstInputObject(kAliHLTDataTypeEOR) ) {
+ //fJetFinder->FinishRun();
}
+
+
for ( iter=GetFirstInputObject(kAliHLTDataTypeESDObject|kAliHLTDataOriginOffline); iter != NULL; iter=GetNextInputObject() ) {
// ADD ESD Object -- Offline
// ** PushBack ** \\
- // ----------- ** \\
- // ** PushBack ** \\
+
return 0;
}
// this is just a tool to switch the logging systems
AliHLTLogging log;
- log.SwitchAliLog(1);
+ log.SwitchAliLog(0);
AliHLTSystem gHLT;
gHLT.SetGlobalLoggingLevel(0x7F);
gHLT.LoadComponentLibraries("libAliHLTUtil.so");
gHLT.LoadComponentLibraries("libAliHLTJET.so");
+ gHLT.LoadComponentLibraries("libESD.so");
+ gHLT.LoadComponentLibraries("libSTEER.so");
+ gHLT.LoadComponentLibraries("libSTEERBase.so");
+ gHLT.LoadComponentLibraries("libAOD.so");
+ gHLT.LoadComponentLibraries("libANALYSIS.so");
+ gHLT.LoadComponentLibraries("libANALYSISalice.so");
+ gHLT.LoadComponentLibraries("libJETAN.so");
+ gHLT.LoadComponentLibraries("libJETANMC.so");
+
+ // ----------------------------//
// - - //
// -- Publisher Components -- //
// - - //
+ // ----------------------------//
arg.Form("-entrytype MCFAST -dataspec 0x0000001F -datapath /home/jthaeder/jet/data/v4-16-Rev-01/FastGen/kPythia6Jets104_125_14TeV/JET-ETA=-0.2,0.2_JET-ET=50,1000_R=0.4_500ev");
if (!writerInput.IsNull()) writerInput+=" ";
writerInput+="ESDMCEventPublisher";
-
+
+ // ----------------------------//
// - - //
// -- Processing Components -- //
// - - //
- AliHLTConfiguration jetFinder("JETConeJet", "JETConeJetFinder",
+ // ----------------------------//
+
+ // -- ConeJetFinder
+ /*
+ AliHLTConfiguration jetFinder("JETConeJet", "JETConeJetFinder",
+ "ESDMCEventPublisher","");
+
+ if (!writerInput.IsNull()) writerInput+=" ";
+ writerInput+="JETConeJet";
+ */
+
+ // -- FasTJetFinder
+ AliHLTConfiguration jetFinder("JETFastJet", "JETFastJetFinder",
"ESDMCEventPublisher","");
if (!writerInput.IsNull()) writerInput+=" ";
- writerInput+="JETConeJet";
+ writerInput+="JETFastJet";
+ // ----------------------------//
// - - //
// -- Sink Components -- //
// - - //
+ // ----------------------------//
// -- The RootFileWriter
AliHLTConfiguration rootWriter("RootWriter", "ROOTFileWriter",