Toot toot, time to use root

Simon At Home: yay - it's created simon.root
Jim: now you need to fire up Root and read the file
Jim: find out what's in it, and try plotting some stuff
Simon At Home: :shudder:
Jim: for example, generator jet et distribution
Simon At Home: ...time to learn root :)
Jim: yep, and to work out the meaning of all the jetmet quantities ;-)

Now to load my root file and plot some graphs.....

Trying again

With advice from Jim:
"- create fresh orca area and check out JetMetAnalysis
- build in domain/test
- extract .orcarc contents from run_db.csh
- run JetMet"

So deleted the Orca directory, redid the check out etc.
- cd JetMetAnalysis/domain/test
- eval `scram runtime -csh`
- scram build JetMet
- rehash

Set up .orcarc (link)
- JetMet

Result: It looks to have worked! It's created simon.root at least.... (log file)

Nope it failed

It really looked like it was doing something but get this:

---

------ Information on L1 jets :
*** isolated em objects : 0
*** central jets : 1
*** forward jets : 3
*** tau objects : 1
*** global objects : 1
------ MEMORY use statistics -------
CobraSignal: SEGV received.
in thread 14331 2051
0x42ad0b8d BackTrace::trace(std::ostream&) const + 0x2d [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libGenUtil.so]
0x42a34156 CobraSignalHandler::TermSignalHandler(int) + 0x246 [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libBaseMonitor.so]
0x43bb0f05 [/lib/i686/libpthread.so.0]
0x43be5188 [/lib/i686/libc.so.6]
0x43c36d24 __libc_free + 0xb4 [/lib/i686/libc.so.6]
0x43b58953 operator delete(void*) + 0x23 [/usr/local/gcc-alt-3.2.3/lib/libstdc++.so.5]
0x4204cdca ConcreteJet::~ConcreteJet() + 0x66 [/afs/cern.ch/cms/Releases/ORCA/ORCA_8_3_0/lib/Linux__2.4/libJetTypes.so]
0x42356fee TReconstructor::~TReconstructor() + 0x4e [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libReco.so]
0x4235542e TRecEvent::~TRecEvent() + 0x8e [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libReco.so]
0x422f8326 TRecEventWP::~TRecEventWP() + 0xa6 [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libEvent.so]
0x4226ffc8 G3EventProxy::reset() + 0x418 [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libG3Event.so]
0x4226f85e G3EventProxy::~G3EventProxy() + 0x3e [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libG3Event.so]
0x4219f64a SimpleRecEventProcessor::clearEvent() + 0xaa [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libRecReader.so]
0x42197698 TEventReader::nextEvent(pool::Ref const&) + 0x48 [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libRecReader.so]
0x421938a5 RecEventProcessor::consume(std::pair > const + 0xa5 [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libRecReader.so]
0x42171eaf CollectionProcessingThread::run() + 0xff [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libBatchRecReader.so]
0x42afd9e3 Thread::operator()() + 0x313 [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libGenUtil.so]
0x42afe3b1 boost::detail::function::void_function_obj_invoker0::invoke(boost::detail::functio + 0x21 [/afs/cern.ch/cms/Releases/COBRA/COBRA_7_8_6/lib/Linux__2.4/libGenUtil.so]
0x43a41b8a [/afs/cern.ch/cms/external/lcg/external/Boost/1.30.2/rh73_gcc32/lib/libboost_thread.so]
0x43badfaf [/lib/i686/libpthread.so.0]
0x43ca392a __clone + 0x3a [/lib/i686/libc.so.6]

Forced KILL,

Killed

---

Balls!

Google is your friend

It seems that all the orca documentation is spread out over the net and not really connected or anything, so google is your friend (cheers Jim for the pointer).
So I was getting: "CobraSignal: SEGV received." googled it, and it turns out its caused by not setting the OutputDataSet in .orcarc. Just set that and it's now running - taking ages mind! - but it is doing stuff.....

Couple of useful links

USCMS JetMet group
The Algorithms

Starting from scratch

I want to do some JetMet analysis. To this end I'm going to check out from CVS the JetMet code and try and run their example jobs on lxplus.
Lets see how it goes (sort of following this):

1. scram project ORCA ORCA_8_3_0
   
2. cd ORCA_8_3_0 /src
   
3. eval `scram runtime -csh`
   
4. cmscvsroot ORCA
   
5. cvs login
   
6. cvs co -r ORCA_8_3_0 JetMetAnalysis

That sets up the workspace and gets all the relevant files....

Now to try and build some executables, and maybe even a root tree!

1. cd JetMetAnalysis/JetMetTest/test
   
2. scram b shared > JetMetTest_scram_b_shared.log & (to store it for posterity)
   
3. scram b bin > JetMetTest_scram_b_bin.log &
   
4. eval `scram runtime -csh`
   
5. rehash

Now to run an exe:

1. JetMetTest (to run the exe - this will probably fail as I have no .orcarc set up)
   
2. Failed with a seg fault - no data to run on time to set the orcarc
   
3. vi .orcarc and enter the following:
   
   PoolCatalogFile = @{xmlcatalog_http://cmsdoc.cern.ch/orca/catalog/PoolFileCatalog_8_3_0.xml}@
   InputCollections=/System/StW830DST2x1033/jets50100/jets50100
   MaxEvents=10
   
4. JetMetTest > JetMetTest.log &
   
5. Takes a while....
   
6. ...but it is doing something...
   
7. ...and then it fails:

Try with a different exe JetMetSmallTest:

1. JetMetSmallTest > JetMetSmallTest.log
2. runs and doesn't crash, but produces no output

First Post

Well I'm struggling to learn ORCA which is bordering on a total nightmare! So going to use this blog as an online work book and keep track of the steps I do, you never know it might come in handy for thesis or for some other poor soul learning the CMS analysis framework.