Test beam "publish"(performance) plot - PS data.

To make progress with the energy scale and non-linearity we need to "publish" a few things from the test beam analysis, to be included in an EMCal analysis note. These are things which must be compared to MC, to be sure that we have the MC under control. We have data for 10 values of beam momenta - 0.5, 0.75, 1, 1.25, 1.5, 2, 3, 4, 5 and 6 GeV.

Cell threshold is 2 adc units (~30 MeV)

The number of cells per cluster at each momenta.

  • pdf file: Presentation for EMCAL phone meeting at Jan 20,2011 by Alexei Pavlinov.

You can read the list of histograms with next peace of code:

    TFile f("PS_th2_NcellsInCluster.root","READ");
    TIter nextkey(f.GetListOfKeys());
    int nKeys = f.GetListOfKeys()->GetSize();
    printf("\t Keys %i \n", nKeys); // Should be one
    TKey *key=(TKey*)nextkey(); // get first key (only one exist) 
    TObject *obj = key->ReadObj();
    TList *l = dynamic_cast<TList *>(obj);
    TH1 *h = dynamic_cast<TH1 *>(l->At(i));  // i is index of hist

-- AlexeiPavlinov - 21-Jan-2011

Edit | Attach | Watch | Print version | History: r4 < r3 < r2 < r1 | Backlinks | Raw View | Raw edit | More topic actions...
Topic revision: r1 - 2011-01-21 - AlexeiPavlinov
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback