Ignore:
Timestamp:
08/02/05 13:41:41 (19 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/datacenter/macros/fillstar.C

    r7228 r7252  
    105105    }
    106106
    107     TProfile *h1 = (TProfile*)arr.FindObjectInCanvas("RingBroadening", "TProfile", "MHMuonPar");
    108     if (!h1)
    109     {
    110         cout << "WARNING - Reading of RingBroadening failed." << endl;
    111         return 0;
    112     }
    113 
    114     Float_t numerator = 0;
    115     Float_t denominator = 0;
    116     Float_t summe = 0;
    117 
    118     for(Int_t i=5; i<15; i++)
    119     {
    120         numerator   += h1->GetBinContent(i)*h1->GetBinEntries(i);
    121         denominator += h1->GetBinEntries(i);
    122     }
    123 
    124     summe = numerator/denominator;
    125 
    126     Float_t psf = (summe - 0.04816)/0.001294;
     107    MHMuonPar *hmuon = (TProfile*)arr.FindObjectInCanvas("MHMuonPar", "MHMuonPar", "MHMuonPar");
     108    if (!hmuon)
     109    {
     110        cout << "WARNING - Reading of MHMuon failed." << endl;
     111        return 0;
     112    }
     113
     114    Float_t psf   = (muon->GetMeanWidth() - 0.04816)/0.001294;
    127115    psf = TMath::Nint(psf*10)/10.;
    128116    TString PSF = Form("%5.1f", psf);
    129     Int_t   num = (int)h1->GetEntries();
    130 
    131     TProfile *h2 = (TProfile*)arr.FindObjectInCanvas("SizeVsRadius", "TProfile", "MHMuonPar");
    132     if (!h1)
    133     {
    134         cout << "WARNING - Reading of SizeVsRadius failed." << endl;
    135         return 0;
    136     }
    137 
    138     numerator = 0;
    139     denominator = 0;
    140     summe = 0;
    141 
    142     for(Int_t i=5; i<15; i++)
    143     {
    144         numerator   += h2->GetBinContent(i)*h2->GetBinEntries(i);
    145         denominator += h2->GetBinEntries(i);
    146     }
    147 
    148     summe = numerator/denominator;
     117    Int_t   num = (int)muon.GetEntries();
    149118
    150119    Float_t integralmc = -1.43*psf + 1035;
    151     Float_t ratiodatamc = (summe/integralmc)*100;
     120    Float_t ratiodatamc = (muon->GetMeanSize()/integralmc)*100;
    152121    TString ratio = Form("%5.1f", ratiodatamc);
    153 
    154122
    155123    TH1 *h = (TH1*)arr.FindObjectInCanvas("Islands", "TH1F", "MHImagePar");
Note: See TracChangeset for help on using the changeset viewer.