Ntuple Plots

distributions

Muons Occupancy

Occupancy dei Muons sui dati 50908Global


Riempimento dentro il ciclo sui muoni:
OccupancyMuons.Fill(treeVars.muonTkAtEcalPhi[MUindex], treeVars.muonTkAtEcalEta[MUindex]);


Occupancy dei Muons sui dati 50908Cosmic in arrivo

Come sono riempite eta e phi per i muoni

nel codice di produzione delle ntuple (EcalCosmicTree.cc)

dentro al ciclo sulla collezione di tracce
for(reco::TrackCollection::const_iterator recoTrack = recoTracks->begin(); recoTrack = recoTracks->end(); ++recoTrack)

TrackDetMatchInfo info = trackAssociator_.associate(iEvent, iSetup, *recoTrack, trackParameters_);
myTreeVariables_.muonTkAtEcalEta[numberOfRecoMuons] = (float)info.trkGlobPosAtEcal.eta();
myTreeVariables_.muonTkAtEcalPhi[numberOfRecoMuons] = (float)info.trkGlobPosAtEcal.phi();

Muons on Inner Ecal Barrel function ect::EtaPhiMuonAtEcal.

Analytical intersection between muon track and Ecal Inner Barrel

No intersections in Ecal Inner Barrel


Two intersections in Ecal Inner Barrel


One intersection in Ecal Inner Barrel: Barrel border

Inner vs outer Hits

Muon Phi intersection with ECAL before and after the use of the appropriate TrackAssociator tag

Distributions refer to cosmic (first plot) and global(second plot) muon reconstruction, without (in red) and with (in black) the use of the TrackAssociator tag V03-01-03-01.

Asymmetry in ECAL intersections (as returned by the TrackAssociator) with the use of globalCosmicMuons still to be understood.

maps

muon length in single crystals

muon length in ECAL

distances most energetic - mostly traversed

(Total Tk in Ecal - SUM over Xtals Tks) / Total Tk in Ecal

Energy vs Length

number of muons and superclusters

The number of recoMuons il plotted against the number of SuperClusters, both in the case of the cosmic muons and of the global muons, for the run number 50908, from data in /castor/cern.ch/user/m/mattia/50908Global and /castor/cern.ch/user/m/mattia/50908Cosmic respectively.

The values superimposed correspond to the number of events and the fraction of events, respectively.

The macro to create the plots:

TFile *_file0 = TFile::Open("testGlobal50908.root")
SCMuons->SetTitle ("global muons run 50908")
SCMuons->GetXaxis ()->SetTitle ("nSuperClusters")
SCMuons->GetYaxis ()->SetTitle ("nRecoMuons")
SCMuons->DrawCopy ("COLZ")
SCMuons->DrawNormalized("textsame")
/* oppure SCMuons->Draw("textsame") */

superclusters distribution map

The number of superclusters found in the (phi,eta) ECAL surface, both for the Cosmic and Global reconstruction respectively, data from run 50908. Plots should change only because of different statistiscs in the two cases.

muons distribution map

The number of muons found in the (phi,eta) ECAL surface, both for the Cosmic and Global reconstruction respectively, data from run 50908. ECAL surface defined by means of EtaPhiMuonAtEcal

matching between superclusters and muons

matched superclusters and muons distribution map

The number of superclusters found in the (phi,eta) ECAL surface, both for the Cosmic and Global reconstruction respectively, data from run 50908. Only superclusters that match with muons, within a cone 0.3, are plotted.

The same plot for the muons inrtersection at ECAL, calculated by means of:

          const GlobalPoint muonInnerHit (treeVars.muonInnerHitX[MUindex],
                                          treeVars.muonInnerHitY[MUindex],
                                          treeVars.muonInnerHitZ[MUindex]);
          const GlobalPoint muonOuterHit (treeVars.muonOuterHitX[MUindex],
                                          treeVars.muonOuterHitY[MUindex],
                                          treeVars.muonOuterHitZ[MUindex]);
    
          std::pair<GlobalPoint, GlobalPoint> muonTkAtEcal = 
                ect::EtaPhiMuonAtEcal(muonInnerHit, muonOuterHit) ;

          if (muonTkAtEcal.first.perp () < 0.001) continue ;
          matchMUdistr.Fill (muonTkAtEcal.first.phi (),
                             muonTkAtEcal.first.eta ()) ;

matching efficiency

with the same conditions as above, distribution as a function of the muons number:

with the same conditions as above, distribution as a function of the superclusters number:

matching distance

with the same conditions as above, DeltaR for cosmics and global:

matching distance

matching distance vs Eta, vs Phi

dE/dx Approx vs Detail

Differenza relativa 2* (dE/dx(Detailed)-dE/dx(Approx)) / (dE/dx(Detailed) + dE/dx(Approx) )

Selections

Angle between muon direction and supercluster axis for matched muon-SC

dEdx with variation of the cut on Angle

cuts on 1 rad, 0.8 rad 0.6 rad 0.4 rad 0.2 rad

Occupancy Xtals

Occupancy all Xtals in matched SC: all events , events angle(MUON-SC) < pi/8

Average Energy vs Eta, Ave Energy vs Phi

All events , events events angle(MUON-SC) < pi/8

Correlation between Angles and SC Energies

Nel grafico seguente in bianco Energia dei SC su tutti gli eventi, in giallo energia dei SC con angolo minore di pi/8. Le statistiche si riferiscono al sample su tutti gli eventi (qualcuno sa come mettere due box di statistiche sullo stesso canvas?).

Nel grafico seguente in bianco Angoli su tutti gli eventi matchati, in verde angoli con cut su 0.2<E<0.5 GeV

CUT su Xtal MAX, runs 51047, 51087, 51020

Correlazione con Angolo Muone-SCDirection

Nero: Tutti eventi dei run 51047, 51087, 51020

CUTS su cristallo più energetico Rosso: 50Mev Blu: 75 Mev Viola: 100 Mev Verde: 125 Mev Azzurro 150 Mev

Occupancy

Totale 50MeV

75Mev 100MeV

125Mev 150MeV

CUTS on Angle: runs 51057, 51087, 51020, 50908, 50911

-- PietroGovoni - 22 Aug 2008