---

VisCuTrackerRecHit Class Reference

#include <VisReco/VisCustomTracker/interface/VisCuTrackerRecHit.h>

Inheritance diagram for VisCuTrackerRecHit:

List of all members.

Public Member Functions
void	drawrechit (VisCuTkModule *mod, int nlay, QPainter *p)
void	setEvent (const edm::Event &event, const edm::EventSetup &eventSetup)
void	setMap (VisCuTkMap2D *m)
void	setTracker (VisCuCmsTracker *tr)
SoSeparator *	update (VisTrackingGeometry::DetectorRep type)
std::string	updatetext (VisCuTkModule *mod)
	VisCuTrackerRecHit ()
Private Attributes
bool	firsteventloaded
edm::ESHandle< TrackerGeometry >	geom
bool	horizontal_view
VisCuTkMap2D *	map
std::string	pixelrechitProducerName
std::vector< std::pair < LocalPoint, const GeomDetUnit * > >	rdigi
std::string	striprechitProducerName
VisCuCmsTracker *	tracker

---

Detailed Description

Definition at line 25 of file VisCuTrackerRecHit.h.

---

Constructor & Destructor Documentation

VisCuTrackerRecHit::VisCuTrackerRecHit

(

)

Definition at line 51 of file VisCuTrackerRecHit.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), firsteventloaded, edm::ParameterSet::getUntrackedParameter(), horizontal_view, p, pixelrechitProducerName, VisConfigurationService::pSet(), python::EventAnalyzer_cfg::pset, striprechitProducerName, and VisConfigurationService_cff::VisConfigurationService.

{
  horizontal_view =true;
  firsteventloaded=false;
  const edm::ParameterSet *pset=VisConfigurationService::pSet();
  if(!pset)
    {
      edm::ParameterSet p;
      VisConfigurationService *visService = new VisConfigurationService(p);
      pset = visService->pSet();
      delete visService;
    }
  striprechitProducerName = pset->getUntrackedParameter<std::string>("striprechitProducer", "siStripMatchedRecHits");
  cout <<"striprechitProducer="<< striprechitProducerName <<  endl;
  pixelrechitProducerName = pset->getUntrackedParameter<std::string>("pixelrechitProducer", "siPixelRecHits");
  cout <<"pixelrechitProducer="<< pixelrechitProducerName <<  endl;
}

---

Member Function Documentation

void VisCuTrackerRecHit::drawrechit	(	VisCuTkModule *	mod,
		int	nlay,
		QPainter *	p
	)

Definition at line 339 of file VisCuTrackerRecHit.cc.

References a, VisCuTkModule::bufvalue, funct::cos(), firsteventloaded, VisCuTkSubLayer::getId(), VisCuTkModule::getId(), VisCuTkModule::getOwner(), cmsRelvalreport::green(), horizontal_view, VisCuTkSubLayer::isStereo(), j, VisCuTkModule::length, phi, VisCuTkMap2D::phival(), VisCuTkModule::posx, VisCuTkModule::posy, VisCuTkModule::posz, pv1, r, rdigi, VisCuTkModule::rechit_size, VisCuTkModule::rechit_start, funct::sin(), sl, funct::sqrt(), VisCuTkModule::value, VisCuTkModule::width, VisCuTkModule::widthAtHalfLength, x, y, and z.

Referenced by VisCuTkMap2D::paintEvent(), and VisCuTkMap2D::print().

