|
|
|
#include <RPCEfficiency.h>
Class for RPC Monitoring using RPCDigi and DT and CSC Segments.
Definition at line 107 of file RPCEfficiency.h.
| RPCEfficiency::RPCEfficiency | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 37 of file RPCEfficiency.cc.
References gather_cfg::cout, cscSegments, dbe, debug, dt4DSegments, dupli, EffRootFileName, EffSaveRootFile, folderPath, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hGlobalResClu1La, hGlobalResClu1R2A, hGlobalResClu1R2B, hGlobalResClu1R2C, hGlobalResClu1R3A, hGlobalResClu1R3B, hGlobalResClu1R3C, hGlobalResClu2La, hGlobalResClu2R2A, hGlobalResClu2R2B, hGlobalResClu2R2C, hGlobalResClu2R3A, hGlobalResClu2R3B, hGlobalResClu2R3C, hGlobalResClu3La, hGlobalResClu3R2A, hGlobalResClu3R2B, hGlobalResClu3R2C, hGlobalResClu3R3A, hGlobalResClu3R3B, hGlobalResClu3R3C, inclcsc, incldt, incldtMB4, MaxD, MaxDrb4, MinCosAng, cmsCodeRules::cppFunctionSkipper::operator, rangestrips, rangestripsRB4, RPCRecHitLabel_, MonitorElement::setBinLabel(), and statistics.
{
incldt=iConfig.getUntrackedParameter<bool>("incldt",true);
incldtMB4=iConfig.getUntrackedParameter<bool>("incldtMB4",true);
inclcsc=iConfig.getUntrackedParameter<bool>("inclcsc",true);
debug=iConfig.getUntrackedParameter<bool>("debug",false);
rangestrips = iConfig.getUntrackedParameter<double>("rangestrips",4.);
rangestripsRB4=iConfig.getUntrackedParameter<double>("rangestripsRB4",4.);
dupli = iConfig.getUntrackedParameter<int>("DuplicationCorrection",2);
MinCosAng=iConfig.getUntrackedParameter<double>("MinCosAng",0.96);
MaxD=iConfig.getUntrackedParameter<double>("MaxD",80.);
MaxDrb4=iConfig.getUntrackedParameter<double>("MaxDrb4",150.);
// muonRPCDigis=iConfig.getUntrackedParameter<std::string>("muonRPCDigis","muonRPCDigis");
cscSegments=iConfig.getParameter<edm::InputTag>("cscSegments");
dt4DSegments=iConfig.getParameter<edm::InputTag>("dt4DSegments");
RPCRecHitLabel_ = iConfig.getParameter<edm::InputTag>("RecHitLabel");
folderPath=iConfig.getUntrackedParameter<std::string>("folderPath","RPC/RPCEfficiency/");
EffSaveRootFile = iConfig.getUntrackedParameter<bool>("EffSaveRootFile", false);
EffRootFileName = iConfig.getUntrackedParameter<std::string>("EffRootFileName", "RPCEfficiency.root");
//Interface
dbe = edm::Service<DQMStore>().operator->();
std::string folder;
dbe->setCurrentFolder(folderPath);
statistics = dbe->book1D("Statistics","All Statistics",33,0.5,33.5);
statistics->setBinLabel(1,"Events ",1);
statistics->setBinLabel(2,"Events with DT seg",1);
statistics->setBinLabel(3,"1 DT seg",1);
statistics->setBinLabel(4,"2 DT seg",1);
statistics->setBinLabel(5,"3 DT seg",1);
statistics->setBinLabel(6,"4 DT seg",1);
statistics->setBinLabel(7,"5 DT seg",1);
statistics->setBinLabel(8,"6 DT seg",1);
statistics->setBinLabel(9,"7 DT seg",1);
statistics->setBinLabel(10,"8 DT seg",1);
statistics->setBinLabel(11,"9 DT seg",1);
statistics->setBinLabel(12,"10 DT seg",1);
statistics->setBinLabel(13,"11 DT seg",1);
statistics->setBinLabel(14,"12 DT seg",1);
statistics->setBinLabel(15,"13 DT seg",1);
statistics->setBinLabel(16,"14 DT seg",1);
statistics->setBinLabel(17,"15 DT seg",1);
statistics->setBinLabel(18,"Events with CSC seg",1);
statistics->setBinLabel(16+3,"1 CSC seg",1);
statistics->setBinLabel(16+4,"2 CSC seg",1);
statistics->setBinLabel(16+5,"3 CSC seg",1);
statistics->setBinLabel(16+6,"4 CSC seg",1);
statistics->setBinLabel(16+7,"5 CSC seg",1);
statistics->setBinLabel(16+8,"6 CSC seg",1);
statistics->setBinLabel(16+9,"7 CSC seg",1);
statistics->setBinLabel(16+10,"8 CSC seg",1);
statistics->setBinLabel(16+11,"9 CSC seg",1);
statistics->setBinLabel(16+12,"10 CSC seg",1);
statistics->setBinLabel(16+13,"11 CSC seg",1);
statistics->setBinLabel(16+14,"12 CSC seg",1);
statistics->setBinLabel(16+15,"13 CSC seg",1);
statistics->setBinLabel(16+16,"14 CSC seg",1);
statistics->setBinLabel(16+17,"15 CSC seg",1);
if(debug) std::cout<<"booking Global histograms with "<<folderPath<<std::endl;
folder = folderPath+"MuonSegEff/"+"Residuals/Barrel";
dbe->setCurrentFolder(folder);
//Barrel
std::stringstream histoName, histoTitle;
for (int layer = 1 ; layer<= 6 ;layer++){
histoName.str("");
histoTitle.str("");
histoName<<"GlobalResidualsClu1La"<<layer;
histoTitle<<"RPC Residuals Layer "<<layer<<" Cluster Size 1";
hGlobalResClu1La[layer-1] = dbe->book1D(histoName.str(), histoTitle.str(),101,-10.,10.);
histoName.str("");
histoTitle.str("");
histoName<<"GlobalResidualsClu2La"<<layer;
histoTitle<<"RPC Residuals Layer "<<layer<<" Cluster Size 2";
hGlobalResClu2La[layer-1] = dbe->book1D(histoName.str(), histoTitle.str(),101,-10.,10.);
histoName.str("");
histoTitle.str("");
histoName<<"GlobalResidualsClu3La"<<layer;
histoTitle<<"RPC Residuals Layer "<<layer<<" Cluster Size 3";
hGlobalResClu3La[layer-1] = dbe->book1D(histoName.str(), histoTitle.str(),101,-10.,10.);
}
if(debug) std::cout<<"Booking Residuals for EndCap"<<std::endl;
folder = folderPath+"MuonSegEff/Residuals/EndCap";
dbe->setCurrentFolder(folder);
//Endcap
hGlobalResClu1R3C = dbe->book1D("GlobalResidualsClu1R3C","RPC Residuals Ring 3 Roll C Cluster Size 1",101,-10.,10.);
hGlobalResClu1R3B = dbe->book1D("GlobalResidualsClu1R3B","RPC Residuals Ring 3 Roll B Cluster Size 1",101,-10.,10.);
hGlobalResClu1R3A = dbe->book1D("GlobalResidualsClu1R3A","RPC Residuals Ring 3 Roll A Cluster Size 1",101,-10.,10.);
hGlobalResClu1R2C = dbe->book1D("GlobalResidualsClu1R2C","RPC Residuals Ring 2 Roll C Cluster Size 1",101,-10.,10.);
hGlobalResClu1R2B = dbe->book1D("GlobalResidualsClu1R2B","RPC Residuals Ring 2 Roll B Cluster Size 1",101,-10.,10.);
hGlobalResClu1R2A = dbe->book1D("GlobalResidualsClu1R2A","RPC Residuals Ring 2 Roll A Cluster Size 1",101,-10.,10.);
hGlobalResClu2R3C = dbe->book1D("GlobalResidualsClu2R3C","RPC Residuals Ring 3 Roll C Cluster Size 2",101,-10.,10.);
hGlobalResClu2R3B = dbe->book1D("GlobalResidualsClu2R3B","RPC Residuals Ring 3 Roll B Cluster Size 2",101,-10.,10.);
hGlobalResClu2R3A = dbe->book1D("GlobalResidualsClu2R3A","RPC Residuals Ring 3 Roll A Cluster Size 2",101,-10.,10.);
hGlobalResClu2R2C = dbe->book1D("GlobalResidualsClu2R2C","RPC Residuals Ring 2 Roll C Cluster Size 2",101,-10.,10.);
hGlobalResClu2R2B = dbe->book1D("GlobalResidualsClu2R2B","RPC Residuals Ring 2 Roll B Cluster Size 2",101,-10.,10.);
hGlobalResClu2R2A = dbe->book1D("GlobalResidualsClu2R2A","RPC Residuals Ring 2 Roll A Cluster Size 2",101,-10.,10.);
hGlobalResClu3R3C = dbe->book1D("GlobalResidualsClu3R3C","RPC Residuals Ring 3 Roll C Cluster Size 3",101,-10.,10.);
hGlobalResClu3R3B = dbe->book1D("GlobalResidualsClu3R3B","RPC Residuals Ring 3 Roll B Cluster Size 3",101,-10.,10.);
hGlobalResClu3R3A = dbe->book1D("GlobalResidualsClu3R3A","RPC Residuals Ring 3 Roll A Cluster Size 3",101,-10.,10.);
hGlobalResClu3R2C = dbe->book1D("GlobalResidualsClu3R2C","RPC Residuals Ring 2 Roll C Cluster Size 3",101,-10.,10.);
hGlobalResClu3R2B = dbe->book1D("GlobalResidualsClu3R2B","RPC Residuals Ring 2 Roll B Cluster Size 3",101,-10.,10.);
hGlobalResClu3R2A = dbe->book1D("GlobalResidualsClu3R2A","RPC Residuals Ring 2 Roll A Cluster Size 3",101,-10.,10.);
}
| RPCEfficiency::~RPCEfficiency | ( | ) |
Definition at line 262 of file RPCEfficiency.cc.
{}
| void RPCEfficiency::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 264 of file RPCEfficiency.cc.
References CSCDetId::chamber(), gather_cfg::cout, cscSegments, funct::D, debug, distsector_tmp(), dt4DSegments, dupli, CSCDetId::endcap(), MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByLabel(), hGlobalResClu1La, hGlobalResClu1R2A, hGlobalResClu1R2B, hGlobalResClu1R2C, hGlobalResClu1R3A, hGlobalResClu1R3B, hGlobalResClu1R3C, hGlobalResClu2La, hGlobalResClu2R2A, hGlobalResClu2R2B, hGlobalResClu2R2C, hGlobalResClu2R3A, hGlobalResClu2R3B, hGlobalResClu2R3C, hGlobalResClu3La, hGlobalResClu3R2A, hGlobalResClu3R2B, hGlobalResClu3R2C, hGlobalResClu3R3A, hGlobalResClu3R3B, hGlobalResClu3R3C, CSCChamber::id(), RPCRoll::id(), inclcsc, incldt, incldtMB4, edm::HandleBase::isValid(), RPCDetId::layer(), MaxD, MaxDrb4, meCollection, MinCosAng, RPCGeomServ::name(), RPCRoll::nstrips(), rangestrips, DetId::rawId(), RPCDetId::ring(), CSCDetId::ring(), RPCDetId::roll(), rollstoreCSC, rollstoreDT, RPCRecHitLabel_, DTChamberId::sector(), mathSSE::sqrt(), RPCDetId::station(), CSCDetId::station(), DTChamberId::station(), statistics, RPCRoll::strip(), GeomDet::surface(), Surface::toGlobal(), GeomDet::toGlobal(), GloballyPositioned< T >::toLocal(), GeomDet::toLocal(), RPCRoll::topology(), DTChamberId::wheel(), PV3DBase< T, PVType, FrameType >::x(), X, PV3DBase< T, PVType, FrameType >::y(), Gflash::Z, and PV3DBase< T, PVType, FrameType >::z().
{
edm::ESHandle<RPCGeometry> rpcGeo;
edm::ESHandle<DTGeometry> dtGeo;
edm::ESHandle<CSCGeometry> cscGeo;
iSetup.get<MuonGeometryRecord>().get(rpcGeo);
iSetup.get<MuonGeometryRecord>().get(dtGeo);
iSetup.get<MuonGeometryRecord>().get(cscGeo);
statistics->Fill(1);
std::stringstream meIdRPC, meIdDT, meIdCSC;
if(debug) std::cout <<"\t Getting the RPC RecHits"<<std::endl;
edm::Handle<RPCRecHitCollection> rpcHits;
iEvent.getByLabel(RPCRecHitLabel_,rpcHits);
if(!rpcHits.isValid()) return;
if(incldt){
if(debug) std::cout<<"\t Getting the DT Segments"<<std::endl;
edm::Handle<DTRecSegment4DCollection> all4DSegments;
iEvent.getByLabel(dt4DSegments, all4DSegments);
if(all4DSegments.isValid()){
if(all4DSegments->size()>0){
if(all4DSegments->size()<=16) statistics->Fill(2);
if(debug) std::cout<<"\t Number of DT Segments in this event = "<<all4DSegments->size()<<std::endl;
std::map<DTChamberId,int> DTSegmentCounter;
DTRecSegment4DCollection::const_iterator segment;
for (segment = all4DSegments->begin();segment!=all4DSegments->end(); ++segment){
DTSegmentCounter[segment->chamberId()]++;
}
statistics->Fill(all4DSegments->size()+2);
if(debug) std::cout<<"\t Loop over all the 4D Segments"<<std::endl;
for (segment = all4DSegments->begin(); segment != all4DSegments->end(); ++segment){
DTChamberId DTId = segment->chamberId();
if(DTSegmentCounter[DTId]==1 && DTId.station()!=4){
int dtWheel = DTId.wheel();
int dtStation = DTId.station();
int dtSector = DTId.sector();
LocalPoint segmentPosition= segment->localPosition();
LocalVector segmentDirection=segment->localDirection();
const GeomDet* gdet=dtGeo->idToDet(segment->geographicalId());
const BoundPlane & DTSurface = gdet->surface();
//check if the dimension of the segment is 4
if(segment->dimension()==4){
float Xo=segmentPosition.x();
float Yo=segmentPosition.y();
float Zo=segmentPosition.z();
float dx=segmentDirection.x();
float dy=segmentDirection.y();
float dz=segmentDirection.z();
std::set<RPCDetId> rollsForThisDT = rollstoreDT[DTStationIndex(0,dtWheel,dtSector,dtStation)];
if(debug) std::cout<<"DT \t \t Loop over all the rolls asociated to this DT"<<std::endl;
for (std::set<RPCDetId>::iterator iteraRoll = rollsForThisDT.begin();iteraRoll != rollsForThisDT.end(); iteraRoll++){
const RPCRoll* rollasociated = rpcGeo->roll(*iteraRoll);
RPCDetId rpcId = rollasociated->id();
const BoundPlane & RPCSurface = rollasociated->surface();
// RPCGeomServ rpcsrv(rpcId);
// std::string nameRoll = rpcsrv.name();
GlobalPoint CenterPointRollGlobal = RPCSurface.toGlobal(LocalPoint(0,0,0));
LocalPoint CenterRollinDTFrame = DTSurface.toLocal(CenterPointRollGlobal);
float D=CenterRollinDTFrame.z();
float X=Xo+dx*D/dz;
float Y=Yo+dy*D/dz;
float Z=D;
const RectangularStripTopology* top_= dynamic_cast<const RectangularStripTopology*> (&(rollasociated->topology()));
LocalPoint xmin = top_->localPosition(0.);
if(debug) std::cout<<"DT \t \t \t xmin of this Roll "<<xmin<<"cm"<<std::endl;
LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
if(debug) std::cout<<"DT \t \t \t xmax of this Roll "<<xmax<<"cm"<<std::endl;
float rsize = fabs( xmax.x()-xmin.x() );
if(debug) std::cout<<"DT \t \t \t Roll Size "<<rsize<<"cm"<<std::endl;
float stripl = top_->stripLength();
float stripw = top_->pitch();
float extrapolatedDistance = sqrt((X-Xo)*(X-Xo)+(Y-Yo)*(Y-Yo)+(Z-Zo)*(Z-Zo));
if(extrapolatedDistance<=MaxD){
GlobalPoint GlobalPointExtrapolated = DTSurface.toGlobal(LocalPoint(X,Y,Z));
LocalPoint PointExtrapolatedRPCFrame = RPCSurface.toLocal(GlobalPointExtrapolated);
if(fabs(PointExtrapolatedRPCFrame.z()) < 10. &&
fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
RPCDetId rollId = rollasociated->id();
RPCGeomServ rpcsrv(rollId);
std::string nameRoll = rpcsrv.name();
if(debug) std::cout<<"DT \t \t \t \t The RPCName is "<<nameRoll<<std::endl;
const float stripPredicted =
rollasociated->strip(LocalPoint(PointExtrapolatedRPCFrame.x(),PointExtrapolatedRPCFrame.y(),0.));
if(debug) std::cout<<"DT \t \t \t \t Candidate (from DT Segment) STRIP---> "<<stripPredicted<< std::endl;
//---- HISTOGRAM STRIP PREDICTED FROM DT ----
std::map<std::string, MonitorElement*> meMap=meCollection[rpcId.rawId()];
meIdDT.str("");
meIdDT<<"ExpectedOccupancyFromDT_"<<rollId.rawId();
meMap[meIdDT.str()]->Fill(stripPredicted);
//-----------------------------------------------------
//-------RecHitPart Just For Residual--------
int countRecHits = 0;
int cluSize = 0;
float minres = 3000.;
typedef std::pair<RPCRecHitCollection::const_iterator, RPCRecHitCollection::const_iterator> rangeRecHits;
rangeRecHits recHitCollection = rpcHits->get(rollasociated->id());
RPCRecHitCollection::const_iterator recHit;
for (recHit = recHitCollection.first; recHit != recHitCollection.second ; recHit++) {
countRecHits++;
// sprintf(meIdRPC,"BXDistribution_%d",rollasociated->id().rawId());
// meMap[meIdRPC]->Fill(recHit->BunchX());
LocalPoint recHitPos=recHit->localPosition();
float res=PointExtrapolatedRPCFrame.x()- recHitPos.x();
if(debug) std::cout<<"DT \t \t \t \t \t Found Rec Hit at "<<res<<"cm of the prediction."<<std::endl;
if(fabs(res)<fabs(minres)){
minres=res;
cluSize = recHit->clusterSize();
if(debug) std::cout<<"DT \t \t \t \t \t \t New Min Res "<<res<<"cm."<<std::endl;
}
}
if(countRecHits==0){
if(debug) std::cout <<"DT \t \t \t \t \t THIS ROLL DOESN'T HAVE ANY RECHIT"<<std::endl;
}else{
assert(minres!=3000);
if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
if(debug) std::cout<<"DT \t \t \t \t \t \t True!"<<std::endl;
// float cosal = dx/sqrt(dx*dx+dz*dz);
if(rollId.station()==1&&rollId.layer()==1) {
if(cluSize==1*dupli) {hGlobalResClu1La[0]->Fill(minres);}
else if(cluSize==2*dupli){ hGlobalResClu2La[0]->Fill(minres);}
else if(cluSize==3*dupli){ hGlobalResClu3La[0]->Fill(minres);}}
else if(rollId.station()==1&&rollId.layer()==2){
if(cluSize==1*dupli) {hGlobalResClu1La[1]->Fill(minres);}
else if(cluSize==2*dupli){ hGlobalResClu2La[1]->Fill(minres);}
else if(cluSize==3*dupli){ hGlobalResClu3La[1]->Fill(minres);}}
else if(rollId.station()==2&&rollId.layer()==1){
if(cluSize==1*dupli) {hGlobalResClu1La[2]->Fill(minres);}
else if(cluSize==2*dupli){ hGlobalResClu2La[2]->Fill(minres);}
else if(cluSize==3*dupli){ hGlobalResClu3La[2]->Fill(minres);}
}
else if(rollId.station()==2&&rollId.layer()==2){
if(cluSize==1*dupli) {hGlobalResClu1La[3]->Fill(minres);}
if(cluSize==2*dupli){ hGlobalResClu2La[3]->Fill(minres);}
else if(cluSize==3*dupli){ hGlobalResClu3La[3]->Fill(minres);}
}
else if(rollId.station()==3){
if(cluSize==1*dupli) {hGlobalResClu1La[4]->Fill(minres);}
else if(cluSize==2*dupli){ hGlobalResClu2La[4]->Fill(minres);}
else if(cluSize==3*dupli){ hGlobalResClu3La[4]->Fill(minres);}
}
meIdRPC.str("");
meIdRPC<<"RPCDataOccupancyFromDT_"<<rollId.rawId();
meMap[meIdRPC.str()]->Fill(stripPredicted);
}
}
}else{
if(debug) std::cout<<"DT \t \t \t \t No the prediction is outside of this roll"<<std::endl;
}//Condition for the right match
}else{
if(debug) std::cout<<"DT \t \t \t No, Exrtrapolation too long!, canceled"<<std::endl;
}//D so big
}//loop over all the rolls asociated
}//Is the segment 4D?
}else {
if(debug) std::cout<<"DT \t \t More than one segment in this chamber, or we are in Station 4"<<std::endl;
}
}
} else {
if(debug) std::cout<<"DT This Event doesn't have any DT4DDSegment"<<std::endl; //is ther more than 1 segment in this event?
}
}
}
if(incldtMB4){
edm::Handle<DTRecSegment4DCollection> all4DSegments;
iEvent.getByLabel(dt4DSegments, all4DSegments);
if(all4DSegments.isValid() && all4DSegments->size()>0){
std::map<DTChamberId,int> DTSegmentCounter;
DTRecSegment4DCollection::const_iterator segment;
for (segment = all4DSegments->begin();segment!=all4DSegments->end(); ++segment){
DTSegmentCounter[segment->chamberId()]++;
}
if(debug) std::cout<<"MB4 \t \t Loop Over all4DSegments"<<std::endl;
for (segment = all4DSegments->begin(); segment != all4DSegments->end(); ++segment){
DTChamberId DTId = segment->chamberId();
if(DTSegmentCounter[DTId] == 1 && DTId.station()==4){
int dtWheel = DTId.wheel();
int dtStation = DTId.station();
int dtSector = DTId.sector();
LocalPoint segmentPosition= segment->localPosition();
LocalVector segmentDirection=segment->localDirection();
//check if the dimension of the segment is 2 and the station is 4
if(segment->dimension()==2){
LocalVector segmentDirectionMB4=segmentDirection;
LocalPoint segmentPositionMB4=segmentPosition;
const BoundPlane& DTSurface4 = dtGeo->idToDet(DTId)->surface();
DTRecSegment4DCollection::const_iterator segMB3;
for(segMB3=all4DSegments->begin();segMB3!=all4DSegments->end();++segMB3){
DTChamberId dtid3 = segMB3->chamberId();
if(distsector_tmp(dtid3.sector(),DTId.sector())<=1
&& dtid3.station()==3
&& dtid3.wheel()==DTId.wheel()
&& DTSegmentCounter[dtid3] == 1
&& segMB3->dimension()==4){
const GeomDet* gdet3=dtGeo->idToDet(segMB3->geographicalId());
const BoundPlane & DTSurface3 = gdet3->surface();
LocalVector segmentDirectionMB3 = segMB3->localDirection();
GlobalPoint segmentPositionMB3inGlobal = DTSurface3.toGlobal(segMB3->localPosition());
GlobalVector segDirMB4inGlobalFrame=DTSurface4.toGlobal(segmentDirectionMB4);
GlobalVector segDirMB3inGlobalFrame=DTSurface3.toGlobal(segmentDirectionMB3);
float dx=segDirMB4inGlobalFrame.x();
float dy=segDirMB4inGlobalFrame.y();
// float dz=segDirMB4inGlobalFrame.z();
float dx3=segDirMB3inGlobalFrame.x();
float dy3=segDirMB3inGlobalFrame.y();
// float dz3=segDirMB3inGlobalFrame.z();
double cosAng=fabs(dx*dx3+dy*dy3/sqrt((dx3*dx3+dy3*dy3)*(dx*dx+dy*dy)));
if(cosAng>MinCosAng){
if(dtSector==13){
dtSector=4;
}
if(dtSector==14){
dtSector=10;
}
std::set<RPCDetId> rollsForThisDT = rollstoreDT[DTStationIndex(0,dtWheel,dtSector,dtStation)]; //It should be always 4
assert(rollsForThisDT.size()>=1);
for (std::set<RPCDetId>::iterator iteraRoll=rollsForThisDT.begin();iteraRoll != rollsForThisDT.end(); iteraRoll++){
const RPCRoll* rollasociated = rpcGeo->roll(*iteraRoll); //roll asociado a MB4
const BoundPlane & RPCSurfaceRB4 = rollasociated->surface(); //surface MB4
// RPCGeomServ rpcsrv(rpcId);
// std::string nameRoll = rpcsrv.name();
GlobalPoint CenterPointRollGlobal=RPCSurfaceRB4.toGlobal(LocalPoint(0,0,0));
LocalPoint CenterRollinMB4Frame = DTSurface4.toLocal(CenterPointRollGlobal); //In MB4
LocalPoint segmentPositionMB3inMB4Frame = DTSurface4.toLocal(segmentPositionMB3inGlobal); //In MB4
LocalVector segmentDirectionMB3inMB4Frame = DTSurface4.toLocal(segDirMB3inGlobalFrame); //In MB4
//The exptrapolation is done in MB4 frame. for local x and z is done from MB4,
float Dxz=CenterRollinMB4Frame.z();
float Xo4=segmentPositionMB4.x();
float dxl=segmentDirectionMB4.x(); //dx local for MB4 segment in MB4 Frame
float dzl=segmentDirectionMB4.z(); //dx local for MB4 segment in MB4 Frame
float X=Xo4+dxl*Dxz/dzl; //In MB4 frame
float Z=Dxz;//In MB4 frame
//for local y is done from MB3
float Yo34=segmentPositionMB3inMB4Frame.y();
float dy34 = segmentDirectionMB3inMB4Frame.y();
float dz34 = segmentDirectionMB3inMB4Frame.z();
float Dy=Dxz-(segmentPositionMB3inMB4Frame.z()); //Distance beetween the segment in MB3 and the RB4 surface
float Y=Yo34+dy34*Dy/dz34;//In MB4 Frame
const RectangularStripTopology* top_
=dynamic_cast<const RectangularStripTopology*>(&(rollasociated->topology())); //Topology roll asociated MB4
LocalPoint xmin = top_->localPosition(0.);
LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
float rsize = fabs( xmax.x()-xmin.x() );
float stripl = top_->stripLength();
float stripw = top_->pitch();
float extrapolatedDistance = sqrt((Y-Yo34)*(Y-Yo34)+Dy*Dy);
if(extrapolatedDistance<=MaxDrb4){
GlobalPoint GlobalPointExtrapolated = DTSurface4.toGlobal(LocalPoint(X,Y,Z));
LocalPoint PointExtrapolatedRPCFrame = RPCSurfaceRB4.toLocal(GlobalPointExtrapolated);
if(fabs(PointExtrapolatedRPCFrame.z()) < 5. &&
fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
RPCDetId rollId = rollasociated->id();
// RPCGeomServ rpcsrv(rollId);
// std::string nameRoll = rpcsrv.name();
// if(debug) std::cout<<"MB4 \t \t \t \t \t The RPCName is "<<nameRoll<<std::endl;
const float stripPredicted=
rollasociated->strip(LocalPoint(PointExtrapolatedRPCFrame.x(),PointExtrapolatedRPCFrame.y(),0.));
if(debug) std::cout<<"MB4 \t \t \t \t Candidate (from DT Segment) STRIP---> "<<stripPredicted<< std::endl;
//--------- HISTOGRAM STRIP PREDICTED FROM DT MB4 -------------------
std::map<std::string, MonitorElement*> meMap=meCollection[rollId.rawId()];
meIdDT.str("");
meIdDT<<"ExpectedOccupancyFromDT_"<<rollId.rawId();
meMap[meIdDT.str()]->Fill(stripPredicted);
//-------------------------------------------------
//-------RecHitPart Just For Residual--------
int countRecHits = 0;
int cluSize = 0;
float minres = 3000.;
if(debug) std::cout<<"MB4 \t \t \t \t Getting RecHits in Roll Asociated"<<std::endl;
typedef std::pair<RPCRecHitCollection::const_iterator, RPCRecHitCollection::const_iterator> rangeRecHits;
rangeRecHits recHitCollection = rpcHits->get(rollasociated->id());
RPCRecHitCollection::const_iterator recHit;
for (recHit = recHitCollection.first; recHit != recHitCollection.second ; recHit++) {
countRecHits++;
LocalPoint recHitPos=recHit->localPosition();
float res=PointExtrapolatedRPCFrame.x()- recHitPos.x();
if(debug) std::cout<<"DT \t \t \t \t \t Found Rec Hit at "<<res<<"cm of the prediction."<<std::endl;
if(fabs(res)<fabs(minres)){
minres=res;
cluSize = recHit->clusterSize();
}
}
if(countRecHits==0){
if(debug) std::cout <<"MB4 \t \t \t \t \t \t THIS ROLL DOESN'T HAVE ANY RECHIT"<<std::endl;
}else{
assert(minres!=3000);
if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
assert(rollId.station()==4);
if(cluSize==1*dupli){ hGlobalResClu1La[5]->Fill(minres);}
else if(cluSize==2*dupli){ hGlobalResClu2La[5]->Fill(minres);}
else if(cluSize==3*dupli){ hGlobalResClu3La[5]->Fill(minres);}
meIdRPC.str("");
meIdRPC<<"RPCDataOccupancyFromDT_"<<rollId.rawId();
meMap[meIdRPC.str()]->Fill(stripPredicted);
}
}
}else{
if(debug) std::cout<<"MB4 \t \t \t \t No the prediction is outside of this roll"<<std::endl;
}
}//Condition for the right match
else{
if(debug) std::cout<<"MB4 \t \t \t No, Exrtrapolation too long!, canceled"<<std::endl;
}
}//loop over all the rollsasociated
}else{
if(debug) std::cout<<"MB4 \t \t \t \t I found segments in MB4 and MB3 adjacent or same wheel and sector but not compatibles Diferent Directions"<<std::endl;
}
}else{//if dtid3.station()==3&&dtid3.sector()==DTId.sector()&&dtid3.wheel()==DTId.wheel()&&segMB3->dim()==4
if(debug) std::cout<<"MB4 \t \t \t No the same station or same wheel or segment dim in mb3 not 4D"<<std::endl;
}
}//loop over all the segments looking for one in MB3
}else{
if(debug) std::cout<<"MB4 \t \t \t Is NOT a 2D Segment"<<std::endl;
}
}else{
if(debug) std::cout<<"MB4 \t \t \t \t There is not just one segment or is not in station 4"<<std::endl;
}//De aca para abajo esta en dtpart.inl
}
}else{
if(debug) std::cout<<"MB4 \t This event doesn't have 4D Segment"<<std::endl;
}
}
if(inclcsc){
if(debug) std::cout <<"\t Getting the CSC Segments"<<std::endl;
edm::Handle<CSCSegmentCollection> allCSCSegments;
iEvent.getByLabel(cscSegments, allCSCSegments);
if(allCSCSegments.isValid()){
if(allCSCSegments->size()>0){
statistics->Fill(18);
if(debug) std::cout<<"CSC \t Number of CSC Segments in this event = "<<allCSCSegments->size()<<std::endl;
std::map<CSCDetId,int> CSCSegmentsCounter;
CSCSegmentCollection::const_iterator segment;
int segmentsInThisEventInTheEndcap=0;
for (segment = allCSCSegments->begin();segment!=allCSCSegments->end(); ++segment){
CSCSegmentsCounter[segment->cscDetId()]++;
segmentsInThisEventInTheEndcap++;
}
statistics->Fill(allCSCSegments->size()+18);
if(debug) std::cout<<"CSC \t loop over all the CSCSegments "<<std::endl;
for (segment = allCSCSegments->begin();segment!=allCSCSegments->end(); ++segment){
CSCDetId CSCId = segment->cscDetId();
if(CSCSegmentsCounter[CSCId]==1 && CSCId.station()!=4 && CSCId.ring()!=1 && allCSCSegments->size()>=2){
if(debug) std::cout<<"CSC \t \t yes"<<std::endl;
int cscEndCap = CSCId.endcap();
int cscStation = CSCId.station();
int cscRing = CSCId.ring();
// int cscChamber = CSCId.chamber();
int rpcRegion = 1; if(cscEndCap==2) rpcRegion= -1;//Relacion entre las endcaps
int rpcRing = cscRing;
if(cscRing==4)rpcRing =1;
int rpcStation = cscStation;
int rpcSegment = CSCId.chamber();
LocalPoint segmentPosition= segment->localPosition();
LocalVector segmentDirection=segment->localDirection();
float dz=segmentDirection.z();
if(debug) std::cout<<"CSC \t \t Is a good Segment? dim = 4, 4 <= nRecHits <= 10 Incident angle int range 45 < "<<acos(dz)*180/3.1415926<<" < 135? "<<std::endl;
if(segment->dimension()==4 && (segment->nRecHits()<=10 && segment->nRecHits()>=4)&& acos(dz)*180/3.1415926 > 45. && acos(dz)*180/3.1415926 < 160. ){
float Xo=segmentPosition.x();
float Yo=segmentPosition.y();
float Zo=segmentPosition.z();
float dx=segmentDirection.x();
float dy=segmentDirection.y();
float dz=segmentDirection.z();
if(debug) std::cout<<"CSC \t \t Getting chamber from Geometry"<<std::endl;
const CSCChamber* TheChamber=cscGeo->chamber(CSCId);
if(debug) std::cout<<"CSC \t \t Getting ID from Chamber"<<std::endl;
const CSCDetId TheId=TheChamber->id();
if(debug) std::cout<<"CSC \t \t Printing The Id"<<TheId<<std::endl;
std::set<RPCDetId> rollsForThisCSC = rollstoreCSC[CSCStationIndex(rpcRegion,rpcStation,rpcRing,rpcSegment)];
if(debug) std::cout<<"CSC \t \t Number of rolls for this CSC = "<<rollsForThisCSC.size()<<std::endl;
if(rpcRing!=1&&rpcStation!=4){
//Loop over all the rolls
for (std::set<RPCDetId>::iterator iteraRoll = rollsForThisCSC.begin();iteraRoll != rollsForThisCSC.end(); iteraRoll++){
const RPCRoll* rollasociated = rpcGeo->roll(*iteraRoll);
RPCDetId rpcId = rollasociated->id();
const BoundPlane & RPCSurface = rollasociated->surface();
GlobalPoint CenterPointRollGlobal = RPCSurface.toGlobal(LocalPoint(0,0,0));
LocalPoint CenterRollinCSCFrame = TheChamber->toLocal(CenterPointRollGlobal);
float D=CenterRollinCSCFrame.z();
float X=Xo+dx*D/dz;
float Y=Yo+dy*D/dz;
float Z=D;
const TrapezoidalStripTopology* top_=dynamic_cast<const TrapezoidalStripTopology*>(&(rollasociated->topology()));
LocalPoint xmin = top_->localPosition(0.);
if(debug) std::cout<<"CSC \t \t \t xmin of this Roll "<<xmin<<"cm"<<std::endl;
LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
if(debug) std::cout<<"CSC \t \t \t xmax of this Roll "<<xmax<<"cm"<<std::endl;
float rsize = fabs( xmax.x()-xmin.x() );
if(debug) std::cout<<"CSC \t \t \t Roll Size "<<rsize<<"cm"<<std::endl;
float stripl = top_->stripLength();
float stripw = top_->pitch();
float extrapolatedDistance = sqrt((X-Xo)*(X-Xo)+(Y-Yo)*(Y-Yo)+(Z-Zo)*(Z-Zo));
if(extrapolatedDistance<=MaxD){
GlobalPoint GlobalPointExtrapolated=TheChamber->toGlobal(LocalPoint(X,Y,Z));
LocalPoint PointExtrapolatedRPCFrame = RPCSurface.toLocal(GlobalPointExtrapolated);
if(fabs(PointExtrapolatedRPCFrame.z()) < 10. &&
fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
RPCDetId rollId = rollasociated->id();
RPCGeomServ rpcsrv(rollId);
std::string nameRoll = rpcsrv.name();
if(debug) std::cout<<"CSC \t \t \t \t The RPCName is "<<nameRoll<<std::endl;
const float stripPredicted =
rollasociated->strip(LocalPoint(PointExtrapolatedRPCFrame.x(),PointExtrapolatedRPCFrame.y(),0.));
if(debug) std::cout<<"CSC \t \t \t \t \t Candidate"<<rollId<<" "<<"(from CSC Segment) STRIP---> "<<stripPredicted<< std::endl;
//--------- HISTOGRAM STRIP PREDICTED FROM CSC -------------------
std::map<std::string, MonitorElement*> meMap=meCollection[rpcId.rawId()];
meIdCSC.str("");
meIdCSC<<"ExpectedOccupancyFromCSC_"<<rollId.rawId();
meMap[meIdCSC.str()]->Fill(stripPredicted);
//--------------------------------------------------------------------
//-------RecHitPart Just For Residual--------
int cluSize = 0;
int countRecHits = 0;
float minres = 3000.;
if(debug) std::cout<<"CSC \t \t \t \t \t Getting RecHits in Roll Asociated"<<std::endl;
typedef std::pair<RPCRecHitCollection::const_iterator, RPCRecHitCollection::const_iterator> rangeRecHits;
rangeRecHits recHitCollection = rpcHits->get(rollasociated->id());
RPCRecHitCollection::const_iterator recHit;
for (recHit = recHitCollection.first; recHit != recHitCollection.second ; recHit++) {
countRecHits++;
LocalPoint recHitPos=recHit->localPosition();
float res=PointExtrapolatedRPCFrame.x()- recHitPos.x();
if(debug) std::cout<<"CSC \t \t \t \t \t \t Found Rec Hit at "<<res<<"cm of the prediction."<<std::endl;
if(fabs(res)<fabs(minres)){
minres=res;
cluSize = recHit->clusterSize();
if(debug) std::cout<<"CSC \t \t \t \t \t \t \t New Min Res "<<res<<"cm."<<std::endl;
}
}
if(countRecHits==0){
if(debug) std::cout <<"CSC \t \t \t \t \t THIS ROLL DOESN'T HAVE ANY RECHIT"<<std::endl;
}else{
assert(minres!=3000);
if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
if(debug) std::cout<<"CSC \t \t \t \t \t \t True!"<<std::endl;
if(rollId.ring()==2&&rollId.roll()==1){
if(cluSize==1*dupli) hGlobalResClu1R2A->Fill(minres);
else if(cluSize==2*dupli) hGlobalResClu2R2A->Fill(minres);
else if(cluSize==3*dupli) hGlobalResClu3R2A->Fill(minres);
}
else if(rollId.ring()==2&&rollId.roll()==2){
if(cluSize==1*dupli) hGlobalResClu1R2B->Fill(minres);
else if(cluSize==2*dupli) hGlobalResClu2R2B->Fill(minres);
else if(cluSize==3*dupli) hGlobalResClu3R2B->Fill(minres);
}
else if(rollId.ring()==2&&rollId.roll()==3){
if(cluSize==1*dupli) hGlobalResClu1R2C->Fill(minres);
else if(cluSize==2*dupli) hGlobalResClu2R2C->Fill(minres);
else if(cluSize==3*dupli) hGlobalResClu3R2C->Fill(minres);
}
else if(rollId.ring()==3&&rollId.roll()==1){
if(cluSize==1*dupli) hGlobalResClu1R3A->Fill(minres);
else if(cluSize==2*dupli) hGlobalResClu2R3A->Fill(minres);
else if(cluSize==3*dupli) hGlobalResClu3R3A->Fill(minres);
}
else if(rollId.ring()==3&&rollId.roll()==2){
if(cluSize==1*dupli) hGlobalResClu1R3B->Fill(minres);
else if(cluSize==2*dupli) hGlobalResClu2R3B->Fill(minres);
else if(cluSize==3*dupli) hGlobalResClu3R3B->Fill(minres);
}
else if(rollId.ring()==3&&rollId.roll()==3){if(cluSize==1*dupli) hGlobalResClu1R3C->Fill(minres); if(cluSize==2*dupli) hGlobalResClu2R3C->Fill(minres); if(cluSize==3*dupli) hGlobalResClu3R3C->Fill(minres);
}
meIdRPC.str("");
meIdRPC<<"RPCDataOccupancyFromCSC_"<<rollId.rawId();
meMap[meIdRPC.str()]->Fill(stripPredicted);
}
}
}else{
if(debug) std::cout<<"CSC \t \t \t \t No the prediction is outside of this roll"<<std::endl;
}//Condition for the right match
}else{//if extrapolation distance D is not too long
if(debug) std::cout<<"CSC \t \t \t No, Exrtrapolation too long!, canceled"<<std::endl;
}//D so big
}//loop over the rolls asociated
}//Condition over the startup geometry!!!!
}//Is the segment 4D?
}else{
if(debug) std::cout<<"CSC \t \t More than one segment in this chamber, or we are in Station Ring 1 or in Station 4"<<std::endl;
}
}
}else{
if(debug) std::cout<<"CSC This Event doesn't have any CSCSegment"<<std::endl;
}
}
}
}
| void RPCEfficiency::beginJob | ( | void | ) | [virtual] |
| void RPCEfficiency::beginRun | ( | const edm::Run & | run, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 162 of file RPCEfficiency.cc.
References bookDetUnitSeg(), gather_cfg::cout, debug, folderPath, edm::EventSetup::get(), inclcsc, incldt, incldtMB4, meCollection, alignCSCRings::r, DetId::rawId(), RPCDetId::region(), RPCDetId::ring(), relativeConstraints::ring, RPCChamber::rolls(), rollstoreCSC, rollstoreDT, RPCDetId::sector(), RPCGeomServ::segment(), relativeConstraints::station, and RPCDetId::station().
{
edm::ESHandle<RPCGeometry> rpcGeo;
iSetup.get<MuonGeometryRecord>().get(rpcGeo);
for (TrackingGeometry::DetContainer::const_iterator it=rpcGeo->dets().begin();it<rpcGeo->dets().end();it++){
if(dynamic_cast< RPCChamber* >( *it ) != 0 ){
RPCChamber* ch = dynamic_cast< RPCChamber* >( *it );
std::vector< const RPCRoll*> roles = (ch->rolls());
for(std::vector<const RPCRoll*>::const_iterator r = roles.begin();r != roles.end(); ++r){
RPCDetId rpcId = (*r)->id();
int region=rpcId.region();
//booking all histograms
// std::string nameRoll = rpcsrv.name();
if(debug) std::cout<<"Booking for "<<rpcId.rawId()<<std::endl;
bookDetUnitSeg(rpcId,(*r)->nstrips(),folderPath+"MuonSegEff/", meCollection[rpcId.rawId()] );
if(region==0&&(incldt||incldtMB4)){
//std::cout<<"--Filling the dtstore"<<rpcId<<std::endl;
int wheel=rpcId.ring();
int sector=rpcId.sector();
int station=rpcId.station();
DTStationIndex ind(region,wheel,sector,station);
std::set<RPCDetId> myrolls;
if (rollstoreDT.find(ind)!=rollstoreDT.end()) myrolls=rollstoreDT[ind];
myrolls.insert(rpcId);
rollstoreDT[ind]=myrolls;
}else if(region!=0 && inclcsc){
int station=rpcId.station();
int ring=rpcId.ring();
int cscring=ring;
int cscstation=station;
RPCGeomServ rpcsrv(rpcId);
int rpcsegment = rpcsrv.segment();
int cscchamber = rpcsegment;
if((station==2||station==3)&&ring==3){
cscring = 2;
}
CSCStationIndex ind(region,cscstation,cscring,cscchamber);
std::set<RPCDetId> myrolls;
if (rollstoreCSC.find(ind)!=rollstoreCSC.end()){
myrolls=rollstoreCSC[ind];
}
myrolls.insert(rpcId);
rollstoreCSC[ind]=myrolls;
}
}
}
}
for (TrackingGeometry::DetContainer::const_iterator it=rpcGeo->dets().begin();it<rpcGeo->dets().end();it++){
if( dynamic_cast< RPCChamber* >( *it ) != 0 ){
RPCChamber* ch = dynamic_cast< RPCChamber* >( *it );
std::vector< const RPCRoll*> roles = (ch->rolls());
for(std::vector<const RPCRoll*>::const_iterator r = roles.begin();r != roles.end(); ++r){
RPCDetId rpcId = (*r)->id();
int region=rpcId.region();
if(region!=0 && inclcsc && (rpcId.ring()==2 || rpcId.ring()==3)){
int region=rpcId.region();
int station=rpcId.station();
int ring=rpcId.ring();
int cscring = ring;
if((station==2||station==3)&&ring==3) cscring = 2; //CSC Ring 2 covers rpc ring 2 & 3
int cscstation=station;
RPCGeomServ rpcsrv(rpcId);
int rpcsegment = rpcsrv.segment();
int cscchamber = rpcsegment+1;
if(cscchamber==37)cscchamber=1;
CSCStationIndex ind(region,cscstation,cscring,cscchamber);
std::set<RPCDetId> myrolls;
if (rollstoreCSC.find(ind)!=rollstoreCSC.end())myrolls=rollstoreCSC[ind];
myrolls.insert(rpcId);
rollstoreCSC[ind]=myrolls;
cscchamber = rpcsegment-1;
if(cscchamber==0)cscchamber=36;
CSCStationIndex indDos(region,cscstation,cscring,cscchamber);
std::set<RPCDetId> myrollsDos;
if (rollstoreCSC.find(indDos)!=rollstoreCSC.end()) myrollsDos=rollstoreCSC[indDos];
myrollsDos.insert(rpcId);
rollstoreCSC[indDos]=myrollsDos;
}
}
}
}
}//beginRun
| void RPCEfficiency::bookDetUnitSeg | ( | RPCDetId & | detId, |
| int | nstrips, | ||
| std::string | folder, | ||
| std::map< std::string, MonitorElement * > & | meMap | ||
| ) |
Definition at line 13 of file RPCEfficiencyBook.cc.
References DQMStore::book1D(), dbe, DetId::rawId(), RPCDetId::region(), and DQMStore::setCurrentFolder().
Referenced by beginRun().
{
//std::map<std::string, MonitorElement*> meMap;
dbe->setCurrentFolder(folder);
char meId [128];
char meTitle [128];
int rawId = detId.rawId();
//Begin booking DT
if(detId.region()==0) {
sprintf(meId,"ExpectedOccupancyFromDT_%d",rawId);
sprintf(meTitle,"ExpectedOccupancyFromDT_for_%d",rawId);
meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
sprintf(meId,"RPCDataOccupancyFromDT_%d",rawId);
sprintf(meTitle,"RPCDataOccupancyFromDT_for_%d",rawId);
meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
// sprintf(meId,"BXDistribution_%d",rawId);
// sprintf(meTitle,"BXDistribution_for_%d",rawId);
// meMap[meId] = dbe->book1D(meId, meTitle, 11,-5.5, 5.5);
}else{
//std::cout<<"Booking for the EndCap"<<detUnitLabel<<std::endl;
sprintf(meId,"ExpectedOccupancyFromCSC_%d",rawId);
sprintf(meTitle,"ExpectedOccupancyFromCSC_for_%d",rawId);
meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
sprintf(meId,"RPCDataOccupancyFromCSC_%d",rawId);
sprintf(meTitle,"RPCDataOccupancyFromCSC_for_%d",rawId);
meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
// sprintf(meId,"BXDistribution_%d",rawId);
// sprintf(meTitle,"BXDistribution_for_%d",rawId);
// meMap[meId] = dbe->book1D(meId, meTitle, 11,-5.5, 5.5);
}
//return meMap;
}
| void RPCEfficiency::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 907 of file RPCEfficiency.cc.
References dbe.
{
dbe =0;
}
| void RPCEfficiency::endRun | ( | const edm::Run & | r, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 900 of file RPCEfficiency.cc.
References dbe, EffRootFileName, EffSaveRootFile, and DQMStore::save().
{
if (EffSaveRootFile){
dbe->save(EffRootFileName);
}
}
std::vector<std::map<RPCDetId, int> > RPCEfficiency::counter [private] |
Definition at line 158 of file RPCEfficiency.h.
edm::InputTag RPCEfficiency::cscSegments [private] |
Definition at line 175 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
DQMStore* RPCEfficiency::dbe [private] |
Definition at line 187 of file RPCEfficiency.h.
Referenced by bookDetUnitSeg(), endJob(), endRun(), and RPCEfficiency().
bool RPCEfficiency::debug [private] |
Definition at line 164 of file RPCEfficiency.h.
Referenced by analyze(), beginRun(), and RPCEfficiency().
edm::InputTag RPCEfficiency::dt4DSegments [private] |
Definition at line 176 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
int RPCEfficiency::dupli [private] |
Definition at line 172 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
std::string RPCEfficiency::EffRootFileName [private] |
Definition at line 186 of file RPCEfficiency.h.
Referenced by endRun(), and RPCEfficiency().
bool RPCEfficiency::EffSaveRootFile [private] |
Definition at line 185 of file RPCEfficiency.h.
Referenced by endRun(), and RPCEfficiency().
std::string RPCEfficiency::folderPath [private] |
Definition at line 179 of file RPCEfficiency.h.
Referenced by beginRun(), and RPCEfficiency().
Definition at line 131 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 141 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 140 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 139 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 138 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 137 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 136 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 132 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 148 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 147 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 146 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 145 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 144 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 143 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 133 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 155 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 154 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 153 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 152 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 151 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 150 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
bool RPCEfficiency::inclcsc [private] |
Definition at line 163 of file RPCEfficiency.h.
Referenced by analyze(), beginRun(), and RPCEfficiency().
bool RPCEfficiency::incldt [private] |
Definition at line 161 of file RPCEfficiency.h.
Referenced by analyze(), beginRun(), and RPCEfficiency().
bool RPCEfficiency::incldtMB4 [private] |
Definition at line 162 of file RPCEfficiency.h.
Referenced by analyze(), beginRun(), and RPCEfficiency().
bool RPCEfficiency::inves [private] |
Definition at line 166 of file RPCEfficiency.h.
double RPCEfficiency::MaxD [private] |
Definition at line 170 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
double RPCEfficiency::MaxDrb4 [private] |
Definition at line 171 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
std::map<int, std::map<std::string, MonitorElement*> > RPCEfficiency::meCollection [private] |
Definition at line 183 of file RPCEfficiency.h.
Referenced by analyze(), and beginRun().
double RPCEfficiency::MinCosAng [private] |
Definition at line 169 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
double RPCEfficiency::rangestrips [private] |
Definition at line 167 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
double RPCEfficiency::rangestripsRB4 [private] |
Definition at line 168 of file RPCEfficiency.h.
Referenced by RPCEfficiency().
std::string RPCEfficiency::rollseff [private] |
Definition at line 181 of file RPCEfficiency.h.
| std::map<CSCStationIndex,std::set<RPCDetId> > RPCEfficiency::rollstoreCSC |
Definition at line 118 of file RPCEfficiency.h.
Referenced by analyze(), and beginRun().
| std::map<DTStationIndex,std::set<RPCDetId> > RPCEfficiency::rollstoreDT |
Definition at line 117 of file RPCEfficiency.h.
Referenced by analyze(), and beginRun().
edm::InputTag RPCEfficiency::RPCRecHitLabel_ [private] |
Definition at line 177 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 124 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
std::vector<int> RPCEfficiency::totalcounter [private] |
Definition at line 159 of file RPCEfficiency.h.
1.7.3