|
| gROOT | Macro ("$VMCWORKDIR/gconfig/rootlogon.C") |
| |
| gROOT | LoadMacro ("$VMCWORKDIR/macro/mvd/Tools.C") |
| |
| | LoadPandaStyle () |
| |
| t | AddFriend ("pndsim","data/digi_sttcombi.root") |
| |
| t | AddFriend ("pndsim","data/reco_sttcombi.root") |
| |
| t | AddFriend ("pndsim","data/lhepidmaker.root") |
| |
| t | AddFriend ("pndsim","data/fit_sttcombi.root") |
| |
| t | SetBranchAddress ("MCTrack",&mctrack_array) |
| |
| t | SetBranchAddress ("MVDPoint",&mc_array) |
| |
| t | SetBranchAddress ("MVDPixelDigis",&digiPixel_array) |
| |
| t | SetBranchAddress ("MVDStripDigis",&digiStrip_array) |
| |
| t | SetBranchAddress ("MVDStripClusterCand",&stripClust_array) |
| |
| t | SetBranchAddress ("MVDClusterCand",&pixelClust_array) |
| |
| t | SetBranchAddress ("MVDHitsStrip",&stripHit_array) |
| |
| t | SetBranchAddress ("MVDHitsPixel",&pixelHit_array) |
| |
| t | SetBranchAddress ("LhePidTrack",&track_array) |
| |
| c1 | Divide (4, 2) |
| |
| c1 | cd (1) |
| |
| t | Draw ("LhePidTrack.GetMomentum().Mag()>>momlhe(100,0.5,1.5)","","") |
| |
| momlhe | SetTitle ("LhePidTrack;Momentum / GeV/c") |
| |
| c1 | cd (2) |
| |
| t | Draw ("LhePidTrack.GetMomentum().Mag()-MCTrack.GetMomentum().Mag()>>momlhedf(100,-2.5,2.5)","","") |
| |
| momlhedf | SetTitle ("LhePidTrack - MC;Momentum / GeV/c") |
| |
| c1 | cd (3) |
| |
| t | Draw ("LheGenTrack.getMom().Mag()>>momfit(100,0.5,1.5)","","") |
| |
| momfit | SetTitle ("genfit Track;Momentum / GeV/c") |
| |
| c1 | cd (4) |
| |
| t | Draw ("LheGenTrack.getMom().Mag()-MCTrack.GetMomentum().Mag()>>momfitdf(100,-2.5,2.5)","","") |
| |
| momfitdf | SetTitle ("genfit Track - MC;Momentum / GeV/c") |
| |
| c1 | cd (5) |
| |
| t | Draw ("MCTrack.GetMomentum().Mag():LhePidTrack.GetMomentum().Mag()-MCTrack.GetMomentum().Mag()>>momlhedfcor(100,-2.5,2.5,100,0.5,1.5)","","colz") |
| |
| momlhedfcor | SetTitle (";LHE - MC Momentum / GeV/c;MC Momentum / GeV/c") |
| |
| c1 | cd (6) |
| |
| t | Draw ("MCTrack.GetMomentum().Mag():LhePidTrack.GetMomentum().Mag()>>momlhecor(100,0.5,1.5,100,0.5,1.5)","","colz") |
| |
| momlhecor | SetTitle (";LHE Momentum / GeV/c;MC Momentum / GeV/c") |
| |
| c1 | cd (7) |
| |
| t | Draw ("MCTrack.GetMomentum().Mag():LheGenTrack.getMom().Mag()-MCTrack.GetMomentum().Mag()>>momfitdfcor(100,-2.5,2.5,100,0.5,1.5)","","colz") |
| |
| momfitdfcor | SetTitle (";GEN - MC Momentum / GeV/c;MC Momentum / GeV/c") |
| |
| c1 | cd (8) |
| |
| t | Draw ("MCTrack.GetMomentum().Mag():LheGenTrack.getMom().Mag()>>momfitcor(100,0.5,1.5,100,0.5,1.5)","","colz") |
| |
| momfitcor | SetTitle (";GEN Momentum / GeV/c;MC Momentum / GeV/c") |
| |
| c2 | Divide (2, 2) |
| |
| t | Draw ("1000*LhePidTrack.GetMvdELoss()/LhePidTrack.GetMvdHitCounts():LhePidTrack.GetMomentum().Mag()>>mvd_pid(100,0,3,100,0,2)","LhePidTrack.GetMvdHitCounts()>0","colz") |
| |
| mvd_pid | SetTitle ("MVD pid;Momentum [GeV/c];MVD dE / # hit [MeV]") |
| |
| t | Draw ("LhePidTrack.GetBeta():LhePidTrack.GetMomentum().Mag()>>tof_pid(100,0,2,100,0,1.2)","LhePidTrack.GetTofIndex()>-1","colz") |
| |
| tof_pid | SetTitle ("TOF pid;Momentum [GeV/c];#beta") |
| |
| t | Draw ("LhePidTrack.GetDrcThetaC()*TMath::RadToDeg():LhePidTrack.GetMomentum().Mag()>>drc_pid(100,0,3,100,0,60)","LhePidTrack.GetDrcIndex()>-1","colz") |
| |
| drc_pid | SetTitle ("DRC pid;Momentum [GeV/c];#theta_{c} [^{0}]") |
| |
| t | Draw ("LhePidTrack.GetEmcELoss()/LhePidTrack.GetMomentum().Mag():LhePidTrack.GetMomentum().Mag()>>emc_pid(100,0,3,100,0,1.2)","LhePidTrack.GetEmcIndex()>-1","colz") |
| |
| emc_pid | SetTitle ("EMC pid;Momentum [GeV/c];EMC Energy / Momentum [c]") |
| |
|
| Int_t | nEvents = 100 |
| |
| TFile * | f = new TFile("data/points_sttcombi.root") |
| |
| TTree * | t =(TTree*)f->Get("pndsim") |
| |
| TClonesArray * | mctrack_array =new TClonesArray("PndMCTrack") |
| |
| TClonesArray * | mc_array =new TClonesArray("PndSdsMCPoint") |
| |
| TClonesArray * | digiPixel_array =new TClonesArray("PndSdsDigiPixel") |
| |
| TClonesArray * | digiStrip_array =new TClonesArray("PndSdsDigiStrip") |
| |
| TClonesArray * | stripClust_array =new TClonesArray("PndSdsCluster") |
| |
| TClonesArray * | pixelClust_array =new TClonesArray("PndSdsCluster") |
| |
| TClonesArray * | stripHit_array =new TClonesArray("PndSdsHit") |
| |
| TClonesArray * | pixelHit_array =new TClonesArray("PndSdsHit") |
| |
| TClonesArray * | track_array =new TClonesArray("PndLhePidTrack") |
| |
| TGeoManager * | geoMan = (TGeoManager*) gDirectory->Get("FAIRGeom") |
| |
| PndMvdGeoHandling * | fGeoH = new PndMvdGeoHandling(geoMan) |
| |
| TH1D * | hMomDiff = new TH1D("hmomdiff","Mom difference track to mctrack;|p-p_{mc}|/ GeV/c;",500,-15.,15.) |
| |
| TCanvas * | c1 = new TCanvas("LHETRACK MOMENTIM", "LHETRACK MOMENTUM", 10, 10, 1200, 600) |
| |
| TCanvas * | c2 = new TCanvas("LHETRACK PID", "LHETRACK PID", 20, 20, 600, 600) |
| |
