DTSegment2DSLPhiQuality Class Reference

#include <DTSegment2DSLPhiQuality.h>

Inheritance diagram for DTSegment2DSLPhiQuality:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &eventSetup)
	Perform the real analysis. More...
	DTSegment2DSLPhiQuality (const edm::ParameterSet &pset)
	Constructor. More...
void	endJob ()
virtual	~DTSegment2DSLPhiQuality ()
	Destructor. More...
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Private Attributes
bool	debug
HEff2DHit *	h2DHitEff_SuperPhi
HRes2DHit *	h2DHitSuperPhi
std::string	rootFileName
std::string	segment4DLabel
double	sigmaResAngle
double	sigmaResPos
std::string	simHitLabel
TFile *	theFile

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Detailed Description

Basic analyzer class which accesses 2D DTSegments reconstructed with both SL Phi and plot resolution comparing reconstructed and simulated quantities

Date:

2007/06/08 15:50:31

Revision:

1.2

Author

S. Bolognesi and G. Cerminara - INFN Torino

Definition at line 28 of file DTSegment2DSLPhiQuality.h.

Constructor & Destructor Documentation

DTSegment2DSLPhiQuality::DTSegment2DSLPhiQuality

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 38 of file DTSegment2DSLPhiQuality.cc.

References debug, DTHitQualityUtils::debug, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), dtT0WireCalibration_cfg::rootFileName, and interactiveExample::theFile.

                                                                          {
   // Get the debug parameter for verbose output
  debug = pset.getUntrackedParameter<bool>("debug");
  DTHitQualityUtils::debug = debug;
 
  rootFileName = pset.getUntrackedParameter<string>("rootFileName");
  // the name of the simhit collection
  simHitLabel = pset.getUntrackedParameter<string>("simHitLabel", "SimG4Object");
  // the name of the 4D rec hit collection
  segment4DLabel = pset.getUntrackedParameter<string>("segment4DLabel");

  //sigma resolution on position
  sigmaResPos = pset.getParameter<double>("sigmaResPos");
  //sigma resolution on angle
  sigmaResAngle = pset.getParameter<double>("sigmaResAngle");

  // Create the root file
  theFile = new TFile(rootFileName.c_str(), "RECREATE");
  theFile->cd();

  // Book the histos
  h2DHitSuperPhi = new HRes2DHit ("SuperPhi");
  h2DHitEff_SuperPhi = new HEff2DHit ("SuperPhi");
}

DTSegment2DSLPhiQuality::~DTSegment2DSLPhiQuality ( )

virtual

Destructor.

Definition at line 64 of file DTSegment2DSLPhiQuality.cc.

                                                 {
}

Member Function Documentation

void DTSegment2DSLPhiQuality::analyze ( const edm::Event &  event, const edm::EventSetup &  eventSetup  )

virtual

Perform the real analysis.

Implements edm::EDAnalyzer.

Definition at line 80 of file DTSegment2DSLPhiQuality.cc.

