CTPPSPixelDQMSource Class Reference

Inheritance diagram for CTPPSPixelDQMSource:

Public Member Functions
	CTPPSPixelDQMSource (const edm::ParameterSet &ps)
	~CTPPSPixelDQMSource () override

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override
void	endRun (edm::Run const &run, edm::EventSetup const &eSetup) override

Private Member Functions
int	getDet (int id)
int	getPixPlane (int id)
int	getPlaneIndex (int arm, int station, int rp, int plane)
int	getRPglobalBin (int arm, int stn)
int	getRPindex (int arm, int station, int rp)
int	getRPInStationBin (int rp)
int	prIndex (int rp, int plane)

Private Attributes
int	cluSizeMax
const int	ClusMultMAX = 10
int	ClusMultPlane [RPotsTotalNumber][NplaneMAX]
MonitorElement *	h2AllPlanesActive
MonitorElement *	h2CluSize [NArms][NStationMAX]
MonitorElement *	h2ClusMultipl [NArms][NStationMAX]
MonitorElement *	h2HitsMultipl [NArms][NStationMAX]
MonitorElement *	h2HitsMultROC [RPotsTotalNumber]
MonitorElement *	h2PlaneActive [NArms][NStationMAX]
MonitorElement *	h2xyHits [RPotsTotalNumber][NplaneMAX]
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 *	hp2HitsMultROC_LS [RPotsTotalNumber]
MonitorElement *	hp2HitsOcc [NArms][NStationMAX]
MonitorElement *	hp2xyADC [RPotsTotalNumber][NplaneMAX]
MonitorElement *	hROCadc [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	hRPotActivBX [RPotsTotalNumber]
MonitorElement *	hRPotActivBXall [RPotsTotalNumber]
MonitorElement *	hRPotActivBXroc [RPotsTotalNumber]
MonitorElement *	hRPotActivPlanes [RPotsTotalNumber]
MonitorElement *	hRPotActivROCs [RPotsTotalNumber]
const int	IndexNotValid = 0
int	multHitsMax
long int	nEvents = 0
int	pixColMAX = 156
int	pixRowMAX = 160
int	ROCSizeInX = pixRowMAX/2
int	ROCSizeInY = pixColMAX/3
int	RPindexValid [RPotsTotalNumber]
int	RPstatus [StationIDMAX][RPotsIDMAX]
unsigned int	rpStatusWord = 0x8008
int	StationStatus [StationIDMAX]
CTPPSPixelIndices	thePixIndices
edm::EDGetTokenT< edm::DetSetVector< CTPPSPixelCluster > >	tokenCluster
edm::EDGetTokenT< edm::DetSetVector< CTPPSPixelDigi > >	tokenDigi
unsigned int	verbosity

Static Private Attributes
static int	ADCMax = 256
static int	NArms =2
static int	NPlaneBins = NplaneMAX*NRPotBinsInStation
static int	NplaneMAX =6
static int	NROCsMAX = 6
static int	NRPglobalBins = 4
static int	NRPotBinsInStation = RPn_last-RPn_first
static int	NRPotsMAX =6
static int	NStationMAX =3
static int	RPn_first = 3
static int	RPn_last = 4
static int	RPotsIDMAX =8
static int	RPotsTotalNumber =NArms*NStationMAX*NRPotsMAX
static int	StationIDMAX =4

Detailed Description

Definition at line 32 of file CTPPSPixelDQMSource.cc.

Constructor & Destructor Documentation

CTPPSPixelDQMSource::CTPPSPixelDQMSource

(

const edm::ParameterSet &

ps

)

Definition at line 146 of file CTPPSPixelDQMSource.cc.

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

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

}

CTPPSPixelDQMSource::~CTPPSPixelDQMSource ( )

override

Definition at line 157 of file CTPPSPixelDQMSource.cc.

{
}

Member Function Documentation

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

overrideprotected

Definition at line 382 of file CTPPSPixelDQMSource.cc.

References ecalMGPA::adc(), CTPPSDetId::arm(), edm::DetSet< T >::begin(), edm::EventBase::bunchCrossing(), cluSizeMax, ClusMultPlane, cuy::col, MonitorElement::Fill(), HcalObjRepresent::Fill(), objects.autophobj::float, getDet(), getPixPlane(), getPlaneIndex(), getRPglobalBin(), getRPindex(), getRPInStationBin(), h2AllPlanesActive, h2CluSize, h2ClusMultipl, h2HitsMultipl, h2HitsMultROC, h2PlaneActive, h2xyHits, h2xyROCHits, hBX, hBXshort, hHitsMult, HitsMultPlane, HitsMultROC, hp2HitsMultROC_LS, hp2HitsOcc, hp2xyADC, hROCadc, hRPotActivBX, hRPotActivBXall, hRPotActivBXroc, hRPotActivPlanes, hRPotActivROCs, edm::HandleBase::isValid(), SiStripPI::max, multHitsMax, nEvents, np, NplaneMAX, NROCsMAX, NRPotsMAX, NStationMAX, AlCaHLTBitMon_ParallelJobs::p, pixColMAX, pixRowMAX, prIndex(), alignCSCRings::r, CTPPSDetId::rp(), RPindexValid, RPotsTotalNumber, RPstatus, mathSSE::sqrt(), CTPPSDetId::station(), relativeConstraints::station, StationStatus, thePixIndices, tokenCluster, tokenDigi, CTPPSPixelIndices::transformToROC(), verbosity, and DetId::VeryForward.

{
  ++nEvents;
  int lumiId = event.getLuminosityBlock().id().luminosityBlock();
  if(lumiId<0) lumiId=0;

  int RPactivity[RPotsTotalNumber], RPdigiSize[RPotsTotalNumber];

  for(int rp=0; rp<RPotsTotalNumber; rp++)
    { RPactivity[rp] = RPdigiSize[rp] = 0;}

  for(int ind=0; ind<RPotsTotalNumber; ind++) { 
    for(int p=0; p<NplaneMAX; p++) {
      HitsMultPlane[ind][p] = 0;
      ClusMultPlane[ind][p] = 0;
    }
  }
  for(int ind=0; ind<RPotsTotalNumber*NplaneMAX; ind++)  {
    for(int roc=0; roc<NROCsMAX; roc++) {
      HitsMultROC[ind][roc] = 0;
    }
  }
  Handle< DetSetVector<CTPPSPixelDigi> > pixDigi;
  event.getByToken(tokenDigi, pixDigi);

  Handle< DetSetVector<CTPPSPixelCluster> > pixClus;
  event.getByToken(tokenCluster, pixClus);

  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;
      int rpInd = getRPindex(arm,station,rpot);
      RPactivity[rpInd] = 1;
      ++RPdigiSize[rpInd];

      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,getRPInStationBin(rpot));
        h2AllPlanesActive->Fill(plane,getRPglobalBin(arm,station));

        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(colROC,rowROC);
                hROCadc[rocHistIndex][trocId]->Fill(adc);
                ++HitsMultROC[rocHistIndex][trocId];
              }
            }
          } //end if(RPindexValid[index]) {
        }
        if(int(ds_digi.data.size()) > multHitsMax) multHitsMax = ds_digi.data.size();

      } // end  if(StationStatus[station]) {
    } // end for(const auto &ds_digi : *pixDigi)
  } //if(pixDigi.isValid()) {

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

      int plane = getPixPlane(ds.id);

      CTPPSDetId theId(ds.id);
      int arm = theId.arm()&0x1;
      int station = theId.station()&0x3;
      int rpot = theId.rp()&0x7;

      if((StationStatus[station]==0) || (RPstatus[station][rpot]==0)) continue;

        int index = getRPindex(arm,station,rpot);
        ++ClusMultPlane[index][plane];

      unsigned int minRow=pixRowMAX, maxRow=0;
      unsigned int minCol=pixColMAX, maxCol=0;
      for (const auto &p : ds) {
        unsigned int max= p.minPixelRow() + p.rowSpan()+1;
        if(minRow > p.minPixelRow()) minRow = p.minPixelRow();
        if(maxRow < max) maxRow = max;
        max= p.minPixelCol() + p.colSpan()+1;
        if(minCol > p.minPixelCol()) minCol = p.minPixelCol();
        if(maxCol < max) maxCol = max;

      }
      int drow = maxRow - minRow;
      int dcol = maxCol - minCol;
      float clusize= sqrt(drow*drow + dcol*dcol);
      if(cluSizeMax < int(clusize)) cluSizeMax = clusize;
      h2CluSize[arm][station]->Fill(prIndex(rpot,plane),clusize);

   } // end if(pixClus.isValid()) for(const auto &ds : *pixClus)

  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[index]==0) continue;

        for(int p=0; p<NplaneMAX; p++)
          h2ClusMultipl[arm][stn]->Fill(prIndex(rp,p),ClusMultPlane[index][p]);

        int np = 0; 
        for(int p=0; p<NplaneMAX; p++) if(HitsMultPlane[index][p]>0) np++;
    for(int p=0; p<=NplaneMAX; p++) {
      if(p == np) hRPotActivPlanes[index]->Fill(p,1.);
      else hRPotActivPlanes[index]->Fill(p,0.);
    }
        if(np>5) hRPotActivBX[index]->Fill(event.bunchCrossing());
        hRPotActivBXall[index]->Fill(event.bunchCrossing(),float(RPdigiSize[index]));

        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]);
            hp2HitsMultROC_LS[index]->Fill(lumiId,p*NROCsMAX+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++) {
      for(int p=0; p<=NplaneMAX; p++)
        if(p==rocf[r]) hRPotActivROCs[index]->Fill(r,rocf[r],1.);
        else hRPotActivROCs[index]->Fill(r,p,0);
    }

        if(max > 4) hRPotActivBXroc[index]->Fill(event.bunchCrossing());
      } //end for(int rp=0; rp<NRPotsMAX; rp++) {
    }
  } //end for(int arm=0; arm<2; arm++) {

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

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

overrideprotected

Definition at line 195 of file CTPPSPixelDQMSource.cc.

References ADCMax, CTPPSDetId::armName(), DQMStore::IBooker::book1D(), DQMStore::IBooker::book1DD(), DQMStore::IBooker::book2D(), DQMStore::IBooker::book2DD(), DQMStore::IBooker::bookProfile(), DQMStore::IBooker::bookProfile2D(), DQMStore::IBooker::cd(), ClusMultMAX, MillePedeFileConverter_cfg::e, JetChargeProducer_cfi::exp, getPlaneIndex(), getRPglobalBin(), CTPPSDetId::getRPId(), getRPindex(), getRPInStationBin(), CTPPSDetId::getStationId(), MonitorElement::getTH2D(), MonitorElement::getTH2F(), MonitorElement::getTProfile2D(), h2AllPlanesActive, h2CluSize, h2ClusMultipl, h2HitsMultipl, h2HitsMultROC, h2PlaneActive, h2xyHits, h2xyROCHits, hBX, hBXshort, hHitsMult, hitMultMAX, hp2HitsMultROC_LS, hp2HitsOcc, hp2xyADC, hROCadc, hRPotActivBX, hRPotActivBXall, hRPotActivBXroc, hRPotActivPlanes, hRPotActivROCs, mps_fire::i, cmsBatch::log, n0, pileupCalc::nbins, CTPPSDetId::nFull, CTPPSDetId::nPath, NPlaneBins, NplaneMAX, NROCsMAX, NRPglobalBins, NRPotBinsInStation, CTPPSDetId::nShort, NStationMAX, AlCaHLTBitMon_ParallelJobs::p, pixRowMAX, alignCSCRings::r, ROCSizeInX, ROCSizeInY, RPindexValid, RPn_first, RPn_last, RPstatus, alignCSCRings::s, sd, CTPPSDetId::sdTrackingPixel, DQMStore::IBooker::setCurrentFolder(), CTPPSDetId::setRP(), CTPPSDetId::setStation(), StationStatus, and AlCaHLTBitMon_QueryRunRegistry::string.

{
  ibooker.cd();
  ibooker.setCurrentFolder("CTPPS/TrackingPixel");
  char s[50];
  string armTitleShort, stnTitleShort;
  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);

  string str1st = "Pixel planes activity";
  h2AllPlanesActive = ibooker.book2DD(str1st,str1st+";Plane #",
                NplaneMAX,0,NplaneMAX, NRPglobalBins, 0.5, NRPglobalBins+0.5);
  TH2D *h1st = h2AllPlanesActive->getTH2D();
  h1st->SetOption("colz");
  TAxis *yah1st = h1st->GetYaxis();

 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);
   ID.armName(armTitleShort, CTPPSDetId::nShort);

   ibooker.setCurrentFolder(sd);

 for(int stn=0; 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);
   CTPPSDetId(ID.getStationId()).stationName(stnTitleShort, CTPPSDetId::nShort);

   ibooker.setCurrentFolder(stnd);

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

   h2PlaneActive[arm][stn] = ibooker.book2DD(st,st+st2+";Plane #",
        NplaneMAX,0,NplaneMAX, NRPotBinsInStation, 0.5, NRPotBinsInStation+0.5);
   TH2D *h = h2PlaneActive[arm][stn]->getTH2D();
   h->SetOption("colz");
   TAxis *yah = h->GetYaxis();
   
   st = "hit average multiplicity in planes";
   int PlaneMultCut = 20;
   hp2HitsOcc[arm][stn]= ibooker.bookProfile2D(st,st+st2+";Plane #",
     NplaneMAX, 0, NplaneMAX, NRPotBinsInStation, RPn_first,RPn_last,-1,PlaneMultCut,"");
   TProfile2D *h2 = hp2HitsOcc[arm][stn]->getTProfile2D();
   h2->SetOption("textcolz");
   TAxis *yah2 = h2->GetYaxis();

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

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

   st = "cluster span in planes";
  const int nyClus = 9; //18;
  float xCh[NPlaneBins+1];
  float yClus[nyClus+1];
  for(int i=0; i<=NPlaneBins; i++) xCh[i]=i;
  double n0 = 1; //1./CluSizeMAX;
  double lnA = log(2.);
  yClus[0] = n0; //yClus[1] = n0;
  for(int j=0; j<nyClus; j++) yClus[j+1] = n0*exp(j*lnA);

  h2CluSize[arm][stn] = ibooker.book2D(st,st+st2+st3+";Cluster size",
    NPlaneBins,xCh,nyClus,yClus);
   h2CluSize[arm][stn]->getTH2F()->SetOption("colz");

//--------- RPots ---
   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(getRPInStationBin(rp), rpTitle.c_str()); // h
       yah2->SetBinLabel(getRPInStationBin(rp), rpTitle.c_str()); //h2
       string rpBinName = armTitleShort + "_" + stnTitleShort+"_"+rpTitle;
       yah1st->SetBinLabel(getRPglobalBin(arm,stn), rpBinName.c_str());

       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.bookProfile("number of fired planes per event", rpTitle+";nPlanes;Probability",
     NplaneMAX+1, -0.5, NplaneMAX+0.5, -0.5,NplaneMAX+0.5,"");
       hRPotActivBX[indexP] = 
        ibooker.book1D("5 fired planes 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,
    0.,ROCSizeInX*ROCSizeInY,"");
       h2HitsMultROC[indexP]->getTProfile2D()->SetOption("colztext");
       h2HitsMultROC[indexP]->getTProfile2D()->SetMinimum(1.e-10);

    
       hp2HitsMultROC_LS[indexP]=ibooker.bookProfile2D("ROCs_hits_multiplicity_per_event vs LS",
     rpTitle+";LumiSection;Plane#___ROC#", 1000,0.,1000.,
     NplaneMAX*NROCsMAX,0.,double(NplaneMAX*NROCsMAX),0.,ROCSizeInX*ROCSizeInY,"");
       hp2HitsMultROC_LS[indexP]->getTProfile2D()->SetOption("colz");
       hp2HitsMultROC_LS[indexP]->getTProfile2D()->SetMinimum(1.0e-10);
       hp2HitsMultROC_LS[indexP]->getTProfile2D()->SetCanExtend(TProfile2D::kXaxis);
       TAxis *yahp2 = hp2HitsMultROC_LS[indexP]->getTProfile2D()->GetYaxis();
       for(int p=0; p<NplaneMAX; p++) {
     sprintf(s,"plane%d_0",p);
     yahp2->SetBinLabel(p*NplaneMAX+1,s); 
     for(int r=1; r<NROCsMAX; r++) {
       sprintf(s,"   %d_%d",p,r);
       yahp2->SetBinLabel(p*NplaneMAX+r+1,s);
     }
       }

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

       ibooker.setCurrentFolder(rpd+"/latency");

       hRPotActivBXroc[indexP] = 
        ibooker.book1D("4 fired ROCs per BX", rpTitle+";Event.BX", 4002, -1.5, 4000.+0.5);

       hRPotActivBXall[indexP] =
        ibooker.book1D("hits per BX", rpTitle+";Event.BX", 4002, -1.5, 4000.+0.5);

       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, 0.,512.,"");
     hp2xyADC[indexP][p]->getTProfile2D()->SetOption("colz");

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

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

         for(int roc=0; roc<NROCsMAX; 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 col;pix row",
        ROCSizeInY,0,ROCSizeInY, ROCSizeInX,0,ROCSizeInX);
           h2xyROCHits[index][roc]->getTH2D()->SetOption("colz");
           hROCadc[index][roc]=ibooker.book1D("adc value",st2+";ADC;number of ROCs",
        ADCMax, 0.,float(ADCMax));
         }
       } // 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 &    )

overrideprotected

Definition at line 163 of file CTPPSPixelDQMSource.cc.

References cluSizeMax, CTPPSPixelIndices::getDefaultColDetSize(), CTPPSPixelIndices::getDefaultRowDetSize(), multHitsMax, nEvents, NRPotsMAX, NStationMAX, pixColMAX, pixRowMAX, ROCSizeInX, ROCSizeInY, RPindexValid, RPotsIDMAX, RPstatus, rpStatusWord, StationIDMAX, 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;

 for(int stn=0; stn<StationIDMAX; stn++) {
  StationStatus[stn]=0;
  for(int rp=0; rp<RPotsIDMAX; rp++) RPstatus[stn][rp]=0;
 } 

 unsigned int rpSts = rpStatusWord<<1;
 for(int stn=0; stn<NStationMAX; 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;

 multHitsMax = cluSizeMax = -1;
}

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

overrideprotected

Definition at line 569 of file CTPPSPixelDQMSource.cc.

References cluSizeMax, DEFINE_FWK_MODULE, multHitsMax, nEvents, edm::RunBase::run(), and verbosity.

{
 if(!verbosity) return; 
 LogPrint("CTPPSPixelDQMSource") 
  <<"end of Run "<<run.run()<<": "<<nEvents<<" events\n"
  <<"multHitsMax= "<<multHitsMax <<"   cluSizeMax= "<<cluSizeMax;
}

int CTPPSPixelDQMSource::getDet ( int  id )

inlineprivate

Definition at line 129 of file CTPPSPixelDQMSource.cc.

References DetId::kDetOffset.

Referenced by analyze().

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

int CTPPSPixelDQMSource::getPixPlane ( int  id )

inlineprivate

Definition at line 131 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 111 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::getRPglobalBin ( int  arm, int  stn  )

inlineprivate

Definition at line 122 of file CTPPSPixelDQMSource.cc.

References constexpr, and NStationMAX.

Referenced by analyze(), and bookHistograms().

                                       { 
   static constexpr int stationBinOrder[NStationMAX] = {0, 4, 1};
   return( arm*2 + stationBinOrder[stn] +1 );
  }

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

inlineprivate

Definition at line 104 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::getRPInStationBin ( int  rp )

inlineprivate

Definition at line 118 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

{ return(rp - RPn_first +1); }

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

inlineprivate

Definition at line 127 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

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

Member Data Documentation

int CTPPSPixelDQMSource::ADCMax = 256

staticprivate

Definition at line 56 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::cluSizeMax

private

