CTPPSPixelDQMSource Class Reference

Inheritance diagram for CTPPSPixelDQMSource:

Public Member Functions
	CTPPSPixelDQMSource (const edm::ParameterSet &ps)
	~CTPPSPixelDQMSource () override
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

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
int	ClusMultPlane [RPotsTotalNumber][NplaneMAX]
MonitorElement *	h2AllPlanesActive
MonitorElement *	h2CluSize [NArms][NStationMAX]
MonitorElement *	h2HitsMultipl [NArms][NStationMAX]
MonitorElement *	h2HitsMultROC [RPotsTotalNumber]
MonitorElement *	h2trackXY0 [RPotsTotalNumber]
MonitorElement *	h2xyHits [RPotsTotalNumber][NplaneMAX]
MonitorElement *	h2xyROCHits [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	hBX
MonitorElement *	hBXshort
MonitorElement *	hHitsMult [RPotsTotalNumber][NplaneMAX]
int	HitsMultPlane [RPotsTotalNumber][NplaneMAX]
int	HitsMultROC [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	hp2HitsMultROC_LS [RPotsTotalNumber]
MonitorElement *	hp2xyADC [RPotsTotalNumber][NplaneMAX]
MonitorElement *	hpixLTrack
MonitorElement *	hpRPactive
MonitorElement *	hROCadc [RPotsTotalNumber *NplaneMAX][NROCsMAX]
MonitorElement *	hRPotActivBX [RPotsTotalNumber]
MonitorElement *	hRPotActivBXall [RPotsTotalNumber]
MonitorElement *	hRPotActivBXroc [RPotsTotalNumber]
MonitorElement *	hRPotActivPlanes [RPotsTotalNumber]
MonitorElement *	htrackHits [RPotsTotalNumber]
MonitorElement *	htrackMult [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
edm::EDGetTokenT< edm::DetSetVector< CTPPSPixelLocalTrack > >	tokenTrack
int	TrackFitDimension = 4
unsigned int	verbosity
float	x0_MAX
float	x0_MIN
float	y0_MAX
float	y0_MIN

Static Private Attributes
static int	ADCMax = 256
static int	ClusMultMAX = 10
static int	ClusterSizeMax = 9
static int	hitMultMAX = 50
static int	NArms =2
static int	NLocalTracksMAX =20
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 33 of file CTPPSPixelDQMSource.cc.

Constructor & Destructor Documentation

CTPPSPixelDQMSource::CTPPSPixelDQMSource

(

const edm::ParameterSet &

ps

)

Definition at line 152 of file CTPPSPixelDQMSource.cc.

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

                                                                  :
  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"));
tokenTrack=consumes<DetSetVector<CTPPSPixelLocalTrack>>(ps.getParameter<edm::InputTag>("tagRPixLTrack"));
}

CTPPSPixelDQMSource::~CTPPSPixelDQMSource ( )

override

Definition at line 163 of file CTPPSPixelDQMSource.cc.

{
}

Member Function Documentation

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

overrideprotected

Definition at line 378 of file CTPPSPixelDQMSource.cc.

References ecalMGPA::adc(), protons_cff::arm, CTPPSDetId::arm(), edm::EventBase::bunchCrossing(), cluSizeMax, ClusMultPlane, ClusterSizeMax, cuy::col, MonitorElement::Fill(), objects.autophobj::float, getDet(), getPixPlane(), getPlaneIndex(), getRPglobalBin(), getRPindex(), h2AllPlanesActive, h2CluSize, h2HitsMultipl, h2HitsMultROC, h2trackXY0, h2xyHits, hBX, hBXshort, hHitsMult, hitMultMAX, HitsMultPlane, HitsMultROC, hp2HitsMultROC_LS, hp2xyADC, hpixLTrack, hpRPactive, hRPotActivBX, hRPotActivBXall, hRPotActivBXroc, hRPotActivPlanes, htrackHits, htrackMult, mps_fire::i, edm::HandleBase::isValid(), SiStripPI::max, multHitsMax, nEvents, NLocalTracksMAX, np, NplaneMAX, NROCsMAX, NRPotsMAX, NStationMAX, AlCaHLTBitMon_ParallelJobs::p, prIndex(), alignCSCRings::r, CTPPSDetId::rp(), RPindexValid, RPotsTotalNumber, RPstatus, CTPPSDetId::station(), relativeConstraints::station, StationStatus, thePixIndices, tokenCluster, tokenDigi, tokenTrack, TrackFitDimension, CTPPSPixelIndices::transformToROC(), verbosity, DetId::VeryForward, x0_MAX, x0_MIN, y0_MAX, and y0_MIN.

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

  int RPactivity[RPotsTotalNumber], RPdigiSize[RPotsTotalNumber];
  int pixRPTracks[RPotsTotalNumber];

  for(int rp=0; rp<RPotsTotalNumber; rp++)
    { RPactivity[rp] = RPdigiSize[rp] = pixRPTracks[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);

  Handle< DetSetVector<CTPPSPixelLocalTrack> > pixTrack;
  event.getByToken(tokenTrack, pixTrack);

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

  if(pixTrack.isValid()) {
   for(const auto &ds_tr : *pixTrack)
   {
     int idet = getDet(ds_tr.id);
     if(idet != DetId::VeryForward) {
      if(verbosity>1) LogPrint("CTPPSPixelDQMSource") <<"not CTPPS: ds_tr.id"<<ds_tr.id;
      continue;
     }
     CTPPSDetId theId(ds_tr.id);
     int arm = theId.arm()&0x1;
     int station = theId.station()&0x3;
     int rpot = theId.rp()&0x7;  
     int rpInd = getRPindex(arm,station,rpot);

     for(DetSet<CTPPSPixelLocalTrack>::const_iterator dit = ds_tr.begin();
            dit != ds_tr.end(); ++dit) {
       ++pixRPTracks[rpInd];
       int nh_tr = (dit->getNDF() + TrackFitDimension)/2;
       for(int i=0; i<=NplaneMAX; i++) {
         if(i==nh_tr) htrackHits[rpInd]->Fill(nh_tr,1.);
         else htrackHits[rpInd]->Fill(i,0.);
       }
       float x0 = dit->getX0();
       float y0 = dit->getY0();
       h2trackXY0[rpInd]->Fill(x0,y0);
       if(x0_MAX < x0) x0_MAX = x0;
       if(y0_MAX < y0) y0_MAX = y0;
       if(x0_MIN > x0) x0_MIN = x0;
       if(y0_MIN > y0) y0_MIN = y0;
     }
   }
  } // end  if(pixTrack.isValid())


  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]) {

        h2HitsMultipl[arm][station]->Fill(prIndex(rpot,plane),ds_digi.data.size());
        h2AllPlanesActive->Fill(plane,getRPglobalBin(arm,station));

        int index = getRPindex(arm,station,rpot);
        HitsMultPlane[index][plane] += ds_digi.data.size();
        if(RPindexValid[index]) {
      int nh = ds_digi.data.size();
          if(nh > hitMultMAX) nh = hitMultMAX;
          hHitsMult[index][plane]->Fill(nh);
        }
        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) {
                ++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];

      for (const auto &p : ds) {
    int clusize = p.size(); 
    h2CluSize[arm][station]->Fill(prIndex(rpot,plane),clusize);
        if(cluSizeMax < clusize) cluSizeMax = clusize;
        if(clusize > ClusterSizeMax) clusize = ClusterSizeMax;
      }
   } // end if(pixClus.isValid()) for(const auto &ds : *pixClus)

  bool allRPactivity = false;
  for(int rp=0; rp<RPotsTotalNumber; rp++) if(RPactivity[rp]>0) allRPactivity=true; 
  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;
        hpRPactive->Fill(getRPglobalBin(arm,stn),RPactivity[index]);
