MuonInsideOutAnalysisPlots.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef MUON_MUONINSIDEOUTANALYSISPLOTS_H
#define MUON_MUONINSIDEOUTANALYSISPLOTS_H
 
#include "MuonRecValidationNtuples/MuonInsideOutValidationNtuple.h"
#include "TString.h"
 
class TH1;
class TH2;
class TDirectory;
 
namespace Muon {
  
  struct ResPlots {
    TH1* res;
    TH1* pull;
    TH1* exerr;
    
    void book( const TString& prefix );
    void fill( float res_, float pull_, float exerr_ );
  };

  struct ResPlotsPt {
    TH2* res_p;
    TH2* pull_p;
    TH2* exerr_p;
    
    void book( const TString& prefix );
    void fill( float res_, float pull_, float exerr_, float pt );
  };

  struct ChamberResPlots {
 
    std::vector<ResPlots>   plots;
    std::vector<ResPlotsPt> plots_p;
 
    void book( TDirectory* dir, const TString& prefix );
    void fill( int chIndex_, float res_, float pull_, float exerr_, float p_ );
 
  };

  struct HitPlots {
    ResPlots         all;
    ResPlotsPt       all_p;
    ChamberResPlots  chamber;
 
    void book( TDirectory* dir, const TString& prefix );
 
    void fill( const MuonValidationHitBlock& hits, int index, float p_ );
 
    void fill( int chIndex_, float res_, float pull_, float exerr_, float p_ );
  };

  struct TimePlots {
    int typeIndex;
    TH1* time;
    TH1* error;
    TH1* timeRes;
    TH1* timePull;
    TH1* beta;
    TH1* betaRes;
    TH1* betaPull;
    TH1* invBetaRes;
    TH1* invBetaPull;
    TH2* timeResScat;
    TH2* timePullScat;
 
    TH1* error_time;
    TH1* errorRT;
    TH1* errorTrack;
    TH2* errorRT_Track;
    TH1* timeResCor;
    TH1* timePullCor;
    TH1* timeRes0_10;
    TH1* timeRes10_70;
    TH1* timeRes70_150;
    TH1* timeRes150_350;
    TH1* timeRes350_800;
 
    TH1* timePull0_10;
    TH1* timePull10_70;
    TH1* timePull70_150;
    TH1* timePull150_350;
    TH1* timePull350_800;
 
    void book( TDirectory* dir, const TString& prefix, int type );
 
    void fill( const MuonValidationTimeBlock& times, int index, float betaTrack );
  };

  struct HoughPlots {
    TH1* max;
    TH2* max_pull;
 
    ResPlots         all;
    ChamberResPlots  chamber;
 
    void book( TDirectory* dir, const TString& prefix );
    void fill( const MuonValidationHoughBlock& hough, int index, float p_ );
    void fill( int chIndex_, float res_, float pull_, float exerr_, float p_, float max_ );
 
  };

  struct SegmentPlots {
    TH1* t0;
    TH1* t0Trig;
    TH1* beta;
    TH2* beta2;
    TH1* betaTrig;
    TH1* t0Res;
    TH1* t0ResTrig;
    TH1* betaRes;
    TH1* betaResTrig;
    TH1* quality;
    TH2* quality_chIndex;
    TH2* quality_sector;
    
    ResPlots        allx;
    ChamberResPlots chamberx;
 
    ResPlots        ally;
    ChamberResPlots chambery;
 
    ResPlots        allxz;
    ChamberResPlots chamberxz;
 
    ResPlots        allyz;
    ChamberResPlots chamberyz;
 
    ResPlots        allcy;
    ChamberResPlots chambercy;
 
    void book( TDirectory* dir, const TString& prefix );
    void fill( int chIndex_, int sector_, int quality_ );
    void fill( const MuonValidationSegmentBlock& segments, int index, float p_, float betaTrack );
 
  };
 

  struct BetaFitPlots {
    TH1* beta;
    TH1* chi2Ndof;
    TH1* ndof;
    TH1* res;
 
    void book( TDirectory* dir, const TString& prefix );
    void fill( float beta_, float betaTruth_, float chi2_, int ndof_ );
  };
 
  struct BetaFitRegionPlots {
    BetaFitPlots mdt;
    BetaFitPlots mdtt;
    BetaFitPlots mdtt_good;
    BetaFitPlots rpc;
    BetaFitPlots rpct;
    BetaFitPlots all;
    BetaFitPlots allt;
    BetaFitPlots allt_good;
    void book( TDirectory* dir, const TString& prefix, bool isBarrel );
  };

  struct CandidatePlots {
    TH1* pt;
    TH1* eta;
    TH1* phi;
    TH1* beta;
    TH1* nseg;
 
    BetaFitPlots betaCandidates;
    BetaFitPlots betaBestCandidate;
 
    void book( TDirectory* dir, const TString& prefix );
    void fill( int ntruth_, int nseg_, int nseg1_, int nhough_ );
    void fill( const MuonValidationTrackParticleBlock& tracks, int index );
 
  };
 
  struct StageSummaryPlots {
    TH1* betaMissed;
    TH1* betaMissedCombined;
    TH1* etaMissed;
    TH2* etaBetaMissed;
    TH1* etaMissedCombined;
    TH1* ncandidates;
    TH1* ncombinedCandidates;
    CandidatePlots allCandidates;
    CandidatePlots tagCandidates;
    CandidatePlots combinedCandidates;
    void book( TDirectory* dir, const TString& prefix );
  };

  struct TrackPlots {
    TH1* pt;
    TH1* eta;
    TH1* phi;
    TH1* beta;
    TH1* pdg;
    TH1* ntruth;
    TH1* nseg;
    TH1* nseg1;
    TH1* nseg2;
    TH1* nseg3;
    TH1* nhough;
    TH2* ntruth_seg;
    TH2* ntruth_seg1;
    TH2* ntruth_hough;
    
    HitPlots     hits;
    HoughPlots   hough;
    SegmentPlots segments;
    SegmentPlots segments1;
    SegmentPlots segments2;
    SegmentPlots segments3;
 
    std::vector<TimePlots> timePlots;
    // get the plot for a given index
    TimePlots* getTimePlots( int typeIndex ) {
      for( auto& plot : timePlots ){
        if( plot.typeIndex == typeIndex ) return &plot;
      }
      return nullptr;
    }
 
    BetaFitRegionPlots barrel;
    BetaFitRegionPlots endcap;
    
    std::vector<StageSummaryPlots> candidateStages;
    
    void book( TDirectory* dir, const TString& prefix );
    void fill( int ntruth_, int nseg_, int nseg1_, int nseg2_, int nseg3_, int nhough_ );
    void fill( const MuonValidationTrackParticleBlock& tracks, int index );
  };
 
  struct MuonInsideOutAnalysisPlots {
    TrackPlots muon;
    TrackPlots stau;
    TrackPlots rest;
    void book( TDirectory* dir, const TString& prefix );
  };
 
}
 
#endif