Definition at line 135 of file CTPPSPixelDQMSource.cc.

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

const int CTPPSPixelDQMSource::ClusMultMAX = 10

private

Definition at line 60 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::ClusMultPlane[RPotsTotalNumber][NplaneMAX]

private

Definition at line 96 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

MonitorElement * CTPPSPixelDQMSource::h2AllPlanesActive

private

Definition at line 71 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 77 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 76 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 74 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2HitsMultROC[RPotsTotalNumber]

private

Definition at line 85 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 75 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 89 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 91 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hBX

private

Definition at line 71 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement * CTPPSPixelDQMSource::hBXshort

private

Definition at line 71 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 88 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

const int CTPPSPixelDQMSource::hitMultMAX = 300

private

Definition at line 59 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::HitsMultPlane[RPotsTotalNumber][NplaneMAX]

private

Definition at line 95 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 94 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

MonitorElement* CTPPSPixelDQMSource::hp2HitsMultROC_LS[RPotsTotalNumber]

private

Definition at line 86 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 73 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 90 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 92 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBX[RPotsTotalNumber]

private

Definition at line 83 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBXall[RPotsTotalNumber]

private

Definition at line 93 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBXroc[RPotsTotalNumber]

private

Definition at line 84 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivPlanes[RPotsTotalNumber]

private

Definition at line 82 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivROCs[RPotsTotalNumber]

private

Definition at line 87 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

const int CTPPSPixelDQMSource::IndexNotValid = 0

private

Definition at line 102 of file CTPPSPixelDQMSource.cc.

int CTPPSPixelDQMSource::multHitsMax

private

Definition at line 135 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::NArms =2

staticprivate

Definition at line 50 of file CTPPSPixelDQMSource.cc.

long int CTPPSPixelDQMSource::nEvents = 0

private

Definition at line 46 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::NPlaneBins = NplaneMAX*NRPotBinsInStation

staticprivate

Definition at line 69 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::NplaneMAX =6

staticprivate

Definition at line 53 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::NROCsMAX = 6

staticprivate

Definition at line 54 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::NRPglobalBins = 4

staticprivate

Definition at line 121 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::NRPotBinsInStation = RPn_last-RPn_first

staticprivate

Definition at line 68 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::NRPotsMAX =6

staticprivate

Definition at line 52 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

int CTPPSPixelDQMSource::NStationMAX =3

staticprivate

Definition at line 51 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::pixColMAX = 156

private

Definition at line 64 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

int CTPPSPixelDQMSource::pixRowMAX = 160

private

Definition at line 63 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::ROCSizeInX = pixRowMAX/2

private

Definition at line 65 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::ROCSizeInY = pixColMAX/3

private

Definition at line 66 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::RPindexValid[RPotsTotalNumber]

private

Definition at line 81 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::RPn_first = 3

staticprivate

Definition at line 55 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::RPn_last = 4

staticprivate

Definition at line 55 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::RPotsIDMAX =8

staticprivate

Definition at line 58 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::RPotsTotalNumber =NArms*NStationMAX*NRPotsMAX

staticprivate

Definition at line 79 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

int CTPPSPixelDQMSource::RPstatus[StationIDMAX][RPotsIDMAX]

private

Definition at line 100 of file CTPPSPixelDQMSource.cc.

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

unsigned int CTPPSPixelDQMSource::rpStatusWord = 0x8008

private

Definition at line 99 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::StationIDMAX =4

staticprivate

Definition at line 57 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::StationStatus[StationIDMAX]

private

Definition at line 101 of file CTPPSPixelDQMSource.cc.

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

CTPPSPixelIndices CTPPSPixelDQMSource::thePixIndices

private

Definition at line 62 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

edm::EDGetTokenT< edm::DetSetVector<CTPPSPixelCluster> > CTPPSPixelDQMSource::tokenCluster

private

Definition at line 48 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and CTPPSPixelDQMSource().

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

private

Definition at line 47 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and CTPPSPixelDQMSource().

unsigned int CTPPSPixelDQMSource::verbosity

private

Definition at line 45 of file CTPPSPixelDQMSource.cc.

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