{
  int x,y;
  double phi,r,dx,dy,dy1;
  double xp[4],yp[4],xp1,yp1;
  double vhbot,vhtop,vhapo;
  double rmedio[]={0.041,0.0701,0.0988,0.255,0.340,0.430,0.520,0.610,0.696,0.782,0.868,0.965,1.080};
  double xt1,yt1,xs1=0.,ys1=0.,xt2,yt2,xs2,ys2,pv1,pv2;
  unsigned int numberOfDigis = 0;
  
  if(firsteventloaded){
    QPointArray a =  QPointArray(4);
    phi = VisCuTkMap2D::phival(mod->posx,mod->posy);
    r = sqrt(mod->posx*mod->posx+mod->posy*mod->posy);
    vhbot = mod->width;           
    vhtop=mod->width;     
    vhapo=mod->length;
    if(nlay < 31){
      vhbot = mod->widthAtHalfLength/2.-(mod->width/2.-mod->widthAtHalfLength/2.);        
      vhtop=mod->width/2.;        
      vhapo=mod->length/2.;
      if(nlay >12 && nlay <19){
        if(map->separated)r = r+r;
        xp[0]=r-vhtop;yp[0]=-vhapo;
        xp[1]=r+vhtop;yp[1]=-vhapo;
        xp[2]=r+vhtop;yp[2]=vhapo;
        xp[3]=r-vhtop;yp[3]=vhapo;
      }else{
        if(map->separated)r = r + r/3.;
        xp[0]=r-vhapo;yp[0]=-vhbot;
        xp[1]=r+vhapo;yp[1]=-vhtop;
        xp[2]=r+vhapo;yp[2]=vhtop;
        xp[3]=r-vhapo;yp[3]=vhbot;
      }
      for(int j=0;j<4;j++){
        xp1 = xp[j]*cos(phi)-yp[j]*sin(phi);
        yp1 = xp[j]*sin(phi)+yp[j]*cos(phi);
        xp[j] = xp1;yp[j]=yp1;
      }
      for(int idigi=mod->rechit_start;idigi!=(mod->rechit_start+mod->rechit_size);idigi++){ 
        if(!map->separated){
          LocalPoint point = rdigi[idigi].first;
          float x = rdigi[idigi].second->surface().toGlobal(point).x()/100;
          float y = rdigi[idigi].second->surface().toGlobal(point).y()/100;
          int xc = map->xpixel(x);
          int yc = map->ypixel(y);
          if(!map->acc){
            p->setPen(QColor(255,0,255));    
            p->setBrush(QColor(255,0,255));
            if(!horizontal_view)p->drawRect(xc,yc,3,3);else p->drawRect(yc,map->x_win-xc,3,3); 
          } else {
              if(!horizontal_view){
              if(xc > (map->ximg -1)) xc =  (map->ximg -1);
              if(yc > (map->yimg -1)) yc =(map->yimg-1);
                } else {
              if(yc > (map->ximg -1)) yc =  (map->ximg -1);
              if(xc > (map->yimg -1)) xc =(map->yimg-1);
                }
              QRgb pixel;
              if(!horizontal_view){
                pixel = map->image->pixel(xc,yc);
              }else pixel = map->image->pixel(yc,map->x_win-xc);
              int green = qGreen(pixel);
              if(green > 1) green = green -20;
              if(green < 0) green = 0;
              if(!horizontal_view)map->image->setPixel(xc,yc,green);else map->image->setPixel(yc,map->x_win-xc,green);
          }
        } 
        numberOfDigis++;
      } //for digis 
    } else { //barrel
      int numod;
      numod=mod->getId();if(numod>100)numod=numod-100;
      int vane = mod->getOwner()->getId();
      if(map->separated){
        dx = vhapo;
        phi=M_PI;
        xt1=rmedio[nlay-31]; yt1=-vhtop/2.;
        xs1 = xt1*cos(phi)-yt1*sin(phi);
        ys1 = xt1*sin(phi)+yt1*cos(phi);
        xt2=rmedio[nlay-31]; yt2=vhtop/2.;
        xs2 = xt2*cos(phi)-yt2*sin(phi);
        ys2 = xt2*sin(phi)+yt2*cos(phi);
        dy=VisCuTkMap2D::phival(xs2,ys2)-VisCuTkMap2D::phival(xs1,ys1);
        dy1 = dy;
        if(nlay==31)dy1=0.39;
        if(nlay==32)dy1=0.23;
        if(nlay==33)dy1=0.16;
        xp[0]=vane*(dx+dx/8.);yp[0]=numod*(dy1);
        xp[1]=vane*(dx+dx/8.)+dx;yp[1]=numod*(dy1);
        xp[2]=vane*(dx+dx/8.)+dx;yp[2]=numod*(dy1)+dy;
        xp[3]=vane*(dx+dx/8.);yp[3]=numod*(dy1)+dy;
      }else{
        xt1=r; yt1=-vhtop/2.;
        xs1 = xt1*cos(phi)-yt1*sin(phi);
        ys1 = xt1*sin(phi)+yt1*cos(phi);
        xt2=r; yt2=vhtop/2.;
        xs2 = xt2*cos(phi)-yt2*sin(phi);
        ys2 = xt2*sin(phi)+yt2*cos(phi);
        pv1=VisCuTkMap2D::phival(xs1,ys1);
        pv2=VisCuTkMap2D::phival(xs2,ys2);
        if(fabs(pv1-pv2)>M_PI && numod==1)pv1=pv1-2.*M_PI;
        if(fabs(pv1-pv2)>M_PI && numod!=1)pv2=pv2+2.*M_PI;
        xp[0]=mod->posz-vhapo/2.;yp[0]=4.2*pv1;
        xp[1]=mod->posz+vhapo/2.;yp[1]=4.2*pv1;
        xp[2]=mod->posz+vhapo/2. ;yp[2]=4.2*pv2;
        xp[3]=mod->posz-vhapo/2.;yp[3]=4.2*pv2;
      }
      for(int idigi=mod->rechit_start;idigi!=(mod->rechit_start+mod->rechit_size);idigi++){ 
        if(!map->separated){
          LocalPoint point = rdigi[idigi].first;
          float x = rdigi[idigi].second->surface().toGlobal(point).x()/100;
          float y = rdigi[idigi].second->surface().toGlobal(point).y()/100;
          float z = rdigi[idigi].second->surface().toGlobal(point).z()/100;
          int xc = map->xpixel(z);
          int yc = map->ypixel(4.2*VisCuTkMap2D::phival(x,y));
          if(!map->acc){
            p->setPen(QColor(255,0,0));//red    
            p->setBrush(QColor(255,0,0));
            if(!horizontal_view)p->drawRect(xc,yc,3,3);else p->drawRect(yc,map->x_win-xc,3,3);
          } else {
              if(!horizontal_view){
              if(xc > (map->ximg -1)) xc =  (map->ximg -1);
              if(yc > (map->yimg -1)) yc =(map->yimg-1);
                } else {
              if(yc > (map->ximg -1)) yc =  (map->ximg -1);
              if(xc > (map->yimg -1)) xc =(map->yimg-1);
                }
              QRgb pixel;
              if(!horizontal_view){
                pixel = map->image->pixel(xc,yc);
              }else pixel = map->image->pixel(yc,xc);
              int green = qGreen(pixel);
              if(green > 1) green = green -20;
              if(green < 0) green = 0;
              if(!horizontal_view)map->image->setPixel(xc,yc,green);else map->image->setPixel(yc,map->x_win-xc,green);
          }
        }//not separated
        numberOfDigis++;        
      }//digis
    }//endcap
    VisCuTkSubLayer* sl = mod->getOwner();
    if(sl->isStereo()&& map->separated)
      {
        if(mod->getId()>100 ){
          for(int j=0;j<3;j++){
            x=map->xpixel(xp[j]);y=map->ypixel(yp[j]);
            if(!horizontal_view)a.setPoint(j,x,y);else a.setPoint(j,y,map->x_win-x);
          }
          if(!horizontal_view)a.setPoint(3,x,y);else a.setPoint(3,y,map->x_win-x);
        }else {
          x=map->xpixel(xp[2]);y=map->ypixel(yp[2]);
          if(!horizontal_view)a.setPoint(0,x,y);else a.setPoint(0,y,map->x_win-x);
          x=map->xpixel(xp[3]);y=map->ypixel(yp[3]);
          if(!horizontal_view)a.setPoint(1,x,y);else a.setPoint(1,y,map->x_win-x);
          x=map->xpixel(xp[0]);y=map->ypixel(yp[0]);
          if(!horizontal_view)a.setPoint(2,x,y);else a.setPoint(2,y,map->x_win-x);
          if(!horizontal_view)a.setPoint(3,x,y);else a.setPoint(3,y,map->x_win-x);
        }
      } else {
        for(int j=0;j<4;j++){
          x=map->xpixel(xp[j]);y=map->ypixel(yp[j]);
          if(!horizontal_view)a.setPoint(j,x,y);else a.setPoint(j,y,map->x_win-x);
        }
      }                 
    if(map->acc){mod->value = mod->value + mod->rechit_size;mod->bufvalue = mod->bufvalue + mod->rechit_size;}
    else mod->value = mod->rechit_size; 
    if(map->separated){
      p->setPen(QColor( 158, 185, 225));
      int green = 255-mod->value; if (green < 0)green = 0;if(green > 255)green=255;
      p->setBrush(QColor(255,green,0));
      if(mod->value==0)p->setBrush(QColor(255,255,255));
      p->drawPolygon(a);
    }
    p->setPen(QColor(0,0,0));    
    p->setBrush(Qt::NoBrush); 
  }//firsteventloaded
}

