RPCDigiValid Class Reference

#include <RPCDigiValid.h>

Inheritance diagram for RPCDigiValid:

Public Member Functions
	RPCDigiValid (const edm::ParameterSet &ps)
	~RPCDigiValid ()
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
bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)
virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &c)
void	beginJob ()
void	beginRun (edm::Run const &, edm::EventSetup const &)
void	endJob (void)
void	endRun (edm::Run const &, edm::EventSetup const &)
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)

Private Attributes
MonitorElement *	BxDisc_4Min
MonitorElement *	BxDisc_4Plus
MonitorElement *	BxDist
MonitorElement *	CLS_Endcap_123_Ring3_A
MonitorElement *	CLS_Endcap_123_Ring3_B
MonitorElement *	CLS_Endcap_123_Ring3_C
MonitorElement *	CLS_Endcap_1_Ring2_A
MonitorElement *	CLS_Endcap_1_Ring2_B
MonitorElement *	CLS_Endcap_1_Ring2_C
MonitorElement *	CLS_Endcap_23_Ring2_A
MonitorElement *	CLS_Endcap_23_Ring2_B
MonitorElement *	CLS_Endcap_23_Ring2_C
MonitorElement *	CLS_Endcap_4
MonitorElement *	clsBarrel
int	countEvent
DQMStore *	dbe_
std::string	digiLabel
std::map< RPCDetId, double >	mapRollArea
std::map< RPCDetId, double >	mapRollCls
std::map< RPCDetId, int >	mapRollFakeCount
std::map< RPCDetId, std::string >	mapRollName
std::map< RPCDetId, std::map < int, double > * >	mapRollNoisyStripRate
std::map< RPCDetId, double >	mapRollStripArea
std::map< RPCDetId, std::map < int, double > * >	mapRollStripRate
std::map< RPCDetId, int >	mapRollTruCount
MonitorElement *	noiseCLS
MonitorElement *	noiseCLSBarrel
MonitorElement *	noiseCLSEndcaps
std::string	outputFile_
MonitorElement *	Res
MonitorElement *	Res_Endcap123_Ring3_A
MonitorElement *	Res_Endcap123_Ring3_B
MonitorElement *	Res_Endcap123_Ring3_C
MonitorElement *	Res_Endcap1_Ring2_A
MonitorElement *	Res_Endcap1_Ring2_B
MonitorElement *	Res_Endcap1_Ring2_C
MonitorElement *	Res_Endcap23_Ring2_A
MonitorElement *	Res_Endcap23_Ring2_B
MonitorElement *	Res_Endcap23_Ring2_C
MonitorElement *	ResDmin1
MonitorElement *	ResDmin2
MonitorElement *	ResDmin3
MonitorElement *	ResDmin4
MonitorElement *	ResDplu1
MonitorElement *	ResDplu2
MonitorElement *	ResDplu3
MonitorElement *	ResDplu4
MonitorElement *	ResLayer1_barrel
MonitorElement *	ResLayer2_barrel
MonitorElement *	ResLayer3_barrel
MonitorElement *	ResLayer4_barrel
MonitorElement *	ResLayer5_barrel
MonitorElement *	ResLayer6_barrel
MonitorElement *	ResWmin1
MonitorElement *	ResWmin2
MonitorElement *	ResWplu1
MonitorElement *	ResWplu2
MonitorElement *	ResWzer0
edm::EDGetTokenT < RPCDigiCollection >	rpcDigiToken
MonitorElement *	rzview
edm::EDGetTokenT < edm::PSimHitContainer >	simHitToken
MonitorElement *	StripProf
MonitorElement *	xyvDmin4
MonitorElement *	xyvDplu4
MonitorElement *	xyview

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 &)

Detailed Description

Definition at line 32 of file RPCDigiValid.h.

Constructor & Destructor Documentation

RPCDigiValid::RPCDigiValid

(

const edm::ParameterSet &

ps

)

