|
|
|
#include <ShallowClustersProducer.h>
Classes | |
| struct | moduleVars |
| struct | NearDigis |
Public Member Functions | |
| ShallowClustersProducer (const edm::ParameterSet &) | |
Private Member Functions | |
| void | produce (edm::Event &, const edm::EventSetup &) |
Private Attributes | |
| std::string | Prefix |
| edm::InputTag | theClustersLabel |
Definition at line 11 of file ShallowClustersProducer.h.
| ShallowClustersProducer::ShallowClustersProducer | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 15 of file ShallowClustersProducer.cc.
References Prefix.
: theClustersLabel(iConfig.getParameter<edm::InputTag>("Clusters")), Prefix(iConfig.getParameter<std::string>("Prefix") ) { produces <std::vector<unsigned> > ( Prefix + "number" ); produces <std::vector<unsigned> > ( Prefix + "width" ); produces <std::vector<float> > ( Prefix + "variance" ); produces <std::vector<float> > ( Prefix + "barystrip" ); produces <std::vector<float> > ( Prefix + "middlestrip" ); produces <std::vector<unsigned> > ( Prefix + "charge" ); produces <std::vector<float> > ( Prefix + "noise" ); produces <std::vector<float> > ( Prefix + "ston" ); produces <std::vector<unsigned> > ( Prefix + "seedstrip" ); produces <std::vector<unsigned> > ( Prefix + "seedindex" ); produces <std::vector<unsigned> > ( Prefix + "seedcharge" ); produces <std::vector<float> > ( Prefix + "seednoise" ); produces <std::vector<float> > ( Prefix + "seedgain" ); produces <std::vector<unsigned> > ( Prefix + "qualityisbad" ); produces <std::vector<float> > ( Prefix + "rawchargeC" ); produces <std::vector<float> > ( Prefix + "rawchargeL" ); produces <std::vector<float> > ( Prefix + "rawchargeR" ); produces <std::vector<float> > ( Prefix + "rawchargeLL" ); produces <std::vector<float> > ( Prefix + "rawchargeRR" ); produces <std::vector<float> > ( Prefix + "eta" ); produces <std::vector<float> > ( Prefix + "foldedeta" ); produces <std::vector<float> > ( Prefix + "etaX" ); produces <std::vector<float> > ( Prefix + "etaasymm" ); produces <std::vector<float> > ( Prefix + "outsideasymm"); produces <std::vector<float> > ( Prefix + "neweta"); produces <std::vector<float> > ( Prefix + "newetaerr"); produces <std::vector<unsigned> > ( Prefix + "detid" ); produces <std::vector<int> > ( Prefix + "subdetid" ); produces <std::vector<int> > ( Prefix + "module" ); produces <std::vector<int> > ( Prefix + "side" ); produces <std::vector<int> > ( Prefix + "layerwheel" ); produces <std::vector<int> > ( Prefix + "stringringrod" ); produces <std::vector<int> > ( Prefix + "petal" ); produces <std::vector<int> > ( Prefix + "stereo" ); }
| void ShallowClustersProducer::produce | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private, virtual] |
Implements edm::EDProducer.
Definition at line 59 of file ShallowClustersProducer.cc.
References edmNew::DetSetVector< T >::begin(), SiStripClusterInfo::charge(), DeDxDiscriminatorTools::charge(), cond::rpcobgas::detid, edmNew::DetSetVector< T >::end(), eta(), ShallowClustersProducer::NearDigis::eta(), ShallowClustersProducer::NearDigis::etaasymm(), ShallowClustersProducer::NearDigis::etaX(), ShallowClustersProducer::NearDigis::first, SiStripClusterInfo::firstStrip(), edm::Event::getByLabel(), edmNew::DetSetVector< T >::id(), info, SiStripClusterInfo::IsAnythingBad(), edm::HandleBase::isValid(), ShallowClustersProducer::NearDigis::last, ShallowClustersProducer::moduleVars::layerwheel, ShallowClustersProducer::NearDigis::left, ShallowClustersProducer::NearDigis::Lleft, ShallowClustersProducer::NearDigis::max, SiStripClusterInfo::maxCharge(), SiStripClusterInfo::maxIndex(), SiStripClusterInfo::maxStrip(), python::rootplot::argparse::module, ShallowClustersProducer::moduleVars::module, SiStripClusterInfo::noiseRescaledByGain(), ShallowClustersProducer::NearDigis::outsideasymm(), ShallowClustersProducer::moduleVars::petal, funct::pow(), Prefix, edm::Event::put(), ShallowClustersProducer::NearDigis::right, ShallowClustersProducer::NearDigis::Rright, ShallowClustersProducer::moduleVars::side, SiStripClusterInfo::signalOverNoise(), mathSSE::sqrt(), ShallowClustersProducer::moduleVars::stereo, ShallowClustersProducer::moduleVars::stringringrod, SiStripClusterInfo::stripGains(), SiStripClusterInfo::stripNoisesRescaledByGain(), ShallowClustersProducer::moduleVars::subdetid, theClustersLabel, SiStripClusterInfo::variance(), tablePrinter::width, and SiStripClusterInfo::width().
{
std::auto_ptr<std::vector<unsigned> > number ( new std::vector<unsigned>(7,0) );
std::auto_ptr<std::vector<unsigned> > width ( new std::vector<unsigned>() );
std::auto_ptr<std::vector<float> > variance ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > barystrip ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > middlestrip ( new std::vector<float>() );
std::auto_ptr<std::vector<unsigned> > charge ( new std::vector<unsigned>() );
std::auto_ptr<std::vector<float> > noise ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > ston ( new std::vector<float>() );
std::auto_ptr<std::vector<unsigned> > seedstrip ( new std::vector<unsigned>() );
std::auto_ptr<std::vector<unsigned> > seedindex ( new std::vector<unsigned>() );
std::auto_ptr<std::vector<unsigned> > seedcharge ( new std::vector<unsigned>() );
std::auto_ptr<std::vector<float> > seednoise ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > seedgain ( new std::vector<float>() );
std::auto_ptr<std::vector<unsigned> > qualityisbad ( new std::vector<unsigned>() );
std::auto_ptr<std::vector<float> > rawchargeC ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > rawchargeL ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > rawchargeR ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > rawchargeLL ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > rawchargeRR ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > etaX ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > eta ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > foldedeta ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > etaasymm ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > outsideasymm ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > neweta ( new std::vector<float>() );
std::auto_ptr<std::vector<float> > newetaerr ( new std::vector<float>() );
std::auto_ptr<std::vector<unsigned> > detid ( new std::vector<unsigned>() );
std::auto_ptr<std::vector<int> > subdetid ( new std::vector<int>() );
std::auto_ptr<std::vector<int> > side ( new std::vector<int>() );
std::auto_ptr<std::vector<int> > module ( new std::vector<int>() );
std::auto_ptr<std::vector<int> > layerwheel ( new std::vector<int>() );
std::auto_ptr<std::vector<int> > stringringrod ( new std::vector<int>() );
std::auto_ptr<std::vector<int> > petal ( new std::vector<int>() );
std::auto_ptr<std::vector<int> > stereo ( new std::vector<int>());
edm::Handle<edmNew::DetSetVector<SiStripCluster> > clusters;
iEvent.getByLabel(theClustersLabel, clusters);
edm::Handle<edm::DetSetVector<SiStripProcessedRawDigi> > rawProcessedDigis;
iEvent.getByLabel("siStripProcessedRawDigis", "", rawProcessedDigis);
edmNew::DetSetVector<SiStripCluster>::const_iterator itClusters=clusters->begin();
for(;itClusters!=clusters->end();++itClusters){
uint32_t id = itClusters->id();
const moduleVars moduleV(id);
for(edmNew::DetSet<SiStripCluster>::const_iterator cluster=itClusters->begin(); cluster!=itClusters->end();++cluster){
const SiStripClusterInfo info(*cluster, iSetup, id);
const NearDigis digis = rawProcessedDigis.isValid() ? NearDigis(info, *rawProcessedDigis) : NearDigis(info);
(number->at(0))++;
(number->at(moduleV.subdetid))++;
width->push_back( cluster->amplitudes().size() );
barystrip->push_back( cluster->barycenter() );
variance->push_back( info.variance() );
middlestrip->push_back( info.firstStrip() + info.width()/2.0 );
charge->push_back( info.charge() );
noise->push_back( info.noiseRescaledByGain() );
ston->push_back( info.signalOverNoise() );
seedstrip->push_back( info.maxStrip() );
seedindex->push_back( info.maxIndex() );
seedcharge->push_back( info.maxCharge() );
seednoise->push_back( info.stripNoisesRescaledByGain().at(info.maxIndex()) );
seedgain->push_back( info.stripGains().at(info.maxIndex()) );
qualityisbad->push_back( info.IsAnythingBad() );
rawchargeC->push_back( digis.max );
rawchargeL->push_back( digis.left );
rawchargeR->push_back( digis.right );
rawchargeLL->push_back( digis.Lleft );
rawchargeRR->push_back( digis.Rright );
etaX->push_back( digis.etaX() );
eta->push_back( digis.eta() );
etaasymm->push_back( digis.etaasymm() );
outsideasymm->push_back( digis.outsideasymm() );
neweta->push_back( (digis.last-digis.first)/info.charge() );
newetaerr->push_back( (sqrt(digis.last+digis.first))/pow(info.charge(),1.5) );
detid->push_back( id );
subdetid->push_back( moduleV.subdetid );
side->push_back( moduleV.side );
module->push_back( moduleV.module );
layerwheel->push_back( moduleV.layerwheel );
stringringrod->push_back( moduleV.stringringrod );
petal->push_back( moduleV.petal );
stereo->push_back( moduleV.stereo );
}
}
iEvent.put( number, Prefix + "number" );
iEvent.put( width, Prefix + "width" );
iEvent.put( variance, Prefix + "variance" );
iEvent.put( barystrip, Prefix + "barystrip" );
iEvent.put( middlestrip, Prefix + "middlestrip" );
iEvent.put( charge, Prefix + "charge" );
iEvent.put( noise, Prefix + "noise" );
iEvent.put( ston, Prefix + "ston" );
iEvent.put( seedstrip, Prefix + "seedstrip" );
iEvent.put( seedindex, Prefix + "seedindex" );
iEvent.put( seedcharge, Prefix + "seedcharge" );
iEvent.put( seednoise, Prefix + "seednoise" );
iEvent.put( seedgain, Prefix + "seedgain" );
iEvent.put( qualityisbad, Prefix + "qualityisbad" );
iEvent.put( rawchargeC, Prefix + "rawchargeC" );
iEvent.put( rawchargeL, Prefix + "rawchargeL" );
iEvent.put( rawchargeR, Prefix + "rawchargeR" );
iEvent.put( rawchargeLL, Prefix + "rawchargeLL" );
iEvent.put( rawchargeRR, Prefix + "rawchargeRR" );
iEvent.put( etaX, Prefix + "etaX" );
iEvent.put( eta, Prefix + "eta" );
iEvent.put( foldedeta, Prefix + "foldedeta" );
iEvent.put( etaasymm, Prefix + "etaasymm" );
iEvent.put( outsideasymm, Prefix + "outsideasymm" );
iEvent.put( neweta, Prefix + "neweta" );
iEvent.put( newetaerr, Prefix + "newetaerr" );
iEvent.put( detid, Prefix + "detid" );
iEvent.put( subdetid, Prefix + "subdetid" );
iEvent.put( module, Prefix + "module" );
iEvent.put( side, Prefix + "side" );
iEvent.put( layerwheel, Prefix + "layerwheel" );
iEvent.put( stringringrod, Prefix + "stringringrod" );
iEvent.put( petal, Prefix + "petal" );
iEvent.put( stereo, Prefix + "stereo" );
}
std::string ShallowClustersProducer::Prefix [private] |
Definition at line 20 of file ShallowClustersProducer.h.
Referenced by produce(), and ShallowClustersProducer().
Definition at line 19 of file ShallowClustersProducer.h.
Referenced by produce().
1.7.3