00001
00002
00003
00004
00005
00006
00007
00008 #include "DQMOffline/Muon/interface/RPCEfficiency.h"
00009 #include <memory>
00010 #include "FWCore/Framework/interface/MakerMacros.h"
00011 #include <DataFormats/RPCDigi/interface/RPCDigiCollection.h>
00012 #include "DataFormats/RPCRecHit/interface/RPCRecHitCollection.h"
00013 #include <DataFormats/MuonDetId/interface/RPCDetId.h>
00014 #include <DataFormats/DTRecHit/interface/DTRecSegment4DCollection.h>
00015 #include <DataFormats/CSCRecHit/interface/CSCSegmentCollection.h>
00016 #include <Geometry/RPCGeometry/interface/RPCGeomServ.h>
00017 #include <Geometry/CommonDetUnit/interface/GeomDet.h>
00018 #include <Geometry/Records/interface/MuonGeometryRecord.h>
00019 #include <Geometry/CommonTopologies/interface/RectangularStripTopology.h>
00020 #include <Geometry/CommonTopologies/interface/TrapezoidalStripTopology.h>
00021
00022 #include <cmath>
00023 #include "TFile.h"
00024 #include "TH1F.h"
00025 #include "TH2F.h"
00026 #include "TCanvas.h"
00027 #include "TAxis.h"
00028 #include "TString.h"
00029
00030 int mySegment(RPCDetId rpcId){
00031 int seg=0;
00032 int nsec=36;
00033 int nsub=6;
00034 if (rpcId.ring()==1 && rpcId.station() > 1) {
00035 nsub=3;
00036 nsec=18;
00037 }
00038 seg =rpcId.subsector()+nsub*(rpcId.sector()-1);
00039 if(seg==nsec+1)seg=1;
00040 return seg;
00041 }
00042
00043 void RPCEfficiency::beginJob(){
00044
00045 }
00046
00047 RPCEfficiency::RPCEfficiency(const edm::ParameterSet& iConfig){
00048 incldt=iConfig.getUntrackedParameter<bool>("incldt",true);
00049 incldtMB4=iConfig.getUntrackedParameter<bool>("incldtMB4",true);
00050 inclcsc=iConfig.getUntrackedParameter<bool>("inclcsc",true);
00051 debug=iConfig.getUntrackedParameter<bool>("debug",false);
00052 inves=iConfig.getUntrackedParameter<bool>("inves");
00053
00054 rangestrips = iConfig.getUntrackedParameter<double>("rangestrips",1.);
00055 rangestripsRB4=iConfig.getUntrackedParameter<double>("rangestripsRB4",4.);
00056 dupli = iConfig.getUntrackedParameter<int>("DuplicationCorrection",2);
00057 MinCosAng=iConfig.getUntrackedParameter<double>("MinCosAng",0.9999);
00058 MaxD=iConfig.getUntrackedParameter<double>("MaxD",20.);
00059 MaxDrb4=iConfig.getUntrackedParameter<double>("MaxDrb4",30.);
00060 muonRPCDigis=iConfig.getUntrackedParameter<std::string>("muonRPCDigis","muonRPCDigis");
00061 cscSegments=iConfig.getUntrackedParameter<std::string>("cscSegments","cscSegments");
00062 dt4DSegments=iConfig.getUntrackedParameter<std::string>("dt4DSegments","dt4DSegments");
00063
00064 nameInLog = iConfig.getUntrackedParameter<std::string>("moduleLogName", "RPC_Eff");
00065 EffSaveRootFile = iConfig.getUntrackedParameter<bool>("EffSaveRootFile", false);
00066 EffRootFileName = iConfig.getUntrackedParameter<std::string>("EffRootFileName", "RPCEfficiency.root");
00067
00068
00069
00070 dbe = edm::Service<DQMStore>().operator->();
00071
00072 std::string folder = "Muons/MuonSegEff/";
00073 dbe->setCurrentFolder(folder);
00074 statistics = dbe->book1D("Statistics","All Statistics",33,0.5,33.5);
00075
00076 if(debug) std::cout<<"booking Global histograms"<<std::endl;
00077
00078 statistics->setBinLabel(1,"Events ",1);
00079 statistics->setBinLabel(2,"Events with DT segments",1);
00080 statistics->setBinLabel(3,"Events with 1 DT segment",1);
00081 statistics->setBinLabel(4,"Events with 2 DT segments",1);
00082 statistics->setBinLabel(5,"Events with 3 DT segments",1);
00083 statistics->setBinLabel(6,"Events with 4 DT segments",1);
00084 statistics->setBinLabel(7,"Events with 5 DT segments",1);
00085 statistics->setBinLabel(8,"Events with 6 DT segments",1);
00086 statistics->setBinLabel(9,"Events with 7 DT segments",1);
00087 statistics->setBinLabel(10,"Events with 8 DT segments",1);
00088 statistics->setBinLabel(11,"Events with 9 DT segments",1);
00089 statistics->setBinLabel(12,"Events with 10 DT segments",1);
00090 statistics->setBinLabel(13,"Events with 11 DT segments",1);
00091 statistics->setBinLabel(14,"Events with 12 DT segments",1);
00092 statistics->setBinLabel(15,"Events with 13 DT segments",1);
00093 statistics->setBinLabel(16,"Events with 14 DT segments",1);
00094 statistics->setBinLabel(17,"Events with 15 DT segments",1);
00095 statistics->setBinLabel(18,"Events with CSC segments",1);
00096 statistics->setBinLabel(16+3,"Events with 1 CSC segment",1);
00097 statistics->setBinLabel(16+4,"Events with 2 CSC segments",1);
00098 statistics->setBinLabel(16+5,"Events with 3 CSC segments",1);
00099 statistics->setBinLabel(16+6,"Events with 4 CSC segments",1);
00100 statistics->setBinLabel(16+7,"Events with 5 CSC segments",1);
00101 statistics->setBinLabel(16+8,"Events with 6 CSC segments",1);
00102 statistics->setBinLabel(16+9,"Events with 7 CSC segments",1);
00103 statistics->setBinLabel(16+10,"Events with 8 CSC segments",1);
00104 statistics->setBinLabel(16+11,"Events with 9 CSC segments",1);
00105 statistics->setBinLabel(16+12,"Events with 10 CSC segments",1);
00106 statistics->setBinLabel(16+13,"Events with 11 CSC segments",1);
00107 statistics->setBinLabel(16+14,"Events with 12 CSC segments",1);
00108 statistics->setBinLabel(16+15,"Events with 13 CSC segments",1);
00109 statistics->setBinLabel(16+16,"Events with 14 CSC segments",1);
00110 statistics->setBinLabel(16+17,"Events with 15 CSC segments",1);
00111
00112 if(debug) std::cout<<"booking Global histograms Change statistics"<<std::endl;
00113
00114 folder = "Muons/MuonSegEff/Residuals/Investigation";
00115 dbe->setCurrentFolder(folder);
00116
00117
00118
00119 DistBorderClu1La1 = dbe->book1D("DistBorderClu1La1","Distance to the Border of the Strip Layer 1 Cluster Size 1",50,-2.,3.);
00120 DistBorderClu1La2 = dbe->book1D("DistBorderClu1La2","Distance to the Border of the Strip Layer 2 Cluster Size 1",50,-2.,3.);
00121 DistBorderClu1La3 = dbe->book1D("DistBorderClu1La3","Distance to the Border of the Strip Layer 3 Cluster Size 1",50,-2.,3.);
00122 DistBorderClu1La4 = dbe->book1D("DistBorderClu1La4","Distance to the Border of the Strip Layer 4 Cluster Size 1",50,-2.,3.);
00123 DistBorderClu1La5 = dbe->book1D("DistBorderClu1La5","Distance to the Border of the Strip Layer 5 Cluster Size 1",50,-2.,3.);
00124 DistBorderClu1La6 = dbe->book1D("DistBorderClu1La6","Distance to the Border of the Strip Layer 6 Cluster Size 1",50,-2.,3.);
00125
00126 DistBorderClu2La1 = dbe->book1D("DistBorderClu2La1","Distance to the Border of the Strip Layer 1 Cluster Size 2",50,-2.,3.);
00127 DistBorderClu2La2 = dbe->book1D("DistBorderClu2La2","Distance to the Border of the Strip Layer 2 Cluster Size 2",50,-2.,3.);
00128 DistBorderClu2La3 = dbe->book1D("DistBorderClu2La3","Distance to the Border of the Strip Layer 3 Cluster Size 2",50,-2.,3.);
00129 DistBorderClu2La4 = dbe->book1D("DistBorderClu2La4","Distance to the Border of the Strip Layer 4 Cluster Size 2",50,-2.,3.);
00130 DistBorderClu2La5 = dbe->book1D("DistBorderClu2La5","Distance to the Border of the Strip Layer 5 Cluster Size 2",50,-2.,3.);
00131 DistBorderClu2La6 = dbe->book1D("DistBorderClu2La6","Distance to the Border of the Strip Layer 6 Cluster Size 2",50,-2.,3.);
00132
00133 DistBorderClu3La1 = dbe->book1D("DistBorderClu3La1","Distance to the Border of the Strip Layer 1 Cluster Size 3",50,-2.,3.);
00134 DistBorderClu3La2 = dbe->book1D("DistBorderClu3La2","Distance to the Border of the Strip Layer 2 Cluster Size 3",50,-2.,3.);
00135 DistBorderClu3La3 = dbe->book1D("DistBorderClu3La3","Distance to the Border of the Strip Layer 3 Cluster Size 3",50,-2.,3.);
00136 DistBorderClu3La4 = dbe->book1D("DistBorderClu3La4","Distance to the Border of the Strip Layer 4 Cluster Size 3",50,-2.,3.);
00137 DistBorderClu3La5 = dbe->book1D("DistBorderClu3La5","Distance to the Border of the Strip Layer 5 Cluster Size 3",50,-2.,3.);
00138 DistBorderClu3La6 = dbe->book1D("DistBorderClu3La6","Distance to the Border of the Strip Layer 6 Cluster Size 3",50,-2.,3.);
00139
00140
00141 float pi = 3.14159265;
00142 AngClu1La1 = dbe->book1D("AngClu1La1","Angle of incident Muon Layer 1 Cluster Size 1",50,0.,pi);
00143 AngClu1La2 = dbe->book1D("AngClu1La2","Angle of incident Muon Layer 2 Cluster Size 1",50,0.,pi);
00144 AngClu1La3 = dbe->book1D("AngClu1La3","Angle of incident Muon Layer 3 Cluster Size 1",50,0.,pi);
00145 AngClu1La4 = dbe->book1D("AngClu1La4","Angle of incident Muon Layer 4 Cluster Size 1",50,0.,pi);
00146 AngClu1La5 = dbe->book1D("AngClu1La5","Angle of incident Muon Layer 5 Cluster Size 1",50,0.,pi);
00147 AngClu1La6 = dbe->book1D("AngClu1La6","Angle of incident Muon Layer 6 Cluster Size 1",50,0.,pi);
00148
00149 AngClu2La1 = dbe->book1D("AngClu2La1","Angle of incident Muon Layer 1 Cluster Size 2",50,0.,pi);
00150 AngClu2La2 = dbe->book1D("AngClu2La2","Angle of incident Muon Layer 2 Cluster Size 2",50,0.,pi);
00151 AngClu2La3 = dbe->book1D("AngClu2La3","Angle of incident Muon Layer 3 Cluster Size 2",50,0.,pi);
00152 AngClu2La4 = dbe->book1D("AngClu2La4","Angle of incident Muon Layer 4 Cluster Size 2",50,0.,pi);
00153 AngClu2La5 = dbe->book1D("AngClu2La5","Angle of incident Muon Layer 5 Cluster Size 2",50,0.,pi);
00154 AngClu2La6 = dbe->book1D("AngClu2La6","Angle of incident Muon Layer 6 Cluster Size 2",50,0.,pi);
00155
00156 AngClu3La1 = dbe->book1D("AngClu3La1","Angle of incident Muon Layer 1 Cluster Size 3",50,0.,pi);
00157 AngClu3La2 = dbe->book1D("AngClu3La2","Angle of incident Muon Layer 2 Cluster Size 3",50,0.,pi);
00158 AngClu3La3 = dbe->book1D("AngClu3La3","Angle of incident Muon Layer 3 Cluster Size 3",50,0.,pi);
00159 AngClu3La4 = dbe->book1D("AngClu3La4","Angle of incident Muon Layer 4 Cluster Size 3",50,0.,pi);
00160 AngClu3La5 = dbe->book1D("AngClu3La5","Angle of incident Muon Layer 5 Cluster Size 3",50,0.,pi);
00161 AngClu3La6 = dbe->book1D("AngClu3La6","Angle of incident Muon Layer 6 Cluster Size 3",50,0.,pi);
00162
00163 folder = "Muons/MuonSegEff/Residuals/Barrel";
00164 dbe->setCurrentFolder(folder);
00165
00166
00167 hGlobalResClu1La1 = dbe->book1D("GlobalResidualsClu1La1","RPC Residuals Layer 1 Cluster Size 1",101,-10.,10.);
00168 hGlobalResClu1La2 = dbe->book1D("GlobalResidualsClu1La2","RPC Residuals Layer 2 Cluster Size 1",101,-10.,10.);
00169 hGlobalResClu1La3 = dbe->book1D("GlobalResidualsClu1La3","RPC Residuals Layer 3 Cluster Size 1",101,-10.,10.);
00170 hGlobalResClu1La4 = dbe->book1D("GlobalResidualsClu1La4","RPC Residuals Layer 4 Cluster Size 1",101,-10.,10.);
00171 hGlobalResClu1La5 = dbe->book1D("GlobalResidualsClu1La5","RPC Residuals Layer 5 Cluster Size 1",101,-10.,10.);
00172 hGlobalResClu1La6 = dbe->book1D("GlobalResidualsClu1La6","RPC Residuals Layer 6 Cluster Size 1",101,-10.,10.);
00173
00174 hGlobalResClu2La1 = dbe->book1D("GlobalResidualsClu2La1","RPC Residuals Layer 1 Cluster Size 2",101,-10.,10.);
00175 hGlobalResClu2La2 = dbe->book1D("GlobalResidualsClu2La2","RPC Residuals Layer 2 Cluster Size 2",101,-10.,10.);
00176 hGlobalResClu2La3 = dbe->book1D("GlobalResidualsClu2La3","RPC Residuals Layer 3 Cluster Size 2",101,-10.,10.);
00177 hGlobalResClu2La4 = dbe->book1D("GlobalResidualsClu2La4","RPC Residuals Layer 4 Cluster Size 2",101,-10.,10.);
00178 hGlobalResClu2La5 = dbe->book1D("GlobalResidualsClu2La5","RPC Residuals Layer 5 Cluster Size 2",101,-10.,10.);
00179 hGlobalResClu2La6 = dbe->book1D("GlobalResidualsClu2La6","RPC Residuals Layer 6 Cluster Size 2",101,-10.,10.);
00180
00181 hGlobalResClu3La1 = dbe->book1D("GlobalResidualsClu3La1","RPC Residuals Layer 1 Cluster Size 3",101,-10.,10.);
00182 hGlobalResClu3La2 = dbe->book1D("GlobalResidualsClu3La2","RPC Residuals Layer 2 Cluster Size 3",101,-10.,10.);
00183 hGlobalResClu3La3 = dbe->book1D("GlobalResidualsClu3La3","RPC Residuals Layer 3 Cluster Size 3",101,-10.,10.);
00184 hGlobalResClu3La4 = dbe->book1D("GlobalResidualsClu3La4","RPC Residuals Layer 4 Cluster Size 3",101,-10.,10.);
00185 hGlobalResClu3La5 = dbe->book1D("GlobalResidualsClu3La5","RPC Residuals Layer 5 Cluster Size 3",101,-10.,10.);
00186 hGlobalResClu3La6 = dbe->book1D("GlobalResidualsClu3La6","RPC Residuals Layer 6 Cluster Size 3",101,-10.,10.);
00187
00188 if(debug) std::cout<<"Booking Residuals for EndCap"<<std::endl;
00189 folder = "Muons/MuonSegEff/Residuals/EndCap";
00190 dbe->setCurrentFolder(folder);
00191
00192
00193 hGlobalResClu1R3C = dbe->book1D("GlobalResidualsClu1R3C","RPC Residuals Ring 3 Roll C Cluster Size 1",101,-10.,10.);
00194 hGlobalResClu1R3B = dbe->book1D("GlobalResidualsClu1R3B","RPC Residuals Ring 3 Roll B Cluster Size 1",101,-10.,10.);
00195 hGlobalResClu1R3A = dbe->book1D("GlobalResidualsClu1R3A","RPC Residuals Ring 3 Roll A Cluster Size 1",101,-10.,10.);
00196 hGlobalResClu1R2C = dbe->book1D("GlobalResidualsClu1R2C","RPC Residuals Ring 2 Roll C Cluster Size 1",101,-10.,10.);
00197 hGlobalResClu1R2B = dbe->book1D("GlobalResidualsClu1R2B","RPC Residuals Ring 2 Roll B Cluster Size 1",101,-10.,10.);
00198 hGlobalResClu1R2A = dbe->book1D("GlobalResidualsClu1R2A","RPC Residuals Ring 2 Roll A Cluster Size 1",101,-10.,10.);
00199
00200 hGlobalResClu2R3C = dbe->book1D("GlobalResidualsClu2R3C","RPC Residuals Ring 3 Roll C Cluster Size 2",101,-10.,10.);
00201 hGlobalResClu2R3B = dbe->book1D("GlobalResidualsClu2R3B","RPC Residuals Ring 3 Roll B Cluster Size 2",101,-10.,10.);
00202 hGlobalResClu2R3A = dbe->book1D("GlobalResidualsClu2R3A","RPC Residuals Ring 3 Roll A Cluster Size 2",101,-10.,10.);
00203 hGlobalResClu2R2C = dbe->book1D("GlobalResidualsClu2R2C","RPC Residuals Ring 2 Roll C Cluster Size 2",101,-10.,10.);
00204 hGlobalResClu2R2B = dbe->book1D("GlobalResidualsClu2R2B","RPC Residuals Ring 2 Roll B Cluster Size 2",101,-10.,10.);
00205 hGlobalResClu2R2A = dbe->book1D("GlobalResidualsClu2R2A","RPC Residuals Ring 2 Roll A Cluster Size 2",101,-10.,10.);
00206
00207 hGlobalResClu3R3C = dbe->book1D("GlobalResidualsClu3R3C","RPC Residuals Ring 3 Roll C Cluster Size 3",101,-10.,10.);
00208 hGlobalResClu3R3B = dbe->book1D("GlobalResidualsClu3R3B","RPC Residuals Ring 3 Roll B Cluster Size 3",101,-10.,10.);
00209 hGlobalResClu3R3A = dbe->book1D("GlobalResidualsClu3R3A","RPC Residuals Ring 3 Roll A Cluster Size 3",101,-10.,10.);
00210 hGlobalResClu3R2C = dbe->book1D("GlobalResidualsClu3R2C","RPC Residuals Ring 2 Roll C Cluster Size 3",101,-10.,10.);
00211 hGlobalResClu3R2B = dbe->book1D("GlobalResidualsClu3R2B","RPC Residuals Ring 2 Roll B Cluster Size 3",101,-10.,10.);
00212 hGlobalResClu3R2A = dbe->book1D("GlobalResidualsClu3R2A","RPC Residuals Ring 2 Roll A Cluster Size 3",101,-10.,10.);
00213
00214
00215 if(debug) ofrej.open("rejected.txt");
00216
00217 if(debug) std::cout<<"Rejected done"<<std::endl;
00218
00219 }
00220
00221 void RPCEfficiency::beginRun(const edm::Run& run, const edm::EventSetup& iSetup){
00222
00223 iSetup.get<MuonGeometryRecord>().get(rpcGeo);
00224 iSetup.get<MuonGeometryRecord>().get(dtGeo);
00225 iSetup.get<MuonGeometryRecord>().get(cscGeo);
00226
00227 for (TrackingGeometry::DetContainer::const_iterator it=rpcGeo->dets().begin();it<rpcGeo->dets().end();it++){
00228 if(dynamic_cast< RPCChamber* >( *it ) != 0 ){
00229 RPCChamber* ch = dynamic_cast< RPCChamber* >( *it );
00230 std::vector< const RPCRoll*> roles = (ch->rolls());
00231 for(std::vector<const RPCRoll*>::const_iterator r = roles.begin();r != roles.end(); ++r){
00232 RPCDetId rpcId = (*r)->id();
00233 int region=rpcId.region();
00234
00235 RPCGeomServ rpcsrv(rpcId);
00236 std::string nameRoll = rpcsrv.name();
00237
00238 meCollection[nameRoll] = bookDetUnitSeg(rpcId,(*r)->nstrips());
00239
00240 if(region==0&&(incldt||incldtMB4)){
00241
00242 int wheel=rpcId.ring();
00243 int sector=rpcId.sector();
00244 int station=rpcId.station();
00245 DTStationIndex ind(region,wheel,sector,station);
00246 std::set<RPCDetId> myrolls;
00247 if (rollstoreDT.find(ind)!=rollstoreDT.end()) myrolls=rollstoreDT[ind];
00248 myrolls.insert(rpcId);
00249 rollstoreDT[ind]=myrolls;
00250
00251 }
00252 if(region!=0 && inclcsc){
00253
00254 int region=rpcId.region();
00255 int station=rpcId.station();
00256 int ring=rpcId.ring();
00257 int cscring=ring;
00258 int cscstation=station;
00259 RPCGeomServ rpcsrv(rpcId);
00260 int rpcsegment = mySegment(rpcId);
00261
00262 int cscchamber = rpcsegment;
00263 if((station==2||station==3)&&ring==3){
00264 cscring = 2;
00265 }
00266
00267 CSCStationIndex ind(region,cscstation,cscring,cscchamber);
00268 std::set<RPCDetId> myrolls;
00269 if (rollstoreCSC.find(ind)!=rollstoreCSC.end()){
00270 myrolls=rollstoreCSC[ind];
00271 }
00272 myrolls.insert(rpcId);
00273 rollstoreCSC[ind]=myrolls;
00274
00275 }
00276 }
00277 }
00278 }
00279
00280
00281
00282
00283
00284
00285
00286
00287
00288
00289
00290
00291
00292
00293
00294
00295
00296
00297
00298
00299
00300
00301
00302
00303
00304
00305
00306
00307
00308
00309
00310
00311
00312
00313
00314
00315
00316
00317
00318
00319
00320
00321
00322
00323
00324
00325
00326
00327
00328
00329
00330
00331
00332
00333
00334
00335
00336
00337
00338
00339
00340
00341
00342
00343
00344
00345
00346
00347
00348
00349
00350
00351
00352
00353
00354
00355
00356
00357
00358 }
00359
00360
00361 RPCEfficiency::~RPCEfficiency()
00362 {
00363
00364 }
00365
00366 void RPCEfficiency::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
00367 {
00368
00369 statistics->Fill(1);
00370 using namespace edm;
00371
00372 char layerLabel[128];
00373 char meIdRPC [128];
00374 char meIdDT [128];
00375 char meIdCSC [128];
00376
00377
00378
00379 if(debug) std::cout <<"Digi Getting the RPC Digis"<<std::endl;
00380 edm::Handle<RPCDigiCollection> rpcDigis;
00381 iEvent.getByLabel(muonRPCDigis, rpcDigis);
00382 char detUnitLabel[128];
00383 for (TrackingGeometry::DetContainer::const_iterator it=rpcGeo->dets().begin();it<rpcGeo->dets().end();it++){
00384 if(dynamic_cast< RPCChamber* >( *it ) != 0 ){
00385 RPCChamber* ch = dynamic_cast< RPCChamber* >( *it );
00386 std::vector< const RPCRoll*> roles = (ch->rolls());
00387 for(std::vector<const RPCRoll*>::const_iterator r = roles.begin();r != roles.end(); ++r){
00388 RPCDetId rpcId = (*r)->id();
00389 RPCGeomServ rpcsrv(rpcId);
00390 std::string nameRoll = rpcsrv.name();
00391 std::map<std::string, MonitorElement*> meMap=meCollection[nameRoll];
00392 sprintf(detUnitLabel ,"%s",nameRoll.c_str());
00393
00394 sprintf(layerLabel ,"%s",nameRoll.c_str());
00395 RPCDigiCollection::Range rpcRangeDigi=rpcDigis->get(rpcId);
00396 for (RPCDigiCollection::const_iterator digiIt = rpcRangeDigi.first;digiIt!=rpcRangeDigi.second;++digiIt){
00397 int stripDetected=digiIt->strip();
00398 sprintf(meIdRPC,"BXDistribution_%s",detUnitLabel);
00399 meMap[meIdRPC]->Fill(digiIt->bx());
00400 if(rpcId.region()==0){
00401 sprintf(meIdRPC,"RealDetectedOccupancyFromDT_%s",detUnitLabel);
00402 }else {
00403 sprintf(meIdRPC,"RealDetectedOccupancyFromCSC_%s",detUnitLabel);
00404 }
00405 meMap[meIdRPC]->Fill(stripDetected);
00406
00407 }
00408 }
00409 }
00410 }
00411
00412
00413
00414 if(debug) std::cout <<"\t Getting the RPC RecHits"<<std::endl;
00415 Handle<RPCRecHitCollection> rpcHits;
00416 iEvent.getByType(rpcHits);
00417
00418 if(incldt){
00419 if(debug) std::cout<<"\t Getting the DT Segments"<<std::endl;
00420 edm::Handle<DTRecSegment4DCollection> all4DSegments;
00421 iEvent.getByLabel(dt4DSegments, all4DSegments);
00422 if(all4DSegments->size()>0){
00423 if(all4DSegments->size()<=16) statistics->Fill(2);
00424
00425 if(debug) std::cout<<"\t Number of DT Segments in this event = "<<all4DSegments->size()<<std::endl;
00426
00427 std::map<DTChamberId,int> DTSegmentCounter;
00428 DTRecSegment4DCollection::const_iterator segment;
00429
00430 for (segment = all4DSegments->begin();segment!=all4DSegments->end(); ++segment){
00431 DTSegmentCounter[segment->chamberId()]++;
00432 }
00433
00434 statistics->Fill(all4DSegments->size()+2);
00435
00436 if(debug) std::cout<<"\t Loop over all the 4D Segments"<<std::endl;
00437 for (segment = all4DSegments->begin(); segment != all4DSegments->end(); ++segment){
00438
00439 DTChamberId DTId = segment->chamberId();
00440
00441
00442 if(debug) std::cout<<"DT \t \t This Segment is in Chamber id: "<<DTId<<std::endl;
00443 if(debug) std::cout<<"DT \t \t Number of segments in this DT = "<<DTSegmentCounter[DTId]<<std::endl;
00444 if(debug) std::cout<<"DT \t \t Is the only one in this DT? and is not in the 4th Station?"<<std::endl;
00445
00446
00447 if(DTSegmentCounter[DTId]==1 && DTId.station()!=4){
00448
00449 int dtWheel = DTId.wheel();
00450 int dtStation = DTId.station();
00451 int dtSector = DTId.sector();
00452
00453 LocalPoint segmentPosition= segment->localPosition();
00454 LocalVector segmentDirection=segment->localDirection();
00455
00456 const GeomDet* gdet=dtGeo->idToDet(segment->geographicalId());
00457 const BoundPlane & DTSurface = gdet->surface();
00458
00459
00460
00461 if(debug) std::cout<<"DT \t \t Is the segment 4D?"<<std::endl;
00462
00463 if(segment->dimension()==4){
00464
00465 if(debug) std::cout<<"DT \t \t yes"<<std::endl;
00466 if(debug) std::cout<<"DT \t \t DT Segment Dimension "<<segment->dimension()<<std::endl;
00467
00468 float Xo=segmentPosition.x();
00469 float Yo=segmentPosition.y();
00470 float Zo=segmentPosition.z();
00471 float dx=segmentDirection.x();
00472 float dy=segmentDirection.y();
00473 float dz=segmentDirection.z();
00474
00475 std::set<RPCDetId> rollsForThisDT = rollstoreDT[DTStationIndex(0,dtWheel,dtSector,dtStation)];
00476
00477 if(debug) std::cout<<"DT \t \t Number of rolls for this DT = "<<rollsForThisDT.size()<<std::endl;
00478
00479 assert(rollsForThisDT.size()>=1);
00480
00481 if(debug) std::cout<<"DT \t \t Loop over all the rolls asociated to this DT"<<std::endl;
00482 for (std::set<RPCDetId>::iterator iteraRoll = rollsForThisDT.begin();iteraRoll != rollsForThisDT.end(); iteraRoll++){
00483 const RPCRoll* rollasociated = rpcGeo->roll(*iteraRoll);
00484 RPCDetId rpcId = rollasociated->id();
00485 const BoundPlane & RPCSurface = rollasociated->surface();
00486
00487 RPCGeomServ rpcsrv(rpcId);
00488 std::string nameRoll = rpcsrv.name();
00489
00490 if(debug) std::cout<<"DT \t \t \t RollName: "<<nameRoll<<std::endl;
00491 if(debug) std::cout<<"DT \t \t \t Doing the extrapolation to this roll"<<std::endl;
00492 if(debug) std::cout<<"DT \t \t \t DT Segment Direction in DTLocal "<<segmentDirection<<std::endl;
00493 if(debug) std::cout<<"DT \t \t \t DT Segment Point in DTLocal "<<segmentPosition<<std::endl;
00494
00495 GlobalPoint CenterPointRollGlobal = RPCSurface.toGlobal(LocalPoint(0,0,0));
00496
00497 LocalPoint CenterRollinDTFrame = DTSurface.toLocal(CenterPointRollGlobal);
00498
00499 if(debug) std::cout<<"DT \t \t \t Center (0,0,0) Roll In DTLocal"<<CenterRollinDTFrame<<std::endl;
00500 if(debug) std::cout<<"DT \t \t \t Center (0,0,0) of the Roll in Global"<<CenterPointRollGlobal<<std::endl;
00501
00502 float D=CenterRollinDTFrame.z();
00503
00504 float X=Xo+dx*D/dz;
00505 float Y=Yo+dy*D/dz;
00506 float Z=D;
00507
00508 const RectangularStripTopology* top_= dynamic_cast<const RectangularStripTopology*> (&(rollasociated->topology()));
00509 LocalPoint xmin = top_->localPosition(0.);
00510 if(debug) std::cout<<"DT \t \t \t xmin of this Roll "<<xmin<<"cm"<<std::endl;
00511 LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
00512 if(debug) std::cout<<"DT \t \t \t xmax of this Roll "<<xmax<<"cm"<<std::endl;
00513 float rsize = fabs( xmax.x()-xmin.x() );
00514 if(debug) std::cout<<"DT \t \t \t Roll Size "<<rsize<<"cm"<<std::endl;
00515 float stripl = top_->stripLength();
00516 float stripw = top_->pitch();
00517
00518 if(debug) std::cout<<"DT \t \t \t Strip Lenght "<<stripl<<"cm"<<std::endl;
00519 if(debug) std::cout<<"DT \t \t \t Strip Width "<<stripw<<"cm"<<std::endl;
00520 if(debug) std::cout<<"DT \t \t \t X Predicted in DTLocal= "<<X<<"cm"<<std::endl;
00521 if(debug) std::cout<<"DT \t \t \t Y Predicted in DTLocal= "<<Y<<"cm"<<std::endl;
00522 if(debug) std::cout<<"DT \t \t \t Z Predicted in DTLocal= "<<Z<<"cm"<<std::endl;
00523
00524 float extrapolatedDistance = sqrt((X-Xo)*(X-Xo)+(Y-Yo)*(Y-Yo)+(Z-Zo)*(Z-Zo));
00525
00526 if(debug) std::cout<<"DT \t \t \t Is the distance of extrapolation less than MaxD? ="<<extrapolatedDistance<<"cm"<<"MaxD="<<MaxD<<"cm"<<std::endl;
00527
00528 if(extrapolatedDistance<=MaxD){
00529 if(debug) std::cout<<"DT \t \t \t yes"<<std::endl;
00530 GlobalPoint GlobalPointExtrapolated = DTSurface.toGlobal(LocalPoint(X,Y,Z));
00531 if(debug) std::cout<<"DT \t \t \t Point ExtraPolated in Global"<<GlobalPointExtrapolated<< std::endl;
00532 LocalPoint PointExtrapolatedRPCFrame = RPCSurface.toLocal(GlobalPointExtrapolated);
00533
00534 if(debug) std::cout<<"DT \t \t \t Point Extrapolated in RPCLocal"<<PointExtrapolatedRPCFrame<< std::endl;
00535 if(debug) std::cout<<"DT \t \t \t Corner of the Roll = ("<<rsize*0.5<<","<<stripl*0.5<<")"<<std::endl;
00536 if(debug) std::cout<<"DT \t \t \t Info About the Point Extrapolated in X Abs ("<<fabs(PointExtrapolatedRPCFrame.x())<<","
00537 <<fabs(PointExtrapolatedRPCFrame.y())<<","<<fabs(PointExtrapolatedRPCFrame.z())<<")"<<std::endl;
00538 if(debug) std::cout<<"DT \t \t \t Does the extrapolation go inside this roll?"<<std::endl;
00539
00540 if(fabs(PointExtrapolatedRPCFrame.z()) < 10. &&
00541 fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
00542 fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
00543
00544 if(debug) std::cout<<"DT \t \t \t \t yes"<<std::endl;
00545
00546 RPCDetId rollId = rollasociated->id();
00547
00548 RPCGeomServ rpcsrv(rollId);
00549 std::string nameRoll = rpcsrv.name();
00550 if(debug) std::cout<<"DT \t \t \t \t The RPCName is "<<nameRoll<<std::endl;
00551 const float stripPredicted =
00552 rollasociated->strip(LocalPoint(PointExtrapolatedRPCFrame.x(),PointExtrapolatedRPCFrame.y(),0.));
00553
00554 if(debug) std::cout<<"DT \t \t \t \t Candidate (from DT Segment) STRIP---> "<<stripPredicted<< std::endl;
00555
00556 char detUnitLabel[128];
00557 sprintf(detUnitLabel ,"%s",nameRoll.c_str());
00558 sprintf(layerLabel ,"%s",nameRoll.c_str());
00559
00560 std::map<std::string, MonitorElement*> meMap=meCollection[nameRoll];
00561
00562 sprintf(meIdDT,"ExpectedOccupancyFromDT_%s",detUnitLabel);
00563 if(debug) std::cout<<"DT \t \t \t \t Filling Expected for "<<meIdDT<<" with "<<stripPredicted<<std::endl;
00564 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;
00565 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;
00566 meMap[meIdDT]->Fill(stripPredicted);
00567
00568
00569
00570
00571 int countRecHits = 0;
00572 int cluSize = 0;
00573 float minres = 3000.;
00574
00575 if(debug) std::cout<<"DT \t \t \t \t Getting RecHits in Roll Asociated"<<std::endl;
00576 typedef std::pair<RPCRecHitCollection::const_iterator, RPCRecHitCollection::const_iterator> rangeRecHits;
00577 rangeRecHits recHitCollection = rpcHits->get(rollasociated->id());
00578 RPCRecHitCollection::const_iterator recHit;
00579
00580 for (recHit = recHitCollection.first; recHit != recHitCollection.second ; recHit++) {
00581 countRecHits++;
00582 LocalPoint recHitPos=recHit->localPosition();
00583 float res=PointExtrapolatedRPCFrame.x()- recHitPos.x();
00584 if(debug) std::cout<<"DT \t \t \t \t \t Found Rec Hit at "<<res<<"cm of the prediction."<<std::endl;
00585 if(fabs(res)<fabs(minres)){
00586 minres=res;
00587 cluSize = recHit->clusterSize();
00588 if(debug) std::cout<<"DT \t \t \t \t \t \t New Min Res "<<res<<"cm."<<std::endl;
00589 }
00590 }
00591
00592 bool anycoincidence=false;
00593
00594 if(countRecHits==0){
00595 if(debug) std::cout <<"DT \t \t \t \t \t THIS ROLL DOESN'T HAVE ANY RECHIT"<<std::endl;
00596 }else{
00597 assert(minres!=3000);
00598
00599 if(debug) std::cout<<"DT \t \t \t \t \t PointExtrapolatedRPCFrame.x="<<PointExtrapolatedRPCFrame.x()<<" Minimal Residual="<<minres<<std::endl;
00600 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;
00601 if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
00602 if(debug) std::cout<<"DT \t \t \t \t \t \t True!"<<std::endl;
00603 anycoincidence=true;
00604 }
00605 }
00606 if(anycoincidence){
00607 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;
00608 if(debug) std::cout<<"DT \t \t \t \t \t Norm of Cosine Directors="<<dx*dx+dy*dy+dz*dz<<"~1?"<<std::endl;
00609
00610
00611 if(inves){
00612 float cosal = dx/sqrt(dx*dx+dz*dz);
00613 if(debug) std::cout<<"DT \t \t \t \t \t Angle="<<acos(cosal)*180/3.1415926<<" degree"<<std::endl;
00614 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;
00615 if(rollId.station()==1&&rollId.layer()==1) { if(cluSize==1*dupli) {hGlobalResClu1La1->Fill(minres); AngClu1La1->Fill(acos(cosal)); DistBorderClu1La1->Fill(minres/stripw+0.5);}if(cluSize==2*dupli){ hGlobalResClu2La1->Fill(minres); AngClu2La1->Fill(acos(cosal));DistBorderClu2La1->Fill(minres/stripw+0.5);} if(cluSize==3*dupli){ hGlobalResClu3La1->Fill(minres); AngClu3La1->Fill(acos(cosal));} DistBorderClu3La1->Fill(minres/stripw+0.5);}
00616 else if(rollId.station()==1&&rollId.layer()==2){ if(cluSize==1*dupli) {hGlobalResClu1La2->Fill(minres); AngClu1La2->Fill(acos(cosal)); DistBorderClu1La2->Fill(minres/stripw+0.5);}if(cluSize==2*dupli){ hGlobalResClu2La2->Fill(minres); AngClu2La2->Fill(acos(cosal));DistBorderClu2La2->Fill(minres/stripw+0.5);} if(cluSize==3*dupli){ hGlobalResClu3La2->Fill(minres); AngClu3La2->Fill(acos(cosal));} DistBorderClu3La2->Fill(minres/stripw+0.5);}
00617 else if(rollId.station()==2&&rollId.layer()==1){ if(cluSize==1*dupli) {hGlobalResClu1La3->Fill(minres); AngClu1La3->Fill(acos(cosal)); DistBorderClu1La3->Fill(minres/stripw+0.5);}if(cluSize==2*dupli){ hGlobalResClu2La3->Fill(minres); AngClu2La3->Fill(acos(cosal));DistBorderClu2La3->Fill(minres/stripw+0.5);} if(cluSize==3*dupli){ hGlobalResClu3La3->Fill(minres); AngClu3La3->Fill(acos(cosal));} DistBorderClu3La3->Fill(minres/stripw+0.5);}
00618 else if(rollId.station()==2&&rollId.layer()==2){ if(cluSize==1*dupli) {hGlobalResClu1La4->Fill(minres); AngClu1La4->Fill(acos(cosal)); DistBorderClu1La4->Fill(minres/stripw+0.5);}if(cluSize==2*dupli){ hGlobalResClu2La4->Fill(minres); AngClu2La4->Fill(acos(cosal));DistBorderClu2La4->Fill(minres/stripw+0.5);} if(cluSize==3*dupli){ hGlobalResClu3La4->Fill(minres); AngClu3La4->Fill(acos(cosal));} DistBorderClu3La4->Fill(minres/stripw+0.5);}
00619 else if(rollId.station()==3) { if(cluSize==1*dupli) {hGlobalResClu1La5->Fill(minres); AngClu1La5->Fill(acos(cosal)); DistBorderClu1La5->Fill(minres/stripw+0.5);}if(cluSize==2*dupli){ hGlobalResClu2La5->Fill(minres); AngClu2La5->Fill(acos(cosal));DistBorderClu2La5->Fill(minres/stripw+0.5);} if(cluSize==3*dupli){ hGlobalResClu3La5->Fill(minres); AngClu3La5->Fill(acos(cosal));} DistBorderClu3La5->Fill(minres/stripw+0.5);}
00620
00621 }
00622 sprintf(meIdRPC,"RPCDataOccupancyFromDT_%s",detUnitLabel);
00623 if(debug) std::cout<<"DT \t \t \t \t \t COINCIDENCE!!! Event="<<iEvent.id()<<" Filling RPC Data Occupancy for "<<meIdRPC<<" with "<<stripPredicted<<std::endl;
00624 meMap[meIdRPC]->Fill(stripPredicted);
00625 }
00626 else{
00627 RPCGeomServ rpcsrv(rollasociated->id());
00628 std::string nameRoll = rpcsrv.name();
00629 if(debug) std::cout<<"DT \t \t \t \t \t A roll was ineficient in event "<<iEvent.id().event()<<std::endl;
00630 if(debug) ofrej<<"DTs \t Wh "<<dtWheel
00631 <<"\t St "<<dtStation
00632 <<"\t Se "<<dtSector
00633 <<"\t Roll "<<nameRoll
00634 <<"\t Event "
00635 <<iEvent.id().event()
00636 <<"\t Run "
00637 <<iEvent.id().run()
00638 <<std::endl;
00639 }
00640 }else {
00641 if(debug) std::cout<<"DT \t \t \t \t No the prediction is outside of this roll"<<std::endl;
00642 }
00643 }else{
00644 if(debug) std::cout<<"DT \t \t \t No, Exrtrapolation too long!, canceled"<<std::endl;
00645 }
00646 }
00647 }
00648 }else {
00649 if(debug) std::cout<<"DT \t \t More than one segment in this chamber, or we are in Station 4"<<std::endl;
00650 }
00651 }
00652 }
00653 else {
00654 if(debug) std::cout<<"DT This Event doesn't have any DT4DDSegment"<<std::endl;
00655 }
00656 }
00657
00658 if(incldtMB4){
00659 if(debug) std::cout <<"MB4 \t Getting ALL the DT Segments"<<std::endl;
00660 edm::Handle<DTRecSegment4DCollection> all4DSegments;
00661 iEvent.getByLabel(dt4DSegments, all4DSegments);
00662
00663 if(all4DSegments->size()>0){
00664
00665 std::map<DTChamberId,int> DTSegmentCounter;
00666 DTRecSegment4DCollection::const_iterator segment;
00667
00668 for (segment = all4DSegments->begin();segment!=all4DSegments->end(); ++segment){
00669 DTSegmentCounter[segment->chamberId()]++;
00670 }
00671
00672 if(debug) std::cout<<"MB4 \t \t Loop Over all4DSegments"<<std::endl;
00673 for (segment = all4DSegments->begin(); segment != all4DSegments->end(); ++segment){
00674
00675 DTChamberId DTId = segment->chamberId();
00676
00677 if(debug) std::cout<<"MB4 \t \t \t Is the only one in the chamber? and is in the Station 4?"<<std::endl;
00678
00679 if(DTSegmentCounter[DTId] == 1 && DTId.station()==4){
00680
00681 if(debug) std::cout<<"MB4 \t \t \t yes"<<std::endl;
00682 int dtWheel = DTId.wheel();
00683 int dtStation = DTId.station();
00684 int dtSector = DTId.sector();
00685
00686 LocalPoint segmentPosition= segment->localPosition();
00687 LocalVector segmentDirection=segment->localDirection();
00688
00689
00690
00691
00692 if(debug) std::cout<<"MB4 \t \t \t \t The Segment in MB4 is 2D?"<<std::endl;
00693 if(segment->dimension()==2){
00694 if(debug) std::cout<<"MB4 \t \t \t \t yes"<<std::endl;
00695 LocalVector segmentDirectionMB4=segmentDirection;
00696 LocalPoint segmentPositionMB4=segmentPosition;
00697
00698 bool compatiblesegments=false;
00699 float dx=segmentDirectionMB4.x();
00700 float dz=segmentDirectionMB4.z();
00701
00702 const BoundPlane& DTSurface4 = dtGeo->idToDet(DTId)->surface();
00703
00704 DTRecSegment4DCollection::const_iterator segMB3;
00705
00706 if(debug) std::cout<<"MB4 \t \t \t \t Loop on segments in =sector && MB3 && adjacent Wheel && y dim=4"<<std::endl;
00707 for(segMB3=all4DSegments->begin();segMB3!=all4DSegments->end();++segMB3){
00708
00709 DTChamberId dtid3 = segMB3->chamberId();
00710
00711 if(dtid3.sector()==DTId.sector()
00712 && dtid3.station()==3
00713 && abs(dtid3.wheel()-DTId.wheel())<2
00714 && DTSegmentCounter[dtid3] == 1
00715 && segMB3->dimension()==4){
00716
00717 const GeomDet* gdet3=dtGeo->idToDet(segMB3->geographicalId());
00718 const BoundPlane & DTSurface3 = gdet3->surface();
00719
00720 float dx3=segMB3->localDirection().x();
00721 float dy3=segMB3->localDirection().y();
00722 float dz3=segMB3->localDirection().z();
00723
00724 LocalVector segDirMB4inMB3Frame=DTSurface3.toLocal(DTSurface4.toGlobal(segmentDirectionMB4));
00725
00726 double cosAng=fabs(dx*dx3+dz*dz3/sqrt((dx3*dx3+dz3*dz3)*(dx*dx+dz*dz)));
00727 if(fabs(cosAng)>1.){
00728 std::cout<<"dx="<<dx<<" dz="<<dz<<std::endl;
00729 std::cout<<"dx3="<<dx3<<" dz3="<<dz<<std::endl;
00730 std::cout<<cosAng<<std::endl;
00731 }
00732
00733 if(cosAng>MinCosAng){
00734 compatiblesegments=true;
00735 if(dtSector==13){
00736 dtSector=4;
00737 }
00738 if(dtSector==14){
00739 dtSector=10;
00740 }
00741
00742 std::set<RPCDetId> rollsForThisDT = rollstoreDT[DTStationIndex(0,dtWheel,dtSector,4)];
00743
00744 assert(rollsForThisDT.size()>=1);
00745
00746 for (std::set<RPCDetId>::iterator iteraRoll=rollsForThisDT.begin();iteraRoll != rollsForThisDT.end(); iteraRoll++){
00747 const RPCRoll* rollasociated = rpcGeo->roll(*iteraRoll);
00748 const BoundPlane & RPCSurfaceRB4 = rollasociated->surface();
00749 const GeomDet* gdet=dtGeo->idToDet(segMB3->geographicalId());
00750 const BoundPlane & DTSurfaceMB3 = gdet->surface();
00751
00752 GlobalPoint CenterPointRollGlobal=RPCSurfaceRB4.toGlobal(LocalPoint(0,0,0));
00753
00754 LocalPoint CenterRollinMB3Frame = DTSurfaceMB3.toLocal(CenterPointRollGlobal);
00755
00756 float D=CenterRollinMB3Frame.z();
00757
00758 float Xo3=segMB3->localPosition().x();
00759 float Yo3=segMB3->localPosition().y();
00760 float Zo3=segMB3->localPosition().z();
00761
00762 float X=Xo3+dx3*D/dz3;
00763 float Y=Yo3+dy3*D/dz3;
00764 float Z=D;
00765
00766
00767 const RectangularStripTopology* top_
00768 =dynamic_cast<const RectangularStripTopology*>(&(rollasociated->topology()));
00769 LocalPoint xmin = top_->localPosition(0.);
00770 LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
00771 float rsize = fabs( xmax.x()-xmin.x() );
00772 float stripl = top_->stripLength();
00773 float stripw = top_->pitch();
00774
00775
00776 if(debug) std::cout<<"MB4 \t \t \t Strip Lenght "<<stripl<<"cm"<<std::endl;
00777 if(debug) std::cout<<"MB4 \t \t \t Strip Width "<<stripw<<"cm"<<std::endl;
00778
00779 if(debug) std::cout<<"MB4 \t \t \t X Predicted in MB3Local= "<<X<<"cm"<<std::endl;
00780 if(debug) std::cout<<"MB4 \t \t \t Y Predicted in MB3DTLocal= "<<Y<<"cm"<<std::endl;
00781 if(debug) std::cout<<"MB4 \t \t \t Z Predicted in MB3DTLocal= "<<Z<<"cm"<<std::endl;
00782
00783 float extrapolatedDistance = sqrt((X-Xo3)*(X-Xo3)+(Y-Yo3)*(Y-Yo3)+(Z-Zo3)*(Z-Zo3));
00784
00785 if(extrapolatedDistance<=MaxDrb4){
00786 if(debug) std::cout<<"MB4 \t \t \t yes"<<std::endl;
00787
00788 GlobalPoint GlobalPointExtrapolated = DTSurfaceMB3.toGlobal(LocalPoint(X,Y,Z));
00789
00790 if(debug) std::cout<<"MB4 \t \t \t Point ExtraPolated in Global"<<GlobalPointExtrapolated<< std::endl;
00791
00792 LocalPoint PointExtrapolatedRPCFrame = RPCSurfaceRB4.toLocal(GlobalPointExtrapolated);
00793
00794 if(debug) std::cout<<"MB4 \t \t \t Point Extrapolated in RPCLocal"<<PointExtrapolatedRPCFrame<< std::endl;
00795 if(debug) std::cout<<"MB4 \t \t \t Corner of the Roll = ("<<rsize*0.5<<","<<stripl*0.5<<")"<<std::endl;
00796 if(debug) std::cout<<"MB4 \t \t \t Info About the Point Extrapolated in X Abs ("<<fabs(PointExtrapolatedRPCFrame.x())<<","
00797 <<fabs(PointExtrapolatedRPCFrame.y())<<","<<fabs(PointExtrapolatedRPCFrame.z())<<")"<<std::endl;
00798
00799 if(debug) std::cout<<"MB4 \t \t \t Does the extrapolation go inside this roll?"<<std::endl;
00800
00801 if(fabs(PointExtrapolatedRPCFrame.z()) < 10. &&
00802 fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
00803 fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
00804
00805 if(debug) std::cout<<"MB4 \t \t \t \t yes"<<std::endl;
00806
00807 RPCDetId rollId = rollasociated->id();
00808
00809 RPCGeomServ rpcsrv(rollId);
00810 std::string nameRoll = rpcsrv.name();
00811 if(debug) std::cout<<"MB4 \t \t \t \t \t The RPCName is "<<nameRoll<<std::endl;
00812 const float stripPredicted=
00813 rollasociated->strip(LocalPoint(PointExtrapolatedRPCFrame.x(),PointExtrapolatedRPCFrame.y(),0.));
00814
00815 if(debug) std::cout<<"MB4 \t \t \t \t Candidate (from DT Segment) STRIP---> "<<stripPredicted<< std::endl;
00816
00817 char detUnitLabel[128];
00818 sprintf(detUnitLabel ,"%s",nameRoll.c_str());
00819 sprintf(layerLabel ,"%s",nameRoll.c_str());
00820
00821 std::map<std::string, MonitorElement*> meMap=meCollection[nameRoll];
00822
00823 if(debug) std::cout<<"MB4 \t \t \t \t \t Filling Expected"<<std::endl;
00824 sprintf(meIdDT,"ExpectedOccupancyFromDT_%s",detUnitLabel);
00825 meMap[meIdDT]->Fill(stripPredicted);
00826
00827
00828
00829
00830 int countRecHits = 0;
00831 int cluSize = 0;
00832 float minres = 3000.;
00833
00834 if(debug) std::cout<<"MB4 \t \t \t \t Getting RecHits in Roll Asociated"<<std::endl;
00835 typedef std::pair<RPCRecHitCollection::const_iterator, RPCRecHitCollection::const_iterator> rangeRecHits;
00836 rangeRecHits recHitCollection = rpcHits->get(rollasociated->id());
00837 RPCRecHitCollection::const_iterator recHit;
00838
00839 for (recHit = recHitCollection.first; recHit != recHitCollection.second ; recHit++) {
00840 countRecHits++;
00841 LocalPoint recHitPos=recHit->localPosition();
00842 float res=PointExtrapolatedRPCFrame.x()- recHitPos.x();
00843 if(debug) std::cout<<"DT \t \t \t \t \t Found Rec Hit at "<<res<<"cm of the prediction."<<std::endl;
00844 if(fabs(res)<fabs(minres)){
00845 minres=res;
00846 cluSize = recHit->clusterSize();
00847 }
00848 }
00849
00850 bool anycoincidence=false;
00851
00852 if(countRecHits==0){
00853 if(debug) std::cout <<"MB4 \t \t \t \t \t \t THIS ROLL DOESN'T HAVE ANY RECHIT"<<std::endl;
00854 }else{
00855 assert(minres!=3000);
00856
00857 if(debug) std::cout<<"MB4 \t \t \t \t \t \t PointExtrapolatedRPCFrame.x="<<PointExtrapolatedRPCFrame.x()<<" Minimal Residual ="<<minres<<std::endl;
00858 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;
00859 if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
00860 if(debug) std::cout<<"MB4 \t \t \t \t \t \t \t True!"<<std::endl;
00861 anycoincidence=true;
00862 }
00863 }
00864 if(anycoincidence){
00865 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;
00866 if(debug) std::cout<<"MB4 \t \t \t \t \t Norm of Cosine Directors="<<dx3*dx3+dy3*dy3+dz3*dz3<<"~1?"<<std::endl;
00867
00868
00869 if(inves){
00870 float cosal = dx/sqrt(dx*dx+dz*dz);
00871 if(debug) std::cout<<"MB4 \t \t \t \t \t Angle="<<acos(cosal)*180/3.1415926<<" degree"<<std::endl;
00872 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;
00873 assert(rollId.station()==4);
00874 if(cluSize==1*dupli){ hGlobalResClu1La6->Fill(minres); if(minres<0.005*stripw) AngClu1La6->Fill(acos(cosal)); if(fabs(minres)<stripw*0.5) DistBorderClu1La6->Fill(minres/stripw+0.5);}
00875 else if(cluSize==2*dupli){ hGlobalResClu2La6->Fill(minres); if(minres<0.005*stripw) AngClu2La6->Fill(acos(cosal));}
00876 else if(cluSize==3*dupli){ hGlobalResClu3La6->Fill(minres); if(minres<0.005*stripw) AngClu3La6->Fill(acos(cosal));}
00877 }
00878
00879
00880 sprintf(meIdRPC,"RPCDataOccupancyFromDT_%s",detUnitLabel);
00881 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;
00882 meMap[meIdRPC]->Fill(stripPredicted);
00883 }
00884 else{
00885 RPCGeomServ rpcsrv(rollasociated->id());
00886 std::string nameRoll = rpcsrv.name();
00887 if(debug) std::cout<<"MB4 \t \t \t \t \t \t A roll was ineficient in event"<<iEvent.id().event()<<std::endl;
00888 if(debug) ofrej<<"MB4 \t Wh "<<dtWheel
00889 <<"\t St "<<dtStation
00890 <<"\t Se "<<dtSector
00891 <<"\t Roll "<<nameRoll
00892 <<"\t Event "
00893 <<iEvent.id().event()
00894 <<"\t Run "
00895 <<iEvent.id().run()
00896 <<std::endl;
00897 }
00898 }else{
00899 if(debug) std::cout<<"MB4 \t \t \t \t No the prediction is outside of this roll"<<std::endl;
00900 }
00901 }
00902 else{
00903 if(debug) std::cout<<"MB4 \t \t \t No, Exrtrapolation too long!, canceled"<<std::endl;
00904 }
00905 }
00906 }else{
00907 compatiblesegments=false;
00908 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;
00909 }
00910 }else{
00911 if(debug) std::cout<<"MB4 \t \t \t No the same station or same wheel or segment dim in mb3 not 4D"<<std::endl;
00912 }
00913 }
00914 }else{
00915 if(debug) std::cout<<"MB4 \t \t \t Is NOT a 2D Segment"<<std::endl;
00916 }
00917 }else{
00918 if(debug) std::cout<<"MB4 \t \t \t \t There is not just one segment or is not in station 4"<<std::endl;
00919 }
00920 }
00921 }else{
00922 if(debug) std::cout<<"MB4 \t This event doesn't have 4D Segment"<<std::endl;
00923 }
00924 }
00925
00926 if(inclcsc){
00927 if(debug) std::cout <<"\t Getting the CSC Segments"<<std::endl;
00928 edm::Handle<CSCSegmentCollection> allCSCSegments;
00929 iEvent.getByLabel(cscSegments, allCSCSegments);
00930
00931 if(allCSCSegments->size()>0){
00932 statistics->Fill(18);
00933
00934 if(debug) std::cout<<"CSC \t Number of CSC Segments in this event = "<<allCSCSegments->size()<<std::endl;
00935
00936 std::map<CSCDetId,int> CSCSegmentsCounter;
00937 CSCSegmentCollection::const_iterator segment;
00938
00939 int segmentsInThisEventInTheEndcap=0;
00940
00941 for (segment = allCSCSegments->begin();segment!=allCSCSegments->end(); ++segment){
00942 CSCSegmentsCounter[segment->cscDetId()]++;
00943 segmentsInThisEventInTheEndcap++;
00944 }
00945
00946 statistics->Fill(allCSCSegments->size()+18);
00947
00948
00949 if(debug) std::cout<<"CSC \t loop over all the CSCSegments "<<std::endl;
00950 for (segment = allCSCSegments->begin();segment!=allCSCSegments->end(); ++segment){
00951 CSCDetId CSCId = segment->cscDetId();
00952
00953 if(debug) std::cout<<"CSC \t \t This Segment is in Chamber id: "<<CSCId<<std::endl;
00954 if(debug) std::cout<<"CSC \t \t Number of segments in this CSC = "<<CSCSegmentsCounter[CSCId]<<std::endl;
00955 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;
00956
00957 if(CSCSegmentsCounter[CSCId]==1 && CSCId.station()!=4 && CSCId.ring()!=1 && allCSCSegments->size()>=2){
00958 if(debug) std::cout<<"CSC \t \t yes"<<std::endl;
00959 int cscEndCap = CSCId.endcap();
00960 int cscStation = CSCId.station();
00961 int cscRing = CSCId.ring();
00962 int cscChamber = CSCId.chamber();
00963 int rpcRegion = 1; if(cscEndCap==2) rpcRegion= -1;
00964 int rpcRing = cscRing;
00965 if(cscRing==4)rpcRing =1;
00966 int rpcStation = cscStation;
00967 int rpcSegment = 0;
00968
00969 if(cscStation!=1&&cscRing==1){
00970 rpcSegment = CSCId.chamber();
00971 }
00972 else{
00973 rpcSegment = (CSCId.chamber()==1) ? 36 : CSCId.chamber()-1;
00974 }
00975
00976 LocalPoint segmentPosition= segment->localPosition();
00977 LocalVector segmentDirection=segment->localDirection();
00978 float dz=segmentDirection.z();
00979
00980 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;
00981
00982 if(segment->dimension()==4 && (segment->nRecHits()<=10 && segment->nRecHits()>=4)&& acos(dz)*180/3.1415926 > 45. && acos(dz)*180/3.1415926 < 160. ){
00983
00984
00985
00986 if(debug) std::cout<<"CSC \t \t yes"<<std::endl;
00987 if(debug) std::cout<<"CSC \t \t CSC Segment Dimension "<<segment->dimension()<<std::endl;
00988
00989 float Xo=segmentPosition.x();
00990 float Yo=segmentPosition.y();
00991 float Zo=segmentPosition.z();
00992 float dx=segmentDirection.x();
00993 float dy=segmentDirection.y();
00994 float dz=segmentDirection.z();
00995
00996
00997 if(debug) std::cout<<"CSC \t \t Getting chamber from Geometry"<<std::endl;
00998 const CSCChamber* TheChamber=cscGeo->chamber(CSCId);
00999 if(debug) std::cout<<"CSC \t \t Getting ID from Chamber"<<std::endl;
01000 const CSCDetId TheId=TheChamber->id();
01001 if(debug) std::cout<<"CSC \t \t Printing The Id"<<TheId<<std::endl;
01002 std::set<RPCDetId> rollsForThisCSC = rollstoreCSC[CSCStationIndex(rpcRegion,rpcStation,rpcRing,rpcSegment)];
01003 if(debug) std::cout<<"CSC \t \t Number of rolls for this CSC = "<<rollsForThisCSC.size()<<std::endl;
01004
01005 if(debug) std::cout<<"CSC \t \t Loop over all the rolls asociated to this CSC"<<std::endl;
01006
01007 if(rpcRing!=1&&rpcStation!=4){
01008
01009 if(rollsForThisCSC.size()==0){
01010 if(debug) std::cout<<"CSC Fail for CSCId="<<TheId<<" rpcRegion="<<rpcRegion<<" rpcStation="<<rpcStation<<" rpcRing="<<rpcRing<<" rpcSegment="<<rpcSegment<<std::endl;
01011 }
01012
01013 assert(rollsForThisCSC.size()>=1);
01014
01015
01016 for (std::set<RPCDetId>::iterator iteraRoll = rollsForThisCSC.begin();iteraRoll != rollsForThisCSC.end(); iteraRoll++){
01017 const RPCRoll* rollasociated = rpcGeo->roll(*iteraRoll);
01018 RPCDetId rpcId = rollasociated->id();
01019
01020 if(debug) std::cout<<"CSC \t \t \t We are in the roll getting the surface"<<rpcId<<std::endl;
01021 const BoundPlane & RPCSurface = rollasociated->surface();
01022
01023 if(debug) std::cout<<"CSC \t \t \t RollID: "<<rpcId<<std::endl;
01024
01025 if(debug) std::cout<<"CSC \t \t \t Doing the extrapolation to this roll"<<std::endl;
01026 if(debug) std::cout<<"CSC \t \t \t CSC Segment Direction in CSCLocal "<<segmentDirection<<std::endl;
01027 if(debug) std::cout<<"CSC \t \t \t CSC Segment Point in CSCLocal "<<segmentPosition<<std::endl;
01028
01029 GlobalPoint CenterPointRollGlobal = RPCSurface.toGlobal(LocalPoint(0,0,0));
01030 if(debug) std::cout<<"CSC \t \t \t Center (0,0,0) of the Roll in Global"<<CenterPointRollGlobal<<std::endl;
01031 GlobalPoint CenterPointCSCGlobal = TheChamber->toGlobal(LocalPoint(0,0,0));
01032 if(debug) std::cout<<"CSC \t \t \t Center (0,0,0) of the CSC in Global"<<CenterPointCSCGlobal<<std::endl;
01033 GlobalPoint segmentPositionInGlobal=TheChamber->toGlobal(segmentPosition);
01034 if(debug) std::cout<<"CSC \t \t \t Segment Position in Global"<<segmentPositionInGlobal<<std::endl;
01035 LocalPoint CenterRollinCSCFrame = TheChamber->toLocal(CenterPointRollGlobal);
01036
01037 if(debug){
01038 float rpcphi=0;
01039 float cscphi=0;
01040
01041 (CenterPointRollGlobal.barePhi()<0)?
01042 rpcphi = 2*3.141592+CenterPointRollGlobal.barePhi():rpcphi=CenterPointRollGlobal.barePhi();
01043
01044 (CenterPointCSCGlobal.barePhi()<0)?
01045 cscphi = 2*3.1415926536+CenterPointCSCGlobal.barePhi():cscphi=CenterPointCSCGlobal.barePhi();
01046
01047 float df=fabs(cscphi-rpcphi);
01048 float dr=fabs(CenterPointRollGlobal.perp()-CenterPointCSCGlobal.perp());
01049 float diffz=CenterPointRollGlobal.z()-CenterPointCSCGlobal.z();
01050 float dfg=df*180./3.14159265;
01051
01052 if(debug) std::cout<<"CSC \t \t \t z of RPC="<<CenterPointRollGlobal.z()<<"z of CSC"<<CenterPointCSCGlobal.z()<<" dfg="<<dfg<<std::endl;
01053
01054
01055 RPCGeomServ rpcsrv(rpcId);
01056
01057
01058 if(dr>200.||fabs(dz)>55.||dfg>1.){
01059
01060 if (debug) std::cout
01061 <<"\t \t \t CSC Station= "<<CSCId.station()
01062 <<" Ring= "<<CSCId.ring()
01063 <<" Chamber= "<<CSCId.chamber()
01064 <<" cscphi="<<cscphi*180/3.14159265
01065 <<"\t RPC Station= "<<rpcId.station()
01066 <<" ring= "<<rpcId.ring()
01067 <<" segment =-> "<<mySegment(rpcId)
01068 <<" rollphi="<<rpcphi*180/3.14159265
01069 <<"\t dfg="<<dfg
01070 <<" dz="<<diffz
01071 <<" dr="<<dr
01072 <<std::endl;
01073
01074 }
01075 }
01076
01077
01078
01079
01080 float D=CenterRollinCSCFrame.z();
01081
01082 float X=Xo+dx*D/dz;
01083 float Y=Yo+dy*D/dz;
01084 float Z=D;
01085
01086 const TrapezoidalStripTopology* top_=dynamic_cast<const TrapezoidalStripTopology*>(&(rollasociated->topology()));
01087 LocalPoint xmin = top_->localPosition(0.);
01088 if(debug) std::cout<<"CSC \t \t \t xmin of this Roll "<<xmin<<"cm"<<std::endl;
01089 LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
01090 if(debug) std::cout<<"CSC \t \t \t xmax of this Roll "<<xmax<<"cm"<<std::endl;
01091 float rsize = fabs( xmax.x()-xmin.x() );
01092 if(debug) std::cout<<"CSC \t \t \t Roll Size "<<rsize<<"cm"<<std::endl;
01093 float stripl = top_->stripLength();
01094 float stripw = top_->pitch();
01095
01096 if(debug) std::cout<<"CSC \t \t \t Strip Lenght "<<stripl<<"cm"<<std::endl;
01097 if(debug) std::cout<<"CSC \t \t \t Strip Width "<<stripw<<"cm"<<std::endl;
01098
01099 if(debug) std::cout<<"CSC \t \t \t X Predicted in CSCLocal= "<<X<<"cm"<<std::endl;
01100 if(debug) std::cout<<"CSC \t \t \t Y Predicted in CSCLocal= "<<Y<<"cm"<<std::endl;
01101 if(debug) std::cout<<"CSC \t \t \t Z Predicted in CSCLocal= "<<Z<<"cm"<<std::endl;
01102
01103 float extrapolatedDistance = sqrt((X-Xo)*(X-Xo)+(Y-Yo)*(Y-Yo)+(Z-Zo)*(Z-Zo));
01104
01105 if(debug) std::cout<<"CSC \t \t \t Is the distance of extrapolation less than MaxD? ="<<extrapolatedDistance<<"cm"<<"MaxD="<<MaxD<<"cm"<<std::endl;
01106
01107 if(extrapolatedDistance<=MaxD){
01108
01109 if(debug) std::cout<<"CSC \t \t \t yes"<<std::endl;
01110
01111 GlobalPoint GlobalPointExtrapolated=TheChamber->toGlobal(LocalPoint(X,Y,Z));
01112 if(debug) std::cout<<"CSC \t \t \t Point ExtraPolated in Global"<<GlobalPointExtrapolated<< std::endl;
01113
01114
01115 LocalPoint PointExtrapolatedRPCFrame = RPCSurface.toLocal(GlobalPointExtrapolated);
01116
01117 if(debug) std::cout<<"CSC \t \t \t Point Extrapolated in RPCLocal"<<PointExtrapolatedRPCFrame<< std::endl;
01118 if(debug) std::cout<<"CSC \t \t \t Corner of the Roll = ("<<rsize*0.5<<","<<stripl*0.5<<")"<<std::endl;
01119 if(debug) std::cout<<"CSC \t \t \t Info About the Point Extrapolated in X Abs ("<<fabs(PointExtrapolatedRPCFrame.x())<<","
01120 <<fabs(PointExtrapolatedRPCFrame.y())<<","<<fabs(PointExtrapolatedRPCFrame.z())<<")"<<std::endl;
01121 if(debug) std::cout<<"CSC \t \t \t dz="
01122 <<fabs(PointExtrapolatedRPCFrame.z())<<" dx="
01123 <<fabs(PointExtrapolatedRPCFrame.x())<<" dy="
01124 <<fabs(PointExtrapolatedRPCFrame.y())<<std::endl;
01125
01126 if(debug) std::cout<<"CSC \t \t \t Does the extrapolation go inside this roll????"<<std::endl;
01127
01128 if(fabs(PointExtrapolatedRPCFrame.z()) < 10. &&
01129 fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
01130 fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
01131
01132 if(debug) std::cout<<"CSC \t \t \t \t yes"<<std::endl;
01133
01134 RPCDetId rollId = rollasociated->id();
01135
01136 RPCGeomServ rpcsrv(rollId);
01137 std::string nameRoll = rpcsrv.name();
01138 if(debug) std::cout<<"CSC \t \t \t \t The RPCName is "<<nameRoll<<std::endl;
01139
01140 const float stripPredicted =
01141 rollasociated->strip(LocalPoint(PointExtrapolatedRPCFrame.x(),PointExtrapolatedRPCFrame.y(),0.));
01142
01143 if(debug) std::cout<<"CSC \t \t \t \t \t Candidate"<<rollId<<" "<<"(from CSC Segment) STRIP---> "<<stripPredicted<< std::endl;
01144
01145
01146 char detUnitLabel[128];
01147 sprintf(detUnitLabel ,"%s",nameRoll.c_str());
01148 sprintf(layerLabel ,"%s",nameRoll.c_str());
01149
01150 std::map<std::string, MonitorElement*> meMap=meCollection[nameRoll];
01151
01152 if(debug) std::cout<<"CSC \t \t \t \t Filling Expected"<<std::endl;
01153 sprintf(meIdCSC,"ExpectedOccupancyFromCSC_%s",detUnitLabel);
01154 meMap[meIdCSC]->Fill(stripPredicted);
01155
01156
01157
01158
01159 int cluSize = 0;
01160 int countRecHits = 0;
01161 float minres = 3000.;
01162
01163 if(debug) std::cout<<"CSC \t \t \t \t \t Getting RecHits in Roll Asociated"<<std::endl;
01164 typedef std::pair<RPCRecHitCollection::const_iterator, RPCRecHitCollection::const_iterator> rangeRecHits;
01165 rangeRecHits recHitCollection = rpcHits->get(rollasociated->id());
01166 RPCRecHitCollection::const_iterator recHit;
01167
01168 for (recHit = recHitCollection.first; recHit != recHitCollection.second ; recHit++) {
01169 countRecHits++;
01170 LocalPoint recHitPos=recHit->localPosition();
01171 float res=PointExtrapolatedRPCFrame.x()- recHitPos.x();
01172 if(debug) std::cout<<"CSC \t \t \t \t \t \t Found Rec Hit at "<<res<<"cm of the prediction."<<std::endl;
01173 if(fabs(res)<fabs(minres)){
01174 minres=res;
01175 cluSize = recHit->clusterSize();
01176 if(debug) std::cout<<"CSC \t \t \t \t \t \t \t New Min Res "<<res<<"cm."<<std::endl;
01177 }
01178 }
01179
01180 bool anycoincidence = false;
01181
01182 if(countRecHits==0){
01183 if(debug) std::cout <<"CSC \t \t \t \t \t THIS ROLL DOESN'T HAVE ANY RECHIT"<<std::endl;
01184 }else{
01185 assert(minres!=3000);
01186
01187 if(debug) std::cout<<"CSC \t \t \t \t \t PointExtrapolatedRPCFrame.x="<<PointExtrapolatedRPCFrame.x()<<" Minimal Residual"<<minres<<std::endl;
01188 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;
01189 if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
01190 if(debug) std::cout<<"CSC \t \t \t \t \t \t True!"<<std::endl;
01191 anycoincidence=true;
01192 }
01193 }
01194 if(anycoincidence){
01195 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;
01196 if(debug) std::cout<<"CSC \t \t \t \t \t Norm of Cosine Directors="<<dx*dx+dy*dy+dz*dz<<"~1?"<<std::endl;
01197
01198
01199 if(inves){
01200 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);}
01201 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);}
01202 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);}
01203 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);}
01204 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);}
01205 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);}
01206 }
01207
01208 sprintf(meIdRPC,"RPCDataOccupancyFromCSC_%s",detUnitLabel);
01209 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;
01210 meMap[meIdRPC]->Fill(stripPredicted);
01211 }
01212 else{
01213 RPCGeomServ rpcsrv(rollasociated->id());
01214 std::string nameRoll = rpcsrv.name();
01215 if(debug) std::cout<<"CSC \t \t \t \t \t \t A roll was ineficient in event"<<iEvent.id().event()<<std::endl;
01216 if(debug) ofrej<<"CSC \t EndCap "<<rpcRegion
01217 <<"\t cscStation "<<cscStation
01218 <<"\t cscRing "<<cscRing
01219 <<"\t cscChamber "<<cscChamber
01220 <<"\t Roll "<<nameRoll
01221 <<"\t Event "<<iEvent.id().event()
01222 <<"\t CSCId "<<CSCId
01223 <<"\t Event "
01224 <<iEvent.id().event()
01225 <<"\t Run "
01226 <<iEvent.id().run()
01227 <<std::endl;
01228 }
01229 }else{
01230 if(debug) std::cout<<"CSC \t \t \t \t No the prediction is outside of this roll"<<std::endl;
01231 }
01232 }else{
01233 if(debug) std::cout<<"CSC \t \t \t No, Exrtrapolation too long!, canceled"<<std::endl;
01234 }
01235 }
01236 }
01237 }
01238 }else{
01239 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;
01240 }
01241 }
01242 }else{
01243 if(debug) std::cout<<"CSC This Event doesn't have any CSCSegment"<<std::endl;
01244 }
01245 }
01246 }
01247
01248 void RPCEfficiency::endRun(const edm::Run& r, const edm::EventSetup& iSetup){
01249 if (EffSaveRootFile){
01250 dbe->save(EffRootFileName);
01251 }
01252 }
01253
01254
01255 void RPCEfficiency::endJob()
01256 {
01257 dbe =0;
01258 }
1.5.4