Difference: GMMsciProjectProgress (1 vs. 3)

Simulations, Reconstruction and Analysis of HZ->bbar e-e+ Using Direct Reconstruction

My analysis looks at Z->e-e+ and H->b bbar using direct reconstruction. Attached files "MarlinProcessor.cc" and "MarlinProcessor.h".

NECESSARY CHANGES TO C++ CODE:

PANDORA and WOLF output different MC particles collections so have to use which is required in the analysis code, PANDORA and WOLF respectively:

if(*name=="MCParticle"){

//if(*name=="MCParticlesSkimmed"){

PANDORA and WOLF call different functions in the analysis code:

PANDORA requires only these two:

if(_HbbMC.size()==2)this->AnalyseFourJetEvent(); //Analyse only when MC b bar from the Higgs was found
this->AnalyseEmEp(); //Z reconstruction and e- e+ energy resolutions

WOLF requires only these two:

// if(_RecoParvec.size()!=0)this->AnalyseRecoParticlesEvent();
// if(_HbbMC.size()==2 && _fourjetvec.size()==4)this->AnalyseWOLFFourJetEvent(); //Analyse only when MC b bar from the Higgs was found

• MyTrackCheater or MyLEPTrackingProcessor
• MyClusterCheater or MyTrackwiseClustering

[ DEBUG1 "Marlin"] hello from ECALGarlicTrainer destructor 0x2b94baf9abe0
[ DEBUG1 "Marlin"] hello from ECALGarlic destructor 0x2b94baf9afa0
**********************************************
A runtime error occured - (uncaught exception):
lcio::Exception: ProcessorMgr::registerProcessor: unknown processor with type TrackwiseClustering !

Marlin will have to be terminated, sorry.
**********************************************

Reconstruction

MarlinReco:

Required MarlinReco processors are outlined in MarlinReco Manual

Setting up Marlin:

My file directory used for reconstruction is ~/public_ppe/ExampleDSTAnalysis where I initialise ILCSoft with initialiseILCSoft.sh, build the Marlin processor using build_the_processor.sh and finally run Marlin_first_run.sh to set up the environment corretly.

Running Marlin:

Marlin steer_ild.xml

steer_ild.xml uses e-gun.slcio and corresponding GearOutput.xml which have been generated in Mokka and copied from ~/public_ppe/Mokka_Test/ParticleGun/

e-gun.slcio:

/generator/generator particleGun
/gun/particle e-
/gun/energy 100 GeV
/gun/momentumSmearing 5 GeV
/gun/position 0 0 0
/gun/direction 1 0 1
/run/beamOn 3000
exit

Looking At SLCIO Files:

Marlin lookAtFilenew.xml

To view what the SLCIO file contains.

Current Status

MarlinReco:

Required processors for reconstruction of e- particle gun:

• MyTPCDigiProcessor
• MyVTXDigiProcessor
• MySimpleCaloDigi
• MyTrackCheater
• MyLEPTrackingProcessor
• MyClusterCheater
• MyTrackwiseClustering
• MyWolf

Errors:

Unknown processors:

• MyTrackCheater
• MyClusterCheater
• MyTrackwiseClustering
• MyWolf

These are required for the Particle Flow Algorithm processor of type "Wolf" (also 'unknown') which creates the Reconstructed Particle Collection called 'RecoParticles' which will then be used to look at the physics of the events.

-- GavinMurdoch - 2013-11-28