Definition at line 18 of file RPCDigiValid.cc.

References dbe_, edm::ParameterSet::getUntrackedParameter(), cppFunctionSkipper::operator, outputFile_, rpcDigiToken, and simHitToken.

                                                 :
  dbe_(0)
{

//  Init the tokens for run data retrieval - stanislav
//  ps.getUntackedParameter<InputTag> retrieves a InputTag from the configuration. The second param is default value
//  module, instance and process labels may be passed in a single string if separated by colon ':'
//  (@see the edm::InputTag constructor documentation)
  simHitToken = consumes<PSimHitContainer>(ps.getUntrackedParameter<edm::InputTag >("simHitTag", edm::InputTag("g4SimHits:MuonRPCHits")));
  rpcDigiToken    = consumes<RPCDigiCollection>(ps.getUntrackedParameter<edm::InputTag>("rpcDigiTag", edm::InputTag("simMuonRPCDigis")));

  outputFile_ = ps.getUntrackedParameter<string> ("outputFile", "rpcDigiValidPlots.root");
  dbe_ = Service<DQMStore> ().operator->();
}

RPCDigiValid::~RPCDigiValid

(

)

Definition at line 33 of file RPCDigiValid.cc.

{
}

Member Function Documentation

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

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 38 of file RPCDigiValid.cc.

References funct::abs(), BxDisc_4Min, BxDisc_4Plus, BxDist, RPCRoll::centreOfStrip(), CLS_Endcap_123_Ring3_A, CLS_Endcap_123_Ring3_B, CLS_Endcap_123_Ring3_C, CLS_Endcap_1_Ring2_A, CLS_Endcap_1_Ring2_B, CLS_Endcap_1_Ring2_C, CLS_Endcap_23_Ring2_A, CLS_Endcap_23_Ring2_B, CLS_Endcap_23_Ring2_C, CLS_Endcap_4, clsBarrel, countEvent, MonitorElement::Fill(), edm::EventSetup::get(), RPCRoll::id(), RPCDetId::layer(), mapRollArea, mapRollCls, mapRollFakeCount, mapRollName, mapRollNoisyStripRate, mapRollStripArea, mapRollStripRate, mapRollTruCount, mergeVDriftHistosByStation::name, noiseCLS, noiseCLSBarrel, noiseCLSEndcaps, RPCRoll::nstrips(), AlCaHLTBitMon_ParallelJobs::p, PV3DBase< T, PVType, FrameType >::perp(), DetId::rawId(), RPCDetId::region(), Res, Res_Endcap123_Ring3_A, Res_Endcap123_Ring3_B, Res_Endcap123_Ring3_C, Res_Endcap1_Ring2_A, Res_Endcap1_Ring2_B, Res_Endcap1_Ring2_C, Res_Endcap23_Ring2_A, Res_Endcap23_Ring2_B, Res_Endcap23_Ring2_C, ResDmin1, ResDmin2, ResDmin3, ResDmin4, ResDplu1, ResDplu2, ResDplu3, ResDplu4, ResLayer1_barrel, ResLayer2_barrel, ResLayer3_barrel, ResLayer4_barrel, ResLayer5_barrel, ResLayer6_barrel, ResWmin1, ResWmin2, ResWplu1, ResWplu2, ResWzer0, RPCDetId::ring(), RPCDetId::roll(), RPCDetId, rpcDigiToken, rzview, simHitToken, RPCDetId::station(), AlCaHLTBitMon_QueryRunRegistry::string, StripProf, GeomDet::toGlobal(), RPCRoll::topology(), x, PV3DBase< T, PVType, FrameType >::x(), SiStripMonitorClusterAlca_cfi::xmax, SiStripMonitorClusterAlca_cfi::xmin, xyvDmin4, xyvDplu4, xyview, PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{

  countEvent++;
  //  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;
  edm::Handle<RPCDigiCollection> rpcDigis;
  event.getByToken(simHitToken, simHit);
  event.getByToken(rpcDigiToken, 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;

      // std::cout << "allsims[Rsid] = "  << std::endl;
    }
    GlobalPoint p = soll->toGlobal(simIt->localPosition());

    double sim_x = p.x();
    double sim_y = p.y();

    xyview->Fill(sim_x, sim_y);

    if (Rsid.region() == (+1))
    {
      if (Rsid.station() == 4)
      {
        xyvDplu4->Fill(sim_x, sim_y);
      }
    }
    else if (Rsid.region() == (-1))
    {
      if (Rsid.station() == 4)
      {
        xyvDmin4->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));

    RPCGeomServ rpcsrv(roll->id());
    std::string name = rpcsrv.name();
    //std:: cout << (roll->id().rawId()) << "\t" << name << std::endl;


    const RPCDigiCollection::Range& range = (*detUnitIt).second;
    std::vector<double> sims;
    if (allsims.find(Rsid) != allsims.end())
    {
      sims = allsims[Rsid];
    }

    //int ndigi=0;
    double ndigi = 0;
    for (RPCDigiCollection::const_iterator digiIt = range.first; digiIt != range.second; ++digiIt)
    {
      StripProf->Fill(digiIt->strip());
      BxDist->Fill(digiIt->bx());
      //bx for 4 endcaps
      if (Rsid.region() == (+1))
      {
        if (Rsid.station() == 4)
          BxDisc_4Plus->Fill(digiIt->bx());
      }
      else if (Rsid.region() == (-1))
      {
        if (Rsid.station() == 4)
          BxDisc_4Min->Fill(digiIt->bx());
      }
      map<int, double>* stripRate = mapRollStripRate[Rsid.rawId()];
      map<int, double>* stripNoisyRate = mapRollNoisyStripRate[Rsid.rawId()];
      //init map strip rate
      if (stripRate == 0)
      {
        stripRate = new map<int, double> ();
        mapRollStripRate[Rsid.rawId()] = stripRate;
      }
      (*stripRate)[digiIt->strip()] += 1;

      //noisy only
      if (sims.size() == 0)
      {
        if (stripNoisyRate == 0)
        {
          stripNoisyRate = new map<int, double> ();
          mapRollNoisyStripRate[Rsid.rawId()] = stripNoisyRate;
        }
        (*stripNoisyRate)[digiIt->strip()] += 1;

      }

      //      std::cout << Rsid.rawId() << "\tstrip = " << digiIt->strip() << std::endl;
      ndigi++;
      //      std::cout << "digis = " <<  ndigi << std::endl;
    }

    double area = 0.0;
    double stripArea = 0.0;

    if (Rsid.region() == 0)
    {
      const RectangularStripTopology* top_ = dynamic_cast<const RectangularStripTopology*> (&(roll->topology()));
      float xmin = (top_->localPosition(0.)).x();
      float xmax = (top_->localPosition((float) roll->nstrips())).x();
      float striplength = (top_->stripLength());
      area = striplength * (xmax - xmin);
      stripArea = area / ((float) roll->nstrips());
    }
    else
    {
      const TrapezoidalStripTopology* top_ = dynamic_cast<const TrapezoidalStripTopology*> (&(roll->topology()));
      float xmin = (top_->localPosition(0.)).x();
      float xmax = (top_->localPosition((float) roll->nstrips())).x();
      float striplength = (top_->stripLength());
      area = striplength * (xmax - xmin);
      stripArea = area / ((float) roll->nstrips());
    }

    mapRollTruCount[Rsid] += 1;

    if (sims.size() == 0)
    {
      noiseCLS->Fill(ndigi);
      if (Rsid.region() == 0)
      {
        noiseCLSBarrel->Fill(ndigi);
      }
      if (Rsid.region() == +1 || Rsid.region() == 1)
      {
        noiseCLSEndcaps->Fill(ndigi);
      }

      mapRollCls[Rsid] += ndigi;
      mapRollFakeCount[Rsid] += 1;
    }
    mapRollArea[Rsid] = area;
    mapRollStripArea[Rsid] = stripArea;
    mapRollName[Rsid] = name;

    //CLS histos
    if (Rsid.region() == 0)
    {
      clsBarrel->Fill(ndigi);
    }
    //endcap 
    if (Rsid.region() != 0)
    {
      if (Rsid.ring() == 2)
      {
        if (abs(Rsid.station()) == 1)
        {
          if (Rsid.roll() == 1)
            CLS_Endcap_1_Ring2_A->Fill(ndigi);
          if (Rsid.roll() == 2)
            CLS_Endcap_1_Ring2_B->Fill(ndigi);
          if (Rsid.roll() == 3)
            CLS_Endcap_1_Ring2_C->Fill(ndigi);
        }
        if (abs(Rsid.station()) == 2 || abs(Rsid.station()) == 3)
        {
          if (Rsid.roll() == 1)
            CLS_Endcap_23_Ring2_A->Fill(ndigi);
          if (Rsid.roll() == 2)
            CLS_Endcap_23_Ring2_B->Fill(ndigi);
          if (Rsid.roll() == 3)
            CLS_Endcap_23_Ring2_C->Fill(ndigi);
        }
      }
      if (Rsid.ring() == 3)
      {
        if (Rsid.roll() == 1)
          CLS_Endcap_123_Ring3_A->Fill(ndigi);
        if (Rsid.roll() == 2)
          CLS_Endcap_123_Ring3_B->Fill(ndigi);
        if (Rsid.roll() == 3)
          CLS_Endcap_123_Ring3_C->Fill(ndigi);
      }
      if (abs(Rsid.station()) == 4)
        CLS_Endcap_4->Fill(ndigi);
    }

    //cls histos
    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);
        //barrel layers
        if (Rsid.station() == 1 && Rsid.layer() == 1)
          ResLayer1_barrel->Fill(dis);
        if (Rsid.station() == 1 && Rsid.layer() == 2)
          ResLayer2_barrel->Fill(dis);
        if (Rsid.station() == 2 && Rsid.layer() == 1)
          ResLayer3_barrel->Fill(dis);
        if (Rsid.station() == 2 && Rsid.layer() == 2)
          ResLayer4_barrel->Fill(dis);
        if (Rsid.station() == 3)
          ResLayer5_barrel->Fill(dis);
        if (Rsid.station() == 4)
          ResLayer6_barrel->Fill(dis);
      }
      //endcap layers residuals
      if (Rsid.region() != 0)
      {
        if (Rsid.ring() == 2)
        {
          if (abs(Rsid.station()) == 1)
          {
            if (Rsid.roll() == 1)
              Res_Endcap1_Ring2_A->Fill(dis);
            if (Rsid.roll() == 2)
              Res_Endcap1_Ring2_B->Fill(dis);
            if (Rsid.roll() == 3)
              Res_Endcap1_Ring2_C->Fill(dis);
          }
          if (abs(Rsid.station()) == 2 || abs(Rsid.station()) == 3)
          {
            if (Rsid.roll() == 1)
              Res_Endcap23_Ring2_A->Fill(dis);
            if (Rsid.roll() == 2)
              Res_Endcap23_Ring2_B->Fill(dis);
            if (Rsid.roll() == 3)
              Res_Endcap23_Ring2_C->Fill(dis);
          }
        }
        if (Rsid.ring() == 3)
        {
          if (Rsid.roll() == 1)
            Res_Endcap123_Ring3_A->Fill(dis);
          if (Rsid.roll() == 2)
            Res_Endcap123_Ring3_B->Fill(dis);
          if (Rsid.roll() == 3)
            Res_Endcap123_Ring3_C->Fill(dis);
        }
      }

      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.station() == 4)
          ResDplu4->Fill(dis);
      }
      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);
        else if (Rsid.station() == 4)
          ResDmin4->Fill(dis);
      }
    }
  }
}

void RPCDigiValid::beginJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 350 of file RPCDigiValid.cc.

{
}

void RPCDigiValid::beginRun ( edm::Run const &  run, edm::EventSetup const &  eSetup  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 361 of file RPCDigiValid.cc.

References DQMStore::book1D(), DQMStore::book2D(), BxDisc_4Min, BxDisc_4Plus, BxDist, CLS_Endcap_123_Ring3_A, CLS_Endcap_123_Ring3_B, CLS_Endcap_123_Ring3_C, CLS_Endcap_1_Ring2_A, CLS_Endcap_1_Ring2_B, CLS_Endcap_1_Ring2_C, CLS_Endcap_23_Ring2_A, CLS_Endcap_23_Ring2_B, CLS_Endcap_23_Ring2_C, CLS_Endcap_4, clsBarrel, countEvent, dbe_, noiseCLS, noiseCLSBarrel, noiseCLSEndcaps, Res, Res_Endcap123_Ring3_A, Res_Endcap123_Ring3_B, Res_Endcap123_Ring3_C, Res_Endcap1_Ring2_A, Res_Endcap1_Ring2_B, Res_Endcap1_Ring2_C, Res_Endcap23_Ring2_A, Res_Endcap23_Ring2_B, Res_Endcap23_Ring2_C, ResDmin1, ResDmin2, ResDmin3, ResDmin4, ResDplu1, ResDplu2, ResDplu3, ResDplu4, ResLayer1_barrel, ResLayer2_barrel, ResLayer3_barrel, ResLayer4_barrel, ResLayer5_barrel, ResLayer6_barrel, ResWmin1, ResWmin2, ResWplu1, ResWplu2, ResWzer0, rzview, DQMStore::setCurrentFolder(), StripProf, xyvDmin4, xyvDplu4, and xyview.

{
  countEvent = 0;

  if (dbe_)
  {
    dbe_->setCurrentFolder("RPCDigisV/RPCDigis");

    xyview = dbe_->book2D("X_Vs_Y_View", "X_Vs_Y_View", 155, -775., 775., 155, -775., 775.);

    xyvDplu4 = dbe_->book2D("Dplu4_XvsY", "Dplu4_XvsY", 155, -775., 775., 155, -775., 775.);
    xyvDmin4 = dbe_->book2D("Dmin4_XvsY", "Dmin4_XvsY", 155, -775., 775., 155, -775., 775.);

    rzview = dbe_->book2D("R_Vs_Z_View", "R_Vs_Z_View", 216, -1080., 1080., 52, 260., 780.);
    Res = dbe_->book1D("Digi_SimHit_difference", "Digi_SimHit_difference", 300, -8, 8);
    ResWmin2 = dbe_->book1D("W_Min2_Residuals", "W_Min2_Residuals", 400, -8, 8);
    ResWmin1 = dbe_->book1D("W_Min1_Residuals", "W_Min1_Residuals", 400, -8, 8);
    ResWzer0 = dbe_->book1D("W_Zer0_Residuals", "W_Zer0_Residuals", 400, -8, 8);
    ResWplu1 = dbe_->book1D("W_Plu1_Residuals", "W_Plu1_Residuals", 400, -8, 8);
    ResWplu2 = dbe_->book1D("W_Plu2_Residuals", "W_Plu2_Residuals", 400, -8, 8);

    ResLayer1_barrel = dbe_->book1D("ResLayer1_barrel", "ResLayer1_barrel", 400, -8, 8);
    ResLayer2_barrel = dbe_->book1D("ResLayer2_barrel", "ResLayer2_barrel", 400, -8, 8);
    ResLayer3_barrel = dbe_->book1D("ResLayer3_barrel", "ResLayer3_barrel", 400, -8, 8);
    ResLayer4_barrel = dbe_->book1D("ResLayer4_barrel", "ResLayer4_barrel", 400, -8, 8);
    ResLayer5_barrel = dbe_->book1D("ResLayer5_barrel", "ResLayer5_barrel", 400, -8, 8);
    ResLayer6_barrel = dbe_->book1D("ResLayer6_barrel", "ResLayer6_barrel", 400, -8, 8);

    BxDist = dbe_->book1D("Bunch_Crossing", "Bunch_Crossing", 20, -10., 10.);
    StripProf = dbe_->book1D("Strip_Profile", "Strip_Profile", 100, 0, 100);

    BxDisc_4Plus = dbe_->book1D("BxDisc_4Plus", "BxDisc_4Plus", 20, -10., 10.);
    BxDisc_4Min = dbe_->book1D("BxDisc_4Min", "BxDisc_4Min", 20, -10., 10.);

    //cls histos
    noiseCLS = dbe_->book1D("noiseCLS", "noiseCLS", 10, 0.5, 10.5);
    noiseCLSBarrel = dbe_->book1D("noiseCLSBarrel", "noiseCLSBarrel", 10, 0.5, 10.5);
    noiseCLSEndcaps = dbe_->book1D("noiseCLSEndcaps", "noiseCLSEndcaps", 10, 0.5, 10.5);

    clsBarrel = dbe_->book1D("clsBarrel", "clsBarrel", 10, 0.5, 10.5);

    //endcap CLS
    CLS_Endcap_1_Ring2_A = dbe_->book1D("CLS_Endcap_1_1Ring2_A", "CLS_Endcap_1_Ring2_A", 10, 0.5, 10.5);
    CLS_Endcap_1_Ring2_B = dbe_->book1D("CLS_Endcap_1_1Ring2_B", "CLS_Endcap_1_Ring2_B", 10, 0.5, 10.5);
    CLS_Endcap_1_Ring2_C = dbe_->book1D("CLS_Endcap_1_1Ring2_C", "CLS_Endcap_1_Ring2_C", 10, 0.5, 10.5);

    CLS_Endcap_23_Ring2_A = dbe_->book1D("CLS_Endcap_23_Ring2_A", "CLS_Endcap_23_Ring2_A", 10, 0.5, 10.5);
    CLS_Endcap_23_Ring2_B = dbe_->book1D("CLS_Endcap_23_Ring2_B", "CLS_Endcap_23_Ring2_B", 10, 0.5, 10.5);
    CLS_Endcap_23_Ring2_C = dbe_->book1D("CLS_Endcap_23_Ring2_C", "CLS_Endcap_23_Ring2_C", 10, 0.5, 10.5);

    CLS_Endcap_123_Ring3_A = dbe_->book1D("CLS_Endcap_123_Ring3_A", "CLS_Endcap_123_Ring3_A", 10, 0.5, 10.5);
    CLS_Endcap_123_Ring3_B = dbe_->book1D("CLS_Endcap_123_Ring3_B", "CLS_Endcap_123_Ring3_B", 10, 0.5, 10.5);
    CLS_Endcap_123_Ring3_C = dbe_->book1D("CLS_Endcap_123_Ring3_C", "CLS_Endcap_123_Ring3_C", 10, 0.5, 10.5);

    CLS_Endcap_4 = dbe_->book1D("CLS_Endcap_4", "CLS_Endcap_4", 10, 0.5, 10.5);

    //endcap residuals
    ResDmin1 = dbe_->book1D("Disk_Min1_Residuals", "Disk_Min1_Residuals", 400, -8, 8);
    ResDmin2 = dbe_->book1D("Disk_Min2_Residuals", "Disk_Min2_Residuals", 400, -8, 8);
    ResDmin3 = dbe_->book1D("Disk_Min3_Residuals", "Disk_Min3_Residuals", 400, -8, 8);
    ResDplu1 = dbe_->book1D("Disk_Plu1_Residuals", "Disk_Plu1_Residuals", 400, -8, 8);
    ResDplu2 = dbe_->book1D("Disk_Plu2_Residuals", "Disk_Plu2_Residuals", 400, -8, 8);
    ResDplu3 = dbe_->book1D("Disk_Plu3_Residuals", "Disk_Plu3_Residuals", 400, -8, 8);

    ResDmin4 = dbe_->book1D("Disk_Min4_Residuals", "Disk_Min4_Residuals", 400, -8, 8);
    ResDplu4 = dbe_->book1D("Disk_Plu4_Residuals", "Disk_Plu4_Residuals", 400, -8, 8);

    Res_Endcap1_Ring2_A = dbe_->book1D("Res_Endcap1_Ring2_A", "Res_Endcap1_Ring2_A", 400, -8, 8);
    Res_Endcap1_Ring2_B = dbe_->book1D("Res_Endcap1_Ring2_B", "Res_Endcap1_Ring2_B", 400, -8, 8);
    Res_Endcap1_Ring2_C = dbe_->book1D("Res_Endcap1_Ring2_C", "Res_Endcap1_Ring2_C", 400, -8, 8);

    Res_Endcap23_Ring2_A = dbe_->book1D("Res_Endcap23_Ring2_A", "Res_Endcap23_Ring2_A", 400, -8, 8);
    Res_Endcap23_Ring2_B = dbe_->book1D("Res_Endcap23_Ring2_B", "Res_Endcap23_Ring2_B", 400, -8, 8);
    Res_Endcap23_Ring2_C = dbe_->book1D("Res_Endcap23_Ring2_C", "Res_Endcap23_Ring2_C", 400, -8, 8);

    Res_Endcap123_Ring3_A = dbe_->book1D("Res_Endcap123_Ring3_A", "Res_Endcap123_Ring3_A", 400, -8, 8);
    Res_Endcap123_Ring3_B = dbe_->book1D("Res_Endcap123_Ring3_B", "Res_Endcap123_Ring3_B", 400, -8, 8);
    Res_Endcap123_Ring3_C = dbe_->book1D("Res_Endcap123_Ring3_C", "Res_Endcap123_Ring3_C", 400, -8, 8);

  }//end dbe
}

void RPCDigiValid::endJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 354 of file RPCDigiValid.cc.

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

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

}

