|
|
|
#include <RPCDigiValid.h>
Public Member Functions | |
| RPCDigiValid (const edm::ParameterSet &ps) | |
| ~RPCDigiValid () | |
Protected Member Functions | |
| void | analyze (const edm::Event &e, const edm::EventSetup &c) |
| void | beginJob () |
| void | endJob (void) |
Private Attributes | |
| MonitorElement * | BxDist |
| DQMStore * | dbe_ |
| std::string | digiLabel |
| std::string | outputFile_ |
| MonitorElement * | Res |
| MonitorElement * | ResDmin1 |
| MonitorElement * | ResDmin2 |
| MonitorElement * | ResDmin3 |
| MonitorElement * | ResDplu1 |
| MonitorElement * | ResDplu2 |
| MonitorElement * | ResDplu3 |
| MonitorElement * | ResWmin1 |
| MonitorElement * | ResWmin2 |
| MonitorElement * | ResWplu1 |
| MonitorElement * | ResWplu2 |
| MonitorElement * | ResWzer0 |
| MonitorElement * | rzview |
| MonitorElement * | StripProf |
| MonitorElement * | xyvDmin1 |
| MonitorElement * | xyvDmin2 |
| MonitorElement * | xyvDmin3 |
| MonitorElement * | xyvDplu1 |
| MonitorElement * | xyvDplu2 |
| MonitorElement * | xyvDplu3 |
| MonitorElement * | xyview |
Definition at line 24 of file RPCDigiValid.h.
| RPCDigiValid::RPCDigiValid | ( | const edm::ParameterSet & | ps | ) |
Definition at line 19 of file RPCDigiValid.cc.
References BxDist, dbe_, digiLabel, edm::ParameterSet::getUntrackedParameter(), cmsCodeRules::cppFunctionSkipper::operator, outputFile_, Res, ResDmin1, ResDmin2, ResDmin3, ResDplu1, ResDplu2, ResDplu3, ResWmin1, ResWmin2, ResWplu1, ResWplu2, ResWzer0, rzview, StripProf, xyvDmin1, xyvDmin2, xyvDmin3, xyvDplu1, xyvDplu2, xyvDplu3, and xyview.
:dbe_(0){ digiLabel = ps.getUntrackedParameter<std::string>("digiLabel"); outputFile_ = ps.getUntrackedParameter<string>("outputFile", "rpcDigiValidPlots.root"); dbe_ = Service<DQMStore>().operator->(); if ( dbe_ ) { dbe_->setCurrentFolder("RPCDigisV/RPCDigis"); xyview = dbe_->book2D("X_Vs_Y_View","X_Vs_Y_View",1000, -760., 760., 1000, -760., 760.); xyvDplu1 = dbe_->book2D("Dplu1_XvsY","Dplu1_XvsY",1000, -760., 760., 1000, -760., 760.); xyvDplu2 = dbe_->book2D("Dplu2_XvsY","Dplu2_XvsY",1000, -760., 760., 1000, -760., 760.); xyvDplu3 = dbe_->book2D("Dplu3_XvsY","Dplu3_XvsY",1000, -760., 760., 1000, -760., 760.); xyvDmin1 = dbe_->book2D("Dmin1_XvsY","Dmin1_XvsY",1000, -760., 760., 1000, -760., 760.); xyvDmin2 = dbe_->book2D("Dmin2_XvsY","Dmin2_XvsY",1000, -760., 760., 1000, -760., 760.); xyvDmin3 = dbe_->book2D("Dmin3_XvsY","Dmin3_XvsY",1000, -760., 760., 1000, -760., 760.); rzview = dbe_->book2D("R_Vs_Z_View","R_Vs_Z_View",1000, -1100., 1100.,1000,0., 800.); Res = dbe_->book1D("Digi_SimHit_difference", "Digi_SimHit_difference", 300, -8, 8); ResWmin2 = dbe_->book1D("W_Min2_Residuals", "W_Min2_Residuals", 300, -8, 8); ResWmin1 = dbe_->book1D("W_Min1_Residuals", "W_Min1_Residuals", 300, -8, 8); ResWzer0 = dbe_->book1D("W_Zer0_Residuals", "W_Zer0_Residuals", 300, -8, 8); ResWplu1 = dbe_->book1D("W_Plu1_Residuals", "W_Plu1_Residuals", 300, -8, 8); ResWplu2 = dbe_->book1D("W_Plu2_Residuals", "W_Plu2_Residuals", 300, -8, 8); BxDist = dbe_->book1D("Bunch_Crossing", "Bunch_Crossing", 20, -10., 10.); StripProf = dbe_->book1D("Strip_Profile", "Strip_Profile", 100, 0, 100); ResDmin1 = dbe_->book1D("Disk_Min1_Residuals", "Disk_Min1_Residuals", 300, -8, 8); ResDmin2 = dbe_->book1D("Disk_Min2_Residuals", "Disk_Min2_Residuals", 300, -8, 8); ResDmin3 = dbe_->book1D("Disk_Min3_Residuals", "Disk_Min3_Residuals", 300, -8, 8); ResDplu1 = dbe_->book1D("Disk_Plu1_Residuals", "Disk_Plu1_Residuals", 300, -8, 8); ResDplu2 = dbe_->book1D("Disk_Plu2_Residuals", "Disk_Plu2_Residuals", 300, -8, 8); ResDplu3 = dbe_->book1D("Disk_Plu3_Residuals", "Disk_Plu3_Residuals", 300, -8, 8); } }
| RPCDigiValid::~RPCDigiValid | ( | ) |
Definition at line 58 of file RPCDigiValid.cc.
{}
| void RPCDigiValid::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [protected, virtual] |
Implements edm::EDAnalyzer.
Definition at line 66 of file RPCDigiValid.cc.
References BxDist, RPCRoll::centreOfStrip(), digiLabel, MonitorElement::Fill(), edm::EventSetup::get(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), RPCDetId::region(), Res, ResDmin1, ResDmin2, ResDmin3, ResDplu1, ResDplu2, ResDplu3, ResWmin1, ResWmin2, ResWplu1, ResWplu2, ResWzer0, RPCDetId::ring(), RPCDetId, rzview, RPCDetId::station(), StripProf, GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), x, xyvDmin1, xyvDmin2, xyvDmin3, xyvDplu1, xyvDplu2, xyvDplu3, xyview, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
// cout << endl <<"--- [RPCDigiQuality] Analysing Event: #Run: " << event.id().run()
// << " #Event: " << event.id().event() << endl;
// Get the RPC Geometry
edm::ESHandle<RPCGeometry> rpcGeom;
eventSetup.get<MuonGeometryRecord>().get(rpcGeom);
edm::Handle<PSimHitContainer> simHit;
event.getByLabel("g4SimHits", "MuonRPCHits", simHit);
edm::Handle<RPCDigiCollection> rpcDigis;
event.getByLabel(digiLabel, rpcDigis);
// Loop on simhits
PSimHitContainer::const_iterator simIt;
//loop over Simhit
std::map<RPCDetId, std::vector<double> > allsims;
for (simIt = simHit->begin(); simIt != simHit->end(); simIt++) {
RPCDetId Rsid = (RPCDetId)(*simIt).detUnitId();
const RPCRoll* soll = dynamic_cast<const RPCRoll* >( rpcGeom->roll(Rsid));
int ptype = simIt->particleType();
// std::cout <<"This is a Simhit with Parent "<<ptype<<std::endl;
if (ptype == 13 || ptype == -13) {
std::vector<double> buff;
if (allsims.find(Rsid) != allsims.end() ){
buff= allsims[Rsid];
}
buff.push_back(simIt->localPosition().x());
allsims[Rsid]=buff;
}
GlobalPoint p=soll->toGlobal(simIt->localPosition());
/*
std::cout <<"Muon Position phi="<<p.phi()
<<" R="<<p.perp()
<<" z="<<p.z()<<std::endl;
*/
double sim_x = p.x();
double sim_y = p.y();
xyview->Fill(sim_x, sim_y);
if (Rsid.region() == 0 ){
}
else if (Rsid.region() == (+1)){
if (Rsid.station() == 1)
xyvDplu1->Fill(sim_x, sim_y);
else if (Rsid.station() == 2)
xyvDplu2->Fill(sim_x, sim_y);
else if (Rsid.station() == 3)
xyvDplu3->Fill(sim_x, sim_y);
}
else if (Rsid.region() == (-1)){
if (Rsid.station() == 1)
xyvDmin1->Fill(sim_x, sim_y);
else if (Rsid.station() == 2)
xyvDmin2->Fill(sim_x, sim_y);
else if (Rsid.station() == 3)
xyvDmin3->Fill(sim_x, sim_y);
}
// xyview->Fill(p.x(),p.y());
rzview->Fill(p.z(),p.perp());
}
//loop over Digis
RPCDigiCollection::DigiRangeIterator detUnitIt;
for (detUnitIt=rpcDigis->begin(); detUnitIt!=rpcDigis->end();++detUnitIt){
const RPCDetId Rsid = (*detUnitIt).first;
const RPCRoll* roll = dynamic_cast<const RPCRoll* >( rpcGeom->roll(Rsid));
const RPCDigiCollection::Range& range = (*detUnitIt).second;
std::vector<double> sims;
if (allsims.find(Rsid) != allsims.end() ){
sims = allsims[Rsid];
}
int ndigi=0;
for (RPCDigiCollection::const_iterator digiIt = range.first;
digiIt != range.second; ++digiIt){
StripProf->Fill(digiIt->strip());
BxDist->Fill(digiIt->bx());
ndigi++;
}
// std::cout<<" Number of Digi "<<ndigi<<" for "<<Rsid<<std::endl;
if (sims.size() == 1 && ndigi == 1){
double dis = roll->centreOfStrip(range.first->strip()).x()-sims[0];
Res->Fill(dis);
if (Rsid.region() == 0 ){
if (Rsid.ring() == -2)
ResWmin2->Fill(dis);
else if (Rsid.ring() == -1)
ResWmin1->Fill(dis);
else if (Rsid.ring() == 0)
ResWzer0->Fill(dis);
else if (Rsid.ring() == 1)
ResWplu1->Fill(dis);
else if (Rsid.ring() == 2)
ResWplu2->Fill(dis);
}
else if (Rsid.region() == (+1)){
if (Rsid.station() == 1)
ResDplu1->Fill(dis);
else if (Rsid.station() == 2)
ResDplu2->Fill(dis);
else if (Rsid.station() == 3)
ResDplu3->Fill(dis);
}
else if (Rsid.region() == (-1)){
if (Rsid.station() == 1)
ResDmin1->Fill(dis);
else if (Rsid.station() == 2)
ResDmin2->Fill(dis);
else if (Rsid.station() == 3)
ResDmin3->Fill(dis);
}
}
}
}
| void RPCDigiValid::beginJob | ( | void | ) | [protected, virtual] |
| void RPCDigiValid::endJob | ( | void | ) | [protected, virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 62 of file RPCDigiValid.cc.
References dbe_, outputFile_, and DQMStore::save().
{
if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
}
MonitorElement* RPCDigiValid::BxDist [private] |
Definition at line 46 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
DQMStore* RPCDigiValid::dbe_ [private] |
Definition at line 66 of file RPCDigiValid.h.
Referenced by endJob(), and RPCDigiValid().
std::string RPCDigiValid::digiLabel [private] |
Definition at line 68 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
std::string RPCDigiValid::outputFile_ [private] |
Definition at line 67 of file RPCDigiValid.h.
Referenced by endJob(), and RPCDigiValid().
MonitorElement* RPCDigiValid::Res [private] |
Definition at line 40 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResDmin1 [private] |
Definition at line 50 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResDmin2 [private] |
Definition at line 51 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResDmin3 [private] |
Definition at line 52 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResDplu1 [private] |
Definition at line 53 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResDplu2 [private] |
Definition at line 54 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResDplu3 [private] |
Definition at line 55 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResWmin1 [private] |
Definition at line 42 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResWmin2 [private] |
Definition at line 41 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResWplu1 [private] |
Definition at line 44 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResWplu2 [private] |
Definition at line 45 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::ResWzer0 [private] |
Definition at line 43 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::rzview [private] |
Definition at line 39 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::StripProf [private] |
Definition at line 47 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::xyvDmin1 [private] |
Definition at line 62 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::xyvDmin2 [private] |
Definition at line 63 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::xyvDmin3 [private] |
Definition at line 64 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::xyvDplu1 [private] |
Definition at line 58 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::xyvDplu2 [private] |
Definition at line 59 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::xyvDplu3 [private] |
Definition at line 60 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
MonitorElement* RPCDigiValid::xyview [private] |
Definition at line 38 of file RPCDigiValid.h.
Referenced by analyze(), and RPCDigiValid().
1.7.3