//        if(RPactivity[index]==0) continue;
          if(!allRPactivity) continue;
        hpixLTrack->Fill(getRPglobalBin(arm,stn),pixRPTracks[index]);
        int ntr = pixRPTracks[index];
        if(ntr > NLocalTracksMAX) ntr = NLocalTracksMAX;
    for(int i=0; i<=NLocalTracksMAX; i++) {
     if(i==ntr) htrackMult[index]->Fill(ntr,1.);
     else htrackMult[index]->Fill(i,0.);
    }

        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]; }
        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 206 of file CTPPSPixelDQMSource.cc.

References protons_cff::arm, CTPPSDetId::armName(), DQMStore::IBooker::book1D(), DQMStore::IBooker::book1DD(), DQMStore::IBooker::book2D(), DQMStore::IBooker::book2DD(), DQMStore::IBooker::bookProfile(), DQMStore::IBooker::bookProfile2D(), DQMStore::IBooker::cd(), ClusterSizeMax, MillePedeFileConverter_cfg::e, getRPglobalBin(), CTPPSDetId::getRPId(), getRPindex(), CTPPSDetId::getStationId(), MonitorElement::getTH2D(), MonitorElement::getTH2F(), MonitorElement::getTProfile(), MonitorElement::getTProfile2D(), h2AllPlanesActive, h2CluSize, h2HitsMultipl, h2HitsMultROC, h2trackXY0, h2xyHits, hBX, hBXshort, hHitsMult, hitMultMAX, hp2HitsMultROC_LS, hp2xyADC, hpixLTrack, hpRPactive, hRPotActivBX, hRPotActivBXall, hRPotActivBXroc, hRPotActivPlanes, htrackHits, htrackMult, pileupCalc::nbins, CTPPSDetId::nFull, NLocalTracksMAX, CTPPSDetId::nPath, NPlaneBins, NplaneMAX, NROCsMAX, NRPglobalBins, 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+"(digi task);Plane #",
                NplaneMAX,0,NplaneMAX, NRPglobalBins, 0.5, NRPglobalBins+0.5);
  TH2D *h1st = h2AllPlanesActive->getTH2D();
  h1st->SetOption("colz");
  TAxis *yah1st = h1st->GetYaxis();

  string str2 = "Pixel RP active";
  hpRPactive = ibooker.bookProfile(str2,str2+" per event(digi task)",
        NRPglobalBins, 0.5,NRPglobalBins+0.5, -0.1, 1.1,"" );
  TAxis *xaRPact = hpRPactive->getTProfile()->GetXaxis();
  hpRPactive->getTProfile()->SetOption("hist");
  hpRPactive->getTProfile()->SetMinimum(0.);
  hpRPactive->getTProfile()->SetMaximum(1.1);

  str2 = "Pixel Local Tracks";
  hpixLTrack = ibooker.bookProfile(str2,str2+" per event",
        NRPglobalBins, 0.5, NRPglobalBins+0.5,-0.1,NLocalTracksMAX,"");

  TAxis *xah1trk = hpixLTrack->getTProfile()->GetXaxis();
  hpixLTrack->getTProfile()->GetYaxis()->SetTitle("average number of tracks per event");
  hpixLTrack->getTProfile()->SetOption("hist");

 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);