void RPCDigiValid::endRun ( edm::Run const &  run, edm::EventSetup const &  eSetup  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 443 of file RPCDigiValid.cc.

{
}

Member Data Documentation

MonitorElement* RPCDigiValid::BxDisc_4Min

private

Definition at line 118 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::BxDisc_4Plus

private

Definition at line 117 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::BxDist

private

Definition at line 57 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_123_Ring3_A

private

Definition at line 108 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_123_Ring3_B

private

Definition at line 109 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_123_Ring3_C

private

Definition at line 110 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_1_Ring2_A

private

Definition at line 98 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_1_Ring2_B

private

Definition at line 99 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_1_Ring2_C

private

Definition at line 100 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_23_Ring2_A

private

Definition at line 103 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_23_Ring2_B

private

Definition at line 104 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_23_Ring2_C

private

Definition at line 105 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::CLS_Endcap_4

private

Definition at line 121 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::clsBarrel

private

Definition at line 94 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

int RPCDigiValid::countEvent

private

Definition at line 132 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

DQMStore* RPCDigiValid::dbe_

private

Definition at line 134 of file RPCDigiValid.h.

Referenced by beginRun(), endJob(), and RPCDigiValid().

std::string RPCDigiValid::digiLabel

private

Definition at line 136 of file RPCDigiValid.h.

std::map<RPCDetId, double> RPCDigiValid::mapRollArea

private

Definition at line 125 of file RPCDigiValid.h.

Referenced by analyze().

std::map<RPCDetId, double> RPCDigiValid::mapRollCls

private

Definition at line 124 of file RPCDigiValid.h.

Referenced by analyze().

std::map<RPCDetId, int> RPCDigiValid::mapRollFakeCount

private

Definition at line 127 of file RPCDigiValid.h.

Referenced by analyze().

std::map<RPCDetId, std::string> RPCDigiValid::mapRollName

private

Definition at line 129 of file RPCDigiValid.h.

Referenced by analyze().

std::map<RPCDetId, std::map<int, double>*> RPCDigiValid::mapRollNoisyStripRate

private

Definition at line 131 of file RPCDigiValid.h.

Referenced by analyze().

std::map<RPCDetId, double> RPCDigiValid::mapRollStripArea

private

Definition at line 126 of file RPCDigiValid.h.

Referenced by analyze().

std::map<RPCDetId, std::map<int, double>*> RPCDigiValid::mapRollStripRate

private

Definition at line 130 of file RPCDigiValid.h.

Referenced by analyze().

std::map<RPCDetId, int> RPCDigiValid::mapRollTruCount

private

Definition at line 128 of file RPCDigiValid.h.

Referenced by analyze().

MonitorElement* RPCDigiValid::noiseCLS

private

Definition at line 90 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::noiseCLSBarrel

private

Definition at line 91 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::noiseCLSEndcaps

private

Definition at line 92 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

std::string RPCDigiValid::outputFile_

private

Definition at line 135 of file RPCDigiValid.h.

Referenced by endJob(), and RPCDigiValid().

MonitorElement* RPCDigiValid::Res

private

Definition at line 51 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::Res_Endcap123_Ring3_A

private

Definition at line 85 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::Res_Endcap123_Ring3_B

private

Definition at line 86 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::Res_Endcap123_Ring3_C

private

Definition at line 87 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::Res_Endcap1_Ring2_A

private

Definition at line 77 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::Res_Endcap1_Ring2_B

private

Definition at line 78 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::Res_Endcap1_Ring2_C

private

Definition at line 79 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::Res_Endcap23_Ring2_A

private

Definition at line 81 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::Res_Endcap23_Ring2_B

private

Definition at line 82 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::Res_Endcap23_Ring2_C

private

Definition at line 83 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResDmin1

private

Definition at line 69 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResDmin2

private

Definition at line 70 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResDmin3

private

Definition at line 71 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResDmin4

private

Definition at line 115 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResDplu1

private

Definition at line 72 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResDplu2

private

Definition at line 73 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResDplu3

private

Definition at line 74 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResDplu4

private

Definition at line 116 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResLayer1_barrel

private

Definition at line 61 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResLayer2_barrel

private

Definition at line 62 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResLayer3_barrel

private

Definition at line 63 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResLayer4_barrel

private

Definition at line 64 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResLayer5_barrel

private

Definition at line 65 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResLayer6_barrel

private

Definition at line 66 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResWmin1

private

Definition at line 53 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResWmin2

private

Definition at line 52 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResWplu1

private

Definition at line 55 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResWplu2

private

Definition at line 56 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::ResWzer0

private

Definition at line 54 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

edm::EDGetTokenT<RPCDigiCollection> RPCDigiValid::rpcDigiToken

private

Definition at line 140 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::rzview

private

Definition at line 50 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

edm::EDGetTokenT<edm::PSimHitContainer> RPCDigiValid::simHitToken

private

Definition at line 139 of file RPCDigiValid.h.

Referenced by analyze(), and RPCDigiValid().

MonitorElement* RPCDigiValid::StripProf

private

Definition at line 58 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::xyvDmin4

private

Definition at line 120 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::xyvDplu4

private

Definition at line 119 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().

MonitorElement* RPCDigiValid::xyview

private

Definition at line 49 of file RPCDigiValid.h.

Referenced by analyze(), and beginRun().