References gather_cfg::cout, debug, PV3DBase< T, PVType, FrameType >::eta(), DTHitQualityUtils::findMuSimSegment(), DTHitQualityUtils::findMuSimSegmentDirAndPos(), DTHitQualityUtils::findSegmentAlphaAndBeta(), edm::EventSetup::get(), DTRecSegment2D::localDirection(), DTRecSegment2D::localDirectionError(), DTRecSegment2D::localPosition(), DTRecSegment2D::localPositionError(), DTHitQualityUtils::mapMuSimHitsPerWire(), DTHitQualityUtils::mapSimHitsPerWire(), PV3DBase< T, PVType, FrameType >::phi(), trackerHits::simHits, mathSSE::sqrt(), interactiveExample::theFile, GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), and LocalError::xx().

                                                                                      {
  theFile->cd();

  // Get the DT Geometry
  ESHandle<DTGeometry> dtGeom;
  eventSetup.get<MuonGeometryRecord>().get(dtGeom);

  // Get the SimHit collection from the event
  edm::Handle<PSimHitContainer> simHits;
  event.getByLabel(simHitLabel, "MuonDTHits", simHits); //FIXME: second string to be removed

  //Map simHits by chamber
  map<DTChamberId, PSimHitContainer > simHitsPerCh;
  for(PSimHitContainer::const_iterator simHit = simHits->begin();
      simHit != simHits->end(); simHit++){
    // Create the id of the chamber (the simHits in the DT known their wireId)
    DTChamberId chamberId = (((DTWireId(simHit->detUnitId())).layerId()).superlayerId()).chamberId();
    // Fill the map
    simHitsPerCh[chamberId].push_back(*simHit);
  }

  // Get the 4D rechits from the event
  Handle<DTRecSegment4DCollection> segment4Ds;
  event.getByLabel(segment4DLabel, segment4Ds);
    
  // Loop over all chambers containing a segment
  DTRecSegment4DCollection::id_iterator chamberId;
  for (chamberId = segment4Ds->id_begin();
       chamberId != segment4Ds->id_end();
       ++chamberId){
    
    //------------------------- simHits ---------------------------//
    //Get simHits of each chamber
    PSimHitContainer simHits =  simHitsPerCh[(*chamberId)];
       
    // Map simhits per wire
    map<DTWireId, PSimHitContainer > simHitsPerWire = DTHitQualityUtils::mapSimHitsPerWire(simHits);
    map<DTWireId, const PSimHit*> muSimHitPerWire = DTHitQualityUtils::mapMuSimHitsPerWire(simHitsPerWire);
    int nMuSimHit = muSimHitPerWire.size();
    if(nMuSimHit == 0 || nMuSimHit == 1) {
      if(debug && nMuSimHit == 1)
        cout << "[DTSegment2DSLPhiQuality] Only " << nMuSimHit << " mu SimHit in this chamber, skipping!" << endl;
      continue; // If no or only one mu SimHit is found skip this chamber
    } 
    if(debug)
      cout << "=== Chamber " << (*chamberId) << " has " << nMuSimHit << " SimHits" << endl;
  
    //Find outer and inner mu SimHit to build a segment
    pair<const PSimHit*, const PSimHit*> inAndOutSimHit = DTHitQualityUtils::findMuSimSegment(muSimHitPerWire); 

    //Find direction and position of the sim Segment in Chamber RF
    pair<LocalVector, LocalPoint> dirAndPosSimSegm = DTHitQualityUtils::findMuSimSegmentDirAndPos(inAndOutSimHit,
                                                  (*chamberId),&(*dtGeom));

    LocalVector simSegmLocalDir = dirAndPosSimSegm.first;
    LocalPoint simSegmLocalPos = dirAndPosSimSegm.second;
    const DTChamber* chamber = dtGeom->chamber(*chamberId);
    GlobalPoint simSegmGlobalPos = chamber->toGlobal(simSegmLocalPos);

    //Atan(x/z) angle and x position in SL RF
    float angleSimSeg = DTHitQualityUtils::findSegmentAlphaAndBeta(simSegmLocalDir).first;
    float posSimSeg = simSegmLocalPos.x(); 
    //Position (in eta,phi coordinates) in lobal RF
    float etaSimSeg = simSegmGlobalPos.eta(); 
    float phiSimSeg = simSegmGlobalPos.phi();

    if(debug)
      cout<<"  Simulated segment:  local direction "<<simSegmLocalDir<<endl
      <<"                      local position  "<<simSegmLocalPos<<endl
      <<"                      angle           "<<angleSimSeg<<endl;
    
    //---------------------------- recHits --------------------------//
    // Get the range of rechit for the corresponding chamberId
    bool recHitFound = false;
    DTRecSegment4DCollection::range range = segment4Ds->get(*chamberId);
    int nsegm = distance(range.first, range.second);
    if(debug)
      cout << "   Chamber: " << *chamberId << " has " << nsegm
       << " 4D segments" << endl;

    if (nsegm!=0) {
      // Find the best RecHit: look for the 4D RecHit with the phi angle closest
      // to that of segment made of SimHits. 
      // RecHits must have delta alpha and delta position within 5 sigma of
      // the residual distribution (we are looking for residuals of segments
      // usefull to the track fit) for efficency purpose
      const DTRecSegment2D* bestRecHit = 0;
      bool bestRecHitFound = false;
      double deltaAlpha = 99999;

      // Loop over the recHits of this chamberId
      for (DTRecSegment4DCollection::const_iterator segment4D = range.first;
           segment4D!=range.second;
           ++segment4D){
        // Check the dimension
        if((*segment4D).dimension() != 4) {
          if(debug) cout << "[DTSegment2DSLPhiQuality]***Error: This is not 4D segment!!!" << endl;
          continue;
        }

        //Get 2D superPhi segments from 4D segments
        const DTChamberRecSegment2D* phiSegment2D = (*segment4D).phiSegment();
        if((*phiSegment2D).dimension() != 2) {
          if(debug) cout << "[DTSegment2DQuality]***Error: This is not 2D segment!!!" << endl;
          abort();
        }

        // Segment Local Direction and position (in Chamber RF)
        LocalVector recSegDirection = (*phiSegment2D).localDirection();

        float recSegAlpha = DTHitQualityUtils::findSegmentAlphaAndBeta(recSegDirection).first;
        if(debug)
          cout << "  RecSegment direction: " << recSegDirection << endl
            << "             position : " <<  (*phiSegment2D).localPosition() << endl
            << "             alpha    : " << recSegAlpha << endl;

        if(fabs(recSegAlpha - angleSimSeg) < deltaAlpha) {
          deltaAlpha = fabs(recSegAlpha - angleSimSeg);
          bestRecHit = &(*phiSegment2D);
          bestRecHitFound = true;
        }
      }  // End of Loop over all 4D RecHits of this chambers

      if(bestRecHitFound) {
        // Best rechit direction and position in Chamber RF
        LocalPoint bestRecHitLocalPos = bestRecHit->localPosition();
        LocalVector bestRecHitLocalDir = bestRecHit->localDirection();

        LocalError bestRecHitLocalPosErr = bestRecHit->localPositionError();
        LocalError bestRecHitLocalDirErr = bestRecHit->localDirectionError();

        float angleBestRHit = DTHitQualityUtils::findSegmentAlphaAndBeta(bestRecHitLocalDir).first;
        if(fabs(angleBestRHit - angleSimSeg) < 5*sigmaResAngle &&
           fabs(bestRecHitLocalPos.x() - posSimSeg) < 5*sigmaResPos) {
          recHitFound = true;
        }

        // Fill Residual histos
        h2DHitSuperPhi->Fill(angleSimSeg,
                            angleBestRHit,
                            posSimSeg,
                            bestRecHitLocalPos.x(),
                            etaSimSeg,
                            phiSimSeg,
                            sqrt(bestRecHitLocalPosErr.xx()),
                            sqrt(bestRecHitLocalDirErr.xx())
                           );
      }
    } //end of if(nsegm!=0)

      // Fill Efficiency plot
    h2DHitEff_SuperPhi->Fill(etaSimSeg,
                            phiSimSeg,
                            posSimSeg,
                            angleSimSeg,
                            recHitFound);
  } // End of loop over chambers
}