//--------- 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);
       string rpBinName = armTitleShort + "_" + stnTitleShort+"_"+rpTitle;
       yah1st->SetBinLabel(getRPglobalBin(arm,stn), rpBinName.c_str());
       xah1trk->SetBinLabel(getRPglobalBin(arm,stn), rpBinName.c_str());
       xaRPact->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);

   string st2 = ": " + stnTitle;

   string 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 size in planes";
  h2CluSize[arm][stn] = ibooker.book2D(st,st+st2+st3+";Cluster size",
    NPlaneBins,0,NPlaneBins, ClusterSizeMax+1,0,ClusterSizeMax+1);  
   h2CluSize[arm][stn]->getTH2F()->SetOption("colz");

   const float x0Maximum = 70.;
   const float y0Maximum =  15.;
   st = "track intercept point";
   h2trackXY0[indexP] = ibooker.book2D(st,st+st2+";x0;y0",
    int(x0Maximum)*2, 0.,x0Maximum,int(y0Maximum)*4,-y0Maximum,y0Maximum);
   h2trackXY0[indexP]->getTH2F()->SetOption("colz");

       st = "number of tracks per event";
       htrackMult[indexP] = ibooker.bookProfile(st,rpTitle+";number of tracks",
     NLocalTracksMAX+1, -0.5,NLocalTracksMAX+0.5, -0.5,NLocalTracksMAX+0.5,"");
       htrackMult[indexP]->getTProfile()->SetOption("hist");

       st = "number of hits per track";
       htrackHits[indexP] = ibooker.bookProfile(st,rpTitle+";number of hits",
     5,1.5,6.5, -0.1,1.1, "");
       htrackHits[indexP]->getTProfile()->SetOption("hist");

       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,"");
       hRPotActivPlanes[indexP]->getTProfile()->SetOption("hist");

       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);
     }
       }

       ibooker.setCurrentFolder(rpd+"/latency");
       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);

       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",
           pixRowMAX,0,pixRowMAX, pixRowMAX,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+1,-0.5,hitMultMAX+0.5);
       } // 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 169 of file CTPPSPixelDQMSource.cc.

References cluSizeMax, CTPPSPixelLocalTrack::dimension, CTPPSPixelIndices::getDefaultColDetSize(), CTPPSPixelIndices::getDefaultRowDetSize(), multHitsMax, nEvents, NRPotsMAX, NStationMAX, pixColMAX, pixRowMAX, ROCSizeInX, ROCSizeInY, RPindexValid, RPotsIDMAX, RPstatus, rpStatusWord, StationIDMAX, StationStatus, thePixIndices, TrackFitDimension, verbosity, x0_MAX, x0_MIN, y0_MAX, and y0_MIN.

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

  CTPPSPixelLocalTrack thePixelLocalTrack;
  TrackFitDimension = thePixelLocalTrack.dimension;

  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;
 x0_MIN = y0_MIN =  1.0e06;
 x0_MAX = y0_MAX = -1.0e06;
}

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

overrideprotected

Definition at line 591 of file CTPPSPixelDQMSource.cc.

References cluSizeMax, DEFINE_FWK_MODULE, multHitsMax, nEvents, edm::RunBase::run(), verbosity, x0_MAX, x0_MIN, y0_MAX, and y0_MIN.

{
 if(!verbosity) return; 
 LogPrint("CTPPSPixelDQMSource") 
  <<"end of Run "<<run.run()<<": "<<nEvents<<" events\n"
  <<"multHitsMax= "<<multHitsMax <<"   cluSizeMax= "<<cluSizeMax
  <<"\nx0: "<< x0_MIN<<"/" << x0_MAX<<"y0: "<<y0_MIN<<"/"<<y0_MAX<<"\n";
}

int CTPPSPixelDQMSource::getDet ( int  id )

inlineprivate

Definition at line 135 of file CTPPSPixelDQMSource.cc.

References DetId::kDetOffset.

Referenced by analyze().

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

int CTPPSPixelDQMSource::getPixPlane ( int  id )

inlineprivate

Definition at line 137 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 117 of file CTPPSPixelDQMSource.cc.

References getRPindex().

