RPCEfficiency Class Reference

#include <RPCEfficiency.h>

Inheritance diagram for RPCEfficiency:

Public Member Functions
	RPCEfficiency (const edm::ParameterSet &)
	~RPCEfficiency () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final
virtual void	analyze (edm::Event const &, edm::EventSetup const &)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) override
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override
	~DQMEDAnalyzer () override=default
Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookDetUnitSeg (DQMStore::IBooker &, RPCDetId &detId, int nstrips, std::string folder, std::map< std::string, MonitorElement * > &)
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Protected Attributes
MonitorElement *	hGlobalResClu1La [6]
MonitorElement *	hGlobalResClu1R2A
MonitorElement *	hGlobalResClu1R2B
MonitorElement *	hGlobalResClu1R2C
MonitorElement *	hGlobalResClu1R3A
MonitorElement *	hGlobalResClu1R3B
MonitorElement *	hGlobalResClu1R3C
MonitorElement *	hGlobalResClu2La [6]
MonitorElement *	hGlobalResClu2R2A
MonitorElement *	hGlobalResClu2R2B
MonitorElement *	hGlobalResClu2R2C
MonitorElement *	hGlobalResClu2R3A
MonitorElement *	hGlobalResClu2R3B
MonitorElement *	hGlobalResClu2R3C
MonitorElement *	hGlobalResClu3La [6]
MonitorElement *	hGlobalResClu3R2A
MonitorElement *	hGlobalResClu3R2B
MonitorElement *	hGlobalResClu3R2C
MonitorElement *	hGlobalResClu3R3A
MonitorElement *	hGlobalResClu3R3B
MonitorElement *	hGlobalResClu3R3C
std::map< CSCStationIndex, std::set< RPCDetId > >	rollstoreCSC
std::map< DTStationIndex, std::set< RPCDetId > >	rollstoreDT
MonitorElement *	statistics
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

Private Attributes
std::vector< std::map< RPCDetId, int > >	counter
edm::EDGetTokenT< CSCSegmentCollection >	cscSegments
bool	debug
edm::EDGetTokenT< DTRecSegment4DCollection >	dt4DSegments
int	dupli
std::string	EffRootFileName
bool	EffSaveRootFile
std::string	folderPath
bool	inclcsc
bool	incldt
bool	incldtMB4
bool	inves
double	MaxD
double	MaxDrb4
std::map< int, std::map< std::string, MonitorElement * > >	meCollection
double	MinCosAng
double	rangestrips
double	rangestripsRB4
std::string	rollseff
edm::EDGetTokenT< RPCRecHitCollection >	RPCRecHitLabel_
std::vector< int >	totalcounter

Additional Inherited Members
Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Class for RPC Monitoring using RPCDigi and DT and CSC Segments. author Camilo Carrillo (Uniandes)

Definition at line 95 of file RPCEfficiency.h.

Constructor & Destructor Documentation

RPCEfficiency::RPCEfficiency ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 33 of file RPCEfficiency.cc.

References cscSegments, debug, dt4DSegments, dupli, EffRootFileName, EffSaveRootFile, folderPath, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), inclcsc, incldt, incldtMB4, MaxD, MaxDrb4, MinCosAng, rangestrips, rangestripsRB4, RPCRecHitLabel_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                          {
  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.);


  cscSegments= consumes<CSCSegmentCollection>(iConfig.getParameter<edm::InputTag>("cscSegments"));
  dt4DSegments =  consumes<DTRecSegment4DCollection>(iConfig.getParameter<edm::InputTag>("dt4DSegments"));
  RPCRecHitLabel_ =  consumes<RPCRecHitCollection>(iConfig.getParameter<edm::InputTag>("RecHitLabel"));

  folderPath=iConfig.getUntrackedParameter<std::string>("folderPath","RPC/RPCEfficiency/");
   
  EffSaveRootFile  = iConfig.getUntrackedParameter<bool>("EffSaveRootFile", false); 
  EffRootFileName  = iConfig.getUntrackedParameter<std::string>("EffRootFileName", "RPCEfficiency.root"); 

}

RPCEfficiency::~RPCEfficiency ( )

override

Definition at line 218 of file RPCEfficiency.cc.

{}

Member Function Documentation

void RPCEfficiency::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprotected

Definition at line 220 of file RPCEfficiency.cc.

References CSCDetId::chamber(), CSCGeometry::chamber(), cscSegments, distsector_tmp(), dt4DSegments, dupli, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, CSCDetId::endcap(), MonitorElement::Fill(), edm::EventSetup::get(), edm::Event::getByToken(), hGlobalResClu1La, hGlobalResClu1R2A, hGlobalResClu1R2B, hGlobalResClu1R2C, hGlobalResClu1R3A, hGlobalResClu1R3B, hGlobalResClu1R3C, hGlobalResClu2La, hGlobalResClu2R2A, hGlobalResClu2R2B, hGlobalResClu2R2C, hGlobalResClu2R3A, hGlobalResClu2R3B, hGlobalResClu2R3C, hGlobalResClu3La, hGlobalResClu3R2A, hGlobalResClu3R2B, hGlobalResClu3R2C, hGlobalResClu3R3A, hGlobalResClu3R3B, hGlobalResClu3R3C, RPCRoll::id(), CSCChamber::id(), DTGeometry::idToDet(), inclcsc, incldt, incldtMB4, edm::HandleBase::isValid(), RPCDetId::layer(), RectangularStripTopology::localPosition(), TrapezoidalStripTopology::localPosition(), LogDebug, MaxD, MaxDrb4, meCollection, MinCosAng, RPCGeomServ::name(), RPCRoll::nstrips(), RectangularStripTopology::pitch(), TrapezoidalStripTopology::pitch(), rangestrips, DetId::rawId(), rpcPointValidation_cfi::recHit, RPCDetId::ring(), CSCDetId::ring(), RPCGeometry::roll(), RPCDetId::roll(), rollstoreCSC, rollstoreDT, RPCRecHitLabel_, DTChamberId::sector(), mathSSE::sqrt(), DTChamberId::station(), CSCDetId::station(), RPCDetId::station(), statistics, AlCaHLTBitMon_QueryRunRegistry::string, RPCRoll::strip(), RectangularStripTopology::stripLength(), TrapezoidalStripTopology::stripLength(), GeomDet::surface(), GeomDet::toGlobal(), GeomDet::toLocal(), RPCRoll::topology(), DTChamberId::wheel(), X, PV3DBase< T, PVType, FrameType >::x(), TrackerOfflineValidation_Dqm_cff::xmax, TrackerOfflineValidation_Dqm_cff::xmin, DOFs::Y, PV3DBase< T, PVType, FrameType >::y(), DOFs::Z, and PV3DBase< T, PVType, FrameType >::z().

                                                                              {
   

  edm::ESHandle<RPCGeometry> rpcGeo;
  edm::ESHandle<DTGeometry> dtGeo;  
  edm::ESHandle<CSCGeometry> cscGeo;
  
  iSetup.get<MuonGeometryRecord>().get(rpcGeo);
  iSetup.get<MuonGeometryRecord>().get(dtGeo);
  iSetup.get<MuonGeometryRecord>().get(cscGeo);
  
  statistics->Fill(1);
  
  std::stringstream  meIdRPC, meIdDT, meIdCSC;
  
   LogDebug("rpcefficiency") <<"\t Getting the RPC RecHits";
  edm::Handle<RPCRecHitCollection> rpcHits;
  iEvent.getByToken(RPCRecHitLabel_,rpcHits);  
  
  if(!rpcHits.isValid()) return;
  
  if(incldt){
     LogDebug("rpcefficiency")<<"\t Getting the DT Segments";
    edm::Handle<DTRecSegment4DCollection> all4DSegments;
    
    iEvent.getByToken(dt4DSegments, all4DSegments);
    
    if(all4DSegments.isValid()){
      
      if(all4DSegments->size()>0){
    
    if(all4DSegments->size()<=16) {statistics->Fill(2);}
    
    LogDebug("rpcefficiency")<<"\t Number of DT Segments in this event = "<<all4DSegments->size();
    
    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);
      
     LogDebug("rpcefficiency")<<"\t Loop over all the 4D Segments";
    for (segment = all4DSegments->begin(); segment != all4DSegments->end(); ++segment){ 
      
      DTChamberId DTId = segment->chamberId();
      
      
      if(DTSegmentCounter[DTId]==1 && DTId.station()!=4){   
        
        int dtWheel = DTId.wheel();
        int dtStation = DTId.station();
        int dtSector = DTId.sector();      
        
        LocalPoint segmentPosition= segment->localPosition();
        LocalVector segmentDirection=segment->localDirection();
        
        const GeomDet* gdet=dtGeo->idToDet(segment->geographicalId());
        const BoundPlane & DTSurface = gdet->surface();
        
        //check if the dimension of the segment is 4 
        
        if(segment->dimension()==4){
          
          float Xo=segmentPosition.x();
          float Yo=segmentPosition.y();
          float Zo=segmentPosition.z();
          float dx=segmentDirection.x();
          float dy=segmentDirection.y();
          float dz=segmentDirection.z();
          
          std::set<RPCDetId> rollsForThisDT = rollstoreDT[DTStationIndex(0,dtWheel,dtSector,dtStation)];
          
           LogDebug("rpcefficiency")<<"DT  \t \t Loop over all the rolls asociated to this DT";
          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(); 
        
        GlobalPoint CenterPointRollGlobal = RPCSurface.toGlobal(LocalPoint(0,0,0));
        
        LocalPoint CenterRollinDTFrame = DTSurface.toLocal(CenterPointRollGlobal);
        
        float D=CenterRollinDTFrame.z();
        
        float X=Xo+dx*D/dz;
        float Y=Yo+dy*D/dz;
        float Z=D;
        
        const RectangularStripTopology* top_= dynamic_cast<const RectangularStripTopology*> (&(rollasociated->topology()));
        LocalPoint xmin = top_->localPosition(0.);
         LogDebug("rpcefficiency")<<"DT  \t \t \t xmin of this  Roll "<<xmin<<"cm";
        LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
         LogDebug("rpcefficiency")<<"DT  \t \t \t xmax of this  Roll "<<xmax<<"cm";
        float rsize = fabs( xmax.x()-xmin.x() );
         LogDebug("rpcefficiency")<<"DT  \t \t \t Roll Size "<<rsize<<"cm";
        float stripl = top_->stripLength();
        float stripw = top_->pitch();
        
        float extrapolatedDistance = sqrt((X-Xo)*(X-Xo)+(Y-Yo)*(Y-Yo)+(Z-Zo)*(Z-Zo));
        
        if(extrapolatedDistance<=MaxD){ 
          
          GlobalPoint GlobalPointExtrapolated = DTSurface.toGlobal(LocalPoint(X,Y,Z));
          LocalPoint PointExtrapolatedRPCFrame = RPCSurface.toLocal(GlobalPointExtrapolated);
          
          if(fabs(PointExtrapolatedRPCFrame.z()) < 10. && 
             fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 && 
             fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
            
            RPCDetId  rollId = rollasociated->id();           
            RPCGeomServ rpcsrv(rollId);
            std::string nameRoll = rpcsrv.name();
             LogDebug("rpcefficiency")<<"DT  \t \t \t \t The RPCName is "<<nameRoll;            
            const float stripPredicted = 
              rollasociated->strip(LocalPoint(PointExtrapolatedRPCFrame.x(),PointExtrapolatedRPCFrame.y(),0.)); 
            
             LogDebug("rpcefficiency")<<"DT  \t \t \t \t Candidate (from DT Segment) STRIP---> "<<stripPredicted<< std::endl;         
            //---- HISTOGRAM STRIP PREDICTED FROM DT ----
            
            std::map<std::string, MonitorElement*> meMap=meCollection[rpcId.rawId()];
            meIdDT.str("");
            meIdDT<<"ExpectedOccupancyFromDT_"<<rollId.rawId();
            meMap[meIdDT.str()]->Fill(stripPredicted);
            //-----------------------------------------------------
              
            
            //-------RecHitPart Just For Residual--------
            int countRecHits = 0;
            int cluSize = 0;
            float minres = 3000.;
            
            typedef std::pair<RPCRecHitCollection::const_iterator, RPCRecHitCollection::const_iterator> rangeRecHits;
            rangeRecHits recHitCollection =  rpcHits->get(rollasociated->id());
            RPCRecHitCollection::const_iterator recHit;
              
            for (recHit = recHitCollection.first; recHit != recHitCollection.second ; recHit++) {
              countRecHits++;
              
              LocalPoint recHitPos=recHit->localPosition();
              float res=PointExtrapolatedRPCFrame.x()- recHitPos.x();       
               LogDebug("rpcefficiency")<<"DT  \t \t \t \t \t Found Rec Hit at "<<res<<"cm of the prediction.";
              if(fabs(res)<fabs(minres)){
            minres=res;
            cluSize = recHit->clusterSize();
             LogDebug("rpcefficiency")<<"DT  \t \t \t \t \t \t New Min Res "<<res<<"cm.";
              }
            }
            
            if(countRecHits==0){
               LogDebug("rpcefficiency") <<"DT \t \t \t \t \t THIS ROLL DOESN'T HAVE ANY RECHIT";
            }else{
              assert(minres!=3000);     
                
              if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
             LogDebug("rpcefficiency")<<"DT  \t \t \t \t \t \t True!";
            
            if(rollId.station()==1&&rollId.layer()==1)     { 
              if(cluSize==1*dupli) {hGlobalResClu1La[0]->Fill(minres);}
              else if(cluSize==2*dupli){ hGlobalResClu2La[0]->Fill(minres);} 
              else if(cluSize==3*dupli){ hGlobalResClu3La[0]->Fill(minres);}
            }else if(rollId.station()==1&&rollId.layer()==2){ 
              if(cluSize==1*dupli) {hGlobalResClu1La[1]->Fill(minres);}
              else if(cluSize==2*dupli){ hGlobalResClu2La[1]->Fill(minres);} 
              else if(cluSize==3*dupli){ hGlobalResClu3La[1]->Fill(minres);}
            }else if(rollId.station()==2&&rollId.layer()==1){ 
              if(cluSize==1*dupli) {hGlobalResClu1La[2]->Fill(minres);}
              else if(cluSize==2*dupli){ hGlobalResClu2La[2]->Fill(minres);} 
              else if(cluSize==3*dupli){ hGlobalResClu3La[2]->Fill(minres);}
            }else if(rollId.station()==2&&rollId.layer()==2){ 
              if(cluSize==1*dupli) {hGlobalResClu1La[3]->Fill(minres);}
              if(cluSize==2*dupli){ hGlobalResClu2La[3]->Fill(minres);} 
              else if(cluSize==3*dupli){ hGlobalResClu3La[3]->Fill(minres);}
            }else if(rollId.station()==3){ 
              if(cluSize==1*dupli) {hGlobalResClu1La[4]->Fill(minres);}
              else if(cluSize==2*dupli){ hGlobalResClu2La[4]->Fill(minres);} 
              else if(cluSize==3*dupli){ hGlobalResClu3La[4]->Fill(minres);}
            }
            meIdRPC.str("");
            meIdRPC<<"RPCDataOccupancyFromDT_"<<rollId.rawId();
            meMap[meIdRPC.str()]->Fill(stripPredicted);
              }
            }
          }else{
             LogDebug("rpcefficiency")<<"DT \t \t \t \t No the prediction is outside of this roll";
          }//Condition for the right match
        }else{
           LogDebug("rpcefficiency")<<"DT \t \t \t No, Exrtrapolation too long!, canceled";
          }//D so big
          }//loop over all the rolls asociated
        }//Is the segment 4D?
      }else {
         LogDebug("rpcefficiency")<<"DT \t \t More than one segment in this chamber, or we are in Station 4";
      }
    }
      } else {  
     LogDebug("rpcefficiency")<<"DT This Event doesn't have any DT4DDSegment"; //is ther more than 1 segment in this event?
      }
    }
  }
  
  if(incldtMB4){
    edm::Handle<DTRecSegment4DCollection> all4DSegments;
    iEvent.getByToken(dt4DSegments, all4DSegments);
      
    if(all4DSegments.isValid() && all4DSegments->size()>0){

    std::map<DTChamberId,int> DTSegmentCounter;
    DTRecSegment4DCollection::const_iterator segment;  
        
    for (segment = all4DSegments->begin();segment!=all4DSegments->end(); ++segment){
      DTSegmentCounter[segment->chamberId()]++;
    }    
    
     LogDebug("rpcefficiency")<<"MB4 \t \t Loop Over all4DSegments";
    for (segment = all4DSegments->begin(); segment != all4DSegments->end(); ++segment){ 
      
      DTChamberId DTId = segment->chamberId();
        
      if(DTSegmentCounter[DTId] == 1 && DTId.station()==4){
              int dtWheel = DTId.wheel();
          int dtStation = DTId.station();
          int dtSector = DTId.sector();
          
          LocalPoint segmentPosition= segment->localPosition();
          LocalVector segmentDirection=segment->localDirection();
          
          //check if the dimension of the segment is 2 and the station is 4

          if(segment->dimension()==2){
            const LocalVector& segmentDirectionMB4=segmentDirection;
            const LocalPoint& segmentPositionMB4=segmentPosition;
            
            
            const BoundPlane& DTSurface4 = dtGeo->idToDet(DTId)->surface();
            
            DTRecSegment4DCollection::const_iterator segMB3;  
            
            for(segMB3=all4DSegments->begin();segMB3!=all4DSegments->end();++segMB3){
              
              DTChamberId dtid3 = segMB3->chamberId();  
            
              if(distsector_tmp(dtid3.sector(),DTId.sector())<=1 
             && dtid3.station()==3
             && dtid3.wheel()==DTId.wheel()
             && DTSegmentCounter[dtid3] == 1
             && segMB3->dimension()==4){
            
            const GeomDet* gdet3=dtGeo->idToDet(segMB3->geographicalId());
            const BoundPlane & DTSurface3 = gdet3->surface();
            
            LocalVector segmentDirectionMB3 =  segMB3->localDirection();
            GlobalPoint segmentPositionMB3inGlobal = DTSurface3.toGlobal(segMB3->localPosition());
            
            GlobalVector segDirMB4inGlobalFrame=DTSurface4.toGlobal(segmentDirectionMB4);
            GlobalVector segDirMB3inGlobalFrame=DTSurface3.toGlobal(segmentDirectionMB3);
            
            float dx=segDirMB4inGlobalFrame.x();
            float dy=segDirMB4inGlobalFrame.y();
            //          float dz=segDirMB4inGlobalFrame.z();
            
            float dx3=segDirMB3inGlobalFrame.x();
            float dy3=segDirMB3inGlobalFrame.y();
            //  float dz3=segDirMB3inGlobalFrame.z();
            
            double cosAng=fabs(dx*dx3+dy*dy3/sqrt((dx3*dx3+dy3*dy3)*(dx*dx+dy*dy)));
            
            if(cosAng>MinCosAng){
              if(dtSector==13){
                dtSector=4;
              }
              if(dtSector==14){
                dtSector=10;
              }
              
              std::set<RPCDetId> rollsForThisDT = rollstoreDT[DTStationIndex(0,dtWheel,dtSector,dtStation)]; //It should be always 4
              
              assert(!rollsForThisDT.empty());
              
              for (std::set<RPCDetId>::iterator iteraRoll=rollsForThisDT.begin();iteraRoll != rollsForThisDT.end(); iteraRoll++){
                const RPCRoll* rollasociated = rpcGeo->roll(*iteraRoll); //roll asociado a MB4
                const BoundPlane & RPCSurfaceRB4 = rollasociated->surface(); //surface MB4
                
                //   RPCGeomServ rpcsrv(rpcId);
                //          std::string nameRoll = rpcsrv.name();
                
                GlobalPoint CenterPointRollGlobal=RPCSurfaceRB4.toGlobal(LocalPoint(0,0,0));
                LocalPoint CenterRollinMB4Frame = DTSurface4.toLocal(CenterPointRollGlobal); //In MB4
                LocalPoint segmentPositionMB3inMB4Frame = DTSurface4.toLocal(segmentPositionMB3inGlobal); //In MB4
                LocalVector segmentDirectionMB3inMB4Frame = DTSurface4.toLocal(segDirMB3inGlobalFrame); //In MB4
                
                //The exptrapolation is done in MB4 frame. for local x and z is done from MB4,
                float Dxz=CenterRollinMB4Frame.z();
                float Xo4=segmentPositionMB4.x();
                float dxl=segmentDirectionMB4.x(); //dx local for MB4 segment in MB4 Frame
                float dzl=segmentDirectionMB4.z(); //dx local for MB4 segment in MB4 Frame
                
                float X=Xo4+dxl*Dxz/dzl; //In MB4 frame
                float Z=Dxz;//In MB4 frame
                
                //for local y is done from MB3
                float Yo34=segmentPositionMB3inMB4Frame.y();
                float dy34 = segmentDirectionMB3inMB4Frame.y();
                float dz34 = segmentDirectionMB3inMB4Frame.z();
                float Dy=Dxz-(segmentPositionMB3inMB4Frame.z()); //Distance beetween the segment in MB3 and the RB4 surface
                
                float Y=Yo34+dy34*Dy/dz34;//In MB4 Frame
                
                const RectangularStripTopology* top_
                  =dynamic_cast<const RectangularStripTopology*>(&(rollasociated->topology())); //Topology roll asociated MB4
                LocalPoint xmin = top_->localPosition(0.);
                LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
                float rsize = fabs( xmax.x()-xmin.x() );
                float stripl = top_->stripLength();
                float stripw = top_->pitch();
                
                float extrapolatedDistance = sqrt((Y-Yo34)*(Y-Yo34)+Dy*Dy);
                
                if(extrapolatedDistance<=MaxDrb4){ 
                  
                  GlobalPoint GlobalPointExtrapolated = DTSurface4.toGlobal(LocalPoint(X,Y,Z));
                  LocalPoint PointExtrapolatedRPCFrame = RPCSurfaceRB4.toLocal(GlobalPointExtrapolated);
                  
                  if(fabs(PointExtrapolatedRPCFrame.z()) < 5.  &&
                 fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 &&
                 fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){
                
                RPCDetId  rollId = rollasociated->id();
                const float stripPredicted=rollasociated->strip(LocalPoint(PointExtrapolatedRPCFrame.x(),PointExtrapolatedRPCFrame.y(),0.)); 
                
                LogDebug("rpcefficiency")<<"MB4 \t \t \t \t Candidate (from DT Segment) STRIP---> "<<stripPredicted<< std::endl;
                //--------- HISTOGRAM STRIP PREDICTED FROM DT  MB4 -------------------
                
                std::map<std::string, MonitorElement*> meMap=meCollection[rollId.rawId()];              
                meIdDT.str("");
                meIdDT<<"ExpectedOccupancyFromDT_"<<rollId.rawId();
                meMap[meIdDT.str()]->Fill(stripPredicted);
                
                //-------RecHitPart Just For Residual--------
                int countRecHits = 0;
                int cluSize = 0;
                float minres = 3000.;
                
                LogDebug("rpcefficiency")<<"MB4 \t \t \t \t Getting RecHits in Roll Asociated";
                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();       
                  LogDebug("rpcefficiency")<<"DT  \t \t \t \t \t Found Rec Hit at "<<res<<"cm of the prediction.";
                  if(fabs(res)<fabs(minres)){
                    minres=res;
                    cluSize = recHit->clusterSize();
                  }
                }       
                
                if(countRecHits==0){
                  LogDebug("rpcefficiency") <<"MB4 \t \t \t \t \t \t THIS ROLL DOESN'T HAVE ANY RECHIT";
                }else{     
                  assert(minres!=3000); 
                  
                  if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
                    assert(rollId.station()==4);
                    if(cluSize==1*dupli){ hGlobalResClu1La[5]->Fill(minres);}
                    else if(cluSize==2*dupli){ hGlobalResClu2La[5]->Fill(minres);}
                    else if(cluSize==3*dupli){ hGlobalResClu3La[5]->Fill(minres);}
                    
                    meIdRPC.str("");
                    meIdRPC<<"RPCDataOccupancyFromDT_"<<rollId.rawId();
                    meMap[meIdRPC.str()]->Fill(stripPredicted);
                  }
                }
                  }else{
                 LogDebug("rpcefficiency")<<"MB4 \t \t \t \t No the prediction is outside of this roll";
                  }
                }//Condition for the right match
                else{
                   LogDebug("rpcefficiency")<<"MB4 \t \t \t No, Exrtrapolation too long!, canceled";
                }
              }//loop over all the rollsasociated
            }else{
               LogDebug("rpcefficiency")<<"MB4 \t \t \t \t I found segments in MB4 and MB3 adjacent or same wheel and sector but not compatibles Diferent Directions";
            }
              }else{//if dtid3.station()==3&&dtid3.sector()==DTId.sector()&&dtid3.wheel()==DTId.wheel()&&segMB3->dim()==4
             LogDebug("rpcefficiency")<<"MB4 \t \t \t No the same station or same wheel or segment dim in mb3 not 4D";
              }
            }//loop over all the segments looking for one in MB3 
          }else{
             LogDebug("rpcefficiency")<<"MB4 \t \t \t Is NOT a 2D Segment";
          }
      }else{
         LogDebug("rpcefficiency")<<"MB4 \t \t \t \t There is not just one segment or is not in station 4";
      }//De aca para abajo esta en dtpart.inl
    }
      }else{
     LogDebug("rpcefficiency")<<"MB4 \t This event doesn't have 4D Segment";
      }
  }
  
  
  if(inclcsc){
    LogDebug("rpcefficiency") <<"\t Getting the CSC Segments";
    edm::Handle<CSCSegmentCollection> allCSCSegments;
    
    iEvent.getByToken(cscSegments, allCSCSegments);
    
    if(allCSCSegments.isValid()){ 
      if(allCSCSegments->size()>0){
    statistics->Fill(18);
    
    LogDebug("rpcefficiency")<<"CSC \t Number of CSC Segments in this event = "<<allCSCSegments->size();
    
    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);
    
     LogDebug("rpcefficiency")<<"CSC \t loop over all the CSCSegments ";
    for (segment = allCSCSegments->begin();segment!=allCSCSegments->end(); ++segment){
      CSCDetId CSCId = segment->cscDetId();
      
      if(CSCSegmentsCounter[CSCId]==1 && CSCId.ring()!=1 && allCSCSegments->size()>=2){
         LogDebug("rpcefficiency")<<"CSC \t \t yes";
        int cscEndCap = CSCId.endcap();
        int cscStation = CSCId.station();
        int cscRing = CSCId.ring();
        //      int cscChamber = CSCId.chamber();
        int rpcRegion = 1; if(cscEndCap==2) rpcRegion= -1;//Relation among the 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();
        
         LogDebug("rpcefficiency")<<"CSC \t \t Is a good Segment? dim = 4, 4 <= nRecHits <= 10 Incident angle int range 45 < "<<acos(dz)*180/3.1415926<<" < 135? ";
        
        if(segment->dimension()==4 && (segment->nRecHits()<=10 && segment->nRecHits()>=4)&& acos(dz)*180/3.1415926 > 45. && acos(dz)*180/3.1415926 < 160. ){ 
          
          float Xo=segmentPosition.x();
          float Yo=segmentPosition.y();
          float Zo=segmentPosition.z();
          float dx=segmentDirection.x();
          float dy=segmentDirection.y();
          float dz=segmentDirection.z();
          
          
           LogDebug("rpcefficiency")<<"CSC \t \t Getting chamber from Geometry";
          const CSCChamber* TheChamber=cscGeo->chamber(CSCId); 
           LogDebug("rpcefficiency")<<"CSC \t \t Getting ID from Chamber";
          const CSCDetId TheId=TheChamber->id();
           LogDebug("rpcefficiency")<<"CSC \t \t Printing The Id"<<TheId;
          std::set<RPCDetId> rollsForThisCSC = rollstoreCSC[CSCStationIndex(rpcRegion,rpcStation,rpcRing,rpcSegment)];
           LogDebug("rpcefficiency")<<"CSC \t \t Number of rolls for this CSC = "<<rollsForThisCSC.size();
          
          if(rpcRing!=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();
          
          const BoundPlane & RPCSurface = rollasociated->surface(); 
          
          GlobalPoint CenterPointRollGlobal = RPCSurface.toGlobal(LocalPoint(0,0,0));
          LocalPoint CenterRollinCSCFrame = TheChamber->toLocal(CenterPointRollGlobal);
          
          float D=CenterRollinCSCFrame.z();
          
          float X=Xo+dx*D/dz;
          float Y=Yo+dy*D/dz;
          float Z=D;
          
          const TrapezoidalStripTopology* top_=dynamic_cast<const TrapezoidalStripTopology*>(&(rollasociated->topology()));
          LocalPoint xmin = top_->localPosition(0.);
           LogDebug("rpcefficiency")<<"CSC \t \t \t xmin of this  Roll "<<xmin<<"cm";
          LocalPoint xmax = top_->localPosition((float)rollasociated->nstrips());
           LogDebug("rpcefficiency")<<"CSC \t \t \t xmax of this  Roll "<<xmax<<"cm";
          float rsize = fabs( xmax.x()-xmin.x() );
           LogDebug("rpcefficiency")<<"CSC \t \t \t Roll Size "<<rsize<<"cm";
          float stripl = top_->stripLength();
          float stripw = top_->pitch();
          
          
          float extrapolatedDistance = sqrt((X-Xo)*(X-Xo)+(Y-Yo)*(Y-Yo)+(Z-Zo)*(Z-Zo));
          
          
          if(extrapolatedDistance<=MaxD){ 
            
            GlobalPoint GlobalPointExtrapolated=TheChamber->toGlobal(LocalPoint(X,Y,Z));
            LocalPoint PointExtrapolatedRPCFrame = RPCSurface.toLocal(GlobalPointExtrapolated);
            
            
            if(fabs(PointExtrapolatedRPCFrame.z()) < 10. && 
               fabs(PointExtrapolatedRPCFrame.x()) < rsize*0.5 && 
               fabs(PointExtrapolatedRPCFrame.y()) < stripl*0.5){ 
              
              RPCDetId  rollId = rollasociated->id();
              RPCGeomServ rpcsrv(rollId);
              std::string nameRoll = rpcsrv.name();
              
               LogDebug("rpcefficiency")<<"CSC \t \t \t \t The RPCName is "<<nameRoll;
              
              const float stripPredicted = 
            rollasociated->strip(LocalPoint(PointExtrapolatedRPCFrame.x(),PointExtrapolatedRPCFrame.y(),0.)); 
              
               LogDebug("rpcefficiency")<<"CSC  \t \t \t \t \t Candidate"<<rollId<<" "<<"(from CSC Segment) STRIP---> "<<stripPredicted<< std::endl;
              //--------- HISTOGRAM STRIP PREDICTED FROM CSC  -------------------
              
              std::map<std::string, MonitorElement*> meMap=meCollection[rpcId.rawId()];
              meIdCSC.str("");
              meIdCSC<<"ExpectedOccupancyFromCSC_"<<rollId.rawId();
              meMap[meIdCSC.str()]->Fill(stripPredicted);
              //--------------------------------------------------------------------
              
              
              //-------RecHitPart Just For Residual--------
              int cluSize = 0;
              int countRecHits = 0;
              float minres = 3000.;
              
               LogDebug("rpcefficiency")<<"CSC  \t \t \t \t \t Getting RecHits in Roll Asociated";
              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();
             LogDebug("rpcefficiency")<<"CSC  \t \t \t \t \t \t Found Rec Hit at "<<res<<"cm of the prediction.";
            if(fabs(res)<fabs(minres)){
              minres=res;
              cluSize = recHit->clusterSize();
               LogDebug("rpcefficiency")<<"CSC  \t \t \t \t \t \t \t New Min Res "<<res<<"cm.";
            }
              }
              
              if(countRecHits==0){
             LogDebug("rpcefficiency") <<"CSC \t \t \t \t \t THIS ROLL DOESN'T HAVE ANY RECHIT";
              }else{  
            assert(minres!=3000); 
            
            if(fabs(minres)<=(rangestrips+cluSize*0.5)*stripw){
               LogDebug("rpcefficiency")<<"CSC  \t \t \t \t \t \t True!";
              
              if(rollId.ring()==2&&rollId.roll()==1){
                if(cluSize==1*dupli) { hGlobalResClu1R2A->Fill(minres);} 
                else if(cluSize==2*dupli) { hGlobalResClu2R2A->Fill(minres); }
                else if(cluSize==3*dupli)  {hGlobalResClu3R2A->Fill(minres);}
              }
              else if(rollId.ring()==2&&rollId.roll()==2){
                if(cluSize==1*dupli) { hGlobalResClu1R2B->Fill(minres);} 
                else if(cluSize==2*dupli) { hGlobalResClu2R2B->Fill(minres); }
                else if(cluSize==3*dupli)  {hGlobalResClu3R2B->Fill(minres);}
              }
              else if(rollId.ring()==2&&rollId.roll()==3){
                if(cluSize==1*dupli) { hGlobalResClu1R2C->Fill(minres);} 
                else if(cluSize==2*dupli)  {hGlobalResClu2R2C->Fill(minres); }
                else if(cluSize==3*dupli) { hGlobalResClu3R2C->Fill(minres);}
              }
              else if(rollId.ring()==3&&rollId.roll()==1){
                if(cluSize==1*dupli)  {hGlobalResClu1R3A->Fill(minres); }
                else if(cluSize==2*dupli) { hGlobalResClu2R3A->Fill(minres); }
                else if(cluSize==3*dupli)  {hGlobalResClu3R3A->Fill(minres);}
              }
              else if(rollId.ring()==3&&rollId.roll()==2){
                if(cluSize==1*dupli) { hGlobalResClu1R3B->Fill(minres); }
                else if(cluSize==2*dupli) { hGlobalResClu2R3B->Fill(minres); }
                else if(cluSize==3*dupli) { hGlobalResClu3R3B->Fill(minres);}
              }
              else if(rollId.ring()==3&&rollId.roll()==3){
                if(cluSize==1*dupli) {hGlobalResClu1R3C->Fill(minres); }
                else if(cluSize==2*dupli) {hGlobalResClu2R3C->Fill(minres); }
                else if(cluSize==3*dupli) {hGlobalResClu3R3C->Fill(minres);}
              }
              meIdRPC.str("");
              meIdRPC<<"RPCDataOccupancyFromCSC_"<<rollId.rawId();
              meMap[meIdRPC.str()]->Fill(stripPredicted);
            }
              }

            }else{
               LogDebug("rpcefficiency")<<"CSC \t \t \t \t No the prediction is outside of this roll";
            }//Condition for the right match
          }else{//if extrapolation distance D is not too long
             LogDebug("rpcefficiency")<<"CSC \t \t \t No, Exrtrapolation too long!, canceled";
          }//D so big
        }//loop over the rolls asociated 
          }//Condition over the LS1 geometry!!!!
        }//Is the segment 4D?
      }else{
         LogDebug("rpcefficiency")<<"CSC \t \t More than one segment in this chamber, or we are in Station Ring 1 or in Station 4";
      }
    }
      }else{
     LogDebug("rpcefficiency")<<"CSC This Event doesn't have any CSCSegment";
      }
    }
    }
    
}

