ヒストグラムを作る

#include <TH1.h>

TH1 *h_electronPtOriginal; //!

TH1 *h_electronPtCorrected; //!

h_electronPtOriginal = new TH1F("h_electronPtOriginal", "h_electronPtOriginal", 100, 0, 100);

wk()->addOutput(h_electronPtOriginal);

h_electronPtCorrected = new TH1F("h_electronPtCorrected", "h_electronPtCorrected", 100, 0, 100);

wk()->addOutput(h_electronPtCorrected);

// Electrons

const xAOD::ElectronContainer* electrons = nullptr;

ANA_CHECK(event->retrieve( electrons, "Electrons"));

xAOD::ElectronContainer::const_iterator electron_itr = electrons->begin();

xAOD::ElectronContainer::const_iterator electron_end = electrons->end();

for ( ; electron_itr != electron_end; ++electron_itr ){

bool LHMediumSel = m_MediumLH->accept((*electron_itr));

if ( !LHMediumSel ) continue;

if( fabs((*electron_itr)->eta()) > 2.47 ) continue;

h_electronPtOriginal->Fill( (*electron_itr)->pt()*0.001 );

Info("execute()", " electron pt = %.2f GeV", ((*electron_itr)->pt()*0.001));

}

std::pair< xAOD::ElectronContainer*, xAOD::ShallowAuxContainer* > electrons_shallowCopy = xAOD::shallowCopyContainer( *electrons );

xAOD::ElectronContainer::iterator electronSC_itr = (electrons_shallowCopy.first)->begin();

xAOD::ElectronContainer::iterator electronSC_end = (electrons_shallowCopy.first)->end();

for( ; electronSC_itr != electronSC_end; ++electronSC_itr ) {

if(m_electronCalibrationAndSmearingTool->applyCorrection(**electronSC_itr) == CP::CorrectionCode::Error){

Error("execute()", "ElectronCalibrationAndSmearingTool returns Error CorrectionCode");

}

bool LHMediumSel = m_MediumLH->accept((*electronSC_itr));

if ( !LHMediumSel ) continue;

if ( fabs((*electronSC_itr)->eta()) > 2.47 ) continue;

h_electronPtCorrected->Fill( (*electronSC_itr)->pt()*0.001 );

Info("execute()", " corrected electron pt = %.2f GeV", ((*electronSC_itr)->pt() * 0.001));

}

delete electrons_shallowCopy.first;

delete electrons_shallowCopy.second;