Referenced by analyze().

                                                             {
    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 128 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 110 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 124 of file CTPPSPixelDQMSource.cc.

{ return(rp - RPn_first +1); }

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

inlineprivate

Definition at line 133 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

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

Member Data Documentation

int CTPPSPixelDQMSource::ADCMax = 256

staticprivate

Definition at line 58 of file CTPPSPixelDQMSource.cc.

int CTPPSPixelDQMSource::cluSizeMax

private

Definition at line 141 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::ClusMultMAX = 10

staticprivate

Definition at line 63 of file CTPPSPixelDQMSource.cc.

int CTPPSPixelDQMSource::ClusMultPlane[RPotsTotalNumber][NplaneMAX]

private

Definition at line 102 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

int CTPPSPixelDQMSource::ClusterSizeMax = 9

staticprivate

Definition at line 64 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement * CTPPSPixelDQMSource::h2AllPlanesActive

private

Definition at line 77 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 81 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 80 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2HitsMultROC[RPotsTotalNumber]

private

Definition at line 92 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::h2trackXY0[RPotsTotalNumber]

private

Definition at line 86 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 95 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 97 of file CTPPSPixelDQMSource.cc.

MonitorElement* CTPPSPixelDQMSource::hBX

private

Definition at line 77 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement * CTPPSPixelDQMSource::hBXshort

private

Definition at line 77 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 94 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::hitMultMAX = 50

staticprivate

Definition at line 62 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::HitsMultPlane[RPotsTotalNumber][NplaneMAX]

private

Definition at line 101 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

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

private

Definition at line 100 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

MonitorElement* CTPPSPixelDQMSource::hp2HitsMultROC_LS[RPotsTotalNumber]

private

Definition at line 93 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 96 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement * CTPPSPixelDQMSource::hpixLTrack

private

Definition at line 77 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hpRPactive

private

Definition at line 78 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 98 of file CTPPSPixelDQMSource.cc.

MonitorElement* CTPPSPixelDQMSource::hRPotActivBX[RPotsTotalNumber]

private

Definition at line 90 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBXall[RPotsTotalNumber]

private

Definition at line 99 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivBXroc[RPotsTotalNumber]

private

Definition at line 91 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::hRPotActivPlanes[RPotsTotalNumber]

private

Definition at line 89 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::htrackHits[RPotsTotalNumber]

private

Definition at line 88 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

MonitorElement* CTPPSPixelDQMSource::htrackMult[RPotsTotalNumber]

private

Definition at line 87 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

const int CTPPSPixelDQMSource::IndexNotValid = 0

private

Definition at line 108 of file CTPPSPixelDQMSource.cc.

int CTPPSPixelDQMSource::multHitsMax

private

Definition at line 141 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::NArms =2

staticprivate

Definition at line 52 of file CTPPSPixelDQMSource.cc.

long int CTPPSPixelDQMSource::nEvents = 0

private

Definition at line 47 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::NLocalTracksMAX =20

staticprivate

Definition at line 61 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::NPlaneBins = NplaneMAX*NRPotBinsInStation

staticprivate

Definition at line 75 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::NplaneMAX =6

staticprivate

Definition at line 55 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::NROCsMAX = 6

staticprivate

Definition at line 56 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and bookHistograms().

int CTPPSPixelDQMSource::NRPglobalBins = 4

staticprivate

Definition at line 127 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::NRPotBinsInStation = RPn_last-RPn_first

staticprivate

Definition at line 74 of file CTPPSPixelDQMSource.cc.

int CTPPSPixelDQMSource::NRPotsMAX =6

staticprivate

Definition at line 54 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

int CTPPSPixelDQMSource::NStationMAX =3

staticprivate

Definition at line 53 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::pixColMAX = 156

private

Definition at line 68 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::pixRowMAX = 160

private

Definition at line 67 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::ROCSizeInX = pixRowMAX/2

private

Definition at line 69 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::ROCSizeInY = pixColMAX/3

private

Definition at line 70 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms(), and dqmBeginRun().

int CTPPSPixelDQMSource::RPindexValid[RPotsTotalNumber]

private

Definition at line 85 of file CTPPSPixelDQMSource.cc.

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

int CTPPSPixelDQMSource::RPn_first = 3

staticprivate

Definition at line 57 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::RPn_last = 4

staticprivate

Definition at line 57 of file CTPPSPixelDQMSource.cc.

Referenced by bookHistograms().

int CTPPSPixelDQMSource::RPotsIDMAX =8

staticprivate

Definition at line 60 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::RPotsTotalNumber =NArms*NStationMAX*NRPotsMAX

staticprivate

Definition at line 83 of file CTPPSPixelDQMSource.cc.

Referenced by analyze().

int CTPPSPixelDQMSource::RPstatus[StationIDMAX][RPotsIDMAX]

private

Definition at line 106 of file CTPPSPixelDQMSource.cc.

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

unsigned int CTPPSPixelDQMSource::rpStatusWord = 0x8008

private

Definition at line 105 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::StationIDMAX =4

staticprivate

Definition at line 59 of file CTPPSPixelDQMSource.cc.

Referenced by dqmBeginRun().

int CTPPSPixelDQMSource::StationStatus[StationIDMAX]

private

Definition at line 107 of file CTPPSPixelDQMSource.cc.

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

CTPPSPixelIndices CTPPSPixelDQMSource::thePixIndices

private

Definition at line 66 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

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

private

Definition at line 49 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and CTPPSPixelDQMSource().

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

private

Definition at line 48 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and CTPPSPixelDQMSource().

edm::EDGetTokenT< edm::DetSetVector<CTPPSPixelLocalTrack> > CTPPSPixelDQMSource::tokenTrack

private

Definition at line 50 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and CTPPSPixelDQMSource().

int CTPPSPixelDQMSource::TrackFitDimension = 4

private

Definition at line 72 of file CTPPSPixelDQMSource.cc.

Referenced by analyze(), and dqmBeginRun().

unsigned int CTPPSPixelDQMSource::verbosity

private

Definition at line 46 of file CTPPSPixelDQMSource.cc.

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

float CTPPSPixelDQMSource::x0_MAX

private

Definition at line 142 of file CTPPSPixelDQMSource.cc.

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

float CTPPSPixelDQMSource::x0_MIN

private

Definition at line 142 of file CTPPSPixelDQMSource.cc.

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

float CTPPSPixelDQMSource::y0_MAX

private

Definition at line 142 of file CTPPSPixelDQMSource.cc.

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

float CTPPSPixelDQMSource::y0_MIN

private

Definition at line 142 of file CTPPSPixelDQMSource.cc.

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