void DTSegment2DSLPhiQuality::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 67 of file DTSegment2DSLPhiQuality.cc.

References interactiveExample::theFile.

                                     {
  // Write the histos to file
  theFile->cd();

  h2DHitSuperPhi->Write();

  h2DHitEff_SuperPhi->ComputeEfficiency();
  h2DHitEff_SuperPhi->Write();

  theFile->Close();
} 

Member Data Documentation

bool DTSegment2DSLPhiQuality::debug

private

Definition at line 50 of file DTSegment2DSLPhiQuality.h.

HEff2DHit* DTSegment2DSLPhiQuality::h2DHitEff_SuperPhi

private

Definition at line 62 of file DTSegment2DSLPhiQuality.h.

HRes2DHit* DTSegment2DSLPhiQuality::h2DHitSuperPhi

private

Definition at line 61 of file DTSegment2DSLPhiQuality.h.

std::string DTSegment2DSLPhiQuality::rootFileName

private

Definition at line 52 of file DTSegment2DSLPhiQuality.h.

std::string DTSegment2DSLPhiQuality::segment4DLabel

private

Definition at line 55 of file DTSegment2DSLPhiQuality.h.

double DTSegment2DSLPhiQuality::sigmaResAngle

private

Definition at line 59 of file DTSegment2DSLPhiQuality.h.

double DTSegment2DSLPhiQuality::sigmaResPos

private

Definition at line 57 of file DTSegment2DSLPhiQuality.h.

std::string DTSegment2DSLPhiQuality::simHitLabel

private

Definition at line 54 of file DTSegment2DSLPhiQuality.h.

TFile* DTSegment2DSLPhiQuality::theFile

private

Definition at line 48 of file DTSegment2DSLPhiQuality.h.