Commit c1d7e289 authored by Mark Hodgkinson's avatar Mark Hodgkinson
Browse files

Add calib hit moments to neutral FlowElements.

Update job options to enable PFNeutralFlowElementCreatorAlgorithm.cxx
Update StreamESD/AOD to unclude all FlowElements from particle flow.
parent 314a765a
......@@ -250,8 +250,12 @@ if jobproperties.eflowRecFlags.usePFEGammaPFOAssoc:
PFEGammaPFOAssoc=PFEGammaPFOAssoc("PFEGammaPFOAssoc")
topSequence += PFEGammaPFOAssoc
#Add new FlowElement creator
#Add new FlowElement creators
from eflowRec.eflowRecConf import PFChargedFlowElementCreatorAlgorithm
PFChargedFlowElementCreatorAlgorithm = PFChargedFlowElementCreatorAlgorithm("PFChargedFlowElementCreatorAlgorithm")
topSequence += PFChargedFlowElementCreatorAlgorithm
from eflowRec.eflowRecConf import PFNeutralFlowElementCreatorAlgorithm
PFNeutralFlowElementCreatorAlgorithm = PFNeutralFlowElementCreatorAlgorithm("PFNeutralFlowElementCreatorAlgorithm")
topSequence += PFNeutralFlowElementCreatorAlgorithm
......@@ -23,4 +23,7 @@ if True == jobproperties.eflowRecFlags.provideShowerSubtractedClusters:
#Add new flow elements
eflowRecAODList += [ "xAOD::FlowElementContainer#JetETMissChargedFlowElements"]
eflowRecAODList += [ "xAOD::FlowElementAuxContainer#JetETMissChargedFlowElementsAux."]
\ No newline at end of file
eflowRecAODList += [ "xAOD::FlowElementAuxContainer#JetETMissChargedFlowElementsAux."]
eflowRecAODList += [ "xAOD::FlowElementContainer#JetETMissNeutralFlowElements"]
eflowRecAODList += [ "xAOD::FlowElementAuxContainer#JetETMissNeutralFlowElementsAux."]
\ No newline at end of file
......@@ -21,3 +21,9 @@ if True == jobproperties.eflowRecFlags.provideShowerSubtractedClusters:
eflowRecESDList += [ "xAOD::CaloClusterContainer#PFCaloCluster"]
eflowRecESDList += [ "xAOD::CaloClusterAuxContainer#PFCaloClusterAux."]
#Add new flow elements
eflowRecAODList += [ "xAOD::FlowElementContainer#JetETMissChargedFlowElements"]
eflowRecAODList += [ "xAOD::FlowElementAuxContainer#JetETMissChargedFlowElementsAux."]
eflowRecAODList += [ "xAOD::FlowElementContainer#JetETMissNeutralFlowElements"]
eflowRecAODList += [ "xAOD::FlowElementAuxContainer#JetETMissNeutralFlowElementsAux."]
\ No newline at end of file
......@@ -120,6 +120,13 @@ StatusCode PFNeutralFlowElementCreatorAlgorithm::createNeutralFlowElement(const
this->addMoment(xAOD::CaloCluster::ISOLATION,"eflowRec_ISOLATION",*cluster,*thisFE);
this->addMoment(xAOD::CaloCluster::SECOND_LAMBDA,"eflowRec_SECOND_LAMBDA",*cluster,*thisFE);
this->addMoment(xAOD::CaloCluster::EM_PROBABILITY,"eflowRec_EM_PROBABILITY",*cluster,*thisFE);
if (m_useCalibHitTruth){
this->addMoment(xAOD::CaloCluster::ENG_CALIB_TOT,"eflowRec_ENG_CALIB_TOT",*cluster,*thisFE);
this->addMoment(xAOD::CaloCluster::ENG_CALIB_FRAC_EM,"eflowRec_ENG_CALIB_FRAC_EM",*cluster,*thisFE);
this->addMoment(xAOD::CaloCluster::ENG_CALIB_FRAC_HAD,"eflowRec_ENG_CALIB_FRAC_HAD",*cluster,*thisFE);
this->addMoment(xAOD::CaloCluster::ENG_CALIB_FRAC_REST,"eflowRec_ENG_CALIB_FRAC_REST",*cluster,*thisFE);
}
}//cluster loop
return StatusCode::SUCCESS;
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment