|
|
|
#include <MEGeom.h>
Public Member Functions | |
| virtual | ~MEGeom () |
Static Public Member Functions | |
| static void | drawGlobalBoundaries (int lineColor) |
| static void | drawHist (int ilmr, int unit, TCanvas *canv=0) |
| static TGraph * | getBoundary (int ilmr, int unit) |
| static TH2 * | getGlobalHist (const char *name=0) |
| static TH2 * | getHist (int ilmr, int unit) |
| static void | setBinGlobalHist (TH2 *h, int ix, int iy, int iz, float val) |
Static Private Attributes | |
| static TH2 * | _h |
| static int | _nbinx = 2*( MEGeom::_nbuf + 100 + 85 ) + 1 |
| static int | _nbiny = 2*( MEGeom::_nbuf + 180 ) |
| static int | _nbuf = 5 |
| static float | _xmax = 0.5 + MEGeom::_nbuf + 100 + 85 |
| static float | _xmin = -( 0.5 + MEGeom::_nbuf + 100 + 85 ) |
| static float | _ymax = 180.5 + ( MEGeom::_nbuf + 180 ) |
| static float | _ymin = 180.5 - ( MEGeom::_nbuf + 180 ) |
| void MEGeom::drawGlobalBoundaries | ( | int | lineColor | ) | [static] |
Definition at line 254 of file MEGeom.cc.
References MEEBGeom::getGraphBoundary(), MEEEGeom::getGraphBoundary(), MEEEGeom::iSector, ecalpyutils::ism(), MEEBGeom::iSuperModule, and n.
{
TGraph* gr(0);
for( int ism=1; ism<=36; ism++ )
{
gr = MEEBGeom::getGraphBoundary( MEEBGeom::iSuperModule, ism, true );
gr->SetLineWidth( 2 );
gr->SetLineColor( lineColor );
gr->Draw("LSame");
}
for( int isec=1; isec<=9; isec++ )
{
gr = MEEEGeom::getGraphBoundary( MEEEGeom::iSector, isec );
TGraph* grm = (TGraph*) gr->Clone();
TGraph* grp = (TGraph*) gr->Clone();
int n = gr->GetN();
// std::cout << std::endl;
for( int ii=0; ii<n; ii++ )
{
double x_, y_;
gr->GetPoint( ii, x_, y_ );
// std::cout << isec << " " << ii << " x=" << x_ << " y=" << y_ << std::endl;
double xx_, yy_;
xx_ = -85-x_;
yy_ = 231-y_;
// std::cout << isec << " " << 0 << " x=" << xx_ << " y=" << yy_ << std::endl;
grm->SetPoint( ii, xx_, yy_ );
xx_ = 85+x_;
yy_ = 231-y_;
// std::cout << isec << " " << 1 << " x=" << xx_ << " y=" << yy_ << std::endl;
grp->SetPoint( ii, xx_, yy_ );
}
grm->SetLineColor( lineColor );
grm->SetLineWidth( 2 );
grm->Draw("LSame");
grp->SetLineColor( lineColor );
grp->SetLineWidth( 2 );
grp->Draw("LSame");
}
}
| void MEGeom::drawHist | ( | int | ilmr, |
| int | unit, | ||
| TCanvas * | canv = 0 |
||
| ) | [static] |
Definition at line 167 of file MEGeom.cc.
References h, ME::iCrystal, ME::iElectronicChannel, ME::iHVChannel, ME::iLMModule, ME::iLMRegion, ME::iLVChannel, ME::iSector, and ME::iSuperCrystal.
{
TH2* h = getHist( ilmr, histtype );
assert( h!=0 );
TString tname = h->GetTitle();
switch( histtype )
{
case ME::iSector: break;
case ME::iLMRegion: tname += " Monitoring Regions"; break;
case ME::iLMModule: tname += " Monitoring Modules"; break;
case ME::iSuperCrystal: tname += " Super Crystals"; break;
case ME::iCrystal: tname += " Crystals"; break;
case ME::iElectronicChannel: tname += " Electronic Channels"; break;
case ME::iHVChannel: tname += " HV Channels"; break;
case ME::iLVChannel: tname += " LV Channels"; break;
}
if( canv==0 )
{
TString cname = tname;
cname.ReplaceAll(" ","_");
canv = new TCanvas( cname, cname, 10, 10, 500, 500 );
}
canv->SetTitle( tname );
canv->cd();
h->Draw("COLZ");
TGraph* gsect = getBoundary( ilmr, ME::iSector );
assert( gsect!=0 );
gsect->SetLineWidth( 1 );
gsect->Draw("LSame");
TGraph* gside = getBoundary( ilmr, ME::iLMRegion );
assert( gside!=0 );
gside->SetLineWidth( 2 );
gside->Draw("LSame");
}
| TGraph * MEGeom::getBoundary | ( | int | ilmr, |
| int | unit | ||
| ) | [static] |
Definition at line 115 of file MEGeom.cc.
References ME::iEBM, ME::iEBP, ME::iEEM, ME::iEEP, ME::iLMRegion, ME::iSector, ecalpyutils::ism(), ME::regionAndSector(), and reco_application_2006_simpleTBanalysis_cfg::rootfile.
{
// for local pictures, only sector or monitoring region
if( histtype!=ME::iSector && histtype!=ME::iLMRegion )
histtype = ME::iSector;
int ireg;
int ism;
int idcc;
int side;
ME::regionAndSector( ilmr, ireg, ism, idcc, side );
//
// to produce these root files, run the runGeom executable
//
TGraph* g_(0);
TFile* rootfile(0);
TString gn_;
if( ireg==ME::iEBM || ireg==ME::iEBP )
{
rootfile = TFile::Open("ebgeom.root");
assert( rootfile!=0 );
switch (histtype)
{
case ME::iSector: gn_ = "SuperModule"; break;
case ME::iLMRegion: gn_ = "Side_"; gn_+=side; break;
}
}
else
{
int isect=ism;
if( ireg==ME::iEEM )
{
isect-=9;
rootfile = TFile::Open("eegeom_1.root");
}
if( ireg==ME::iEEP ) rootfile = TFile::Open("eegeom_2.root");
assert( rootfile!=0 );
int lmr_= ilmr;
if( ireg==ME::iEEP ) lmr_-=72;
else if( ireg==ME::iEEM ) lmr_-=82;
switch (histtype)
{
case ME::iSector: gn_ = "Sector_"; gn_+=isect; break;
case ME::iLMRegion: gn_ = "LMRegion_"; gn_+=lmr_; break;
}
}
g_ = (TGraph*) rootfile->Get( gn_ );
return g_;
}
| TH2 * MEGeom::getGlobalHist | ( | const char * | name = 0 | ) | [static] |
| TH2 * MEGeom::getHist | ( | int | ilmr, |
| int | unit | ||
| ) | [static] |
Definition at line 34 of file MEGeom.cc.
References gather_cfg::cout, ME::iCrystal, ME::iEBM, ME::iEBP, ME::iEEM, ME::iEEP, ME::iElectronicChannel, ME::iHVChannel, ME::iLMModule, ME::iLMRegion, ME::iLVChannel, ME::iSector, ecalpyutils::ism(), ME::iSuperCrystal, ME::regionAndSector(), reco_application_2006_simpleTBanalysis_cfg::rootfile, MEEEGeom::smName(), and MEEBGeom::smName().
{
int ireg;
int ism;
int idcc;
int side;
ME::regionAndSector( ilmr, ireg, ism, idcc, side );
std::cout << "LM region=" << ilmr << " reg/sm/dcc/side " << ireg << "/" << ism << "/" << idcc << "/" << side << std::endl;
TString hname = "LMR="; hname += ilmr; hname += " ";
if( ireg==ME::iEEM || ireg==ME::iEEP )
{
hname += MEEEGeom::smName( ism );
}
else if( ireg==ME::iEBM || ireg==ME::iEBP )
{
hname += MEEBGeom::smName( ism );
}
else
abort();
hname += " ";
hname += "DCC="; hname += idcc; hname += "/" ; hname += side;
//
// to produce these root files, run the runGeom executable
//
TH2* h_(0);
TFile* rootfile(0);
TString hn_;
if( ireg==ME::iEBM || ireg==ME::iEBP )
{
rootfile = TFile::Open("ebgeom.root");
assert( rootfile!=0 );
hn_="eb_loc";
switch ( unit )
{
case ME::iSector: break;
case ME::iLMRegion: hn_ += "_side"; break;
case ME::iLMModule: hn_ += "_lmmod"; break;
case ME::iSuperCrystal: hn_ += "_tt"; break;
case ME::iCrystal: hn_ += "_cr"; break;
case ME::iElectronicChannel: hn_ += "_elecr"; break;
case ME::iHVChannel: hn_ += "_hv"; break;
case ME::iLVChannel: hn_ += "_lv"; break;
}
}
else
{
int isect=ism;
if( ireg==ME::iEEM )
{
isect-=9;
rootfile = TFile::Open("eegeom_1.root");
}
if( ireg==ME::iEEP ) rootfile = TFile::Open("eegeom_2.root");
assert( rootfile!=0 );
hn_="eem_S"; hn_+= isect;
switch (unit)
{
case ME::iSector: break;
case ME::iLMRegion: break;
case ME::iLMModule: hn_ += "_lmmod"; break;
case ME::iSuperCrystal: hn_ += "_sc"; break;
case ME::iCrystal: hn_ += "_cr"; break;
case ME::iElectronicChannel: hn_ += "_cr"; break;
case ME::iHVChannel: break;
case ME::iLVChannel: break;
}
}
h_ = (TH2*) rootfile->Get( hn_ );
h_->SetTitle( hname );
h_->GetXaxis()->SetTitle("ix");
h_->GetXaxis()->CenterTitle();
h_->GetYaxis()->SetTitle("iy");
h_->GetYaxis()->CenterTitle();
return h_;
}
| void MEGeom::setBinGlobalHist | ( | TH2 * | h, |
| int | ix, | ||
| int | iy, | ||
| int | iz, | ||
| float | val | ||
| ) | [static] |
Definition at line 215 of file MEGeom.cc.
References abs.
{
// make sure it's a global hist
TAxis* ax = h->GetXaxis();
TAxis* ay = h->GetYaxis();
assert( ax->GetNbins()==_nbinx && ay->GetNbins()==_nbiny );
assert( ax->GetXmax()==_xmax && ax->GetXmin()==_xmin );
assert( ay->GetXmax()==_ymax && ay->GetXmin()==_ymin );
int ibinx(0);
int ibiny(0);
if( iz==0 )
{
int ieta = ix;
int iphi = iy;
assert( abs(ieta)>=1 && abs(ieta)<=85 );
assert( iphi>=1 && iphi<=360 );
ibinx = ax->FindBin( ieta );
ibiny = ay->FindBin( iphi );
}
else if( iz==-1 )
{
assert( ix>=1 && ix<=100 && iy>=1 && iy<=100 );
ibinx = (_nbuf+100+1)-ix;
ibiny = (_nbuf+180+50+1)-iy;
}
else if( iz==1 )
{
assert( ix>=1 && ix<=100 && iy>=1 && iy<=100 );
ibinx = _nbinx+1-((_nbuf+100+1)-ix);
ibiny = (_nbuf+180+50+1)-iy;
}
h->SetBinContent( ibinx, ibiny, val );
}
TH2 * MEGeom::_h [static, private] |
new TH2F( "globalEcal", "Global representation of ECAL", MEGeom::_nbinx, MEGeom::_xmin, MEGeom::_xmax, MEGeom::_nbiny, MEGeom::_ymin, MEGeom::_ymax )
int MEGeom::_nbinx = 2*( MEGeom::_nbuf + 100 + 85 ) + 1 [static, private] |
int MEGeom::_nbiny = 2*( MEGeom::_nbuf + 180 ) [static, private] |
int MEGeom::_nbuf = 5 [static, private] |
float MEGeom::_xmax = 0.5 + MEGeom::_nbuf + 100 + 85 [static, private] |
float MEGeom::_xmin = -( 0.5 + MEGeom::_nbuf + 100 + 85 ) [static, private] |
float MEGeom::_ymax = 180.5 + ( MEGeom::_nbuf + 180 ) [static, private] |
float MEGeom::_ymin = 180.5 - ( MEGeom::_nbuf + 180 ) [static, private] |
1.7.3