|
|
|
#include <SiPixelTrackResidualSource.h>
Definition at line 41 of file SiPixelTrackResidualSource.h.
| SiPixelTrackResidualSource::SiPixelTrackResidualSource | ( | const edm::ParameterSet & | pSet | ) | [explicit] |
Definition at line 64 of file SiPixelTrackResidualSource.cc.
References bladeOn, clustersrc_, dbe_, debug_, diskOn, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ladOn, layOn, modOn, cppFunctionSkipper::operator, phiOn, pSet_, ptminres_, ringOn, src_, AlCaHLTBitMon_QueryRunRegistry::string, tracksrc_, and ttrhbuilder_.
: pSet_(pSet), modOn( pSet.getUntrackedParameter<bool>("modOn",true) ), reducedSet( pSet.getUntrackedParameter<bool>("reducedSet",true) ), ladOn( pSet.getUntrackedParameter<bool>("ladOn",false) ), layOn( pSet.getUntrackedParameter<bool>("layOn",false) ), phiOn( pSet.getUntrackedParameter<bool>("phiOn",false) ), ringOn( pSet.getUntrackedParameter<bool>("ringOn",false) ), bladeOn( pSet.getUntrackedParameter<bool>("bladeOn",false) ), diskOn( pSet.getUntrackedParameter<bool>("diskOn",false) ) { pSet_ = pSet; debug_ = pSet_.getUntrackedParameter<bool>("debug", false); src_ = pSet_.getParameter<edm::InputTag>("src"); clustersrc_ = pSet_.getParameter<edm::InputTag>("clustersrc"); tracksrc_ = pSet_.getParameter<edm::InputTag>("trajectoryInput"); dbe_ = edm::Service<DQMStore>().operator->(); ttrhbuilder_ = pSet_.getParameter<std::string>("TTRHBuilder"); ptminres_= pSet.getUntrackedParameter<double>("PtMinRes",4.0) ; LogInfo("PixelDQM") << "SiPixelTrackResidualSource constructor" << endl; LogInfo ("PixelDQM") << "Mod/Lad/Lay/Phi " << modOn << "/" << ladOn << "/" << layOn << "/" << phiOn << std::endl; LogInfo ("PixelDQM") << "Blade/Disk/Ring" << bladeOn << "/" << diskOn << "/" << ringOn << std::endl; }
| SiPixelTrackResidualSource::~SiPixelTrackResidualSource | ( | ) |
Definition at line 93 of file SiPixelTrackResidualSource.cc.
References theSiPixelStructure.
{
LogInfo("PixelDQM") << "SiPixelTrackResidualSource destructor" << endl;
std::map<uint32_t,SiPixelTrackResidualModule*>::iterator struct_iter;
for (struct_iter = theSiPixelStructure.begin() ; struct_iter != theSiPixelStructure.end() ; struct_iter++){
delete struct_iter->second;
struct_iter->second = 0;
}
}
| void SiPixelTrackResidualSource::analyze | ( | const edm::Event & | iEvent, |
| const edm::EventSetup & | iSetup | ||
| ) | [virtual] |
Implements edm::EDAnalyzer.
Definition at line 515 of file SiPixelTrackResidualSource.cc.
References abs, Reference_intrackfit_cff::barrel, edmNew::DetSetVector< T >::begin(), edm::AssociationMap< Tag >::begin(), bladeOn, DeDxDiscriminatorTools::charge(), clustersrc_, gather_cfg::cout, debug_, DetId::det(), PixelEndcapName::diskName(), diskOn, edmNew::DetSetVector< T >::end(), edm::AssociationMap< Tag >::end(), Reference_intrackfit_cff::endcap, edm::HandleBase::failedToGet(), MonitorElement::Fill(), HcalObjRepresent::Fill(), edmNew::DetSetVector< T >::find(), edm::EventSetup::get(), edm::Event::getByLabel(), corrVsCorr::gX, corrVsCorr::gY, i, TrackerGeometry::idToDet(), reco::TransientTrack::initialFreeState(), reco::TransientTrack::innermostMeasurementState(), edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), edm::ESHandleBase::isValid(), kap, edm::Ref< C, T, F >::key(), ladOn, PixelBarrelName::layerName(), layOn, TrajectoryStateOnSurface::localParameters(), Topology::localPosition(), PV3DBase< T, PVType, FrameType >::mag(), match(), meClChargeNotOnTrack_all, meClChargeNotOnTrack_bpix, meClChargeNotOnTrack_diskm1, meClChargeNotOnTrack_diskm2, meClChargeNotOnTrack_diskp1, meClChargeNotOnTrack_diskp2, meClChargeNotOnTrack_fpix, meClChargeNotOnTrack_layer1, meClChargeNotOnTrack_layer2, meClChargeNotOnTrack_layer3, meClChargeOnTrack_all, meClChargeOnTrack_bpix, meClChargeOnTrack_diskm1, meClChargeOnTrack_diskm2, meClChargeOnTrack_diskp1, meClChargeOnTrack_diskp2, meClChargeOnTrack_fpix, meClChargeOnTrack_layer1, meClChargeOnTrack_layer2, meClChargeOnTrack_layer3, meClPosDisk1mzNotOnTrack, meClPosDisk1mzOnTrack, meClPosDisk1pzNotOnTrack, meClPosDisk1pzOnTrack, meClPosDisk2mzNotOnTrack, meClPosDisk2mzOnTrack, meClPosDisk2pzNotOnTrack, meClPosDisk2pzOnTrack, meClPosLayer1NotOnTrack, meClPosLayer1OnTrack, meClPosLayer2NotOnTrack, meClPosLayer2OnTrack, meClPosLayer3NotOnTrack, meClPosLayer3OnTrack, meClSizeNotOnTrack_all, meClSizeNotOnTrack_bpix, meClSizeNotOnTrack_diskm1, meClSizeNotOnTrack_diskm2, meClSizeNotOnTrack_diskp1, meClSizeNotOnTrack_diskp2, meClSizeNotOnTrack_fpix, meClSizeNotOnTrack_layer1, meClSizeNotOnTrack_layer2, meClSizeNotOnTrack_layer3, meClSizeOnTrack_all, meClSizeOnTrack_bpix, meClSizeOnTrack_diskm1, meClSizeOnTrack_diskm2, meClSizeOnTrack_diskp1, meClSizeOnTrack_diskp2, meClSizeOnTrack_fpix, meClSizeOnTrack_layer1, meClSizeOnTrack_layer2, meClSizeOnTrack_layer3, meClSizeXNotOnTrack_all, meClSizeXNotOnTrack_bpix, meClSizeXNotOnTrack_diskm1, meClSizeXNotOnTrack_diskm2, meClSizeXNotOnTrack_diskp1, meClSizeXNotOnTrack_diskp2, meClSizeXNotOnTrack_fpix, meClSizeXNotOnTrack_layer1, meClSizeXNotOnTrack_layer2, meClSizeXNotOnTrack_layer3, meClSizeXOnTrack_all, meClSizeXOnTrack_bpix, meClSizeXOnTrack_diskm1, meClSizeXOnTrack_diskm2, meClSizeXOnTrack_diskp1, meClSizeXOnTrack_diskp2, meClSizeXOnTrack_fpix, meClSizeXOnTrack_layer1, meClSizeXOnTrack_layer2, meClSizeXOnTrack_layer3, meClSizeYNotOnTrack_all, meClSizeYNotOnTrack_bpix, meClSizeYNotOnTrack_diskm1, meClSizeYNotOnTrack_diskm2, meClSizeYNotOnTrack_diskp1, meClSizeYNotOnTrack_diskp2, meClSizeYNotOnTrack_fpix, meClSizeYNotOnTrack_layer1, meClSizeYNotOnTrack_layer2, meClSizeYNotOnTrack_layer3, meClSizeYOnTrack_all, meClSizeYOnTrack_bpix, meClSizeYOnTrack_diskm1, meClSizeYOnTrack_diskm2, meClSizeYOnTrack_diskp1, meClSizeYOnTrack_diskp2, meClSizeYOnTrack_fpix, meClSizeYOnTrack_layer1, meClSizeYOnTrack_layer2, meClSizeYOnTrack_layer3, meHitProbability, meNClustersNotOnTrack_all, meNClustersNotOnTrack_bpix, meNClustersNotOnTrack_diskm1, meNClustersNotOnTrack_diskm2, meNClustersNotOnTrack_diskp1, meNClustersNotOnTrack_diskp2, meNClustersNotOnTrack_fpix, meNClustersNotOnTrack_layer1, meNClustersNotOnTrack_layer2, meNClustersNotOnTrack_layer3, meNClustersOnTrack_all, meNClustersOnTrack_bpix, meNClustersOnTrack_diskm1, meNClustersOnTrack_diskm2, meNClustersOnTrack_diskp1, meNClustersOnTrack_diskp2, meNClustersOnTrack_fpix, meNClustersOnTrack_layer1, meNClustersOnTrack_layer2, meNClustersOnTrack_layer3, meNofClustersNotOnTrack_, meNofClustersOnTrack_, meNofTracks_, meNofTracksInPixVol_, modOn, LocalTrajectoryParameters::momentum(), reco::Vertex::ndof(), NLowProb, NTotal, PV3DBase< T, PVType, FrameType >::phi(), phi, phiOn, PixelSubdetector::PixelBarrel, GeomDetEnumerators::PixelBarrel, PixelSubdetector::PixelEndcap, GeomDetEnumerators::PixelEndcap, funct::pow(), edm::ESHandle< T >::product(), edm::Handle< T >::product(), ptminres_, TrackerTopology::pxbLayer(), reducedSet, ringOn, edm::AssociationMap< Tag >::size(), edmNew::DetSetVector< T >::size(), PixelGeomDetUnit::specificTopology(), mathSSE::sqrt(), DetId::subdetId(), GeomDet::surface(), funct::tan(), groupFilesInBlocks::temp, theSiPixelStructure, Surface::toGlobal(), align::Tracker, testEve_cfg::tracks, tracksrc_, FreeTrajectoryState::transverseCurvature(), triplets(), ttrhbuilder_, PV3DBase< T, PVType, FrameType >::x(), x, PV3DBase< T, PVType, FrameType >::y(), detailsBasic3DVector::y, PV3DBase< T, PVType, FrameType >::z(), and z.
{
//Retrieve tracker topology from geometry
edm::ESHandle<TrackerTopology> tTopoHandle;
iSetup.get<IdealGeometryRecord>().get(tTopoHandle);
const TrackerTopology* const tTopo = tTopoHandle.product();
// retrieve TrackerGeometry again and MagneticField for use in transforming
// a TrackCandidate's P(ersistent)TrajectoryStateoOnDet (PTSoD) to a TrajectoryStateOnSurface (TSoS)
ESHandle<TrackerGeometry> TG;
iSetup.get<TrackerDigiGeometryRecord>().get(TG);
const TrackerGeometry* theTrackerGeometry = TG.product();
//analytic triplet method to calculate the track residuals in the pixe barrel detector
//--------------------------------------------------------------------
// beam spot:
//
edm::Handle<reco::BeamSpot> rbs;
iEvent.getByLabel( "offlineBeamSpot", rbs );
math::XYZPoint bsP = math::XYZPoint(0,0,0);
if( !rbs.failedToGet() && rbs.isValid() )
{
bsP = math::XYZPoint( rbs->x0(), rbs->y0(), rbs->z0() );
}
//--------------------------------------------------------------------
// primary vertices:
//
edm::Handle<reco::VertexCollection> vertices;
iEvent.getByLabel("offlinePrimaryVertices", vertices );
if( vertices.failedToGet() ) return;
if( !vertices.isValid() ) return;
math::XYZPoint vtxN = math::XYZPoint(0,0,0);
math::XYZPoint vtxP = math::XYZPoint(0,0,0);
double bestNdof = 0;
double maxSumPt = 0;
reco::Vertex bestPvx;
for(reco::VertexCollection::const_iterator iVertex = vertices->begin();
iVertex != vertices->end(); ++iVertex ) {
if( iVertex->ndof() > bestNdof ) {
bestNdof = iVertex->ndof();
vtxN = math::XYZPoint( iVertex->x(), iVertex->y(), iVertex->z() );
}//ndof
if( iVertex->p4().pt() > maxSumPt ) {
maxSumPt = iVertex->p4().pt();
vtxP = math::XYZPoint( iVertex->x(), iVertex->y(), iVertex->z() );
bestPvx = *iVertex;
}//sumpt
}//vertex
if( maxSumPt < 1 ) return;
if( maxSumPt < 1 ) vtxP = vtxN;
//---------------------------------------------
//get Tracks
//
edm:: Handle<reco::TrackCollection> TracksForRes;
iEvent.getByLabel( "generalTracks", TracksForRes );
//
// transient track builder, needs B-field from data base (global tag in .py)
//
edm::ESHandle<TransientTrackBuilder> theB;
iSetup.get<TransientTrackRecord>().get( "TransientTrackBuilder", theB );
//get the TransienTrackingRecHitBuilder needed for extracting the global position of the hits in the pixel
edm::ESHandle<TransientTrackingRecHitBuilder> theTrackerRecHitBuilder;
iSetup.get<TransientRecHitRecord>().get(ttrhbuilder_,theTrackerRecHitBuilder);
//check that tracks are valid
if( TracksForRes.failedToGet() ) return;
if( !TracksForRes.isValid() ) return;
//get tracker geometry
edm::ESHandle<TrackerGeometry> pDD;
iSetup.get<TrackerDigiGeometryRecord>().get(pDD);
if( !pDD.isValid() ) {
cout << "Unable to find TrackerDigiGeometry. Return\n";
return;
}
int kk = -1;
//----------------------------------------------------------------------------
// Residuals:
//
for( reco::TrackCollection::const_iterator iTrack = TracksForRes->begin();
iTrack != TracksForRes->end(); ++iTrack ) {
//count
kk++;
//Calculate minimal track pt before curling
// cpt = cqRB = 0.3*R[m]*B[T] = 1.14*R[m] for B=3.8T
// D = 2R = 2*pt/1.14
// calo: D = 1.3 m => pt = 0.74 GeV/c
double pt = iTrack->pt();
if( pt < 0.75 ) continue;// curls up
if( abs( iTrack->dxy(vtxP) ) > 5*iTrack->dxyError() ) continue; // not prompt
double charge = iTrack->charge();
reco::TransientTrack tTrack = theB->build(*iTrack);
//get curvature of the track, needed for the residuals
double kap = tTrack.initialFreeState().transverseCurvature();
//needed for the TransienTrackingRecHitBuilder
TrajectoryStateOnSurface initialTSOS = tTrack.innermostMeasurementState();
if( iTrack->extra().isNonnull() &&iTrack->extra().isAvailable() ){
double x1 = 0;
double y1 = 0;
double z1 = 0;
double x2 = 0;
double y2 = 0;
double z2 = 0;
double x3 = 0;
double y3 = 0;
double z3 = 0;
int n1 = 0;
int n2 = 0;
int n3 = 0;
//for saving the pixel barrel hits
vector<TransientTrackingRecHit::RecHitPointer> GoodPixBarrelHits;
//looping through the RecHits of the track
for( trackingRecHit_iterator irecHit = iTrack->recHitsBegin();
irecHit != iTrack->recHitsEnd(); ++irecHit){
if( (*irecHit)->isValid() ){
DetId detId = (*irecHit)->geographicalId();
// enum Detector { Tracker=1, Muon=2, Ecal=3, Hcal=4, Calo=5 };
if( detId.det() != 1 ){
if(debug_){
cout << "rec hit ID = " << detId.det() << " not in tracker!?!?\n";
}
continue;
}
uint32_t subDet = detId.subdetId();
// enum SubDetector{ PixelBarrel=1, PixelEndcap=2 };
// enum SubDetector{ TIB=3, TID=4, TOB=5, TEC=6 };
TransientTrackingRecHit::RecHitPointer trecHit = theTrackerRecHitBuilder->build( &*(*irecHit), initialTSOS);
double gX = trecHit->globalPosition().x();
double gY = trecHit->globalPosition().y();
double gZ = trecHit->globalPosition().z();
if( subDet == PixelSubdetector::PixelBarrel ) {
int ilay = tTopo->pxbLayer(detId);
if( ilay == 1 ){
n1++;
x1 = gX;
y1 = gY;
z1 = gZ;
GoodPixBarrelHits.push_back((trecHit));
}//PXB1
if( ilay == 2 ){
n2++;
x2 = gX;
y2 = gY;
z2 = gZ;
GoodPixBarrelHits.push_back((trecHit));
}//PXB2
if( ilay == 3 ){
n3++;
x3 = gX;
y3 = gY;
z3 = gZ;
GoodPixBarrelHits.push_back((trecHit));
}
}//PXB
}//valid
}//loop rechits
//CS extra plots
if( n1+n2+n3 == 3 && n1*n2*n3 > 0) {
for( unsigned int i = 0; i < GoodPixBarrelHits.size(); i++){
if( GoodPixBarrelHits[i]->isValid() ){
DetId detId = GoodPixBarrelHits[i]->geographicalId().rawId();
int ilay = tTopo->pxbLayer(detId);
if(pt > ptminres_){
double dca2 = 0.0, dz2=0.0;
double ptsig = pt;
if(charge<0.) ptsig = -pt;
//Filling the histograms in modules
MeasurementPoint Test;
MeasurementPoint Test2;
Test=MeasurementPoint(0,0);
Test2=MeasurementPoint(0,0);
Measurement2DVector residual;
if( ilay == 1 ){
triplets(x2,y2,z2,x1,y1,z1,x3,y3,z3,ptsig,dca2,dz2, kap);
Test=MeasurementPoint(dca2*1E4,dz2*1E4);
residual=Test-Test2;
}
if( ilay == 2 ){
triplets(x1,y1,z1,x2,y2,z2,x3,y3,z3,ptsig,dca2,dz2, kap);
Test=MeasurementPoint(dca2*1E4,dz2*1E4);
residual=Test-Test2;
}
if( ilay == 3 ){
triplets(x1,y1,z1,x3,y3,z3,x2,y2,z2,ptsig,dca2,dz2, kap);
Test=MeasurementPoint(dca2*1E4,dz2*1E4);
residual=Test-Test2;
}
// fill the residual histograms
std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.find(detId);
if (pxd!=theSiPixelStructure.end()) (*pxd).second->fill(residual, reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);
}//three hits
}//is valid
}//rechits loop
}//pt 4
}
}//-----Tracks
//get trajectories
edm::Handle<std::vector<Trajectory> > trajCollectionHandle;
iEvent.getByLabel(tracksrc_,trajCollectionHandle);
const std::vector<Trajectory> trajColl = *(trajCollectionHandle.product());
//get tracks
edm::Handle<std::vector<reco::Track> > trackCollectionHandle;
iEvent.getByLabel(tracksrc_,trackCollectionHandle);
const std::vector<reco::Track> trackColl = *(trackCollectionHandle.product());
//get the map
edm::Handle<TrajTrackAssociationCollection> match;
iEvent.getByLabel(tracksrc_,match);
const TrajTrackAssociationCollection ttac = *(match.product());
// get clusters
edm::Handle< edmNew::DetSetVector<SiPixelCluster> > clusterColl;
iEvent.getByLabel( clustersrc_, clusterColl );
const edmNew::DetSetVector<SiPixelCluster> clustColl = *(clusterColl.product());
if(debug_){
std::cout << "Trajectories\t : " << trajColl.size() << std::endl;
std::cout << "recoTracks \t : " << trackColl.size() << std::endl;
std::cout << "Map entries \t : " << ttac.size() << std::endl;
}
std::set<SiPixelCluster> clusterSet;
TrajectoryStateCombiner tsoscomb;
int tracks=0, pixeltracks=0, bpixtracks=0, fpixtracks=0;
int trackclusters=0, barreltrackclusters=0, endcaptrackclusters=0;
int otherclusters=0, barrelotherclusters=0, endcapotherclusters=0;
//Loop over map entries
for(TrajTrackAssociationCollection::const_iterator it = ttac.begin();it != ttac.end(); ++it){
const edm::Ref<std::vector<Trajectory> > traj_iterator = it->key;
// Trajectory Map, extract Trajectory for this track
reco::TrackRef trackref = it->val;
tracks++;
bool isBpixtrack = false, isFpixtrack = false, crossesPixVol=false;
//find out whether track crosses pixel fiducial volume (for cosmic tracks)
double d0 = (*trackref).d0(), dz = (*trackref).dz();
if(abs(d0)<15 && abs(dz)<50) crossesPixVol = true;
std::vector<TrajectoryMeasurement> tmeasColl =traj_iterator->measurements();
std::vector<TrajectoryMeasurement>::const_iterator tmeasIt;
//loop on measurements to find out whether there are bpix and/or fpix hits
for(tmeasIt = tmeasColl.begin();tmeasIt!=tmeasColl.end();tmeasIt++){
if(! tmeasIt->updatedState().isValid()) continue;
TransientTrackingRecHit::ConstRecHitPointer testhit = tmeasIt->recHit();
if(! testhit->isValid() || testhit->geographicalId().det() != DetId::Tracker) continue;
uint testSubDetID = (testhit->geographicalId().subdetId());
if(testSubDetID==PixelSubdetector::PixelBarrel) isBpixtrack = true;
if(testSubDetID==PixelSubdetector::PixelEndcap) isFpixtrack = true;
}//end loop on measurements
if(isBpixtrack) {
bpixtracks++;
if(debug_) std::cout << "bpixtrack\n";
}
if(isFpixtrack) {
fpixtracks++;
if(debug_) std::cout << "fpixtrack\n";
}
if(isBpixtrack || isFpixtrack){
pixeltracks++;
if(crossesPixVol) meNofTracksInPixVol_->Fill(0,1);
std::vector<TrajectoryMeasurement> tmeasColl = traj_iterator->measurements();
for(std::vector<TrajectoryMeasurement>::const_iterator tmeasIt = tmeasColl.begin(); tmeasIt!=tmeasColl.end(); tmeasIt++){
if(! tmeasIt->updatedState().isValid()) continue;
TrajectoryStateOnSurface tsos = tsoscomb( tmeasIt->forwardPredictedState(), tmeasIt->backwardPredictedState() );
TransientTrackingRecHit::ConstRecHitPointer hit = tmeasIt->recHit();
if(! hit->isValid() || hit->geographicalId().det() != DetId::Tracker ) {
continue;
} else {
// //residual
const DetId & hit_detId = hit->geographicalId();
//uint IntRawDetID = (hit_detId.rawId());
uint IntSubDetID = (hit_detId.subdetId());
if(IntSubDetID == 0 ) continue; // don't look at SiStrip hits!
// get the enclosed persistent hit
const TrackingRecHit *persistentHit = hit->hit();
// check if it's not null, and if it's a valid pixel hit
if ((persistentHit != 0) && (typeid(*persistentHit) == typeid(SiPixelRecHit))) {
// tell the C++ compiler that the hit is a pixel hit
const SiPixelRecHit* pixhit = static_cast<const SiPixelRecHit*>( hit->hit() );
//Hit probability:
float hit_prob = -1.;
if(pixhit->hasFilledProb()){
hit_prob = pixhit->clusterProbability(0);
//std::cout<<"HITPROB= "<<hit_prob<<std::endl;
if(hit_prob<pow(10.,-15.)) NLowProb++;
NTotal++;
if(NTotal>0) meHitProbability->Fill(float(NLowProb/NTotal));
}
// get the edm::Ref to the cluster
edm::Ref<edmNew::DetSetVector<SiPixelCluster>, SiPixelCluster> const& clust = (*pixhit).cluster();
// check if the ref is not null
if (clust.isNonnull()) {
//define tracker and pixel geometry and topology
const TrackerGeometry& theTracker(*theTrackerGeometry);
const PixelGeomDetUnit* theGeomDet = static_cast<const PixelGeomDetUnit*> (theTracker.idToDet(hit_detId) );
//test if PixelGeomDetUnit exists
if(theGeomDet == 0) {
if(debug_) std::cout << "NO THEGEOMDET\n";
continue;
}
const PixelTopology * topol = &(theGeomDet->specificTopology());
//fill histograms for clusters on tracks
//correct SiPixelTrackResidualModule
std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.find((*hit).geographicalId().rawId());
//CHARGE CORRECTION (for track impact angle)
// calculate alpha and beta from cluster position
LocalTrajectoryParameters ltp = tsos.localParameters();
LocalVector localDir = ltp.momentum()/ltp.momentum().mag();
float clust_alpha = atan2(localDir.z(), localDir.x());
float clust_beta = atan2(localDir.z(), localDir.y());
double corrCharge = clust->charge() * sqrt( 1.0 / ( 1.0/pow( tan(clust_alpha), 2 ) +
1.0/pow( tan(clust_beta ), 2 ) +
1.0 )
)/1000.;
if (pxd!=theSiPixelStructure.end()) (*pxd).second->fill((*clust), true, corrCharge, reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);
trackclusters++;
//CORR CHARGE
meClChargeOnTrack_all->Fill(corrCharge);
meClSizeOnTrack_all->Fill((*clust).size());
meClSizeXOnTrack_all->Fill((*clust).sizeX());
meClSizeYOnTrack_all->Fill((*clust).sizeY());
clusterSet.insert(*clust);
//find cluster global position (rphi, z)
// get cluster center of gravity (of charge)
float xcenter = clust->x();
float ycenter = clust->y();
// get the cluster position in local coordinates (cm)
LocalPoint clustlp = topol->localPosition( MeasurementPoint(xcenter, ycenter) );
// get the cluster position in global coordinates (cm)
GlobalPoint clustgp = theGeomDet->surface().toGlobal( clustlp );
//find location of hit (barrel or endcap, same for cluster)
bool barrel = DetId((*hit).geographicalId()).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
bool endcap = DetId((*hit).geographicalId()).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
if(barrel) {
barreltrackclusters++;
//CORR CHARGE
meClChargeOnTrack_bpix->Fill(corrCharge);
meClSizeOnTrack_bpix->Fill((*clust).size());
meClSizeXOnTrack_bpix->Fill((*clust).sizeX());
meClSizeYOnTrack_bpix->Fill((*clust).sizeY());
uint32_t DBlayer = PixelBarrelName(DetId((*hit).geographicalId())).layerName();
float phi = clustgp.phi();
float z = clustgp.z();
switch(DBlayer){
case 1: {
meClPosLayer1OnTrack->Fill(z,phi);
meClChargeOnTrack_layer1->Fill(corrCharge);
meClSizeOnTrack_layer1->Fill((*clust).size());
meClSizeXOnTrack_layer1->Fill((*clust).sizeX());
meClSizeYOnTrack_layer1->Fill((*clust).sizeY());
break;
}
case 2: {
meClPosLayer2OnTrack->Fill(z,phi);
meClChargeOnTrack_layer2->Fill(corrCharge);
meClSizeOnTrack_layer2->Fill((*clust).size());
meClSizeXOnTrack_layer2->Fill((*clust).sizeX());
meClSizeYOnTrack_layer2->Fill((*clust).sizeY());
break;
}
case 3: {
meClPosLayer3OnTrack->Fill(z,phi);
meClChargeOnTrack_layer3->Fill(corrCharge);
meClSizeOnTrack_layer3->Fill((*clust).size());
meClSizeXOnTrack_layer3->Fill((*clust).sizeX());
meClSizeYOnTrack_layer3->Fill((*clust).sizeY());
break;
}
}
}
if(endcap) {
endcaptrackclusters++;
//CORR CHARGE
meClChargeOnTrack_fpix->Fill(corrCharge);
meClSizeOnTrack_fpix->Fill((*clust).size());
meClSizeXOnTrack_fpix->Fill((*clust).sizeX());
meClSizeYOnTrack_fpix->Fill((*clust).sizeY());
uint32_t DBdisk = PixelEndcapName(DetId((*hit).geographicalId())).diskName();
float x = clustgp.x();
float y = clustgp.y();
float z = clustgp.z();
if(z>0){
if(DBdisk==1) {
meClPosDisk1pzOnTrack->Fill(x,y);
meClChargeOnTrack_diskp1->Fill(corrCharge);
meClSizeOnTrack_diskp1->Fill((*clust).size());
meClSizeXOnTrack_diskp1->Fill((*clust).sizeX());
meClSizeYOnTrack_diskp1->Fill((*clust).sizeY());
}
if(DBdisk==2) {
meClPosDisk2pzOnTrack->Fill(x,y);
meClChargeOnTrack_diskp2->Fill(corrCharge);
meClSizeOnTrack_diskp2->Fill((*clust).size());
meClSizeXOnTrack_diskp2->Fill((*clust).sizeX());
meClSizeYOnTrack_diskp2->Fill((*clust).sizeY());
}
}
else{
if(DBdisk==1) {
meClPosDisk1mzOnTrack->Fill(x,y);
meClChargeOnTrack_diskm1->Fill(corrCharge);
meClSizeOnTrack_diskm1->Fill((*clust).size());
meClSizeXOnTrack_diskm1->Fill((*clust).sizeX());
meClSizeYOnTrack_diskm1->Fill((*clust).sizeY());
}
if(DBdisk==2) {
meClPosDisk2mzOnTrack->Fill(x,y);
meClChargeOnTrack_diskm2->Fill(corrCharge);
meClSizeOnTrack_diskm2->Fill((*clust).size());
meClSizeXOnTrack_diskm2->Fill((*clust).sizeX());
meClSizeYOnTrack_diskm2->Fill((*clust).sizeY());
}
}
}
}//end if (cluster exists)
}//end if (persistent hit exists and is pixel hit)
}//end of else
}//end for (all traj measurements of pixeltrack)
}//end if (is pixeltrack)
else {
if(debug_) std::cout << "no pixeltrack:\n";
if(crossesPixVol) meNofTracksInPixVol_->Fill(1,1);
}
}//end loop on map entries
//find clusters that are NOT on track
//edmNew::DetSet<SiPixelCluster>::const_iterator di;
if(debug_) std::cout << "clusters not on track: (size " << clustColl.size() << ") ";
for(TrackerGeometry::DetContainer::const_iterator it = TG->dets().begin(); it != TG->dets().end(); it++){
//if(dynamic_cast<PixelGeomDetUnit*>((*it))!=0){
DetId detId = (*it)->geographicalId();
if(detId>=302055684 && detId<=352477708){ // make sure it's a Pixel module WITHOUT using dynamic_cast!
int nofclOnTrack = 0, nofclOffTrack=0;
uint32_t DBlayer=10, DBdisk=10;
float z=0.;
//set layer/disk
if(DetId(detId).subdetId() == 1) { // Barrel module
DBlayer = PixelBarrelName(DetId(detId)).layerName();
}
if(DetId(detId).subdetId() == 2){ // Endcap module
DBdisk = PixelEndcapName(DetId(detId )).diskName();
}
edmNew::DetSetVector<SiPixelCluster>::const_iterator isearch = clustColl.find(detId);
if( isearch != clustColl.end() ) { // Not an empty iterator
edmNew::DetSet<SiPixelCluster>::const_iterator di;
for(di=isearch->begin(); di!=isearch->end(); di++){
unsigned int temp = clusterSet.size();
clusterSet.insert(*di);
//check if cluster is off track
if(clusterSet.size()>temp) {
otherclusters++;
nofclOffTrack++;
//fill histograms for clusters off tracks
//correct SiPixelTrackResidualModule
std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.find((*it)->geographicalId().rawId());
if (pxd!=theSiPixelStructure.end()) (*pxd).second->fill((*di), false, -1., reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);
meClSizeNotOnTrack_all->Fill((*di).size());
meClSizeXNotOnTrack_all->Fill((*di).sizeX());
meClSizeYNotOnTrack_all->Fill((*di).sizeY());
meClChargeNotOnTrack_all->Fill((*di).charge()/1000);
//find cluster global position (rphi, z) get cluster
//define tracker and pixel geometry and topology
const TrackerGeometry& theTracker(*theTrackerGeometry);
const PixelGeomDetUnit* theGeomDet = static_cast<const PixelGeomDetUnit*> (theTracker.idToDet(detId) );
//test if PixelGeomDetUnit exists
if(theGeomDet == 0) {
if(debug_) std::cout << "NO THEGEOMDET\n";
continue;
}
const PixelTopology * topol = &(theGeomDet->specificTopology());
//center of gravity (of charge)
float xcenter = di->x();
float ycenter = di->y();
// get the cluster position in local coordinates (cm)
LocalPoint clustlp = topol->localPosition( MeasurementPoint(xcenter, ycenter) );
// get the cluster position in global coordinates (cm)
GlobalPoint clustgp = theGeomDet->surface().toGlobal( clustlp );
//barrel
if(DetId(detId).subdetId() == 1) {
meClSizeNotOnTrack_bpix->Fill((*di).size());
meClSizeXNotOnTrack_bpix->Fill((*di).sizeX());
meClSizeYNotOnTrack_bpix->Fill((*di).sizeY());
meClChargeNotOnTrack_bpix->Fill((*di).charge()/1000);
barrelotherclusters++;
//DBlayer = PixelBarrelName(DetId(detId)).layerName();
float phi = clustgp.phi();
//float r = clustgp.perp();
z = clustgp.z();
switch(DBlayer){
case 1: {
meClPosLayer1NotOnTrack->Fill(z,phi);
meClSizeNotOnTrack_layer1->Fill((*di).size());
meClSizeXNotOnTrack_layer1->Fill((*di).sizeX());
meClSizeYNotOnTrack_layer1->Fill((*di).sizeY());
meClChargeNotOnTrack_layer1->Fill((*di).charge()/1000);
break;
}
case 2: {
meClPosLayer2NotOnTrack->Fill(z,phi);
meClSizeNotOnTrack_layer2->Fill((*di).size());
meClSizeXNotOnTrack_layer2->Fill((*di).sizeX());
meClSizeYNotOnTrack_layer2->Fill((*di).sizeY());
meClChargeNotOnTrack_layer2->Fill((*di).charge()/1000);
break;
}
case 3: {
meClPosLayer3NotOnTrack->Fill(z,phi);
meClSizeNotOnTrack_layer3->Fill((*di).size());
meClSizeXNotOnTrack_layer3->Fill((*di).sizeX());
meClSizeYNotOnTrack_layer3->Fill((*di).sizeY());
meClChargeNotOnTrack_layer3->Fill((*di).charge()/1000);
break;
}
}
}
//endcap
if(DetId(detId).subdetId() == 2) {
meClSizeNotOnTrack_fpix->Fill((*di).size());
meClSizeXNotOnTrack_fpix->Fill((*di).sizeX());
meClSizeYNotOnTrack_fpix->Fill((*di).sizeY());
meClChargeNotOnTrack_fpix->Fill((*di).charge()/1000);
endcapotherclusters++;
//DBdisk = PixelEndcapName(DetId(detId )).diskName();
float x = clustgp.x();
float y = clustgp.y();
z = clustgp.z();
if(z>0){
if(DBdisk==1) {
meClPosDisk1pzNotOnTrack->Fill(x,y);
meClSizeNotOnTrack_diskp1->Fill((*di).size());
meClSizeXNotOnTrack_diskp1->Fill((*di).sizeX());
meClSizeYNotOnTrack_diskp1->Fill((*di).sizeY());
meClChargeNotOnTrack_diskp1->Fill((*di).charge()/1000);
}
if(DBdisk==2) {
meClPosDisk2pzNotOnTrack->Fill(x,y);
meClSizeNotOnTrack_diskp2->Fill((*di).size());
meClSizeXNotOnTrack_diskp2->Fill((*di).sizeX());
meClSizeYNotOnTrack_diskp2->Fill((*di).sizeY());
meClChargeNotOnTrack_diskp2->Fill((*di).charge()/1000);
}
}
else{
if(DBdisk==1) {
meClPosDisk1mzNotOnTrack->Fill(x,y);
meClSizeNotOnTrack_diskm1->Fill((*di).size());
meClSizeXNotOnTrack_diskm1->Fill((*di).sizeX());
meClSizeYNotOnTrack_diskm1->Fill((*di).sizeY());
meClChargeNotOnTrack_diskm1->Fill((*di).charge()/1000);
}
if(DBdisk==2) {
meClPosDisk2mzNotOnTrack->Fill(x,y);
meClSizeNotOnTrack_diskm2->Fill((*di).size());
meClSizeXNotOnTrack_diskm2->Fill((*di).sizeX());
meClSizeYNotOnTrack_diskm2->Fill((*di).sizeY());
meClChargeNotOnTrack_diskm2->Fill((*di).charge()/1000);
}
}
}
}// end "if cluster off track"
else {
nofclOnTrack++;
if(z == 0 && DBdisk != 10){
//find cluster global position (rphi, z) get cluster
//define tracker and pixel geometry and topology
const TrackerGeometry& theTracker(*theTrackerGeometry);
const PixelGeomDetUnit* theGeomDet = static_cast<const PixelGeomDetUnit*> (theTracker.idToDet(detId) );
//test if PixelGeomDetUnit exists
if(theGeomDet == 0) {
if(debug_) std::cout << "NO THEGEOMDET\n";
continue;
}
const PixelTopology * topol = &(theGeomDet->specificTopology());
//center of gravity (of charge)
float xcenter = di->x();
float ycenter = di->y();
// get the cluster position in local coordinates (cm)
LocalPoint clustlp = topol->localPosition( MeasurementPoint(xcenter, ycenter) );
// get the cluster position in global coordinates (cm)
GlobalPoint clustgp = theGeomDet->surface().toGlobal( clustlp );
z = clustgp.z();
}
}
}
}
//++ fill the number of clusters on a module
std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.find((*it)->geographicalId().rawId());
if (pxd!=theSiPixelStructure.end()) (*pxd).second->nfill(nofclOnTrack, nofclOffTrack, reducedSet, modOn, ladOn, layOn, phiOn, bladeOn, diskOn, ringOn);
if(nofclOnTrack!=0) meNClustersOnTrack_all->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_all->Fill(nofclOffTrack);
//barrel
if(DetId(detId).subdetId() == 1){
if(nofclOnTrack!=0) meNClustersOnTrack_bpix->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_bpix->Fill(nofclOffTrack);
//DBlayer = PixelBarrelName(DetId(detId)).layerName();
switch(DBlayer){
case 1: {
if(nofclOnTrack!=0) meNClustersOnTrack_layer1->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_layer1->Fill(nofclOffTrack); break;
}
case 2: {
if(nofclOnTrack!=0) meNClustersOnTrack_layer2->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_layer2->Fill(nofclOffTrack); break;
}
case 3: {
if(nofclOnTrack!=0) meNClustersOnTrack_layer3->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_layer3->Fill(nofclOffTrack); break;
}
}
}//end barrel
//endcap
if(DetId(detId).subdetId() == 2) {
//DBdisk = PixelEndcapName(DetId(detId )).diskName();
//z = clustgp.z();
if(nofclOnTrack!=0) meNClustersOnTrack_fpix->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_fpix->Fill(nofclOffTrack);
if(z>0){
if(DBdisk==1) {
if(nofclOnTrack!=0) meNClustersOnTrack_diskp1->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_diskp1->Fill(nofclOffTrack);
}
if(DBdisk==2) {
if(nofclOnTrack!=0) meNClustersOnTrack_diskp2->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_diskp2->Fill(nofclOffTrack);
}
}
if(z<0){
if(DBdisk==1) {
if(nofclOnTrack!=0) meNClustersOnTrack_diskm1->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_diskm1->Fill(nofclOffTrack);
}
if(DBdisk==2) {
if(nofclOnTrack!=0) meNClustersOnTrack_diskm2->Fill(nofclOnTrack);
if(nofclOffTrack!=0) meNClustersNotOnTrack_diskm2->Fill(nofclOffTrack);
}
}
}
}//end if it's a Pixel module
}//end for loop over tracker detector geometry modules
if(trackclusters>0) (meNofClustersOnTrack_)->Fill(0,trackclusters);
if(barreltrackclusters>0)(meNofClustersOnTrack_)->Fill(1,barreltrackclusters);
if(endcaptrackclusters>0)(meNofClustersOnTrack_)->Fill(2,endcaptrackclusters);
if(otherclusters>0)(meNofClustersNotOnTrack_)->Fill(0,otherclusters);
if(barrelotherclusters>0)(meNofClustersNotOnTrack_)->Fill(1,barrelotherclusters);
if(endcapotherclusters>0)(meNofClustersNotOnTrack_)->Fill(2,endcapotherclusters);
if(tracks>0)(meNofTracks_)->Fill(0,tracks);
if(pixeltracks>0)(meNofTracks_)->Fill(1,pixeltracks);
if(bpixtracks>0)(meNofTracks_)->Fill(2,bpixtracks);
if(fpixtracks>0)(meNofTracks_)->Fill(3,fpixtracks);
}
| void SiPixelTrackResidualSource::beginJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 103 of file SiPixelTrackResidualSource.cc.
| void SiPixelTrackResidualSource::beginRun | ( | const edm::Run & | r, |
| edm::EventSetup const & | iSetup | ||
| ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 111 of file SiPixelTrackResidualSource.cc.
References bladeOn, DQMStore::book1D(), DQMStore::book2D(), clustersrc_, dbe_, debug_, diskOn, Exception, firstRun, edm::EventSetup::get(), edm::InputTag::label(), ladOn, layOn, meClChargeNotOnTrack_all, meClChargeNotOnTrack_bpix, meClChargeNotOnTrack_diskm1, meClChargeNotOnTrack_diskm2, meClChargeNotOnTrack_diskp1, meClChargeNotOnTrack_diskp2, meClChargeNotOnTrack_fpix, meClChargeNotOnTrack_layer1, meClChargeNotOnTrack_layer2, meClChargeNotOnTrack_layer3, meClChargeOnTrack_all, meClChargeOnTrack_bpix, meClChargeOnTrack_diskm1, meClChargeOnTrack_diskm2, meClChargeOnTrack_diskp1, meClChargeOnTrack_diskp2, meClChargeOnTrack_fpix, meClChargeOnTrack_layer1, meClChargeOnTrack_layer2, meClChargeOnTrack_layer3, meClPosDisk1mzNotOnTrack, meClPosDisk1mzOnTrack, meClPosDisk1pzNotOnTrack, meClPosDisk1pzOnTrack, meClPosDisk2mzNotOnTrack, meClPosDisk2mzOnTrack, meClPosDisk2pzNotOnTrack, meClPosDisk2pzOnTrack, meClPosLayer1NotOnTrack, meClPosLayer1OnTrack, meClPosLayer2NotOnTrack, meClPosLayer2OnTrack, meClPosLayer3NotOnTrack, meClPosLayer3OnTrack, meClSizeNotOnTrack_all, meClSizeNotOnTrack_bpix, meClSizeNotOnTrack_diskm1, meClSizeNotOnTrack_diskm2, meClSizeNotOnTrack_diskp1, meClSizeNotOnTrack_diskp2, meClSizeNotOnTrack_fpix, meClSizeNotOnTrack_layer1, meClSizeNotOnTrack_layer2, meClSizeNotOnTrack_layer3, meClSizeOnTrack_all, meClSizeOnTrack_bpix, meClSizeOnTrack_diskm1, meClSizeOnTrack_diskm2, meClSizeOnTrack_diskp1, meClSizeOnTrack_diskp2, meClSizeOnTrack_fpix, meClSizeOnTrack_layer1, meClSizeOnTrack_layer2, meClSizeOnTrack_layer3, meClSizeXNotOnTrack_all, meClSizeXNotOnTrack_bpix, meClSizeXNotOnTrack_diskm1, meClSizeXNotOnTrack_diskm2, meClSizeXNotOnTrack_diskp1, meClSizeXNotOnTrack_diskp2, meClSizeXNotOnTrack_fpix, meClSizeXNotOnTrack_layer1, meClSizeXNotOnTrack_layer2, meClSizeXNotOnTrack_layer3, meClSizeXOnTrack_all, meClSizeXOnTrack_bpix, meClSizeXOnTrack_diskm1, meClSizeXOnTrack_diskm2, meClSizeXOnTrack_diskp1, meClSizeXOnTrack_diskp2, meClSizeXOnTrack_fpix, meClSizeXOnTrack_layer1, meClSizeXOnTrack_layer2, meClSizeXOnTrack_layer3, meClSizeYNotOnTrack_all, meClSizeYNotOnTrack_bpix, meClSizeYNotOnTrack_diskm1, meClSizeYNotOnTrack_diskm2, meClSizeYNotOnTrack_diskp1, meClSizeYNotOnTrack_diskp2, meClSizeYNotOnTrack_fpix, meClSizeYNotOnTrack_layer1, meClSizeYNotOnTrack_layer2, meClSizeYNotOnTrack_layer3, meClSizeYOnTrack_all, meClSizeYOnTrack_bpix, meClSizeYOnTrack_diskm1, meClSizeYOnTrack_diskm2, meClSizeYOnTrack_diskp1, meClSizeYOnTrack_diskp2, meClSizeYOnTrack_fpix, meClSizeYOnTrack_layer1, meClSizeYOnTrack_layer2, meClSizeYOnTrack_layer3, meHitProbability, meNClustersNotOnTrack_all, meNClustersNotOnTrack_bpix, meNClustersNotOnTrack_diskm1, meNClustersNotOnTrack_diskm2, meNClustersNotOnTrack_diskp1, meNClustersNotOnTrack_diskp2, meNClustersNotOnTrack_fpix, meNClustersNotOnTrack_layer1, meNClustersNotOnTrack_layer2, meNClustersNotOnTrack_layer3, meNClustersOnTrack_all, meNClustersOnTrack_bpix, meNClustersOnTrack_diskm1, meNClustersOnTrack_diskm2, meNClustersOnTrack_diskp1, meNClustersOnTrack_diskp2, meNClustersOnTrack_fpix, meNClustersOnTrack_layer1, meNClustersOnTrack_layer2, meNClustersOnTrack_layer3, meNofClustersNotOnTrack_, meNofClustersOnTrack_, meNofTracks_, meNofTracksInPixVol_, meSubdetResidualX, meSubdetResidualY, modOn, python::rootplot::argparse::module, phiOn, pSet_, reducedSet, ringOn, alignCSCRings::s, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), DQMStore::setCurrentFolder(), SiPixelFolderOrganizer::setModuleFolder(), theSiPixelStructure, and tracksrc_.
{
LogInfo("PixelDQM") << "SiPixelTrackResidualSource beginRun()" << endl;
if(firstRun){
// retrieve TrackerGeometry for pixel dets
edm::ESHandle<TrackerGeometry> TG;
iSetup.get<TrackerDigiGeometryRecord>().get(TG);
if (debug_) LogVerbatim("PixelDQM") << "TrackerGeometry "<< &(*TG) <<" size is "<< TG->dets().size() << endl;
// build theSiPixelStructure with the pixel barrel and endcap dets from TrackerGeometry
for (TrackerGeometry::DetContainer::const_iterator pxb = TG->detsPXB().begin();
pxb!=TG->detsPXB().end(); pxb++) {
if (dynamic_cast<PixelGeomDetUnit*>((*pxb))!=0) {
SiPixelTrackResidualModule* module = new SiPixelTrackResidualModule((*pxb)->geographicalId().rawId());
theSiPixelStructure.insert(pair<uint32_t, SiPixelTrackResidualModule*>((*pxb)->geographicalId().rawId(), module));
}
}
for (TrackerGeometry::DetContainer::const_iterator pxf = TG->detsPXF().begin();
pxf!=TG->detsPXF().end(); pxf++) {
if (dynamic_cast<PixelGeomDetUnit*>((*pxf))!=0) {
SiPixelTrackResidualModule* module = new SiPixelTrackResidualModule((*pxf)->geographicalId().rawId());
theSiPixelStructure.insert(pair<uint32_t, SiPixelTrackResidualModule*>((*pxf)->geographicalId().rawId(), module));
}
}
LogInfo("PixelDQM") << "SiPixelStructure size is " << theSiPixelStructure.size() << endl;
// book residual histograms in theSiPixelFolder - one (x,y) pair of histograms per det
SiPixelFolderOrganizer theSiPixelFolder;
for (std::map<uint32_t, SiPixelTrackResidualModule*>::iterator pxd = theSiPixelStructure.begin();
pxd!=theSiPixelStructure.end(); pxd++) {
if(modOn){
if (theSiPixelFolder.setModuleFolder((*pxd).first)) (*pxd).second->book(pSet_,reducedSet);
else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Folder Creation Failed! ";
}
if(ladOn){
if (theSiPixelFolder.setModuleFolder((*pxd).first,1)) {
(*pxd).second->book(pSet_,reducedSet,1);
}
else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource ladder Folder Creation Failed! ";
}
if(layOn){
if (theSiPixelFolder.setModuleFolder((*pxd).first,2)) (*pxd).second->book(pSet_,reducedSet,2);
else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource layer Folder Creation Failed! ";
}
if(phiOn){
if (theSiPixelFolder.setModuleFolder((*pxd).first,3)) (*pxd).second->book(pSet_,reducedSet,3);
else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource phi Folder Creation Failed! ";
}
if(bladeOn){
if (theSiPixelFolder.setModuleFolder((*pxd).first,4)) (*pxd).second->book(pSet_,reducedSet,4);
else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Blade Folder Creation Failed! ";
}
if(diskOn){
if (theSiPixelFolder.setModuleFolder((*pxd).first,5)) (*pxd).second->book(pSet_,reducedSet,5);
else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Disk Folder Creation Failed! ";
}
if(ringOn){
if (theSiPixelFolder.setModuleFolder((*pxd).first,6)) (*pxd).second->book(pSet_,reducedSet,6);
else throw cms::Exception("LogicError") << "SiPixelTrackResidualSource Ring Folder Creation Failed! ";
}
}
// edm::InputTag tracksrc = pSet_.getParameter<edm::InputTag>("trajectoryInput");
// edm::InputTag clustersrc = pSet_.getParameter<edm::InputTag>("clustersrc");
//number of tracks
dbe_->setCurrentFolder("Pixel/Tracks");
meNofTracks_ = dbe_->book1D("ntracks_" + tracksrc_.label(),"Number of Tracks",4,0,4);
meNofTracks_->setAxisTitle("Number of Tracks",1);
meNofTracks_->setBinLabel(1,"All");
meNofTracks_->setBinLabel(2,"Pixel");
meNofTracks_->setBinLabel(3,"BPix");
meNofTracks_->setBinLabel(4,"FPix");
//number of tracks in pixel fiducial volume
dbe_->setCurrentFolder("Pixel/Tracks");
meNofTracksInPixVol_ = dbe_->book1D("ntracksInPixVol_" + tracksrc_.label(),"Number of Tracks crossing Pixel fiducial Volume",2,0,2);
meNofTracksInPixVol_->setAxisTitle("Number of Tracks",1);
meNofTracksInPixVol_->setBinLabel(1,"With Hits");
meNofTracksInPixVol_->setBinLabel(2,"Without Hits");
//number of clusters (associated to track / not associated)
dbe_->setCurrentFolder("Pixel/Clusters/OnTrack");
meNofClustersOnTrack_ = dbe_->book1D("nclusters_" + clustersrc_.label() + "_tot","Number of Clusters (on track)",3,0,3);
meNofClustersOnTrack_->setAxisTitle("Number of Clusters on Track",1);
meNofClustersOnTrack_->setBinLabel(1,"All");
meNofClustersOnTrack_->setBinLabel(2,"BPix");
meNofClustersOnTrack_->setBinLabel(3,"FPix");
dbe_->setCurrentFolder("Pixel/Clusters/OffTrack");
meNofClustersNotOnTrack_ = dbe_->book1D("nclusters_" + clustersrc_.label() + "_tot","Number of Clusters (off track)",3,0,3);
meNofClustersNotOnTrack_->setAxisTitle("Number of Clusters off Track",1);
meNofClustersNotOnTrack_->setBinLabel(1,"All");
meNofClustersNotOnTrack_->setBinLabel(2,"BPix");
meNofClustersNotOnTrack_->setBinLabel(3,"FPix");
//cluster charge and size
//charge
//on track
dbe_->setCurrentFolder("Pixel/Clusters/OnTrack");
meClChargeOnTrack_all = dbe_->book1D("charge_" + clustersrc_.label(),"Charge (on track)",500,0.,500.);
meClChargeOnTrack_all->setAxisTitle("Charge size (in ke)",1);
meClChargeOnTrack_bpix = dbe_->book1D("charge_" + clustersrc_.label() + "_Barrel","Charge (on track, barrel)",500,0.,500.);
meClChargeOnTrack_bpix->setAxisTitle("Charge size (in ke)",1);
meClChargeOnTrack_fpix = dbe_->book1D("charge_" + clustersrc_.label() + "_Endcap","Charge (on track, endcap)",500,0.,500.);
meClChargeOnTrack_fpix->setAxisTitle("Charge size (in ke)",1);
meClChargeOnTrack_layer1 = dbe_->book1D("charge_" + clustersrc_.label() + "_Layer_1","Charge (on track, layer1)",500,0.,500.);
meClChargeOnTrack_layer1->setAxisTitle("Charge size (in ke)",1);
meClChargeOnTrack_layer2 = dbe_->book1D("charge_" + clustersrc_.label() + "_Layer_2","Charge (on track, layer2)",500,0.,500.);
meClChargeOnTrack_layer2->setAxisTitle("Charge size (in ke)",1);
meClChargeOnTrack_layer3 = dbe_->book1D("charge_" + clustersrc_.label() + "_Layer_3","Charge (on track, layer3)",500,0.,500.);
meClChargeOnTrack_layer3->setAxisTitle("Charge size (in ke)",1);
meClChargeOnTrack_diskp1 = dbe_->book1D("charge_" + clustersrc_.label() + "_Disk_p1","Charge (on track, diskp1)",500,0.,500.);
meClChargeOnTrack_diskp1->setAxisTitle("Charge size (in ke)",1);
meClChargeOnTrack_diskp2 = dbe_->book1D("charge_" + clustersrc_.label() + "_Disk_p2","Charge (on track, diskp2)",500,0.,500.);
meClChargeOnTrack_diskp2->setAxisTitle("Charge size (in ke)",1);
meClChargeOnTrack_diskm1 = dbe_->book1D("charge_" + clustersrc_.label() + "_Disk_m1","Charge (on track, diskm1)",500,0.,500.);
meClChargeOnTrack_diskm1->setAxisTitle("Charge size (in ke)",1);
meClChargeOnTrack_diskm2 = dbe_->book1D("charge_" + clustersrc_.label() + "_Disk_m2","Charge (on track, diskm2)",500,0.,500.);
meClChargeOnTrack_diskm2->setAxisTitle("Charge size (in ke)",1);
//off track
dbe_->setCurrentFolder("Pixel/Clusters/OffTrack");
meClChargeNotOnTrack_all = dbe_->book1D("charge_" + clustersrc_.label(),"Charge (off track)",500,0.,500.);
meClChargeNotOnTrack_all->setAxisTitle("Charge size (in ke)",1);
meClChargeNotOnTrack_bpix = dbe_->book1D("charge_" + clustersrc_.label() + "_Barrel","Charge (off track, barrel)",500,0.,500.);
meClChargeNotOnTrack_bpix->setAxisTitle("Charge size (in ke)",1);
meClChargeNotOnTrack_fpix = dbe_->book1D("charge_" + clustersrc_.label() + "_Endcap","Charge (off track, endcap)",500,0.,500.);
meClChargeNotOnTrack_fpix->setAxisTitle("Charge size (in ke)",1);
meClChargeNotOnTrack_layer1 = dbe_->book1D("charge_" + clustersrc_.label() + "_Layer_1","Charge (off track, layer1)",500,0.,500.);
meClChargeNotOnTrack_layer1->setAxisTitle("Charge size (in ke)",1);
meClChargeNotOnTrack_layer2 = dbe_->book1D("charge_" + clustersrc_.label() + "_Layer_2","Charge (off track, layer2)",500,0.,500.);
meClChargeNotOnTrack_layer2->setAxisTitle("Charge size (in ke)",1);
meClChargeNotOnTrack_layer3 = dbe_->book1D("charge_" + clustersrc_.label() + "_Layer_3","Charge (off track, layer3)",500,0.,500.);
meClChargeNotOnTrack_layer3->setAxisTitle("Charge size (in ke)",1);
meClChargeNotOnTrack_diskp1 = dbe_->book1D("charge_" + clustersrc_.label() + "_Disk_p1","Charge (off track, diskp1)",500,0.,500.);
meClChargeNotOnTrack_diskp1->setAxisTitle("Charge size (in ke)",1);
meClChargeNotOnTrack_diskp2 = dbe_->book1D("charge_" + clustersrc_.label() + "_Disk_p2","Charge (off track, diskp2)",500,0.,500.);
meClChargeNotOnTrack_diskp2->setAxisTitle("Charge size (in ke)",1);
meClChargeNotOnTrack_diskm1 = dbe_->book1D("charge_" + clustersrc_.label() + "_Disk_m1","Charge (off track, diskm1)",500,0.,500.);
meClChargeNotOnTrack_diskm1->setAxisTitle("Charge size (in ke)",1);
meClChargeNotOnTrack_diskm2 = dbe_->book1D("charge_" + clustersrc_.label() + "_Disk_m2","Charge (off track, diskm2)",500,0.,500.);
meClChargeNotOnTrack_diskm2->setAxisTitle("Charge size (in ke)",1);
//size
//on track
dbe_->setCurrentFolder("Pixel/Clusters/OnTrack");
meClSizeOnTrack_all = dbe_->book1D("size_" + clustersrc_.label(),"Size (on track)",100,0.,100.);
meClSizeOnTrack_all->setAxisTitle("Cluster size (in pixels)",1);
meClSizeOnTrack_bpix = dbe_->book1D("size_" + clustersrc_.label() + "_Barrel","Size (on track, barrel)",100,0.,100.);
meClSizeOnTrack_bpix->setAxisTitle("Cluster size (in pixels)",1);
meClSizeOnTrack_fpix = dbe_->book1D("size_" + clustersrc_.label() + "_Endcap","Size (on track, endcap)",100,0.,100.);
meClSizeOnTrack_fpix->setAxisTitle("Cluster size (in pixels)",1);
meClSizeOnTrack_layer1 = dbe_->book1D("size_" + clustersrc_.label() + "_Layer_1","Size (on track, layer1)",100,0.,100.);
meClSizeOnTrack_layer1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeOnTrack_layer2 = dbe_->book1D("size_" + clustersrc_.label() + "_Layer_2","Size (on track, layer2)",100,0.,100.);
meClSizeOnTrack_layer2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeOnTrack_layer3 = dbe_->book1D("size_" + clustersrc_.label() + "_Layer_3","Size (on track, layer3)",100,0.,100.);
meClSizeOnTrack_layer3->setAxisTitle("Cluster size (in pixels)",1);
meClSizeOnTrack_diskp1 = dbe_->book1D("size_" + clustersrc_.label() + "_Disk_p1","Size (on track, diskp1)",100,0.,100.);
meClSizeOnTrack_diskp1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeOnTrack_diskp2 = dbe_->book1D("size_" + clustersrc_.label() + "_Disk_p2","Size (on track, diskp2)",100,0.,100.);
meClSizeOnTrack_diskp2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeOnTrack_diskm1 = dbe_->book1D("size_" + clustersrc_.label() + "_Disk_m1","Size (on track, diskm1)",100,0.,100.);
meClSizeOnTrack_diskm1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeOnTrack_diskm2 = dbe_->book1D("size_" + clustersrc_.label() + "_Disk_m2","Size (on track, diskm2)",100,0.,100.);
meClSizeOnTrack_diskm2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXOnTrack_all = dbe_->book1D("sizeX_" + clustersrc_.label(),"SizeX (on track)",100,0.,100.);
meClSizeXOnTrack_all->setAxisTitle("Cluster sizeX (in pixels)",1);
meClSizeXOnTrack_bpix = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Barrel","SizeX (on track, barrel)",100,0.,100.);
meClSizeXOnTrack_bpix->setAxisTitle("Cluster sizeX (in pixels)",1);
meClSizeXOnTrack_fpix = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Endcap","SizeX (on track, endcap)",100,0.,100.);
meClSizeXOnTrack_fpix->setAxisTitle("Cluster sizeX (in pixels)",1);
meClSizeXOnTrack_layer1 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Layer_1","SizeX (on track, layer1)",100,0.,100.);
meClSizeXOnTrack_layer1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXOnTrack_layer2 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Layer_2","SizeX (on track, layer2)",100,0.,100.);
meClSizeXOnTrack_layer2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXOnTrack_layer3 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Layer_3","SizeX (on track, layer3)",100,0.,100.);
meClSizeXOnTrack_layer3->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXOnTrack_diskp1 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Disk_p1","SizeX (on track, diskp1)",100,0.,100.);
meClSizeXOnTrack_diskp1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXOnTrack_diskp2 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Disk_p2","SizeX (on track, diskp2)",100,0.,100.);
meClSizeXOnTrack_diskp2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXOnTrack_diskm1 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Disk_m1","SizeX (on track, diskm1)",100,0.,100.);
meClSizeXOnTrack_diskm1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXOnTrack_diskm2 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Disk_m2","SizeX (on track, diskm2)",100,0.,100.);
meClSizeXOnTrack_diskm2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYOnTrack_all = dbe_->book1D("sizeY_" + clustersrc_.label(),"SizeY (on track)",100,0.,100.);
meClSizeYOnTrack_all->setAxisTitle("Cluster sizeY (in pixels)",1);
meClSizeYOnTrack_bpix = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Barrel","SizeY (on track, barrel)",100,0.,100.);
meClSizeYOnTrack_bpix->setAxisTitle("Cluster sizeY (in pixels)",1);
meClSizeYOnTrack_fpix = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Endcap","SizeY (on track, endcap)",100,0.,100.);
meClSizeYOnTrack_fpix->setAxisTitle("Cluster sizeY (in pixels)",1);
meClSizeYOnTrack_layer1 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Layer_1","SizeY (on track, layer1)",100,0.,100.);
meClSizeYOnTrack_layer1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYOnTrack_layer2 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Layer_2","SizeY (on track, layer2)",100,0.,100.);
meClSizeYOnTrack_layer2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYOnTrack_layer3 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Layer_3","SizeY (on track, layer3)",100,0.,100.);
meClSizeYOnTrack_layer3->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYOnTrack_diskp1 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Disk_p1","SizeY (on track, diskp1)",100,0.,100.);
meClSizeYOnTrack_diskp1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYOnTrack_diskp2 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Disk_p2","SizeY (on track, diskp2)",100,0.,100.);
meClSizeYOnTrack_diskp2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYOnTrack_diskm1 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Disk_m1","SizeY (on track, diskm1)",100,0.,100.);
meClSizeYOnTrack_diskm1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYOnTrack_diskm2 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Disk_m2","SizeY (on track, diskm2)",100,0.,100.);
meClSizeYOnTrack_diskm2->setAxisTitle("Cluster size (in pixels)",1);
//off track
dbe_->setCurrentFolder("Pixel/Clusters/OffTrack");
meClSizeNotOnTrack_all = dbe_->book1D("size_" + clustersrc_.label(),"Size (off track)",100,0.,100.);
meClSizeNotOnTrack_all->setAxisTitle("Cluster size (in pixels)",1);
meClSizeNotOnTrack_bpix = dbe_->book1D("size_" + clustersrc_.label() + "_Barrel","Size (off track, barrel)",100,0.,100.);
meClSizeNotOnTrack_bpix->setAxisTitle("Cluster size (in pixels)",1);
meClSizeNotOnTrack_fpix = dbe_->book1D("size_" + clustersrc_.label() + "_Endcap","Size (off track, endcap)",100,0.,100.);
meClSizeNotOnTrack_fpix->setAxisTitle("Cluster size (in pixels)",1);
meClSizeNotOnTrack_layer1 = dbe_->book1D("size_" + clustersrc_.label() + "_Layer_1","Size (off track, layer1)",100,0.,100.);
meClSizeNotOnTrack_layer1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeNotOnTrack_layer2 = dbe_->book1D("size_" + clustersrc_.label() + "_Layer_2","Size (off track, layer2)",100,0.,100.);
meClSizeNotOnTrack_layer2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeNotOnTrack_layer3 = dbe_->book1D("size_" + clustersrc_.label() + "_Layer_3","Size (off track, layer3)",100,0.,100.);
meClSizeNotOnTrack_layer3->setAxisTitle("Cluster size (in pixels)",1);
meClSizeNotOnTrack_diskp1 = dbe_->book1D("size_" + clustersrc_.label() + "_Disk_p1","Size (off track, diskp1)",100,0.,100.);
meClSizeNotOnTrack_diskp1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeNotOnTrack_diskp2 = dbe_->book1D("size_" + clustersrc_.label() + "_Disk_p2","Size (off track, diskp2)",100,0.,100.);
meClSizeNotOnTrack_diskp2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeNotOnTrack_diskm1 = dbe_->book1D("size_" + clustersrc_.label() + "_Disk_m1","Size (off track, diskm1)",100,0.,100.);
meClSizeNotOnTrack_diskm1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeNotOnTrack_diskm2 = dbe_->book1D("size_" + clustersrc_.label() + "_Disk_m2","Size (off track, diskm2)",100,0.,100.);
meClSizeNotOnTrack_diskm2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXNotOnTrack_all = dbe_->book1D("sizeX_" + clustersrc_.label(),"SizeX (off track)",100,0.,100.);
meClSizeXNotOnTrack_all->setAxisTitle("Cluster sizeX (in pixels)",1);
meClSizeXNotOnTrack_bpix = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Barrel","SizeX (off track, barrel)",100,0.,100.);
meClSizeXNotOnTrack_bpix->setAxisTitle("Cluster sizeX (in pixels)",1);
meClSizeXNotOnTrack_fpix = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Endcap","SizeX (off track, endcap)",100,0.,100.);
meClSizeXNotOnTrack_fpix->setAxisTitle("Cluster sizeX (in pixels)",1);
meClSizeXNotOnTrack_layer1 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Layer_1","SizeX (off track, layer1)",100,0.,100.);
meClSizeXNotOnTrack_layer1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXNotOnTrack_layer2 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Layer_2","SizeX (off track, layer2)",100,0.,100.);
meClSizeXNotOnTrack_layer2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXNotOnTrack_layer3 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Layer_3","SizeX (off track, layer3)",100,0.,100.);
meClSizeXNotOnTrack_layer3->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXNotOnTrack_diskp1 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Disk_p1","SizeX (off track, diskp1)",100,0.,100.);
meClSizeXNotOnTrack_diskp1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXNotOnTrack_diskp2 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Disk_p2","SizeX (off track, diskp2)",100,0.,100.);
meClSizeXNotOnTrack_diskp2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXNotOnTrack_diskm1 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Disk_m1","SizeX (off track, diskm1)",100,0.,100.);
meClSizeXNotOnTrack_diskm1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeXNotOnTrack_diskm2 = dbe_->book1D("sizeX_" + clustersrc_.label() + "_Disk_m2","SizeX (off track, diskm2)",100,0.,100.);
meClSizeXNotOnTrack_diskm2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYNotOnTrack_all = dbe_->book1D("sizeY_" + clustersrc_.label(),"SizeY (off track)",100,0.,100.);
meClSizeYNotOnTrack_all->setAxisTitle("Cluster sizeY (in pixels)",1);
meClSizeYNotOnTrack_bpix = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Barrel","SizeY (off track, barrel)",100,0.,100.);
meClSizeYNotOnTrack_bpix->setAxisTitle("Cluster sizeY (in pixels)",1);
meClSizeYNotOnTrack_fpix = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Endcap","SizeY (off track, endcap)",100,0.,100.);
meClSizeYNotOnTrack_fpix->setAxisTitle("Cluster sizeY (in pixels)",1);
meClSizeYNotOnTrack_layer1 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Layer_1","SizeY (off track, layer1)",100,0.,100.);
meClSizeYNotOnTrack_layer1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYNotOnTrack_layer2 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Layer_2","SizeY (off track, layer2)",100,0.,100.);
meClSizeYNotOnTrack_layer2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYNotOnTrack_layer3 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Layer_3","SizeY (off track, layer3)",100,0.,100.);
meClSizeYNotOnTrack_layer3->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYNotOnTrack_diskp1 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Disk_p1","SizeY (off track, diskp1)",100,0.,100.);
meClSizeYNotOnTrack_diskp1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYNotOnTrack_diskp2 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Disk_p2","SizeY (off track, diskp2)",100,0.,100.);
meClSizeYNotOnTrack_diskp2->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYNotOnTrack_diskm1 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Disk_m1","SizeY (off track, diskm1)",100,0.,100.);
meClSizeYNotOnTrack_diskm1->setAxisTitle("Cluster size (in pixels)",1);
meClSizeYNotOnTrack_diskm2 = dbe_->book1D("sizeY_" + clustersrc_.label() + "_Disk_m2","SizeY (off track, diskm2)",100,0.,100.);
meClSizeYNotOnTrack_diskm2->setAxisTitle("Cluster size (in pixels)",1);
//cluster global position
//on track
dbe_->setCurrentFolder("Pixel/Clusters/OnTrack");
//bpix
meClPosLayer1OnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_Layer_1","Clusters Layer1 (on track)",200,-30.,30.,128,-3.2,3.2);
meClPosLayer1OnTrack->setAxisTitle("Global Z (cm)",1);
meClPosLayer1OnTrack->setAxisTitle("Global #phi",2);
meClPosLayer2OnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_Layer_2","Clusters Layer2 (on track)",200,-30.,30.,128,-3.2,3.2);
meClPosLayer2OnTrack->setAxisTitle("Global Z (cm)",1);
meClPosLayer2OnTrack->setAxisTitle("Global #phi",2);
meClPosLayer3OnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_Layer_3","Clusters Layer3 (on track)",200,-30.,30.,128,-3.2,3.2);
meClPosLayer3OnTrack->setAxisTitle("Global Z (cm)",1);
meClPosLayer3OnTrack->setAxisTitle("Global #phi",2);
//fpix
meClPosDisk1pzOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_pz_Disk_1","Clusters +Z Disk1 (on track)",80,-20.,20.,80,-20.,20.);
meClPosDisk1pzOnTrack->setAxisTitle("Global X (cm)",1);
meClPosDisk1pzOnTrack->setAxisTitle("Global Y (cm)",2);
meClPosDisk2pzOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_pz_Disk_2","Clusters +Z Disk2 (on track)",80,-20.,20.,80,-20.,20.);
meClPosDisk2pzOnTrack->setAxisTitle("Global X (cm)",1);
meClPosDisk2pzOnTrack->setAxisTitle("Global Y (cm)",2);
meClPosDisk1mzOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_mz_Disk_1","Clusters -Z Disk1 (on track)",80,-20.,20.,80,-20.,20.);
meClPosDisk1mzOnTrack->setAxisTitle("Global X (cm)",1);
meClPosDisk1mzOnTrack->setAxisTitle("Global Y (cm)",2);
meClPosDisk2mzOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_mz_Disk_2","Clusters -Z Disk2 (on track)",80,-20.,20.,80,-20.,20.);
meClPosDisk2mzOnTrack->setAxisTitle("Global X (cm)",1);
meClPosDisk2mzOnTrack->setAxisTitle("Global Y (cm)",2);
meNClustersOnTrack_all = dbe_->book1D("nclusters_" + clustersrc_.label(),"Number of Clusters (on Track)",50,0.,50.);
meNClustersOnTrack_all->setAxisTitle("Number of Clusters",1);
meNClustersOnTrack_bpix = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Barrel","Number of Clusters (on track, barrel)",50,0.,50.);
meNClustersOnTrack_bpix->setAxisTitle("Number of Clusters",1);
meNClustersOnTrack_fpix = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Endcap","Number of Clusters (on track, endcap)",50,0.,50.);
meNClustersOnTrack_fpix->setAxisTitle("Number of Clusters",1);
meNClustersOnTrack_layer1 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Layer_1","Number of Clusters (on track, layer1)",50,0.,50.);
meNClustersOnTrack_layer1->setAxisTitle("Number of Clusters",1);
meNClustersOnTrack_layer2 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Layer_2","Number of Clusters (on track, layer2)",50,0.,50.);
meNClustersOnTrack_layer2->setAxisTitle("Number of Clusters",1);
meNClustersOnTrack_layer3 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Layer_3","Number of Clusters (on track, layer3)",50,0.,50.);
meNClustersOnTrack_layer3->setAxisTitle("Number of Clusters",1);
meNClustersOnTrack_diskp1 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Disk_p1","Number of Clusters (on track, diskp1)",50,0.,50.);
meNClustersOnTrack_diskp1->setAxisTitle("Number of Clusters",1);
meNClustersOnTrack_diskp2 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Disk_p2","Number of Clusters (on track, diskp2)",50,0.,50.);
meNClustersOnTrack_diskp2->setAxisTitle("Number of Clusters",1);
meNClustersOnTrack_diskm1 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Disk_m1","Number of Clusters (on track, diskm1)",50,0.,50.);
meNClustersOnTrack_diskm1->setAxisTitle("Number of Clusters",1);
meNClustersOnTrack_diskm2 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Disk_m2","Number of Clusters (on track, diskm2)",50,0.,50.);
meNClustersOnTrack_diskm2->setAxisTitle("Number of Clusters",1);
//not on track
dbe_->setCurrentFolder("Pixel/Clusters/OffTrack");
//bpix
meClPosLayer1NotOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_Layer_1","Clusters Layer1 (off track)",200,-30.,30.,128,-3.2,3.2);
meClPosLayer1NotOnTrack->setAxisTitle("Global Z (cm)",1);
meClPosLayer1NotOnTrack->setAxisTitle("Global #phi",2);
meClPosLayer2NotOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_Layer_2","Clusters Layer2 (off track)",200,-30.,30.,128,-3.2,3.2);
meClPosLayer2NotOnTrack->setAxisTitle("Global Z (cm)",1);
meClPosLayer2NotOnTrack->setAxisTitle("Global #phi",2);
meClPosLayer3NotOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_Layer_3","Clusters Layer3 (off track)",200,-30.,30.,128,-3.2,3.2);
meClPosLayer3NotOnTrack->setAxisTitle("Global Z (cm)",1);
meClPosLayer3NotOnTrack->setAxisTitle("Global #phi",2);
//fpix
meClPosDisk1pzNotOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_pz_Disk_1","Clusters +Z Disk1 (off track)",80,-20.,20.,80,-20.,20.);
meClPosDisk1pzNotOnTrack->setAxisTitle("Global X (cm)",1);
meClPosDisk1pzNotOnTrack->setAxisTitle("Global Y (cm)",2);
meClPosDisk2pzNotOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_pz_Disk_2","Clusters +Z Disk2 (off track)",80,-20.,20.,80,-20.,20.);
meClPosDisk2pzNotOnTrack->setAxisTitle("Global X (cm)",1);
meClPosDisk2pzNotOnTrack->setAxisTitle("Global Y (cm)",2);
meClPosDisk1mzNotOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_mz_Disk_1","Clusters -Z Disk1 (off track)",80,-20.,20.,80,-20.,20.);
meClPosDisk1mzNotOnTrack->setAxisTitle("Global X (cm)",1);
meClPosDisk1mzNotOnTrack->setAxisTitle("Global Y (cm)",2);
meClPosDisk2mzNotOnTrack = dbe_->book2D("position_" + clustersrc_.label() + "_mz_Disk_2","Clusters -Z Disk2 (off track)",80,-20.,20.,80,-20.,20.);
meClPosDisk2mzNotOnTrack->setAxisTitle("Global X (cm)",1);
meClPosDisk2mzNotOnTrack->setAxisTitle("Global Y (cm)",2);
meNClustersNotOnTrack_all = dbe_->book1D("nclusters_" + clustersrc_.label(),"Number of Clusters (off Track)",50,0.,50.);
meNClustersNotOnTrack_all->setAxisTitle("Number of Clusters",1);
meNClustersNotOnTrack_bpix = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Barrel","Number of Clusters (off track, barrel)",50,0.,50.);
meNClustersNotOnTrack_bpix->setAxisTitle("Number of Clusters",1);
meNClustersNotOnTrack_fpix = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Endcap","Number of Clusters (off track, endcap)",50,0.,50.);
meNClustersNotOnTrack_fpix->setAxisTitle("Number of Clusters",1);
meNClustersNotOnTrack_layer1 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Layer_1","Number of Clusters (off track, layer1)",50,0.,50.);
meNClustersNotOnTrack_layer1->setAxisTitle("Number of Clusters",1);
meNClustersNotOnTrack_layer2 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Layer_2","Number of Clusters (off track, layer2)",50,0.,50.);
meNClustersNotOnTrack_layer2->setAxisTitle("Number of Clusters",1);
meNClustersNotOnTrack_layer3 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Layer_3","Number of Clusters (off track, layer3)",50,0.,50.);
meNClustersNotOnTrack_layer3->setAxisTitle("Number of Clusters",1);
meNClustersNotOnTrack_diskp1 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Disk_p1","Number of Clusters (off track, diskp1)",50,0.,50.);
meNClustersNotOnTrack_diskp1->setAxisTitle("Number of Clusters",1);
meNClustersNotOnTrack_diskp2 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Disk_p2","Number of Clusters (off track, diskp2)",50,0.,50.);
meNClustersNotOnTrack_diskp2->setAxisTitle("Number of Clusters",1);
meNClustersNotOnTrack_diskm1 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Disk_m1","Number of Clusters (off track, diskm1)",50,0.,50.);
meNClustersNotOnTrack_diskm1->setAxisTitle("Number of Clusters",1);
meNClustersNotOnTrack_diskm2 = dbe_->book1D("nclusters_" + clustersrc_.label() + "_Disk_m2","Number of Clusters (off track, diskm2)",50,0.,50.);
meNClustersNotOnTrack_diskm2->setAxisTitle("Number of Clusters",1);
//HitProbability
//on track
dbe_->setCurrentFolder("Pixel/Clusters/OnTrack");
meHitProbability = dbe_->book1D("FractionLowProb","Fraction of hits with low probability;FractionLowProb;#HitsOnTrack",100,0.,1.);
if (debug_) {
// book summary residual histograms in a debugging folder - one (x,y) pair of histograms per subdetector
dbe_->setCurrentFolder("debugging");
char hisID[80];
for (int s=0; s<3; s++) {
sprintf(hisID,"residual_x_subdet_%i",s);
meSubdetResidualX[s] = dbe_->book1D(hisID,"Pixel Hit-to-Track Residual in X",500,-5.,5.);
sprintf(hisID,"residual_y_subdet_%i",s);
meSubdetResidualY[s] = dbe_->book1D(hisID,"Pixel Hit-to-Track Residual in Y",500,-5.,5.);
}
}
firstRun = false;
}
}
| void SiPixelTrackResidualSource::endJob | ( | void | ) | [virtual] |
Reimplemented from edm::EDAnalyzer.
Definition at line 501 of file SiPixelTrackResidualSource.cc.
References dbe_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), download_sqlite_cfg::outputFile, pSet_, DQMStore::save(), and AlCaHLTBitMon_QueryRunRegistry::string.
{
LogInfo("PixelDQM") << "SiPixelTrackResidualSource endJob()";
// save the residual histograms to an output root file
bool saveFile = pSet_.getUntrackedParameter<bool>("saveFile", true);
if (saveFile) {
std::string outputFile = pSet_.getParameter<std::string>("outputFile");
LogInfo("PixelDQM") << " - saving histograms to "<< outputFile.data();
dbe_->save(outputFile);
}
LogInfo("PixelDQM") << endl; // dbe_->showDirStructure();
}
| void SiPixelTrackResidualSource::triplets | ( | double | x1, |
| double | y1, | ||
| double | z1, | ||
| double | x2, | ||
| double | y2, | ||
| double | z2, | ||
| double | x3, | ||
| double | y3, | ||
| double | z3, | ||
| double | ptsig, | ||
| double & | dc, | ||
| double & | dz, | ||
| double | kap | ||
| ) |
Definition at line 1263 of file SiPixelTrackResidualSource.cc.
References gather_cfg::cout, kap, dttmaxenums::L, lam, diffTwoXMLs::r1, rho, and mathSSE::sqrt().
Referenced by analyze().
{
//Define some constants
using namespace std;
//Curvature kap from global Track
//inverse of the curvature is the radius in the transverse plane
double rho = 1/kap;
//Check that the hits are in the correct layers
double r1 = sqrt( x1*x1 + y1*y1 );
double r3 = sqrt( x3*x3 + y3*y3 );
if( r3-r1 < 2.0 ) cout << "warn r1 = " << r1 << ", r3 = " << r3 << endl;
// Calculate the centre of the helix in xy-projection with radius rho from the track.
//start with a line (sekante) connecting the two points (x1,y1) and (x3,y3) vec_L = vec_x3-vec_x1
//with L being the length of that vector.
double L=sqrt((x3-x1)*(x3-x1)+(y3-y1)*(y3-y1));
//lam is the line from the middel point of vec_q towards the center of the circle X0,Y0
// we already have kap and rho = 1/kap
double lam = sqrt(rho*rho - L*L/4);
// There are two solutions, the sign of kap gives the information
// which of them is correct.
//
if( kap > 0 ) lam = -lam;
//
// ( X0, Y0 ) is the centre of the circle that describes the helix in xy-projection.
//
double x0 = 0.5*( x1 + x3 ) + lam/L * ( -y1 + y3 );
double y0 = 0.5*( y1 + y3 ) + lam/L * ( x1 - x3 );
// Calculate the dipangle in z direction (needed later for z residual) :
//Starting from the heliz equation whihc has to hold for both points z1,z3
double num = ( y3 - y0 ) * ( x1 - x0 ) - ( x3 - x0 ) * ( y1 - y0 );
double den = ( x1 - x0 ) * ( x3 - x0 ) + ( y1 - y0 ) * ( y3 - y0 );
double tandip = kap * ( z3 - z1 ) / atan( num / den );
// angle from first hit to dca point:
//
double dphi = atan( ( ( x1 - x0 ) * y0 - ( y1 - y0 ) * x0 )
/ ( ( x1 - x0 ) * x0 + ( y1 - y0 ) * y0 ) );
//z position of the track based on the middle of the circle
//track equation for the z component
double uz0 = z1 + tandip * dphi * rho;
//RESIDUAL IN R-PHI
//Calculate distance dca2 from point (x2,y2) to the circle which is given by
//the distance of the point to the middlepoint dcM = sqrt((x0-x2)^2+(y0-y2)) and rho
//dca = rho +- dcM
if(kap>0) dca2=rho-sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2));
else dca2=rho+sqrt((-x0+x2)*(-x0+x2)+(-y0+y2)*(-y0+y2));
//RESIDUAL IN Z
double xx =0 ;
double yy =0 ;
//sign of kappa determines the calculation
//xx and yy are the new coordinates starting from x2, y2 that are on the track itself
//vec_X2+-dca2*vec(X0-X2)/|(X0-X2)|
if(kap<0){
xx = x2+(dca2*((x0-x2))/sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2)));
yy = y2+(dca2*((y0-y2))/sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2)));
}
else if(kap>=0){
xx = x2-(dca2*((x0-x2))/sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2)));
yy = y2-(dca2*((y0-y2))/sqrt((x0-x2)*(x0-x2)+(y0-y2)*(y0-y2)));
}
//to get residual in z start with calculating the new uz2 position if one has moved to xx, yy
//on the track. First calculate the change in phi2 with respect to the center X0, Y0
double dphi2 = atan( ( ( xx - x0 ) * y0 - ( yy - y0 ) * x0 )
/ ( ( xx - x0 ) * x0 + ( yy - y0 ) * y0 ) );
//Solve track equation for this new z depending on the dip angle of the track (see above
//calculated based on X1, X3 and X0, use uz0 as reference point again.
double uz2= uz0 - dphi2*tandip*rho;
//subtract new z position from the old one
dz2=z2-uz2;
//if we are interested in the arclength this is unsigned though
// double cosphi2 = (x2*xx+y2*yy)/(sqrt(x2*x2+y2*y2)*sqrt(xx*xx+yy*yy));
//double arcdca2=sqrt(x2*x2+y2*y2)*acos(cosphi2);
}
bool SiPixelTrackResidualSource::bladeOn [private] |
Definition at line 66 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
Definition at line 55 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
DQMStore* SiPixelTrackResidualSource::dbe_ [private] |
Definition at line 58 of file SiPixelTrackResidualSource.h.
Referenced by beginRun(), endJob(), and SiPixelTrackResidualSource().
bool SiPixelTrackResidualSource::debug_ [private] |
Definition at line 60 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
bool SiPixelTrackResidualSource::diskOn [private] |
Definition at line 66 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
bool SiPixelTrackResidualSource::firstRun [private] |
Definition at line 68 of file SiPixelTrackResidualSource.h.
Referenced by beginJob(), and beginRun().
bool SiPixelTrackResidualSource::ladOn [private] |
Definition at line 64 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
bool SiPixelTrackResidualSource::layOn [private] |
Definition at line 64 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
Definition at line 91 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 92 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 99 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 100 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 97 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 98 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 93 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 94 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 95 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 96 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 81 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 82 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 89 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 90 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 87 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 88 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 83 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 84 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 85 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 86 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 198 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 194 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 196 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 192 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 199 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 195 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 197 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 193 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 188 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 185 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 189 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 186 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 190 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 187 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 111 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 112 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 119 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 120 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 117 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 118 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 113 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 114 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 115 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 116 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 101 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 102 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 109 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 110 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 107 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 108 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 103 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 104 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 105 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 106 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 131 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 132 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 139 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 140 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 137 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 138 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 133 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 134 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 135 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 136 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 121 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 122 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 129 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 130 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 127 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 128 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 123 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 124 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 125 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 126 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 151 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 152 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 159 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 160 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 157 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 158 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 153 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 154 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 155 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 156 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 141 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 142 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 149 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 150 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 147 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 148 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 143 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 144 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 145 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 146 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 201 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 173 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 174 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 181 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 182 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 179 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 180 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 175 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 176 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 177 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 178 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 163 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 164 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 171 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 172 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 169 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 170 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 165 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 166 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 167 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 168 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 80 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 79 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 77 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
Definition at line 78 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
MonitorElement* SiPixelTrackResidualSource::meSubdetResidualX[3] [private] |
Definition at line 74 of file SiPixelTrackResidualSource.h.
Referenced by beginRun().
MonitorElement* SiPixelTrackResidualSource::meSubdetResidualY[3] [private] |
Definition at line 75 of file SiPixelTrackResidualSource.h.
Referenced by beginRun().
bool SiPixelTrackResidualSource::modOn [private] |
Definition at line 61 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
int SiPixelTrackResidualSource::NLowProb [private] |
Definition at line 70 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginJob().
int SiPixelTrackResidualSource::NTotal [private] |
Definition at line 69 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginJob().
bool SiPixelTrackResidualSource::phiOn [private] |
Definition at line 64 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
Definition at line 53 of file SiPixelTrackResidualSource.h.
Referenced by beginRun(), endJob(), and SiPixelTrackResidualSource().
double SiPixelTrackResidualSource::ptminres_ [private] |
Definition at line 67 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and SiPixelTrackResidualSource().
bool SiPixelTrackResidualSource::reducedSet [private] |
Definition at line 62 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and beginRun().
bool SiPixelTrackResidualSource::ringOn [private] |
Definition at line 66 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
Definition at line 54 of file SiPixelTrackResidualSource.h.
Referenced by SiPixelTrackResidualSource().
std::map<uint32_t, SiPixelTrackResidualModule*> SiPixelTrackResidualSource::theSiPixelStructure [private] |
Definition at line 72 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and ~SiPixelTrackResidualSource().
Definition at line 56 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), beginRun(), and SiPixelTrackResidualSource().
std::string SiPixelTrackResidualSource::ttrhbuilder_ [private] |
Definition at line 57 of file SiPixelTrackResidualSource.h.
Referenced by analyze(), and SiPixelTrackResidualSource().
1.7.3