|
|
|
#include <DisplayManager.h>
Public Member Functions | |
| void | display (int ientry) |
| void | displayAll (bool noRedraw=true) |
| void | displayEvent (int run, int lumi, int event) |
| DisplayManager (PFRootEventManager *em, const char *optfile) | |
| void | displayNext () |
| void | displayNextInteresting (int ientry) |
| void | displayPFBlock (int blockNb) |
| void | displayPrevious () |
| void | drawWithNewGraphicAttributes () |
| void | enableDrawBrem (bool state) |
| void | enableDrawPFBlock (bool state) |
| void | findAndDraw (int ident) |
| void | findBlock (int ident) |
| void | lookForGenParticle (unsigned barcode) |
| look for particle with index i in MC truth. | |
| void | lookForMaxRecHit (bool ecal) |
| look for rechit with max energy in ecal or hcal. | |
| void | printDisplay (const char *directory="") const |
| void | printGenParticleInfo (std::string name, int barcode, int barcodeMother) |
| void | readOptions (const char *file) |
| void | reset () |
| void | setNewAttrToSimParticles () |
| void | unZoom () |
| void | updateDisplay () |
| virtual | ~DisplayManager () |
Public Attributes | |
| std::multimap< int, std::pair < int, int > > | blockIdentsMap_ |
| std::vector< float > | bremAttributes_ |
| TAttLine * | bremPatternL_ |
| TAttMarker * | bremPatternM_ |
| double | clusEnMin_ |
| TAttMarker * | clusPattern_ |
| TAttMarker * | clusPatternecal_ |
| TAttMarker * | clusPatternhfem_ |
| TAttMarker * | clusPatternho_ |
| TAttMarker * | clusPSPattern_ |
| std::vector< float > | clusterAttributes_ |
| bool | drawBrems_ |
| bool | drawClus_ |
| bool | drawClusterL_ |
| bool | drawGenParticles_ |
| bool | drawGsfTracks_ |
| bool | drawHits_ |
| bool | drawHO_ |
| bool | drawParticles_ |
| bool | drawPFBlocks_ |
| bool | drawTracks_ |
| double | genParticlePtMin_ |
| TAttMarker * | genPartPattern_ |
| std::multimap< int, GPFBase * > | graphicMap_ |
| std::vector< float > | gsfAttributes_ |
| TAttLine * | gsfPatternL_ |
| TAttMarker * | gsfPatternM_ |
| double | gsfPtMin_ |
| double | hitEnMin_ |
| double | particlePtMin_ |
| TAttMarker * | simPartPatternDefault_ |
| TAttMarker * | simPartPatternElec_ |
| TAttMarker * | simPartPatternK_ |
| TAttLine * | simPartPatternL_ |
| std::vector< TAttMarker * > | simPartPatternM_ |
| TAttMarker * | simPartPatternMuon_ |
| TAttMarker * | simPartPatternNeutron_ |
| TAttMarker * | simPartPatternPhoton_ |
| TAttMarker * | simPartPatternPi_ |
| TAttMarker * | simPartPatternProton_ |
| std::vector< float > | trackAttributes_ |
| TAttLine * | trackPatternL_ |
| TAttMarker * | trackPatternM_ |
| double | trackPtMin_ |
Private Member Functions | |
| void | createCanvas () |
| void | createGCluster (const reco::PFCluster &cluster, int ident, double phi0=0.) |
| void | createGGenParticle (HepMC::GenParticle *p) |
| void | createGPart (const reco::PFSimParticle &ptc, const std::vector< reco::PFTrajectoryPoint > &points, int ident, double pt, double phi0, double sign, bool displayInitial, int markerIndex) |
| void | createGRecHit (reco::PFRecHit &rh, int ident, double maxe, double phi0=0., int color=4) |
| void | createGTrack (reco::PFRecTrack &tr, const std::vector< reco::PFTrajectoryPoint > &points, int ident, double pt, double phi0, double sign, bool displayInitial, int linestyle, int kfgsfbrem=0) |
| void | displayCanvas () |
| void | drawGObject (int ident, int color, bool toInitialColor) |
| bool | findBadBremsId (int ident) |
| void | getDisplayOptions () |
| double | getMaxE (int layer) const |
| double | getMaxEEcal () |
| double | getMaxEHcal () |
| double | getMaxEHo () |
| void | loadGClusters () |
| void | loadGGenParticles () |
| void | loadGGsfRecTracks () |
| void | loadGPFBlocks () |
| void | loadGraphicObjects () |
| void | loadGRecHits () |
| void | loadGRecTracks () |
| void | loadGSimParticles () |
| void | retrieveBadBrems () |
| void | rubOutGPFBlock () |
Private Attributes | |
| std::vector< int > | badBremsId_ |
| std::vector< TH2F * > | displayHist_ |
| support histogram for x/y or r/z display. | |
| std::vector< TCanvas * > | displayView_ |
| vector of canvas for x/y or r/z display | |
| PFRootEventManager * | em_ |
| int | eventNumber_ |
| TBox | frontFaceECALRZ_ |
| ECAL in RZ view. | |
| TEllipse | frontFaceECALXY_ |
| ECAL in XY view. | |
| TEllipse | frontFaceHCALXY_ |
| HCAL in XY view. | |
| TEllipse | frontFaceHOXY_ |
| HO in XY view. | |
| bool | isGraphicLoaded_ |
| double | maxERecHitEcal_ |
| double | maxERecHitHcal_ |
| double | maxERecHitHo_ |
| int | maxEvents_ |
| IO * | options_ |
| options file parser | |
| std::vector< int > | selectedGObj_ |
| graphic object containers | |
| int | shiftId_ |
| std::vector< int > | viewSize_ |
| display pad xy size for (x,y) or (r,z) display | |
| std::vector< int > | viewSizeEtaPhi_ |
| display pad xy size for eta/phi view | |
| double | zoomFactor_ |
Definition at line 34 of file DisplayManager.h.
| DisplayManager::DisplayManager | ( | PFRootEventManager * | em, |
| const char * | optfile | ||
| ) |
Definition at line 41 of file DisplayManager.cc.
References createCanvas(), em_, PFRootEventManager::ev_, PFRootEventManager::eventNumber(), eventNumber_, maxEvents_, readOptions(), and fwlite::ChainEvent::size().
: em_(em), options_(0), maxERecHitEcal_(-1), maxERecHitHcal_(-1), maxERecHitHo_(-1), isGraphicLoaded_(false), shiftId_(SHIFTID) { readOptions( optfile ); eventNumber_ = em_->eventNumber(); //TODOLIST: re_initialize if new option file- better in em maxEvents_= em_->ev_->size(); createCanvas(); }
| DisplayManager::~DisplayManager | ( | ) | [virtual] |
| void DisplayManager::createCanvas | ( | ) | [private] |
Definition at line 278 of file DisplayManager.cc.
References gather_cfg::cout, displayHist_, displayView_, drawHO_, PFGeometry::ECALBarrel, PFGeometry::ECALEndcap, EHO, EPE, EPH, frontFaceECALRZ_, frontFaceECALXY_, frontFaceHCALXY_, frontFaceHOXY_, PFGeometry::HCALBarrel, PFGeometry::HOBarrel, PFGeometry::innerRadius(), PFGeometry::innerZ(), NViews, PFGeometry::outerRadius(), RZ, viewSize_, and XY.
Referenced by DisplayManager().
{
//TODOLIST: better TCanvas *displayView_[4]
displayView_.resize(NViews);
displayHist_.resize(NViews);
// TODOLIST:Canvases size
// Add menu to mofify canvas size
// Add menu of views to be drawn
displayView_[XY] = new TCanvas("displayXY_", "XY view",viewSize_[0], viewSize_[1]);
displayView_[RZ] = new TCanvas("displayRZ_", "RZ view",viewSize_[0], viewSize_[1]);
displayView_[EPE] = new TCanvas("displayEPE_", "eta/phi view, ECAL",viewSize_[0], viewSize_[1]);
displayView_[EPH] = new TCanvas("displayEPH_", "eta/phi view, HCAL",viewSize_[0], viewSize_[1]);
if (drawHO_) displayView_[EHO] = new TCanvas("displayEHO_", "eta/phi view, HO",viewSize_[0], viewSize_[1]);
for (int viewType=0;viewType<NViews;++viewType) {
if (!drawHO_ && viewType==EHO) continue;
displayView_[viewType]->SetGrid(0, 0);
displayView_[viewType]->SetBottomMargin(0.14);
displayView_[viewType]->SetLeftMargin(0.15);
displayView_[viewType]->SetRightMargin(0.05);
displayView_[viewType]->ToggleToolBar();
}
// Draw support histogram
double zLow = -650.;
double zUp = +650.;
double rLow = -450.;
double rUp = +450.;
if (!drawHO_) {zLow = -500.; zUp = +500.; rLow = -300.; rUp = +300.;}
displayHist_[XY] = new TH2F("hdisplayHist_XY", "", 500, rLow, rUp,
500, rLow, rUp);
displayHist_[XY]->SetXTitle("X [cm]");
displayHist_[XY]->SetYTitle("Y [cm]");
displayHist_[RZ] = new TH2F("hdisplayHist_RZ", "",500, zLow, zUp,
500, rLow, rUp);
displayHist_[RZ]->SetXTitle("Z [cm]");
displayHist_[RZ]->SetYTitle("R [cm]");
displayHist_[EPE] = new TH2F("hdisplayHist_EP", "", 500, -5, 5,
500, -3.5, 3.5);
displayHist_[EPE]->SetXTitle("#eta");
displayHist_[EPE]->SetYTitle("#phi [rad]");
displayHist_[EPH] = displayHist_[EPE];
// displayHist_[EHO] = new TH2F("hdisplayHist_HO", "", 150, -1.5, 1.5,
// 500, -3.5, 3.5);
// displayHist_[EHO]->SetXTitle("#eta");
// displayHist_[EHO]->SetYTitle("#phi [rad]");
if (drawHO_) displayHist_[EHO] = displayHist_[EPE];
for (int viewType=0;viewType<NViews;++viewType){
if (!drawHO_ && viewType==EHO) continue;
displayHist_[viewType]->SetStats(kFALSE);
displayHist_[viewType]->GetYaxis()->SetTitleSize(0.06);
displayHist_[viewType]->GetYaxis()->SetTitleOffset(1.2);
displayHist_[viewType]->GetXaxis()->SetTitleSize(0.06);
displayHist_[viewType]->GetYaxis()->SetLabelSize(0.045);
displayHist_[viewType]->GetXaxis()->SetLabelSize(0.045);
}
// Draw ECAL front face
frontFaceECALXY_.SetX1(0);
frontFaceECALXY_.SetY1(0);
frontFaceECALXY_.SetR1(PFGeometry::innerRadius(PFGeometry::ECALBarrel));
frontFaceECALXY_.SetR2(PFGeometry::innerRadius(PFGeometry::ECALBarrel));
frontFaceECALXY_.SetFillStyle(0);
// Draw HCAL front face
frontFaceHCALXY_.SetX1(0);
frontFaceHCALXY_.SetY1(0);
frontFaceHCALXY_.SetR1(PFGeometry::innerRadius(PFGeometry::HCALBarrel));
frontFaceHCALXY_.SetR2(PFGeometry::innerRadius(PFGeometry::HCALBarrel));
frontFaceHCALXY_.SetFillStyle(0);
if (drawHO_) {
// Draw HO front face
frontFaceHOXY_.SetX1(0);
frontFaceHOXY_.SetY1(0);
frontFaceHOXY_.SetR1(PFGeometry::innerRadius(PFGeometry::HOBarrel));
frontFaceHOXY_.SetR2(PFGeometry::outerRadius(PFGeometry::HOBarrel));
frontFaceHOXY_.SetFillStyle(0);
}
// Draw ECAL side
frontFaceECALRZ_.SetX1(-1.*PFGeometry::innerZ(PFGeometry::ECALEndcap));
frontFaceECALRZ_.SetY1(-1.*PFGeometry::innerRadius(PFGeometry::ECALBarrel));
frontFaceECALRZ_.SetX2(PFGeometry::innerZ(PFGeometry::ECALEndcap));
frontFaceECALRZ_.SetY2(PFGeometry::innerRadius(PFGeometry::ECALBarrel));
frontFaceECALRZ_.SetFillStyle(0);
cout <<"End of DisplayManager::createCanvas()"<<endl;
}
| void DisplayManager::createGCluster | ( | const reco::PFCluster & | cluster, |
| int | ident, | ||
| double | phi0 = 0. |
||
| ) | [private] |
Definition at line 378 of file DisplayManager.cc.
References clusPattern_, clusPatternecal_, clusPatternhfem_, clusPatternho_, clusPSPattern_, funct::cos(), drawHO_, PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, EHO, EPE, EPH, eta(), graphicMap_, PFLayer::HCAL_BARREL2, PFLayer::HF_EM, reco::PFCluster::layer(), NViews, phi, reco::CaloCluster::position(), PFLayer::PS1, PFLayer::PS2, RZ, XY, and reco::CaloCluster::z().
Referenced by loadGClusters().
{
double eta = cluster.position().Eta();
double phi = cluster.position().Phi();
// int type = cluster.type();
// if(algosToDisplay_.find(type) == algosToDisplay_.end() )
// return;
// TCutG* cutg = (TCutG*) gROOT->FindObject("CUTG");
// if( cutg && !cutg->IsInside( eta, phi ) ) return;
//int color = clusterAttributes_[2];
//if ( cluster.layer()==PFLayer::PS1 || cluster.layer()==PFLayer::PS2 )
// color = clusterAttributes_[3];
//int markerSize = clusterAttributes_[1];
//int markerStyle = clusterAttributes_[0];
int clusType=0;
if ( cluster.layer()==PFLayer::PS1 || cluster.layer()==PFLayer::PS2 )
clusType=1;
if ( cluster.layer()==PFLayer::ECAL_BARREL || cluster.layer()==PFLayer::ECAL_ENDCAP)
clusType=2;
if ( cluster.layer()==PFLayer::HF_EM) clusType=3;
if ( cluster.layer()==PFLayer::HCAL_BARREL2) clusType=4;
const math::XYZPoint& xyzPos = cluster.position();
GPFCluster *gc;
for (int viewType=0;viewType<NViews;viewType++){
if (!drawHO_ && viewType==EHO) continue;
switch(viewType) {
case XY:
{
if (clusType==0) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.X(), xyzPos.Y(), clusPattern_);
} else if (clusType==1) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.X(), xyzPos.Y(), clusPSPattern_);
} else if (clusType==2) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.X(), xyzPos.Y(), clusPatternecal_);
} else if (clusType==3) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.X(), xyzPos.Y(), clusPatternhfem_);
} else {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.X(), xyzPos.Y(), clusPatternho_);
}
graphicMap_.insert(pair<int,GPFBase *> (ident, gc));
}
break;
case RZ:
{
double sign = 1.;
if (cos(phi0 - phi) < 0.)
sign = -1.;
if ( clusType==0) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.z(),sign*xyzPos.Rho(),
clusPattern_);
} else if (clusType==1) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.z(),sign*xyzPos.Rho(),
clusPSPattern_);
} else if (clusType==2) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.z(),sign*xyzPos.Rho(),
clusPatternecal_);
} else if (clusType==3) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.z(),sign*xyzPos.Rho(),
clusPatternhfem_);
} else {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
xyzPos.z(),sign*xyzPos.Rho(),
clusPatternho_);
}
graphicMap_.insert(pair<int,GPFBase *> (ident, gc));
}
break;
case EPE:
{
if( cluster.layer()<0 || cluster.layer()==PFLayer::HF_EM) {
if (clusType==2) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
eta,phi,
clusPatternecal_);
}
else if (clusType==3) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
eta,phi,
clusPatternhfem_);
}
else {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
eta,phi,
clusPSPattern_);
}
graphicMap_.insert(pair<int,GPFBase *> (ident, gc));
}
}
break;
case EPH:
{
if( cluster.layer()>0 && cluster.layer()!=PFLayer::HF_EM && cluster.layer()!=PFLayer::HCAL_BARREL2) {
if (clusType==0) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
eta,phi,clusPattern_);
}
else {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
eta,phi,clusPSPattern_);
}
graphicMap_.insert(pair<int,GPFBase *> (ident, gc));
}
}
break;
case EHO:
{
if( cluster.layer()>0 && cluster.layer()==PFLayer::HCAL_BARREL2) {
if (clusType==4) {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
eta,phi,clusPatternho_);
}
else {
gc = new GPFCluster(this,
viewType,ident,
&cluster,
eta,phi,clusPSPattern_);
}
graphicMap_.insert(pair<int,GPFBase *> (ident, gc));
}
}
break;
default :break;
}
}
}
| void DisplayManager::createGGenParticle | ( | HepMC::GenParticle * | p | ) | [private] |
Definition at line 1517 of file DisplayManager.cc.
References drawHO_, alignCSCRings::e, em_, eta(), configurableAnalysis::GenParticle, GENPARTICLEID, genPartPattern_, PFRootEventManager::getGenParticleName(), graphicMap_, mergeVDriftHistosByStation::name, NViews, phi, shiftId_, x, and detailsBasic3DVector::y.
Referenced by loadGGenParticles().
{
// these are the beam protons
if ( !p->production_vertex() && p->pdg_id() == 2212 ) return;
int partId = p->pdg_id();
std::string name;
std::string latexStringName;
name = em_->getGenParticleName(partId,latexStringName);
int barcode = p->barcode();
int genPartId=(GENPARTICLEID<<shiftId_) | barcode;
// int vertexId1 = 0;
// vertexId1 = p->production_vertex()->barcode();
// math::XYZVector vertex1 (p->production_vertex()->position().x(),
// p->production_vertex()->position().y(),
// p->production_vertex()->position().z());
math::XYZTLorentzVector momentum1(p->momentum().px(),
p->momentum().py(),
p->momentum().pz(),
p->momentum().e());
double eta = momentum1.eta();
if ( eta > +10. ) eta = +10.;
if ( eta < -10. ) eta = -10.;
double phi = momentum1.phi();
double pt = momentum1.pt();
double e = momentum1.e();
//mother ?
// Colin: the following line gives a segmentation fault when there is
// no particles entering the production vertex.
// const HepMC::GenParticle* mother =
// *(p->production_vertex()->particles_in_const_begin());
// protecting against this in the following way:
const HepMC::GenParticle* mother = 0;
if(p->production_vertex() &&
p->production_vertex()->particles_in_size() ) {
mother =
*(p->production_vertex()->particles_in_const_begin());
}
// Colin: no need to declare this pointer in this context.
// the declaration can be more local.
// GPFGenParticle *gp;
int modnViews = (drawHO_) ? NViews : NViews-1;
if ( mother ) {
int barcodeMother = mother->barcode();
math::XYZTLorentzVector momentumMother(mother->momentum().px(),
mother->momentum().py(),
mother->momentum().pz(),
mother->momentum().e());
double etaMother = momentumMother.eta();
if ( etaMother > +10. ) etaMother = +10.;
if ( etaMother < -10. ) etaMother = -10.;
double phiMother = momentumMother.phi();
double x[2],y[2];
x[0]=etaMother;x[1]=eta;
y[0]=phiMother;y[1]=phi;
for (int view = 2; view< modnViews; view++) {
GPFGenParticle* gp = new GPFGenParticle(this,
view, genPartId,
x, y, //double *, double *
e,pt,barcode,barcodeMother,
genPartPattern_,
name,latexStringName);
graphicMap_.insert(pair<int,GPFBase *> (genPartId, gp));
}
}
else { //no Mother
for (int view = 2; view< modnViews; view++) {
GPFGenParticle* gp = new GPFGenParticle(this,
view, genPartId,
eta, phi, //double double
e,pt,barcode,
genPartPattern_,
name, latexStringName);
graphicMap_.insert(pair<int,GPFBase *> (genPartId, gp));
}
}
}
| void DisplayManager::createGPart | ( | const reco::PFSimParticle & | ptc, |
| const std::vector< reco::PFTrajectoryPoint > & | points, | ||
| int | ident, | ||
| double | pt, | ||
| double | phi0, | ||
| double | sign, | ||
| bool | displayInitial, | ||
| int | markerIndex | ||
| ) | [private] |
no point inside graphical cut.
Definition at line 573 of file DisplayManager.cc.
References reco::PFTrajectoryPoint::ClosestApproach, gather_cfg::cout, drawHO_, EHO, EPE, EPH, eta(), graphicMap_, i, NViews, phi, RZ, simPartPatternL_, simPartPatternM_, xPos, XY, and yPos.
Referenced by loadGSimParticles().
{
//bool inside = false;
//TCutG* cutg = (TCutG*) gROOT->FindObject("CUTG");
bool debug_createGpart = false;
// bool debug_createGpart = true;
for (int viewType=0;viewType<NViews;++viewType) {
if (!drawHO_ && viewType==EHO) continue;
// for (int viewType=0;viewType<4;++viewType) {
// reserving space. nb not all trajectory points are valid
vector<double> xPos;
xPos.reserve( points.size() );
vector<double> yPos;
yPos.reserve( points.size() );
for(unsigned i=0; i<points.size(); i++) {
if( !points[i].isValid() ) continue;
const math::XYZPoint& xyzPos = points[i].position();
double eta = xyzPos.Eta();
double phi = xyzPos.Phi();
// always take momentum for the first point
// otherwise decay products from resonances will have wrong eta, phi
//calculated from position (0,0,0) (MDN 11 June 2008)
if( points[i].layer() == reco::PFTrajectoryPoint::ClosestApproach ) {
// if( !displayInitial &&
// points[i].layer() == reco::PFTrajectoryPoint::ClosestApproach ) {
const math::XYZTLorentzVector& mom = points[i].momentum();
eta = mom.Eta();
phi = mom.Phi();
}
//if( !cutg || cutg->IsInside( eta, phi ) )
// inside = true;
switch(viewType) {
case XY:
xPos.push_back(xyzPos.X());
yPos.push_back(xyzPos.Y());
break;
case RZ:
xPos.push_back(xyzPos.Z());
yPos.push_back(sign*xyzPos.Rho());
break;
case EPE:
case EPH:
case EHO:
xPos.push_back( eta );
yPos.push_back( phi );
break;
default:break;
}
}
if (viewType == EPE && debug_createGpart) {
cout << " display PFsim eta/phi view ECAL" << endl;
cout << " nb of points for display " << xPos.size() << endl;
for(unsigned i=0; i<xPos.size(); i++) {
cout << " point " << i << " x/y " << xPos[i] <<"/" << yPos[i]<< endl;
}
}
//if( !inside ) return;
//int color = 4;
GPFSimParticle *gc = new GPFSimParticle(this,
viewType, ident,
&ptc,
xPos.size(),&xPos[0],&yPos[0],
pt,
simPartPatternM_[markerIndex],
simPartPatternL_,
"pl");
graphicMap_.insert(pair<int,GPFBase *> (ident, gc));
//graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFSimParticle(this,viewType,ident,&ptc,xPos.size(),&xPos[0],&yPos[0],pt,simPartPattern_[indexMarker], "pl")));
}
}
| void DisplayManager::createGRecHit | ( | reco::PFRecHit & | rh, |
| int | ident, | ||
| double | maxe, | ||
| double | phi0 = 0., |
||
| int | color = 4 |
||
| ) | [private] |
Definition at line 658 of file DisplayManager.cc.
References dtNoiseDBValidation_cfg::cerr, PFRootEventManager::clusterAlgoECAL_, PFRootEventManager::clusterAlgoHCAL_, PFRootEventManager::clusterAlgoHFEM_, PFRootEventManager::clusterAlgoHFHAD_, PFRootEventManager::clusterAlgoPS_, funct::cos(), drawHO_, PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, EHO, em_, reco::PFRecHit::energy(), EPE, EPH, eta(), reco::PFRecHit::getCornersXYZ(), graphicMap_, PFLayer::HCAL_BARREL1, PFLayer::HCAL_BARREL2, PFLayer::HCAL_ENDCAP, PFLayer::HF_EM, PFLayer::HF_HAD, reco::PFRecHit::layer(), create_public_lumi_plots::log, NViews, phi, Pi, reco::PFRecHit::position(), PFLayer::PS1, PFLayer::PS2, alignCSCRings::r, RZ, GOODCOLL_filter_cfg::thresh, PFClusterAlgo::threshBarrel(), PFClusterAlgo::threshEndcap(), x, XY, detailsBasic3DVector::y, and z.
Referenced by loadGRecHits().
{
double me = maxe;
double thresh = 0;
int layer = rh.layer();
switch(layer) {
case PFLayer::ECAL_BARREL:
thresh = em_->clusterAlgoECAL_.threshBarrel();
break;
case PFLayer::ECAL_ENDCAP:
thresh = em_->clusterAlgoECAL_.threshEndcap();
break;
case PFLayer::HCAL_BARREL1:
case PFLayer::HCAL_BARREL2:
thresh = em_->clusterAlgoHCAL_.threshBarrel();
break;
case PFLayer::HCAL_ENDCAP:
thresh = em_->clusterAlgoHCAL_.threshEndcap();
break;
case PFLayer::HF_HAD:
// how to handle a different threshold for HF?
thresh = em_->clusterAlgoHFHAD_.threshEndcap();
break;
case PFLayer::HF_EM:
// how to handle a different threshold for HF?
thresh = em_->clusterAlgoHFEM_.threshEndcap();
break;
case PFLayer::PS1:
case PFLayer::PS2:
me = -1;
thresh = em_->clusterAlgoPS_.threshBarrel();
break;
default:
{
cerr<<"DisplayManager::createGRecHit : manage other layers."
<<" GRechit notcreated."<<endl;
return;
}
}
if( rh.energy() < thresh ) return;
//loop on all views
for(int viewType=0;viewType<NViews;++viewType) {
if (!drawHO_ && viewType==EHO) continue;
bool isHCAL = (layer == PFLayer::HCAL_BARREL1 ||
// layer == PFLayer::HCAL_BARREL2 ||
layer == PFLayer::HCAL_ENDCAP ||
layer == PFLayer::HF_HAD);
if( viewType == EPH &&
! isHCAL) {
continue;
}
// on EPE view, draw only HCAL and preshower
if( viewType == EPE && (isHCAL || layer ==PFLayer::HCAL_BARREL2)) {
continue;
}
if (viewType == EHO && layer !=PFLayer::HCAL_BARREL2) continue;
double rheta = rh.position().Eta();
double rhphi = rh.position().Phi();
double sign = 1.;
if (cos(phi0 - rhphi) < 0.) sign = -1.;
double etaSize[4];
double phiSize[4];
double x[5];
double y[5];
double z[5];
double r[5];
double eta[5];
double phi[5];
double xprop[5];
double yprop[5];
double etaprop[5];
double phiprop[5];
const std::vector< math::XYZPoint >& corners = rh.getCornersXYZ();
assert(corners.size() == 4);
double propfact = 0.95; // so that the cells don't overlap ?
double ampl=0;
if(me>0) ampl = (log(rh.energy() + 1.)/log(me + 1.));
for ( unsigned jc=0; jc<4; ++jc ) {
phiSize[jc] = rhphi-corners[jc].Phi();
etaSize[jc] = rheta-corners[jc].Eta();
if ( phiSize[jc] > 1. ) phiSize[jc] -= 2.*TMath::Pi();
if ( phiSize[jc] < -1. ) phiSize[jc]+= 2.*TMath::Pi();
phiSize[jc] *= propfact;
etaSize[jc] *= propfact;
math::XYZPoint cornerposxyz = corners[jc];
x[jc] = cornerposxyz.X();
y[jc] = cornerposxyz.Y();
z[jc] = cornerposxyz.Z();
r[jc] = sign*cornerposxyz.Rho();
eta[jc] = rheta - etaSize[jc];
phi[jc] = rhphi - phiSize[jc];
// cell area is prop to log(E)
// not drawn for preshower.
// otherwise, drawn for eta/phi view, and for endcaps in xy view
if( layer != PFLayer::PS1 &&
layer != PFLayer::PS2 &&
( viewType == EPE ||
viewType == EPH ||
viewType == EHO ||
( viewType == XY &&
( layer == PFLayer::ECAL_ENDCAP ||
layer == PFLayer::HCAL_ENDCAP ||
layer == PFLayer::HF_HAD
) ) ) ) {
math::XYZPoint centreXYZrot = rh.position();
math::XYZPoint centertocorner(x[jc] - centreXYZrot.X(),
y[jc] - centreXYZrot.Y(),
0 );
math::XYZPoint centertocornerep(eta[jc] - centreXYZrot.Eta(),
phi[jc] - centreXYZrot.Phi(),
0 );
// centertocorner -= centreXYZrot;
xprop[jc] = centreXYZrot.X() + centertocorner.X()*ampl;
yprop[jc] = centreXYZrot.Y() + centertocorner.Y()*ampl;
etaprop[jc] = centreXYZrot.Eta() + centertocornerep.X()*ampl;
phiprop[jc] = centreXYZrot.Phi() + centertocornerep.Y()*ampl;
}
}//loop on jc
if(layer == PFLayer::ECAL_BARREL ||
layer == PFLayer::HCAL_BARREL1 ||
layer == PFLayer::HCAL_BARREL2 || viewType == RZ) {
// we are in the barrel. Determining which corners to shift
// away from the center to represent the cell energy
int i1 = -1;
int i2 = -1;
if(fabs(phiSize[1]-phiSize[0]) > 0.0001) {
if (viewType == XY) {
i1 = 2;
i2 = 3;
} else if (viewType == RZ) {
i1 = 1;
i2 = 2;
}
} else {
if (viewType == XY) {
i1 = 1;
i2 = 2;
} else if (viewType == RZ) {
i1 = 2;
i2 = 3;
}
}
double signx=1.;
if(layer == PFLayer::HCAL_BARREL2) { signx=-1.;}
x[i1] *= 1+signx*ampl/2.;
x[i2] *= 1+signx*ampl/2.;
y[i1] *= 1+signx*ampl/2.;
y[i2] *= 1+signx*ampl/2.;
z[i1] *= 1+signx*ampl/2.;
z[i2] *= 1+signx*ampl/2.;
r[i1] *= 1+signx*ampl/2.;
r[i2] *= 1+signx*ampl/2.;
}
x[4]=x[0];
y[4]=y[0]; // closing the polycell
z[4]=z[0];
r[4]=r[0]; // closing the polycell
eta[4]=eta[0];
phi[4]=phi[0]; // closing the polycell
int npoints=5;
switch( viewType ) {
case XY:
{
if(layer == PFLayer::ECAL_BARREL ||
layer == PFLayer::HCAL_BARREL1 ||
layer == PFLayer::HCAL_BARREL2) {
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,x,y,color,"f")));
} else {
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,x,y,color,"l")));
if( ampl>0 ) { // not for preshower
xprop[4]=xprop[0];
yprop[4]=yprop[0]; // closing the polycell
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,xprop,yprop,color,"f")));
}
}
}
break;
case RZ:
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,z,r,color,"f")));
break;
case EPE:
{
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,eta,phi,color,"l")));
if( ampl>0 ) { // not for preshower
etaprop[4]=etaprop[0];
phiprop[4]=phiprop[0]; // closing the polycell
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,etaprop,phiprop,color,"f")));
}
}
break;
case EPH:
{
if(layer == PFLayer::ECAL_BARREL ||
layer == PFLayer::HCAL_BARREL1) {
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,eta,phi,color,"l")));
if( ampl>0 ) { // not for preshower
etaprop[4]=etaprop[0];
phiprop[4]=phiprop[0]; // closing the polycell
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,etaprop,phiprop,color,"f")));
}
}
}
break;
case EHO:
{
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,eta,phi,color,"l")));
if(layer == PFLayer::HCAL_BARREL2) {
if( ampl>0 ) { // not for preshower
etaprop[4]=etaprop[0];
phiprop[4]=phiprop[0]; // closing the polycell
graphicMap_.insert(pair<int,GPFBase *> (ident,new GPFRecHit(this, viewType,ident,&rh,npoints,etaprop,phiprop,color,"f")));
}
}
}
break;
default: break;
}//switch end
} //loop on views
}
| void DisplayManager::createGTrack | ( | reco::PFRecTrack & | tr, |
| const std::vector< reco::PFTrajectoryPoint > & | points, | ||
| int | ident, | ||
| double | pt, | ||
| double | phi0, | ||
| double | sign, | ||
| bool | displayInitial, | ||
| int | linestyle, | ||
| int | kfgsfbrem = 0 |
||
| ) | [private] |
no point inside graphical cut.
Definition at line 926 of file DisplayManager.cc.
References bremPatternL_, bremPatternM_, reco::PFTrajectoryPoint::ClosestApproach, drawHO_, EHO, EPE, EPH, eta(), graphicMap_, gsfPatternL_, gsfPatternM_, gt, i, NViews, phi, RZ, trackPatternL_, trackPatternM_, xPos, XY, and yPos.
Referenced by loadGGsfRecTracks(), and loadGRecTracks().
{
// bool inside = false;
//TCutG* cutg = (TCutG*) gROOT->FindObject("CUTG");
for (int viewType=0;viewType<NViews;++viewType) {
if (!drawHO_ && viewType==EHO) continue;
// reserving space. nb not all trajectory points are valid
vector<double> xPos;
xPos.reserve( points.size() );
vector<double> yPos;
yPos.reserve( points.size() );
for(unsigned i=0; i<points.size(); i++) {
if( !points[i].isValid() ) continue;
const math::XYZPoint& xyzPos = points[i].position();
//muriel
//if(kfgsfbrem)
// std::cout << xyzPos << std::endl;
double eta = xyzPos.Eta();
double phi = xyzPos.Phi();
if( !displayInitial &&
points[i].layer() == reco::PFTrajectoryPoint::ClosestApproach ) {
const math::XYZTLorentzVector& mom = points[i].momentum();
eta = mom.Eta();
phi = mom.Phi();
}
//if( !cutg || cutg->IsInside( eta, phi ) )
// inside = true;
switch(viewType) {
case XY:
xPos.push_back(xyzPos.X());
yPos.push_back(xyzPos.Y());
break;
case RZ:
xPos.push_back(xyzPos.Z());
yPos.push_back(sign*xyzPos.Rho());
break;
case EPE:
case EPH:
case EHO:
xPos.push_back( eta );
yPos.push_back( phi );
break;
}
}
//if( !inside ) return;
//fill map with graphic objects
GPFTrack *gt=0;
if(kfgsfbrem==0) {
gt = new GPFTrack(this,
viewType,ident,
&tr,
xPos.size(),&xPos[0],&yPos[0],pt,
trackPatternM_,trackPatternL_,"pl");
}
else if (kfgsfbrem==1) {
// std::cout << " Creating the GSF track " << std::endl;
gt = new GPFTrack(this,
viewType,ident,
&tr,
xPos.size(),&xPos[0],&yPos[0],pt,
gsfPatternM_,gsfPatternL_,"pl");
}
else if (kfgsfbrem==2) {
// std::cout << " Creating the Brem " << std::endl;
// std::cout<<"adr brem dans create:"<<&tr<<std::flush<<std::endl;
gt = new GPFTrack(this,
viewType,ident,
&tr,
xPos.size(),&xPos[0],&yPos[0],pt,
bremPatternM_,bremPatternL_,"pl");
}
graphicMap_.insert(pair<int,GPFBase *> (ident, gt));
}
}
| void DisplayManager::display | ( | int | ientry | ) |
Definition at line 1026 of file DisplayManager.cc.
References dtNoiseDBValidation_cfg::cerr, displayAll(), em_, PFRootEventManager::eventNumber(), eventNumber_, isGraphicLoaded_, loadGraphicObjects(), maxEvents_, PFRootEventManager::processEntry(), and reset().
Referenced by displayNext(), displayPrevious(), and DialogFrame::doReProcessEvent().
{
if (ientry<0 || ientry>maxEvents_) {
std::cerr<<"DisplayManager::no event matching criteria"<<std::endl;
return;
}
reset();
em_->processEntry(ientry);
eventNumber_= em_->eventNumber();
loadGraphicObjects();
isGraphicLoaded_= true;
displayAll();
}
| void DisplayManager::displayAll | ( | bool | noRedraw = true | ) |
Definition at line 1042 of file DisplayManager.cc.
References BREMID, clusEnMin_, CLUSTERECALID, CLUSTERHCALID, CLUSTERHFEMID, CLUSTERHFHADID, CLUSTERHOID, CLUSTERIBID, CLUSTERPSID, gather_cfg::cout, displayCanvas(), displayView_, drawBrems_, drawClus_, drawGenParticles_, drawGsfTracks_, drawHits_, drawHO_, drawParticles_, drawTracks_, EHO, EPE, EPH, GENPARTICLEID, genParticlePtMin_, graphicMap_, gsfPtMin_, GSFRECTRACKID, hitEnMin_, i, isGraphicLoaded_, NViews, AlCaHLTBitMon_ParallelJobs::p, particlePtMin_, RECHITECALID, RECHITHCALID, RECHITHFEMID, RECHITHFHADID, RECHITHOID, RECHITPSID, RECTRACKID, shiftId_, SIMPARTICLEID, and trackPtMin_.
Referenced by DialogFrame::areBremVisible(), display(), displayEvent(), displayNextInteresting(), DialogFrame::doModifyOptions(), DialogFrame::doModifyPtThreshold(), drawWithNewGraphicAttributes(), and findAndDraw().
{
if (!isGraphicLoaded_) {
std::cout<<" no Graphic Objects to draw"<<std::endl;
return;
}
if (noRedraw) {
for (int viewType=0;viewType<NViews;++viewType) {
if (!drawHO_ && viewType==EHO) continue;
displayView_[viewType]->cd();
gPad->Clear();
}
//TODOLIST: add test on view to draw
displayCanvas();
}
std::multimap<int,GPFBase *>::iterator p;
for (p=graphicMap_.begin();p!=graphicMap_.end();p++) {
int ident=p->first;
int type=ident >> shiftId_;
int view = p->second->getView();
if (!drawHO_ && (view ==EHO || type==CLUSTERHOID || type==RECHITHOID)) continue;
switch (type) {
case CLUSTERECALID:
case CLUSTERHCALID:
case CLUSTERHOID:
case CLUSTERHFEMID:
case CLUSTERHFHADID:
case CLUSTERPSID:
case CLUSTERIBID:
{
// cout<<"displaying "<<type<<" "<<p->second->getEnergy()<<" "<<clusEnMin_<<" on view "<<view<<endl;
if (drawClus_)
if (p->second->getEnergy() > clusEnMin_) {
displayView_[view]->cd();
p->second->draw();
}
}
break;
case RECHITECALID:
case RECHITHCALID:
case RECHITHOID:
case RECHITHFEMID:
case RECHITHFHADID:
case RECHITPSID:
{
if (!noRedraw) break;
if (drawHits_)
if(p->second->getEnergy() > hitEnMin_) {
displayView_[view]->cd();
p->second->draw();
}
break;
}
case RECTRACKID:
{
if (drawTracks_)
if (p->second->getPt() > trackPtMin_) {
displayView_[view]->cd();
p->second->draw();
}
}
break;
case GSFRECTRACKID:
{
if (drawGsfTracks_)
if (p->second->getPt() > gsfPtMin_) {
displayView_[view]->cd();
p->second->draw();
}
}
break;
case BREMID:
{
if (drawBrems_)
{
displayView_[view]->cd();
p->second->draw();
}
}
break;
case SIMPARTICLEID:
{
if (drawParticles_)
if (p->second->getPt() > particlePtMin_) {
displayView_[view]->cd();
p->second->draw();
}
}
break;
case GENPARTICLEID:
{
if (drawGenParticles_)
if (p->second->getPt() > genParticlePtMin_)
if (view == EPH || view ==EPE || view ==EHO) {
displayView_[view]->cd();
p->second->draw();
}
}
break;
default : std::cout<<"DisplayManager::displayAll()-- unknown object "<<std::endl;
} //switch end
} //for end
for (int i=0;i<NViews;i++) {
if (!drawHO_ && i==EHO) continue;
displayView_[i]->cd();
gPad->Modified();
displayView_[i]->Update();
}
}
| void DisplayManager::displayCanvas | ( | ) | [private] |
Definition at line 1191 of file DisplayManager.cc.
References displayHist_, displayView_, drawHO_, EHO, eta(), jptDQMConfig_cff::etaMax, jptDQMConfig_cff::etaMin, frontFaceECALRZ_, frontFaceECALXY_, frontFaceHCALXY_, frontFaceHOXY_, prof2calltree::l, NViews, alignCSCRings::r, RZ, funct::tan(), and XY.
Referenced by displayAll().
{
double zLow = -500.;
double zUp = +500.;
double rUp = +400.;
if (!drawHO_) rUp=300.0;
double etarng=(drawHO_) ? 1.19 : 1.39;
double ratio =(drawHO_) ? 1.15 : 0.90;
double scal = (drawHO_) ? 1.00 : 0.85;
//TODOLIST : test wether view on/off
//if (!displayView_[viewType] || !gROOT->GetListOfCanvases()->FindObject(displayView_[viewType]) ) {
// assert(viewSize_.size() == 2);
for (int viewType=0;viewType<NViews;++viewType) {
if (!drawHO_ && viewType==EHO) continue;
displayView_[viewType]->cd();
displayHist_[viewType]->Draw();
switch(viewType) {
case XY:
frontFaceECALXY_.Draw();
frontFaceHCALXY_.Draw();
frontFaceHOXY_.Draw();
break;
case RZ:
{// Draw lines at different etas
TLine l;
l.SetLineColor(1);
l.SetLineStyle(3);
TLatex etaLeg;
etaLeg.SetTextSize( (drawHO_) ? 0.03 : 0.02);
float etaMin = -3.;
float etaMax = +3.;
float etaBin = 0.2;
int nEtas = int((etaMax - etaMin)/0.2) + 1;
for (int iEta = 0; iEta <= nEtas; iEta++) {
float eta = etaMin + iEta*etaBin;
float r = ratio*rUp;
TVector3 etaImpact;
etaImpact.SetPtEtaPhi(r, eta, 0.);
etaLeg.SetTextAlign(21);
if (eta <= -etarng) {
etaImpact.SetXYZ(0.,scal*zLow*tan(etaImpact.Theta()),scal*zLow);
etaLeg.SetTextAlign(31);
} else if (eta >= etarng) {
etaImpact.SetXYZ(0.,scal*zUp*tan(etaImpact.Theta()),scal*zUp);
etaLeg.SetTextAlign(11);
}
if (drawHO_) {
if (fabs(eta)<1.29) {
etaLeg.SetTextSize(0.030);
} else {
etaLeg.SetTextSize(0.022);
}
}
l.DrawLine(0., 0., etaImpact.Z(), etaImpact.Perp());
etaLeg.DrawLatex(1.0*etaImpact.Z(), etaImpact.Perp(), Form("%2.1f", eta));
}
frontFaceECALRZ_.Draw();
}
break;
default: break;
} //end switch
}
}
| void DisplayManager::displayEvent | ( | int | run, |
| int | lumi, | ||
| int | event | ||
| ) |
Definition at line 1015 of file DisplayManager.cc.
References displayAll(), em_, PFRootEventManager::eventNumber(), eventNumber_, isGraphicLoaded_, loadGraphicObjects(), PFRootEventManager::processEvent(), and reset().
{
reset();
em_->processEvent(run, lumi, event);
eventNumber_= em_->eventNumber();
loadGraphicObjects();
isGraphicLoaded_= true;
displayAll();
}
| void DisplayManager::displayNext | ( | ) |
Definition at line 1259 of file DisplayManager.cc.
References display(), em_, PFRootEventManager::eventNumber(), and eventNumber_.
Referenced by DialogFrame::doNextEvent().
{
int eventNumber_=em_->eventNumber();
display(++eventNumber_);
}
| void DisplayManager::displayNextInteresting | ( | int | ientry | ) |
Definition at line 1265 of file DisplayManager.cc.
References dtNoiseDBValidation_cfg::cerr, displayAll(), em_, PFRootEventManager::eventNumber(), eventNumber_, isGraphicLoaded_, loadGraphicObjects(), convertSQLiteXML::ok, PFRootEventManager::processEntry(), and reset().
{
bool ok=false;
while (!ok && ientry<em_->ev_->size() ) {
ok = em_->processEntry(ientry);
ientry++;
}
eventNumber_ = em_->eventNumber();
if (ok) {
reset();
loadGraphicObjects();
isGraphicLoaded_= true;
displayAll();
}
else
std::cerr<<"DisplayManager::dislayNextInteresting : no event matching criteria"<<std::endl;
}
| void DisplayManager::displayPFBlock | ( | int | blockNb | ) |
Definition at line 1298 of file DisplayManager.cc.
References blockIdentsMap_, create_public_lumi_plots::color, gather_cfg::cout, drawGObject(), drawPFBlocks_, first, AlCaHLTBitMon_ParallelJobs::p, rubOutGPFBlock(), and selectedGObj_.
Referenced by findBlock().
{
rubOutGPFBlock();
selectedGObj_.clear();
if (!drawPFBlocks_) return;
int color=1;
multimap<int,pair <int,int > >::const_iterator p;
p= blockIdentsMap_.find(blockNb);
if (p !=blockIdentsMap_.end()) {
do {
int ident=(p->second).first;
drawGObject(ident,color,false);
p++;
} while (p!=blockIdentsMap_.upper_bound(blockNb));
}
else
cout<<"DisplayManager::displayPFBlock :not found"<<endl;
}
| void DisplayManager::displayPrevious | ( | ) |
Definition at line 1283 of file DisplayManager.cc.
References display(), em_, PFRootEventManager::eventNumber(), and eventNumber_.
Referenced by DialogFrame::doPreviousEvent().
{
int eventNumber_=em_->eventNumber();
display(--eventNumber_);
}
| void DisplayManager::drawGObject | ( | int | ident, |
| int | color, | ||
| bool | toInitialColor | ||
| ) | [private] |
Definition at line 1317 of file DisplayManager.cc.
References gather_cfg::cout, displayView_, drawHO_, EHO, graphicMap_, AlCaHLTBitMon_ParallelJobs::p, query::result, and selectedGObj_.
Referenced by displayPFBlock(), findAndDraw(), and rubOutGPFBlock().
{
typedef std::multimap<int,GPFBase *>::const_iterator iter;
iter p;
std::pair<iter, iter > result = graphicMap_.equal_range(ident);
if(result.first == graphicMap_.end()) {
std::cout<<"pas d'objet avec cet ident: "<<ident<<std::flush<<std::endl;
return;
}
p=result.first;
while (p != result.second) {
int view=p->second->getView();
if (!drawHO_ && view ==EHO) continue;
displayView_[view]->cd();
if (toInitial) p->second->setInitialColor();
else p->second->setColor(color);
p->second->draw();
gPad->Modified();
// displayView_[view]->Update();
if (!toInitial) selectedGObj_.push_back(ident);
p++;
}
}
| void DisplayManager::drawWithNewGraphicAttributes | ( | ) |
Definition at line 1158 of file DisplayManager.cc.
References CLUSTERECALID, CLUSTERHCALID, CLUSTERHOID, CLUSTERIBID, CLUSTERPSID, displayAll(), drawHO_, graphicMap_, AlCaHLTBitMon_ParallelJobs::p, RECTRACKID, shiftId_, and SIMPARTICLEID.
Referenced by DialogFrame::modifyAttr().
{
std::multimap<int,GPFBase *>::iterator p;
for (p=graphicMap_.begin();p!=graphicMap_.end();p++) {
int ident=p->first;
int type=ident >> shiftId_;
if (!drawHO_ && (type==CLUSTERHOID)) continue;
switch (type) {
case CLUSTERECALID: case CLUSTERHCALID: case CLUSTERHOID: case CLUSTERPSID: case CLUSTERIBID:
{
p->second->setNewStyle();
p->second->setNewSize();
p->second->setColor();
}
break;
case RECTRACKID:
{
p->second->setColor();
p->second->setNewStyle();
p->second->setNewSize();
}
break;
case SIMPARTICLEID:
{
}
break;
default : break;
} //switch end
} //for end
displayAll(false);
}
| void DisplayManager::enableDrawBrem | ( | bool | state | ) |
Definition at line 1346 of file DisplayManager.cc.
References drawBrems_, and evf::utils::state.
Referenced by DialogFrame::areBremVisible().
{
drawBrems_=state;
}
| void DisplayManager::enableDrawPFBlock | ( | bool | state | ) |
Definition at line 1341 of file DisplayManager.cc.
References drawPFBlocks_, and evf::utils::state.
Referenced by DialogFrame::isPFBlockVisible().
{
drawPFBlocks_=state;
}
| void DisplayManager::findAndDraw | ( | int | ident | ) |
Definition at line 1351 of file DisplayManager.cc.
References create_public_lumi_plots::color, gather_cfg::cout, displayAll(), drawGObject(), drawHO_, drawPFBlocks_, HITTYPES, RECHITECALID, RECHITHCALID, RECHITHFEMID, RECHITHFHADID, RECHITHOID, RECHITPSID, rubOutGPFBlock(), selectedGObj_, shiftId_, SIMPARTICLEID, and updateDisplay().
Referenced by GPFCluster::ExecuteEvent(), GPFRecHit::ExecuteEvent(), GPFSimParticle::ExecuteEvent(), and GPFTrack::ExecuteEvent().
{
int type=ident >> shiftId_; // defined in DisplayCommon.h
if (drawHO_ || type!=RECHITHOID) {
int color=1;
if (type>15) {
std ::cout<<"DisplayManager::findAndDraw :object Type unknown"<<std::endl;
return;
}
if (drawPFBlocks_==0 ||
type==RECHITECALID || type==RECHITHCALID || type==RECHITHOID ||
type==RECHITHFEMID || type==RECHITHFHADID ||
type==RECHITPSID || type==SIMPARTICLEID) {
rubOutGPFBlock();
selectedGObj_.clear();
bool toInitial=false;
drawGObject(ident,color,toInitial);
if (type<HITTYPES) {
//redrawWithoutHits_=true;
displayAll(false);
//redrawWithoutHits_=false;
}
}
updateDisplay();
}
}
| bool DisplayManager::findBadBremsId | ( | int | ident | ) | [private] |
Definition at line 1401 of file DisplayManager.cc.
References badBremsId_, and i.
Referenced by loadGGsfRecTracks(), and loadGPFBlocks().
{
for (unsigned i=0;i<badBremsId_.size();i++)
if (badBremsId_[i]==ident) return true;
return false;
}
| void DisplayManager::findBlock | ( | int | ident | ) |
Definition at line 1378 of file DisplayManager.cc.
References Association::block, blockIdentsMap_, PFRootEventManager::blocks(), gather_cfg::cout, displayPFBlock(), em_, first, AlCaHLTBitMon_ParallelJobs::p, edm::second(), and updateDisplay().
Referenced by GPFCluster::ExecuteEvent(), and GPFTrack::ExecuteEvent().
{
int blockNb=-1;
int elemNb=-1;
multimap<int, pair <int,int > >::const_iterator p;
for (p=blockIdentsMap_.begin();p!=blockIdentsMap_.end();p++) {
int id=(p->second).first;
if (id == ident) {
blockNb=p->first;
elemNb=(p->second).second;
break;
}
}
if (blockNb > -1) {
std::cout<<"this object is element "<<elemNb<<" of PFblock nb "<<blockNb<<std::endl;
assert( blockNb < static_cast<int>(em_->blocks().size()) );
const reco::PFBlock& block = em_->blocks()[blockNb];
std::cout<<block<<std::endl;
displayPFBlock(blockNb);
}
updateDisplay();
}
| void DisplayManager::getDisplayOptions | ( | ) | [private] |
| double DisplayManager::getMaxE | ( | int | layer | ) | const [private] |
Definition at line 1421 of file DisplayManager.cc.
References dtNoiseDBValidation_cfg::cerr, PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, em_, PFLayer::HCAL_BARREL1, PFLayer::HCAL_BARREL2, PFLayer::HCAL_ENDCAP, PFLayer::HF_EM, PFLayer::HF_HAD, i, PFLayer::PS1, PFLayer::PS2, PFRootEventManager::rechitsECAL_, PFRootEventManager::rechitsHCAL_, PFRootEventManager::rechitsHFEM_, PFRootEventManager::rechitsHFHAD_, PFRootEventManager::rechitsHO_, and PFRootEventManager::rechitsPS_.
Referenced by getMaxEEcal(), getMaxEHcal(), and getMaxEHo().
{
double maxe = -9999;
// in which vector should we look for these rechits ?
const reco::PFRecHitCollection* vec = 0;
switch(layer) {
case PFLayer::ECAL_ENDCAP:
case PFLayer::ECAL_BARREL:
vec = &(em_->rechitsECAL_);
break;
case PFLayer::HCAL_ENDCAP:
case PFLayer::HCAL_BARREL1:
// case PFLayer::HCAL_BARREL2:
vec = &(em_->rechitsHCAL_);
break;
case PFLayer::HCAL_BARREL2:
vec = &(em_->rechitsHO_);
break;
case PFLayer::HF_EM:
vec = &(em_->rechitsHFEM_);
break;
case PFLayer::HF_HAD:
vec = &(em_->rechitsHFHAD_);
break;
case PFLayer::PS1:
case PFLayer::PS2:
vec = &(em_->rechitsPS_);
break;
default:
cerr<<"DisplayManager::getMaxE : manage other layers"<<endl;
return maxe;
}
for( unsigned i=0; i<vec->size(); i++) {
if( (*vec)[i].layer() != layer ) continue;
if( (*vec)[i].energy() > maxe)
maxe = (*vec)[i].energy();
}
return maxe;
}
| double DisplayManager::getMaxEEcal | ( | ) | [private] |
Definition at line 1467 of file DisplayManager.cc.
References PFLayer::ECAL_BARREL, PFLayer::ECAL_ENDCAP, getMaxE(), PFLayer::HF_EM, and maxERecHitEcal_.
Referenced by loadGRecHits().
{
if( maxERecHitEcal_<0 ) {
double maxeec = getMaxE( PFLayer::ECAL_ENDCAP );
double maxeb = getMaxE( PFLayer::ECAL_BARREL );
double maxehf = getMaxE( PFLayer::HF_EM );
maxERecHitEcal_ = maxeec>maxeb ? maxeec:maxeb;
maxERecHitEcal_ = maxERecHitEcal_>maxehf ? maxERecHitEcal_:maxehf;
// max of both barrel and endcap
}
return maxERecHitEcal_;
}
| double DisplayManager::getMaxEHcal | ( | ) | [private] |
Definition at line 1480 of file DisplayManager.cc.
References getMaxE(), PFLayer::HCAL_BARREL1, PFLayer::HCAL_ENDCAP, PFLayer::HF_HAD, and maxERecHitHcal_.
Referenced by loadGRecHits().
{
if(maxERecHitHcal_ < 0) {
double maxehf = getMaxE( PFLayer::HF_HAD );
double maxeec = getMaxE( PFLayer::HCAL_ENDCAP );
double maxeb = getMaxE( PFLayer::HCAL_BARREL1 );
maxERecHitHcal_ = maxeec>maxeb ? maxeec:maxeb;
maxERecHitHcal_ = maxERecHitHcal_>maxehf ? maxERecHitHcal_:maxehf;
}
return maxERecHitHcal_;
}
| double DisplayManager::getMaxEHo | ( | ) | [private] |
Definition at line 1493 of file DisplayManager.cc.
References getMaxE(), PFLayer::HCAL_BARREL2, and maxERecHitHo_.
{
if(maxERecHitHo_ < 0) {
maxERecHitHo_ = getMaxE( PFLayer::HCAL_BARREL2 );
}
return maxERecHitHo_;
}
| void DisplayManager::loadGClusters | ( | ) | [private] |
Definition at line 1615 of file DisplayManager.cc.
References CLUSTERECALID, CLUSTERHCALID, CLUSTERHFEMID, CLUSTERHFHADID, CLUSTERHOID, CLUSTERPSID, PFRootEventManager::clustersECAL_, PFRootEventManager::clustersHCAL_, PFRootEventManager::clustersHFEM_, PFRootEventManager::clustersHFHAD_, PFRootEventManager::clustersHO_, PFRootEventManager::clustersPS_, createGCluster(), drawHO_, em_, i, and shiftId_.
Referenced by loadGraphicObjects().
{
double phi0=0;
for(unsigned i=0; i<em_->clustersECAL_->size(); i++){
//int clusId=(i<<shiftId_) | CLUSTERECALID;
int clusId=(CLUSTERECALID<<shiftId_) | i;
createGCluster( (*(em_->clustersECAL_))[i],clusId, phi0);
}
for(unsigned i=0; i<em_->clustersHCAL_->size(); i++) {
//int clusId=(i<<shiftId_) | CLUSTERHCALID;
int clusId=(CLUSTERHCALID<<shiftId_) | i;
createGCluster( (*(em_->clustersHCAL_))[i],clusId, phi0);
}
if (drawHO_) {
for(unsigned i=0; i<em_->clustersHO_->size(); i++) {
//int clusId=(i<<shiftId_) | CLUSTERHOID;
int clusId=(CLUSTERHOID<<shiftId_) | i;
createGCluster( (*(em_->clustersHO_))[i],clusId, phi0);
}
}
for(unsigned i=0; i<em_->clustersHFEM_->size(); i++) {
//int clusId=(i<<shiftId_) | CLUSTERHFEMID;
int clusId=(CLUSTERHFEMID<<shiftId_) | i;
createGCluster( (*(em_->clustersHFEM_))[i],clusId, phi0);
}
for(unsigned i=0; i<em_->clustersHFHAD_->size(); i++) {
//int clusId=(i<<shiftId_) | CLUSTERHFHADID;
int clusId=(CLUSTERHFHADID<<shiftId_) | i;
createGCluster( (*(em_->clustersHFHAD_))[i],clusId, phi0);
}
for(unsigned i=0; i<em_->clustersPS_->size(); i++){
//int clusId=(i<<shiftId_) | CLUSTERPSID;
int clusId=(CLUSTERPSID<<shiftId_) | i;
createGCluster( (*(em_->clustersPS_))[i],clusId,phi0);
}
// for(unsigned i=0; i<em_->clustersIslandBarrel_.size(); i++) {
// PFLayer::Layer layer = PFLayer::ECAL_BARREL;
// //int clusId=(i<<shiftId_) | CLUSTERIBID;
// int clusId=(CLUSTERIBID<<shiftId_) | i;
// reco::PFCluster cluster( layer,
// em_->clustersIslandBarrel_[i].energy(),
// em_->clustersIslandBarrel_[i].x(),
// em_->clustersIslandBarrel_[i].y(),
// em_->clustersIslandBarrel_[i].z() );
// createGCluster( cluster,clusId, phi0);
// }
}
| void DisplayManager::loadGGenParticles | ( | ) | [private] |
Definition at line 1503 of file DisplayManager.cc.
References createGGenParticle(), em_, configurableAnalysis::GenParticle, edm::HepMCProduct::GetEvent(), PFRootEventManager::MCTruth_, and AlCaHLTBitMon_ParallelJobs::p.
Referenced by loadGraphicObjects().
{
const HepMC::GenEvent* myGenEvent = em_->MCTruth_.GetEvent();
if(!myGenEvent) return;
for ( HepMC::GenEvent::particle_const_iterator piter = myGenEvent->particles_begin();
piter != myGenEvent->particles_end();
++piter ) {
HepMC::GenParticle* p = *piter;
createGGenParticle(p);
}
}
| void DisplayManager::loadGGsfRecTracks | ( | ) | [private] |
Definition at line 1980 of file DisplayManager.cc.
References BREMID, reco::PFTrajectoryPoint::ClosestApproach, funct::cos(), createGTrack(), reco::PFTrajectoryPoint::ECALEntrance, em_, findBadBremsId(), GSFRECTRACKID, PFRootEventManager::gsfrecTracks_, create_public_lumi_plots::linestyle, reco::PFTrajectoryPoint::momentum(), retrieveBadBrems(), and shiftId_.
Referenced by loadGraphicObjects().
{
double phi0=0;
int ind=-1;
int indbrem=-1;
// allows not to draw Brems with no informations in PFBlocks
retrieveBadBrems();
std::vector<reco::GsfPFRecTrack>::iterator itRecTrack;
for (itRecTrack = em_->gsfrecTracks_.begin(); itRecTrack != em_->gsfrecTracks_.end();itRecTrack++) {
double sign = 1.;
const reco::PFTrajectoryPoint& tpinitial
= itRecTrack->extrapolatedPoint(reco::PFTrajectoryPoint::ClosestApproach);
double pt = tpinitial.momentum().Pt();
//if( pt<em_->displayRecTracksPtMin_ ) continue;
const reco::PFTrajectoryPoint& tpatecal
= itRecTrack->trajectoryPoint(itRecTrack->nTrajectoryMeasurements() +
reco::PFTrajectoryPoint::ECALEntrance );
if ( cos(phi0 - tpatecal.momentum().Phi()) < 0.)
sign = -1.;
const std::vector<reco::PFTrajectoryPoint>& points =
itRecTrack->trajectoryPoints();
int linestyle = itRecTrack->algoType();
ind++;
int recTrackId=(GSFRECTRACKID <<shiftId_) | ind;
createGTrack(*itRecTrack,points,recTrackId, pt, phi0, sign, false,linestyle,1);
// now the Brems - bad to copy, but problems otherwise
std::vector<reco::PFBrem> brems=itRecTrack->PFRecBrem();
unsigned nbrems=brems.size();
for(unsigned ibrem=0;ibrem<nbrems;++ibrem)
{
unsigned indTrajPoint=brems[ibrem].indTrajPoint();
if(indTrajPoint==99) continue;
double signBrem = 1. ; // this is not the charge
int linestyleBrem = brems[ibrem].algoType();
indbrem++;
// build the index by hand. Assume that there are less 40 Brems per GSF track (ultrasafe!)
// make it from the GSF index and the brem index
unsigned indexBrem= ind*40+brems[ibrem].indTrajPoint();
int recTrackIdBrem=(BREMID << shiftId_ ) | indexBrem;
//check if there is information on this brem in the PFBlock
//before creating a graphic object
if (!findBadBremsId(recTrackIdBrem)) {
// the vertex is not stored, need to make it by hand
std::vector<reco::PFTrajectoryPoint> pointsBrem;
//the first two trajectory points are dummy; copy them
pointsBrem.push_back(brems[ibrem].trajectoryPoints()[0]);
pointsBrem.push_back(brems[ibrem].trajectoryPoints()[1]);
// then get the vertex from the GSF track
pointsBrem.push_back(itRecTrack->trajectoryPoint(indTrajPoint));
unsigned ntp=brems[ibrem].trajectoryPoints().size();
for(unsigned itp=2;itp<ntp;++itp)
{
pointsBrem.push_back(brems[ibrem].trajectoryPoints()[itp]);
}
double deltaP=brems[ibrem].DeltaP();
const reco::PFTrajectoryPoint& tpatecalbrem
= brems[ibrem].trajectoryPoint(brems[ibrem].nTrajectoryMeasurements() +
reco::PFTrajectoryPoint::ECALEntrance );
if ( cos(phi0 - tpatecalbrem.momentum().Phi()) < 0.)
signBrem = -1.; // again, not the charge
createGTrack(brems[ibrem],pointsBrem,recTrackIdBrem,deltaP,phi0,signBrem,false,linestyleBrem,2);
}
}
}
}
| void DisplayManager::loadGPFBlocks | ( | ) | [private] |
Definition at line 1701 of file DisplayManager.cc.
References begin, blockIdentsMap_, reco::PFBlockElement::BREM, BREMID, CLUSTERECALID, CLUSTERHCALID, CLUSTERHFEMID, CLUSTERHFHADID, CLUSTERHOID, CLUSTERIBID, CLUSTERPSID, gather_cfg::cout, drawHO_, reco::PFBlockElement::ECAL, em_, end, findBadBremsId(), reco::PFBlockElement::GSF, GSFRECTRACKID, reco::PFBlockElementBrem::GsftrackRefPF(), reco::PFBlockElementGsfTrack::GsftrackRefPF(), reco::PFBlockElement::HCAL, reco::PFBlockElement::HFEM, reco::PFBlockElement::HFHAD, reco::PFBlockElement::HO, getHLTprescales::index, reco::PFBlockElementBrem::indTrajPoint(), edm::Ref< C, T, F >::isNull(), edm::Ref< C, T, F >::key(), reco::PFBlockElement::NONE, PFRootEventManager::pfBlocks_, reco::PFBlockElement::PS1, reco::PFBlockElement::PS2, RECTRACKID, reco::PFBlockElement::SC, shiftId_, findQualityFiles::size, reco::PFBlockElementSuperCluster::superClusterRef(), and reco::PFBlockElement::TRACK.
Referenced by loadGraphicObjects().
{
int size = em_->pfBlocks_->size();
for (int ibl=0;ibl<size;ibl++) {
//int elemNb=((*(em_->pfBlocks_))[ibl].elements()).size();
//std::cout<<"block "<<ibl<<":"<<elemNb<<" elements"<<std::flush<<std::endl;
edm::OwnVector< reco::PFBlockElement >::const_iterator iter;
for( iter =((*(em_->pfBlocks_))[ibl].elements()).begin();
iter != ((*(em_->pfBlocks_))[ibl].elements()).end();iter++) {
//COLIN
// std::cout<<"elem index "<<(*iter).index()<<"-type:"
// <<(*iter).type()<<std::flush<<std::endl;
int ident=-1;
reco::PFBlockElement::Type type = (*iter).type();
if (!drawHO_ && type==reco::PFBlockElement::HO) continue;
switch (type) {
case reco::PFBlockElement::NONE :
assert(0);
break;
case reco::PFBlockElement::TRACK:
{
reco::PFRecTrackRef trackref =(*iter).trackRefPF();
assert( !trackref.isNull() );
// std::cout<<" - key "<<trackref.key()<<std::flush<<std::endl<<std::endl;
ident=(RECTRACKID <<shiftId_) | trackref.key();
}
break;
case reco::PFBlockElement::PS1:
{
reco::PFClusterRef clusref=(*iter).clusterRef();
assert( !clusref.isNull() );
//std::cout<<"- key "<<clusref.key()<<std::flush<<std::endl<<std::endl;
ident=(CLUSTERPSID <<shiftId_) |clusref.key();
}
break;
case reco::PFBlockElement::PS2:
{
reco::PFClusterRef clusref=(*iter).clusterRef();
assert( !clusref.isNull() );
//std::cout<<"key "<<clusref.key()<<std::flush<<std::endl<<std::endl;
ident=(CLUSTERPSID <<shiftId_) |clusref.key();
}
break;
case reco::PFBlockElement::ECAL:
{
reco::PFClusterRef clusref=(*iter).clusterRef();
assert( !clusref.isNull() );
//std::cout<<"key "<<clusref.key()<<std::flush<<std::endl<<std::endl;
ident=(CLUSTERECALID <<shiftId_) |clusref.key();
}
break;
case reco::PFBlockElement::HCAL:
{
reco::PFClusterRef clusref=(*iter).clusterRef();
assert( !clusref.isNull() );
//std::cout<<"key "<<clusref.key()<<std::flush<<std::endl<<std::endl;
ident=(CLUSTERHCALID <<shiftId_) |clusref.key();
}
break;
case reco::PFBlockElement::HO:
{
reco::PFClusterRef clusref=(*iter).clusterRef();
assert( !clusref.isNull() );
//std::cout<<"key "<<clusref.key()<<std::flush<<std::endl<<std::endl;
ident=(CLUSTERHOID <<shiftId_) |clusref.key();
}
break;
case reco::PFBlockElement::HFEM:
{
reco::PFClusterRef clusref=(*iter).clusterRef();
assert( !clusref.isNull() );
//std::cout<<"key "<<clusref.key()<<std::flush<<std::endl<<std::endl;
ident=(CLUSTERHFEMID <<shiftId_) |clusref.key();
}
break;
case reco::PFBlockElement::HFHAD:
{
reco::PFClusterRef clusref=(*iter).clusterRef();
assert( !clusref.isNull() );
//std::cout<<"key "<<clusref.key()<<std::flush<<std::endl<<std::endl;
ident=(CLUSTERHFHADID <<shiftId_) |clusref.key();
}
break;
case reco::PFBlockElement::GSF:
{
const reco::PFBlockElementGsfTrack * GsfEl =
dynamic_cast<const reco::PFBlockElementGsfTrack*>(&(*iter));
reco::GsfPFRecTrackRef trackref=GsfEl->GsftrackRefPF();
assert( !trackref.isNull() );
ident=(GSFRECTRACKID << shiftId_) | trackref.key();
}
break;
case reco::PFBlockElement::BREM:
{
const reco::PFBlockElementBrem * Brem = dynamic_cast<const reco::PFBlockElementBrem*>(&(*iter));
reco::GsfPFRecTrackRef trackref = Brem->GsftrackRefPF();
unsigned index=trackref.key()*40+Brem->indTrajPoint();
ident = (BREMID << shiftId_ ) | index ;
if (findBadBremsId(ident)) ident=-1;
}
break;
case reco::PFBlockElement::SC:
{
const reco::PFBlockElementSuperCluster * sc = dynamic_cast<const reco::PFBlockElementSuperCluster*>(&(*iter));
reco::SuperClusterRef scref = sc->superClusterRef();
assert( !scref.isNull() );
ident = (CLUSTERIBID << shiftId_ ) | scref.key();
}
break;
default:
std::cout<<"unknown PFBlock element of type "<<type<<std::endl;
break;
} //end switch
pair <int, int> idElem;
idElem.first=ident;
idElem.second=(*iter).index();
if (ident != -1) blockIdentsMap_.insert(pair<int,pair <int,int> > (ibl,idElem));
} //end for elements
} //end for blocks
}
| void DisplayManager::loadGraphicObjects | ( | ) | [private] |
Definition at line 1830 of file DisplayManager.cc.
References loadGClusters(), loadGGenParticles(), loadGGsfRecTracks(), loadGPFBlocks(), loadGRecHits(), loadGRecTracks(), and loadGSimParticles().
Referenced by display(), displayEvent(), and displayNextInteresting().
{
loadGClusters();
loadGRecHits();
loadGRecTracks();
loadGGsfRecTracks();
loadGSimParticles();
loadGPFBlocks();
loadGGenParticles();
}
| void DisplayManager::loadGRecHits | ( | ) | [private] |
Definition at line 1841 of file DisplayManager.cc.
References dtNoiseDBValidation_cfg::cerr, PFRootEventManager::clusterAlgoECAL_, PFRootEventManager::clusterAlgoHCAL_, PFRootEventManager::clusterAlgoHFEM_, PFRootEventManager::clusterAlgoHFHAD_, PFRootEventManager::clusterAlgoHO_, PFRootEventManager::clusterAlgoPS_, create_public_lumi_plots::color, PFClusterAlgo::color(), createGRecHit(), drawHO_, em_, getMaxEEcal(), getMaxEHcal(), i, RECHITHCALID, RECHITHFEMID, RECHITHFHADID, RECHITHOID, RECHITPSID, PFRootEventManager::rechitsECAL_, PFRootEventManager::rechitsHCAL_, PFRootEventManager::rechitsHFEM_, PFRootEventManager::rechitsHFHAD_, PFRootEventManager::rechitsHO_, PFRootEventManager::rechitsPS_, PFClusterAlgo::SEED, shiftId_, and PFClusterAlgo::SPECIAL.
Referenced by loadGraphicObjects().
{
double phi0=0;
double maxee = getMaxEEcal();
double maxeh = getMaxEHcal();
double maxe = maxee>maxeh ? maxee : maxeh;
int color = TColor::GetColor(210,210,210);
int seedcolor = TColor::GetColor(145,145,145);
int specialcolor = TColor::GetColor(255,140,0);
for(unsigned i=0; i<em_->rechitsECAL_.size(); i++) {
int rhcolor = color;
if( unsigned col = em_->clusterAlgoECAL_.color(i) ) {
switch(col) {
case PFClusterAlgo::SEED: rhcolor = seedcolor; break;
case PFClusterAlgo::SPECIAL: rhcolor = specialcolor; break;
default:
cerr<<"DisplayManager::loadGRecHits: unknown color"<<endl;
}
}
//int recHitId=(i<<shiftId_) | RECHITECALID;
int recHitId=i;
createGRecHit(em_->rechitsECAL_[i],recHitId, maxe, phi0, rhcolor);
}
for(unsigned i=0; i<em_->rechitsHCAL_.size(); i++) {
int rhcolor = color;
if(unsigned col = em_->clusterAlgoHCAL_.color(i) ) {
switch(col) {
case PFClusterAlgo::SEED: rhcolor = seedcolor; break;
case PFClusterAlgo::SPECIAL: rhcolor = specialcolor; break;
default:
cerr<<"DisplayManager::loadGRecHits: unknown color"<<endl;
}
}
//int recHitId=(i<<shiftId_) | RECHITHCALID;
int recHitId=(RECHITHCALID <<shiftId_) | i;
createGRecHit(em_->rechitsHCAL_[i],recHitId, maxe, phi0, rhcolor);
}
if (drawHO_) {
for(unsigned i=0; i<em_->rechitsHO_.size(); i++) {
int rhcolor = color;
if(unsigned col = em_->clusterAlgoHO_.color(i) ) {
switch(col) {
case PFClusterAlgo::SEED: rhcolor = seedcolor; break;
case PFClusterAlgo::SPECIAL: rhcolor = specialcolor; break;
default:
cerr<<"DisplayManager::loadGRecHits: unknown color"<<endl;
}
}
//int recHitId=(i<<shiftId_) | RECHITHOID;
int recHitId=(RECHITHOID <<shiftId_) | i;
createGRecHit(em_->rechitsHO_[i],recHitId, 2*maxe, phi0, rhcolor);
}
}
for(unsigned i=0; i<em_->rechitsHFEM_.size(); i++) {
int rhcolor = color;
if(unsigned col = em_->clusterAlgoHFEM_.color(i) ) {
switch(col) {
case PFClusterAlgo::SEED: rhcolor = seedcolor; break;
case PFClusterAlgo::SPECIAL: rhcolor = specialcolor; break;
default:
cerr<<"DisplayManager::loadGRecHits: unknown color"<<endl;
}
}
//int recHitId=(i<<shiftId_) | RECHITHFEMID;
int recHitId=(RECHITHFEMID <<shiftId_) | i;
createGRecHit(em_->rechitsHFEM_[i],recHitId, maxe, phi0, rhcolor);
}
for(unsigned i=0; i<em_->rechitsHFHAD_.size(); i++) {
int rhcolor = color;
if(unsigned col = em_->clusterAlgoHFHAD_.color(i) ) {
switch(col) {
case PFClusterAlgo::SEED: rhcolor = seedcolor; break;
case PFClusterAlgo::SPECIAL: rhcolor = specialcolor; break;
default:
cerr<<"DisplayManager::loadGRecHits: unknown color"<<endl;
}
}
//int recHitId=(i<<shiftId_) | RECHITHFHADID;
int recHitId=(RECHITHFHADID <<shiftId_) | i;
createGRecHit(em_->rechitsHFHAD_[i],recHitId, maxe, phi0, rhcolor);
}
for(unsigned i=0; i<em_->rechitsPS_.size(); i++) {
int rhcolor = color;
if( unsigned col = em_->clusterAlgoPS_.color(i) ) {
switch(col) {
case PFClusterAlgo::SEED: rhcolor = seedcolor; break;
case PFClusterAlgo::SPECIAL: rhcolor = specialcolor; break;
default:
cerr<<"DisplayManager::loadGRecHits: unknown color"<<endl;
}
}
//int recHitId=(i<<shiftId_) | RECHITPSID;
int recHitId=(RECHITPSID<<shiftId_) | i;
createGRecHit(em_->rechitsPS_[i],recHitId, maxe, phi0, rhcolor);
}
}
| void DisplayManager::loadGRecTracks | ( | ) | [private] |
Definition at line 1947 of file DisplayManager.cc.
References reco::PFTrajectoryPoint::ClosestApproach, funct::cos(), createGTrack(), reco::PFTrajectoryPoint::ECALEntrance, em_, create_public_lumi_plots::linestyle, reco::PFTrajectoryPoint::momentum(), RECTRACKID, PFRootEventManager::recTracks_, and shiftId_.
Referenced by loadGraphicObjects().
{
double phi0=0;
int ind=-1;
std::vector<reco::PFRecTrack>::iterator itRecTrack;
for (itRecTrack = em_->recTracks_.begin(); itRecTrack != em_->recTracks_.end();itRecTrack++) {
double sign = 1.;
const reco::PFTrajectoryPoint& tpinitial
= itRecTrack->extrapolatedPoint(reco::PFTrajectoryPoint::ClosestApproach);
double pt = tpinitial.momentum().Pt();
//if( pt<em_->displayRecTracksPtMin_ ) continue;
const reco::PFTrajectoryPoint& tpatecal
= itRecTrack->trajectoryPoint(itRecTrack->nTrajectoryMeasurements() +
reco::PFTrajectoryPoint::ECALEntrance );
if ( cos(phi0 - tpatecal.momentum().Phi()) < 0.)
sign = -1.;
const std::vector<reco::PFTrajectoryPoint>& points =
itRecTrack->trajectoryPoints();
int linestyle = itRecTrack->algoType();
ind++;
//int recTrackId=(ind<<shiftId_) | RECTRACKID;
int recTrackId=(RECTRACKID <<shiftId_) | ind;
createGTrack(*itRecTrack,points,recTrackId, pt, phi0, sign, false,linestyle);
}
}
| void DisplayManager::loadGSimParticles | ( | ) | [private] |
Definition at line 2064 of file DisplayManager.cc.
References abs, reco::PFTrajectoryPoint::ClosestApproach, createGPart(), em_, reco::PFTrack::extrapolatedPoint(), i, reco::PFTrajectoryPoint::momentum(), reco::PFSimParticle::motherId(), reco::PFSimParticle::pdgCode(), reco::PFTrajectoryPoint::position(), shiftId_, SIMPARTICLEID, reco::PFTrack::trajectoryPoint(), reco::PFTrack::trajectoryPoints(), and PFRootEventManager::trueParticles_.
Referenced by loadGraphicObjects().
{
double phi0=0;
// bool debug_loadGSim = true;
unsigned simParticlesVSize = em_->trueParticles_.size();
for(unsigned i=0; i<simParticlesVSize; i++) {
const reco::PFSimParticle& ptc = em_->trueParticles_[i];
const reco::PFTrajectoryPoint& tpinitial
= ptc.extrapolatedPoint( reco::PFTrajectoryPoint::ClosestApproach );
double pt = tpinitial.momentum().Pt();
//if( pt<em_->getDisplayTrueParticlesPtMin()) continue;
double sign = 1.;
const reco::PFTrajectoryPoint& tpFirst = ptc.trajectoryPoint(0);
if ( tpFirst.position().X() < 0. )
sign = -1.;
//double sign2 = 1.;
// if position vector is 0,0,0: X component is undefined (MDN 11 June 2008)
// const reco::PFTrajectoryPoint& tpFirst = ptc.trajectoryPoint(0);
// if ( tpFirst.positionXYZ().X() < 0. )
// sign2 = -1.;
const std::vector<reco::PFTrajectoryPoint>& points =
ptc.trajectoryPoints();
int indexMarker;
switch( std::abs(ptc.pdgCode() ) ) {
case 22: indexMarker=0; break; // photons
case 11: indexMarker=1; break; // electrons
case 13: indexMarker=2; break; // muons
case 130:
case 321: indexMarker=3; break; // K
case 211: indexMarker=4; break; // pi+/pi-
case 2212: indexMarker=5; break; // protons
case 2112: indexMarker=6; break; // neutrons
default: indexMarker=7; break;
}
bool displayInitial=true;
if( ptc.motherId() < 0 ) displayInitial=false;
int partId=(SIMPARTICLEID << shiftId_) | i;
createGPart(ptc, points,partId, pt, phi0, sign, displayInitial,indexMarker);
//cout << " sign " << sign << " sign2 " << sign2 << endl;
//if ( sign*sign2 <0 ) cout << " ++++++Warning sign*sign2 <0 ++++++++++++++++++ " << endl;
}
}
| void DisplayManager::lookForGenParticle | ( | unsigned | barcode | ) |
look for particle with index i in MC truth.
Definition at line 2180 of file DisplayManager.cc.
References dtNoiseDBValidation_cfg::cerr, displayHist_, displayView_, drawHO_, EHO, em_, EPE, EPH, eta(), event(), configurableAnalysis::GenParticle, edm::HepMCProduct::GetEvent(), PFRootEventManager::MCTruth_, phi, and zoomFactor_.
Referenced by DialogFrame::doLookForGenParticle().
{
const HepMC::GenEvent* event = em_->MCTruth_.GetEvent();
if(!event) {
cerr<<"no GenEvent"<<endl;
return;
}
const HepMC::GenParticle* particle = event->barcode_to_particle(barcode);
if(!particle) {
cerr<<"no particle with barcode "<<barcode<<endl;
return;
}
math::XYZTLorentzVector momentum(particle->momentum().px(),
particle->momentum().py(),
particle->momentum().pz(),
particle->momentum().e());
double eta = momentum.Eta();
double phi = momentum.phi();
double phisize = 0.05;
double etasize = 0.05;
double etagate = zoomFactor_ * etasize;
double phigate = zoomFactor_ * phisize;
if(displayHist_[EPE]) {
displayHist_[EPE]->GetXaxis()->SetRangeUser(eta-etagate, eta+etagate);
displayHist_[EPE]->GetYaxis()->SetRangeUser(phi-phigate, phi+phigate);
displayView_[EPE]->Modified();
displayView_[EPE]->Update();
}
if(displayHist_[EPH]) {
displayHist_[EPH]->GetXaxis()->SetRangeUser(eta-etagate, eta+etagate);
displayHist_[EPH]->GetYaxis()->SetRangeUser(phi-phigate, phi+phigate);
displayView_[EPH]->Modified();
displayView_[EPH]->Update();
}
if(drawHO_ && displayHist_[EHO]) {
displayHist_[EHO]->GetXaxis()->SetRangeUser(eta-etagate, eta+etagate);
displayHist_[EHO]->GetYaxis()->SetRangeUser(phi-phigate, phi+phigate);
displayView_[EHO]->Modified();
displayView_[EHO]->Update();
}
}
| void DisplayManager::lookForMaxRecHit | ( | bool | ecal | ) |
look for rechit with max energy in ecal or hcal.
Definition at line 2117 of file DisplayManager.cc.
References displayHist_, displayView_, drawHO_, EHO, em_, relval_parameters_module::energy, EPE, EPH, eta(), i, phi, reco::PFRecHit::position(), HI_PhotonSkim_cff::rechits, PFRootEventManager::rechitsECAL_, PFRootEventManager::rechitsHCAL_, reco::PFRecHit::size(), and zoomFactor_.
{
// look for the rechit with max e in ecal or hcal
double maxe = -999;
reco::PFRecHit* maxrh = 0;
reco::PFRecHitCollection* rechits = 0;
if(ecal) rechits = &(em_->rechitsECAL_);
else rechits = &(em_->rechitsHCAL_);
assert(rechits);
for(unsigned i=0; i<(*rechits).size(); i++) {
double energy = (*rechits)[i].energy();
if(energy > maxe ) {
maxe = energy;
maxrh = &((*rechits)[i]);
}
}
if(!maxrh) return;
// center view on this rechit
// get the cell size to set the eta and phi width
// of the display window from one of the cells
double phisize = -1;
double etasize = -1;
maxrh->size(phisize, etasize);
double etagate = zoomFactor_ * etasize;
double phigate = zoomFactor_ * phisize;
double eta = maxrh->position().Eta();
double phi = maxrh->position().Phi();
if(displayHist_[EPE]) {
displayHist_[EPE]->GetXaxis()->SetRangeUser(eta-etagate, eta+etagate);
displayHist_[EPE]->GetYaxis()->SetRangeUser(phi-phigate, phi+phigate);
displayView_[EPE]->Modified();
displayView_[EPE]->Update();
}
if(displayHist_[EPH]) {
displayHist_[EPH]->GetXaxis()->SetRangeUser(eta-etagate, eta+etagate);
displayHist_[EPH]->GetYaxis()->SetRangeUser(phi-phigate, phi+phigate);
displayView_[EPH]->Modified();
displayView_[EPH]->Update();
}
if(drawHO_ && displayHist_[EHO]) {
displayHist_[EHO]->GetXaxis()->SetRangeUser(eta-etagate, eta+etagate);
displayHist_[EHO]->GetYaxis()->SetRangeUser(phi-phigate, phi+phigate);
displayView_[EHO]->Modified();
displayView_[EHO]->Update();
}
}
| void DisplayManager::printDisplay | ( | const char * | directory = "" | ) | const |
Definition at line 2231 of file DisplayManager.cc.
References dtNoiseDBValidation_cfg::cerr, gather_cfg::cout, createBeamHaloJobs::directory, displayView_, em_, eventNumber_, mergeVDriftHistosByStation::name, dbtoconf::out, and PFRootEventManager::print().
{
string directory = sdirectory;
if( directory.empty() ) {
directory = "Event_";
}
char num[10];
sprintf(num,"%d", eventNumber_);
directory += num;
string mkdir = "mkdir "; mkdir += directory;
int code = system( mkdir.c_str() );
if( code ) {
cerr<<"cannot create directory "<<directory<<endl;
return;
}
cout<<"Event display printed in directory "<<directory<<endl;
directory += "/";
for(unsigned iView=0; iView<displayView_.size(); iView++) {
if( !displayView_[iView] ) continue;
string name = directory;
name += displayView_[iView]->GetName();
cout<<displayView_[iView]->GetName()<<endl;
string eps = name; eps += ".eps";
displayView_[iView]->SaveAs( eps.c_str() );
string png = name; png += ".png";
displayView_[iView]->SaveAs( png.c_str() );
}
string txt = directory;
txt += "event.txt";
ofstream out( txt.c_str() );
if( !out )
cerr<<"cannot open "<<txt<<endl;
em_->print( out );
}
| void DisplayManager::printGenParticleInfo | ( | std::string | name, |
| int | barcode, | ||
| int | barcodeMother | ||
| ) |
Definition at line 2320 of file DisplayManager.cc.
References gather_cfg::cout, em_, configurableAnalysis::GenParticle, edm::HepMCProduct::GetEvent(), PFRootEventManager::MCTruth_, and AlCaHLTBitMon_ParallelJobs::p.
Referenced by GPFGenParticle::Print().
{
const HepMC::GenEvent* myGenEvent = em_->MCTruth_.GetEvent();
HepMC::GenParticle *p = myGenEvent->barcode_to_particle(barcode);
std::cout<<"genParticle "<<name<<" with barcode "<<barcode<<std::flush<<std::endl;
p->print();
if (barcodeMother) {
HepMC:: GenParticle *mother = myGenEvent->barcode_to_particle(barcodeMother);
std::cout<<"mother particle with barcode "<<barcodeMother<<std::flush<<std::endl;
mother->print();
}
}
| void DisplayManager::readOptions | ( | const char * | file | ) |
Definition at line 65 of file DisplayManager.cc.
References bremAttributes_, bremPatternL_, bremPatternM_, dtNoiseDBValidation_cfg::cerr, clusEnMin_, clusPattern_, clusPatternecal_, clusPatternhfem_, clusPatternho_, clusPSPattern_, clusterAttributes_, gather_cfg::cout, drawBrems_, drawClus_, drawClusterL_, drawGenParticles_, drawGsfTracks_, drawHits_, drawHO_, drawParticles_, drawPFBlocks_, drawTracks_, genParticlePtMin_, genPartPattern_, IO::GetOpt(), gsfAttributes_, gsfPatternL_, gsfPatternM_, gsfPtMin_, hitEnMin_, options_, particlePtMin_, setNewAttrToSimParticles(), simPartPatternDefault_, simPartPatternElec_, simPartPatternK_, simPartPatternL_, simPartPatternM_, simPartPatternMuon_, simPartPatternNeutron_, simPartPatternPhoton_, simPartPatternPi_, simPartPatternProton_, trackAttributes_, trackPatternL_, trackPatternM_, trackPtMin_, viewSize_, viewSizeEtaPhi_, and zoomFactor_.
Referenced by DisplayManager().
{
try {
delete options_;
options_ = new IO(optfile);
}
catch( const string& err ) {
cout<<err<<endl;
return;
}
drawHO_=true;
options_->GetOpt("display", "drawHO", drawHO_);
viewSizeEtaPhi_.clear();
options_->GetOpt("display", "viewsize_etaphi", viewSizeEtaPhi_);
if(viewSizeEtaPhi_.size() != 2) {
cerr<<"PFRootEventManager::ReadOptions, bad display/viewsize_etaphi tag...using 700/350"
<<endl;
viewSizeEtaPhi_.clear();
viewSizeEtaPhi_.push_back(700);
viewSizeEtaPhi_.push_back(350);
}
viewSize_.clear();
options_->GetOpt("display", "viewsize_xy", viewSize_);
if(viewSize_.size() != 2) {
cerr<<"PFRootEventManager::ReadOptions, bad display/viewsize_xy tag...using 700/350"
<<endl;
viewSize_.clear();
viewSize_.push_back(600);
viewSize_.push_back(600);
}
clusterAttributes_.clear();
options_->GetOpt("display", "cluster_attributes", clusterAttributes_);
if(clusterAttributes_.size() != 7) {
cerr<<"PFRootEventManager::::ReadOptions, bad display/cluster_attributes tag...using 20 10 2 5"
<<endl;
clusterAttributes_.clear();
clusterAttributes_.push_back(2); //color
clusterAttributes_.push_back(5); // color if clusterPS
clusterAttributes_.push_back(20); // marker style
clusterAttributes_.push_back(1.0);
clusterAttributes_.push_back(6); //For ECAL
clusterAttributes_.push_back(9); //For HF_EM
clusterAttributes_.push_back(46); //For HO
}
trackAttributes_.clear();
options_->GetOpt("display", "track_attributes", trackAttributes_);
if(trackAttributes_.size() != 4) {
cerr<<"PFRootEventManager::::ReadOptions, bad display/track_attributes tag...using 103 1 8 8"
<<endl;
trackAttributes_.clear();
trackAttributes_.push_back(103); //color Line and Marker
trackAttributes_.push_back(1); //line style
trackAttributes_.push_back(8); //Marker style
trackAttributes_.push_back(1.0); //Marker size
}
gsfAttributes_.clear();
options_->GetOpt("display", "gsf_attributes", gsfAttributes_);
if(gsfAttributes_.size() != 4) {
cerr<<"PFRootEventManager::::ReadOptions, bad display/gsf_attributes tag...using 105 1 8 8"
<<endl;
gsfAttributes_.clear();
gsfAttributes_.push_back(105); //color Line and Marker
gsfAttributes_.push_back(1); //line style
gsfAttributes_.push_back(8); //Marker style
gsfAttributes_.push_back(1.0); //Marker size
}
bremAttributes_.clear();
options_->GetOpt("display", "brem_attributes", bremAttributes_);
if(bremAttributes_.size() != 4) {
cerr<<"PFRootEventManager::::ReadOptions, bad display/gsf_attributes tag...using 106 1 8 8"
<<endl;
bremAttributes_.clear();
bremAttributes_.push_back(106); //color Line and Marker
bremAttributes_.push_back(1); //line style
bremAttributes_.push_back(8); //Marker style
bremAttributes_.push_back(1.0); //Marker size
}
double attrScale = (drawHO_) ? 0.8 : 1.0;
clusPattern_ = new TAttMarker( (int)clusterAttributes_[0],
(int)clusterAttributes_[2],
attrScale*clusterAttributes_[3]);
clusPatternecal_ = new TAttMarker( (int)clusterAttributes_[4],
(int)clusterAttributes_[2],
attrScale*clusterAttributes_[3]);
clusPatternhfem_ = new TAttMarker( (int)clusterAttributes_[5],
(int)clusterAttributes_[2],
attrScale*clusterAttributes_[3]);
clusPatternho_ = new TAttMarker( (int)clusterAttributes_[6],
(int)clusterAttributes_[2],
attrScale*clusterAttributes_[3]);
clusPSPattern_ = new TAttMarker( (int)clusterAttributes_[1],
(int)clusterAttributes_[2],
attrScale*clusterAttributes_[3]);
trackPatternL_ = new TAttLine( (int)trackAttributes_[0],
(int)trackAttributes_[1],
1);
trackPatternM_ = new TAttMarker( (int)trackAttributes_[0],
(int)trackAttributes_[2],
attrScale*trackAttributes_[3]);
gsfPatternL_ = new TAttLine( (int)gsfAttributes_[0],
(int)gsfAttributes_[1],
1);
gsfPatternM_ = new TAttMarker( (int)gsfAttributes_[0],
(int)gsfAttributes_[2],
attrScale*gsfAttributes_[3]);
bremPatternL_ = new TAttLine( (int)bremAttributes_[0],
(int)bremAttributes_[1],
1);
bremPatternM_ = new TAttMarker( (int)bremAttributes_[0],
(int)bremAttributes_[2],
attrScale*bremAttributes_[3]);
genPartPattern_= new TAttMarker(kGreen-1,22,1.);
std::vector<float> simPartAttributes;
simPartAttributes.clear();
options_->GetOpt("display", "simPart_attributes", simPartAttributes);
if(simPartAttributes.size() != 3) {
cerr<<"PFRootEventManager::::ReadOptions, bad display/simPart_attributes tag...using 103 1 8 8"
<<endl;
simPartAttributes.clear();
simPartAttributes.push_back(3); //color Line and Marker
simPartAttributes.push_back(2); //line style
simPartAttributes.push_back(0.6);
}
int simColor = (int)simPartAttributes[0];
int simLStyle = (int)simPartAttributes[1];
float simMSize = attrScale*simPartAttributes[2];
simPartPatternPhoton_ = new TAttMarker(simColor,3,simMSize);
simPartPatternElec_ = new TAttMarker(simColor,5,simMSize);
simPartPatternMuon_ = new TAttMarker(simColor,2,simMSize);
simPartPatternK_ = new TAttMarker(simColor,24,simMSize);
simPartPatternPi_ = new TAttMarker(simColor,25,simMSize);
simPartPatternProton_ = new TAttMarker(simColor,26,simMSize);
simPartPatternNeutron_= new TAttMarker(simColor,27,simMSize);
simPartPatternDefault_= new TAttMarker(simColor,30,simMSize);
simPartPatternL_ = new TAttLine(simColor,simLStyle,1);
simPartPatternM_.resize(8);
setNewAttrToSimParticles();
drawHits_= true;
options_->GetOpt("display", "rechits",drawHits_);
drawClus_ = true;
options_->GetOpt("display", "clusters",drawClus_);
drawClusterL_ = false;
options_->GetOpt("display", "cluster_lines", drawClusterL_);
drawTracks_ = true;
options_->GetOpt("display", "rectracks", drawTracks_);
drawGsfTracks_ = true;
options_->GetOpt("display", "gsftracks", drawGsfTracks_);
drawBrems_ = false;
options_->GetOpt("display", "brems", drawBrems_);
drawParticles_ = true;
options_->GetOpt("display", "particles", drawParticles_);
particlePtMin_ = -1;
options_->GetOpt("display", "particles_ptmin", particlePtMin_);
drawGenParticles_=false;
genParticlePtMin_ = 0;
trackPtMin_ = -1;
options_->GetOpt("display", "rectracks_ptmin", trackPtMin_);
gsfPtMin_ = -1;
options_->GetOpt("display", "gsfrectracks_ptmin", gsfPtMin_);
hitEnMin_ = -1;
options_->GetOpt("display","rechits_enmin",hitEnMin_);
clusEnMin_ = -1;
options_->GetOpt("display","clusters_enmin",clusEnMin_);
drawPFBlocks_ = false;
options_->GetOpt("display","drawPFBlock",drawPFBlocks_);
//redrawWithoutHits_=false;
zoomFactor_ = 10;
options_->GetOpt("display", "zoom_factor", zoomFactor_);
}
| void DisplayManager::reset | ( | void | ) |
Definition at line 2276 of file DisplayManager.cc.
References badBremsId_, blockIdentsMap_, graphicMap_, isGraphicLoaded_, maxERecHitEcal_, maxERecHitHcal_, maxERecHitHo_, AlCaHLTBitMon_ParallelJobs::p, and selectedGObj_.
Referenced by display(), displayEvent(), displayNextInteresting(), and ~DisplayManager().
{
maxERecHitEcal_=-1;
maxERecHitHcal_=-1;
maxERecHitHo_=-1;
isGraphicLoaded_= false;
std::multimap<int,GPFBase *>::iterator p;
for (p=graphicMap_.begin();p!=graphicMap_.end();p++)
delete p->second;
graphicMap_.clear();
blockIdentsMap_.clear();
selectedGObj_.clear();
badBremsId_.clear();
}
| void DisplayManager::retrieveBadBrems | ( | ) | [private] |
Definition at line 1667 of file DisplayManager.cc.
References badBremsId_, begin, reco::PFBlockElement::BREM, BREMID, ECAL, em_, end, reco::PFBlockElementBrem::GsftrackRefPF(), reco::PFBlockElementBrem::indTrajPoint(), edm::Ref< C, T, F >::key(), reco::PFBlock::LINKTEST_ALL, PFRootEventManager::pfBlocks_, shiftId_, and findQualityFiles::size.
Referenced by loadGGsfRecTracks().
{
//selects Brems with no information in PFBlock.Those selected Brems are not displayed
int size = em_->pfBlocks_->size();
for (int ibl=0;ibl<size;ibl++) {
edm::OwnVector< reco::PFBlockElement >::const_iterator iter;
for( iter =((*(em_->pfBlocks_))[ibl].elements()).begin();
iter != ((*(em_->pfBlocks_))[ibl].elements()).end();iter++) {
int ident=-1;
reco::PFBlockElement::Type type = (*iter).type();
if (type == reco::PFBlockElement::BREM) {
std::multimap<double, unsigned> ecalElems;
(*(em_->pfBlocks_))[ibl].associatedElements( (*iter).index(),(*(em_->pfBlocks_))[ibl].linkData(),
ecalElems ,
reco::PFBlockElement::ECAL,
reco::PFBlock::LINKTEST_ALL );
if (ecalElems.size()==0) {
// std::cout<<" PfBlock Nb "<<ibl<<" -- brem elem "<<(*iter).index()<<"-type "<<(*iter).type()<<" not drawn"<<std::flush<<std::endl;
const reco::PFBlockElementBrem * Brem = dynamic_cast<const reco::PFBlockElementBrem*>(&(*iter));
reco::GsfPFRecTrackRef trackref = Brem->GsftrackRefPF();
unsigned ind=trackref.key()*40+Brem->indTrajPoint();
ident = (BREMID << shiftId_ ) | ind ;
badBremsId_.push_back(ident);
}
}
}
}
}
| void DisplayManager::rubOutGPFBlock | ( | ) | [private] |
Definition at line 1289 of file DisplayManager.cc.
References create_public_lumi_plots::color, drawGObject(), i, selectedGObj_, and findQualityFiles::size.
Referenced by displayPFBlock(), and findAndDraw().
{
int size = selectedGObj_.size();
bool toInitial=true;
int color=0;
for (int i=0;i<size;i++)
drawGObject(selectedGObj_[i],color,toInitial);
}
| void DisplayManager::setNewAttrToSimParticles | ( | ) |
Definition at line 2306 of file DisplayManager.cc.
References simPartPatternDefault_, simPartPatternElec_, simPartPatternK_, simPartPatternM_, simPartPatternMuon_, simPartPatternNeutron_, simPartPatternPhoton_, simPartPatternPi_, and simPartPatternProton_.
Referenced by readOptions().
{
simPartPatternM_.clear();
simPartPatternM_.push_back(simPartPatternPhoton_);
simPartPatternM_.push_back(simPartPatternElec_);
simPartPatternM_.push_back(simPartPatternMuon_);
simPartPatternM_.push_back(simPartPatternK_);
simPartPatternM_.push_back(simPartPatternPi_);
simPartPatternM_.push_back(simPartPatternProton_);
simPartPatternM_.push_back(simPartPatternNeutron_);
simPartPatternM_.push_back(simPartPatternDefault_);
}
| void DisplayManager::unZoom | ( | ) |
Definition at line 2295 of file DisplayManager.cc.
References displayHist_, i, and updateDisplay().
Referenced by DialogFrame::unZoom().
{
for( unsigned i=0; i<displayHist_.size(); i++) {
// the corresponding view was not requested
if( ! displayHist_[i] ) continue;
displayHist_[i]->GetXaxis()->UnZoom();
displayHist_[i]->GetYaxis()->UnZoom();
}
updateDisplay();
}
| void DisplayManager::updateDisplay | ( | ) |
Definition at line 1409 of file DisplayManager.cc.
References displayView_, drawHO_, EHO, i, and AlCaHLTBitMon_ParallelJobs::p.
Referenced by findAndDraw(), findBlock(), unZoom(), and DialogFrame::updateDisplay().
{
for(unsigned i=0; i<displayView_.size(); i++) {
if (!drawHO_ && i==EHO) continue;
TPad* p = displayView_[i];
assert( p );
p->Modified();
p->Update();
}
}
std::vector<int> DisplayManager::badBremsId_ [private] |
Definition at line 185 of file DisplayManager.h.
Referenced by findBadBremsId(), reset(), and retrieveBadBrems().
| std::multimap<int , std::pair <int,int> > DisplayManager::blockIdentsMap_ |
Definition at line 95 of file DisplayManager.h.
Referenced by displayPFBlock(), findBlock(), loadGPFBlocks(), and reset().
| std::vector<float> DisplayManager::bremAttributes_ |
Definition at line 100 of file DisplayManager.h.
Referenced by readOptions().
| TAttLine* DisplayManager::bremPatternL_ |
Definition at line 116 of file DisplayManager.h.
Referenced by createGTrack(), and readOptions().
| TAttMarker* DisplayManager::bremPatternM_ |
Definition at line 115 of file DisplayManager.h.
Referenced by createGTrack(), and readOptions().
| double DisplayManager::clusEnMin_ |
Definition at line 72 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyPtThreshold(), and readOptions().
| TAttMarker* DisplayManager::clusPattern_ |
Definition at line 104 of file DisplayManager.h.
Referenced by DialogFrame::createCanvasAttr(), createGCluster(), DialogFrame::modifyAttr(), and readOptions().
| TAttMarker* DisplayManager::clusPatternecal_ |
Definition at line 107 of file DisplayManager.h.
Referenced by createGCluster(), and readOptions().
| TAttMarker* DisplayManager::clusPatternhfem_ |
Definition at line 109 of file DisplayManager.h.
Referenced by createGCluster(), and readOptions().
| TAttMarker* DisplayManager::clusPatternho_ |
Definition at line 108 of file DisplayManager.h.
Referenced by createGCluster(), and readOptions().
| TAttMarker* DisplayManager::clusPSPattern_ |
Definition at line 105 of file DisplayManager.h.
Referenced by createGCluster(), and readOptions().
| std::vector<float> DisplayManager::clusterAttributes_ |
Definition at line 101 of file DisplayManager.h.
Referenced by readOptions().
std::vector<TH2F*> DisplayManager::displayHist_ [private] |
support histogram for x/y or r/z display.
Definition at line 163 of file DisplayManager.h.
Referenced by createCanvas(), displayCanvas(), lookForGenParticle(), lookForMaxRecHit(), and unZoom().
std::vector<TCanvas*> DisplayManager::displayView_ [private] |
vector of canvas for x/y or r/z display
Definition at line 154 of file DisplayManager.h.
Referenced by createCanvas(), displayAll(), displayCanvas(), drawGObject(), lookForGenParticle(), lookForMaxRecHit(), printDisplay(), and updateDisplay().
Definition at line 82 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), enableDrawBrem(), and readOptions().
Definition at line 83 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyOptions(), and readOptions().
Definition at line 84 of file DisplayManager.h.
Referenced by readOptions().
Definition at line 86 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyOptions(), and readOptions().
Definition at line 81 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyOptions(), and readOptions().
Definition at line 79 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyOptions(), and readOptions().
Definition at line 88 of file DisplayManager.h.
Referenced by createCanvas(), createGCluster(), createGGenParticle(), createGPart(), createGRecHit(), createGTrack(), displayAll(), displayCanvas(), drawGObject(), drawWithNewGraphicAttributes(), findAndDraw(), loadGClusters(), loadGPFBlocks(), loadGRecHits(), lookForGenParticle(), lookForMaxRecHit(), readOptions(), and updateDisplay().
Definition at line 85 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyOptions(), and readOptions().
Definition at line 87 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayPFBlock(), enableDrawPFBlock(), findAndDraw(), and readOptions().
Definition at line 80 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyOptions(), and readOptions().
PFRootEventManager* DisplayManager::em_ [private] |
Definition at line 134 of file DisplayManager.h.
Referenced by createGGenParticle(), createGRecHit(), display(), displayEvent(), DisplayManager(), displayNext(), displayNextInteresting(), displayPrevious(), findBlock(), getMaxE(), loadGClusters(), loadGGenParticles(), loadGGsfRecTracks(), loadGPFBlocks(), loadGRecHits(), loadGRecTracks(), loadGSimParticles(), lookForGenParticle(), lookForMaxRecHit(), printDisplay(), printGenParticleInfo(), and retrieveBadBrems().
int DisplayManager::eventNumber_ [private] |
Definition at line 143 of file DisplayManager.h.
Referenced by display(), displayEvent(), DisplayManager(), displayNext(), displayNextInteresting(), displayPrevious(), and printDisplay().
TBox DisplayManager::frontFaceECALRZ_ [private] |
ECAL in RZ view.
Definition at line 169 of file DisplayManager.h.
Referenced by createCanvas(), and displayCanvas().
TEllipse DisplayManager::frontFaceECALXY_ [private] |
ECAL in XY view.
Definition at line 166 of file DisplayManager.h.
Referenced by createCanvas(), and displayCanvas().
TEllipse DisplayManager::frontFaceHCALXY_ [private] |
HCAL in XY view.
Definition at line 172 of file DisplayManager.h.
Referenced by createCanvas(), and displayCanvas().
TEllipse DisplayManager::frontFaceHOXY_ [private] |
HO in XY view.
Definition at line 175 of file DisplayManager.h.
Referenced by createCanvas(), and displayCanvas().
Definition at line 77 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyPtThreshold(), and readOptions().
| TAttMarker* DisplayManager::genPartPattern_ |
Definition at line 117 of file DisplayManager.h.
Referenced by createGGenParticle(), and readOptions().
| std::multimap<int,GPFBase *> DisplayManager::graphicMap_ |
Definition at line 93 of file DisplayManager.h.
Referenced by createGCluster(), createGGenParticle(), createGPart(), createGRecHit(), createGTrack(), displayAll(), drawGObject(), drawWithNewGraphicAttributes(), and reset().
| std::vector<float> DisplayManager::gsfAttributes_ |
Definition at line 99 of file DisplayManager.h.
Referenced by readOptions().
| TAttLine* DisplayManager::gsfPatternL_ |
Definition at line 114 of file DisplayManager.h.
Referenced by createGTrack(), and readOptions().
| TAttMarker* DisplayManager::gsfPatternM_ |
Definition at line 113 of file DisplayManager.h.
Referenced by createGTrack(), and readOptions().
| double DisplayManager::gsfPtMin_ |
Definition at line 75 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyPtThreshold(), and readOptions().
| double DisplayManager::hitEnMin_ |
Definition at line 73 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyPtThreshold(), and readOptions().
bool DisplayManager::isGraphicLoaded_ [private] |
Definition at line 142 of file DisplayManager.h.
Referenced by display(), displayAll(), displayEvent(), displayNextInteresting(), and reset().
double DisplayManager::maxERecHitEcal_ [private] |
Definition at line 139 of file DisplayManager.h.
Referenced by getMaxEEcal(), and reset().
double DisplayManager::maxERecHitHcal_ [private] |
Definition at line 140 of file DisplayManager.h.
Referenced by getMaxEHcal(), and reset().
double DisplayManager::maxERecHitHo_ [private] |
Definition at line 141 of file DisplayManager.h.
Referenced by getMaxEHo(), and reset().
int DisplayManager::maxEvents_ [private] |
Definition at line 144 of file DisplayManager.h.
Referenced by display(), and DisplayManager().
IO* DisplayManager::options_ [private] |
Definition at line 76 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyPtThreshold(), and readOptions().
std::vector<int> DisplayManager::selectedGObj_ [private] |
graphic object containers
Definition at line 184 of file DisplayManager.h.
Referenced by displayPFBlock(), drawGObject(), findAndDraw(), reset(), and rubOutGPFBlock().
int DisplayManager::shiftId_ [private] |
shiftId_ is the number of low bits indicating the object nb in the map key the object type is obtained by shifting the key number of shiftId_ bits ShiftId_ == SHIFTID == 26
Definition at line 150 of file DisplayManager.h.
Referenced by createGGenParticle(), displayAll(), drawWithNewGraphicAttributes(), findAndDraw(), loadGClusters(), loadGGsfRecTracks(), loadGPFBlocks(), loadGRecHits(), loadGRecTracks(), loadGSimParticles(), and retrieveBadBrems().
| TAttMarker* DisplayManager::simPartPatternDefault_ |
Definition at line 127 of file DisplayManager.h.
Referenced by readOptions(), and setNewAttrToSimParticles().
| TAttMarker* DisplayManager::simPartPatternElec_ |
Definition at line 121 of file DisplayManager.h.
Referenced by readOptions(), and setNewAttrToSimParticles().
| TAttMarker* DisplayManager::simPartPatternK_ |
Definition at line 123 of file DisplayManager.h.
Referenced by readOptions(), and setNewAttrToSimParticles().
| TAttLine* DisplayManager::simPartPatternL_ |
Definition at line 119 of file DisplayManager.h.
Referenced by createGPart(), and readOptions().
| std::vector<TAttMarker *> DisplayManager::simPartPatternM_ |
Definition at line 129 of file DisplayManager.h.
Referenced by createGPart(), readOptions(), and setNewAttrToSimParticles().
| TAttMarker* DisplayManager::simPartPatternMuon_ |
Definition at line 122 of file DisplayManager.h.
Referenced by readOptions(), and setNewAttrToSimParticles().
| TAttMarker* DisplayManager::simPartPatternNeutron_ |
Definition at line 126 of file DisplayManager.h.
Referenced by readOptions(), and setNewAttrToSimParticles().
| TAttMarker* DisplayManager::simPartPatternPhoton_ |
Definition at line 120 of file DisplayManager.h.
Referenced by readOptions(), and setNewAttrToSimParticles().
| TAttMarker* DisplayManager::simPartPatternPi_ |
Definition at line 124 of file DisplayManager.h.
Referenced by readOptions(), and setNewAttrToSimParticles().
| TAttMarker* DisplayManager::simPartPatternProton_ |
Definition at line 125 of file DisplayManager.h.
Referenced by readOptions(), and setNewAttrToSimParticles().
| std::vector<float> DisplayManager::trackAttributes_ |
Definition at line 98 of file DisplayManager.h.
Referenced by readOptions().
| TAttLine* DisplayManager::trackPatternL_ |
Definition at line 112 of file DisplayManager.h.
Referenced by DialogFrame::createCanvasAttr(), createGTrack(), DialogFrame::modifyAttr(), and readOptions().
| TAttMarker* DisplayManager::trackPatternM_ |
Definition at line 111 of file DisplayManager.h.
Referenced by DialogFrame::createCanvasAttr(), createGTrack(), DialogFrame::modifyAttr(), and readOptions().
| double DisplayManager::trackPtMin_ |
Definition at line 74 of file DisplayManager.h.
Referenced by DialogFrame::createCmdFrame(), displayAll(), DialogFrame::doModifyPtThreshold(), and readOptions().
std::vector<int> DisplayManager::viewSize_ [private] |
display pad xy size for (x,y) or (r,z) display
Definition at line 157 of file DisplayManager.h.
Referenced by createCanvas(), and readOptions().
std::vector<int> DisplayManager::viewSizeEtaPhi_ [private] |
display pad xy size for eta/phi view
Definition at line 160 of file DisplayManager.h.
Referenced by readOptions().
double DisplayManager::zoomFactor_ [private] |
Definition at line 145 of file DisplayManager.h.
Referenced by lookForGenParticle(), lookForMaxRecHit(), and readOptions().
1.7.3