#include <RPCEfficiencySecond.h>

Inheritance diagram for RPCEfficiencySecond:

Public Member Functions
int	rollY (std::string shortname, const std::vector< std::string > &rollNames)

	RPCEfficiencySecond (const edm::ParameterSet &)

	~RPCEfficiencySecond ()

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Public Attributes
MonitorElement *	DiskSummary [10]

MonitorElement *	EffDistroD [10]

MonitorElement *	EffDistroW [5]

MonitorElement *	EffGlobD [10]

MonitorElement *	EffGlobW [5]

MonitorElement *	ExGregR2D [10]

MonitorElement *	ExGregR3D [10]

MonitorElement *	ExpLayerW [5]

MonitorElement *	ExsectorEffW [5]

MonitorElement *	GregR2D [10]

MonitorElement *	GregR3D [10]

MonitorElement *	histoCSC

MonitorElement *	histoCSC2

MonitorElement *	histoDT

MonitorElement *	histoDT2

MonitorElement *	histoeffIdRPC_CSC

MonitorElement *	histoeffIdRPC_DT

MonitorElement *	histoPRO

MonitorElement *	histoRealRPC

MonitorElement *	histoRealRPC2

MonitorElement *	histoRPC

MonitorElement *	histoRPC2

MonitorElement *	ObsLayerW [5]

MonitorElement *	OcGregR2D [10]

MonitorElement *	OcGregR3D [10]

MonitorElement *	OcsectorEffW [5]

MonitorElement *	sectorEffW [5]

MonitorElement *	WheelSummary [5]

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginRun (const edm::Run &, const edm::EventSetup &iSetup)

std::map< std::string, MonitorElement * >	bookDetUnitSeg (RPCDetId &detId, int nstrips, std::string folder)

virtual void	endJob ()

virtual void	endRun (const edm::Run &, const edm::EventSetup &)

Private Attributes
DQMStore *	dbe

bool	debug

std::string	folderPath

int	innermostRings_

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

std::string	NameFile

int	numberOfDisks_

bool	SaveFile

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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)

Detailed Description

author Camilo Carrillo (Uniandes)

Definition at line 23 of file RPCEfficiencySecond.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 22 of file RPCEfficiencySecond.cc.

References debug, folderPath, edm::ParameterSet::getUntrackedParameter(), innermostRings_, NameFile, numberOfDisks_, SaveFile, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                       {
   SaveFile  = iConfig.getUntrackedParameter<bool>("SaveFile", false); 
   NameFile  = iConfig.getUntrackedParameter<std::string>("NameFile","RPCEfficiency.root"); 
   folderPath  = iConfig.getUntrackedParameter<std::string>("folderPath","RPC/RPCEfficiency/"); 
   debug = iConfig.getUntrackedParameter<bool>("debug",false); 
   numberOfDisks_ =   iConfig.getUntrackedParameter<int>("NumberOfEndcapDisks", 4);
   innermostRings_ = iConfig.getUntrackedParameter<int>("NumberOfInnermostEndcapRings", 2);
 
 }

RPCEfficiencySecond::~RPCEfficiencySecond ( )

Definition at line 32 of file RPCEfficiencySecond.cc.

32 {}

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 175 of file RPCEfficiencySecond.cc.

175 { }

void RPCEfficiencySecond::beginRun	(	const edm::Run &	,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 34 of file RPCEfficiencySecond.cc.

References DQMStore::book1D(), DQMStore::book2D(), dbe, DiskSummary, EffDistroD, EffDistroW, EffGlobD, EffGlobW, ExGregR2D, ExGregR3D, ExpLayerW, ExsectorEffW, folderPath, GregR2D, GregR3D, cmsHarvester::index, innermostRings_, rpcdqm::utils::labelXAxisSector(), rpcdqm::utils::labelXAxisSegment(), rpcdqm::utils::labelYAxisRing(), rpcdqm::utils::labelYAxisRoll(), numberOfDisks_, ObsLayerW, OcGregR2D, OcGregR3D, OcsectorEffW, cppFunctionSkipper::operator, sectorEffW, DQMStore::setCurrentFolder(), w, and WheelSummary.

                                                                             {
   
   dbe = edm::Service<DQMStore>().operator->();
   rpcdqm::utils rpcUtils;
   //Barrel 
   std::stringstream folderName;
   std::stringstream MeName;   std::stringstream MeTitle;
   for (int w =-2; w<=2; w++){  
 
     folderName.str("");
     folderName<<folderPath<<"Wheel_"<<w;
 
     dbe->setCurrentFolder(folderName.str());
 
     MeName.str("");
     MeName<<"EffDistroWheel_"<<w;
     MeTitle.str("");
     MeTitle<<"Efficiency Distribution for Wheel "<<w;
     EffDistroW[w+2]=dbe->book1D(MeName.str(),MeTitle.str(),20,0.5,100.5);
 
     MeName.str("");
     MeName<<"GlobEfficiencyWheel_"<<w;
     MeTitle.str("");
     MeTitle<<"Efficiency Wheel "<<w;
     EffGlobW[w+2]=dbe->book1D(MeName.str(),MeTitle.str(),206,0.5,206.5);
 
     dbe->setCurrentFolder(folderPath);
     MeName.str("");
     MeName<<"Efficiency_Roll_vs_Sector_Wheel_"<<w;                               
     WheelSummary[w+2] = dbe->book2D(MeName.str(),MeName.str(),12, 0.5,12.5, 21, 0.5, 21.5);
     rpcUtils.labelXAxisSector( WheelSummary[w+2] );
     rpcUtils.labelYAxisRoll( WheelSummary[w+2], 0, w, true);
 
     dbe->setCurrentFolder(folderPath+"Azimutal/");
     MeName.str("");
     MeName<<"AzimutalDistroW"<<w;
     MeTitle.str("");
     MeTitle<<"Efficiency per Sector Wheel "<<w;
     sectorEffW[w+2]= dbe->book1D(MeName.str(),MeTitle.str(),12,0.5,12.5);
 
     MeName.str("");
     MeName<<"AzimutalDistroW"<<w<<"Ex";
     MeTitle.str("");
     MeTitle<<"Expected per Sector Wheel "<<w;
     ExsectorEffW[w+2]= dbe->book1D(MeName.str(),MeTitle.str(),12,0.5,12.5);
    
     MeName.str("");
     MeName<<"AzimutalDistroW"<<w<<"Oc";
     MeTitle.str("");
     MeTitle<<"Occupancy per Sector Wheel "<<w;
     OcsectorEffW[w+2]= dbe->book1D(MeName.str(),MeTitle.str(),12,0.5,12.5);
 
     dbe->setCurrentFolder(folderPath+"BarrelPerLayer/");
     MeName.str("");
     MeName<<"ExpLayerW"<<w;
     MeTitle.str("");
     MeTitle<<"Expected Wheel "<<w;
     ExpLayerW[w+2]= dbe->book1D(MeName.str(),MeTitle.str(),6,0.5,6.5);
  
     MeName.str("");
     MeName<<"ObsLayerW"<<w;
     MeTitle.str("");
     MeTitle<<"Observed Wheel "<<w;
     ObsLayerW[w+2]= dbe->book1D(MeName.str(),MeTitle.str(),6,0.5,6.5);
   }
 
   //EndCap
   int index = 0;
   for (int d = (-1 *numberOfDisks_); d<=numberOfDisks_; d++){
     
     if (d==0) {continue;}
 
     folderName.str("");
     folderName<<folderPath<<"Disk_"<<d;
     dbe->setCurrentFolder(folderName.str());
 
     MeName.str("");
     MeName<<"EffDistroDisk_"<<d;
     MeTitle.str("");
     MeTitle<<"Efficiency Distribution Disk "<<d;
     EffDistroD[index]=dbe->book1D(MeName.str(),MeTitle.str(),20,0.5,100.5);
     
     MeName.str("");
     MeName<<"GlobEfficiencyDisk_"<<d;
     MeTitle.str("");
     MeTitle<<"Efficiency Disk "<<d;
     EffGlobD[index] = dbe->book1D(MeName.str(),MeTitle.str(),218,0.5,218.5);
 
 
     dbe->setCurrentFolder(folderPath);
     MeName.str("");
     MeName<<"Efficiency_Roll_vs_Segment_Disk_"<<d;
     DiskSummary[index] = dbe->book2D(MeName.str(),MeName.str(),36,0.5,36.5,6,0.5,6.5);
  
     rpcUtils.labelXAxisSegment(DiskSummary[index]);
     rpcUtils.labelYAxisRing(DiskSummary[index], innermostRings_ , true);
     
 
 
     dbe->setCurrentFolder(folderPath+"Azimutal/");
     MeName.str("");
     MeName<<"GregDistroR2D"<<d;
     MeTitle.str("");
     MeTitle<<"Efficiency for Station "<<d<<" Ring 2";
     GregR2D[index]= dbe->book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
 
     MeName.str("");
     MeName<<"GregDistroR3D"<<d;
     MeTitle.str("");
     MeTitle<<"Efficiency for Station "<<d<<" Ring 3";
     GregR3D[index]= dbe->book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
 
     MeName.str("");
     MeName<<"OcGregDistroR2D"<<d;
     MeTitle.str("");
     MeTitle<<"Occupancy Distribution for Station "<<d<<" Ring 2";
     OcGregR2D[index]= dbe->book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
 
     MeName.str("");
     MeName<<"OcGregDistroR3D"<<d;
     MeTitle.str("");
     MeTitle<<"Occupancy Distribution for Station "<<d<<" Ring 3";
     OcGregR3D[index]= dbe->book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
 
     MeName.str("");
     MeName<<"ExGregDistroR2D"<<d;
     MeTitle.str("");
     MeTitle<<"Expected Distribution for Station "<<d<<" Ring 2";
     ExGregR2D[index]= dbe->book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
 
 
     MeName.str("");
     MeName<<"ExGregDistroR3D"<<d;
     MeTitle.str("");
     MeTitle<<"Expected Distribution for Station "<<d<<" Ring 3";
     ExGregR3D[index]= dbe->book1D(MeName.str(),MeTitle.str(),36,0.5,36.5);
 
     index++;
   }
 }

std::map< std::string, MonitorElement * > RPCEfficiencySecond::bookDetUnitSeg	(	RPCDetId &	detId,
		int	nstrips,
		std::string	folder
	)

private

Definition at line 12 of file RPCEfficiencyBookSecondStep.cc.

References DQMStore::book1D(), dbe, RPCBookFolderStructure::folderStructure(), RPCGeomServ::name(), RPCDetId::region(), DQMStore::setCurrentFolder(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                                          {
   
   std::map<std::string, MonitorElement*> meMap;
    
   RPCBookFolderStructure *  folderStr = new RPCBookFolderStructure(); 
 
   std::string folder = folderPath+ "RollByRoll/" +  folderStr->folderStructure(detId);
 
   delete folderStr;
 
   dbe->setCurrentFolder(folder);
 
   RPCGeomServ RPCname(detId);
   std::string nameRoll = RPCname.name();
   char detUnitLabel[128];
   char layerLabel[128];
 
   sprintf(detUnitLabel ,"%s",nameRoll.c_str());
   sprintf(layerLabel ,"%s",nameRoll.c_str());
 
   char meId [128];
   char meTitle [128];
   
   //Begin booking DT
   if(detId.region()==0) {
 
     //std::cout<<"Booking "<<folder<<meId<<std::endl;
     sprintf(meId,"ExpectedOccupancyFromDT_%s",detUnitLabel);
     sprintf(meTitle,"ExpectedOccupancyFromDT_for_%s",layerLabel);
     meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
     //std::cout<<"Booking "<<meId<<std::endl;
  
     sprintf(meId,"RPCDataOccupancyFromDT_%s",detUnitLabel);
     sprintf(meTitle,"RPCDataOccupancyFromDT_for_%s",layerLabel);
     meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
 
     sprintf(meId,"Profile_%s",detUnitLabel);
     sprintf(meTitle,"Profile_for_%s",layerLabel);
     meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
     //std::cout<<"Booking "<<folder<<meId<<std::endl;
 
 //     sprintf(meId,"BXDistribution_%s",detUnitLabel);
 //     sprintf(meTitle,"BXDistribution_for_%s",layerLabel);
 //     meMap[meId] = dbe->book1D(meId, meTitle, 11,-5.5, 5.5);
     
   }else{
     //std::cout<<"Booking for the EndCap"<<detUnitLabel<<std::endl;
 
     //std::cout<<"Booking "<<meId<<std::endl;
     sprintf(meId,"ExpectedOccupancyFromCSC_%s",detUnitLabel);
     sprintf(meTitle,"ExpectedOccupancyFromCSC_for_%s",layerLabel);
     meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
        
     //std::cout<<"Booking "<<meId<<std::endl;
     sprintf(meId,"RPCDataOccupancyFromCSC_%s",detUnitLabel);
     sprintf(meTitle,"RPCDataOccupancyFromCSC_for_%s",layerLabel);
     meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
     
     //std::cout<<"Booking "<<meId<<std::endl;
     sprintf(meId,"Profile_%s",detUnitLabel);
     sprintf(meTitle,"Profile_for_%s",layerLabel);
     meMap[meId] = dbe->book1D(meId, meTitle, nstrips, 0.5, nstrips+0.5);
     
 //     //std::cout<<"Booking "<<meId<<std::endl;
 //     sprintf(meId,"BXDistribution_%s",detUnitLabel);
 //     sprintf(meTitle,"BXDistribution_for_%s",layerLabel);
 //     meMap[meId] = dbe->book1D(meId, meTitle, 11,-5.5, 5.5);
   }
   return meMap;
 }

void RPCEfficiencySecond::endJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 495 of file RPCEfficiencySecond.cc.

495 {}

void RPCEfficiencySecond::endRun	(	const edm::Run &	r,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 177 of file RPCEfficiencySecond.cc.

References dbe, rpcdqm::utils::detId2RollNr(), DiskSummary, EffDistroD, EffDistroW, EffGlobD, EffGlobW, ExGregR2D, ExGregR3D, ExpLayerW, ExsectorEffW, MonitorElement::Fill(), folderPath, edm::EventSetup::get(), DQMStore::get(), MonitorElement::getBinContent(), MonitorElement::getTH1F(), GregR2D, GregR3D, histoCSC, histoDT, histoRPC, i, j, gen::k, RPCDetId::layer(), LogDebug, N, RPCGeomServ::name(), NameFile, numberOfDisks_, python.connectstrParser::o, ObsLayerW, OcGregR2D, OcGregR3D, OcsectorEffW, AlCaHLTBitMon_ParallelJobs::p, alignCSCRings::r, DetId::rawId(), RPCDetId::region(), relativeConstraints::ring, RPCDetId::ring(), RPCChamber::rolls(), DQMStore::save(), SaveFile, RPCDetId::sector(), sectorEffW, RPCGeomServ::segment(), MonitorElement::setAxisRange(), MonitorElement::setAxisTitle(), MonitorElement::setBinContent(), MonitorElement::setBinError(), MonitorElement::setBinLabel(), RPCGeomServ::shortname(), mathSSE::sqrt(), relativeConstraints::station, RPCDetId::station(), AlCaHLTBitMon_QueryRunRegistry::string, w, WheelSummary, x, detailsBasic3DVector::y, and cuy::yBin.

                                                                             {
   
   LogDebug("rpcefficiencysecond")<<"Getting the RPC Geometry";
   
   edm::ESHandle<RPCGeometry> rpcGeo;
   iSetup.get<MuonGeometryRecord>().get(rpcGeo);
   
   
   //Setting Labels in Summary Label.
   std::stringstream binLabel;
   for (int w = -2; w<=2 ;w++){
     for(int x = 1;x<=12;x++){
       for(int y = 1;y<=21;y++){
         WheelSummary[w+2]->setBinContent(x,y,-1);
       }
     }
   }
  
   
   for (int d = 0 ; d<(numberOfDisks_*2); d++){
     for(int x = 1;x<=36;x++){
       for(int y = 1;y<=6;y++){
     DiskSummary[d]->setBinContent(x,y,-1);
       }
     }
   }
   
   
   binLabel.str("");
   int indexWheel[5];
   for(int j=0;j<5;j++){
     indexWheel[j]=0;
   }
    
   int indexDisk[10];
   for(int j=0;j<10;j++){
     indexDisk[j]=0;
   }
    
   
   for(TrackingGeometry::DetContainer::const_iterator it=rpcGeo->dets().begin(); it!=rpcGeo->dets().end();it++){
    
     
 
     if(dynamic_cast< const RPCChamber* >( *it ) != 0 ){
       const RPCChamber* ch = dynamic_cast< const RPCChamber* >( *it ); 
       std::vector< const RPCRoll*> roles = (ch->rolls());
       for(std::vector<const RPCRoll*>::const_iterator r = roles.begin();r != roles.end(); ++r){
     RPCDetId rpcId = (*r)->id();
     RPCGeomServ rpcsrv(rpcId);
     //  std::string nameRoll = rpcsrv.name();
     std::string camera = rpcsrv.name();
     
     //Breaking down the geometry
     int region = rpcId.region();
     int wheel = rpcId.ring();  int ring = rpcId.ring();  
     int sector = rpcId.sector();
     int station = rpcId.station();
     int geolayer = rpcId.layer(); 
       
     
     if(region==0){//Barrel
       
       std::stringstream meIdRPC,  meIdDT; //,  bxDistroId;    
       
       meIdRPC<<folderPath<<"MuonSegEff/RPCDataOccupancyFromDT_"<<rpcId.rawId();
       meIdDT<<folderPath<<"MuonSegEff/ExpectedOccupancyFromDT_"<<rpcId.rawId();
       
       histoRPC = dbe->get(meIdRPC.str());
       histoDT = dbe->get(meIdDT.str());
          
       int NumberWithOutPrediction=0;
       double p = 0.;
       double o = 0.;
       double ef = 0.;
       double er = 0.;
       double buffef = 0.;
       double buffer = 0.;
       double sumbuffef = 0.;
       double sumbuffer = 0.;
       double averageeff = 0.;
       int NumberStripsPointed = 0;
      
       
       if(histoRPC && histoDT){ 
 
         for(int i=1;i<=int((*r)->nstrips());++i){
       
           if(histoDT->getBinContent(i)!=0){
         LogDebug("rpcefficiencysecond")<<"Inside the If";
         
         buffef = double(histoRPC->getBinContent(i))/double(histoDT->getBinContent(i));
         buffer = sqrt(buffef*(1.-buffef)/double(histoDT->getBinContent(i)));
         
         sumbuffef=sumbuffef+buffef;
         sumbuffer = sumbuffer + buffer*buffer;
         NumberStripsPointed++;
           }else{
         NumberWithOutPrediction++;
           }
           LogDebug("rpcefficiencysecond")<<"Strip="<<i<<" RPC="<<histoRPC->getBinContent(i)<<" DT="<<histoDT->getBinContent(i)<<" buffef="<<buffef<<" buffer="<<buffer<<" sumbuffef="<<sumbuffef<<" sumbuffer="<<sumbuffer<<" NumberStripsPointed="<<NumberStripsPointed<<" NumberWithOutPrediction"<<NumberWithOutPrediction;
         }
 
         p=histoDT->getTH1F()->Integral();
         o=histoRPC->getTH1F()->Integral();
 
         if(NumberStripsPointed!=0){
           averageeff = (sumbuffef/double(NumberStripsPointed))*100.;
           //averageerr = sqrt(sumbuffer/double(NumberStripsPointed))*100.;
         }
 
       } 
       
       //      int Ring = rpcId.ring();
       
       if(p!=0){
         ef = double(o)/double(p); 
         er = sqrt(ef*(1.-ef)/double(p));
       }
       
       ef=ef*100;
       er=er*100;
 
       
        //Filling azimutal Wheel Histograms
           
       int layer = 0;
       
         if(station==1&&geolayer==1) layer = 1;
         else if(station==1&&geolayer==2) layer = 2;
         else if(station==2&&geolayer==1) layer = 3;
         else if(station==2&&geolayer==2)  layer = 4;
         else if(station==3) layer = 5;
         else if(station==4) layer = 6;
 
 
         ExsectorEffW[wheel+2]->Fill(sector,p); OcsectorEffW[wheel+2]->Fill(sector,o); 
         ExpLayerW[wheel+2]->Fill(layer, p); ObsLayerW[wheel+2]->Fill(layer, o);
        
                 
         
         LogDebug("rpcefficiencysecond")<<"Pigi "<<camera<<" "<<rpcsrv.shortname()<<" "<<(*r)->id();
         
         if(p > 100){//We need at least 100 predictions to fill the summary plot
           
           int xBin,yBin;
           xBin= (*r)->id().sector();
           rpcdqm::utils rollNumber;
           yBin = rollNumber.detId2RollNr((*r)->id());
         
         
           WheelSummary[wheel+2]->setBinContent(xBin,yBin,averageeff);
         }
         
         
         EffDistroW[wheel+2]->Fill(averageeff);
         indexWheel[wheel+2]++;
         EffGlobW[wheel+2]->setBinContent(indexWheel[wheel+2],ef);  
         EffGlobW[wheel+2]->setBinError(indexWheel[wheel+2],er);  
         EffGlobW[wheel+2]->setBinLabel(indexWheel[wheel+2],camera,1);
 
     }else{//EndCap
 
       
       std::stringstream meIdRPC,meIdCSC; //, bxDistroId;
       std::string      meIdPRO;
       
       
       meIdRPC<<folderPath<<"MuonSegEff/RPCDataOccupancyFromCSC_"<<rpcId.rawId();
       meIdCSC<<folderPath<<"MuonSegEff/ExpectedOccupancyFromCSC_"<<rpcId.rawId();
       
       meIdPRO = "Profile_"+ rpcsrv.name();
       
       histoRPC= dbe->get(meIdRPC.str());
       histoCSC= dbe->get(meIdCSC.str());
       //BXDistribution = dbe->get(bxDistroId.str());
       
       int NumberWithOutPrediction=0;
       double p = 0;
       double o = 0;
       double ef =0;
       double er =0;
       double buffef = 0;
       double buffer = 0;
       double sumbuffef = 0;
       double sumbuffer = 0;
       double averageeff = 0;
       int NumberStripsPointed = 0;
         
       
       if(histoRPC && histoCSC) {
         LogDebug("rpcefficiencysecond")<<rpcsrv.name();
         
         for(int i=1;i<=int((*r)->nstrips());++i){
           if(histoCSC->getBinContent(i)!=0){
         LogDebug("rpcefficiencysecond")<<"Inside the If";
         buffef = double(histoRPC->getBinContent(i))/double(histoCSC->getBinContent(i));  
         //  meMap[meIdPRO]->setBinContent(i,buffef); 
         buffer = sqrt(buffef*(1.-buffef)/double(histoCSC->getBinContent(i)));
         //  meMap[meIdPRO]->setBinError(i,buffer);
         sumbuffef=sumbuffef+buffef;   
         sumbuffer = sumbuffer + buffer*buffer;
         NumberStripsPointed++;
           }else{
         NumberWithOutPrediction++;
           }
           
           LogDebug("rpcefficiencysecond")<<"Strip="<<i<<" RPC="<<histoRPC->getBinContent(i)<<" CSC="<<histoCSC->getBinContent(i)<<" buffef="<<buffef<<" buffer="<<buffer<<" sumbuffef="<<sumbuffef<<" sumbuffer="<<sumbuffer<<" NumberStripsPointed="<<NumberStripsPointed<<" NumberWithOutPrediction"<<NumberWithOutPrediction;
         }
         
         p=histoCSC->getTH1F()->Integral();
         o=histoRPC->getTH1F()->Integral();
         
         if(NumberStripsPointed!=0){
           averageeff = (sumbuffef/double(NumberStripsPointed))*100.;
         }
         
       }
       
       int Disk = station*region;
 
       if(Disk > numberOfDisks_ || Disk<-numberOfDisks_){continue;} //remove strange disk numbers!!
       
       int dIndex = -1;
       if(region == -1){
         dIndex = Disk+numberOfDisks_;
       }else if(region == 1){
         dIndex = Disk+numberOfDisks_-1;
       }
 
       if( dIndex<0 || dIndex>= numberOfDisks_*2){continue;} //extra check on disk numering
       
       if(p!=0){
         ef = double(o)/double(p); 
         er = sqrt(ef*(1.-ef)/double(p));
       }
       
       ef=ef*100;
       er=er*100;
       
       if(ring==2){ 
         if (ExGregR2D[dIndex] && OcGregR2D[dIndex] ){
           ExGregR2D[dIndex]->Fill(rpcsrv.segment(),p);
           OcGregR2D[dIndex]->Fill(rpcsrv.segment(),o);
         }
       } else if(ring==3){
         ExGregR3D[dIndex]->Fill(rpcsrv.segment(),p);
         OcGregR3D[dIndex]->Fill(rpcsrv.segment(),o);
       }
           
       if(p > 100){ //We need at least 100 predictions to fill the summary plot
         RPCGeomServ RPCServ((*r)->id());
         int xBin = RPCServ.segment();
         int yBin= (ring-1)*3-(*r)->id().roll()+1;
         
         DiskSummary[dIndex]->setBinContent(xBin, yBin, averageeff);
       }
       
         
       if (  EffDistroD[dIndex]){EffDistroD[dIndex]->Fill(averageeff);}
       indexDisk[dIndex]++;     
       if ( EffGlobD[dIndex]){
         EffGlobD[dIndex]->setBinContent(indexDisk[dIndex],ef);     
         EffGlobD[dIndex]->setBinError(indexDisk[dIndex],er);  
         EffGlobD[dIndex]->setBinLabel(indexDisk[dIndex],camera,1);     
         
       }
     }
       }
     }
   }
   
   double eff,N,err;
   int k;
      
   for (int d = 0; d<(numberOfDisks_*2); d++){
     for(k=1;k<=36;k++){
          
       err=0; eff=0; N=ExGregR2D[d]->getBinContent(k);
       if(N!=0.){ eff = OcGregR2D[d]->getBinContent(k)/N; err=sqrt(eff*(1-eff)/N);}
       GregR2D[d]->setBinContent(k,eff); GregR2D[d]->setBinError(k,err);
       
       err=0; eff=0; N=ExGregR3D[d]->getBinContent(k);
       if(N!=0.){ eff = OcGregR3D[d]->getBinContent(k)/N; err=sqrt(eff*(1-eff)/N);}
       GregR3D[d]->setBinContent(k,eff); GregR3D[d]->setBinError(k,err);
     }
   }
      
   
   for (int w =-2; w<=2; w++){
     for(k=1;k<=12;k++){
       err=0; eff=0; N=ExsectorEffW[w+2]->getBinContent(k);
       if(N!=0.){ eff = OcsectorEffW[w+2]->getBinContent(k)/N;err=sqrt(eff*(1-eff)/N);}
       sectorEffW[w+2]->setBinContent(k,eff); sectorEffW[w+2]->setBinError(k,err);
     }
   }
      
   
   //Ranges for Both
   //Barrel
      
   for (int w=-2; w<=2; w++){
     EffGlobW[w+2]->setAxisRange(-4.,100.,2);
     EffGlobW[w+2]->setAxisTitle("%",2);
   }   
   
   for (int d=0; d<(numberOfDisks_*2); d++){
     EffGlobD[d]->setAxisRange(-4.,100.,2);
     EffGlobD[d]->setAxisTitle("%",2);
   }
   
      
   if(SaveFile){
     dbe->save(NameFile);
   }
   
 }

int RPCEfficiencySecond::rollY	(	std::string	shortname,
		const std::vector< std::string > &	rollNames
	)

Member Data Documentation

DQMStore* RPCEfficiencySecond::dbe

private

Definition at line 99 of file RPCEfficiencySecond.h.

Referenced by beginRun(), bookDetUnitSeg(), and endRun().

bool RPCEfficiencySecond::debug

private

Definition at line 93 of file RPCEfficiencySecond.h.

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

MonitorElement* RPCEfficiencySecond::DiskSummary[10]

Definition at line 61 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::EffDistroD[10]

Definition at line 57 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::EffDistroW[5]

Definition at line 51 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::EffGlobD[10]

Definition at line 54 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::EffGlobW[5]

Definition at line 48 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::ExGregR2D[10]

Definition at line 75 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::ExGregR3D[10]

Definition at line 76 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::ExpLayerW[5]

Definition at line 78 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::ExsectorEffW[5]

Definition at line 67 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

std::string RPCEfficiencySecond::folderPath

private

Definition at line 96 of file RPCEfficiencySecond.h.

Referenced by beginRun(), endRun(), and RPCEfficiencySecond().

MonitorElement* RPCEfficiencySecond::GregR2D[10]

Definition at line 69 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::GregR3D[10]

Definition at line 70 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::histoCSC

Definition at line 33 of file RPCEfficiencySecond.h.

Referenced by endRun().

MonitorElement* RPCEfficiencySecond::histoCSC2

Definition at line 43 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoDT

Definition at line 31 of file RPCEfficiencySecond.h.

Referenced by endRun().

MonitorElement* RPCEfficiencySecond::histoDT2

Definition at line 41 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoeffIdRPC_CSC

Definition at line 37 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoeffIdRPC_DT

Definition at line 36 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoPRO

Definition at line 35 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoRealRPC

Definition at line 32 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoRealRPC2

Definition at line 42 of file RPCEfficiencySecond.h.

MonitorElement* RPCEfficiencySecond::histoRPC

Definition at line 30 of file RPCEfficiencySecond.h.

Referenced by endRun().

MonitorElement* RPCEfficiencySecond::histoRPC2

Definition at line 40 of file RPCEfficiencySecond.h.

int RPCEfficiencySecond::innermostRings_

private

Definition at line 98 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and RPCEfficiencySecond().

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

private

Definition at line 91 of file RPCEfficiencySecond.h.

std::string RPCEfficiencySecond::NameFile

private

Definition at line 95 of file RPCEfficiencySecond.h.

Referenced by endRun(), and RPCEfficiencySecond().

int RPCEfficiencySecond::numberOfDisks_

private

Definition at line 97 of file RPCEfficiencySecond.h.

Referenced by beginRun(), endRun(), and RPCEfficiencySecond().

MonitorElement* RPCEfficiencySecond::ObsLayerW[5]

Definition at line 79 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::OcGregR2D[10]

Definition at line 72 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::OcGregR3D[10]

Definition at line 73 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::OcsectorEffW[5]

Definition at line 66 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

bool RPCEfficiencySecond::SaveFile

private

Definition at line 94 of file RPCEfficiencySecond.h.

Referenced by endRun(), and RPCEfficiencySecond().

MonitorElement* RPCEfficiencySecond::sectorEffW[5]

Definition at line 65 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

MonitorElement* RPCEfficiencySecond::WheelSummary[5]

Definition at line 60 of file RPCEfficiencySecond.h.

Referenced by beginRun(), and endRun().

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation