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

All using ILCSoft version 01-14, with Marlin Processor build details outlined here.

Final report for project is attached as "Murdoch_Gavin_0900886_ProjectReport_40Credit.pdf".

Simulations:

Firstly, simulations were ran for stdhep files obtained from .... Process and working script available in automating job submission section of RunMokkaLocally.

Reconstruction:

Reconstruction utilises the following particle flow algorithms (PFAs): PANDORA and WOLF.

Some MarlinReco processors are outlined in MarlinReco Manual however, the attached files are working with the most current processors!! More info on UsingMarlinReco.

PANDORA:

Attached file "CLIC_ILD_CDR500_STEERING_POST_PANDORA_AND_TAGGING_def_10000.xml".

Input all simulated events and be sure to use one of the corresponding gear files (N.B. all will have the same properties). This xml steering file will use the PANDORA PFA and create an output slcio file ready for analysis with MyMarlinProcessor c++ code outlined later.

WOLF:

Attached file "steer_ldc_WOLF_AND_TAGGING_def_10000.xml".

Same procedure as for PANDORA.

ANALYSIS:

XML STEERING FILE

Using the xml steering files for either PANDORA or WOLF, comment out all processors and uncomment the "MyMarlinProcessor"(specific naming for my analysis processor and c++ code). Be sure to change the input slcio file to the output slcio created when utilising the PFAs.

C++ ANALYSIS CODE

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