|
|
|
#include <SiPixelErrorEstimation.h>
Definition at line 60 of file SiPixelErrorEstimation.h.
| SiPixelErrorEstimation::SiPixelErrorEstimation | ( | const edm::ParameterSet & | ps | ) | [explicit] |
Definition at line 30 of file SiPixelErrorEstimation.cc.
References checkType_, genType_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), include_trk_hits_, outputFile_, and src_.
:tfile_(0), ttree_all_hits_(0), ttree_track_hits_(0) { //Read config file outputFile_ = ps.getUntrackedParameter<string>( "outputFile", "SiPixelErrorEstimation_Ntuple.root" ); // Replace "ctfWithMaterialTracks" with "generalTracks" //src_ = ps.getUntrackedParameter<std::string>( "src", "ctfWithMaterialTracks" ); src_ = ps.getUntrackedParameter<std::string>( "src", "generalTracks" ); checkType_ = ps.getParameter<bool>( "checkType" ); genType_ = ps.getParameter<int>( "genType" ); include_trk_hits_ = ps.getParameter<bool>( "include_trk_hits" ); }
| SiPixelErrorEstimation::~SiPixelErrorEstimation | ( | ) | [virtual] |
Definition at line 44 of file SiPixelErrorEstimation.cc.
{}
| void SiPixelErrorEstimation::analyze | ( | const edm::Event & | e, |
| const edm::EventSetup & | es | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 261 of file SiPixelErrorEstimation.cc.
References abs, SiPixelCluster::Pixel::adc, all_alpha, all_beta, all_bigx, all_bigy, all_blade, all_clust_alpha, all_clust_beta, all_clust_col, all_clust_geoid, all_clust_maxpixcol, all_clust_maxpixrow, all_clust_minpixcol, all_clust_minpixrow, all_clust_q, all_clust_row, all_clust_x, all_clust_y, all_col1, all_col2, all_cols, all_disk, all_edgex, all_edgey, all_eloss, all_flipped, all_gx1, all_gx2, all_gy1, all_gy2, all_gz1, all_gz2, all_half, all_hit_cprob0, all_hit_cprob1, all_hit_cprob2, all_hit_probx, all_hit_proby, all_ladder, all_layer, all_mod, all_npix, all_nsimhit, all_nxpix, all_nypix, all_panel, all_pidhit, all_pixadc, all_pixcol, all_pixgx, all_pixgy, all_pixgz, all_pixrow, all_pixx, all_pixy, all_plaq, all_rechiterrx, all_rechiterry, all_rechitpullx, all_rechitpully, all_rechitresx, all_rechitresy, all_rechitx, all_rechity, all_rechitz, all_row1, all_row2, all_rows, all_side, all_simhitx, all_simhity, all_simphi, all_simproc, all_simpx, all_simpy, all_simpz, all_simtheta, all_simtrketa, all_simtrkphi, all_subdetid, all_trkid, all_vtxr, all_vtxz, all_x1, all_x2, all_y1, all_y2, all_z1, all_z2, alpha, TrackerHitAssociator::associateHit(), edmNew::DetSet< T >::begin(), beta, bigx, bigy, PXFDetId::blade(), blade, charge, checkType_, PixelTopology::containsBigPixelInX(), PixelTopology::containsBigPixelInY(), gather_cfg::cout, PXFDetId::disk(), disk, edgex, edgey, edmNew::DetSet< T >::end(), error, eta, edm::EventID::event(), evt, flipped, genType_, edm::EventSetup::get(), edm::Event::getByLabel(), half, hit_cprob0, hit_cprob1, hit_cprob2, hit_probx, hit_proby, i, edm::EventBase::id(), include_trk_hits_, PixelTopology::isItEdgePixelInX(), PixelTopology::isItEdgePixelInY(), ladder, PXBDetId::ladder(), PXBDetId::layer(), layer, asciidump::le, Topology::localPosition(), m, math_pi, Topology::measurementPosition(), mod, PXFDetId::module(), PXBDetId::module(), PixelTopology::ncolumns(), npix, PixelTopology::nrows(), nsimhit, nxpix, nypix, PXFDetId::panel(), panel, PV3DBase< T, PVType, FrameType >::perp(), phi, evf::utils::pid, pidhit, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, plaq, position, GloballyPositioned< T >::position(), edm::Handle< T >::product(), rechiterrx, rechiterry, rechitpullx, rechitpully, rechitresx, rechitresy, rechitx, rechity, rechitz, edm::EventID::run(), run, PXFDetId::side(), side, simhitx, simhity, simproc, PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), src_, subdetId, DetId::subdetId(), GeomDet::surface(), Surface::toGlobal(), patCandidatesForDimuonsSequences_cff::tracker, testEve_cfg::tracks, trk_alpha, trk_beta, ttree_all_hits_, ttree_track_hits_, SiPixelCluster::Pixel::x, PV3DBase< T, PVType, FrameType >::x(), PV2DBase< T, PVType, FrameType >::x(), LocalError::xx(), SiPixelCluster::Pixel::y, PV3DBase< T, PVType, FrameType >::y(), PV2DBase< T, PVType, FrameType >::y(), LocalError::yy(), and PV3DBase< T, PVType, FrameType >::z().
{
using namespace edm;
run = e.id().run();
evt = e.id().event();
if ( evt%1000 == 0 )
cout << "evt = " << evt << endl;
float math_pi = 3.14159265;
float radtodeg = 180.0 / math_pi;
DetId detId;
LocalPoint position;
LocalError error;
float mindist = 999999.9;
std::vector<PSimHit> matched;
TrackerHitAssociator associate(e);
edm::ESHandle<TrackerGeometry> pDD;
es.get<TrackerDigiGeometryRecord> ().get (pDD);
const TrackerGeometry* tracker = &(* pDD);
// --------------------------------------- all hits -----------------------------------------------------------
edm::Handle<SiPixelRecHitCollection> recHitColl;
e.getByLabel( "siPixelRecHits", recHitColl);
Handle<edm::SimTrackContainer> simtracks;
e.getByLabel("g4SimHits", simtracks);
//-----Iterate over detunits
for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++)
{
DetId detId = ((*it)->geographicalId());
SiPixelRecHitCollection::const_iterator dsmatch = recHitColl->find(detId);
if (dsmatch == recHitColl->end()) continue;
SiPixelRecHitCollection::DetSet pixelrechitRange = *dsmatch;
SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorBegin = pixelrechitRange.begin();
SiPixelRecHitCollection::DetSet::const_iterator pixelrechitRangeIteratorEnd = pixelrechitRange.end();
SiPixelRecHitCollection::DetSet::const_iterator pixeliter = pixelrechitRangeIteratorBegin;
std::vector<PSimHit> matched;
//----Loop over rechits for this detId
for ( ; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter)
{
matched.clear();
matched = associate.associateHit(*pixeliter);
// only consider rechits that have associated simhit
// otherwise cannot determine residiual
if ( matched.empty() )
{
cout << "SiPixelErrorEstimation::analyze: rechits without associated simhit !!!!!!!"
<< endl;
continue;
}
all_subdetid = -9999;
all_layer = -9999;
all_ladder = -9999;
all_mod = -9999;
all_side = -9999;
all_disk = -9999;
all_blade = -9999;
all_panel = -9999;
all_plaq = -9999;
all_half = -9999;
all_flipped = -9999;
all_cols = -9999;
all_rows = -9999;
all_rechitx = -9999;
all_rechity = -9999;
all_rechitz = -9999;
all_simhitx = -9999;
all_simhity = -9999;
all_rechiterrx = -9999;
all_rechiterry = -9999;
all_rechitresx = -9999;
all_rechitresy = -9999;
all_rechitpullx = -9999;
all_rechitpully = -9999;
all_npix = -9999;
all_nxpix = -9999;
all_nypix = -9999;
all_edgex = -9999;
all_edgey = -9999;
all_bigx = -9999;
all_bigy = -9999;
all_alpha = -9999;
all_beta = -9999;
all_simphi = -9999;
all_simtheta = -9999;
all_simhitx = -9999;
all_simhity = -9999;
all_nsimhit = -9999;
all_pidhit = -9999;
all_simproc = -9999;
all_vtxr = -9999;
all_vtxz = -9999;
all_simpx = -9999;
all_simpy = -9999;
all_simpz = -9999;
all_eloss = -9999;
all_trkid = -9999;
all_x1 = -9999;
all_x2 = -9999;
all_y1 = -9999;
all_y2 = -9999;
all_z1 = -9999;
all_z2 = -9999;
all_row1 = -9999;
all_row2 = -9999;
all_col1 = -9999;
all_col2 = -9999;
all_gx1 = -9999;
all_gx2 = -9999;
all_gy1 = -9999;
all_gy2 = -9999;
all_gz1 = -9999;
all_gz2 = -9999;
all_simtrketa = -9999;
all_simtrkphi = -9999;
all_clust_row = -9999;
all_clust_col = -9999;
all_clust_x = -9999;
all_clust_y = -9999;
all_clust_q = -9999;
all_clust_maxpixcol = -9999;
all_clust_maxpixrow = -9999;
all_clust_minpixcol = -9999;
all_clust_minpixrow = -9999;
all_clust_geoid = -9999;
all_clust_alpha = -9999;
all_clust_beta = -9999;
/*
for (int i=0; i<all_npix; ++i)
{
all_pixrow[i] = -9999;
all_pixcol[i] = -9999;
all_pixadc[i] = -9999;
all_pixx[i] = -9999;
all_pixy[i] = -9999;
all_pixgx[i] = -9999;
all_pixgy[i] = -9999;
all_pixgz[i] = -9999;
}
*/
all_hit_probx = -9999;
all_hit_proby = -9999;
all_hit_cprob0 = -9999;
all_hit_cprob1 = -9999;
all_hit_cprob2 = -9999;
all_nsimhit = (int)matched.size();
all_subdetid = (int)detId.subdetId();
// only consider rechits in pixel barrel and pixel forward
if ( !(all_subdetid==1 || all_subdetid==2) )
{
cout << "SiPixelErrorEstimation::analyze: Not in a pixel detector !!!!!" << endl;
continue;
}
const PixelGeomDetUnit* theGeomDet
= dynamic_cast<const PixelGeomDetUnit*> ( tracker->idToDet(detId) );
const PixelTopology* topol = &(theGeomDet->specificTopology());
const int maxPixelCol = pixeliter->cluster()->maxPixelCol();
const int maxPixelRow = pixeliter->cluster()->maxPixelRow();
const int minPixelCol = pixeliter->cluster()->minPixelCol();
const int minPixelRow = pixeliter->cluster()->minPixelRow();
//all_hit_probx = (float)pixeliter->probabilityX();
//all_hit_proby = (float)pixeliter->probabilityY();
all_hit_cprob0 = (float)pixeliter->clusterProbability(0);
all_hit_cprob1 = (float)pixeliter->clusterProbability(1);
all_hit_cprob2 = (float)pixeliter->clusterProbability(2);
// check whether the cluster is at the module edge
if ( topol->isItEdgePixelInX( minPixelRow ) ||
topol->isItEdgePixelInX( maxPixelRow ) )
all_edgex = 1;
else
all_edgex = 0;
if ( topol->isItEdgePixelInY( minPixelCol ) ||
topol->isItEdgePixelInY( maxPixelCol ) )
all_edgey = 1;
else
all_edgey = 0;
// check whether this rechit contains big pixels
if ( topol->containsBigPixelInX(minPixelRow, maxPixelRow) )
all_bigx = 1;
else
all_bigx = 0;
if ( topol->containsBigPixelInY(minPixelCol, maxPixelCol) )
all_bigy = 1;
else
all_bigy = 0;
if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel )
{
PXBDetId bdetid(detId);
all_layer = bdetid.layer();
all_ladder = bdetid.ladder();
all_mod = bdetid.module();
int tmp_nrows = theGeomDet->specificTopology().nrows();
if ( tmp_nrows == 80 )
all_half = 1;
else if ( tmp_nrows == 160 )
all_half = 0;
else
cout << "-------------------------------------------------- Wrong module size !!!" << endl;
float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
if ( tmp2<tmp1 )
all_flipped = 1;
else
all_flipped = 0;
}
else if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelEndcap )
{
PXFDetId fdetid(detId);
all_side = fdetid.side();
all_disk = fdetid.disk();
all_blade = fdetid.blade();
all_panel = fdetid.panel();
all_plaq = fdetid.module(); // also known as plaquette
} // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
else std::cout << "We are not in the pixel detector" << (int)detId.subdetId() << endl;
all_cols = theGeomDet->specificTopology().ncolumns();
all_rows = theGeomDet->specificTopology().nrows();
LocalPoint lp = pixeliter->localPosition();
// gavril: change this name
all_rechitx = lp.x();
all_rechity = lp.y();
all_rechitz = lp.z();
LocalError le = pixeliter->localPositionError();
all_rechiterrx = sqrt( le.xx() );
all_rechiterry = sqrt( le.yy() );
bool found_hit_from_generated_particle = false;
//---Loop over sim hits, fill closest
float closest_dist = 99999.9;
std::vector<PSimHit>::const_iterator closest_simhit = matched.begin();
for (std::vector<PSimHit>::const_iterator m = matched.begin(); m < matched.end(); m++)
{
if ( checkType_ )
{
int pid = (*m).particleType();
if ( abs(pid) != genType_ )
continue;
}
float simhitx = 0.5 * ( (*m).entryPoint().x() + (*m).exitPoint().x() );
float simhity = 0.5 * ( (*m).entryPoint().y() + (*m).exitPoint().y() );
float x_res = simhitx - rechitx;
float y_res = simhity - rechity;
float dist = sqrt( x_res*x_res + y_res*y_res );
if ( dist < closest_dist )
{
closest_dist = dist;
closest_simhit = m;
found_hit_from_generated_particle = true;
}
} // end sim hit loop
// If this recHit does not have any simHit with the same particleType as the particles generated
// ignore it as most probably comes from delta rays.
if ( checkType_ && !found_hit_from_generated_particle )
continue;
all_x1 = (*closest_simhit).entryPoint().x(); // width (row index, in col direction)
all_y1 = (*closest_simhit).entryPoint().y(); // length (col index, in row direction)
all_z1 = (*closest_simhit).entryPoint().z();
all_x2 = (*closest_simhit).exitPoint().x();
all_y2 = (*closest_simhit).exitPoint().y();
all_z2 = (*closest_simhit).exitPoint().z();
GlobalPoint GP1 =
theGeomDet->surface().toGlobal( Local3DPoint( (*closest_simhit).entryPoint().x(),
(*closest_simhit).entryPoint().y(),
(*closest_simhit).entryPoint().z() ) );
GlobalPoint GP2 =
theGeomDet->surface().toGlobal (Local3DPoint( (*closest_simhit).exitPoint().x(),
(*closest_simhit).exitPoint().y(),
(*closest_simhit).exitPoint().z() ) );
all_gx1 = GP1.x();
all_gx2 = GP2.x();
all_gy1 = GP1.y();
all_gy2 = GP2.y();
all_gz1 = GP1.z();
all_gz2 = GP2.z();
MeasurementPoint mp1 =
topol->measurementPosition( LocalPoint( (*closest_simhit).entryPoint().x(),
(*closest_simhit).entryPoint().y(),
(*closest_simhit).entryPoint().z() ) );
MeasurementPoint mp2 =
topol->measurementPosition( LocalPoint( (*closest_simhit).exitPoint().x(),
(*closest_simhit).exitPoint().y(),
(*closest_simhit).exitPoint().z() ) );
all_row1 = mp1.x();
all_col1 = mp1.y();
all_row2 = mp2.x();
all_col2 = mp2.y();
all_simhitx = 0.5*(all_x1+all_x2);
all_simhity = 0.5*(all_y1+all_y2);
all_rechitresx = all_rechitx - all_simhitx;
all_rechitresy = all_rechity - all_simhity;
all_rechitpullx = all_rechitresx / all_rechiterrx;
all_rechitpully = all_rechitresy / all_rechiterry;
SiPixelRecHit::ClusterRef const& clust = pixeliter->cluster();
all_npix = clust->size();
all_nxpix = clust->sizeX();
all_nypix = clust->sizeY();
all_clust_row = clust->x();
all_clust_col = clust->y();
LocalPoint lp2 = topol->localPosition( MeasurementPoint( all_clust_row, all_clust_col ) );
all_clust_x = lp2.x();
all_clust_y = lp2.y();
all_clust_q = clust->charge();
all_clust_maxpixcol = clust->maxPixelCol();
all_clust_maxpixrow = clust->maxPixelRow();
all_clust_minpixcol = clust->minPixelCol();
all_clust_minpixrow = clust->minPixelRow();
all_clust_geoid = clust->geographicalId();
all_simpx = (*closest_simhit).momentumAtEntry().x();
all_simpy = (*closest_simhit).momentumAtEntry().y();
all_simpz = (*closest_simhit).momentumAtEntry().z();
all_eloss = (*closest_simhit).energyLoss();
all_simphi = (*closest_simhit).phiAtEntry();
all_simtheta = (*closest_simhit).thetaAtEntry();
all_pidhit = (*closest_simhit).particleType();
all_trkid = (*closest_simhit).trackId();
//--- Fill alpha and beta -- more useful for exploring the residuals...
all_beta = atan2(all_simpz, all_simpy);
all_alpha = atan2(all_simpz, all_simpx);
all_simproc = (int)closest_simhit->processType();
const edm::SimTrackContainer& trks = *(simtracks.product());
SimTrackContainer::const_iterator trksiter;
for (trksiter = trks.begin(); trksiter != trks.end(); trksiter++)
if ( (int)trksiter->trackId() == all_trkid )
{
all_simtrketa = trksiter->momentum().eta();
all_simtrkphi = trksiter->momentum().phi();
}
all_vtxz = theGeomDet->surface().position().z();
all_vtxr = theGeomDet->surface().position().perp();
//computeAnglesFromDetPosition(clust,
// theGeomDet,
// all_clust_alpha, all_clust_beta )
const std::vector<SiPixelCluster::Pixel>& pixvector = clust->pixels();
for ( int i=0; i<(int)pixvector.size(); ++i)
{
SiPixelCluster::Pixel holdpix = pixvector[i];
all_pixrow[i] = holdpix.x;
all_pixcol[i] = holdpix.y;
all_pixadc[i] = holdpix.adc;
LocalPoint lp = topol->localPosition(MeasurementPoint(holdpix.x, holdpix.y));
all_pixx[i] = lp.x();
all_pixy[i]= lp.y();
GlobalPoint GP = theGeomDet->surface().toGlobal(Local3DPoint(lp.x(),lp.y(),lp.z()));
all_pixgx[i] = GP.x();
all_pixgy[i]= GP.y();
all_pixgz[i]= GP.z();
}
ttree_all_hits_->Fill();
} // for ( ; pixeliter != pixelrechitRangeIteratorEnd; ++pixeliter)
} // for (TrackerGeometry::DetContainer::const_iterator it = pDD->dets().begin(); it != pDD->dets().end(); it++)
// ------------------------------------------------ all hits ---------------------------------------------------------------
// ------------------------------------------------ track hits only --------------------------------------------------------
if ( include_trk_hits_ )
{
// Get tracks
edm::Handle<reco::TrackCollection> trackCollection;
e.getByLabel(src_, trackCollection);
const reco::TrackCollection *tracks = trackCollection.product();
reco::TrackCollection::const_iterator tciter;
if ( tracks->size() > 0 )
{
// Loop on tracks
for ( tciter=tracks->begin(); tciter!=tracks->end(); ++tciter)
{
// First loop on hits: find matched hits
for ( trackingRecHit_iterator it = tciter->recHitsBegin(); it != tciter->recHitsEnd(); ++it)
{
const TrackingRecHit &thit = **it;
// Is it a matched hit?
const SiPixelRecHit* matchedhit = dynamic_cast<const SiPixelRecHit*>(&thit);
if ( matchedhit )
{
rechitx = -9999.9;
rechity = -9999.9;
rechitz = -9999.9;
rechiterrx = -9999.9;
rechiterry = -9999.9;
rechitresx = -9999.9;
rechitresy = -9999.9;
rechitpullx = -9999.9;
rechitpully = -9999.9;
npix = -9999;
nxpix = -9999;
nypix = -9999;
charge = -9999.9;
edgex = -9999;
edgey = -9999;
bigx = -9999;
bigy = -9999;
alpha = -9999.9;
beta = -9999.9;
phi = -9999.9;
eta = -9999.9;
subdetId = -9999;
layer = -9999;
ladder = -9999;
mod = -9999;
side = -9999;
disk = -9999;
blade = -9999;
panel = -9999;
plaq = -9999;
half = -9999;
flipped = -9999;
nsimhit = -9999;
pidhit = -9999;
simproc = -9999;
simhitx = -9999.9;
simhity = -9999.9;
hit_probx = -9999.9;
hit_proby = -9999.9;
hit_cprob0 = -9999.9;
hit_cprob1 = -9999.9;
hit_cprob2 = -9999.9;
position = (*it)->localPosition();
error = (*it)->localPositionError();
rechitx = position.x();
rechity = position.y();
rechitz = position.z();
rechiterrx = sqrt(error.xx());
rechiterry = sqrt(error.yy());
npix = matchedhit->cluster()->size();
nxpix = matchedhit->cluster()->sizeX();
nypix = matchedhit->cluster()->sizeY();
charge = matchedhit->cluster()->charge();
//hit_probx = (float)matchedhit->probabilityX();
//hit_proby = (float)matchedhit->probabilityY();
hit_cprob0 = (float)matchedhit->clusterProbability(0);
hit_cprob1 = (float)matchedhit->clusterProbability(1);
hit_cprob2 = (float)matchedhit->clusterProbability(2);
//Association of the rechit to the simhit
matched.clear();
matched = associate.associateHit(*matchedhit);
nsimhit = (int)matched.size();
if ( !matched.empty() )
{
mindist = 999999.9;
float distx, disty, dist;
bool found_hit_from_generated_particle = false;
int n_assoc_muon = 0;
vector<PSimHit>::const_iterator closestit = matched.begin();
for (vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); ++m)
{
if ( checkType_ )
{ // only consider associated simhits with the generated pid (muons)
int pid = (*m).particleType();
if ( abs(pid) != genType_ )
continue;
}
float simhitx = 0.5 * ( (*m).entryPoint().x() + (*m).exitPoint().x() );
float simhity = 0.5 * ( (*m).entryPoint().y() + (*m).exitPoint().y() );
distx = fabs(rechitx - simhitx);
disty = fabs(rechity - simhity);
dist = sqrt( distx*distx + disty*disty );
if ( dist < mindist )
{
n_assoc_muon++;
mindist = dist;
closestit = m;
found_hit_from_generated_particle = true;
}
} // for (vector<PSimHit>::const_iterator m=matched.begin(); m<matched.end(); m++)
// This recHit does not have any simHit with the same particleType as the particles generated
// Ignore it as most probably come from delta rays.
if ( checkType_ && !found_hit_from_generated_particle )
continue;
//if ( n_assoc_muon > 1 )
//{
// // cout << " ----- This is not good: n_assoc_muon = " << n_assoc_muon << endl;
// // cout << "evt = " << evt << endl;
//}
detId = (*it)->geographicalId();
const PixelGeomDetUnit* theGeomDet =
dynamic_cast<const PixelGeomDetUnit*> ((*tracker).idToDet(detId) );
const PixelTopology* theTopol = &(theGeomDet->specificTopology());
pidhit = (*closestit).particleType();
simproc = (int)(*closestit).processType();
simhitx = 0.5*( (*closestit).entryPoint().x() + (*closestit).exitPoint().x() );
simhity = 0.5*( (*closestit).entryPoint().y() + (*closestit).exitPoint().y() );
rechitresx = rechitx - simhitx;
rechitresy = rechity - simhity;
rechitpullx = ( rechitx - simhitx ) / sqrt(error.xx());
rechitpully = ( rechity - simhity ) / sqrt(error.yy());
float simhitpx = (*closestit).momentumAtEntry().x();
float simhitpy = (*closestit).momentumAtEntry().y();
float simhitpz = (*closestit).momentumAtEntry().z();
beta = atan2(simhitpz, simhitpy) * radtodeg;
alpha = atan2(simhitpz, simhitpx) * radtodeg;
//beta = fabs(atan2(simhitpz, simhitpy)) * radtodeg;
//alpha = fabs(atan2(simhitpz, simhitpx)) * radtodeg;
// calculate alpha and beta exactly as in PixelCPEBase.cc
float locx = simhitpx;
float locy = simhitpy;
float locz = simhitpz;
bool isFlipped = false;
float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
if ( tmp2<tmp1 )
isFlipped = true;
else
isFlipped = false;
trk_alpha = acos(locx/sqrt(locx*locx+locz*locz)) * radtodeg;
if ( isFlipped ) // &&& check for FPIX !!!
trk_alpha = 180.0 - trk_alpha ;
trk_beta = acos(locy/sqrt(locy*locy+locz*locz)) * radtodeg;
phi = tciter->momentum().phi() / math_pi*180.0;
eta = tciter->momentum().eta();
const int maxPixelCol = (*matchedhit).cluster()->maxPixelCol();
const int maxPixelRow = (*matchedhit).cluster()->maxPixelRow();
const int minPixelCol = (*matchedhit).cluster()->minPixelCol();
const int minPixelRow = (*matchedhit).cluster()->minPixelRow();
// check whether the cluster is at the module edge
if ( theTopol->isItEdgePixelInX( minPixelRow ) ||
theTopol->isItEdgePixelInX( maxPixelRow ) )
edgex = 1;
else
edgex = 0;
if ( theTopol->isItEdgePixelInY( minPixelCol ) ||
theTopol->isItEdgePixelInY( maxPixelCol ) )
edgey = 1;
else
edgey = 0;
// check whether this rechit contains big pixels
if ( theTopol->containsBigPixelInX(minPixelRow, maxPixelRow) )
bigx = 1;
else
bigx = 0;
if ( theTopol->containsBigPixelInY(minPixelCol, maxPixelCol) )
bigy = 1;
else
bigy = 0;
subdetId = (int)detId.subdetId();
if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelBarrel )
{
int tmp_nrows = theGeomDet->specificTopology().nrows();
if ( tmp_nrows == 80 )
half = 1;
else if ( tmp_nrows == 160 )
half = 0;
else
cout << "-------------------------------------------------- Wrong module size !!!" << endl;
float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
if ( tmp2<tmp1 )
flipped = 1;
else
flipped = 0;
PXBDetId bdetid(detId);
layer = bdetid.layer(); // Layer: 1,2,3.
ladder = bdetid.ladder(); // Ladder: 1-20, 32, 44.
mod = bdetid.module(); // Mod: 1-8.
}
else if ( (int)detId.subdetId() == (int)PixelSubdetector::PixelEndcap )
{
PXFDetId fdetid(detId);
side = fdetid.side();
disk = fdetid.disk();
blade = fdetid.blade();
panel = fdetid.panel();
plaq = fdetid.module(); // also known as plaquette
float tmp1 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,0.)).perp();
float tmp2 = theGeomDet->surface().toGlobal(Local3DPoint(0.,0.,1.)).perp();
if ( tmp2<tmp1 )
flipped = 1;
else
flipped = 0;
} // else if ( detId.subdetId()==PixelSubdetector::PixelEndcap )
//else std::// cout << "We are not in the pixel detector. detId.subdetId() = " << (int)detId.subdetId() << endl;
ttree_track_hits_->Fill();
} // if ( !matched.empty() )
else
cout << "---------------- RecHit with no associated SimHit !!! -------------------------- " << endl;
} // if ( matchedhit )
} // end of loop on hits
} //end of loop on track
} // tracks > 0.
} // if ( include_trk_hits_ )
// ----------------------------------------------- track hits only -----------------------------------------------------------
}
| void SiPixelErrorEstimation::beginJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 47 of file SiPixelErrorEstimation.cc.
References all_alpha, all_beta, all_bigx, all_bigy, all_blade, all_clust_alpha, all_clust_beta, all_clust_col, all_clust_geoid, all_clust_maxpixcol, all_clust_maxpixrow, all_clust_minpixcol, all_clust_minpixrow, all_clust_q, all_clust_row, all_clust_x, all_clust_y, all_col1, all_col2, all_cols, all_disk, all_edgex, all_edgey, all_eloss, all_flipped, all_gx1, all_gx2, all_half, all_hit_cprob0, all_hit_cprob1, all_hit_cprob2, all_hit_probx, all_hit_proby, all_ladder, all_layer, all_mod, all_npix, all_nsimhit, all_nxpix, all_nypix, all_panel, all_pidhit, all_pixadc, all_pixcol, all_pixgx, all_pixgy, all_pixgz, all_pixrow, all_pixx, all_pixy, all_plaq, all_rechiterrx, all_rechiterry, all_rechitpullx, all_rechitpully, all_rechitresx, all_rechitresy, all_rechitx, all_rechity, all_rechitz, all_row1, all_row2, all_rows, all_side, all_simhitx, all_simhity, all_simphi, all_simproc, all_simpx, all_simpy, all_simpz, all_simtheta, all_simtrketa, all_simtrkphi, all_subdetid, all_trkid, all_vtxr, all_vtxz, all_x1, all_x2, alpha, beta, bigx, bigy, blade, charge, disk, edgex, edgey, eta, evt, flipped, half, hit_cprob0, hit_cprob1, hit_cprob2, hit_probx, hit_proby, include_trk_hits_, ladder, layer, mod, npix, nsimhit, nxpix, nypix, outputFile_, panel, phi, pidhit, plaq, rechiterrx, rechiterry, rechitpullx, rechitpully, rechitresx, rechitresy, rechitx, rechity, rechitz, run, side, simhitx, simhity, simproc, subdetId, tfile_, trk_alpha, trk_beta, ttree_all_hits_, and ttree_track_hits_.
{
int bufsize = 64000;
if ( include_trk_hits_ )
{
//tfile_ = new TFile ("SiPixelErrorEstimation_Ntuple.root" , "RECREATE");
//const char* tmp_name = outputFile_.c_str();
tfile_ = new TFile ( outputFile_.c_str() , "RECREATE");
ttree_track_hits_ = new TTree("TrackHitNtuple", "TrackHitNtuple");
ttree_track_hits_->Branch("evt", &evt, "evt/I", bufsize);
ttree_track_hits_->Branch("run", &run, "run/I", bufsize);
ttree_track_hits_->Branch("subdetId", &subdetId, "subdetId/I", bufsize);
ttree_track_hits_->Branch("layer" , &layer , "layer/I" , bufsize);
ttree_track_hits_->Branch("ladder", &ladder, "ladder/I", bufsize);
ttree_track_hits_->Branch("mod" , &mod , "mod/I" , bufsize);
ttree_track_hits_->Branch("side" , &side , "side/I" , bufsize);
ttree_track_hits_->Branch("disk" , &disk , "disk/I" , bufsize);
ttree_track_hits_->Branch("blade" , &blade , "blade/I" , bufsize);
ttree_track_hits_->Branch("panel" , &panel , "panel/I" , bufsize);
ttree_track_hits_->Branch("plaq" , &plaq , "plaq/I" , bufsize);
ttree_track_hits_->Branch("half" , &half , "half/I" , bufsize);
ttree_track_hits_->Branch("flipped", &flipped, "flipped/I", bufsize);
ttree_track_hits_->Branch("rechitx", &rechitx, "rechitx/F" , bufsize);
ttree_track_hits_->Branch("rechity", &rechity, "rechity/F" , bufsize);
ttree_track_hits_->Branch("rechitz", &rechitz, "rechitz/F" , bufsize);
ttree_track_hits_->Branch("rechiterrx", &rechiterrx, "rechiterrx/F" , bufsize);
ttree_track_hits_->Branch("rechiterry", &rechiterry, "rechiterry/F" , bufsize);
ttree_track_hits_->Branch("rechitresx", &rechitresx, "rechitresx/F" , bufsize);
ttree_track_hits_->Branch("rechitresy", &rechitresy, "rechitresy/F" , bufsize);
ttree_track_hits_->Branch("rechitpullx", &rechitpullx, "rechitpullx/F", bufsize);
ttree_track_hits_->Branch("rechitpully", &rechitpully, "rechitpully/F", bufsize);
ttree_track_hits_->Branch("npix" , &npix , "npix/I" , bufsize);
ttree_track_hits_->Branch("nxpix" , &nxpix , "nxpix/I" , bufsize);
ttree_track_hits_->Branch("nypix" , &nypix , "nypix/I" , bufsize);
ttree_track_hits_->Branch("charge", &charge, "charge/F", bufsize);
ttree_track_hits_->Branch("edgex", &edgex, "edgex/I", bufsize);
ttree_track_hits_->Branch("edgey", &edgey, "edgey/I", bufsize);
ttree_track_hits_->Branch("bigx", &bigx, "bigx/I", bufsize);
ttree_track_hits_->Branch("bigy", &bigy, "bigy/I", bufsize);
ttree_track_hits_->Branch("alpha", &alpha, "alpha/F", bufsize);
ttree_track_hits_->Branch("beta" , &beta , "beta/F" , bufsize);
ttree_track_hits_->Branch("trk_alpha", &trk_alpha, "trk_alpha/F", bufsize);
ttree_track_hits_->Branch("trk_beta" , &trk_beta , "trk_beta/F" , bufsize);
ttree_track_hits_->Branch("phi", &phi, "phi/F", bufsize);
ttree_track_hits_->Branch("eta", &eta, "eta/F", bufsize);
ttree_track_hits_->Branch("simhitx", &simhitx, "simhitx/F", bufsize);
ttree_track_hits_->Branch("simhity", &simhity, "simhity/F", bufsize);
ttree_track_hits_->Branch("nsimhit", &nsimhit, "nsimhit/I", bufsize);
ttree_track_hits_->Branch("pidhit" , &pidhit , "pidhit/I" , bufsize);
ttree_track_hits_->Branch("simproc", &simproc, "simproc/I", bufsize);
ttree_track_hits_->Branch("hit_probx" , &hit_probx , "hit_probx/F" , bufsize);
ttree_track_hits_->Branch("hit_proby" , &hit_proby , "hit_proby/F" , bufsize);
ttree_track_hits_->Branch("hit_cprob0", &hit_cprob0, "hit_cprob0/F", bufsize);
ttree_track_hits_->Branch("hit_cprob1", &hit_cprob1, "hit_cprob1/F", bufsize);
ttree_track_hits_->Branch("hit_cprob2", &hit_cprob2, "hit_cprob2/F", bufsize);
} // if ( include_trk_hits_ )
// ----------------------------------------------------------------------
ttree_all_hits_ = new TTree("AllHitNtuple", "AllHitNtuple");
ttree_all_hits_->Branch("evt", &evt, "evt/I", bufsize);
ttree_all_hits_->Branch("run", &run, "run/I", bufsize);
ttree_all_hits_->Branch("subdetid", &all_subdetid, "subdetid/I", bufsize);
ttree_all_hits_->Branch("layer" , &all_layer , "layer/I" , bufsize);
ttree_all_hits_->Branch("ladder", &all_ladder, "ladder/I", bufsize);
ttree_all_hits_->Branch("mod" , &all_mod , "mod/I" , bufsize);
ttree_all_hits_->Branch("side" , &all_side , "side/I" , bufsize);
ttree_all_hits_->Branch("disk" , &all_disk , "disk/I" , bufsize);
ttree_all_hits_->Branch("blade" , &all_blade , "blade/I" , bufsize);
ttree_all_hits_->Branch("panel" , &all_panel , "panel/I" , bufsize);
ttree_all_hits_->Branch("plaq" , &all_plaq , "plaq/I" , bufsize);
ttree_all_hits_->Branch("half" , &all_half , "half/I" , bufsize);
ttree_all_hits_->Branch("flipped", &all_flipped, "flipped/I", bufsize);
ttree_all_hits_->Branch("cols", &all_cols, "cols/I", bufsize);
ttree_all_hits_->Branch("rows", &all_rows, "rows/I", bufsize);
ttree_all_hits_->Branch("rechitx" , &all_rechitx , "rechitx/F" , bufsize);
ttree_all_hits_->Branch("rechity" , &all_rechity , "rechity/F" , bufsize);
ttree_all_hits_->Branch("rechitz" , &all_rechitz , "rechitz/F" , bufsize);
ttree_all_hits_->Branch("rechiterrx" , &all_rechiterrx , "rechiterrx/F" , bufsize);
ttree_all_hits_->Branch("rechiterry" , &all_rechiterry , "rechiterry/F" , bufsize);
ttree_all_hits_->Branch("rechitresx" , &all_rechitresx , "rechitresx/F" , bufsize);
ttree_all_hits_->Branch("rechitresy" , &all_rechitresy , "rechitresy/F" , bufsize);
ttree_all_hits_->Branch("rechitpullx", &all_rechitpullx, "rechitpullx/F", bufsize);
ttree_all_hits_->Branch("rechitpully", &all_rechitpully, "rechitpully/F", bufsize);
ttree_all_hits_->Branch("npix" , &all_npix , "npix/I" , bufsize);
ttree_all_hits_->Branch("nxpix" , &all_nxpix , "nxpix/I" , bufsize);
ttree_all_hits_->Branch("nypix" , &all_nypix , "nypix/I" , bufsize);
ttree_all_hits_->Branch("edgex", &all_edgex, "edgex/I", bufsize);
ttree_all_hits_->Branch("edgey", &all_edgey, "edgey/I", bufsize);
ttree_all_hits_->Branch("bigx", &all_bigx, "bigx/I", bufsize);
ttree_all_hits_->Branch("bigy", &all_bigy, "bigy/I", bufsize);
ttree_all_hits_->Branch("alpha", &all_alpha, "alpha/F", bufsize);
ttree_all_hits_->Branch("beta" , &all_beta , "beta/F" , bufsize);
ttree_all_hits_->Branch("simhitx", &all_simhitx, "simhitx/F", bufsize);
ttree_all_hits_->Branch("simhity", &all_simhity, "simhity/F", bufsize);
ttree_all_hits_->Branch("nsimhit", &all_nsimhit, "nsimhit/I", bufsize);
ttree_all_hits_->Branch("pidhit" , &all_pidhit , "pidhit/I" , bufsize);
ttree_all_hits_->Branch("simproc", &all_simproc, "simproc/I", bufsize);
ttree_all_hits_->Branch("vtxr", &all_vtxr, "vtxr/F", bufsize);
ttree_all_hits_->Branch("vtxz", &all_vtxz, "vtxz/F", bufsize);
ttree_all_hits_->Branch("simpx", &all_simpx, "simpx/F", bufsize);
ttree_all_hits_->Branch("simpy", &all_simpy, "simpy/F", bufsize);
ttree_all_hits_->Branch("simpz", &all_simpz, "simpz/F", bufsize);
ttree_all_hits_->Branch("eloss", &all_eloss, "eloss/F", bufsize);
ttree_all_hits_->Branch("simphi", &all_simphi, "simphi/F", bufsize);
ttree_all_hits_->Branch("simtheta", &all_simtheta, "simtheta/F", bufsize);
ttree_all_hits_->Branch("trkid", &all_trkid, "trkid/I", bufsize);
ttree_all_hits_->Branch("x1", &all_x1, "x1/F", bufsize);
ttree_all_hits_->Branch("x2", &all_x2, "x2/F", bufsize);
ttree_all_hits_->Branch("y1", &all_x1, "y1/F", bufsize);
ttree_all_hits_->Branch("y2", &all_x2, "y2/F", bufsize);
ttree_all_hits_->Branch("z1", &all_x1, "z1/F", bufsize);
ttree_all_hits_->Branch("z2", &all_x2, "z2/F", bufsize);
ttree_all_hits_->Branch("row1", &all_row1, "row1/F", bufsize);
ttree_all_hits_->Branch("row2", &all_row2, "row2/F", bufsize);
ttree_all_hits_->Branch("col1", &all_col1, "col1/F", bufsize);
ttree_all_hits_->Branch("col2", &all_col2, "col2/F", bufsize);
ttree_all_hits_->Branch("gx1", &all_gx1, "gx1/F", bufsize);
ttree_all_hits_->Branch("gx2", &all_gx2, "gx2/F", bufsize);
ttree_all_hits_->Branch("gy1", &all_gx1, "gy1/F", bufsize);
ttree_all_hits_->Branch("gy2", &all_gx2, "gy2/F", bufsize);
ttree_all_hits_->Branch("gz1", &all_gx1, "gz1/F", bufsize);
ttree_all_hits_->Branch("gz2", &all_gx2, "gz2/F", bufsize);
ttree_all_hits_->Branch("simtrketa", &all_simtrketa, "simtrketa/F", bufsize);
ttree_all_hits_->Branch("simtrkphi", &all_simtrkphi, "simtrkphi/F", bufsize);
ttree_all_hits_->Branch("clust_row", &all_clust_row, "clust_row/F", bufsize);
ttree_all_hits_->Branch("clust_col", &all_clust_col, "clust_col/F", bufsize);
ttree_all_hits_->Branch("clust_x", &all_clust_x, "clust_x/F", bufsize);
ttree_all_hits_->Branch("clust_y", &all_clust_y, "clust_y/F", bufsize);
ttree_all_hits_->Branch("clust_q", &all_clust_q, "clust_q/F", bufsize);
ttree_all_hits_->Branch("clust_maxpixcol", &all_clust_maxpixcol, "clust_maxpixcol/I", bufsize);
ttree_all_hits_->Branch("clust_maxpixrow", &all_clust_maxpixrow, "clust_maxpixrow/I", bufsize);
ttree_all_hits_->Branch("clust_minpixcol", &all_clust_minpixcol, "clust_minpixcol/I", bufsize);
ttree_all_hits_->Branch("clust_minpixrow", &all_clust_minpixrow, "clust_minpixrow/I", bufsize);
ttree_all_hits_->Branch("clust_geoid", &all_clust_geoid, "clust_geoid/I", bufsize);
ttree_all_hits_->Branch("clust_alpha", &all_clust_alpha, "clust_alpha/F", bufsize);
ttree_all_hits_->Branch("clust_beta" , &all_clust_beta , "clust_beta/F" , bufsize);
ttree_all_hits_->Branch("rowpix", all_pixrow, "row[npix]/F", bufsize);
ttree_all_hits_->Branch("colpix", all_pixcol, "col[npix]/F", bufsize);
ttree_all_hits_->Branch("adc", all_pixadc, "adc[npix]/F", bufsize);
ttree_all_hits_->Branch("xpix", all_pixx, "x[npix]/F", bufsize);
ttree_all_hits_->Branch("ypix", all_pixy, "y[npix]/F", bufsize);
ttree_all_hits_->Branch("gxpix", all_pixgx, "gx[npix]/F", bufsize);
ttree_all_hits_->Branch("gypix", all_pixgy, "gy[npix]/F", bufsize);
ttree_all_hits_->Branch("gzpix", all_pixgz, "gz[npix]/F", bufsize);
ttree_all_hits_->Branch("all_hit_probx", &all_hit_probx, "all_hit_probx/F" , bufsize);
ttree_all_hits_->Branch("all_hit_proby", &all_hit_proby, "all_hit_proby/F" , bufsize);
ttree_all_hits_->Branch("all_hit_cprob0", &all_hit_cprob0, "all_hit_cprob0/F", bufsize);
ttree_all_hits_->Branch("all_hit_cprob1", &all_hit_cprob1, "all_hit_cprob1/F", bufsize);
ttree_all_hits_->Branch("all_hit_cprob2", &all_hit_cprob2, "all_hit_cprob2/F", bufsize);
}
| void SiPixelErrorEstimation::computeAnglesFromDetPosition | ( | const SiPixelCluster & | cl, |
| const GeomDetUnit & | det, | ||
| float & | alpha, | ||
| float & | beta | ||
| ) |
Definition at line 1004 of file SiPixelErrorEstimation.cc.
References gather_cfg::cout, Topology::localPosition(), PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), GeomDet::surface(), Surface::toGlobal(), SiPixelCluster::x(), PV3DBase< T, PVType, FrameType >::x(), SiPixelCluster::y(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().
{
//--- This is a new det unit, so cache it
const PixelGeomDetUnit* theDet = dynamic_cast<const PixelGeomDetUnit*>( &det );
if (! theDet)
{
cout << "---------------------------------------------- Not a pixel detector !!!!!!!!!!!!!!" << endl;
assert(0);
}
const PixelTopology* theTopol = &(theDet->specificTopology());
// get cluster center of gravity (of charge)
float xcenter = cl.x();
float ycenter = cl.y();
// get the cluster position in local coordinates (cm)
LocalPoint lp = theTopol->localPosition( MeasurementPoint(xcenter, ycenter) );
//float lp_mod = sqrt( lp.x()*lp.x() + lp.y()*lp.y() + lp.z()*lp.z() );
// get the cluster position in global coordinates (cm)
GlobalPoint gp = theDet->surface().toGlobal( lp );
float gp_mod = sqrt( gp.x()*gp.x() + gp.y()*gp.y() + gp.z()*gp.z() );
// normalize
float gpx = gp.x()/gp_mod;
float gpy = gp.y()/gp_mod;
float gpz = gp.z()/gp_mod;
// make a global vector out of the global point; this vector will point from the
// origin of the detector to the cluster
GlobalVector gv(gpx, gpy, gpz);
// make local unit vector along local X axis
const Local3DVector lvx(1.0, 0.0, 0.0);
// get the unit X vector in global coordinates/
GlobalVector gvx = theDet->surface().toGlobal( lvx );
// make local unit vector along local Y axis
const Local3DVector lvy(0.0, 1.0, 0.0);
// get the unit Y vector in global coordinates
GlobalVector gvy = theDet->surface().toGlobal( lvy );
// make local unit vector along local Z axis
const Local3DVector lvz(0.0, 0.0, 1.0);
// get the unit Z vector in global coordinates
GlobalVector gvz = theDet->surface().toGlobal( lvz );
// calculate the components of gv (the unit vector pointing to the cluster)
// in the local coordinate system given by the basis {gvx, gvy, gvz}
// note that both gv and the basis {gvx, gvy, gvz} are given in global coordinates
float gv_dot_gvx = gv.x()*gvx.x() + gv.y()*gvx.y() + gv.z()*gvx.z();
float gv_dot_gvy = gv.x()*gvy.x() + gv.y()*gvy.y() + gv.z()*gvy.z();
float gv_dot_gvz = gv.x()*gvz.x() + gv.y()*gvz.y() + gv.z()*gvz.z();
// calculate angles
alpha = atan2( gv_dot_gvz, gv_dot_gvx );
beta = atan2( gv_dot_gvz, gv_dot_gvy );
// calculate cotalpha and cotbeta
// cotalpha_ = 1.0/tan(alpha_);
// cotbeta_ = 1.0/tan(beta_ );
// or like this
//cotalpha_ = gv_dot_gvx / gv_dot_gvz;
//cotbeta_ = gv_dot_gvy / gv_dot_gvz;
}
| void SiPixelErrorEstimation::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 254 of file SiPixelErrorEstimation.cc.
References tfile_.
float SiPixelErrorEstimation::all_alpha [private] |
Definition at line 190 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_beta [private] |
Definition at line 191 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_bigx [private] |
Definition at line 187 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_bigy [private] |
Definition at line 188 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_blade [private] |
Definition at line 154 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_clust_alpha [private] |
Definition at line 248 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_clust_beta [private] |
Definition at line 249 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_clust_col [private] |
Definition at line 234 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_clust_geoid [private] |
Definition at line 246 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_clust_maxpixcol [private] |
Definition at line 241 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_clust_maxpixrow [private] |
Definition at line 242 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_clust_minpixcol [private] |
Definition at line 243 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_clust_minpixrow [private] |
Definition at line 244 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_clust_q [private] |
Definition at line 239 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_clust_row [private] |
Definition at line 233 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_clust_x [private] |
Definition at line 236 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_clust_y [private] |
Definition at line 237 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_col1 [private] |
Definition at line 220 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_col2 [private] |
Definition at line 221 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_cols [private] |
Definition at line 161 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_disk [private] |
Definition at line 153 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_edgex [private] |
Definition at line 184 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_edgey [private] |
Definition at line 185 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_eloss [private] |
Definition at line 207 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_flipped [private] |
Definition at line 159 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_gx1 [private] |
Definition at line 223 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_gx2 [private] |
Definition at line 224 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_gy1 [private] |
Definition at line 225 of file SiPixelErrorEstimation.h.
Referenced by analyze().
float SiPixelErrorEstimation::all_gy2 [private] |
Definition at line 226 of file SiPixelErrorEstimation.h.
Referenced by analyze().
float SiPixelErrorEstimation::all_gz1 [private] |
Definition at line 227 of file SiPixelErrorEstimation.h.
Referenced by analyze().
float SiPixelErrorEstimation::all_gz2 [private] |
Definition at line 228 of file SiPixelErrorEstimation.h.
Referenced by analyze().
int SiPixelErrorEstimation::all_half [private] |
Definition at line 158 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_hit_cprob0 [private] |
Definition at line 264 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_hit_cprob1 [private] |
Definition at line 265 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_hit_cprob2 [private] |
Definition at line 266 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_hit_probx [private] |
Definition at line 262 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_hit_proby [private] |
Definition at line 263 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_ladder [private] |
Definition at line 149 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_layer [private] |
Definition at line 148 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_mod [private] |
Definition at line 150 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_npix [private] |
Definition at line 180 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_nsimhit [private] |
Definition at line 196 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_nxpix [private] |
Definition at line 181 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_nypix [private] |
Definition at line 182 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_panel [private] |
Definition at line 155 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_pidhit [private] |
Definition at line 197 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_pixadc[maxpix] [private] |
Definition at line 254 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_pixcol[maxpix] [private] |
Definition at line 253 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_pixgx[maxpix] [private] |
Definition at line 258 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_pixgy[maxpix] [private] |
Definition at line 259 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_pixgz[maxpix] [private] |
Definition at line 260 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_pixrow[maxpix] [private] |
Definition at line 252 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_pixx[maxpix] [private] |
Definition at line 256 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_pixy[maxpix] [private] |
Definition at line 257 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_plaq [private] |
Definition at line 156 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_rechiterrx [private] |
Definition at line 171 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_rechiterry [private] |
Definition at line 172 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_rechitpullx [private] |
Definition at line 177 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_rechitpully [private] |
Definition at line 178 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_rechitresx [private] |
Definition at line 174 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_rechitresy [private] |
Definition at line 175 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_rechitx [private] |
Definition at line 164 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_rechity [private] |
Definition at line 165 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_rechitz [private] |
Definition at line 166 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_row1 [private] |
Definition at line 218 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_row2 [private] |
Definition at line 219 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_rows [private] |
Definition at line 162 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_side [private] |
Definition at line 152 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_simhitx [private] |
Definition at line 168 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_simhity [private] |
Definition at line 169 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_simphi [private] |
Definition at line 193 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_simproc [private] |
Definition at line 198 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_simpx [private] |
Definition at line 203 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_simpy [private] |
Definition at line 204 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_simpz [private] |
Definition at line 205 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_simtheta [private] |
Definition at line 194 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_simtrketa [private] |
Definition at line 230 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_simtrkphi [private] |
Definition at line 231 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_subdetid [private] |
Definition at line 146 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::all_trkid [private] |
Definition at line 209 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_vtxr [private] |
Definition at line 200 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_vtxz [private] |
Definition at line 201 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_x1 [private] |
Definition at line 211 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_x2 [private] |
Definition at line 212 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::all_y1 [private] |
Definition at line 213 of file SiPixelErrorEstimation.h.
Referenced by analyze().
float SiPixelErrorEstimation::all_y2 [private] |
Definition at line 214 of file SiPixelErrorEstimation.h.
Referenced by analyze().
float SiPixelErrorEstimation::all_z1 [private] |
Definition at line 215 of file SiPixelErrorEstimation.h.
Referenced by analyze().
float SiPixelErrorEstimation::all_z2 [private] |
Definition at line 216 of file SiPixelErrorEstimation.h.
Referenced by analyze().
float SiPixelErrorEstimation::alpha [private] |
Definition at line 106 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::beta [private] |
Definition at line 107 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::bigx [private] |
Definition at line 103 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::bigy [private] |
Definition at line 104 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::blade [private] |
Definition at line 121 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::charge [private] |
Definition at line 98 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
bool SiPixelErrorEstimation::checkType_ [private] |
Definition at line 80 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and SiPixelErrorEstimation().
Definition at line 77 of file SiPixelErrorEstimation.h.
int SiPixelErrorEstimation::disk [private] |
Definition at line 120 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::edgex [private] |
Definition at line 100 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::edgey [private] |
Definition at line 101 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::eta [private] |
Definition at line 113 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::evt [private] |
Definition at line 135 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::flipped [private] |
Definition at line 126 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::genType_ [private] |
Definition at line 81 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and SiPixelErrorEstimation().
int SiPixelErrorEstimation::half [private] |
Definition at line 125 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::hit_cprob0 [private] |
Definition at line 140 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::hit_cprob1 [private] |
Definition at line 141 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::hit_cprob2 [private] |
Definition at line 142 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::hit_probx [private] |
Definition at line 138 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::hit_proby [private] |
Definition at line 139 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
bool SiPixelErrorEstimation::include_trk_hits_ [private] |
Definition at line 82 of file SiPixelErrorEstimation.h.
Referenced by analyze(), beginJob(), and SiPixelErrorEstimation().
int SiPixelErrorEstimation::ladder [private] |
Definition at line 117 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::layer [private] |
Definition at line 116 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
const int SiPixelErrorEstimation::maxpix = 10000 [static, private] |
Definition at line 251 of file SiPixelErrorEstimation.h.
int SiPixelErrorEstimation::mod [private] |
Definition at line 118 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::npix [private] |
Definition at line 95 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::nsimhit [private] |
Definition at line 128 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::nxpix [private] |
Definition at line 96 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::nypix [private] |
Definition at line 97 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
std::string SiPixelErrorEstimation::outputFile_ [private] |
Definition at line 78 of file SiPixelErrorEstimation.h.
Referenced by beginJob(), and SiPixelErrorEstimation().
int SiPixelErrorEstimation::panel [private] |
Definition at line 122 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::phi [private] |
Definition at line 112 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::pidhit [private] |
Definition at line 129 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::plaq [private] |
Definition at line 123 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::rechiterrx [private] |
Definition at line 88 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::rechiterry [private] |
Definition at line 89 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::rechitpullx [private] |
Definition at line 92 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::rechitpully [private] |
Definition at line 93 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::rechitresx [private] |
Definition at line 90 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::rechitresy [private] |
Definition at line 91 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::rechitx [private] |
Definition at line 85 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::rechity [private] |
Definition at line 86 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::rechitz [private] |
Definition at line 87 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::run [private] |
Definition at line 136 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::side [private] |
Definition at line 119 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::simhitx [private] |
Definition at line 132 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::simhity [private] |
Definition at line 133 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
int SiPixelErrorEstimation::simproc [private] |
Definition at line 130 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
std::string SiPixelErrorEstimation::src_ [private] |
Definition at line 79 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and SiPixelErrorEstimation().
int SiPixelErrorEstimation::subdetId [private] |
Definition at line 115 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
TFile* SiPixelErrorEstimation::tfile_ [private] |
Definition at line 270 of file SiPixelErrorEstimation.h.
Referenced by beginJob(), and endJob().
float SiPixelErrorEstimation::trk_alpha [private] |
Definition at line 109 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
float SiPixelErrorEstimation::trk_beta [private] |
Definition at line 110 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
TTree* SiPixelErrorEstimation::ttree_all_hits_ [private] |
Definition at line 271 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
TTree* SiPixelErrorEstimation::ttree_track_hits_ [private] |
Definition at line 272 of file SiPixelErrorEstimation.h.
Referenced by analyze(), and beginJob().
1.7.3