|
|
|
#include <GlobalRecHitsProducer.h>
Definition at line 158 of file GlobalRecHitsProducer.h.
| typedef std::vector<double> GlobalRecHitsProducer::DoubleVector |
Definition at line 164 of file GlobalRecHitsProducer.h.
| typedef std::vector<float> GlobalRecHitsProducer::FloatVector |
Definition at line 163 of file GlobalRecHitsProducer.h.
| typedef std::vector<int> GlobalRecHitsProducer::IntVector |
Definition at line 165 of file GlobalRecHitsProducer.h.
| typedef std::map<uint32_t,float,std::less<uint32_t> > GlobalRecHitsProducer::MapType |
Definition at line 166 of file GlobalRecHitsProducer.h.
| GlobalRecHitsProducer::GlobalRecHitsProducer | ( | const edm::ParameterSet & | iPSet | ) | [explicit] |
Definition at line 13 of file GlobalRecHitsProducer.cc.
References conf_, ECalEBSrc_, ECalEESrc_, ECalESSrc_, ECalUncalEBSrc_, ECalUncalEESrc_, fName, frequency, getAllProvenances, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HCalSrc_, edm::InputTag::instance(), edm::InputTag::label(), label, MuCSCSrc_, MuDTSimSrc_, MuDTSrc_, MuRPCSimSrc_, MuRPCSrc_, printProvenanceInfo, SiPxlSrc_, SiStripSrc_, and verbosity.
: fName(""), verbosity(0), frequency(0), label(""), getAllProvenances(false), printProvenanceInfo(false), count(0) { std::string MsgLoggerCat = "GlobalRecHitsProducer_GlobalRecHitsProducer"; // get information from parameter set fName = iPSet.getUntrackedParameter<std::string>("Name"); verbosity = iPSet.getUntrackedParameter<int>("Verbosity"); frequency = iPSet.getUntrackedParameter<int>("Frequency"); label = iPSet.getParameter<std::string>("Label"); edm::ParameterSet m_Prov = iPSet.getParameter<edm::ParameterSet>("ProvenanceLookup"); getAllProvenances = m_Prov.getUntrackedParameter<bool>("GetAllProvenances"); printProvenanceInfo = m_Prov.getUntrackedParameter<bool>("PrintProvenanceInfo"); //get Labels to use to extract information ECalEBSrc_ = iPSet.getParameter<edm::InputTag>("ECalEBSrc"); ECalUncalEBSrc_ = iPSet.getParameter<edm::InputTag>("ECalUncalEBSrc"); ECalEESrc_ = iPSet.getParameter<edm::InputTag>("ECalEESrc"); ECalUncalEESrc_ = iPSet.getParameter<edm::InputTag>("ECalUncalEESrc"); ECalESSrc_ = iPSet.getParameter<edm::InputTag>("ECalESSrc"); HCalSrc_ = iPSet.getParameter<edm::InputTag>("HCalSrc"); SiStripSrc_ = iPSet.getParameter<edm::InputTag>("SiStripSrc"); SiPxlSrc_ = iPSet.getParameter<edm::InputTag>("SiPxlSrc"); MuDTSrc_ = iPSet.getParameter<edm::InputTag>("MuDTSrc"); MuDTSimSrc_ = iPSet.getParameter<edm::InputTag>("MuDTSimSrc"); MuCSCSrc_ = iPSet.getParameter<edm::InputTag>("MuCSCSrc"); MuRPCSrc_ = iPSet.getParameter<edm::InputTag>("MuRPCSrc"); MuRPCSimSrc_ = iPSet.getParameter<edm::InputTag>("MuRPCSimSrc"); conf_ = iPSet; // use value of first digit to determine default output level (inclusive) // 0 is none, 1 is basic, 2 is fill output, 3 is gather output verbosity %= 10; // create persistent object produces<PGlobalRecHit>(label); // print out Parameter Set information being used if (verbosity >= 0) { edm::LogInfo(MsgLoggerCat) << "\n===============================\n" << "Initialized as EDProducer with parameter values:\n" << " Name = " << fName << "\n" << " Verbosity = " << verbosity << "\n" << " Frequency = " << frequency << "\n" << " Label = " << label << "\n" << " GetProv = " << getAllProvenances << "\n" << " PrintProv = " << printProvenanceInfo << "\n" << " ECalEBSrc = " << ECalEBSrc_.label() << ":" << ECalEBSrc_.instance() << "\n" << " ECalUncalEBSrc = " << ECalUncalEBSrc_.label() << ":" << ECalUncalEBSrc_.instance() << "\n" << " ECalEESrc = " << ECalEESrc_.label() << ":" << ECalUncalEESrc_.instance() << "\n" << " ECalUncalEESrc = " << ECalUncalEESrc_.label() << ":" << ECalEESrc_.instance() << "\n" << " ECalESSrc = " << ECalESSrc_.label() << ":" << ECalESSrc_.instance() << "\n" << " HCalSrc = " << HCalSrc_.label() << ":" << HCalSrc_.instance() << "\n" << " SiStripSrc = " << SiStripSrc_.label() << ":" << SiStripSrc_.instance() << "\n" << " SiPixelSrc = " << SiPxlSrc_.label() << ":" << SiPxlSrc_.instance() << "\n" << " MuDTSrc = " << MuDTSrc_.label() << ":" << MuDTSrc_.instance() << "\n" << " MuDTSimSrc = " << MuDTSimSrc_.label() << ":" << MuDTSimSrc_.instance() << "\n" << " MuCSCSrc = " << MuCSCSrc_.label() << ":" << MuCSCSrc_.instance() << "\n" << " MuRPCSrc = " << MuRPCSrc_.label() << ":" << MuRPCSrc_.instance() << "\n" << " MuRPCSimSrc = " << MuRPCSimSrc_.label() << ":" << MuRPCSimSrc_.instance() << "\n" << "===============================\n"; } }
| GlobalRecHitsProducer::~GlobalRecHitsProducer | ( | ) | [virtual] |
Definition at line 96 of file GlobalRecHitsProducer.cc.
{
}
| void GlobalRecHitsProducer::beginJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDProducer.
Definition at line 100 of file GlobalRecHitsProducer.cc.
References clear().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_beginJob";
// clear storage vectors
clear();
return;
}
| void GlobalRecHitsProducer::clear | ( | void | ) | [private] |
Definition at line 1926 of file GlobalRecHitsProducer.cc.
References BRL1RX, BRL1RY, BRL1SX, BRL1SY, BRL2RX, BRL2RY, BRL2SX, BRL2SY, BRL3RX, BRL3RY, BRL3SX, BRL3SY, CSCRHPERP, CSCRHPHI, CSCSHPHI, DTRHD, DTSHD, EBRE, EBSHE, EERE, EESHE, ESRE, ESSHE, FWD1nRX, FWD1nRY, FWD1nSX, FWD1nSY, FWD1pRX, FWD1pRY, FWD1pSX, FWD1pSY, FWD2nRX, FWD2nRY, FWD2nSX, FWD2nSY, FWD2pRX, FWD2pRY, FWD2pSX, FWD2pSY, HBCalR, HBCalREC, HBCalSHE, HECalR, HECalREC, HECalSHE, HFCalR, HFCalREC, HFCalSHE, HOCalR, HOCalREC, HOCalSHE, RPCRHX, RPCSHX, TECW1RX, TECW1RY, TECW1SX, TECW1SY, TECW2RX, TECW2RY, TECW2SX, TECW2SY, TECW3RX, TECW3RY, TECW3SX, TECW3SY, TECW4RX, TECW4RY, TECW4SX, TECW4SY, TECW5RX, TECW5RY, TECW5SX, TECW5SY, TECW6RX, TECW6RY, TECW6SX, TECW6SY, TECW7RX, TECW7RY, TECW7SX, TECW7SY, TECW8RX, TECW8RY, TECW8SX, TECW8SY, TIBL1RX, TIBL1RY, TIBL1SX, TIBL1SY, TIBL2RX, TIBL2RY, TIBL2SX, TIBL2SY, TIBL3RX, TIBL3RY, TIBL3SX, TIBL3SY, TIBL4RX, TIBL4RY, TIBL4SX, TIBL4SY, TIDW1RX, TIDW1RY, TIDW1SX, TIDW1SY, TIDW2RX, TIDW2RY, TIDW2SX, TIDW2SY, TIDW3RX, TIDW3RY, TIDW3SX, TIDW3SY, TOBL1RX, TOBL1RY, TOBL1SX, TOBL1SY, TOBL2RX, TOBL2RY, TOBL2SX, TOBL2SY, TOBL3RX, TOBL3RY, TOBL3SX, TOBL3SY, TOBL4RX, TOBL4RY, TOBL4SX, TOBL4SY, and verbosity.
Referenced by beginJob(), and produce().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_clear";
if (verbosity > 0)
edm::LogInfo(MsgLoggerCat)
<< "Clearing event holders";
// reset electromagnetic info
// EE info
EERE.clear();
EESHE.clear();
// EB info
EBRE.clear();
EBSHE.clear();
// ES info
ESRE.clear();
ESSHE.clear();
// reset HCal Info
HBCalREC.clear();
HBCalR.clear();
HBCalSHE.clear();
HECalREC.clear();
HECalR.clear();
HECalSHE.clear();
HOCalREC.clear();
HOCalR.clear();
HOCalSHE.clear();
HFCalREC.clear();
HFCalR.clear();
HFCalSHE.clear();
// reset Track Info
TIBL1RX.clear();
TIBL2RX.clear();
TIBL3RX.clear();
TIBL4RX.clear();
TIBL1RY.clear();
TIBL2RY.clear();
TIBL3RY.clear();
TIBL4RY.clear();
TIBL1SX.clear();
TIBL2SX.clear();
TIBL3SX.clear();
TIBL4SX.clear();
TIBL1SY.clear();
TIBL2SY.clear();
TIBL3SY.clear();
TIBL4SY.clear();
TOBL1RX.clear();
TOBL2RX.clear();
TOBL3RX.clear();
TOBL4RX.clear();
TOBL1RY.clear();
TOBL2RY.clear();
TOBL3RY.clear();
TOBL4RY.clear();
TOBL1SX.clear();
TOBL2SX.clear();
TOBL3SX.clear();
TOBL4SX.clear();
TOBL1SY.clear();
TOBL2SY.clear();
TOBL3SY.clear();
TOBL4SY.clear();
TIDW1RX.clear();
TIDW2RX.clear();
TIDW3RX.clear();
TIDW1RY.clear();
TIDW2RY.clear();
TIDW3RY.clear();
TIDW1SX.clear();
TIDW2SX.clear();
TIDW3SX.clear();
TIDW1SY.clear();
TIDW2SY.clear();
TIDW3SY.clear();
TECW1RX.clear();
TECW2RX.clear();
TECW3RX.clear();
TECW4RX.clear();
TECW5RX.clear();
TECW6RX.clear();
TECW7RX.clear();
TECW8RX.clear();
TECW1RY.clear();
TECW2RY.clear();
TECW3RY.clear();
TECW4RY.clear();
TECW5RY.clear();
TECW6RY.clear();
TECW7RY.clear();
TECW8RY.clear();
TECW1SX.clear();
TECW2SX.clear();
TECW3SX.clear();
TECW4SX.clear();
TECW5SX.clear();
TECW6SX.clear();
TECW7SX.clear();
TECW8SX.clear();
TECW1SY.clear();
TECW2SY.clear();
TECW3SY.clear();
TECW4SY.clear();
TECW5SY.clear();
TECW6SY.clear();
TECW7SY.clear();
TECW8SY.clear();
BRL1RX.clear();
BRL1RY.clear();
BRL1SX.clear();
BRL1SY.clear();
BRL2RX.clear();
BRL2RY.clear();
BRL2SX.clear();
BRL2SY.clear();
BRL3RX.clear();
BRL3RY.clear();
BRL3SX.clear();
BRL3SY.clear();
FWD1pRX.clear();
FWD1pRY.clear();
FWD1pSX.clear();
FWD1pSY.clear();
FWD1nRX.clear();
FWD1nRY.clear();
FWD1nSX.clear();
FWD1nSY.clear();
FWD2pRX.clear();
FWD2pRY.clear();
FWD2pSX.clear();
FWD2pSY.clear();
FWD2nRX.clear();
FWD2nRY.clear();
FWD2nSX.clear();
FWD2nSY.clear();
//muon clear
DTRHD.clear();
DTSHD.clear();
CSCRHPHI.clear();
CSCRHPERP.clear();
CSCSHPHI.clear();
RPCRHX.clear();
RPCSHX.clear();
return;
}
| int GlobalRecHitsProducer::compute | ( | const DTGeometry * | dtGeom, |
| std::map< DTWireId, std::vector< PSimHit > > | simHitsPerWire, | ||
| std::map< DTWireId, std::vector< type > > | recHitsPerWire, | ||
| int | step | ||
| ) | [private] |
Definition at line 2184 of file GlobalRecHitsProducer.cc.
References DTRHD, DTSHD, findBestRecHit(), DTHitQualityUtils::findMuSimHit(), DTGeometry::layer(), Association::map, recHitDistFromWire(), and simHitDistFromWire().
Referenced by fillMuon().
{
int nDt = 0;
// Loop over cells with a muon SimHit
for(std::map<DTWireId, std::vector<PSimHit> >::const_iterator wireAndSHits =
simHitsPerWire.begin();
wireAndSHits != simHitsPerWire.end();
wireAndSHits++) {
DTWireId wireId = (*wireAndSHits).first;
std::vector<PSimHit> simHitsInCell = (*wireAndSHits).second;
// Get the layer
const DTLayer* layer = dtGeom->layer(wireId);
// Look for a mu hit in the cell
const PSimHit* muSimHit = DTHitQualityUtils::findMuSimHit(simHitsInCell);
if (muSimHit==0) {
continue; // Skip this cell
}
// Find the distance of the simhit from the wire
float simHitWireDist = simHitDistFromWire(layer, wireId, *muSimHit);
// Skip simhits out of the cell
if(simHitWireDist>2.1) {
continue; // Skip this cell
}
//GlobalPoint simHitGlobalPos = layer->toGlobal(muSimHit->localPosition());
// Look for RecHits in the same cell
if(recHitsPerWire.find(wireId) == recHitsPerWire.end()) {
continue; // No RecHit found in this cell
} else {
// vector<type> recHits = (*wireAndRecHits).second;
std::vector<type> recHits = recHitsPerWire[wireId];
// Find the best RecHit
const type* theBestRecHit =
findBestRecHit(layer, wireId, recHits, simHitWireDist);
float recHitWireDist = recHitDistFromWire(*theBestRecHit, layer);
++nDt;
DTRHD.push_back(recHitWireDist);
DTSHD.push_back(simHitWireDist);
} // find rechits
} // loop over simhits
return nDt;
}
| void GlobalRecHitsProducer::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDProducer.
Definition at line 109 of file GlobalRecHitsProducer.cc.
References count, and verbosity.
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_endJob";
if (verbosity >= 0)
edm::LogInfo(MsgLoggerCat)
<< "Terminating having processed " << count << " events.";
return;
}
| void GlobalRecHitsProducer::fillECal | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private] |
Definition at line 220 of file GlobalRecHitsProducer.cc.
References edm::SortedCollection< T, SORT >::begin(), EBRE, EBSHE, ECalEBSrc_, ECalEESrc_, ECalESSrc_, ECalUncalEBSrc_, ECalUncalEESrc_, EERE, EESHE, edm::SortedCollection< T, SORT >::end(), ESRE, ESSHE, edm::SortedCollection< T, SORT >::find(), edm::Event::getByLabel(), edm::HandleBase::isValid(), edm::Handle< T >::product(), DetId::rawId(), and verbosity.
Referenced by produce().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_fillECal";
TString eventout;
if (verbosity > 0)
eventout = "\nGathering info:";
// extract crossing frame from event
//edm::Handle<CrossingFrame> crossingFrame;
edm::Handle<CrossingFrame<PCaloHit> > crossingFrame;
//iEvent.getByType(crossingFrame);
//if (!crossingFrame.isValid()) {
// edm::LogWarning(MsgLoggerCat)
// << "Unable to find crossingFrame in event!";
// return;
//}
//extract EB information
edm::Handle<EBUncalibratedRecHitCollection> EcalUncalibRecHitEB;
iEvent.getByLabel(ECalUncalEBSrc_, EcalUncalibRecHitEB);
if (!EcalUncalibRecHitEB.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find EcalUncalRecHitEB in event!";
return;
}
edm::Handle<EBRecHitCollection> EcalRecHitEB;
iEvent.getByLabel(ECalEBSrc_, EcalRecHitEB);
if (!EcalRecHitEB.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find EcalRecHitEB in event!";
return;
}
// loop over simhits
const std::string barrelHitsName("EcalHitsEB");
iEvent.getByLabel("mix",barrelHitsName,crossingFrame);
if (!crossingFrame.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find cal barrel crossingFrame in event!";
return;
}
//std::auto_ptr<MixCollection<PCaloHit> >
// barrelHits(new MixCollection<PCaloHit>
// (crossingFrame.product(), barrelHitsName));
std::auto_ptr<MixCollection<PCaloHit> >
barrelHits(new MixCollection<PCaloHit>(crossingFrame.product()));
// keep track of sum of simhit energy in each crystal
MapType ebSimMap;
for (MixCollection<PCaloHit>::MixItr hitItr
= barrelHits->begin();
hitItr != barrelHits->end();
++hitItr) {
EBDetId ebid = EBDetId(hitItr->id());
uint32_t crystid = ebid.rawId();
ebSimMap[crystid] += hitItr->energy();
}
int nEBRecHits = 0;
// loop over RecHits
const EBUncalibratedRecHitCollection *EBUncalibRecHit =
EcalUncalibRecHitEB.product();
const EBRecHitCollection *EBRecHit = EcalRecHitEB.product();
for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit =
EBUncalibRecHit->begin();
uncalibRecHit != EBUncalibRecHit->end();
++uncalibRecHit) {
EBDetId EBid = EBDetId(uncalibRecHit->id());
EcalRecHitCollection::const_iterator myRecHit = EBRecHit->find(EBid);
if (myRecHit != EBRecHit->end()) {
++nEBRecHits;
EBRE.push_back(myRecHit->energy());
EBSHE.push_back(ebSimMap[EBid.rawId()]);
}
}
if (verbosity > 1) {
eventout += "\n Number of EBRecHits collected:............ ";
eventout += nEBRecHits;
}
//extract EE information
edm::Handle<EEUncalibratedRecHitCollection> EcalUncalibRecHitEE;
iEvent.getByLabel(ECalUncalEESrc_, EcalUncalibRecHitEE);
if (!EcalUncalibRecHitEE.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find EcalUncalRecHitEE in event!";
return;
}
edm::Handle<EERecHitCollection> EcalRecHitEE;
iEvent.getByLabel(ECalEESrc_, EcalRecHitEE);
if (!EcalRecHitEE.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find EcalRecHitEE in event!";
return;
}
// loop over simhits
const std::string endcapHitsName("EcalHitsEE");
iEvent.getByLabel("mix",endcapHitsName,crossingFrame);
if (!crossingFrame.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find cal endcap crossingFrame in event!";
return;
}
//std::auto_ptr<MixCollection<PCaloHit> >
// endcapHits(new MixCollection<PCaloHit>
// (crossingFrame.product(), endcapHitsName));
std::auto_ptr<MixCollection<PCaloHit> >
endcapHits(new MixCollection<PCaloHit>(crossingFrame.product()));
// keep track of sum of simhit energy in each crystal
MapType eeSimMap;
for (MixCollection<PCaloHit>::MixItr hitItr
= endcapHits->begin();
hitItr != endcapHits->end();
++hitItr) {
EEDetId eeid = EEDetId(hitItr->id());
uint32_t crystid = eeid.rawId();
eeSimMap[crystid] += hitItr->energy();
}
int nEERecHits = 0;
// loop over RecHits
const EEUncalibratedRecHitCollection *EEUncalibRecHit =
EcalUncalibRecHitEE.product();
const EERecHitCollection *EERecHit = EcalRecHitEE.product();
for (EcalUncalibratedRecHitCollection::const_iterator uncalibRecHit =
EEUncalibRecHit->begin();
uncalibRecHit != EEUncalibRecHit->end();
++uncalibRecHit) {
EEDetId EEid = EEDetId(uncalibRecHit->id());
EcalRecHitCollection::const_iterator myRecHit = EERecHit->find(EEid);
if (myRecHit != EERecHit->end()) {
++nEERecHits;
EERE.push_back(myRecHit->energy());
EESHE.push_back(eeSimMap[EEid.rawId()]);
}
}
if (verbosity > 1) {
eventout += "\n Number of EERecHits collected:............ ";
eventout += nEERecHits;
}
//extract ES information
edm::Handle<ESRecHitCollection> EcalRecHitES;
iEvent.getByLabel(ECalESSrc_, EcalRecHitES);
if (!EcalRecHitES.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find EcalRecHitES in event!";
return;
}
// loop over simhits
const std::string preshowerHitsName("EcalHitsES");
iEvent.getByLabel("mix",preshowerHitsName,crossingFrame);
if (!crossingFrame.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find cal preshower crossingFrame in event!";
return;
}
//std::auto_ptr<MixCollection<PCaloHit> >
// preshowerHits(new MixCollection<PCaloHit>
// (crossingFrame.product(), preshowerHitsName));
std::auto_ptr<MixCollection<PCaloHit> >
preshowerHits(new MixCollection<PCaloHit>(crossingFrame.product()));
// keep track of sum of simhit energy in each crystal
MapType esSimMap;
for (MixCollection<PCaloHit>::MixItr hitItr
= preshowerHits->begin();
hitItr != preshowerHits->end();
++hitItr) {
ESDetId esid = ESDetId(hitItr->id());
uint32_t crystid = esid.rawId();
esSimMap[crystid] += hitItr->energy();
}
int nESRecHits = 0;
// loop over RecHits
const ESRecHitCollection *ESRecHit = EcalRecHitES.product();
for (EcalRecHitCollection::const_iterator recHit =
ESRecHit->begin();
recHit != ESRecHit->end();
++recHit) {
ESDetId ESid = ESDetId(recHit->id());
++nESRecHits;
ESRE.push_back(recHit->energy());
ESSHE.push_back(esSimMap[ESid.rawId()]);
}
if (verbosity > 1) {
eventout += "\n Number of ESRecHits collected:............ ";
eventout += nESRecHits;
}
if (verbosity > 0)
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
return;
}
| void GlobalRecHitsProducer::fillHCal | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private] |
Definition at line 491 of file GlobalRecHitsProducer.cc.
References PV3DBase< T, PVType, FrameType >::eta(), geometry, edm::EventSetup::get(), edm::Event::getByLabel(), edm::Event::getManyByType(), CaloCellGeometry::getPosition(), HBCalR, HBCalREC, HBCalSHE, HCalSrc_, HECalR, HECalREC, HECalSHE, HFCalR, HFCalREC, HFCalSHE, HOCalR, HOCalREC, HOCalSHE, edm::HandleBase::isValid(), edm::ESHandleBase::isValid(), PV3DBase< T, PVType, FrameType >::phi(), PI, edm::Handle< T >::product(), csvReporter::r, DetId::rawId(), sdHcalBrl, sdHcalEC, sdHcalFwd, sdHcalOut, mathSSE::sqrt(), and verbosity.
Referenced by produce().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_fillHCal";
TString eventout;
if (verbosity > 0)
eventout = "\nGathering info:";
// get geometry
edm::ESHandle<CaloGeometry> geometry;
iSetup.get<CaloGeometryRecord>().get(geometry);
if (!geometry.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find CaloGeometry in event!";
return;
}
// extract simhit info
edm::Handle<edm::PCaloHitContainer> hcalHits;
iEvent.getByLabel(HCalSrc_,hcalHits);
if (!hcalHits.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find hcalHits in event!";
return;
}
const edm::PCaloHitContainer *simhitResult = hcalHits.product();
MapType fHBEnergySimHits;
MapType fHEEnergySimHits;
MapType fHOEnergySimHits;
MapType fHFEnergySimHits;
for (std::vector<PCaloHit>::const_iterator simhits = simhitResult->begin();
simhits != simhitResult->end();
++simhits) {
HcalDetId detId(simhits->id());
uint32_t cellid = detId.rawId();
if (detId.subdet() == sdHcalBrl){
fHBEnergySimHits[cellid] += simhits->energy();
}
if (detId.subdet() == sdHcalEC){
fHEEnergySimHits[cellid] += simhits->energy();
}
if (detId.subdet() == sdHcalOut){
fHOEnergySimHits[cellid] += simhits->energy();
}
if (detId.subdet() == sdHcalFwd){
fHFEnergySimHits[cellid] += simhits->energy();
}
}
// max values to be used (HO is found in HB)
Double_t maxHBEnergy = 0.;
Double_t maxHEEnergy = 0.;
Double_t maxHOEnergy = 0.;
Double_t maxHFEnergy = 0.;
Double_t maxHBPhi = -1000.;
Double_t maxHEPhi = -1000.;
Double_t maxHOPhi = -1000.;
Double_t maxHFPhi = -1000.;
Double_t maxHBEta = -1000.;
Double_t maxHEEta = -1000.;
Double_t maxHOEta = -1000.;
Double_t maxHFEta = -1000.;
Double_t PI = 3.141592653589;
// get HBHE information
std::vector<edm::Handle<HBHERecHitCollection> > hbhe;
iEvent.getManyByType(hbhe);
if (!hbhe[0].isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find any HBHERecHitCollections in event!";
return;
}
std::vector<edm::Handle<HBHERecHitCollection> >::iterator ihbhe;
int iHB = 0;
int iHE = 0;
for (ihbhe = hbhe.begin(); ihbhe != hbhe.end(); ++ihbhe) {
// find max values
for (HBHERecHitCollection::const_iterator jhbhe = (*ihbhe)->begin();
jhbhe != (*ihbhe)->end(); ++jhbhe) {
HcalDetId cell(jhbhe->id());
if (cell.subdet() == sdHcalBrl) {
const CaloCellGeometry* cellGeometry =
geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
double fEta = cellGeometry->getPosition().eta () ;
double fPhi = cellGeometry->getPosition().phi () ;
if ( (jhbhe->energy()) > maxHBEnergy ) {
maxHBEnergy = jhbhe->energy();
maxHOEnergy = maxHBEnergy;
maxHBPhi = fPhi;
maxHOPhi = maxHBPhi;
maxHBEta = fEta;
maxHOEta = maxHBEta;
}
}
if (cell.subdet() == sdHcalEC) {
const CaloCellGeometry* cellGeometry =
geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
double fEta = cellGeometry->getPosition().eta () ;
double fPhi = cellGeometry->getPosition().phi () ;
if ( (jhbhe->energy()) > maxHEEnergy ) {
maxHEEnergy = jhbhe->energy();
maxHEPhi = fPhi;
maxHEEta = fEta;
}
}
} // end find max values
for (HBHERecHitCollection::const_iterator jhbhe = (*ihbhe)->begin();
jhbhe != (*ihbhe)->end(); ++jhbhe) {
HcalDetId cell(jhbhe->id());
if (cell.subdet() == sdHcalBrl) {
++iHB;
const CaloCellGeometry* cellGeometry =
geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
double fEta = cellGeometry->getPosition().eta () ;
double fPhi = cellGeometry->getPosition().phi () ;
float deltaphi = maxHBPhi - fPhi;
if (fPhi > maxHBPhi) { deltaphi = fPhi - maxHBPhi;}
if (deltaphi > PI) { deltaphi = 2.0 * PI - deltaphi;}
float deltaeta = fEta - maxHBEta;
Double_t r = sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
HBCalREC.push_back(jhbhe->energy());
HBCalR.push_back(r);
HBCalSHE.push_back(fHBEnergySimHits[cell.rawId()]);
}
if (cell.subdet() == sdHcalEC) {
++iHE;
const CaloCellGeometry* cellGeometry =
geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
double fEta = cellGeometry->getPosition().eta () ;
double fPhi = cellGeometry->getPosition().phi () ;
float deltaphi = maxHEPhi - fPhi;
if (fPhi > maxHEPhi) { deltaphi = fPhi - maxHEPhi;}
if (deltaphi > PI) { deltaphi = 2.0 * PI - deltaphi;}
float deltaeta = fEta - maxHEEta;
Double_t r = sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
HECalREC.push_back(jhbhe->energy());
HECalR.push_back(r);
HECalSHE.push_back(fHEEnergySimHits[cell.rawId()]);
}
}
} // end loop through collection
if (verbosity > 1) {
eventout += "\n Number of HBRecHits collected:............ ";
eventout += iHB;
}
if (verbosity > 1) {
eventout += "\n Number of HERecHits collected:............ ";
eventout += iHE;
}
// get HF information
std::vector<edm::Handle<HFRecHitCollection> > hf;
iEvent.getManyByType(hf);
if (!hf[0].isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find any HFRecHitCollections in event!";
return;
}
std::vector<edm::Handle<HFRecHitCollection> >::iterator ihf;
int iHF = 0;
for (ihf = hf.begin(); ihf != hf.end(); ++ihf) {
// find max values
for (HFRecHitCollection::const_iterator jhf = (*ihf)->begin();
jhf != (*ihf)->end(); ++jhf) {
HcalDetId cell(jhf->id());
if (cell.subdet() == sdHcalFwd) {
const CaloCellGeometry* cellGeometry =
geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
double fEta = cellGeometry->getPosition().eta () ;
double fPhi = cellGeometry->getPosition().phi () ;
if ( (jhf->energy()) > maxHFEnergy ) {
maxHFEnergy = jhf->energy();
maxHFPhi = fPhi;
maxHFEta = fEta;
}
}
} // end find max values
for (HFRecHitCollection::const_iterator jhf = (*ihf)->begin();
jhf != (*ihf)->end(); ++jhf) {
HcalDetId cell(jhf->id());
if (cell.subdet() == sdHcalFwd) {
++iHF;
const CaloCellGeometry* cellGeometry =
geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
double fEta = cellGeometry->getPosition().eta () ;
double fPhi = cellGeometry->getPosition().phi () ;
float deltaphi = maxHBPhi - fPhi;
if (fPhi > maxHFPhi) { deltaphi = fPhi - maxHFPhi;}
if (deltaphi > PI) { deltaphi = 2.0 * PI - deltaphi;}
float deltaeta = fEta - maxHFEta;
Double_t r = sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
HFCalREC.push_back(jhf->energy());
HFCalR.push_back(r);
HFCalSHE.push_back(fHFEnergySimHits[cell.rawId()]);
}
}
} // end loop through collection
if (verbosity > 1) {
eventout += "\n Number of HFDigis collected:.............. ";
eventout += iHF;
}
// get HO information
std::vector<edm::Handle<HORecHitCollection> > ho;
iEvent.getManyByType(ho);
if (!ho[0].isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find any HORecHitCollections in event!";
return;
}
std::vector<edm::Handle<HORecHitCollection> >::iterator iho;
int iHO = 0;
for (iho = ho.begin(); iho != ho.end(); ++iho) {
for (HORecHitCollection::const_iterator jho = (*iho)->begin();
jho != (*iho)->end(); ++jho) {
HcalDetId cell(jho->id());
if (cell.subdet() == sdHcalOut) {
++iHO;
const CaloCellGeometry* cellGeometry =
geometry->getSubdetectorGeometry (cell)->getGeometry (cell) ;
double fEta = cellGeometry->getPosition().eta () ;
double fPhi = cellGeometry->getPosition().phi () ;
float deltaphi = maxHOPhi - fPhi;
if (fPhi > maxHOPhi) { deltaphi = fPhi - maxHOPhi;}
if (deltaphi > PI) { deltaphi = 2.0 * PI - deltaphi;}
float deltaeta = fEta - maxHOEta;
Double_t r = sqrt(deltaeta * deltaeta + deltaphi * deltaphi);
HOCalREC.push_back(jho->energy());
HOCalR.push_back(r);
HOCalSHE.push_back(fHOEnergySimHits[cell.rawId()]);
}
}
} // end loop through collection
if (verbosity > 1) {
eventout += "\n Number of HODigis collected:.............. ";
eventout += iHO;
}
if (verbosity > 0)
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
return;
}
| void GlobalRecHitsProducer::fillMuon | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private] |
Definition at line 1651 of file GlobalRecHitsProducer.cc.
References CSCLayer::chamber(), CSCChamberSpecs::chamberType(), compute(), CSCDetId, edm::EventSetup::get(), edm::Event::getByLabel(), i, CSCGeometry::idToDetUnit(), RPCRoll::isForward(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), map1DRecHitsPerWire(), DTHitQualityUtils::mapSimHitsPerWire(), MuCSCSrc_, MuDTSimSrc_, MuDTSrc_, MuRPCSimSrc_, MuRPCSrc_, plotResolution(), edm::ESHandle< T >::product(), edm::Handle< T >::product(), csvReporter::r, RPCDetId, RPCRHX, RPCSHX, trackerHits::simHits, CSCChamber::specs(), theMap, verbosity, and PV3DBase< T, PVType, FrameType >::x().
Referenced by produce().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_fillMuon";
TString eventout;
if (verbosity > 0)
eventout = "\nGathering info:";
// get DT information
edm::ESHandle<DTGeometry> dtGeom;
iSetup.get<MuonGeometryRecord>().get(dtGeom);
if (!dtGeom.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find DTMuonGeometryRecord in event!";
return;
}
edm::Handle<edm::PSimHitContainer> dtsimHits;
iEvent.getByLabel(MuDTSimSrc_, dtsimHits);
if (!dtsimHits.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find dtsimHits in event!";
return;
}
std::map<DTWireId, edm::PSimHitContainer> simHitsPerWire =
DTHitQualityUtils::mapSimHitsPerWire(*(dtsimHits.product()));
edm::Handle<DTRecHitCollection> dtRecHits;
iEvent.getByLabel(MuDTSrc_, dtRecHits);
if (!dtRecHits.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find dtRecHits in event!";
return;
}
std::map<DTWireId, std::vector<DTRecHit1DPair> > recHitsPerWire =
map1DRecHitsPerWire(dtRecHits.product());
int nDt = compute(dtGeom.product(), simHitsPerWire, recHitsPerWire, 1);
if (verbosity > 1) {
eventout += "\n Number of DtMuonRecHits collected:........ ";
eventout += nDt;
}
// get CSC Strip information
// get map of sim hits
theMap.clear();
//edm::Handle<CrossingFrame> cf;
edm::Handle<CrossingFrame<PSimHit> > cf;
//iEvent.getByType(cf);
//if (!cf.isValid()) {
// edm::LogWarning(MsgLoggerCat)
// << "Unable to find CrossingFrame in event!";
// return;
//}
//MixCollection<PSimHit> simHits(cf.product(), "MuonCSCHits");
iEvent.getByLabel("mix","MuonCSCHits",cf);
if (!cf.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find muo CSC crossingFrame in event!";
return;
}
MixCollection<PSimHit> simHits(cf.product());
// arrange the hits by detUnit
for(MixCollection<PSimHit>::MixItr hitItr = simHits.begin();
hitItr != simHits.end(); ++hitItr)
{
theMap[hitItr->detUnitId()].push_back(*hitItr);
}
// get geometry
edm::ESHandle<CSCGeometry> hGeom;
iSetup.get<MuonGeometryRecord>().get(hGeom);
if (!hGeom.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find CSCMuonGeometryRecord in event!";
return;
}
const CSCGeometry *theCSCGeometry = &*hGeom;
// get rechits
edm::Handle<CSCRecHit2DCollection> hRecHits;
iEvent.getByLabel(MuCSCSrc_, hRecHits);
if (!hRecHits.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find CSC RecHits in event!";
return;
}
const CSCRecHit2DCollection *cscRecHits = hRecHits.product();
int nCSC = 0;
for (CSCRecHit2DCollection::const_iterator recHitItr = cscRecHits->begin();
recHitItr != cscRecHits->end(); ++recHitItr) {
int detId = (*recHitItr).cscDetId().rawId();
edm::PSimHitContainer simHits;
std::map<int, edm::PSimHitContainer>::const_iterator mapItr =
theMap.find(detId);
if (mapItr != theMap.end()) {
simHits = mapItr->second;
}
if (simHits.size() == 1) {
++nCSC;
const GeomDetUnit* detUnit =
theCSCGeometry->idToDetUnit(CSCDetId(detId));
const CSCLayer *layer = dynamic_cast<const CSCLayer *>(detUnit);
int chamberType = layer->chamber()->specs()->chamberType();
plotResolution(simHits[0], *recHitItr, layer, chamberType);
}
}
if (verbosity > 1) {
eventout += "\n Number of CSCRecHits collected:........... ";
eventout += nCSC;
}
// get RPC information
std::map<double, int> mapsim, maprec;
std::map<int, double> nmapsim, nmaprec;
edm::ESHandle<RPCGeometry> rpcGeom;
iSetup.get<MuonGeometryRecord>().get(rpcGeom);
if (!rpcGeom.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find RPCMuonGeometryRecord in event!";
return;
}
edm::Handle<edm::PSimHitContainer> simHit;
iEvent.getByLabel(MuRPCSimSrc_, simHit);
if (!simHit.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find RPCSimHit in event!";
return;
}
edm::Handle<RPCRecHitCollection> recHit;
iEvent.getByLabel(MuRPCSrc_, recHit);
if (!simHit.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find RPCRecHit in event!";
return;
}
int nRPC = 0;
RPCRecHitCollection::const_iterator recIt;
int nrec = 0;
for (recIt = recHit->begin(); recIt != recHit->end(); ++recIt) {
RPCDetId Rid = (RPCDetId)(*recIt).rpcId();
const RPCRoll *roll = dynamic_cast<const RPCRoll*>(rpcGeom->roll(Rid));
if (roll->isForward()) {
if (verbosity > 1) {
eventout += "\n Number of RPCRecHits collected:........... ";
eventout += nRPC;
}
if (verbosity > 0)
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
return;
}
nrec = nrec + 1;
LocalPoint rhitlocal = (*recIt).localPosition();
double rhitlocalx = rhitlocal.x();
maprec[rhitlocalx] = nrec;
}
int i = 0;
for (std::map<double,int>::iterator iter = maprec.begin();
iter != maprec.end(); ++iter) {
i = i + 1;
nmaprec[i] = (*iter).first;
}
edm::PSimHitContainer::const_iterator simIt;
int nsim = 0;
for (simIt = simHit->begin(); simIt != simHit->end(); simIt++) {
int ptype = (*simIt).particleType();
//RPCDetId Rsid = (RPCDetId)(*simIt).detUnitId();
if (ptype == 13 || ptype == -13) {
nsim = nsim + 1;
LocalPoint shitlocal = (*simIt).localPosition();
double shitlocalx = shitlocal.x();
mapsim[shitlocalx] = nsim;
}
}
i = 0;
for (std::map<double,int>::iterator iter = mapsim.begin();
iter != mapsim.end(); ++iter) {
i = i + 1;
nmapsim[i] = (*iter).first;
}
if (nsim == nrec) {
for (int r = 0; r < nsim; r++) {
++nRPC;
RPCRHX.push_back(nmaprec[r+1]);
RPCSHX.push_back(nmapsim[r+1]);
}
}
if (verbosity > 1) {
eventout += "\n Number of RPCRecHits collected:........... ";
eventout += nRPC;
}
if (verbosity > 0)
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
return;
}
| void GlobalRecHitsProducer::fillTrk | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [private] |
Definition at line 855 of file GlobalRecHitsProducer.cc.
References TrackerHitAssociator::associateHit(), edmNew::DetSet< T >::begin(), BRL1RX, BRL1RY, BRL1SX, BRL1SY, BRL2RX, BRL2RY, BRL2SX, BRL2SY, BRL3RX, BRL3RY, BRL3SX, BRL3SY, conf_, cond::rpcobgas::detid, PXFDetId::disk(), edmNew::DetSet< T >::end(), FWD1nRX, FWD1nRY, FWD1nSX, FWD1nSY, FWD1pRX, FWD1pRY, FWD1pSX, FWD1pSY, FWD2nRX, FWD2nRY, FWD2nSX, FWD2nSY, FWD2pRX, FWD2pRY, FWD2pSX, FWD2pSY, TrackingRecHit::geographicalId(), relativeConstraints::geom, edm::EventSetup::get(), edm::Event::getByLabel(), TrackerGeometry::idToDet(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), TIBDetId::layer(), TOBDetId::layer(), PXBDetId::layer(), BaseSiTrackerRecHit2DLocalPos::localPosition(), m, matched, position, projectHit(), sdPxlBrl, sdPxlFwd, sdSiTEC, sdSiTIB, sdSiTID, sdSiTOB, PXFDetId::side(), SiPxlSrc_, SiStripSrc_, mathSSE::sqrt(), GluedGeomDet::stereoDet(), DetId::subdetId(), GeomDet::surface(), TECW1RX, TECW1RY, TECW1SX, TECW1SY, TECW2RX, TECW2RY, TECW2SX, TECW2SY, TECW3RX, TECW3RY, TECW3SX, TECW3SY, TECW4RX, TECW4RY, TECW4SX, TECW4SY, TECW5RX, TECW5RY, TECW5SX, TECW5SY, TECW6RX, TECW6RY, TECW6SX, TECW6SY, TECW7RX, TECW7RY, TECW7SX, TECW7SY, TECW8RX, TECW8RY, TECW8SX, TECW8SY, TIBL1RX, TIBL1RY, TIBL1SX, TIBL1SY, TIBL2RX, TIBL2RY, TIBL2SX, TIBL2SY, TIBL3RX, TIBL3RY, TIBL3SX, TIBL3SY, TIBL4RX, TIBL4RY, TIBL4SX, TIBL4SY, TIDW1RX, TIDW1RY, TIDW1SX, TIDW1SY, TIDW2RX, TIDW2RY, TIDW2SX, TIDW2SY, TIDW3RX, TIDW3RY, TIDW3SX, TIDW3SY, TOBL1RX, TOBL1RY, TOBL1SX, TOBL1SY, TOBL2RX, TOBL2RY, TOBL2SX, TOBL2SY, TOBL3RX, TOBL3RY, TOBL3SX, TOBL3SY, TOBL4RX, TOBL4RY, TOBL4SX, TOBL4SY, patCandidatesForDimuonsSequences_cff::tracker, verbosity, TIDDetId::wheel(), TECDetId::wheel(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().
Referenced by produce().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_fillTrk";
TString eventout;
if (verbosity > 0)
eventout = "\nGathering info:";
// get strip information
edm::Handle<SiStripMatchedRecHit2DCollection> rechitsmatched;
iEvent.getByLabel(SiStripSrc_, rechitsmatched);
if (!rechitsmatched.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find stripmatchedrechits in event!";
return;
}
TrackerHitAssociator associate(iEvent,conf_);
edm::ESHandle<TrackerGeometry> pDD;
iSetup.get<TrackerDigiGeometryRecord>().get(pDD);
if (!pDD.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find TrackerDigiGeometry in event!";
return;
}
const TrackerGeometry &tracker(*pDD);
int nStripBrl = 0, nStripFwd = 0;
// loop over det units
for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin();
it != pDD->dets().end(); ++it) {
uint32_t myid = ((*it)->geographicalId()).rawId();
DetId detid = ((*it)->geographicalId());
//loop over rechits-matched in the same subdetector
SiStripMatchedRecHit2DCollection::const_iterator rechitmatchedMatch = rechitsmatched->find(detid);
if (rechitmatchedMatch != rechitsmatched->end()) {
SiStripMatchedRecHit2DCollection::DetSet rechitmatchedRange = *rechitmatchedMatch;
SiStripMatchedRecHit2DCollection::DetSet::const_iterator rechitmatchedRangeIteratorBegin = rechitmatchedRange.begin();
SiStripMatchedRecHit2DCollection::DetSet::const_iterator rechitmatchedRangeIteratorEnd = rechitmatchedRange.end();
SiStripMatchedRecHit2DCollection::DetSet::const_iterator itermatched = rechitmatchedRangeIteratorBegin;
for ( itermatched = rechitmatchedRangeIteratorBegin;
itermatched != rechitmatchedRangeIteratorEnd;
++itermatched) {
SiStripMatchedRecHit2D const rechit = *itermatched;
LocalPoint position = rechit.localPosition();
float mindist = 999999.;
float distx = 999999.;
float disty = 999999.;
float dist = 999999.;
std::pair<LocalPoint,LocalVector> closestPair;
matched.clear();
float rechitmatchedx = position.x();
float rechitmatchedy = position.y();
matched = associate.associateHit(rechit);
if (!matched.empty()) {
//project simhit;
const GluedGeomDet* gluedDet =
(const GluedGeomDet*)tracker.idToDet(rechit.geographicalId());
const StripGeomDetUnit* partnerstripdet =
(StripGeomDetUnit*) gluedDet->stereoDet();
std::pair<LocalPoint,LocalVector> hitPair;
for(std::vector<PSimHit>::const_iterator m = matched.begin();
m != matched.end(); m++){
//project simhit;
hitPair = projectHit((*m),partnerstripdet,gluedDet->surface());
distx = fabs(rechitmatchedx - hitPair.first.x());
disty = fabs(rechitmatchedy - hitPair.first.y());
dist = sqrt(distx*distx+disty*disty);
if(dist < mindist){
mindist = dist;
closestPair = hitPair;
}
}
// get TIB
if (detid.subdetId() == sdSiTIB) {
TIBDetId tibid(myid);
++nStripBrl;
if (tibid.layer() == 1) {
TIBL1RX.push_back(rechitmatchedx);
TIBL1RY.push_back(rechitmatchedy);
TIBL1SX.push_back(closestPair.first.x());
TIBL1SY.push_back(closestPair.first.y());
}
if (tibid.layer() == 2) {
TIBL2RX.push_back(rechitmatchedx);
TIBL2RY.push_back(rechitmatchedy);
TIBL2SX.push_back(closestPair.first.x());
TIBL2SY.push_back(closestPair.first.y());
}
if (tibid.layer() == 3) {
TIBL3RX.push_back(rechitmatchedx);
TIBL3RY.push_back(rechitmatchedy);
TIBL3SX.push_back(closestPair.first.x());
TIBL3SY.push_back(closestPair.first.y());
}
if (tibid.layer() == 4) {
TIBL4RX.push_back(rechitmatchedx);
TIBL4RY.push_back(rechitmatchedy);
TIBL4SX.push_back(closestPair.first.x());
TIBL4SY.push_back(closestPair.first.y());
}
}
// get TOB
if (detid.subdetId() == sdSiTOB) {
TOBDetId tobid(myid);
++nStripBrl;
if (tobid.layer() == 1) {
TOBL1RX.push_back(rechitmatchedx);
TOBL1RY.push_back(rechitmatchedy);
TOBL1SX.push_back(closestPair.first.x());
TOBL1SY.push_back(closestPair.first.y());
}
if (tobid.layer() == 2) {
TOBL2RX.push_back(rechitmatchedx);
TOBL2RY.push_back(rechitmatchedy);
TOBL2SX.push_back(closestPair.first.x());
TOBL2SY.push_back(closestPair.first.y());
}
if (tobid.layer() == 3) {
TOBL3RX.push_back(rechitmatchedx);
TOBL3RY.push_back(rechitmatchedy);
TOBL3SX.push_back(closestPair.first.x());
TOBL3SY.push_back(closestPair.first.y());
}
if (tobid.layer() == 4) {
TOBL4RX.push_back(rechitmatchedx);
TOBL4RY.push_back(rechitmatchedy);
TOBL4SX.push_back(closestPair.first.x());
TOBL4SY.push_back(closestPair.first.y());
}
}
// get TID
if (detid.subdetId() == sdSiTID) {
TIDDetId tidid(myid);
++nStripFwd;
if (tidid.wheel() == 1) {
TIDW1RX.push_back(rechitmatchedx);
TIDW1RY.push_back(rechitmatchedy);
TIDW1SX.push_back(closestPair.first.x());
TIDW1SY.push_back(closestPair.first.y());
}
if (tidid.wheel() == 2) {
TIDW2RX.push_back(rechitmatchedx);
TIDW2RY.push_back(rechitmatchedy);
TIDW2SX.push_back(closestPair.first.x());
TIDW2SY.push_back(closestPair.first.y());
}
if (tidid.wheel() == 3) {
TIDW3RX.push_back(rechitmatchedx);
TIDW3RY.push_back(rechitmatchedy);
TIDW3SX.push_back(closestPair.first.x());
TIDW3SY.push_back(closestPair.first.y());
}
}
// get TEC
if (detid.subdetId() == sdSiTEC) {
TECDetId tecid(myid);
++nStripFwd;
if (tecid.wheel() == 1) {
TECW1RX.push_back(rechitmatchedx);
TECW1RY.push_back(rechitmatchedy);
TECW1SX.push_back(closestPair.first.x());
TECW1SY.push_back(closestPair.first.y());
}
if (tecid.wheel() == 2) {
TECW2RX.push_back(rechitmatchedx);
TECW2RY.push_back(rechitmatchedy);
TECW2SX.push_back(closestPair.first.x());
TECW2SY.push_back(closestPair.first.y());
}
if (tecid.wheel() == 3) {
TECW3RX.push_back(rechitmatchedx);
TECW3RY.push_back(rechitmatchedy);
TECW3SX.push_back(closestPair.first.x());
TECW3SY.push_back(closestPair.first.y());
}
if (tecid.wheel() == 4) {
TECW4RX.push_back(rechitmatchedx);
TECW4RY.push_back(rechitmatchedy);
TECW4SX.push_back(closestPair.first.x());
TECW4SY.push_back(closestPair.first.y());
}
if (tecid.wheel() == 5) {
TECW5RX.push_back(rechitmatchedx);
TECW5RY.push_back(rechitmatchedy);
TECW5SX.push_back(closestPair.first.x());
TECW5SY.push_back(closestPair.first.y());
}
if (tecid.wheel() == 6) {
TECW6RX.push_back(rechitmatchedx);
TECW6RY.push_back(rechitmatchedy);
TECW6SX.push_back(closestPair.first.x());
TECW6SY.push_back(closestPair.first.y());
}
if (tecid.wheel() == 7) {
TECW7RX.push_back(rechitmatchedx);
TECW7RY.push_back(rechitmatchedy);
TECW7SX.push_back(closestPair.first.x());
TECW7SY.push_back(closestPair.first.y());
}
if (tecid.wheel() == 8) {
TECW8RX.push_back(rechitmatchedx);
TECW8RY.push_back(rechitmatchedy);
TECW8SX.push_back(closestPair.first.x());
TECW8SY.push_back(closestPair.first.y());
}
}
} // end if matched empty
}
}
} // end loop over det units
if (verbosity > 1) {
eventout += "\n Number of BrlStripRecHits collected:...... ";
eventout += nStripBrl;
}
if (verbosity > 1) {
eventout += "\n Number of FrwdStripRecHits collected:..... ";
eventout += nStripFwd;
}
// get pixel information
//Get RecHits
edm::Handle<SiPixelRecHitCollection> recHitColl;
iEvent.getByLabel(SiPxlSrc_, recHitColl);
if (!recHitColl.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find SiPixelRecHitCollection in event!";
return;
}
//Get event setup
edm::ESHandle<TrackerGeometry> geom;
iSetup.get<TrackerDigiGeometryRecord>().get(geom);
if (!geom.isValid()) {
edm::LogWarning(MsgLoggerCat)
<< "Unable to find TrackerDigiGeometry in event!";
return;
}
//const TrackerGeometry& theTracker(*geom);
int nPxlBrl = 0, nPxlFwd = 0;
//iterate over detunits
for (TrackerGeometry::DetContainer::const_iterator it = geom->dets().begin();
it != geom->dets().end(); ++it) {
uint32_t myid = ((*it)->geographicalId()).rawId();
DetId detId = ((*it)->geographicalId());
int subid = detId.subdetId();
if (! ((subid == sdPxlBrl) || (subid == sdPxlFwd))) continue;
//const PixelGeomDetUnit * theGeomDet =
// dynamic_cast<const PixelGeomDetUnit*>(theTracker.idToDet(detId) );
SiPixelRecHitCollection::const_iterator pixeldet = recHitColl->find(detId);
if (pixeldet == recHitColl->end()) continue;
SiPixelRecHitCollection::DetSet pixelrechitRange = *pixeldet;
SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;
std::vector<PSimHit> matched;
//----Loop over rechits for this detId
for ( ; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter) {
matched.clear();
matched = associate.associateHit(*pixeliter);
if ( !matched.empty() ) {
float closest = 9999.9;
//std::vector<PSimHit>::const_iterator closestit = matched.begin();
LocalPoint lp = pixeliter->localPosition();
float rechit_x = lp.x();
float rechit_y = lp.y();
float sim_x = 0.;
float sim_y = 0.;
//loop over sim hits and fill closet
for (std::vector<PSimHit>::const_iterator m = matched.begin();
m != matched.end(); ++m) {
float sim_x1 = (*m).entryPoint().x();
float sim_x2 = (*m).exitPoint().x();
float sim_xpos = 0.5*(sim_x1+sim_x2);
float sim_y1 = (*m).entryPoint().y();
float sim_y2 = (*m).exitPoint().y();
float sim_ypos = 0.5*(sim_y1+sim_y2);
float x_res = fabs(sim_xpos - rechit_x);
float y_res = fabs(sim_ypos - rechit_y);
float dist = sqrt(x_res*x_res + y_res*y_res);
if ( dist < closest ) {
closest = dist;
sim_x = sim_xpos;
sim_y = sim_ypos;
}
} // end sim hit loop
// get Barrel pixels
if (subid == sdPxlBrl) {
PXBDetId bdetid(myid);
++nPxlBrl;
if (bdetid.layer() == 1) {
BRL1RX.push_back(rechit_x);
BRL1RY.push_back(rechit_y);
BRL1SX.push_back(sim_x);
BRL1SY.push_back(sim_y);
}
if (bdetid.layer() == 2) {
BRL2RX.push_back(rechit_x);
BRL2RY.push_back(rechit_y);
BRL2SX.push_back(sim_x);
BRL2SY.push_back(sim_y);
}
if (bdetid.layer() == 3) {
BRL3RX.push_back(rechit_x);
BRL3RY.push_back(rechit_y);
BRL3SX.push_back(sim_x);
BRL3SY.push_back(sim_y);
}
}
// get Forward pixels
if (subid == sdPxlFwd) {
PXFDetId fdetid(myid);
++nPxlFwd;
if (fdetid.disk() == 1) {
if (fdetid.side() == 1) {
FWD1nRX.push_back(rechit_x);
FWD1nRY.push_back(rechit_y);
FWD1nSX.push_back(sim_x);
FWD1nSY.push_back(sim_y);
}
if (fdetid.side() == 2) {
FWD1pRX.push_back(rechit_x);
FWD1pRY.push_back(rechit_y);
FWD1pSX.push_back(sim_x);
FWD1pSY.push_back(sim_y);
}
}
if (fdetid.disk() == 2) {
if (fdetid.side() == 1) {
FWD2nRX.push_back(rechit_x);
FWD2nRY.push_back(rechit_y);
FWD2nSX.push_back(sim_x);
FWD2nSY.push_back(sim_y);
}
if (fdetid.side() == 2) {
FWD2pRX.push_back(rechit_x);
FWD2pRY.push_back(rechit_y);
FWD2pSX.push_back(sim_x);
FWD2pSY.push_back(sim_y);
}
}
}
} // end matched emtpy
} // <-----end rechit loop
} // <------ end detunit loop
if (verbosity > 1) {
eventout += "\n Number of BrlPixelRecHits collected:...... ";
eventout += nPxlBrl;
}
if (verbosity > 1) {
eventout += "\n Number of FrwdPixelRecHits collected:..... ";
eventout += nPxlFwd;
}
if (verbosity > 0)
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
return;
}
| const type * GlobalRecHitsProducer::findBestRecHit | ( | const DTLayer * | layer, |
| DTWireId | wireId, | ||
| const std::vector< type > & | recHits, | ||
| const float | simHitDist | ||
| ) | [private] |
Definition at line 2146 of file GlobalRecHitsProducer.cc.
References recHitDistFromWire().
Referenced by compute().
{
float res = 99999;
const type* theBestRecHit = 0;
// Loop over RecHits within the cell
for(typename std::vector<type>::const_iterator recHit = recHits.begin();
recHit != recHits.end();
recHit++) {
float distTmp = recHitDistFromWire(*recHit, layer);
if(fabs(distTmp-simHitDist) < res) {
res = fabs(distTmp-simHitDist);
theBestRecHit = &(*recHit);
}
} // End of loop over RecHits within the cell
return theBestRecHit;
}
| std::map< DTWireId, std::vector< DTRecHit1DPair > > GlobalRecHitsProducer::map1DRecHitsPerWire | ( | const DTRecHitCollection * | dt1DRecHitPairs | ) | [private] |
Definition at line 2117 of file GlobalRecHitsProducer.cc.
References runTheMatrix::ret.
Referenced by fillMuon().
{
std::map<DTWireId, std::vector<DTRecHit1DPair> > ret;
for(DTRecHitCollection::const_iterator rechit = dt1DRecHitPairs->begin();
rechit != dt1DRecHitPairs->end(); rechit++) {
ret[(*rechit).wireId()].push_back(*rechit);
}
return ret;
}
| void GlobalRecHitsProducer::plotResolution | ( | const PSimHit & | simHit, |
| const CSCRecHit2D & | recHit, | ||
| const CSCLayer * | layer, | ||
| int | chamberType | ||
| ) | [private] |
Definition at line 2243 of file GlobalRecHitsProducer.cc.
References CSCRHPERP, CSCRHPHI, CSCSHPHI, PSimHit::localPosition(), CSCRecHit2D::localPosition(), PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), and GeomDet::toGlobal().
Referenced by fillMuon().
{
GlobalPoint simHitPos = layer->toGlobal(simHit.localPosition());
GlobalPoint recHitPos = layer->toGlobal(recHit.localPosition());
CSCRHPHI.push_back(recHitPos.phi());
CSCRHPERP.push_back(recHitPos.perp());
CSCSHPHI.push_back(simHitPos.phi());
}
| void GlobalRecHitsProducer::produce | ( | edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDProducer.
Definition at line 118 of file GlobalRecHitsProducer.cc.
References clear(), count, edm::EventID::event(), fillECal(), fillHCal(), fillMuon(), fillTrk(), frequency, edm::Event::getAllProvenance(), getAllProvenances, i, edm::EventBase::id(), label, printProvenanceInfo, edm::Event::put(), edm::EventID::run(), storeECal(), storeHCal(), storeMuon(), storeTrk(), and verbosity.
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_produce";
// keep track of number of events processed
++count;
// get event id information
int nrun = iEvent.id().run();
int nevt = iEvent.id().event();
if (verbosity > 0) {
edm::LogInfo(MsgLoggerCat)
<< "Processing run " << nrun << ", event " << nevt
<< " (" << count << " events total)";
} else if (verbosity == 0) {
if (nevt%frequency == 0 || nevt == 1) {
edm::LogInfo(MsgLoggerCat)
<< "Processing run " << nrun << ", event " << nevt
<< " (" << count << " events total)";
}
}
// clear event holders
clear();
// look at information available in the event
if (getAllProvenances) {
std::vector<const edm::Provenance*> AllProv;
iEvent.getAllProvenance(AllProv);
if (verbosity >= 0)
edm::LogInfo(MsgLoggerCat)
<< "Number of Provenances = " << AllProv.size();
if (printProvenanceInfo && (verbosity >= 0)) {
TString eventout("\nProvenance info:\n");
for (unsigned int i = 0; i < AllProv.size(); ++i) {
eventout += "\n ******************************";
eventout += "\n Module : ";
//eventout += (AllProv[i]->product).moduleLabel();
eventout += AllProv[i]->moduleLabel();
eventout += "\n ProductID : ";
//eventout += (AllProv[i]->product).productID_.id_;
eventout += AllProv[i]->productID().id();
eventout += "\n ClassName : ";
//eventout += (AllProv[i]->product).fullClassName_;
eventout += AllProv[i]->className();
eventout += "\n InstanceName : ";
//eventout += (AllProv[i]->product).productInstanceName_;
eventout += AllProv[i]->productInstanceName();
eventout += "\n BranchName : ";
//eventout += (AllProv[i]->product).branchName_;
eventout += AllProv[i]->branchName();
}
eventout += "\n ******************************\n";
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
printProvenanceInfo = false;
}
getAllProvenances = false;
}
// call fill functions
// gather Ecal information from event
fillECal(iEvent, iSetup);
// gather Hcal information from event
fillHCal(iEvent, iSetup);
// gather Track information from event
fillTrk(iEvent, iSetup);
// gather Muon information from event
fillMuon(iEvent, iSetup);
if (verbosity > 0)
edm::LogInfo (MsgLoggerCat)
<< "Done gathering data from event.";
// produce object to put into event
std::auto_ptr<PGlobalRecHit> pOut(new PGlobalRecHit);
if (verbosity > 2)
edm::LogInfo (MsgLoggerCat)
<< "Saving event contents:";
// call store functions
// store ECal information in produce
storeECal(*pOut);
// store HCal information in produce
storeHCal(*pOut);
// store Track information in produce
storeTrk(*pOut);
// store Muon information in produce
storeMuon(*pOut);
// store information in event
iEvent.put(pOut,label);
return;
}
| std::pair< LocalPoint, LocalVector > GlobalRecHitsProducer::projectHit | ( | const PSimHit & | hit, |
| const StripGeomDetUnit * | stripDet, | ||
| const BoundPlane & | plane | ||
| ) | [private] |
Definition at line 2086 of file GlobalRecHitsProducer.cc.
References funct::cos(), dir, PSimHit::localDirection(), PSimHit::localPosition(), funct::sin(), StripGeomDetUnit::specificTopology(), StripTopology::strip(), StripTopology::stripAngle(), GeomDet::surface(), Surface::toGlobal(), GloballyPositioned< T >::toLocal(), and PV3DBase< T, PVType, FrameType >::z().
Referenced by fillTrk().
{
const StripTopology& topol = stripDet->specificTopology();
GlobalPoint globalpos= stripDet->surface().toGlobal(hit.localPosition());
LocalPoint localHit = plane.toLocal(globalpos);
//track direction
LocalVector locdir=hit.localDirection();
//rotate track in new frame
GlobalVector globaldir= stripDet->surface().toGlobal(locdir);
LocalVector dir=plane.toLocal(globaldir);
float scale = -localHit.z() / dir.z();
LocalPoint projectedPos = localHit + scale*dir;
float selfAngle = topol.stripAngle( topol.strip( hit.localPosition()));
// vector along strip in hit frame
LocalVector stripDir( sin(selfAngle), cos(selfAngle), 0);
LocalVector
localStripDir(plane.toLocal(stripDet->surface().toGlobal(stripDir)));
return std::pair<LocalPoint,LocalVector>( projectedPos, localStripDir);
}
| float GlobalRecHitsProducer::recHitDistFromWire | ( | const DTRecHit1D & | recHit, |
| const DTLayer * | layer | ||
| ) | [private] |
Definition at line 2177 of file GlobalRecHitsProducer.cc.
References DTRecHit1D::localPosition(), DTLayer::specificTopology(), DTWireId::wire(), DTRecHit1D::wireId(), DTTopology::wirePosition(), and PV3DBase< T, PVType, FrameType >::x().
{
return fabs(recHit.localPosition().x() -
layer->specificTopology().wirePosition(recHit.wireId().wire()));
}
| float GlobalRecHitsProducer::recHitDistFromWire | ( | const DTRecHit1DPair & | hitPair, |
| const DTLayer * | layer | ||
| ) | [private] |
Definition at line 2168 of file GlobalRecHitsProducer.cc.
References DTEnums::Left, DTRecHit1DPair::localPosition(), DTEnums::Right, and PV3DBase< T, PVType, FrameType >::x().
Referenced by compute(), and findBestRecHit().
{
// Compute the rechit distance from wire
return fabs(hitPair.localPosition(DTEnums::Left).x() -
hitPair.localPosition(DTEnums::Right).x())/2.;
}
| float GlobalRecHitsProducer::simHitDistFromWire | ( | const DTLayer * | layer, |
| DTWireId | wireId, | ||
| const PSimHit & | hit | ||
| ) | [private] |
Definition at line 2130 of file GlobalRecHitsProducer.cc.
References PSimHit::entryPoint(), PSimHit::exitPoint(), DTLayer::specificTopology(), DTWireId::wire(), DTTopology::wirePosition(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().
Referenced by compute().
{
float xwire = layer->specificTopology().wirePosition(wireId.wire());
LocalPoint entryP = hit.entryPoint();
LocalPoint exitP = hit.exitPoint();
float xEntry = entryP.x()-xwire;
float xExit = exitP.x()-xwire;
//FIXME: check...
return fabs(xEntry - (entryP.z()*(xExit-xEntry))/(exitP.z()-entryP.z()));
}
| void GlobalRecHitsProducer::storeECal | ( | PGlobalRecHit & | product | ) | [private] |
Definition at line 449 of file GlobalRecHitsProducer.cc.
References EBRE, EBSHE, EERE, EESHE, ESRE, ESSHE, i, PGlobalRecHit::putEBCalRecHits(), PGlobalRecHit::putEECalRecHits(), PGlobalRecHit::putESCalRecHits(), and verbosity.
Referenced by produce().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_storeECal";
if (verbosity > 2) {
TString eventout("\n nEBRecHits = ");
eventout += EBRE.size();
for (unsigned int i = 0; i < EBRE.size(); ++i) {
eventout += "\n (RE, SHE) = (";
eventout += EBRE[i];
eventout += ", ";
eventout += EBSHE[i];
eventout += ")";
}
eventout += "\n nEERecHits = ";
eventout += EERE.size();
for (unsigned int i = 0; i < EERE.size(); ++i) {
eventout += "\n (RE, SHE) = (";
eventout += EERE[i];
eventout += ", ";
eventout += EESHE[i];
eventout += ")";
}
eventout += "\n nESRecHits = ";
eventout += ESRE.size();
for (unsigned int i = 0; i < ESRE.size(); ++i) {
eventout += "\n (RE, SHE) = (";
eventout += ESRE[i];
eventout += ", ";
eventout += ESSHE[i];
eventout += ")";
}
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
}
product.putEBCalRecHits(EBRE,EBSHE);
product.putEECalRecHits(EERE,EESHE);
product.putESCalRecHits(ESRE,ESSHE);
return;
}
| void GlobalRecHitsProducer::storeHCal | ( | PGlobalRecHit & | product | ) | [private] |
Definition at line 794 of file GlobalRecHitsProducer.cc.
References HBCalR, HBCalREC, HBCalSHE, HECalR, HECalREC, HECalSHE, HFCalR, HFCalREC, HFCalSHE, HOCalR, HOCalREC, HOCalSHE, i, PGlobalRecHit::putHBCalRecHits(), PGlobalRecHit::putHECalRecHits(), PGlobalRecHit::putHFCalRecHits(), PGlobalRecHit::putHOCalRecHits(), and verbosity.
Referenced by produce().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_storeHCal";
if (verbosity > 2) {
TString eventout("\n nHBRecHits = ");
eventout += HBCalREC.size();
for (unsigned int i = 0; i < HBCalREC.size(); ++i) {
eventout += "\n (REC, R, SHE) = (";
eventout += HBCalREC[i];
eventout += ", ";
eventout += HBCalR[i];
eventout += ", ";
eventout += HBCalSHE[i];
eventout += ")";
}
eventout += "\n nHERecHits = ";
eventout += HECalREC.size();
for (unsigned int i = 0; i < HECalREC.size(); ++i) {
eventout += "\n (REC, R, SHE) = (";
eventout += HECalREC[i];
eventout += ", ";
eventout += HECalR[i];
eventout += ", ";
eventout += HECalSHE[i];
eventout += ")";
}
eventout += "\n nHFRecHits = ";
eventout += HFCalREC.size();
for (unsigned int i = 0; i < HFCalREC.size(); ++i) {
eventout += "\n (REC, R, SHE) = (";
eventout += HFCalREC[i];
eventout += ", ";
eventout += HFCalR[i];
eventout += ", ";
eventout += HFCalSHE[i];
eventout += ")";
}
eventout += "\n nHORecHits = ";
eventout += HOCalREC.size();
for (unsigned int i = 0; i < HOCalREC.size(); ++i) {
eventout += "\n (REC, R, SHE) = (";
eventout += HOCalREC[i];
eventout += ", ";
eventout += HOCalR[i];
eventout += ", ";
eventout += HOCalSHE[i];
eventout += ")";
}
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
}
product.putHBCalRecHits(HBCalREC,HBCalR,HBCalSHE);
product.putHECalRecHits(HECalREC,HECalR,HECalSHE);
product.putHOCalRecHits(HOCalREC,HOCalR,HOCalSHE);
product.putHFCalRecHits(HFCalREC,HFCalR,HFCalSHE);
return;
}
| void GlobalRecHitsProducer::storeMuon | ( | PGlobalRecHit & | product | ) | [private] |
Definition at line 1873 of file GlobalRecHitsProducer.cc.
References CSCRHPERP, CSCRHPHI, CSCSHPHI, DTRHD, DTSHD, i, PGlobalRecHit::putCSCRecHits(), PGlobalRecHit::putDTRecHits(), PGlobalRecHit::putRPCRecHits(), RPCRHX, RPCSHX, and verbosity.
Referenced by produce().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_storeMuon";
if (verbosity > 2) {
// dt output
TString eventout("\n nDT = ");
eventout += DTRHD.size();
for (unsigned int i = 0; i < DTRHD.size(); ++i) {
eventout += "\n (RHD, SHD) = (";
eventout += DTRHD[i];
eventout += ", ";
eventout += DTSHD[i];
eventout += ")";
}
// CSC Strip
eventout += "\n nCSC = ";
eventout += CSCRHPHI.size();
for (unsigned int i = 0; i < CSCRHPHI.size(); ++i) {
eventout += "\n (rhphi, rhperp, shphi) = (";
eventout += CSCRHPHI[i];
eventout += ", ";
eventout += CSCRHPERP[i];
eventout += ", ";
eventout += CSCSHPHI[i];
eventout += ")";
}
// RPC
eventout += "\n nRPC = ";
eventout += RPCRHX.size();
for (unsigned int i = 0; i < RPCRHX.size(); ++i) {
eventout += "\n (rhx, shx) = (";
eventout += RPCRHX[i];
eventout += ", ";
eventout += RPCSHX[i];
eventout += ")";
}
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
}
product.putDTRecHits(DTRHD,DTSHD);
product.putCSCRecHits(CSCRHPHI,CSCRHPERP,CSCSHPHI);
product.putRPCRecHits(RPCRHX,RPCSHX);
return;
}
| void GlobalRecHitsProducer::storeTrk | ( | PGlobalRecHit & | product | ) | [private] |
Definition at line 1267 of file GlobalRecHitsProducer.cc.
References BRL1RX, BRL1RY, BRL1SX, BRL1SY, BRL2RX, BRL2RY, BRL2SX, BRL2SY, BRL3RX, BRL3RY, BRL3SX, BRL3SY, FWD1nRX, FWD1nRY, FWD1nSX, FWD1nSY, FWD1pRX, FWD1pRY, FWD1pSX, FWD1pSY, FWD2nRX, FWD2nRY, FWD2nSX, FWD2nSY, FWD2pRX, FWD2pRY, FWD2pSX, FWD2pSY, i, PGlobalRecHit::putBRL1RecHits(), PGlobalRecHit::putBRL2RecHits(), PGlobalRecHit::putBRL3RecHits(), PGlobalRecHit::putFWD1nRecHits(), PGlobalRecHit::putFWD1pRecHits(), PGlobalRecHit::putFWD2nRecHits(), PGlobalRecHit::putFWD2pRecHits(), PGlobalRecHit::putTECW1RecHits(), PGlobalRecHit::putTECW2RecHits(), PGlobalRecHit::putTECW3RecHits(), PGlobalRecHit::putTECW4RecHits(), PGlobalRecHit::putTECW5RecHits(), PGlobalRecHit::putTECW6RecHits(), PGlobalRecHit::putTECW7RecHits(), PGlobalRecHit::putTECW8RecHits(), PGlobalRecHit::putTIBL1RecHits(), PGlobalRecHit::putTIBL2RecHits(), PGlobalRecHit::putTIBL3RecHits(), PGlobalRecHit::putTIBL4RecHits(), PGlobalRecHit::putTIDW1RecHits(), PGlobalRecHit::putTIDW2RecHits(), PGlobalRecHit::putTIDW3RecHits(), PGlobalRecHit::putTOBL1RecHits(), PGlobalRecHit::putTOBL2RecHits(), PGlobalRecHit::putTOBL3RecHits(), PGlobalRecHit::putTOBL4RecHits(), TECW1RX, TECW1RY, TECW1SX, TECW1SY, TECW2RX, TECW2RY, TECW2SX, TECW2SY, TECW3RX, TECW3RY, TECW3SX, TECW3SY, TECW4RX, TECW4RY, TECW4SX, TECW4SY, TECW5RX, TECW5RY, TECW5SX, TECW5SY, TECW6RX, TECW6RY, TECW6SX, TECW6SY, TECW7RX, TECW7RY, TECW7SX, TECW7SY, TECW8RX, TECW8RY, TECW8SX, TECW8SY, TIBL1RX, TIBL1RY, TIBL1SX, TIBL1SY, TIBL2RX, TIBL2RY, TIBL2SX, TIBL2SY, TIBL3RX, TIBL3RY, TIBL3SX, TIBL3SY, TIBL4RX, TIBL4RY, TIBL4SX, TIBL4SY, TIDW1RX, TIDW1RY, TIDW1SX, TIDW1SY, TIDW2RX, TIDW2RY, TIDW2SX, TIDW2SY, TIDW3RX, TIDW3RY, TIDW3SX, TIDW3SY, TOBL1RX, TOBL1RY, TOBL1SX, TOBL1SY, TOBL2RX, TOBL2RY, TOBL2SX, TOBL2SY, TOBL3RX, TOBL3RY, TOBL3SX, TOBL3SY, TOBL4RX, TOBL4RY, TOBL4SX, TOBL4SY, and verbosity.
Referenced by produce().
{
std::string MsgLoggerCat = "GlobalRecHitsProducer_storeTrk";
if (verbosity > 2) {
// strip output
TString eventout("\n nTIBL1 = ");
eventout += TIBL1RX.size();
for (unsigned int i = 0; i < TIBL1RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TIBL1RX[i];
eventout += ", ";
eventout += TIBL1RY[i];
eventout += ", ";
eventout += TIBL1SX[i];
eventout += ", ";
eventout += TIBL1SY[i];
eventout += ")";
}
eventout += "\n nTIBL2 = ";
eventout += TIBL2RX.size();
for (unsigned int i = 0; i < TIBL2RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TIBL2RX[i];
eventout += ", ";
eventout += TIBL2RY[i];
eventout += ", ";
eventout += TIBL2SX[i];
eventout += ", ";
eventout += TIBL2SY[i];
eventout += ")";
}
eventout += "\n nTIBL3 = ";
eventout += TIBL3RX.size();
for (unsigned int i = 0; i < TIBL3RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TIBL3RX[i];
eventout += ", ";
eventout += TIBL3RY[i];
eventout += ", ";
eventout += TIBL3SX[i];
eventout += ", ";
eventout += TIBL3SY[i];
eventout += ")";
}
eventout += "\n nTIBL4 = ";
eventout += TIBL4RX.size();
for (unsigned int i = 0; i < TIBL4RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TIBL4RX[i];
eventout += ", ";
eventout += TIBL4RY[i];
eventout += ", ";
eventout += TIBL4SX[i];
eventout += ", ";
eventout += TIBL4SY[i];
eventout += ")";
}
eventout += "\n nTOBL1 = ";
eventout += TOBL1RX.size();
for (unsigned int i = 0; i < TOBL1RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TOBL1RX[i];
eventout += ", ";
eventout += TOBL1RY[i];
eventout += ", ";
eventout += TOBL1SX[i];
eventout += ", ";
eventout += TOBL1SY[i];
eventout += ")";
}
eventout += "\n nTOBL2 = ";
eventout += TOBL2RX.size();
for (unsigned int i = 0; i < TOBL2RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TOBL2RX[i];
eventout += ", ";
eventout += TOBL2RY[i];
eventout += ", ";
eventout += TOBL2SX[i];
eventout += ", ";
eventout += TOBL2SY[i];
eventout += ")";
}
eventout += "\n nTOBL3 = ";
eventout += TOBL3RX.size();
for (unsigned int i = 0; i < TOBL3RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TOBL3RX[i];
eventout += ", ";
eventout += TOBL3RY[i];
eventout += ", ";
eventout += TOBL3SX[i];
eventout += ", ";
eventout += TOBL3SY[i];
eventout += ")";
}
eventout += "\n nTOBL4 = ";
eventout += TOBL4RX.size();
for (unsigned int i = 0; i < TOBL4RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TOBL4RX[i];
eventout += ", ";
eventout += TOBL4RY[i];
eventout += ", ";
eventout += TOBL4SX[i];
eventout += ", ";
eventout += TOBL4SY[i];
eventout += ")";
}
eventout += "\n nTIDW1 = ";
eventout += TIDW1RX.size();
for (unsigned int i = 0; i < TIDW1RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TIDW1RX[i];
eventout += ", ";
eventout += TIDW1RY[i];
eventout += ", ";
eventout += TIDW1SX[i];
eventout += ", ";
eventout += TIDW1SY[i];
eventout += ")";
}
eventout += "\n nTIDW2 = ";
eventout += TIDW2RX.size();
for (unsigned int i = 0; i < TIDW2RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TIDW2RX[i];
eventout += ", ";
eventout += TIDW2RY[i];
eventout += ", ";
eventout += TIDW2SX[i];
eventout += ", ";
eventout += TIDW2SY[i];
eventout += ")";
}
eventout += "\n nTIDW3 = ";
eventout += TIDW3RX.size();
for (unsigned int i = 0; i < TIDW3RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TIDW3RX[i];
eventout += ", ";
eventout += TIDW3RY[i];
eventout += ", ";
eventout += TIDW3SX[i];
eventout += ", ";
eventout += TIDW3SY[i];
eventout += ")";
}
eventout += "\n nTECW1 = ";
eventout += TECW1RX.size();
for (unsigned int i = 0; i < TECW1RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TECW1RX[i];
eventout += ", ";
eventout += TECW1RY[i];
eventout += ", ";
eventout += TECW1SX[i];
eventout += ", ";
eventout += TECW1SY[i];
eventout += ")";
}
eventout += "\n nTECW2 = ";
eventout += TECW2RX.size();
for (unsigned int i = 0; i < TECW2RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TECW2RX[i];
eventout += ", ";
eventout += TECW2RY[i];
eventout += ", ";
eventout += TECW2SX[i];
eventout += ", ";
eventout += TECW2SY[i];
eventout += ")";
}
eventout += "\n nTECW3 = ";
eventout += TECW3RX.size();
for (unsigned int i = 0; i < TECW3RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TECW3RX[i];
eventout += ", ";
eventout += TECW3RY[i];
eventout += ", ";
eventout += TECW3SX[i];
eventout += ", ";
eventout += TECW3SY[i];
eventout += ")";
}
eventout += "\n nTECW4 = ";
eventout += TECW4RX.size();
for (unsigned int i = 0; i < TECW4RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TECW4RX[i];
eventout += ", ";
eventout += TECW4RY[i];
eventout += ", ";
eventout += TECW4SX[i];
eventout += ", ";
eventout += TECW4SY[i];
eventout += ")";
}
eventout += "\n nTECW5 = ";
eventout += TECW5RX.size();
for (unsigned int i = 0; i < TECW5RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TECW5RX[i];
eventout += ", ";
eventout += TECW5RY[i];
eventout += ", ";
eventout += TECW5SX[i];
eventout += ", ";
eventout += TECW5SY[i];
eventout += ")";
}
eventout += "\n nTECW6 = ";
eventout += TECW6RX.size();
for (unsigned int i = 0; i < TECW6RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TECW6RX[i];
eventout += ", ";
eventout += TECW6RY[i];
eventout += ", ";
eventout += TECW6SX[i];
eventout += ", ";
eventout += TECW6SY[i];
eventout += ")";
}
eventout += "\n nTECW7 = ";
eventout += TECW7RX.size();
for (unsigned int i = 0; i < TECW7RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TECW7RX[i];
eventout += ", ";
eventout += TECW7RY[i];
eventout += ", ";
eventout += TECW7SX[i];
eventout += ", ";
eventout += TECW7SY[i];
eventout += ")";
}
eventout += "\n nTECW8 = ";
eventout += TECW8RX.size();
for (unsigned int i = 0; i < TECW8RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += TECW8RX[i];
eventout += ", ";
eventout += TECW8RY[i];
eventout += ", ";
eventout += TECW8SX[i];
eventout += ", ";
eventout += TECW8SY[i];
eventout += ")";
}
// pixel output
eventout += "\n nBRL1 = ";
eventout += BRL1RX.size();
for (unsigned int i = 0; i < BRL1RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += BRL1RX[i];
eventout += ", ";
eventout += BRL1RY[i];
eventout += ", ";
eventout += BRL1SX[i];
eventout += ", ";
eventout += BRL1SY[i];
eventout += ")";
}
eventout += "\n nBRL2 = ";
eventout += BRL2RX.size();
for (unsigned int i = 0; i < BRL2RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += BRL2RX[i];
eventout += ", ";
eventout += BRL2RY[i];
eventout += ", ";
eventout += BRL2SX[i];
eventout += ", ";
eventout += BRL2SY[i];
eventout += ")";
}
eventout += "\n nBRL3 = ";
eventout += BRL3RX.size();
for (unsigned int i = 0; i < BRL3RX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += BRL3RX[i];
eventout += ", ";
eventout += BRL3RY[i];
eventout += ", ";
eventout += BRL3SX[i];
eventout += ", ";
eventout += BRL3SY[i];
eventout += ")";
}
eventout += "\n nFWD1p = ";
eventout += FWD1pRX.size();
for (unsigned int i = 0; i < FWD1pRX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += FWD1pRX[i];
eventout += ", ";
eventout += FWD1pRY[i];
eventout += ", ";
eventout += FWD1pSX[i];
eventout += ", ";
eventout += FWD1pSY[i];
eventout += ")";
}
eventout += "\n nFWD1n = ";
eventout += FWD1nRX.size();
for (unsigned int i = 0; i < FWD1nRX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += FWD1nRX[i];
eventout += ", ";
eventout += FWD1nRY[i];
eventout += ", ";
eventout += FWD1nSX[i];
eventout += ", ";
eventout += FWD1nSY[i];
eventout += ")";
}
eventout += "\n nFWD2p = ";
eventout += FWD2pRX.size();
for (unsigned int i = 0; i < FWD2pRX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += FWD2pRX[i];
eventout += ", ";
eventout += FWD2pRY[i];
eventout += ", ";
eventout += FWD2pSX[i];
eventout += ", ";
eventout += FWD2pSY[i];
eventout += ")";
}
eventout += "\n nFWD2p = ";
eventout += FWD2nRX.size();
for (unsigned int i = 0; i < FWD2nRX.size(); ++i) {
eventout += "\n (RX, RY, SX, SY) = (";
eventout += FWD2nRX[i];
eventout += ", ";
eventout += FWD2nRY[i];
eventout += ", ";
eventout += FWD2nSX[i];
eventout += ", ";
eventout += FWD2nSY[i];
eventout += ")";
}
edm::LogInfo(MsgLoggerCat) << eventout << "\n";
}
// strip output
product.putTIBL1RecHits(TIBL1RX,TIBL1RY,TIBL1SX,TIBL1SY);
product.putTIBL2RecHits(TIBL2RX,TIBL2RY,TIBL2SX,TIBL2SY);
product.putTIBL3RecHits(TIBL3RX,TIBL3RY,TIBL3SX,TIBL3SY);
product.putTIBL4RecHits(TIBL4RX,TIBL4RY,TIBL4SX,TIBL4SY);
product.putTOBL1RecHits(TOBL1RX,TOBL1RY,TOBL1SX,TOBL1SY);
product.putTOBL2RecHits(TOBL2RX,TOBL2RY,TOBL2SX,TOBL2SY);
product.putTOBL3RecHits(TOBL3RX,TOBL3RY,TOBL3SX,TOBL3SY);
product.putTOBL4RecHits(TOBL4RX,TOBL4RY,TOBL4SX,TOBL4SY);
product.putTIDW1RecHits(TIDW1RX,TIDW1RY,TIDW1SX,TIDW1SY);
product.putTIDW2RecHits(TIDW2RX,TIDW2RY,TIDW2SX,TIDW2SY);
product.putTIDW3RecHits(TIDW3RX,TIDW3RY,TIDW3SX,TIDW3SY);
product.putTECW1RecHits(TECW1RX,TECW1RY,TECW1SX,TECW1SY);
product.putTECW2RecHits(TECW2RX,TECW2RY,TECW2SX,TECW2SY);
product.putTECW3RecHits(TECW3RX,TECW3RY,TECW3SX,TECW3SY);
product.putTECW4RecHits(TECW4RX,TECW4RY,TECW4SX,TECW4SY);
product.putTECW5RecHits(TECW5RX,TECW5RY,TECW5SX,TECW5SY);
product.putTECW6RecHits(TECW6RX,TECW6RY,TECW6SX,TECW6SY);
product.putTECW7RecHits(TECW7RX,TECW7RY,TECW7SX,TECW7SY);
product.putTECW8RecHits(TECW8RX,TECW8RY,TECW8SX,TECW8SY);
// pixel output
product.putBRL1RecHits(BRL1RX,BRL1RY,BRL1SX,BRL1SY);
product.putBRL2RecHits(BRL2RX,BRL2RY,BRL2SX,BRL2SY);
product.putBRL3RecHits(BRL3RX,BRL3RY,BRL3SX,BRL3SY);
product.putFWD1pRecHits(FWD1pRX,FWD1pRY,FWD1pSX,FWD1pSY);
product.putFWD1nRecHits(FWD1nRX,FWD1nRY,FWD1nSX,FWD1nSY);
product.putFWD2pRecHits(FWD2pRX,FWD2pRY,FWD2pSX,FWD2pSY);
product.putFWD2nRecHits(FWD2nRX,FWD2nRY,FWD2nSX,FWD2nSY);
return;
}
FloatVector GlobalRecHitsProducer::BRL1RX [private] |
Definition at line 274 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL1RY [private] |
Definition at line 275 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL1SX [private] |
Definition at line 276 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL1SY [private] |
Definition at line 277 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL2RX [private] |
Definition at line 274 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL2RY [private] |
Definition at line 275 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL2SX [private] |
Definition at line 276 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL2SY [private] |
Definition at line 277 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL3RX [private] |
Definition at line 274 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL3RY [private] |
Definition at line 275 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL3SX [private] |
Definition at line 276 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::BRL3SY [private] |
Definition at line 277 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
Definition at line 270 of file GlobalRecHitsProducer.h.
Referenced by fillTrk(), and GlobalRecHitsProducer().
unsigned int GlobalRecHitsProducer::count [private] |
Definition at line 346 of file GlobalRecHitsProducer.h.
FloatVector GlobalRecHitsProducer::CSCRHPERP [private] |
Definition at line 328 of file GlobalRecHitsProducer.h.
Referenced by clear(), plotResolution(), and storeMuon().
FloatVector GlobalRecHitsProducer::CSCRHPHI [private] |
Definition at line 327 of file GlobalRecHitsProducer.h.
Referenced by clear(), plotResolution(), and storeMuon().
FloatVector GlobalRecHitsProducer::CSCSHPHI [private] |
Definition at line 329 of file GlobalRecHitsProducer.h.
Referenced by clear(), plotResolution(), and storeMuon().
FloatVector GlobalRecHitsProducer::DTRHD [private] |
Definition at line 289 of file GlobalRecHitsProducer.h.
Referenced by clear(), compute(), and storeMuon().
FloatVector GlobalRecHitsProducer::DTSHD [private] |
Definition at line 290 of file GlobalRecHitsProducer.h.
Referenced by clear(), compute(), and storeMuon().
FloatVector GlobalRecHitsProducer::EBRE [private] |
Definition at line 201 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillECal(), and storeECal().
FloatVector GlobalRecHitsProducer::EBSHE [private] |
Definition at line 202 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillECal(), and storeECal().
Definition at line 210 of file GlobalRecHitsProducer.h.
Referenced by fillECal(), and GlobalRecHitsProducer().
Definition at line 212 of file GlobalRecHitsProducer.h.
Referenced by fillECal(), and GlobalRecHitsProducer().
Definition at line 214 of file GlobalRecHitsProducer.h.
Referenced by fillECal(), and GlobalRecHitsProducer().
Definition at line 211 of file GlobalRecHitsProducer.h.
Referenced by fillECal(), and GlobalRecHitsProducer().
Definition at line 213 of file GlobalRecHitsProducer.h.
Referenced by fillECal(), and GlobalRecHitsProducer().
FloatVector GlobalRecHitsProducer::EERE [private] |
Definition at line 204 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillECal(), and storeECal().
FloatVector GlobalRecHitsProducer::EESHE [private] |
Definition at line 205 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillECal(), and storeECal().
FloatVector GlobalRecHitsProducer::ESRE [private] |
Definition at line 207 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillECal(), and storeECal().
FloatVector GlobalRecHitsProducer::ESSHE [private] |
Definition at line 208 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillECal(), and storeECal().
std::string GlobalRecHitsProducer::fName [private] |
Definition at line 191 of file GlobalRecHitsProducer.h.
Referenced by GlobalRecHitsProducer().
int GlobalRecHitsProducer::frequency [private] |
Definition at line 193 of file GlobalRecHitsProducer.h.
Referenced by GlobalRecHitsProducer(), and produce().
FloatVector GlobalRecHitsProducer::FWD1nRX [private] |
Definition at line 279 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD1nRY [private] |
Definition at line 280 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD1nSX [private] |
Definition at line 281 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD1nSY [private] |
Definition at line 282 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD1pRX [private] |
Definition at line 279 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD1pRY [private] |
Definition at line 280 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD1pSX [private] |
Definition at line 281 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD1pSY [private] |
Definition at line 282 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD2nRX [private] |
Definition at line 279 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD2nRY [private] |
Definition at line 280 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD2nSX [private] |
Definition at line 281 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD2nSY [private] |
Definition at line 282 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD2pRX [private] |
Definition at line 279 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD2pRY [private] |
Definition at line 280 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD2pSX [private] |
Definition at line 281 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::FWD2pSY [private] |
Definition at line 282 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
bool GlobalRecHitsProducer::getAllProvenances [private] |
Definition at line 195 of file GlobalRecHitsProducer.h.
Referenced by GlobalRecHitsProducer(), and produce().
FloatVector GlobalRecHitsProducer::HBCalR [private] |
Definition at line 219 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HBCalREC [private] |
Definition at line 218 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HBCalSHE [private] |
Definition at line 220 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
edm::InputTag GlobalRecHitsProducer::HCalSrc_ [private] |
Definition at line 234 of file GlobalRecHitsProducer.h.
Referenced by fillHCal(), and GlobalRecHitsProducer().
FloatVector GlobalRecHitsProducer::HECalR [private] |
Definition at line 223 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HECalREC [private] |
Definition at line 222 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HECalSHE [private] |
Definition at line 224 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HFCalR [private] |
Definition at line 231 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HFCalREC [private] |
Definition at line 230 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HFCalSHE [private] |
Definition at line 232 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HOCalR [private] |
Definition at line 227 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HOCalREC [private] |
Definition at line 226 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
FloatVector GlobalRecHitsProducer::HOCalSHE [private] |
Definition at line 228 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillHCal(), and storeHCal().
std::string GlobalRecHitsProducer::label [private] |
Definition at line 194 of file GlobalRecHitsProducer.h.
Referenced by GlobalRecHitsProducer(), and produce().
std::vector<PSimHit> GlobalRecHitsProducer::matched [private] |
Definition at line 265 of file GlobalRecHitsProducer.h.
Referenced by fillTrk().
Definition at line 331 of file GlobalRecHitsProducer.h.
Referenced by fillMuon(), and GlobalRecHitsProducer().
Definition at line 293 of file GlobalRecHitsProducer.h.
Referenced by fillMuon(), and GlobalRecHitsProducer().
edm::InputTag GlobalRecHitsProducer::MuDTSrc_ [private] |
Definition at line 292 of file GlobalRecHitsProducer.h.
Referenced by fillMuon(), and GlobalRecHitsProducer().
Definition at line 343 of file GlobalRecHitsProducer.h.
Referenced by fillMuon(), and GlobalRecHitsProducer().
Definition at line 342 of file GlobalRecHitsProducer.h.
Referenced by fillMuon(), and GlobalRecHitsProducer().
bool GlobalRecHitsProducer::printProvenanceInfo [private] |
Definition at line 196 of file GlobalRecHitsProducer.h.
Referenced by GlobalRecHitsProducer(), and produce().
FloatVector GlobalRecHitsProducer::RPCRHX [private] |
Definition at line 339 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillMuon(), and storeMuon().
FloatVector GlobalRecHitsProducer::RPCSHX [private] |
Definition at line 340 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillMuon(), and storeMuon().
Definition at line 284 of file GlobalRecHitsProducer.h.
Referenced by fillTrk(), and GlobalRecHitsProducer().
Definition at line 263 of file GlobalRecHitsProducer.h.
Referenced by fillTrk(), and GlobalRecHitsProducer().
FloatVector GlobalRecHitsProducer::TECW1RX [private] |
Definition at line 254 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW1RY [private] |
Definition at line 256 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW1SX [private] |
Definition at line 258 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW1SY [private] |
Definition at line 260 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW2RX [private] |
Definition at line 254 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW2RY [private] |
Definition at line 256 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW2SX [private] |
Definition at line 258 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW2SY [private] |
Definition at line 260 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW3RX [private] |
Definition at line 254 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW3RY [private] |
Definition at line 256 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW3SX [private] |
Definition at line 258 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW3SY [private] |
Definition at line 260 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW4RX [private] |
Definition at line 254 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW4RY [private] |
Definition at line 256 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW4SX [private] |
Definition at line 258 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW4SY [private] |
Definition at line 260 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW5RX [private] |
Definition at line 254 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW5RY [private] |
Definition at line 256 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW5SX [private] |
Definition at line 258 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW5SY [private] |
Definition at line 260 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW6RX [private] |
Definition at line 254 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW6RY [private] |
Definition at line 256 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW6SX [private] |
Definition at line 258 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW6SY [private] |
Definition at line 260 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW7RX [private] |
Definition at line 254 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW7RY [private] |
Definition at line 256 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW7SX [private] |
Definition at line 258 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW7SY [private] |
Definition at line 260 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW8RX [private] |
Definition at line 254 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW8RY [private] |
Definition at line 256 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW8SX [private] |
Definition at line 258 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TECW8SY [private] |
Definition at line 260 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
std::map<int, edm::PSimHitContainer> GlobalRecHitsProducer::theMap [private] |
Definition at line 333 of file GlobalRecHitsProducer.h.
Referenced by fillMuon().
FloatVector GlobalRecHitsProducer::TIBL1RX [private] |
Definition at line 239 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL1RY [private] |
Definition at line 240 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL1SX [private] |
Definition at line 241 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL1SY [private] |
Definition at line 242 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL2RX [private] |
Definition at line 239 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL2RY [private] |
Definition at line 240 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL2SX [private] |
Definition at line 241 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL2SY [private] |
Definition at line 242 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL3RX [private] |
Definition at line 239 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL3RY [private] |
Definition at line 240 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL3SX [private] |
Definition at line 241 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL3SY [private] |
Definition at line 242 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL4RX [private] |
Definition at line 239 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL4RY [private] |
Definition at line 240 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL4SX [private] |
Definition at line 241 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIBL4SY [private] |
Definition at line 242 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW1RX [private] |
Definition at line 249 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW1RY [private] |
Definition at line 250 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW1SX [private] |
Definition at line 251 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW1SY [private] |
Definition at line 252 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW2RX [private] |
Definition at line 249 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW2RY [private] |
Definition at line 250 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW2SX [private] |
Definition at line 251 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW2SY [private] |
Definition at line 252 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW3RX [private] |
Definition at line 249 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW3RY [private] |
Definition at line 250 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW3SX [private] |
Definition at line 251 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TIDW3SY [private] |
Definition at line 252 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL1RX [private] |
Definition at line 244 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL1RY [private] |
Definition at line 245 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL1SX [private] |
Definition at line 246 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL1SY [private] |
Definition at line 247 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL2RX [private] |
Definition at line 244 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL2RY [private] |
Definition at line 245 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL2SX [private] |
Definition at line 246 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL2SY [private] |
Definition at line 247 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL3RX [private] |
Definition at line 244 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL3RY [private] |
Definition at line 245 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL3SX [private] |
Definition at line 246 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL3SY [private] |
Definition at line 247 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL4RX [private] |
Definition at line 244 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL4RY [private] |
Definition at line 245 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL4SX [private] |
Definition at line 246 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
FloatVector GlobalRecHitsProducer::TOBL4SY [private] |
Definition at line 247 of file GlobalRecHitsProducer.h.
Referenced by clear(), fillTrk(), and storeTrk().
int GlobalRecHitsProducer::verbosity [private] |
Definition at line 192 of file GlobalRecHitsProducer.h.
Referenced by clear(), endJob(), fillECal(), fillHCal(), fillMuon(), fillTrk(), GlobalRecHitsProducer(), produce(), storeECal(), storeHCal(), storeMuon(), and storeTrk().
1.7.3