void VisCuTrackerRecHit::setEvent	(	const edm::Event &	event,
		const edm::EventSetup &	eventSetup
	)

Definition at line 70 of file VisCuTrackerRecHit.cc.

References begin, GenMuonPlsPt100GeV_cfg::cout, VisCuTkModule::detUnit, lat::endl(), firsteventloaded, GeomDet::geographicalId(), geom, edm::EventSetup::get(), edm::RangeMap< ID, C, P >::get(), i, edm::RangeMap< ID, C, P >::ids(), iter, BaseSiTrackerRecHit2DLocalPos::localPosition(), mod(), VisCuTkModuleMap::moduleMap, pixelrechitProducerName, rdigi, VisCuTkModule::rechit_size, VisCuTkModule::rechit_start, EgammaValidation_Zee_cff::sel, edm::RangeMap< ID, C, P >::size(), striprechitProducerName, and DetId::subdetId().

Referenced by VisCuTkGeometryTwig::onNewEvent().

{
  eventSetup.get<TrackerDigiGeometryRecord> ().get (geom);
  
  edm::Handle<SiStripMatchedRecHit2DCollection> rechitsmatched;
  edm::Handle<SiStripRecHit2DCollection> rechitsrphi;
  edm::Handle<SiStripRecHit2DCollection> rechitsstereo;
  rdigi.clear();
  int rd_pointer=0;
  bool norechits=true;
  firsteventloaded=true; 
if(striprechitProducerName!="ignore"){
 std::vector< edm::Handle<SiStripMatchedRecHit2DCollection> > rechit2d_matched_collections;
    ModuleLabelSelector sel(striprechitProducerName);
    event.getMany (sel, rechit2d_matched_collections);
 std::vector< edm::Handle<SiStripRecHit2DCollection> > rechit2d_collections;
    event.getMany (sel, rechit2d_collections);
      if (! rechit2d_matched_collections.empty () && !rechit2d_collections.empty ())norechits=false;
      if(norechits) cout << "no striprechits type "<<striprechitProducerName<<" in event" << endl;
      if(!norechits){
      event.getByLabel(striprechitProducerName,"matchedRecHit", rechitsmatched);
      event.getByLabel(striprechitProducerName,"rphiRecHit", rechitsrphi);
      event.getByLabel(striprechitProducerName,"stereoRecHit", rechitsstereo);
    if(rechitsmatched->size()>0||rechitsrphi->size()>0||rechitsstereo->size()>0){
      std::vector<std::pair<const GeomDetUnit *,std::pair<LocalPoint, const GeomDetUnit *> > > rtemp;
      if(rechitsmatched->size()>0){ 
      const std::vector<DetId> detIDs = rechitsmatched->ids();
      for ( std::vector<DetId>::const_iterator detunit_iterator = detIDs.begin(); detunit_iterator != detIDs.end(); detunit_iterator++ )
        {
          unsigned int id = (*detunit_iterator).rawId();
          if(id!=999999999){ //if is valid detector
            VisCuTkModule * mod = 0, *mod1 = 0,*mod2 = 0;
            DetId detIdObject(id);
            DetId detIdObject1(id+1);
            DetId detIdObject2(id+2);
            const GeomDetUnit * geoUnit,*geoUnit1,*geoUnit2;
            int isStereo=id&0x3;
            if(isStereo!=0){
              geoUnit = geom->idToDetUnit( detIdObject );
              mod = VisCuTkModuleMap::moduleMap[geoUnit];
            }else{
              geoUnit1 = geom->idToDetUnit( detIdObject1 );
              geoUnit2 = geom->idToDetUnit( detIdObject2 );
              mod1 = VisCuTkModuleMap::moduleMap[geoUnit1];
              mod2 = VisCuTkModuleMap::moduleMap[geoUnit2];
            }
            SiStripMatchedRecHit2DCollection::range rechitRange = rechitsmatched->get((*detunit_iterator));
            SiStripMatchedRecHit2DCollection::const_iterator rechitRangeIteratorBegin = rechitRange.first;
            SiStripMatchedRecHit2DCollection::const_iterator rechitRangeIteratorEnd   = rechitRange.second;
            SiStripMatchedRecHit2DCollection::const_iterator iter=rechitRangeIteratorBegin;
            
            for(iter=rechitRangeIteratorBegin;iter!=rechitRangeIteratorEnd;++iter){//loop on the rechit
              SiStripMatchedRecHit2D const rechit=*iter;
              LocalPoint position=rechit.localPosition();
              if(isStereo!=0){
                std::pair<const GeomDetUnit*, std::pair<LocalPoint, const GeomDetUnit* > > pxd;
                pxd.first = geoUnit;
                pxd.second.first = position;
                pxd.second.second = geoUnit;
                rtemp.push_back(pxd);
                if(mod){ mod->rechit_size++; } else cout << "module not found " << endl; 
              }else{
                std::pair<const GeomDetUnit*, std::pair<LocalPoint, const GeomDetUnit* > > pxd1;
                pxd1.first = geoUnit1;
                pxd1.second.first = position;
                pxd1.second.second = geoUnit2;
                rtemp.push_back(pxd1);
                if(mod1){ mod1->rechit_size++; } else cout << "module not found " << endl; 
                std::pair<const GeomDetUnit*, std::pair<LocalPoint, const GeomDetUnit*  > > pxd2;
                pxd2.first = geoUnit2;
                pxd2.second.first = position;
                pxd2.second.second = geoUnit2;
                rtemp.push_back(pxd2);
                if(mod2){ mod2->rechit_size++; } else cout << "module not found " << endl; 
              }
            }
          }
        } 
      } 
      if(rechitsrphi->size()>0){ 
        const std::vector<DetId> detIDs1 = rechitsrphi->ids();
        // loop over detunits
        for ( std::vector<DetId>::const_iterator detunit_iterator = detIDs1.begin(); detunit_iterator != detIDs1.end(); detunit_iterator++ )
          {
            unsigned int id = (*detunit_iterator).rawId();
            if(id!=999999999){ //if is valid detector
              DetId detIdObject(id);
              const GeomDetUnit      * geoUnit = geom->idToDetUnit( detIdObject );
              VisCuTkModule * mod = VisCuTkModuleMap::moduleMap[geoUnit];
              SiStripRecHit2DCollection::range rechitRange = rechitsrphi->get((*detunit_iterator));
              SiStripRecHit2DCollection::const_iterator rechitRangeIteratorBegin = rechitRange.first;
              SiStripRecHit2DCollection::const_iterator rechitRangeIteratorEnd   = rechitRange.second;
              SiStripRecHit2DCollection::const_iterator iter=rechitRangeIteratorBegin;
              for(iter=rechitRangeIteratorBegin;iter!=rechitRangeIteratorEnd;++iter){//loop on the rechit
                SiStripRecHit2D const rechit=*iter;
                LocalPoint position=rechit.localPosition();
                std::pair<const GeomDetUnit*, std::pair<LocalPoint, const GeomDetUnit* > > pxd;
                pxd.first = geoUnit;
                pxd.second.first = position;
                pxd.second.second = geoUnit;
                rtemp.push_back(pxd);
                if(mod){ mod->rechit_size++; } else cout << "module not found " << endl; 
              }
            }
          } 
      } 
      if(rechitsstereo->size()>0){ 
        const std::vector<DetId> detIDs2 = rechitsstereo->ids();
        // loop over detunits
        for ( std::vector<DetId>::const_iterator detunit_iterator = detIDs2.begin(); detunit_iterator != detIDs2.end(); detunit_iterator++ )
          {
            unsigned int id = (*detunit_iterator).rawId();
            if(id!=999999999){ //if is valid detector
              DetId detIdObject(id);
              const GeomDetUnit      * geoUnit = geom->idToDetUnit( detIdObject );
              VisCuTkModule * mod = VisCuTkModuleMap::moduleMap[geoUnit];
              SiStripRecHit2DCollection::range rechitRange = rechitsstereo->get((*detunit_iterator));
              SiStripRecHit2DCollection::const_iterator rechitRangeIteratorBegin = rechitRange.first;
              SiStripRecHit2DCollection::const_iterator rechitRangeIteratorEnd   = rechitRange.second;
              SiStripRecHit2DCollection::const_iterator iter=rechitRangeIteratorBegin;
              for(iter=rechitRangeIteratorBegin;iter!=rechitRangeIteratorEnd;++iter){//loop on the rechit
                SiStripRecHit2D const rechit=*iter;
                LocalPoint position=rechit.localPosition();
                std::pair<const GeomDetUnit*, std::pair<LocalPoint, const GeomDetUnit* > > pxd;
                pxd.first = geoUnit;
                pxd.second.first = position;
                pxd.second.second = geoUnit;
                rtemp.push_back(pxd);
                if(mod){ mod->rechit_size++; } else cout << "module not found " << endl; 
              }
            }
          } 
      } 
      cout << rtemp.size() << " rechits in sistrip " << endl; 
      
      std::map<const GeomDetUnit * , VisCuTkModule *>::iterator imod;
      for (imod=VisCuTkModuleMap::moduleMap.begin();imod != VisCuTkModuleMap::moduleMap.end(); imod++){
        VisCuTkModule* mod = imod->second;
        if(mod->rechit_size>0){
          mod->rechit_start=rd_pointer; rd_pointer = rd_pointer + mod->rechit_size;
          int rt = rtemp.size(); 
          for(int i =0;i<rt;i++){
            std::pair<const GeomDetUnit*, std::pair<LocalPoint, const GeomDetUnit* > > pxd = rtemp[i];
            const GeomDetUnit * geomdet = pxd.first;
            if (geomdet==mod->detUnit)rdigi.push_back(pxd.second);
          }
        } 
      }
      }
    }
      }
    norechits=true;
    edm::Handle<SiPixelRecHitCollection> rechitspixel;
if(pixelrechitProducerName!="ignore"){
std::vector< edm::Handle<SiPixelRecHitCollection> > pixelrechit_collections;
    ModuleLabelSelector sel(pixelrechitProducerName);
    event.getMany (sel, pixelrechit_collections);
      if (! pixelrechit_collections.empty ())norechits=false;
      if(norechits) cout << "no pixelrechits type "<<pixelrechitProducerName<<" in event" << endl;
      if(!norechits){
            std::vector< edm::Handle<SiPixelRecHitCollection> >::iterator i;
            i = pixelrechit_collections.begin ();
            const SiPixelRecHitCollection&rechitspixel = *(*i);
    if(rechitspixel.size()>0){
      std::vector<std::pair<const GeomDetUnit *,std::pair<LocalPoint, const GeomDetUnit* > > > rtemp;
      const std::vector<DetId> detIDs3 = rechitspixel.ids();
      for ( std::vector<DetId>::const_iterator detunit_iterator = detIDs3.begin(); detunit_iterator != detIDs3.end(); detunit_iterator++ )
        {
          unsigned int id = (*detunit_iterator).rawId();
          if(id!=999999999){ //if is valid detector
            DetId detIdObject(id);
            const GeomDetUnit      * geoUnit = geom->idToDetUnit( detIdObject );
            VisCuTkModule * mod = VisCuTkModuleMap::moduleMap[geoUnit];
            SiPixelRecHitCollection::range rechitRange = rechitspixel.get((*detunit_iterator));
            SiPixelRecHitCollection::const_iterator rechitRangeIteratorBegin = rechitRange.first;
            SiPixelRecHitCollection::const_iterator rechitRangeIteratorEnd   = rechitRange.second;
            SiPixelRecHitCollection::const_iterator iter=rechitRangeIteratorBegin;
            
            for(iter=rechitRangeIteratorBegin;iter!=rechitRangeIteratorEnd;++iter){//loop on the rechit
              LocalPoint position=iter->localPosition();
              std::pair<const GeomDetUnit*, std::pair<LocalPoint, const GeomDetUnit* > > pxd;
                pxd.first = geoUnit;
                pxd.second.first = position;
                pxd.second.second = geoUnit;
              rtemp.push_back(pxd);
              if(mod){ mod->rechit_size++; } else cout << "module not found " << endl; 
            }
          }
        } 
      cout << rtemp.size() << " pixel rechits  " << endl; 
      std::map<const GeomDetUnit * , VisCuTkModule *>::iterator imod;
      for (imod=VisCuTkModuleMap::moduleMap.begin();imod != VisCuTkModuleMap::moduleMap.end(); imod++){
        VisCuTkModule* mod = imod->second;
        unsigned int detType = mod->detUnit->geographicalId().subdetId();
        if(detType<3 && mod->rechit_size>0){
          mod->rechit_start=rd_pointer; rd_pointer = rd_pointer + mod->rechit_size;
          int rt = rtemp.size(); 
          for(int i =0;i<rt;i++){
            std::pair<const GeomDetUnit*, std::pair<LocalPoint, const GeomDetUnit* > > pxd = rtemp[i];
            const GeomDetUnit * geomdet = pxd.first;
            if (geomdet==mod->detUnit)rdigi.push_back(pxd.second);
          }
        } 
      }
      }
    }
    }
}

void VisCuTrackerRecHit::setMap

(

VisCuTkMap2D *

m

)

[inline]

Definition at line 33 of file VisCuTrackerRecHit.h.

Referenced by VisCuTkMapWindow::VisCuTkMapWindow().

{map=m;};

void VisCuTrackerRecHit::setTracker

(

VisCuCmsTracker *

tr

)

[inline]

Definition at line 31 of file VisCuTrackerRecHit.h.

References tracker.

Referenced by VisCuTkGeometryTwig::onNewEvent().

{tracker=tr;};

SoSeparator * VisCuTrackerRecHit::update

(

VisTrackingGeometry::DetectorRep

type

)

Definition at line 280 of file VisCuTrackerRecHit.cc.

References begin, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), VisCuTkModule::isVisible(), mod(), VisCuTkModuleMap::moduleMap, rdigi, VisCuTkModule::rechit_size, VisCuTkModule::rechit_start, VisTrackingGeometry::RZDet, edm::second(), funct::sqrt(), x, y, and z.

Referenced by VisCuTkGeometryTwig::update().

{

   SoSeparator *sep = new SoSeparator;
  if (! rdigi.empty ())
    { 



      SoMaterial *mat = new SoMaterial;
      mat->diffuseColor.setValue (1.0, 0.0, 0.0);//red
      mat->emissiveColor.setValue(1.0, 0.0, 0.0);
      
      SoDrawStyle *drawStyle = new SoDrawStyle;
      drawStyle->pointSize = 4.0;
      
       // global
      sep->addChild (mat);
      sep->addChild (drawStyle);
      
      try 
        {
          SoVertexProperty *vertices = new SoVertexProperty;
          int nVrtx = 0;
          std::map<const GeomDetUnit * , VisCuTkModule *>::iterator imod;
          for (imod=VisCuTkModuleMap::moduleMap.begin();imod != VisCuTkModuleMap::moduleMap.end(); imod++){
            VisCuTkModule* mod = imod->second;
            if(mod->rechit_size>0){
              if(mod->isVisible()){
                for(int idigi=mod->rechit_start;idigi!=(mod->rechit_start+mod->rechit_size);idigi++){
                  LocalPoint point = rdigi[idigi].first;
                  float x = rdigi[idigi].second->surface().toGlobal(point).x()/100;
                  float y;
                  if(type==RZDet){
                    x=0.0;
                   y = sqrt(rdigi[idigi].second->surface().toGlobal(point).x()*rdigi[idigi].second->surface().toGlobal(point).x()+rdigi[idigi].second->surface().toGlobal(point).y()*rdigi[idigi].second->surface().toGlobal(point).y()) / 100.0; if(rdigi[idigi].second->surface().toGlobal(point).y()< 0. )y = - y;} else   y = rdigi[idigi].second->surface().toGlobal(point).y()/100;
                  float z = rdigi[idigi].second->surface().toGlobal(point).z()/100;
                  vertices->vertex.set1Value (nVrtx++, SbVec3f (x, y, z));
                }
              }//if mod is visible
            }// if rechits in module
          }//loop on modules
          
          vertices->vertex.setNum (nVrtx);
          SoPointSet *points = new SoPointSet;
          points->vertexProperty.setValue (vertices);
          points->numPoints.setValue (nVrtx);
          
          sep->addChild (points);
        }
      catch (...) 
        {
          std::cout << "No Rechits" << std::endl;
        }
    }
      return sep;
 
}

