#include <RPCDigiValid.h>

Inheritance diagram for RPCDigiValid:

Public Member Functions
	RPCDigiValid (const edm::ParameterSet &ps)

	~RPCDigiValid ()

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

std::string	workerType () const

virtual	~EDAnalyzer ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)

void	beginJob ()

void	endJob (void)

Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Attributes
MonitorElement *	BxDist

DQMStore *	dbe_

std::string	digiLabel

std::string	outputFile_

MonitorElement *	Res

MonitorElement *	ResDmin1

MonitorElement *	ResDmin2

MonitorElement *	ResDmin3

MonitorElement *	ResDplu1

MonitorElement *	ResDplu2

MonitorElement *	ResDplu3

MonitorElement *	ResWmin1

MonitorElement *	ResWmin2

MonitorElement *	ResWplu1

MonitorElement *	ResWplu2

MonitorElement *	ResWzer0

MonitorElement *	rzview

MonitorElement *	StripProf

MonitorElement *	xyvDmin1

MonitorElement *	xyvDmin2

MonitorElement *	xyvDmin3

MonitorElement *	xyvDplu1

MonitorElement *	xyvDplu2

MonitorElement *	xyvDplu3

MonitorElement *	xyview

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

typedef WorkerT< EDAnalyzer >	WorkerType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 24 of file RPCDigiValid.h.

Constructor & Destructor Documentation

RPCDigiValid::RPCDigiValid ( const edm::ParameterSet & ps )

Definition at line 19 of file RPCDigiValid.cc.

References BxDist, dbe_, digiLabel, edm::ParameterSet::getUntrackedParameter(), cmsCodeRules.cppFunctionSkipper::operator, outputFile_, Res, ResDmin1, ResDmin2, ResDmin3, ResDplu1, ResDplu2, ResDplu3, ResWmin1, ResWmin2, ResWplu1, ResWplu2, ResWzer0, rzview, StripProf, xyvDmin1, xyvDmin2, xyvDmin3, xyvDplu1, xyvDplu2, xyvDplu3, and xyview.

                                                 :dbe_(0){
     
   digiLabel = ps.getUntrackedParameter<std::string>("digiLabel");
   outputFile_ = ps.getUntrackedParameter<string>("outputFile", "rpcDigiValidPlots.root");
   dbe_ = Service<DQMStore>().operator->();
   
   if ( dbe_ ) {
     dbe_->setCurrentFolder("RPCDigisV/RPCDigis");
     
     xyview = dbe_->book2D("X_Vs_Y_View","X_Vs_Y_View",1000, -760., 760., 1000, -760., 760.);
 
     xyvDplu1 = dbe_->book2D("Dplu1_XvsY","Dplu1_XvsY",1000, -760., 760., 1000, -760., 760.);
     xyvDplu2 = dbe_->book2D("Dplu2_XvsY","Dplu2_XvsY",1000, -760., 760., 1000, -760., 760.);
     xyvDplu3 = dbe_->book2D("Dplu3_XvsY","Dplu3_XvsY",1000, -760., 760., 1000, -760., 760.);
     xyvDmin1 = dbe_->book2D("Dmin1_XvsY","Dmin1_XvsY",1000, -760., 760., 1000, -760., 760.);
     xyvDmin2 = dbe_->book2D("Dmin2_XvsY","Dmin2_XvsY",1000, -760., 760., 1000, -760., 760.);
     xyvDmin3 = dbe_->book2D("Dmin3_XvsY","Dmin3_XvsY",1000, -760., 760., 1000, -760., 760.);
 
     rzview = dbe_->book2D("R_Vs_Z_View","R_Vs_Z_View",1000, -1100., 1100.,1000,0., 800.);
     Res  = dbe_->book1D("Digi_SimHit_difference", "Digi_SimHit_difference", 300, -8, 8);
     ResWmin2 = dbe_->book1D("W_Min2_Residuals", "W_Min2_Residuals", 300, -8, 8);
     ResWmin1 = dbe_->book1D("W_Min1_Residuals", "W_Min1_Residuals", 300, -8, 8);
     ResWzer0 = dbe_->book1D("W_Zer0_Residuals", "W_Zer0_Residuals", 300, -8, 8);
     ResWplu1 = dbe_->book1D("W_Plu1_Residuals", "W_Plu1_Residuals", 300, -8, 8);
     ResWplu2 = dbe_->book1D("W_Plu2_Residuals", "W_Plu2_Residuals", 300, -8, 8);
 
     BxDist = dbe_->book1D("Bunch_Crossing", "Bunch_Crossing", 20, -10., 10.);
     StripProf = dbe_->book1D("Strip_Profile", "Strip_Profile", 100, 0, 100);
 
     ResDmin1 = dbe_->book1D("Disk_Min1_Residuals", "Disk_Min1_Residuals", 300, -8, 8);
     ResDmin2 = dbe_->book1D("Disk_Min2_Residuals", "Disk_Min2_Residuals", 300, -8, 8);
     ResDmin3 = dbe_->book1D("Disk_Min3_Residuals", "Disk_Min3_Residuals", 300, -8, 8);
     ResDplu1 = dbe_->book1D("Disk_Plu1_Residuals", "Disk_Plu1_Residuals", 300, -8, 8);
     ResDplu2 = dbe_->book1D("Disk_Plu2_Residuals", "Disk_Plu2_Residuals", 300, -8, 8);
     ResDplu3 = dbe_->book1D("Disk_Plu3_Residuals", "Disk_Plu3_Residuals", 300, -8, 8);   
 
   }
 }

RPCDigiValid::~RPCDigiValid ( )

Definition at line 58 of file RPCDigiValid.cc.

58 {}

Member Function Documentation

void RPCDigiValid::analyze	(	const edm::Event &	e,
		const edm::EventSetup &	c
	)

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 66 of file RPCDigiValid.cc.

References BxDist, RPCRoll::centreOfStrip(), digiLabel, MonitorElement::Fill(), edm::EventSetup::get(), L1TEmulatorMonitor_cff::p, PV3DBase< T, PVType, FrameType >::perp(), RPCDetId::region(), Res, ResDmin1, ResDmin2, ResDmin3, ResDplu1, ResDplu2, ResDplu3, ResWmin1, ResWmin2, ResWplu1, ResWplu2, ResWzer0, RPCDetId::ring(), RPCDetId, rzview, RPCDetId::station(), StripProf, GeomDet::toGlobal(), x, PV3DBase< T, PVType, FrameType >::x(), xyvDmin1, xyvDmin2, xyvDmin3, xyvDplu1, xyvDplu2, xyvDplu3, xyview, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

                                                                           {
 
   //  cout << endl <<"--- [RPCDigiQuality] Analysing Event: #Run: " << event.id().run()
   //       << " #Event: " << event.id().event() << endl;
   
   // Get the RPC Geometry
   edm::ESHandle<RPCGeometry> rpcGeom;
   eventSetup.get<MuonGeometryRecord>().get(rpcGeom);
   
   edm::Handle<PSimHitContainer> simHit;
   event.getByLabel("g4SimHits", "MuonRPCHits", simHit);
   
   edm::Handle<RPCDigiCollection> rpcDigis;
   event.getByLabel(digiLabel, rpcDigis);
 
   // Loop on simhits
   PSimHitContainer::const_iterator simIt;
 
   //loop over Simhit
   std::map<RPCDetId, std::vector<double> > allsims;
 
   for (simIt = simHit->begin(); simIt != simHit->end(); simIt++) {
     RPCDetId Rsid = (RPCDetId)(*simIt).detUnitId();
     const RPCRoll* soll = dynamic_cast<const RPCRoll* >( rpcGeom->roll(Rsid));
     int ptype = simIt->particleType();
 
     //    std::cout <<"This is a Simhit with Parent "<<ptype<<std::endl;
     if (ptype == 13 || ptype == -13) {
       std::vector<double> buff;
       if (allsims.find(Rsid) != allsims.end() ){
     buff= allsims[Rsid];
       }
       buff.push_back(simIt->localPosition().x());
       allsims[Rsid]=buff;
     }
     GlobalPoint p=soll->toGlobal(simIt->localPosition());
     /*
     std::cout <<"Muon Position phi="<<p.phi()
         <<" R="<<p.perp()
           <<" z="<<p.z()<<std::endl;
     */
 
     double sim_x = p.x();
     double sim_y = p.y();
 
     xyview->Fill(sim_x, sim_y);
 
       if (Rsid.region() == 0 ){
      }
      else if (Rsid.region() == (+1)){
         if (Rsid.station() == 1)
           xyvDplu1->Fill(sim_x, sim_y);
         else if (Rsid.station() == 2)
           xyvDplu2->Fill(sim_x, sim_y);
         else if (Rsid.station() == 3)
           xyvDplu3->Fill(sim_x, sim_y);
       }
      else if (Rsid.region() == (-1)){
         if (Rsid.station() == 1)
           xyvDmin1->Fill(sim_x, sim_y);
         else if (Rsid.station() == 2)
           xyvDmin2->Fill(sim_x, sim_y);
         else if (Rsid.station() == 3)
           xyvDmin3->Fill(sim_x, sim_y);
       }
 
 
 //    xyview->Fill(p.x(),p.y());
     rzview->Fill(p.z(),p.perp());
 
   }
   //loop over Digis
   RPCDigiCollection::DigiRangeIterator detUnitIt;
   for (detUnitIt=rpcDigis->begin(); detUnitIt!=rpcDigis->end();++detUnitIt){
     const RPCDetId Rsid = (*detUnitIt).first;
     const RPCRoll* roll = dynamic_cast<const RPCRoll* >( rpcGeom->roll(Rsid));
     const RPCDigiCollection::Range& range = (*detUnitIt).second;
     std::vector<double> sims;
     if (allsims.find(Rsid) != allsims.end() ){
       sims = allsims[Rsid];
     }
     int ndigi=0;
     for (RPCDigiCollection::const_iterator digiIt = range.first;
      digiIt != range.second; ++digiIt){
       StripProf->Fill(digiIt->strip());
       BxDist->Fill(digiIt->bx());
       ndigi++;
     }
     
 
     //    std::cout<<" Number of Digi "<<ndigi<<" for "<<Rsid<<std::endl;
 
     if (sims.size() == 1 &&  ndigi == 1){
       double dis = roll->centreOfStrip(range.first->strip()).x()-sims[0];
       Res->Fill(dis);   
       
       if (Rsid.region() == 0 ){
     if (Rsid.ring() == -2)
       ResWmin2->Fill(dis);
     else if (Rsid.ring() == -1)
       ResWmin1->Fill(dis);
     else if (Rsid.ring() == 0)
       ResWzer0->Fill(dis);
     else if (Rsid.ring() == 1)
       ResWplu1->Fill(dis);
     else if (Rsid.ring() == 2)
       ResWplu2->Fill(dis);
       }
      else if (Rsid.region() == (+1)){
         if (Rsid.station() == 1)
           ResDplu1->Fill(dis);
         else if (Rsid.station() == 2)
           ResDplu2->Fill(dis);
         else if (Rsid.station() == 3)
           ResDplu3->Fill(dis);
       }
      else if (Rsid.region() == (-1)){
         if (Rsid.station() == 1)
           ResDmin1->Fill(dis);
         else if (Rsid.station() == 2)
           ResDmin2->Fill(dis);
         else if (Rsid.station() == 3)
           ResDmin3->Fill(dis);
       }
     }
   }
 }

void RPCDigiValid::beginJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 60 of file RPCDigiValid.cc.

60 {}

void RPCDigiValid::endJob ( void )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 62 of file RPCDigiValid.cc.

References dbe_, outputFile_, and DQMStore::save().

                           {
  if ( outputFile_.size() != 0 && dbe_ ) dbe_->save(outputFile_);
 }

Member Data Documentation

MonitorElement* RPCDigiValid::BxDist

private

Definition at line 46 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

DQMStore* RPCDigiValid::dbe_

private

Definition at line 66 of file RPCDigiValid.h.

Referenced by endJob(), and RPCDigiValid().

std::string RPCDigiValid::digiLabel

private

Definition at line 68 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

std::string RPCDigiValid::outputFile_

private

Definition at line 67 of file RPCDigiValid.h.

Referenced by endJob(), and RPCDigiValid().

MonitorElement* RPCDigiValid::Res

private

Definition at line 40 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResDmin1

private

Definition at line 50 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResDmin2

private

Definition at line 51 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResDmin3

private

Definition at line 52 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResDplu1

private

Definition at line 53 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResDplu2

private

Definition at line 54 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResDplu3

private

Definition at line 55 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResWmin1

private

Definition at line 42 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResWmin2

private

Definition at line 41 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResWplu1

private

Definition at line 44 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResWplu2

private

Definition at line 45 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::ResWzer0

private

Definition at line 43 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::rzview

private

Definition at line 39 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::StripProf

private

Definition at line 47 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::xyvDmin1

private

Definition at line 62 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::xyvDmin2

private

Definition at line 63 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::xyvDmin3

private

Definition at line 64 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::xyvDplu1

private

Definition at line 58 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::xyvDplu2

private

Definition at line 59 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::xyvDplu3

private

Definition at line 60 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::xyview

private

Definition at line 38 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation