Ignore:
Timestamp:
10/24/01 15:12:21 (23 years ago)
Author:
tbretz
Message:
*** empty log message ***
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/mmain/MMonteCarlo.cc

    r954 r988  
    219219    //
    220220    MReadTree reader("Events", fInputFile);
     221    reader.UseLeaf("fImpact");
     222    reader.UseLeaf("fEnergy");
     223    reader.UseLeaf("fNumFirstLevel");
     224
    221225    tlist.AddToList(&reader);
    222226
     
    288292    // Check if the list really contains the right number of histograms
    289293    //
    290     if (hists.GetEntriesFast() != dim)
     294    if (hists.GetEntriesFast() != dim && dim)
    291295        return;
    292296
     
    304308    //
    305309    MReadTree reader("Events", fInputFile);
     310    reader.UseLeaf("fImpact");
     311    reader.UseLeaf("fEnergy");
     312    reader.UseLeaf("fPhi");
     313    reader.UseLeaf("fTheta");
     314    reader.UseLeaf("fNumFirstLevel");
     315    reader.UseLeaf("fPhotElfromShower");
     316
    306317    tlist.AddToList(&reader);
    307318
    308     Float_t BgR[10]={660,4,0,0,0,0,0,0,0,0};
    309 
    310     MMcTriggerRateCalc crate(dim, 14, BgR, 100000, 2.75, 10.91e-2);
     319    Float_t BgR[10]={660, 4, 0, 0, 0, 0, 0, 0, 0, 0};
     320
     321    MMcTriggerRateCalc crate(dim, 14, BgR, 100000);
    311322    tlist.AddToList(&crate);
    312323
     
    323334        return;
    324335
    325     TIter Next(&hists);
    326     MHMcRate *rate=NULL;
    327     while ((rate=(MHMcRate*)Next()))
    328         rate->Print();
     336    hists.Print();
    329337}
    330338
     
    378386    //      like one dimension MMcThresholdCalc
    379387    //
    380     MReadTree        read("Events", fInputFile);
     388    MReadTree read("Events", fInputFile);
     389    read.VetoBranch("MRawEvtData");
     390    read.VetoBranch("MRawEvtHeader");
     391
    381392    MMcThresholdCalc calc(dim);
    382 
    383393    tlist.AddToList(&read);
    384394    tlist.AddToList(&calc);
     
    398408    //
    399409    TIter Next(&hists);
    400     MHMcRate *hist=NULL;
    401     while ((hist=(MHMcRate*)Next()))
    402         hist->DrawClone();
     410    TObject *obj;
     411    while ((obj=Next()))
     412        obj->DrawClone();
    403413}
    404414
Note: See TracChangeset for help on using the changeset viewer.