void RPCEfficiency::bookDetUnitSeg ( DQMStore::IBooker &  ibooker, RPCDetId &  detId, int  nstrips, std::string  folder, std::map< std::string, MonitorElement * > &  meMap  )

protected

Definition at line 5 of file RPCEfficiencyBook.cc.

References DQMStore::IBooker::book1D(), printsummarytable::folder, RPCBookFolderStructure::folderStructure(), DetId::rawId(), RPCDetId::region(), DQMStore::IBooker::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by bookHistograms().

                                                                                                                                                           {
  
  
  RPCBookFolderStructure folderStr;

  std::string folder = folderPath +  folderStr.folderStructure(detId);



  ibooker.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] = ibooker.book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
    
    sprintf(meId,"RPCDataOccupancyFromDT_%d",rawId);
    sprintf(meTitle,"RPCDataOccupancyFromDT_for_%d",rawId);
    meMap[meId] = ibooker.book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
    
  }else{  //Begin booking CSC

    sprintf(meId,"ExpectedOccupancyFromCSC_%d",rawId);
    sprintf(meTitle,"ExpectedOccupancyFromCSC_for_%d",rawId);
    meMap[meId] = ibooker.book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
    
    sprintf(meId,"RPCDataOccupancyFromCSC_%d",rawId);
    sprintf(meTitle,"RPCDataOccupancyFromCSC_for_%d",rawId);
    meMap[meId] = ibooker.book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
    
  }

}

