come leggere il preshower

i link:

http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/RecoEcal/EgammaClusterProducers/interface/PreshowerAnalyzer.h?revision=1.3&view=markup
http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/RecoEcal/EgammaClusterProducers/src/PreshowerAnalyzer.cc?revision=1.4&view=markup

(forse) con questa riga nel BuildFile:



con gli include:

#include "DataFormats/EgammaReco/interface/PreshowerCluster.h"
#include "RecoEcal/EgammaClusterProducers/interface/PreshowerClusterProducer.h"

con i parametri:

preshClusterCollectionX_ = ps.getParameter("preshClusterCollectionX");
preshClusterCollectionY_ = ps.getParameter("preshClusterCollectionY");
preshClusterProducer_ = ps.getParameter("preshClusterProducer");

con li codice:

// Get ES clusters in X plane
Handle pPreshowerClustersX;
evt.getByLabel(preshClusterProducer_, preshClusterCollectionX_, pPreshowerClustersX);
const reco::PreshowerClusterCollection *clustersX = pPreshowerClustersX.product();
h1_nclu_x->Fill( clustersX->size() );
//std::cout << "\n pPreshowerClustersX->size() = " <<>size() <<> pPreshowerClustersY;
evt.getByLabel(preshClusterProducer_, preshClusterCollectionY_, pPreshowerClustersY);
const reco::PreshowerClusterCollection *clustersY = pPreshowerClustersY.product();
h1_nclu_y->Fill( clustersY->size() );
//std::cout << "\n pPreshowerClustersY->size() = " <<>size() << e1 =" 0;" esclus =" clustersX-">begin();
esClus !=clustersX->end(); esClus++) {
e1 += esClus->energy();
h1_esE_x->Fill( esClus->energy() );
h1_esEta_x->Fill( esClus->eta() );
h1_esPhi_x->Fill( esClus->phi() );
h1_esNhits_x->Fill( esClus->nhits() );
}

float e2 = 0;
for(reco::PreshowerClusterCollection::const_iterator esClus = clustersY->begin();
esClus !=clustersY->end(); esClus++) {
e2 += esClus->energy();
h1_esE_y->Fill( esClus->energy() );
h1_esEta_y->Fill( esClus->eta() );
h1_esPhi_y->Fill( esClus->phi() );
h1_esNhits_y->Fill( esClus->nhits() );
}

float deltaE = 0;
if(e1+e2 > 1.0e-10) {
// GeV to #MIPs
e1 = e1 / mip_;
e2 = e2 / mip_;
deltaE = gamma_*(calib_planeX_*e1+calib_planeY_*e2);
}

h1_esDeltaE->Fill(deltaE);

OPPURE, PIU' SEMPLICEMENTE:

SuperCluster::preshowerEnergy ()
http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/CMSSW/DataFormats/EgammaReco/interface/SuperCluster.h?revision=1.10&view=markup