std::string VisCuTrackerRecHit::updatetext

(

VisCuTkModule *

mod

)

Definition at line 518 of file VisCuTrackerRecHit.cc.

References rdigi, VisCuTkModule::rechit_size, VisCuTkModule::rechit_start, x, y, and z.

Referenced by VisCuTkGeometryTwig::update().

{
  std::ostringstream  text;
  if(mod->rechit_size>0){
    int nrechit=0;
    text << mod->rechit_size << " rechit : --------------------------------------- <p>" ;
    for(int idigi=mod->rechit_start;idigi!=(mod->rechit_start+mod->rechit_size);idigi++){
      LocalPoint point = rdigi[idigi].first;
       float x = rdigi[idigi].second->surface().toGlobal(point).x()/100;
       float y = rdigi[idigi].second->surface().toGlobal(point).y()/100;
       float z = rdigi[idigi].second->surface().toGlobal(point).z()/100;
      text << " Number=" << nrechit++ 
           << " Position=" << x  << ", " << y  << ", " << z << "<p>";
    }
    text <<  " -----------------------------" ;
  }
  return text.str();
}

---

Member Data Documentation

bool VisCuTrackerRecHit::firsteventloaded [private]

Definition at line 39 of file VisCuTrackerRecHit.h.

Referenced by drawrechit(), setEvent(), and VisCuTrackerRecHit().

