|
|
|
#include <RPCEfficiency.h>
Class for RPC Monitoring: uses RPCDigi and DT and CSC Segments.
Definition at line 104 of file RPCEfficiency.h.
| RPCEfficiency::RPCEfficiency | ( | const edm::ParameterSet & | iConfig | ) | [explicit] |
Definition at line 49 of file RPCEfficiency.cc.
References gather_cfg::cout, cscSegments, dbe, debug, dt4DSegments, dupli, EffRootFileName, EffSaveRootFile, folderPath, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hGlobalResClu1La1, hGlobalResClu1La2, hGlobalResClu1La3, hGlobalResClu1La4, hGlobalResClu1La5, hGlobalResClu1La6, hGlobalResClu1R2A, hGlobalResClu1R2B, hGlobalResClu1R2C, hGlobalResClu1R3A, hGlobalResClu1R3B, hGlobalResClu1R3C, hGlobalResClu2La1, hGlobalResClu2La2, hGlobalResClu2La3, hGlobalResClu2La4, hGlobalResClu2La5, hGlobalResClu2La6, hGlobalResClu2R2A, hGlobalResClu2R2B, hGlobalResClu2R2C, hGlobalResClu2R3A, hGlobalResClu2R3B, hGlobalResClu2R3C, hGlobalResClu3La1, hGlobalResClu3La2, hGlobalResClu3La3, hGlobalResClu3La4, hGlobalResClu3La5, hGlobalResClu3La6, hGlobalResClu3R2A, hGlobalResClu3R2B, hGlobalResClu3R2C, hGlobalResClu3R3A, hGlobalResClu3R3B, hGlobalResClu3R3C, inclcsc, incldt, incldtMB4, MaxD, MaxDrb4, MinCosAng, muonRPCDigis, nameInLog, ofrej, cmsCodeRules::cppFunctionSkipper::operator, rangestrips, rangestripsRB4, RPCRecHitLabel_, 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");
RPCRecHitLabel_ = iConfig.getParameter<edm::InputTag>("RecHitLabel");
cscSegments=iConfig.getUntrackedParameter<std::string>("cscSegments","cscSegments");
dt4DSegments=iConfig.getUntrackedParameter<std::string>("dt4DSegments","dt4DSegments");
folderPath=iConfig.getUntrackedParameter<std::string>("folderPath","RPC/RPCEfficiency/");
nameInLog = iConfig.getUntrackedParameter<std::string>("moduleLogName", "RPC_Eff");
EffSaveRootFile = iConfig.getUntrackedParameter<bool>("EffSaveRootFile", false);
EffRootFileName = iConfig.getUntrackedParameter<std::string>("EffRootFileName", "RPCEfficiency.root");
//Interface
dbe = edm::Service<DQMStore>().operator->();
std::string folder = folderPath+"MuonSegEff";
dbe->setCurrentFolder(folder);
statistics = dbe->book1D("Statistics","All Statistics",33,0.5,33.5);
if(debug) std::cout<<"booking Global histograms with "<<folderPath<<std::endl;
folder = folderPath+"MuonSegEff/"+"Residuals/Barrel";
dbe->setCurrentFolder(folder);
//Barrel
hGlobalResClu1La1 = dbe->book1D("GlobalResidualsClu1La1","RPC Residuals Layer 1 Cluster Size 1",101,-10.,10.);
hGlobalResClu1La2 = dbe->book1D("GlobalResidualsClu1La2","RPC Residuals Layer 2 Cluster Size 1",101,-10.,10.);
hGlobalResClu1La3 = dbe->book1D("GlobalResidualsClu1La3","RPC Residuals Layer 3 Cluster Size 1",101,-10.,10.);
hGlobalResClu1La4 = dbe->book1D("GlobalResidualsClu1La4","RPC Residuals Layer 4 Cluster Size 1",101,-10.,10.);
hGlobalResClu1La5 = dbe->book1D("GlobalResidualsClu1La5","RPC Residuals Layer 5 Cluster Size 1",101,-10.,10.);
hGlobalResClu1La6 = dbe->book1D("GlobalResidualsClu1La6","RPC Residuals Layer 6 Cluster Size 1",101,-10.,10.);
hGlobalResClu2La1 = dbe->book1D("GlobalResidualsClu2La1","RPC Residuals Layer 1 Cluster Size 2",101,-10.,10.);
hGlobalResClu2La2 = dbe->book1D("GlobalResidualsClu2La2","RPC Residuals Layer 2 Cluster Size 2",101,-10.,10.);
hGlobalResClu2La3 = dbe->book1D("GlobalResidualsClu2La3","RPC Residuals Layer 3 Cluster Size 2",101,-10.,10.);
hGlobalResClu2La4 = dbe->book1D("GlobalResidualsClu2La4","RPC Residuals Layer 4 Cluster Size 2",101,-10.,10.);
hGlobalResClu2La5 = dbe->book1D("GlobalResidualsClu2La5","RPC Residuals Layer 5 Cluster Size 2",101,-10.,10.);
hGlobalResClu2La6 = dbe->book1D("GlobalResidualsClu2La6","RPC Residuals Layer 6 Cluster Size 2",101,-10.,10.);
hGlobalResClu3La1 = dbe->book1D("GlobalResidualsClu3La1","RPC Residuals Layer 1 Cluster Size 3",101,-10.,10.);
hGlobalResClu3La2 = dbe->book1D("GlobalResidualsClu3La2","RPC Residuals Layer 2 Cluster Size 3",101,-10.,10.);
hGlobalResClu3La3 = dbe->book1D("GlobalResidualsClu3La3","RPC Residuals Layer 3 Cluster Size 3",101,-10.,10.);
hGlobalResClu3La4 = dbe->book1D("GlobalResidualsClu3La4","RPC Residuals Layer 4 Cluster Size 3",101,-10.,10.);
hGlobalResClu3La5 = dbe->book1D("GlobalResidualsClu3La5","RPC Residuals Layer 5 Cluster Size 3",101,-10.,10.);
hGlobalResClu3La6 = dbe->book1D("GlobalResidualsClu3La6","RPC Residuals Layer 6 Cluster Size 3",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.);
if(debug) ofrej.open("rejected.txt");
if(debug) std::cout<<"Rejected done"<<std::endl;
}
| RPCEfficiency::~RPCEfficiency | ( | ) |
Definition at line 267 of file RPCEfficiency.cc.
{
}
| void RPCEfficiency::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 272 of file RPCEfficiency.cc.
References PV3DBase< T, PVType, FrameType >::barePhi(), CSCDetId::chamber(), gather_cfg::cout, cscGeo, cscSegments, funct::D, debug, distsector(), dt4DSegments, dtGeo, dupli, CSCDetId::endcap(), edm::EventID::event(), MonitorElement::Fill(), edm::Event::getByLabel(), hGlobalResClu1La1, hGlobalResClu1La2, hGlobalResClu1La3, hGlobalResClu1La4, hGlobalResClu1La5, hGlobalResClu1La6, hGlobalResClu1R2A, hGlobalResClu1R2B, hGlobalResClu1R2C, hGlobalResClu1R3A, hGlobalResClu1R3B, hGlobalResClu1R3C, hGlobalResClu2La1, hGlobalResClu2La2, hGlobalResClu2La3, hGlobalResClu2La4, hGlobalResClu2La5, hGlobalResClu2La6, hGlobalResClu2R2A, hGlobalResClu2R2B, hGlobalResClu2R2C, hGlobalResClu2R3A, hGlobalResClu2R3B, hGlobalResClu2R3C, hGlobalResClu3La1, hGlobalResClu3La2, hGlobalResClu3La3, hGlobalResClu3La4, hGlobalResClu3La5, hGlobalResClu3La6, hGlobalResClu3R2A, hGlobalResClu3R2B, hGlobalResClu3R2C, hGlobalResClu3R3A, hGlobalResClu3R3B, hGlobalResClu3R3C, CSCChamber::id(), edm::EventBase::id(), RPCRoll::id(), inclcsc, incldt, incldtMB4, edm::HandleBase::isValid(), RPCDetId::layer(), TrapezoidalStripTopology::localPosition(), RectangularStripTopology::localPosition(), MaxD, MaxDrb4, meCollection, MinCosAng, RPCGeomServ::name(), RPCRoll::nstrips(), ofrej, PV3DBase< T, PVType, FrameType >::perp(), RectangularStripTopology::pitch(), TrapezoidalStripTopology::pitch(), rangestrips, DetId::rawId(), RPCDetId::ring(), CSCDetId::ring(), RPCDetId::roll(), rollstoreCSC, rollstoreDT, rpcGeo, RPCRecHitLabel_, edm::EventID::run(), DTChamberId::sector(), mathSSE::sqrt(), RPCDetId::station(), CSCDetId::station(), DTChamberId::station(), statistics, RPCRoll::strip(), RectangularStripTopology::stripLength(), TrapezoidalStripTopology::stripLength(), 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().
{
statistics->Fill(1);
char meIdRPC [128];
char meIdDT [128];
char meIdCSC [128];
if(debug) std::cout <<"\t Getting the RPC RecHits"<<std::endl;
edm::Handle<RPCRecHitCollection> rpcHits;
iEvent.getByLabel(RPCRecHitLabel_,rpcHits);
if(rpcHits.isValid()){
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(debug) std::cout<<"DT \t \t This Segment is in Chamber id: "<<DTId<<std::endl;
if(debug) std::cout<<"DT \t \t Number of segments in this DT = "<<DTSegmentCounter[DTId]<<std::endl;
if(debug) std::cout<<"DT \t \t Is the only one in this DT? and is not in the 4th Station?"<<std::endl;
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(debug) std::cout<<"DT \t \t Is the segment 4D?"<<std::endl;
if(segment->dimension()==4){
if(debug) std::cout<<"DT \t \t yes"<<std::endl;
if(debug) std::cout<<"DT \t \t DT Segment Dimension "<<segment->dimension()<<std::endl;
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 Number of rolls for this DT = "<<rollsForThisDT.size()<<std::endl;
assert(rollsForThisDT.size()>=1);
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();
if(debug) std::cout<<"DT \t \t \t RollName: "<<nameRoll<<std::endl;
if(debug) std::cout<<"DT \t \t \t Doing the extrapolation to this roll"<<std::endl;
if(debug) std::cout<<"DT \t \t \t DT Segment Direction in DTLocal "<<segmentDirection<<std::endl;
if(debug) std::cout<<"DT \t \t \t DT Segment Point in DTLocal "<<segmentPosition<<std::endl;
GlobalPoint CenterPointRollGlobal = RPCSurface.toGlobal(LocalPoint(0,0,0));
LocalPoint CenterRollinDTFrame = DTSurface.toLocal(CenterPointRollGlobal);
if(debug) std::cout<<"DT \t \t \t Center (0,0,0) Roll In DTLocal"<<CenterRollinDTFrame<<std::endl;
if(debug) std::cout<<"DT \t \t \t Center (0,0,0) of the Roll in Global"<<CenterPointRollGlobal<<std::endl;
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();
if(debug) std::cout<<"DT \t \t \t Strip Lenght "<<stripl<<"cm"<<std::endl;
if(debug) std::cout<<"DT \t \t \t Strip Width "<<stripw<<"cm"<<std::endl;
if(debug) std::cout<<"DT \t \t \t X Predicted in DTLocal= "<<X<<"cm"<<std::endl;
if(debug) std::cout<<"DT \t \t \t Y Predicted in DTLocal= "<<Y<<"cm"<<std::endl;
if(debug) std::cout<<"DT \t \t \t Z Predicted in DTLocal= "<<Z<<"cm"<<std::endl;
float extrapolatedDistance = sqrt((X-Xo)*(X-Xo)+(Y-Yo)*(Y-Yo)+(Z-Zo)*(Z-Zo));
if(debug) std::cout<<"DT \t \t \t Is the distance of extrapolation less than MaxD? ="<<extrapolatedDistance<<"cm"<<"MaxD="<<MaxD<<"cm"<<std::endl;
if(extrapolatedDistance<=MaxD){
if(debug) std::cout<<"DT \t \t \t yes"<<std::endl;
GlobalPoint GlobalPointExtrapolated = DTSurface.toGlobal(LocalPoint(X,Y,Z));
if(debug) std::cout<<"DT \t \t \t Point ExtraPolated in Global"<<GlobalPointExtrapolated<< std::endl;
LocalPoint PointExtrapolatedRPCFrame = RPCSurface.toLocal(GlobalPointExtrapolated);
if(debug) std::cout<<"DT \t \t \t Point Extrapolated in RPCLocal"<<PointExtrapolatedRPCFrame<< std::endl;
if(debug) std::cout<<"DT \t \t \t Corner of the Roll = ("<<rsize*0.5<<","<<stripl*0.5<<")"<<std::endl;
if(debug) std::cout<<"DT \t \t \t Info About the Point Extrapolated in X Abs ("<<fabs(PointExtrapolatedRPCFrame.x())<<","
<<fabs(PointExtrapolatedRPCFrame.y())<<","<<fabs(PointExtrapolatedRPCFrame.z())<<")"<<std::endl;
if(debug) std::cout<<"DT \t \t \t Does the extrapolation go inside this roll?"<<std::endl;
if(fabs(PointExtrapolatedRPCFrame.z()) < 10. &&
fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
if(debug) std::cout<<"DT \t \t \t \t yes"<<std::endl;
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()];
sprintf(meIdDT,"ExpectedOccupancyFromDT_%d",rollId.rawId());
if(debug) std::cout<<"DT \t \t \t \t Filling Expected for "<<meIdDT<<" with "<<stripPredicted<<std::endl;
if(fabs(stripPredicted-rollasociated->nstrips())<1.) if(debug) std::cout<<"DT \t \t \t \t Extrapolating near last strip, Event"<<iEvent.id()<<" stripPredicted="<<stripPredicted<<" Number of strips="<<rollasociated->nstrips()<<std::endl;
if(fabs(stripPredicted)<1.) if(debug) std::cout<<"DT \t \t \t \t Extrapolating near first strip, Event"<<iEvent.id()<<" stripPredicted="<<stripPredicted<<" Number of strips="<<rollasociated->nstrips()<<std::endl;
meMap[meIdDT]->Fill(stripPredicted);
//-----------------------------------------------------
//-------RecHitPart Just For Residual--------
int countRecHits = 0;
int cluSize = 0;
float minres = 3000.;
if(debug) std::cout<<"DT \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++;
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;
}
}
bool anycoincidence=false;
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(debug) std::cout<<"DT \t \t \t \t \t PointExtrapolatedRPCFrame.x="<<PointExtrapolatedRPCFrame.x()<<" Minimal Residual="<<minres<<std::endl;
if(debug) std::cout<<"DT \t \t \t \t \t Minimal Residual less than stripw*rangestrips? minres="<<minres<<" range="<<rangestrips<<" stripw="<<stripw<<" cluSize"<<cluSize<<" <=compare minres with"<<(rangestrips+cluSize*0.5)*stripw<<std::endl;
if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
if(debug) std::cout<<"DT \t \t \t \t \t \t True!"<<std::endl;
anycoincidence=true;
}
}
if(anycoincidence){
if(debug) std::cout<<"DT \t \t \t \t \t At least one RecHit inside the range, Predicted="<<stripPredicted<<" minres="<<minres<<"cm range="<<rangestrips<<"strips stripw="<<stripw<<"cm"<<std::endl;
if(debug) std::cout<<"DT \t \t \t \t \t Norm of Cosine Directors="<<dx*dx+dy*dy+dz*dz<<"~1?"<<std::endl;
float cosal = dx/sqrt(dx*dx+dz*dz);
if(debug) std::cout<<"DT \t \t \t \t \t Angle="<<acos(cosal)*180/3.1415926<<" degree"<<std::endl;
if(debug) std::cout<<"DT \t \t \t \t \t Filling the Residuals Histogram for globals with "<<minres<<"And the angular incidence with Cos Alpha="<<cosal<<std::endl;
if(rollId.station()==1&&rollId.layer()==1) { if(cluSize==1*dupli) {hGlobalResClu1La1->Fill(minres);}if(cluSize==2*dupli){ hGlobalResClu2La1->Fill(minres);} else if(cluSize==3*dupli){ hGlobalResClu3La1->Fill(minres);}}
else if(rollId.station()==1&&rollId.layer()==2){ if(cluSize==1*dupli) {hGlobalResClu1La2->Fill(minres);}if(cluSize==2*dupli){ hGlobalResClu2La2->Fill(minres);} else if(cluSize==3*dupli){ hGlobalResClu3La2->Fill(minres);}}
else if(rollId.station()==2&&rollId.layer()==1){ if(cluSize==1*dupli) {hGlobalResClu1La3->Fill(minres);}if(cluSize==2*dupli){ hGlobalResClu2La3->Fill(minres);} else if(cluSize==3*dupli){ hGlobalResClu3La3->Fill(minres);}}
else if(rollId.station()==2&&rollId.layer()==2){ if(cluSize==1*dupli) {hGlobalResClu1La4->Fill(minres);}if(cluSize==2*dupli){ hGlobalResClu2La4->Fill(minres);} else if(cluSize==3*dupli){ hGlobalResClu3La4->Fill(minres);}}
else if(rollId.station()==3) { if(cluSize==1*dupli) {hGlobalResClu1La5->Fill(minres);}if(cluSize==2*dupli){ hGlobalResClu2La5->Fill(minres);} else if(cluSize==3*dupli){ hGlobalResClu3La5->Fill(minres);}}
sprintf(meIdRPC,"RPCDataOccupancyFromDT_%d",rollId.rawId());
if(debug) std::cout<<"DT \t \t \t \t \t COINCIDENCE!!! Event="<<iEvent.id()<<" Filling RPC Data Occupancy for "<<meIdRPC<<" with "<<stripPredicted<<std::endl;
meMap[meIdRPC]->Fill(stripPredicted);
}
else{
RPCGeomServ rpcsrv(rollasociated->id());
std::string nameRoll = rpcsrv.name();
if(debug) std::cout<<"DT \t \t \t \t \t A roll was ineficient in event "<<iEvent.id().event()<<std::endl;
if(debug) ofrej<<"DTs \t Wh "<<dtWheel
<<"\t St "<<dtStation
<<"\t Se "<<dtSector
<<"\t Roll "<<nameRoll
<<"\t Event "
<<iEvent.id().event()
<<"\t Run "
<<iEvent.id().run()
<<std::endl;
}
}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){
if(debug) std::cout <<"MB4 \t Getting ALL the DT Segments"<<std::endl;
edm::Handle<DTRecSegment4DCollection> all4DSegments;
iEvent.getByLabel(dt4DSegments, all4DSegments);
iEvent.getByLabel(dt4DSegments, all4DSegments);
if(all4DSegments.isValid()){
if(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(debug) std::cout<<"MB4 \t \t This Segment is in Chamber id: "<<DTId<<std::endl;
if(debug) std::cout<<"MB4 \t \t Number of segments in this DT = "<<DTSegmentCounter[DTId]<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t Is the only one in this DT? and is in the Station 4?"<<std::endl;
if(DTSegmentCounter[DTId] == 1 && DTId.station()==4){
if(debug) std::cout<<"MB4 \t \t \t yes"<<std::endl;
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(debug) std::cout<<"MB4 \t \t \t \t The Segment in MB4 is 2D?"<<std::endl;
if(segment->dimension()==2){
if(debug) std::cout<<"MB4 \t \t \t \t yes"<<std::endl;
LocalVector segmentDirectionMB4=segmentDirection;
LocalPoint segmentPositionMB4=segmentPosition;
bool compatiblesegments=false;
const BoundPlane& DTSurface4 = dtGeo->idToDet(DTId)->surface();
DTRecSegment4DCollection::const_iterator segMB3;
if(debug) std::cout<<"MB4 \t \t \t \t Loop on segments in =sector && MB3 && adjacent sectors && y dim=4"<<std::endl;
for(segMB3=all4DSegments->begin();segMB3!=all4DSegments->end();++segMB3){
DTChamberId dtid3 = segMB3->chamberId();
if(distsector(dtid3.sector(),DTId.sector())<=1
&& dtid3.station()==3
&& dtid3.wheel()==DTId.wheel()
&& DTSegmentCounter[dtid3] == 1
&& segMB3->dimension()==4){
if(debug) std::cout<<"MB4 \t \t \t \t distsector ="<<distsector(dtid3.sector(),DTId.sector())<<std::endl;
const GeomDet* gdet3=dtGeo->idToDet(segMB3->geographicalId());
const BoundPlane & DTSurface3 = gdet3->surface();
LocalVector segmentDirectionMB3 = segMB3->localDirection();
GlobalPoint segmentPositionMB3inGlobal = DTSurface3.toGlobal(segMB3->localPosition());
LocalVector segDirMB4inMB3Frame=DTSurface3.toLocal(DTSurface4.toGlobal(segmentDirectionMB4));
LocalVector segDirMB3inMB4Frame=DTSurface4.toLocal(DTSurface3.toGlobal(segmentDirectionMB3));
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(debug) std::cout<<"MB4 \t \t \t \t cosAng"<<cosAng<<"Beetween "<<dtid3<<" and "<<DTId<<std::endl;
if(fabs(cosAng)>1.){
if(debug) std::cout<<"dx="<<dx<<" dz="<<dz<<std::endl;
if(debug) std::cout<<"dx3="<<dx3<<" dz3="<<dz<<std::endl;
if(debug) std::cout<<cosAng<<std::endl;
}
if(cosAng>MinCosAng){
compatiblesegments=true;
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
if(debug) std::cout<<"MB4 \t \t Number of rolls for this DT = "<<rollsForThisDT.size()<<std::endl;
assert(rollsForThisDT.size()>=1);
if(debug) std::cout<<"MB4 \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); //roll asociado a MB4
RPCDetId rpcId = rollasociated->id();
const BoundPlane & RPCSurfaceRB4 = rollasociated->surface(); //surface MB4
RPCGeomServ rpcsrv(rpcId);
std::string nameRoll = rpcsrv.name();
if(debug) std::cout<<"MB4 \t \t \t RollName: "<<nameRoll<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t Doing the extrapolation to this roll"<<std::endl;
GlobalPoint CenterPointRollGlobal=RPCSurfaceRB4.toGlobal(LocalPoint(0,0,0));
LocalPoint CenterRollinMB4Frame = DTSurface4.toLocal(CenterPointRollGlobal); //In MB4
LocalPoint segmentPositionMB3inMB4Frame = DTSurface4.toLocal(segmentPositionMB3inGlobal); //In MB4
LocalPoint segmentPositionMB3inRB4Frame = RPCSurfaceRB4.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
if(debug) std::cout<<"MB4 \t \t \t The distance to extrapolate in Y from MB3 is "<<Dy<<"cm"<<std::endl;
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();
if(debug) std::cout<<"MB4 \t \t \t Strip Lenght "<<stripl<<"cm"<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t Strip Width "<<stripw<<"cm"<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t X Predicted in MB4DTLocal= "<<X<<"cm"<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t Y Predicted in MB4DTLocal= "<<Y<<"cm"<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t Z Predicted in MB4DTLocal= "<<Z<<"cm"<<std::endl;
float extrapolatedDistance = sqrt((Y-Yo34)*(Y-Yo34)+Dy*Dy);
if(debug) std::cout<<"MB4 \t \t \t segmentPositionMB3inMB4Frame"<<segmentPositionMB3inMB4Frame<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t segmentPositionMB4inMB4Frame"<<segmentPosition<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t segmentDirMB3inMB4Frame"<<segDirMB3inMB4Frame<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t segmentDirMB4inMB4Frame"<<segmentDirectionMB4<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t CenterRB4PositioninMB4Frame"<<CenterRollinMB4Frame<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t Is the extrapolation distance ="<<extrapolatedDistance<<"less than "<<MaxDrb4<<std::endl;
if(extrapolatedDistance<=MaxDrb4){
if(debug) std::cout<<"MB4 \t \t \t yes"<<std::endl;
GlobalPoint GlobalPointExtrapolated = DTSurface4.toGlobal(LocalPoint(X,Y,Z));
if(debug) std::cout<<"MB4 \t \t \t Point ExtraPolated in Global"<<GlobalPointExtrapolated<< std::endl;
LocalPoint PointExtrapolatedRPCFrame = RPCSurfaceRB4.toLocal(GlobalPointExtrapolated);
if(debug) std::cout<<"MB4 \t \t \t Point Extrapolated in RPCLocal"<<PointExtrapolatedRPCFrame<< std::endl;
if(debug) std::cout<<"MB4 \t \t \t Corner of the Roll = ("<<rsize*0.5<<","<<stripl*0.5<<")"<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t Info About the Point Extrapolated in X Abs ("<<fabs(PointExtrapolatedRPCFrame.x())<<","
<<fabs(PointExtrapolatedRPCFrame.y())<<","<<fabs(PointExtrapolatedRPCFrame.z())<<")"<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t Does the extrapolation go inside this roll?"<<std::endl;
if(fabs(PointExtrapolatedRPCFrame.z()) < 5. &&
fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
if(debug) std::cout<<"MB4 \t \t \t \t yes"<<std::endl;
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()];
if(debug) std::cout<<"MB4 \t \t \t \t \t Filling Expected"<<std::endl;
sprintf(meIdDT,"ExpectedOccupancyFromDT_%d",rollId.rawId());
meMap[meIdDT]->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();
}
}
bool anycoincidence=false;
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(debug) std::cout<<"MB4 \t \t \t \t \t \t PointExtrapolatedRPCFrame.x="<<PointExtrapolatedRPCFrame.x()<<" Minimal Residual ="<<minres<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t \t \t \t Minimal Residual less than stripw*rangestrips? minres="<<minres<<" range="<<rangestrips<<" stripw="<<stripw<<" cluSize"<<cluSize<<" <=compare minres with"<<(rangestrips+cluSize*0.5)*stripw<<std::endl;
if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
if(debug) std::cout<<"MB4 \t \t \t \t \t \t \t True!"<<std::endl;
anycoincidence=true;
}
}
if(anycoincidence){
if(debug) std::cout<<"MB4 \t \t \t \t \t At least one RecHit inside the range, Predicted="<<stripPredicted<<" minres="<<minres<<"cm range="<<rangestrips<<"strips stripw="<<stripw<<"cm"<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t \t \t Norm of Cosine Directors="<<dx3*dx3+dy3*dy3+dz3*dz3<<"~1?"<<std::endl;
float cosal = dx/sqrt(dx*dx+dz*dz);
if(debug) std::cout<<"MB4 \t \t \t \t \t Angle="<<acos(cosal)*180/3.1415926<<" degree"<<std::endl;
if(debug) std::cout<<"MB4 \t \t \t \t \t Filling the Residuals Histogram for globals with "<<minres<<"And the angular incidence with Cos Theta="<<-1*dz<<std::endl;
assert(rollId.station()==4);
if(cluSize==1*dupli){ hGlobalResClu1La6->Fill(minres);}
else if(cluSize==2*dupli){ hGlobalResClu2La6->Fill(minres);}
else if(cluSize==3*dupli){ hGlobalResClu3La6->Fill(minres);}
sprintf(meIdRPC,"RPCDataOccupancyFromDT_%d",rollId.rawId());
if(debug) std::cout<<"MB4 \t \t \t \t \t \t COINCIDENCE!!! Event="<<iEvent.id()<<"Filling RPC Data Occupancy for "<<meIdRPC<<" with "<<stripPredicted<<std::endl;
meMap[meIdRPC]->Fill(stripPredicted);
}
else{
RPCGeomServ rpcsrv(rollasociated->id());
std::string nameRoll = rpcsrv.name();
if(debug) std::cout<<"MB4 \t \t \t \t \t \t A roll was ineficient in event"<<iEvent.id().event()<<std::endl;
if(debug) ofrej<<"MB4 \t Wh "<<dtWheel
<<"\t St "<<dtStation
<<"\t Se "<<dtSector
<<"\t Roll "<<nameRoll
<<"\t Event "
<<iEvent.id().event()
<<"\t Run "
<<iEvent.id().run()
<<std::endl;
}
}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{
compatiblesegments=false;
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(debug) std::cout<<"CSC \t \t This Segment is in Chamber id: "<<CSCId<<std::endl;
if(debug) std::cout<<"CSC \t \t Number of segments in this CSC = "<<CSCSegmentsCounter[CSCId]<<std::endl;
if(debug) std::cout<<"CSC \t \t Is the only one in this CSC? is not ind the ring 1 or station 4? Are there more than 2 segments in the event?"<<std::endl;
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. ){
//&& segment->chi2()< ??)Add 3 segmentes in the endcaps???
if(debug) std::cout<<"CSC \t \t yes"<<std::endl;
if(debug) std::cout<<"CSC \t \t CSC Segment Dimension "<<segment->dimension()<<std::endl;
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(debug) std::cout<<"CSC \t \t Loop over all the rolls asociated to this CSC"<<std::endl;
if(rpcRing!=1&&rpcStation!=4){
if(rollsForThisCSC.size()==0){
if(debug) std::cout<<"CSC Fail for CSCId="<<TheId<<" rpcRegion="<<rpcRegion<<" rpcStation="<<rpcStation<<" rpcRing="<<rpcRing<<" rpcSegment="<<rpcSegment<<std::endl;
}
assert(rollsForThisCSC.size()>=1);
//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();
if(debug) std::cout<<"CSC \t \t \t We are in the roll getting the surface"<<rpcId<<std::endl;
const BoundPlane & RPCSurface = rollasociated->surface();
if(debug) std::cout<<"CSC \t \t \t RollID: "<<rpcId<<std::endl;
if(debug) std::cout<<"CSC \t \t \t Doing the extrapolation to this roll"<<std::endl;
if(debug) std::cout<<"CSC \t \t \t CSC Segment Direction in CSCLocal "<<segmentDirection<<std::endl;
if(debug) std::cout<<"CSC \t \t \t CSC Segment Point in CSCLocal "<<segmentPosition<<std::endl;
GlobalPoint CenterPointRollGlobal = RPCSurface.toGlobal(LocalPoint(0,0,0));
if(debug) std::cout<<"CSC \t \t \t Center (0,0,0) of the Roll in Global"<<CenterPointRollGlobal<<std::endl;
GlobalPoint CenterPointCSCGlobal = TheChamber->toGlobal(LocalPoint(0,0,0));
if(debug) std::cout<<"CSC \t \t \t Center (0,0,0) of the CSC in Global"<<CenterPointCSCGlobal<<std::endl;
GlobalPoint segmentPositionInGlobal=TheChamber->toGlobal(segmentPosition); //new way to convert to global
if(debug) std::cout<<"CSC \t \t \t Segment Position in Global"<<segmentPositionInGlobal<<std::endl;
LocalPoint CenterRollinCSCFrame = TheChamber->toLocal(CenterPointRollGlobal);
if(debug){//to check CSC RPC phi relation!
float rpcphi=0;
float cscphi=0;
(CenterPointRollGlobal.barePhi()<0)?
rpcphi = 2*3.141592+CenterPointRollGlobal.barePhi():rpcphi=CenterPointRollGlobal.barePhi();
(CenterPointCSCGlobal.barePhi()<0)?
cscphi = 2*3.1415926536+CenterPointCSCGlobal.barePhi():cscphi=CenterPointCSCGlobal.barePhi();
float df=fabs(cscphi-rpcphi);
float dr=fabs(CenterPointRollGlobal.perp()-CenterPointCSCGlobal.perp());
float diffz=CenterPointRollGlobal.z()-CenterPointCSCGlobal.z();
float dfg=df*180./3.14159265;
if(debug) std::cout<<"CSC \t \t \t z of RPC="<<CenterPointRollGlobal.z()<<"z of CSC"<<CenterPointCSCGlobal.z()<<" dfg="<<dfg<<std::endl;
RPCGeomServ rpcsrv(rpcId);
if(dr>200.||fabs(dz)>55.||dfg>1.){
//if(rpcRegion==1&&dfg>1.&&dr>100.){
if (debug) std::cout
<<"\t \t \t CSC Station= "<<CSCId.station()
<<" Ring= "<<CSCId.ring()
<<" Chamber= "<<CSCId.chamber()
<<" cscphi="<<cscphi*180/3.14159265
<<"\t RPC Station= "<<rpcId.station()
<<" ring= "<<rpcId.ring()
<<" segment =-> "<<rpcsrv.name()
<<" rollphi="<<rpcphi*180/3.14159265
<<"\t dfg="<<dfg
<<" dz="<<diffz
<<" dr="<<dr
<<std::endl;
}
}
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();
if(debug) std::cout<<"CSC \t \t \t Strip Lenght "<<stripl<<"cm"<<std::endl;
if(debug) std::cout<<"CSC \t \t \t Strip Width "<<stripw<<"cm"<<std::endl;
if(debug) std::cout<<"CSC \t \t \t X Predicted in CSCLocal= "<<X<<"cm"<<std::endl;
if(debug) std::cout<<"CSC \t \t \t Y Predicted in CSCLocal= "<<Y<<"cm"<<std::endl;
if(debug) std::cout<<"CSC \t \t \t Z Predicted in CSCLocal= "<<Z<<"cm"<<std::endl;
float extrapolatedDistance = sqrt((X-Xo)*(X-Xo)+(Y-Yo)*(Y-Yo)+(Z-Zo)*(Z-Zo));
if(debug) std::cout<<"CSC \t \t \t Is the distance of extrapolation less than MaxD? ="<<extrapolatedDistance<<"cm"<<"MaxD="<<MaxD<<"cm"<<std::endl;
if(extrapolatedDistance<=MaxD){
if(debug) std::cout<<"CSC \t \t \t yes"<<std::endl;
GlobalPoint GlobalPointExtrapolated=TheChamber->toGlobal(LocalPoint(X,Y,Z));
if(debug) std::cout<<"CSC \t \t \t Point ExtraPolated in Global"<<GlobalPointExtrapolated<< std::endl;
LocalPoint PointExtrapolatedRPCFrame = RPCSurface.toLocal(GlobalPointExtrapolated);
if(debug) std::cout<<"CSC \t \t \t Point Extrapolated in RPCLocal"<<PointExtrapolatedRPCFrame<< std::endl;
if(debug) std::cout<<"CSC \t \t \t Corner of the Roll = ("<<rsize*0.5<<","<<stripl*0.5<<")"<<std::endl;
if(debug) std::cout<<"CSC \t \t \t Info About the Point Extrapolated in X Abs ("<<fabs(PointExtrapolatedRPCFrame.x())<<","
<<fabs(PointExtrapolatedRPCFrame.y())<<","<<fabs(PointExtrapolatedRPCFrame.z())<<")"<<std::endl;
if(debug) std::cout<<"CSC \t \t \t dz="
<<fabs(PointExtrapolatedRPCFrame.z())<<" dx="
<<fabs(PointExtrapolatedRPCFrame.x())<<" dy="
<<fabs(PointExtrapolatedRPCFrame.y())<<std::endl;
if(debug) std::cout<<"CSC \t \t \t Does the extrapolation go inside this roll????"<<std::endl;
if(fabs(PointExtrapolatedRPCFrame.z()) < 10. &&
fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
if(debug) std::cout<<"CSC \t \t \t \t yes"<<std::endl;
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()];
if(debug) std::cout<<"CSC \t \t \t \t Filling Expected"<<std::endl;
sprintf(meIdCSC,"ExpectedOccupancyFromCSC_%d",rollId.rawId());
meMap[meIdCSC]->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++) {
sprintf(meIdRPC,"BXDistribution_%d",rollasociated->id().rawId());
meMap[meIdRPC]->Fill(recHit->BunchX());
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;
}
}
bool anycoincidence = false;
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(debug) std::cout<<"CSC \t \t \t \t \t PointExtrapolatedRPCFrame.x="<<PointExtrapolatedRPCFrame.x()<<" Minimal Residual"<<minres<<std::endl;
if(debug) std::cout<<"CSC \t \t \t \t \t Minimal Residual less than stripw*rangestrips? minres="<<minres<<" range="<<rangestrips<<" stripw="<<stripw<<" cluSize"<<cluSize<<" <=compare minres with"<<(rangestrips+cluSize*0.5)*stripw<<std::endl;
if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
if(debug) std::cout<<"CSC \t \t \t \t \t \t True!"<<std::endl;
anycoincidence=true;
}
}
if(anycoincidence){
if(debug) std::cout<<"CSC \t \t \t \t \t At least one RecHit inside the range, Predicted="<<stripPredicted<<" minres="<<minres<<"cm range="<<rangestrips<<"strips stripw="<<stripw<<"cm"<<std::endl;
if(debug) std::cout<<"CSC \t \t \t \t \t Norm of Cosine Directors="<<dx*dx+dy*dy+dz*dz<<"~1?"<<std::endl;
float cosal = dx/sqrt(dx*dx+dz*dz);
if(debug) std::cout<<"CSC \t \t \t \t \t Angle="<<acos(cosal)*180/3.1415926<<" degree"<<std::endl;
if(debug) std::cout<<"CSC \t \t \t \t \t Filling the Residuals Histogram for globals with "<<minres<<"And the angular incidence with Cos Theta="<<-1*dz<<std::endl;
if(rollId.ring()==2&&rollId.roll()==1){if(cluSize==1*dupli) hGlobalResClu1R2A->Fill(minres); if(cluSize==2*dupli) hGlobalResClu2R2A->Fill(minres); if(cluSize==3*dupli) hGlobalResClu3R2A->Fill(minres);}
if(rollId.ring()==2&&rollId.roll()==2){if(cluSize==1*dupli) hGlobalResClu1R2B->Fill(minres); if(cluSize==2*dupli) hGlobalResClu2R2B->Fill(minres); if(cluSize==3*dupli) hGlobalResClu3R2B->Fill(minres);}
if(rollId.ring()==2&&rollId.roll()==3){if(cluSize==1*dupli) hGlobalResClu1R2C->Fill(minres); if(cluSize==2*dupli) hGlobalResClu2R2C->Fill(minres); if(cluSize==3*dupli) hGlobalResClu3R2C->Fill(minres);}
if(rollId.ring()==3&&rollId.roll()==1){if(cluSize==1*dupli) hGlobalResClu1R3A->Fill(minres); if(cluSize==2*dupli) hGlobalResClu2R3A->Fill(minres); if(cluSize==3*dupli) hGlobalResClu3R3A->Fill(minres);}
if(rollId.ring()==3&&rollId.roll()==2){if(cluSize==1*dupli) hGlobalResClu1R3B->Fill(minres); if(cluSize==2*dupli) hGlobalResClu2R3B->Fill(minres); if(cluSize==3*dupli) hGlobalResClu3R3B->Fill(minres);}
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);}
sprintf(meIdRPC,"RPCDataOccupancyFromCSC_%d",rollId.rawId());
if(debug) std::cout <<"CSC \t \t \t \t \t \t COINCEDENCE!!! Event="<<iEvent.id()<<"Filling RPC Data Occupancy for "<<meIdRPC<<" with "<<stripPredicted<<std::endl;
meMap[meIdRPC]->Fill(stripPredicted);
}
else{
RPCGeomServ rpcsrv(rollasociated->id());
std::string nameRoll = rpcsrv.name();
if(debug) std::cout<<"CSC \t \t \t \t \t \t A roll was ineficient in event"<<iEvent.id().event()<<std::endl;
if(debug) ofrej<<"CSC \t EndCap "<<rpcRegion
<<"\t cscStation "<<cscStation
<<"\t cscRing "<<cscRing
<<"\t cscChamber "<<cscChamber
<<"\t Roll "<<nameRoll
<<"\t Event "<<iEvent.id().event()
<<"\t CSCId "<<CSCId
<<"\t Event "
<<iEvent.id().event()
<<"\t Run "
<<iEvent.id().run()
<<std::endl;
}
}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 140 of file RPCEfficiency.cc.
References bookDetUnitSeg(), gather_cfg::cout, cscGeo, debug, dtGeo, folderPath, edm::EventSetup::get(), inclcsc, incldt, incldtMB4, meCollection, RPCGeomServ::name(), csvReporter::r, DetId::rawId(), RPCDetId::region(), RPCDetId::ring(), relativeConstraints::ring, RPCChamber::rolls(), rollstoreCSC, rollstoreDT, rpcGeo, RPCDetId::sector(), RPCGeomServ::segment(), relativeConstraints::station, and RPCDetId::station().
{
iSetup.get<MuonGeometryRecord>().get(rpcGeo);
iSetup.get<MuonGeometryRecord>().get(dtGeo);
iSetup.get<MuonGeometryRecord>().get(cscGeo);
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
RPCGeomServ rpcsrv(rpcId);
std::string nameRoll = rpcsrv.name();
if(debug) std::cout<<"Booking for "<<rpcId.rawId()<<std::endl;
meCollection[rpcId.rawId()] = bookDetUnitSeg(rpcId,(*r)->nstrips(),folderPath+"MuonSegEff/");
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;
}
if(region!=0 && inclcsc){
int region=rpcId.region();
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&&(incldt||incldtMB4)&&rpcId.ring()!=0&&rpcId.station()!=4){
//std::cout<<"--Filling the dtstore for statistics"<<rpcId<<std::endl;
int sidewheel = 0;
if(rpcId.ring()==-2){
sidewheel=-1;
}
else if(rpcId.ring()==-1){
sidewheel=0;
}
else if(rpcId.ring()==1){
sidewheel=0;
}
else if(rpcId.ring()==2){
sidewheel=1;
}
int wheel= sidewheel;
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;
}*/
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
| std::map< std::string, MonitorElement * > RPCEfficiency::bookDetUnitSeg | ( | RPCDetId & | detId, |
| int | nstrips, | ||
| std::string | folder | ||
| ) |
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 1166 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 1159 of file RPCEfficiency.cc.
References dbe, EffRootFileName, EffSaveRootFile, and DQMStore::save().
{
if (EffSaveRootFile){
dbe->save(EffRootFileName);
}
}
Definition at line 145 of file RPCEfficiency.h.
Definition at line 146 of file RPCEfficiency.h.
Definition at line 147 of file RPCEfficiency.h.
Definition at line 148 of file RPCEfficiency.h.
Definition at line 149 of file RPCEfficiency.h.
Definition at line 150 of file RPCEfficiency.h.
Definition at line 152 of file RPCEfficiency.h.
Definition at line 153 of file RPCEfficiency.h.
Definition at line 154 of file RPCEfficiency.h.
Definition at line 155 of file RPCEfficiency.h.
Definition at line 156 of file RPCEfficiency.h.
Definition at line 157 of file RPCEfficiency.h.
Definition at line 159 of file RPCEfficiency.h.
Definition at line 160 of file RPCEfficiency.h.
Definition at line 161 of file RPCEfficiency.h.
Definition at line 162 of file RPCEfficiency.h.
Definition at line 163 of file RPCEfficiency.h.
Definition at line 164 of file RPCEfficiency.h.
std::vector<std::map<RPCDetId, int> > RPCEfficiency::counter [private] |
Definition at line 212 of file RPCEfficiency.h.
Definition at line 118 of file RPCEfficiency.h.
Referenced by analyze(), and beginRun().
std::string RPCEfficiency::cscSegments [private] |
Definition at line 228 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
DQMStore* RPCEfficiency::dbe [private] |
Definition at line 240 of file RPCEfficiency.h.
Referenced by bookDetUnitSeg(), endJob(), endRun(), and RPCEfficiency().
bool RPCEfficiency::debug [private] |
Definition at line 218 of file RPCEfficiency.h.
Referenced by analyze(), beginRun(), and RPCEfficiency().
Definition at line 123 of file RPCEfficiency.h.
Definition at line 124 of file RPCEfficiency.h.
Definition at line 125 of file RPCEfficiency.h.
Definition at line 126 of file RPCEfficiency.h.
Definition at line 127 of file RPCEfficiency.h.
Definition at line 128 of file RPCEfficiency.h.
Definition at line 130 of file RPCEfficiency.h.
Definition at line 131 of file RPCEfficiency.h.
Definition at line 132 of file RPCEfficiency.h.
Definition at line 133 of file RPCEfficiency.h.
Definition at line 134 of file RPCEfficiency.h.
Definition at line 135 of file RPCEfficiency.h.
Definition at line 137 of file RPCEfficiency.h.
Definition at line 138 of file RPCEfficiency.h.
Definition at line 139 of file RPCEfficiency.h.
Definition at line 140 of file RPCEfficiency.h.
Definition at line 141 of file RPCEfficiency.h.
Definition at line 142 of file RPCEfficiency.h.
std::string RPCEfficiency::dt4DSegments [private] |
Definition at line 229 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 117 of file RPCEfficiency.h.
Referenced by analyze(), and beginRun().
int RPCEfficiency::dupli [private] |
Definition at line 226 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
std::string RPCEfficiency::EffRootFileName [private] |
Definition at line 237 of file RPCEfficiency.h.
Referenced by endRun(), and RPCEfficiency().
bool RPCEfficiency::EffSaveRootFile [private] |
Definition at line 236 of file RPCEfficiency.h.
Referenced by endRun(), and RPCEfficiency().
std::string RPCEfficiency::folderPath [private] |
Definition at line 230 of file RPCEfficiency.h.
Referenced by beginRun(), and RPCEfficiency().
Definition at line 168 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 169 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 170 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 171 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 172 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 173 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 195 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 194 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 193 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 192 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 191 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 190 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 175 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 176 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 177 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 178 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 179 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 180 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 202 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 201 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 200 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 199 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 198 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 197 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 182 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 183 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 184 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 185 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 186 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 187 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 209 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 208 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 207 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 206 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 205 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 204 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
bool RPCEfficiency::inclcsc [private] |
Definition at line 217 of file RPCEfficiency.h.
Referenced by analyze(), beginRun(), and RPCEfficiency().
bool RPCEfficiency::incldt [private] |
Definition at line 215 of file RPCEfficiency.h.
Referenced by analyze(), beginRun(), and RPCEfficiency().
bool RPCEfficiency::incldtMB4 [private] |
Definition at line 216 of file RPCEfficiency.h.
Referenced by analyze(), beginRun(), and RPCEfficiency().
bool RPCEfficiency::inves [private] |
Definition at line 220 of file RPCEfficiency.h.
double RPCEfficiency::MaxD [private] |
Definition at line 224 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
double RPCEfficiency::MaxDrb4 [private] |
Definition at line 225 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
std::map<int, std::map<std::string, MonitorElement*> > RPCEfficiency::meCollection [private] |
Definition at line 234 of file RPCEfficiency.h.
Referenced by analyze(), and beginRun().
double RPCEfficiency::MinCosAng [private] |
Definition at line 223 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
std::string RPCEfficiency::muonRPCDigis [private] |
Definition at line 227 of file RPCEfficiency.h.
Referenced by RPCEfficiency().
std::string RPCEfficiency::nameInLog [private] |
Definition at line 238 of file RPCEfficiency.h.
Referenced by RPCEfficiency().
std::ofstream RPCEfficiency::ofrej [private] |
Definition at line 214 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
bool RPCEfficiency::paper [private] |
Definition at line 219 of file RPCEfficiency.h.
double RPCEfficiency::rangestrips [private] |
Definition at line 221 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
double RPCEfficiency::rangestripsRB4 [private] |
Definition at line 222 of file RPCEfficiency.h.
Referenced by RPCEfficiency().
std::string RPCEfficiency::rejected [private] |
Definition at line 231 of file RPCEfficiency.h.
std::string RPCEfficiency::rollseff [private] |
Definition at line 232 of file RPCEfficiency.h.
| std::map<CSCStationIndex,std::set<RPCDetId> > RPCEfficiency::rollstoreCSC |
Definition at line 115 of file RPCEfficiency.h.
Referenced by analyze(), and beginRun().
| std::map<DTStationIndex,std::set<RPCDetId> > RPCEfficiency::rollstoreDT |
Definition at line 114 of file RPCEfficiency.h.
Referenced by analyze(), and beginRun().
Definition at line 116 of file RPCEfficiency.h.
Referenced by analyze(), and beginRun().
Definition at line 121 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
Definition at line 120 of file RPCEfficiency.h.
Referenced by analyze(), and RPCEfficiency().
std::vector<int> RPCEfficiency::totalcounter [private] |
Definition at line 213 of file RPCEfficiency.h.
1.7.3