void RPCEfficiency::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  run, edm::EventSetup const &  iSetup  )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 59 of file RPCEfficiency.cc.

References DQMStore::IBooker::book1D(), bookDetUnitSeg(), DQMStore::IBooker::cd(), RPCGeometry::dets(), printsummarytable::folder, folderPath, edm::EventSetup::get(), hGlobalResClu1La, hGlobalResClu1R2A, hGlobalResClu1R2B, hGlobalResClu1R2C, hGlobalResClu1R3A, hGlobalResClu1R3B, hGlobalResClu1R3C, hGlobalResClu2La, hGlobalResClu2R2A, hGlobalResClu2R2B, hGlobalResClu2R2C, hGlobalResClu2R3A, hGlobalResClu2R3B, hGlobalResClu2R3C, hGlobalResClu3La, hGlobalResClu3R2A, hGlobalResClu3R2B, hGlobalResClu3R2C, hGlobalResClu3R3A, hGlobalResClu3R3B, hGlobalResClu3R3C, mps_fire::i, inclcsc, incldt, incldtMB4, LogDebug, meCollection, nullptr, alignCSCRings::r, DetId::rawId(), RPCDetId::region(), relativeConstraints::ring, RPCDetId::ring(), RPCChamber::rolls(), rollstoreCSC, rollstoreDT, RPCDetId::sector(), RPCGeomServ::segment(), MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), relativeConstraints::station, RPCDetId::station(), statistics, AlCaHLTBitMon_QueryRunRegistry::string, and makeMuonMisalignmentScenario::wheel.

                                                                                                               {

  std::string folder;
  ibooker.cd();
  ibooker.setCurrentFolder(folderPath);
  statistics = ibooker.book1D("Statistics","All Statistics",33,0.5,33.5);

   
  statistics->setBinLabel(1,"Events ",1);
  statistics->setBinLabel(2,"Events with DT seg",1);
  
  std::stringstream sstr;
  for( int i = 1; i<=15; i++ ){ //DT form bin 3 to bin 17
    sstr.str("");
    sstr<<i<<" DT seg";
    statistics->setBinLabel(i+2, sstr.str() ,1);
  }

  statistics->setBinLabel(18,"Events with CSC seg",1);
  for( int i = 1; i<=15; i++ ){ //CSC form bin 19 to bin 33
    sstr.str("");
    sstr<<i<<" CSC seg";
    statistics->setBinLabel(i+18, sstr.str(),1);
  }
  

   LogDebug("rpcefficiency")<<"booking Global histograms with "<<folderPath;
   
  folder = folderPath+"MuonSegEff/Residuals/Barrel";
  ibooker.setCurrentFolder(folder);
 
  //Barrel
  std::stringstream histoName, histoTitle;

  for (int layer = 1 ; layer<= 6 ;layer++){
    histoName.str("");
    histoTitle.str("");
    histoName<<"GlobalResidualsClu1La"<<layer;
    histoTitle<<"RPC Residuals Layer "<<layer<<" Cluster Size 1"; 
    hGlobalResClu1La[layer-1] = ibooker.book1D(histoName.str(), histoTitle.str(),101,-10.,10.);
 
    histoName.str("");
    histoTitle.str("");
    histoName<<"GlobalResidualsClu2La"<<layer;
    histoTitle<<"RPC Residuals Layer "<<layer<<" Cluster Size 2"; 
    hGlobalResClu2La[layer-1] = ibooker.book1D(histoName.str(), histoTitle.str(),101,-10.,10.);
    
    histoName.str("");
    histoTitle.str("");
    histoName<<"GlobalResidualsClu3La"<<layer;
    histoTitle<<"RPC Residuals Layer "<<layer<<" Cluster Size 3"; 
    hGlobalResClu3La[layer-1] = ibooker.book1D(histoName.str(), histoTitle.str(),101,-10.,10.);
    
  }
  
  LogDebug("rpcefficiency")<<"Booking Residuals for EndCap";
  folder = folderPath+"MuonSegEff/Residuals/Endcap";
  ibooker.setCurrentFolder(folder);
  
  //Endcap   
  hGlobalResClu1R3C = ibooker.book1D("GlobalResidualsClu1R3C","RPC Residuals Ring 3 Roll C Cluster Size 1",101,-10.,10.);
  hGlobalResClu1R3B = ibooker.book1D("GlobalResidualsClu1R3B","RPC Residuals Ring 3 Roll B Cluster Size 1",101,-10.,10.);
  hGlobalResClu1R3A = ibooker.book1D("GlobalResidualsClu1R3A","RPC Residuals Ring 3 Roll A Cluster Size 1",101,-10.,10.);
  hGlobalResClu1R2C = ibooker.book1D("GlobalResidualsClu1R2C","RPC Residuals Ring 2 Roll C Cluster Size 1",101,-10.,10.);
  hGlobalResClu1R2B = ibooker.book1D("GlobalResidualsClu1R2B","RPC Residuals Ring 2 Roll B Cluster Size 1",101,-10.,10.);
  hGlobalResClu1R2A = ibooker.book1D("GlobalResidualsClu1R2A","RPC Residuals Ring 2 Roll A Cluster Size 1",101,-10.,10.);

  hGlobalResClu2R3C = ibooker.book1D("GlobalResidualsClu2R3C","RPC Residuals Ring 3 Roll C Cluster Size 2",101,-10.,10.);
  hGlobalResClu2R3B = ibooker.book1D("GlobalResidualsClu2R3B","RPC Residuals Ring 3 Roll B Cluster Size 2",101,-10.,10.);
  hGlobalResClu2R3A = ibooker.book1D("GlobalResidualsClu2R3A","RPC Residuals Ring 3 Roll A Cluster Size 2",101,-10.,10.);
  hGlobalResClu2R2C = ibooker.book1D("GlobalResidualsClu2R2C","RPC Residuals Ring 2 Roll C Cluster Size 2",101,-10.,10.);
  hGlobalResClu2R2B = ibooker.book1D("GlobalResidualsClu2R2B","RPC Residuals Ring 2 Roll B Cluster Size 2",101,-10.,10.);
  hGlobalResClu2R2A = ibooker.book1D("GlobalResidualsClu2R2A","RPC Residuals Ring 2 Roll A Cluster Size 2",101,-10.,10.);

  hGlobalResClu3R3C = ibooker.book1D("GlobalResidualsClu3R3C","RPC Residuals Ring 3 Roll C Cluster Size 3",101,-10.,10.);
  hGlobalResClu3R3B = ibooker.book1D("GlobalResidualsClu3R3B","RPC Residuals Ring 3 Roll B Cluster Size 3",101,-10.,10.);
  hGlobalResClu3R3A = ibooker.book1D("GlobalResidualsClu3R3A","RPC Residuals Ring 3 Roll A Cluster Size 3",101,-10.,10.);
  hGlobalResClu3R2C = ibooker.book1D("GlobalResidualsClu3R2C","RPC Residuals Ring 2 Roll C Cluster Size 3",101,-10.,10.);
  hGlobalResClu3R2B = ibooker.book1D("GlobalResidualsClu3R2B","RPC Residuals Ring 2 Roll B Cluster Size 3",101,-10.,10.);
  hGlobalResClu3R2A = ibooker.book1D("GlobalResidualsClu3R2A","RPC Residuals Ring 2 Roll A Cluster Size 3",101,-10.,10.);


  edm::ESHandle<RPCGeometry> rpcGeo;
  iSetup.get<MuonGeometryRecord>().get(rpcGeo);
  
  
  for (TrackingGeometry::DetContainer::const_iterator it=rpcGeo->dets().begin();it<rpcGeo->dets().end();it++){ //Loop on all detector units

    if(dynamic_cast< const RPCChamber* >( *it ) != nullptr ){ // check if chamber exists
      const RPCChamber* ch = dynamic_cast< const RPCChamber* >( *it ); 
      std::vector< const RPCRoll*> roles = (ch->rolls()); //get all rolls in a chambers

      for(std::vector<const RPCRoll*>::const_iterator r = roles.begin();r != roles.end(); ++r){ //Loop on all rolls
    
    RPCDetId rpcId = (*r)->id();
    int region=rpcId.region();
    
    LogDebug("rpcefficiency")<<"Booking for "<<rpcId.rawId();

    bookDetUnitSeg(ibooker, rpcId,(*r)->nstrips(),  folderPath+"MuonSegEff/", meCollection[rpcId.rawId()] ); // book histograms

    if(region==0&&(incldt||incldtMB4)){
      //LogDebug("rpcefficiency")<<"--Filling the dtstore"<<rpcId;
      int wheel=rpcId.ring();
      int sector=rpcId.sector();
      int station=rpcId.station();
      DTStationIndex ind(region,wheel,sector,station);
      std::set<RPCDetId> myrolls;
      if (rollstoreDT.find(ind)!=rollstoreDT.end()) {myrolls=rollstoreDT[ind];}
      myrolls.insert(rpcId);
      rollstoreDT[ind]=myrolls;
 
    }else if(region!=0 && inclcsc){
      int station=rpcId.station();
          int ring=rpcId.ring();
          int cscring=ring;
          int cscstation=station;
      RPCGeomServ rpcsrv(rpcId);
      int rpcsegment = rpcsrv.segment();
      int cscchamber = rpcsegment;

          if((station==2||station==3||station==4)&&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;

      if(rpcId.ring()==2 || rpcId.ring()==3){
        
        cscchamber = rpcsegment+1;                                                                                     
        if(cscchamber==37){cscchamber=1;}                                                                                    
        CSCStationIndex ind2(region,cscstation,cscring,cscchamber);                                                         
        std::set<RPCDetId> myrolls2;                                                                                        
        if (rollstoreCSC.find(ind2)!=rollstoreCSC.end()){myrolls2=rollstoreCSC[ind2];}                                          
        myrolls2.insert(rpcId);                                                                                             
        rollstoreCSC[ind2]=myrolls2;                                                                                         
        
        cscchamber = rpcsegment-1;                                                                                         
        if(cscchamber==0){cscchamber=36;}                                                                                    
        CSCStationIndex ind3(region,cscstation,cscring,cscchamber);                                                      
        std::set<RPCDetId> myrolls3;                                                                                     
        if (rollstoreCSC.find(ind3)!=rollstoreCSC.end()) {myrolls3=rollstoreCSC[ind3]; }                                
        myrolls3.insert(rpcId);                                                                                          
        rollstoreCSC[ind3]=myrolls3;    
        
      }
      
    }

      }
    }
  }
}//beginRun

Member Data Documentation

std::vector<std::map<RPCDetId, int> > RPCEfficiency::counter

private

Definition at line 139 of file RPCEfficiency.h.

Referenced by counter.Counters::__getitem__().

edm::EDGetTokenT<CSCSegmentCollection> RPCEfficiency::cscSegments

private

Definition at line 155 of file RPCEfficiency.h.

Referenced by analyze(), and RPCEfficiency().

bool RPCEfficiency::debug

private

Definition at line 145 of file RPCEfficiency.h.

Referenced by util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), rrapi.RRApi::get(), util.rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), RPCEfficiency(), and pkg.AbstractPkg::write().