edm::ESHandle<TrackerGeometry> VisCuTrackerRecHit::geom [private]

Definition at line 41 of file VisCuTrackerRecHit.h.

Referenced by setEvent().

bool VisCuTrackerRecHit::horizontal_view [private]

Definition at line 39 of file VisCuTrackerRecHit.h.

Referenced by drawrechit(), and VisCuTrackerRecHit().

VisCuTkMap2D* VisCuTrackerRecHit::map [private]

Definition at line 38 of file VisCuTrackerRecHit.h.

std::string VisCuTrackerRecHit::pixelrechitProducerName [private]

Definition at line 43 of file VisCuTrackerRecHit.h.

Referenced by setEvent(), and VisCuTrackerRecHit().

std::vector<std::pair<LocalPoint, const GeomDetUnit* > > VisCuTrackerRecHit::rdigi [private]

Definition at line 40 of file VisCuTrackerRecHit.h.

Referenced by drawrechit(), setEvent(), update(), and updatetext().

std::string VisCuTrackerRecHit::striprechitProducerName [private]

Definition at line 42 of file VisCuTrackerRecHit.h.

Referenced by setEvent(), and VisCuTrackerRecHit().

VisCuCmsTracker* VisCuTrackerRecHit::tracker [private]

Definition at line 37 of file VisCuTrackerRecHit.h.

Referenced by setTracker().

---

The documentation for this class was generated from the following files:

• VisReco/VisCustomTracker/interface/VisCuTrackerRecHit.h
• VisReco/VisCustomTracker/src/VisCuTrackerRecHit.cc

---