CTPPSPixelDQMSource Class Reference

Inheritance diagram for CTPPSPixelDQMSource:

Public Member Functions
	CTPPSPixelDQMSource (const edm::ParameterSet &ps)
virtual	~CTPPSPixelDQMSource ()
Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final
virtual void	beginStream (edm::StreamID id) final
	DQMEDAnalyzer (void)
virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
uint32_t	streamId () const
Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default
Public Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
virtual	~EDAnalyzerBase ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	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 &	itemsToGetFromEvent () 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
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup)
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &eSetup)
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &eSetup)
void	endRun (edm::Run const &run, edm::EventSetup const &eSetup)
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 Member Functions
int	getDet (int id)
int	getPixPlane (int id)
int	getPlaneIndex (int arm, int station, int rp, int plane)
int	getRPindex (int arm, int station, int rp)
int	prIndex (int rp, int plane)

Private Attributes
const int	CluSizeMAX = 25
int	cluSizeMaxData
const int	ClusMultMAX = 90
MonitorElement *	h2CluSize [NArms][NStationMAX]
MonitorElement *	h2ClusMultipl [NArms][NStationMAX]
MonitorElement *	h2HitsMultipl [NArms][NStationMAX]
MonitorElement *	h2HitsMultROC [RPotsTotalNumber]
MonitorElement *	h2PlaneActive [NArms][NStationMAX]
MonitorElement *	h2xyHits [RPotsTotalNumber][NplaneMAX]
MonitorElement *	h2xyROCadc [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	h2xyROCHits [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	hBX
MonitorElement *	hBXshort
MonitorElement *	hHitsMult [RPotsTotalNumber][NplaneMAX]
const int	hitMultMAX = 300
int	HitsMultPlane [RPotsTotalNumber][NplaneMAX]
int	HitsMultROC [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	hp2HitsOcc [NArms][NStationMAX]
MonitorElement *	hp2xyADC [RPotsTotalNumber][NplaneMAX]
MonitorElement *	hRPotActivBX [RPotsTotalNumber]
MonitorElement *	hRPotActivBXroc [RPotsTotalNumber]
MonitorElement *	hRPotActivPlanes [RPotsTotalNumber]
MonitorElement *	hRPotActivROCs [RPotsTotalNumber]
MonitorElement *	hRPotActivROCsMax [RPotsTotalNumber]
const int	IndexNotValid = 0
int	multClus
int	multHits
long int	nEvents = 0
int	pixColMAX = 156
int	pixRowMAX = 160
int	ROCSizeInX = pixRowMAX/2
int	ROCSizeInY = pixColMAX/3
int	RPindexValid [RPotsTotalNumber]
const int	RPn_first = 3
const int	RPn_last = 4
int	RPstatus [NStationMAX][NRPotsMAX]
unsigned int	rpStatusWord = 0x8000
int	StationStatus [NStationMAX]
CTPPSPixelIndices	thePixIndices
edm::EDGetTokenT< edm::DetSetVector< CTPPSPixelDigi > >	tokenDigi
unsigned int	verbosity

Static Private Attributes
static int	NArms =2
static int	NplaneMAX =6
static int	NROCsMAX = 6
static int	NRPotsMAX =6
static int	NStationMAX =3
static int	RPotsTotalNumber =NArms*NStationMAX*NRPotsMAX

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
static std::shared_ptr< dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 31 of file CTPPSPixelDQMSource.cc.

Constructor & Destructor Documentation

CTPPSPixelDQMSource::CTPPSPixelDQMSource

(

const edm::ParameterSet &

ps

)

Definition at line 132 of file CTPPSPixelDQMSource.cc.

References edm::ParameterSet::getParameter(), and tokenDigi.

                                                                  :
  verbosity(ps.getUntrackedParameter<unsigned int>("verbosity", 0)),
  rpStatusWord(ps.getUntrackedParameter<unsigned int>("RPStatusWord",0x8000))
{
 tokenDigi = consumes<DetSetVector<CTPPSPixelDigi> >(ps.getParameter<edm::InputTag>("tagRPixDigi"));
// tokenCluster=consumes<DetSetVector<CTPPSPixelCluster>>(ps.getParameter<edm::InputTag>("tagRPixCluster"));

}

CTPPSPixelDQMSource::~CTPPSPixelDQMSource ( )

virtual

Definition at line 143 of file CTPPSPixelDQMSource.cc.

{
}

Member Function Documentation

void CTPPSPixelDQMSource::analyze ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

protectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 334 of file CTPPSPixelDQMSource.cc.

References ecalMGPA::adc(), CTPPSDetId::arm(), edm::DetSet< T >::begin(), edm::EventBase::bunchCrossing(), cuy::col, MonitorElement::Fill(), HcalObjRepresent::Fill(), getDet(), getPixPlane(), getPlaneIndex(), getRPindex(), h2HitsMultipl, h2HitsMultROC, h2PlaneActive, h2xyHits, h2xyROCadc, h2xyROCHits, hBX, hBXshort, hHitsMult, HitsMultPlane, HitsMultROC, hp2HitsOcc, hp2xyADC, hRPotActivBX, hRPotActivBXroc, hRPotActivPlanes, hRPotActivROCs, hRPotActivROCsMax, diffTreeTool::index, edm::HandleBase::isValid(), hpstanc_transforms::max, multHits, nEvents, np, NplaneMAX, NROCsMAX, NRPotsMAX, NStationMAX, AlCaHLTBitMon_ParallelJobs::p, prIndex(), alignCSCRings::r, CTPPSDetId::rp(), RPindexValid, RPstatus, CTPPSDetId::station(), relativeConstraints::station, StationStatus, thePixIndices, tokenDigi, CTPPSPixelIndices::transformToROC(), verbosity, and DetId::VeryForward.

{
  ++nEvents;
 int RPactivity[2][NRPotsMAX];
 for(int arm = 0; arm <2; arm++) for(int rp=0; rp<NRPotsMAX; rp++)
   RPactivity[arm][rp] = 0;
 for(int ind=0; ind<2*3*NRPotsMAX; ind++) 
   for(int p=0; p<NplaneMAX; p++) HitsMultPlane[ind][p] = 0;

 for(int ind=0; ind<2*3*NRPotsMAX*NplaneMAX; ind++) 
   for(int rp=0; rp<NROCsMAX; rp++) HitsMultROC[ind][rp] = 0;

  Handle< DetSetVector<CTPPSPixelDigi> > pixDigi;
  event.getByToken(tokenDigi, pixDigi);

  hBX->Fill(event.bunchCrossing());
  hBXshort->Fill(event.bunchCrossing());

  bool valid = false;
  valid |= pixDigi.isValid();
//  valid |= Clus.isValid();

 if(!valid && verbosity) LogPrint("CTPPSPixelDQMSource") <<"No valid data in Event "<<nEvents;

if(pixDigi.isValid())
  for(const auto &ds_digi : *pixDigi)
  {
   int idet = getDet(ds_digi.id);
   if(idet != DetId::VeryForward) {
     if(verbosity>1) LogPrint("CTPPSPixelDQMSource") <<"not CTPPS: ds_digi.id"<<ds_digi.id;
     continue;
   }
//   int subdet = getSubdet(ds_digi.id);

   int plane = getPixPlane(ds_digi.id);
   
   CTPPSDetId theId(ds_digi.id);
   int arm = theId.arm()&0x1;
   int station = theId.station()&0x3;
   int rpot = theId.rp()&0x7;
   RPactivity[arm][rpot] = 1;

 if(StationStatus[station] && RPstatus[station][rpot]) {

   hp2HitsOcc[arm][station]->Fill(plane,rpot,(int)ds_digi.data.size());
   h2HitsMultipl[arm][station]->Fill(prIndex(rpot,plane),ds_digi.data.size());
   h2PlaneActive[arm][station]->Fill(plane,rpot);

    int index = getRPindex(arm,station,rpot);
    HitsMultPlane[index][plane] += ds_digi.data.size();
    if(RPindexValid[index])
      hHitsMult[index][plane]->Fill(ds_digi.data.size());
    int rocHistIndex = getPlaneIndex(arm,station,rpot,plane);

    for(DetSet<CTPPSPixelDigi>::const_iterator dit = ds_digi.begin();
     dit != ds_digi.end(); ++dit) {
      int row = dit->row();
      int col = dit->column();
      int adc = dit->adc();

      if(RPindexValid[index]) {
        h2xyHits[index][plane]->Fill(col,row);
        hp2xyADC[index][plane]->Fill(col,row,adc);
        int colROC, rowROC;
        int trocId;
        if(!thePixIndices.transformToROC(col,row, trocId, colROC, rowROC)) {
     if(trocId>=0 && trocId<NROCsMAX) {
       h2xyROCHits[rocHistIndex][trocId]->Fill(rowROC,colROC);
       h2xyROCadc[rocHistIndex][trocId]->Fill(rowROC,colROC,adc);
       ++HitsMultROC[rocHistIndex][trocId];
     }
        }
      } //end if(RPindexValid[index]) {
     }

   if(int(ds_digi.data.size()) > multHits) multHits = ds_digi.data.size();
  } // end  if(StationStatus[station]) {
 } // end for(const auto &ds_digi : *pixDigi)

 for(int arm=0; arm<2; arm++) for(int stn=0; stn<NStationMAX; stn++) 
  for(int rp=0; rp<NRPotsMAX; rp++) {
    int index = getRPindex(arm,stn,rp);
    if(RPindexValid[index]==0) continue;
    if(RPactivity[arm][rp]==0) continue;

    int np = 0; 
    for(int p=0; p<NplaneMAX; p++) if(HitsMultPlane[index][p]>0) np++;
    hRPotActivPlanes[index]->Fill(np);
    if(np>5) hRPotActivBX[index]->Fill(event.bunchCrossing());

    int rocf[NplaneMAX];
    for(int r=0; r<NROCsMAX; r++) rocf[r]=0;
    for(int p=0; p<NplaneMAX; p++) {
      int indp = getPlaneIndex(arm,stn,rp,p);
      for(int r=0; r<NROCsMAX; r++) if(HitsMultROC[indp][r] > 0) ++rocf[r];
      for(int r=0; r<NROCsMAX; r++) h2HitsMultROC[index]->Fill(p,r,HitsMultROC[indp][r]);
    }
    int max = 0;
    for(int r=0; r<NROCsMAX; r++) if(max < rocf[r]) 
    { max = rocf[r]; }

    for(int r=0; r<NROCsMAX; r++) hRPotActivROCs[index]->Fill(r,rocf[r]);
    for(int r=0; r<NROCsMAX; r++) if(rocf[r] == max)
    hRPotActivROCsMax[index]->Fill(r,max);

    if(max > 4) hRPotActivBXroc[index]->Fill(event.bunchCrossing());
  }

  if((nEvents % 100)) return;
  if(verbosity) LogPrint("CTPPSPixelDQMSource")<<"analyze event "<<nEvents;
}

void CTPPSPixelDQMSource::beginLuminosityBlock ( edm::LuminosityBlock const &  lumi, edm::EventSetup const &  eSetup  )

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 328 of file CTPPSPixelDQMSource.cc.

{
}

void CTPPSPixelDQMSource::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 176 of file CTPPSPixelDQMSource.cc.

References CTPPSDetId::armName(), DQMStore::IBooker::book1D(), DQMStore::IBooker::book1DD(), DQMStore::IBooker::book2D(), DQMStore::IBooker::book2DD(), DQMStore::IBooker::bookProfile2D(), DQMStore::IBooker::cd(), CluSizeMAX, ClusMultMAX, MillePedeFileConverter_cfg::e, getPlaneIndex(), CTPPSDetId::getRPId(), getRPindex(), CTPPSDetId::getStationId(), MonitorElement::getTH2D(), MonitorElement::getTH2F(), MonitorElement::getTProfile2D(), h2CluSize, h2ClusMultipl, h2HitsMultipl, h2HitsMultROC, h2PlaneActive, h2xyHits, h2xyROCadc, h2xyROCHits, hBX, hBXshort, hHitsMult, hitMultMAX, hp2HitsOcc, hp2xyADC, hRPotActivBX, hRPotActivBXroc, hRPotActivPlanes, hRPotActivROCs, hRPotActivROCsMax, diffTreeTool::index, pileupCalc::nbins, CTPPSDetId::nFull, CTPPSDetId::nPath, NplaneMAX, NROCsMAX, CTPPSDetId::nShort, NStationMAX, AlCaHLTBitMon_ParallelJobs::p, pixRowMAX, ROCSizeInX, ROCSizeInY, RPindexValid, RPn_first, RPn_last, RPstatus, alignCSCRings::s, sd, CTPPSDetId::sdTrackingPixel, DQMStore::IBooker::setCurrentFolder(), CTPPSDetId::setRP(), CTPPSDetId::setStation(), StationStatus, AlCaHLTBitMon_QueryRunRegistry::string, and TrackerOfflineValidation_Dqm_cff::xmax.

{
  ibooker.cd();
  ibooker.setCurrentFolder("CTPPS/TrackingPixel");
  char s[50];
  hBX = ibooker.book1D("events per BX", "ctpps_pixel;Event.BX", 4002, -1.5, 4000. +0.5);
  hBXshort = ibooker.book1D("events per BX(short)", "ctpps_pixel;Event.BX", 102, -1.5, 100. + 0.5);

 for(int arm=0; arm<2; arm++) {
   CTPPSDetId ID(CTPPSDetId::sdTrackingPixel,arm,0);
   string sd, armTitle;
   ID.armName(sd, CTPPSDetId::nPath);
   ID.armName(armTitle, CTPPSDetId::nFull);

   ibooker.setCurrentFolder(sd);

 for(int stn=2; stn<NStationMAX; stn++) {
   if(StationStatus[stn]==0) continue;
   ID.setStation(stn);
   string stnd, stnTitle;

   CTPPSDetId(ID.getStationId()).stationName(stnd, CTPPSDetId::nPath);
   CTPPSDetId(ID.getStationId()).stationName(stnTitle, CTPPSDetId::nFull);

   ibooker.setCurrentFolder(stnd);

   string st = "planes activity";
   string st2 = ": " + stnTitle;

   int rpnbins = RPn_last-RPn_first; 
  
   h2PlaneActive[arm][stn] = ibooker.book2DD(st,st+st2+";Plane #",
        NplaneMAX,0,NplaneMAX, rpnbins, RPn_first,RPn_last);
   TH2D *h = h2PlaneActive[arm][stn]->getTH2D();
   h->SetOption("colz");
   TAxis *yah = h->GetYaxis();
   
   st = "hit average multiplicity in planes";

   hp2HitsOcc[arm][stn]= ibooker.bookProfile2D(st,st+st2+";Plane #",
     NplaneMAX, 0, NplaneMAX, rpnbins, RPn_first,RPn_last,0,20000);
   TProfile2D *h2 = hp2HitsOcc[arm][stn]->getTProfile2D();
   h2->SetOption("textcolz");
   TAxis *yah2 = h2->GetYaxis();

   int xmax = NplaneMAX*rpnbins;

   st = "hit multiplicity in planes";
   string st3 = ";PlaneIndex(=pixelPot*PlaneMAX + plane)";
   h2HitsMultipl[arm][stn]= ibooker.book2DD(st,st+st2+st3+";multiplicity",
    xmax,0,xmax,hitMultMAX,0,hitMultMAX);
   h2HitsMultipl[arm][stn]->getTH2D()->SetOption("colz");

   st = "cluster multiplicity in planes";
   h2ClusMultipl[arm][stn] = ibooker.book2DD(st,st+st2+st3+";multiplicity",
    xmax,0,xmax, ClusMultMAX,0,ClusMultMAX);
   h2ClusMultipl[arm][stn]->getTH2D()->SetOption("colz");

   st = "cluster size in planes";
  h2CluSize[arm][stn] = ibooker.book2D(st,st+st2+st3+";Cluster size",
    xmax,0,xmax, CluSizeMAX,0,CluSizeMAX);
   h2CluSize[arm][stn]->getTH2F()->SetOption("colz");

//--------- Hits ---
   int pixBinW = 4;
     for(int rp=RPn_first; rp<RPn_last; rp++) { // only installed pixel pots
       ID.setRP(rp);
       string rpd, rpTitle;
       CTPPSDetId(ID.getRPId()).rpName(rpTitle, CTPPSDetId::nShort);
    yah->SetBinLabel(rp - RPn_first +1, rpTitle.c_str()); // h
       yah2->SetBinLabel(rp - RPn_first +1, rpTitle.c_str()); //h2

       if(RPstatus[stn][rp]==0) continue;
       int indexP = getRPindex(arm,stn,rp);
       RPindexValid[indexP] = 1;

       CTPPSDetId(ID.getRPId()).rpName(rpTitle, CTPPSDetId::nFull);
       CTPPSDetId(ID.getRPId()).rpName(rpd, CTPPSDetId::nPath);

       ibooker.setCurrentFolder(rpd);

       hRPotActivPlanes[indexP] = 
        ibooker.book1D("number of fired planes per event", rpTitle+";nPlanes",
     NplaneMAX, -0.5, NplaneMAX+0.5);
       hRPotActivBX[indexP] = 
        ibooker.book1D("5 fired planes per BX", rpTitle+";Event.BX", 4002, -1.5, 4000.+0.5);
       hRPotActivBXroc[indexP] = 
        ibooker.book1D("4 fired ROCs per BX", rpTitle+";Event.BX", 4002, -1.5, 4000.+0.5);

       h2HitsMultROC[indexP] = ibooker.bookProfile2D("ROCs hits multiplicity per event",
       rpTitle+";plane # ;ROC #", NplaneMAX,-0.5,NplaneMAX-0.5, NROCsMAX,-0.5,NROCsMAX-0.5, 100,0,1.e3);

       h2HitsMultROC[indexP]->getTProfile2D()->SetOption("colztext");
       h2HitsMultROC[indexP]->getTProfile2D()->SetMinimum(1.e-10);

       hRPotActivROCs[indexP] = ibooker.book2D("number of fired aligned_ROCs per event",
     rpTitle+";ROC ID;number of fired ROCs", NROCsMAX,-0.5,NROCsMAX-0.5, 7,-0.5,6.5);
       hRPotActivROCs[indexP]->getTH2F()->SetOption("colz");

       hRPotActivROCsMax[indexP]= ibooker.book2D("max number of fired aligned_ROCs per event",
     rpTitle+";ROC ID;number of fired ROCs", NROCsMAX,-0.5,NROCsMAX-0.5, 7,-0.5,6.5);
       hRPotActivROCsMax[indexP]->getTH2F()->SetOption("colz");

       int nbins = pixRowMAX/pixBinW;

       for(int p=0; p<NplaneMAX; p++) {
         sprintf(s,"plane_%d",p);
         string pd = rpd+"/"+string(s);
         ibooker.setCurrentFolder(pd);
         string st1 = ": "+rpTitle+"_"+string(s);
         st = "hits position";
         h2xyHits[indexP][p]=ibooker.book2DD(st,st1+";pix col;pix row", 
       nbins,0,pixRowMAX,nbins,0,pixRowMAX);
     h2xyHits[indexP][p]->getTH2D()->SetOption("colz");

     st = "adc average value";
     hp2xyADC[indexP][p]=ibooker.bookProfile2D(st,st1+";pix col;pix row",
       nbins,0,pixRowMAX,nbins,0,pixRowMAX,100,0,1.e10,"");
     hp2xyADC[indexP][p]->getTProfile2D()->SetOption("colz");

         st = "hits multiplicity";
         hHitsMult[indexP][p]=ibooker.book1DD(st,st1+";number of hits;N / 1 hit",
       hitMultMAX,0,hitMultMAX);

         ibooker.setCurrentFolder(pd + "/ROCs");
         int index = getPlaneIndex(arm,stn,rp,p);

         for(int roc=0; roc<6; roc++) {
       sprintf(s,"ROC_%d",roc);
       string st2 = st1 + "_" + string(s);
           ibooker.setCurrentFolder(pd + "/ROCs/" + string(s));

           h2xyROCHits[index][roc]=ibooker.book2DD("hits",st2+";pix row;pix col",
        ROCSizeInX,0,ROCSizeInX,ROCSizeInY,0,ROCSizeInY);
           h2xyROCHits[index][roc]->getTH2D()->SetOption("colz");

           string st = "adc average value";
           h2xyROCadc[index][roc]=ibooker.bookProfile2D(st,st2+";pix row;pix col",
        ROCSizeInX,0,ROCSizeInX,ROCSizeInY,0,ROCSizeInY, 0,512);
           h2xyROCadc[index][roc]->getTProfile2D()->SetOption("colz");
         }
       } // end of for(int p=0; p<NplaneMAX;..

     } // end for(int rp=0; rp<NRPotsMAX;...
   } // end of for(int stn=0; stn<
  } //end of for(int arm=0; arm<2;...

 return;
}

void CTPPSPixelDQMSource::dqmBeginRun ( edm::Run const &  run, edm::EventSetup const &    )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 149 of file CTPPSPixelDQMSource.cc.

References cluSizeMaxData, CTPPSPixelIndices::getDefaultColDetSize(), CTPPSPixelIndices::getDefaultRowDetSize(), multClus, multHits, nEvents, NRPotsMAX, pixColMAX, pixRowMAX, ROCSizeInX, ROCSizeInY, RPindexValid, RPstatus, rpStatusWord, StationStatus, thePixIndices, and verbosity.

{
  if(verbosity) LogPrint("CTPPSPixelDQMSource") <<"RPstatusWord= "<<rpStatusWord;
  nEvents = 0;

  pixRowMAX = thePixIndices.getDefaultRowDetSize();
  pixColMAX = thePixIndices.getDefaultColDetSize();
  ROCSizeInX = pixRowMAX/2;  // ROC row size in pixels = 80
  ROCSizeInY = pixColMAX/3;

 unsigned int rpSts = rpStatusWord<<1;
 for(int stn=0; stn<3; stn++) {
   int stns = 0;
   for(int rp=0; rp<NRPotsMAX; rp++) {
     rpSts = (rpSts >> 1); RPstatus[stn][rp] = rpSts&1;
     if(RPstatus[stn][rp] > 0) stns = 1;
   }
   StationStatus[stn]=stns;
 }
 
 for(int ind=0; ind<2*3*NRPotsMAX; ind++) RPindexValid[ind] = 0;

 multHits = multClus = cluSizeMaxData = -1;
}

void CTPPSPixelDQMSource::endLuminosityBlock ( edm::LuminosityBlock const &  lumi, edm::EventSetup const &  eSetup  )

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 447 of file CTPPSPixelDQMSource.cc.

{
}

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

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 452 of file CTPPSPixelDQMSource.cc.

References cluSizeMaxData, DEFINE_FWK_MODULE, multClus, multHits, nEvents, edm::RunBase::run(), and verbosity.

{
 if(!verbosity) return; 
 LogPrint("CTPPSPixelDQMSource") 
  <<"end of Run "<<run.run()<<": "<<nEvents<<" events\n"
  <<"mult Hits/Clus: "<<multHits<<" / "<<multClus
  <<"   cluSizeMaxData= "<<cluSizeMaxData;
}

int CTPPSPixelDQMSource::getDet ( int  id )

inlineprivate

Definition at line 115 of file CTPPSPixelDQMSource.cc.

References DetId::kDetOffset.

Referenced by analyze().

    { return (id>>DetId::kDetOffset)&0xF; }

int CTPPSPixelDQMSource::getPixPlane ( int  id )

inlineprivate

Definition at line 117 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

    { return ((id>>16)&0x7); }

int CTPPSPixelDQMSource::getPlaneIndex ( int  arm, int  station, int  rp, int  plane  )

inlineprivate

Definition at line 106 of file CTPPSPixelDQMSource.cc.

References getRPindex().

Referenced by analyze(), and bookHistograms().

                                                             {
    if(plane<0 || plane>=NplaneMAX) return(IndexNotValid);
    int rc = getRPindex(arm, station, rp);
    if(rc == IndexNotValid) return(IndexNotValid);
    return(rc*NplaneMAX + plane);
  }

int CTPPSPixelDQMSource::getRPindex ( int  arm, int  station, int  rp  )

inlineprivate

Definition at line 99 of file CTPPSPixelDQMSource.cc.

References relativeConstraints::station.

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

                                               {
    if(arm<0 || station<0 || rp<0) return(IndexNotValid);
    if(arm>1 || station>=NStationMAX || rp>=NRPotsMAX) return(IndexNotValid);
    int rc = (arm*NStationMAX+station)*NRPotsMAX + rp;
    return(rc);
  }

int CTPPSPixelDQMSource::prIndex ( int  rp, int  plane  )

inlineprivate

Definition at line 113 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

    {return((rp - RPn_first)*NplaneMAX + plane);}

Member Data Documentation

const int CTPPSPixelDQMSource::CluSizeMAX = 25

private

Definition at line 59 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::cluSizeMaxData

private

Definition at line 121 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun(), and endRun().

const int CTPPSPixelDQMSource::ClusMultMAX = 90

private

Definition at line 58 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2CluSize[NArms][NStationMAX]

private

Definition at line 74 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2ClusMultipl[NArms][NStationMAX]

private

Definition at line 73 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2HitsMultipl[NArms][NStationMAX]

private

Definition at line 71 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2HitsMultROC[RPotsTotalNumber]

private

Definition at line 82 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2PlaneActive[NArms][NStationMAX]

private

Definition at line 72 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2xyHits[RPotsTotalNumber][NplaneMAX]

private

Definition at line 86 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2xyROCadc[RPotsTotalNumber *NplaneMAX][NROCsMAX]

private

Definition at line 89 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2xyROCHits[RPotsTotalNumber *NplaneMAX][NROCsMAX]

private

Definition at line 88 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hBX

private

Definition at line 68 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement * CTPPSPixelDQMSource::hBXshort

private

Definition at line 68 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hHitsMult[RPotsTotalNumber][NplaneMAX]

private

Definition at line 85 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

const int CTPPSPixelDQMSource::hitMultMAX = 300

private

Definition at line 57 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::HitsMultPlane[RPotsTotalNumber][NplaneMAX]

private

Definition at line 91 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

int CTPPSPixelDQMSource::HitsMultROC[RPotsTotalNumber *NplaneMAX][NROCsMAX]

private

Definition at line 90 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

MonitorElement* CTPPSPixelDQMSource::hp2HitsOcc[NArms][NStationMAX]

private

Definition at line 70 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hp2xyADC[RPotsTotalNumber][NplaneMAX]

private

Definition at line 87 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBX[RPotsTotalNumber]

private

Definition at line 80 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBXroc[RPotsTotalNumber]

private

Definition at line 81 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivPlanes[RPotsTotalNumber]

private

Definition at line 79 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivROCs[RPotsTotalNumber]

private

Definition at line 83 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivROCsMax[RPotsTotalNumber]

private

Definition at line 84 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

const int CTPPSPixelDQMSource::IndexNotValid = 0

private

Definition at line 97 of file CTPPSPixelDQMSource.cc.

int CTPPSPixelDQMSource::multClus

private

Definition at line 121 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun(), and endRun().

int CTPPSPixelDQMSource::multHits

private

Definition at line 121 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), dqmBeginRun(), and endRun().

int CTPPSPixelDQMSource::NArms =2

staticprivate

Definition at line 51 of file CTPPSPixelDQMSource.cc.

long int CTPPSPixelDQMSource::nEvents = 0

private

Definition at line 66 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), dqmBeginRun(), endRun(), and looper.Looper::loop().

int CTPPSPixelDQMSource::NplaneMAX =6

staticprivate

Definition at line 54 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::NROCsMAX = 6

staticprivate

Definition at line 55 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::NRPotsMAX =6

staticprivate

Definition at line 53 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

int CTPPSPixelDQMSource::NStationMAX =3

staticprivate

Definition at line 52 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::pixColMAX = 156

private

Definition at line 63 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::pixRowMAX = 160

private

Definition at line 62 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::ROCSizeInX = pixRowMAX/2

private

Definition at line 64 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::ROCSizeInY = pixColMAX/3

private

Definition at line 65 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::RPindexValid[RPotsTotalNumber]

private

Definition at line 78 of file CTPPSPixelDQMSource.cc.

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

const int CTPPSPixelDQMSource::RPn_first = 3

private

Definition at line 56 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

const int CTPPSPixelDQMSource::RPn_last = 4

private

Definition at line 56 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::RPotsTotalNumber =NArms*NStationMAX*NRPotsMAX

staticprivate

Definition at line 76 of file CTPPSPixelDQMSource.cc.

int CTPPSPixelDQMSource::RPstatus[NStationMAX][NRPotsMAX]

private

Definition at line 95 of file CTPPSPixelDQMSource.cc.

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

unsigned int CTPPSPixelDQMSource::rpStatusWord = 0x8000

private

Definition at line 94 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::StationStatus[NStationMAX]

private

Definition at line 96 of file CTPPSPixelDQMSource.cc.

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

CTPPSPixelIndices CTPPSPixelDQMSource::thePixIndices

private

Definition at line 61 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

edm::EDGetTokenT<edm::DetSetVector<CTPPSPixelDigi> > CTPPSPixelDQMSource::tokenDigi

private

Definition at line 48 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and CTPPSPixelDQMSource().

unsigned int CTPPSPixelDQMSource::verbosity

private

Definition at line 46 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), dqmBeginRun(), and endRun().