edm::EDGetTokenT<DTRecSegment4DCollection> RPCEfficiency::dt4DSegments

private

Definition at line 156 of file RPCEfficiency.h.

Referenced by analyze(), and RPCEfficiency().

int RPCEfficiency::dupli

private

Definition at line 153 of file RPCEfficiency.h.

Referenced by analyze(), and RPCEfficiency().

std::string RPCEfficiency::EffRootFileName

private

Definition at line 165 of file RPCEfficiency.h.

Referenced by RPCEfficiency().

bool RPCEfficiency::EffSaveRootFile

private

Definition at line 164 of file RPCEfficiency.h.

Referenced by RPCEfficiency().

std::string RPCEfficiency::folderPath

private

Definition at line 159 of file RPCEfficiency.h.

Referenced by bookHistograms(), and RPCEfficiency().

MonitorElement* RPCEfficiency::hGlobalResClu1La[6]

protected

Definition at line 112 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu1R2A

protected

Definition at line 122 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu1R2B

protected

Definition at line 121 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu1R2C

protected

Definition at line 120 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu1R3A

protected

Definition at line 119 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu1R3B

protected

Definition at line 118 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu1R3C

protected

Definition at line 117 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu2La[6]

protected

Definition at line 113 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu2R2A

protected

Definition at line 129 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu2R2B

protected

Definition at line 128 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu2R2C

protected

Definition at line 127 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu2R3A

protected

Definition at line 126 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu2R3B

protected

Definition at line 125 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu2R3C

protected

Definition at line 124 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu3La[6]

protected

Definition at line 114 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu3R2A

protected

Definition at line 136 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu3R2B

protected

Definition at line 135 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu3R2C

protected

Definition at line 134 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu3R3A

protected

Definition at line 133 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu3R3B

protected

Definition at line 132 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

MonitorElement* RPCEfficiency::hGlobalResClu3R3C

protected

Definition at line 131 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

bool RPCEfficiency::inclcsc

private

Definition at line 144 of file RPCEfficiency.h.

Referenced by analyze(), bookHistograms(), and RPCEfficiency().

bool RPCEfficiency::incldt

private

Definition at line 142 of file RPCEfficiency.h.

Referenced by analyze(), bookHistograms(), and RPCEfficiency().

bool RPCEfficiency::incldtMB4

private

Definition at line 143 of file RPCEfficiency.h.

Referenced by analyze(), bookHistograms(), and RPCEfficiency().

bool RPCEfficiency::inves

private

Definition at line 147 of file RPCEfficiency.h.

double RPCEfficiency::MaxD

private

Definition at line 151 of file RPCEfficiency.h.

Referenced by analyze(), and RPCEfficiency().

double RPCEfficiency::MaxDrb4

private

Definition at line 152 of file RPCEfficiency.h.

Referenced by analyze(), and RPCEfficiency().

std::map<int, std::map<std::string, MonitorElement*> > RPCEfficiency::meCollection

private

Definition at line 162 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

double RPCEfficiency::MinCosAng

private

Definition at line 150 of file RPCEfficiency.h.

Referenced by analyze(), and RPCEfficiency().

double RPCEfficiency::rangestrips

private

Definition at line 148 of file RPCEfficiency.h.

Referenced by analyze(), and RPCEfficiency().

double RPCEfficiency::rangestripsRB4

private

Definition at line 149 of file RPCEfficiency.h.

Referenced by RPCEfficiency().

std::string RPCEfficiency::rollseff

private

Definition at line 160 of file RPCEfficiency.h.

std::map<CSCStationIndex,std::set<RPCDetId> > RPCEfficiency::rollstoreCSC

protected

Definition at line 107 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

std::map<DTStationIndex,std::set<RPCDetId> > RPCEfficiency::rollstoreDT

protected

Definition at line 106 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

edm::EDGetTokenT<RPCRecHitCollection> RPCEfficiency::RPCRecHitLabel_

private

Definition at line 157 of file RPCEfficiency.h.

Referenced by analyze(), and RPCEfficiency().

MonitorElement* RPCEfficiency::statistics

protected

Definition at line 109 of file RPCEfficiency.h.

Referenced by analyze(), and bookHistograms().

std::vector<int> RPCEfficiency::totalcounter

private

Definition at line 140 of file RPCEfficiency.h.