Changeset 1837


Ignore:
Timestamp:
03/19/03 12:41:28 (22 years ago)
Author:
moralejo
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/MagicSoft/Mars/Changelog

    r1836 r1837  
    11                                                 -*-*- END OF LINE -*-*-
     2
     3 2003/03/19: Abelardo Moralejo
     4
     5    * macros/CT1collarea.C:
     6      - Added filter to cut in hadronness (now available in new version of
     7        root CT1 Monte Carlo file from Wolfgang).
    28
    39 2003/03/19: Robert Wagner
     
    612      - changed MH3::New() such that constructor according to Thomas'
    713        previous changes is supported
    8 
    914
    1015
  • trunk/MagicSoft/Mars/macros/CT1collarea.C

    r1823 r1837  
    2424
    2525
    26 void CT1collarea(TString filename="MC_ON1.root", TString outname="")
     26void CT1collarea(TString filename="MC_ON2.root", TString outname="")
    2727{
    2828    //
     
    4040    //
    4141    MReadMarsFile reader("Events", filename);
     42    reader.DisableAutoScheme();
     43    //    reader.EnableBranch("fHadronness");
     44
    4245    MMcCT1CollectionAreaCalc effi;
     46    MHMcCT1CollectionArea* collarea = new MHMcCT1CollectionArea("MHMcCT1CollectionArea","",30,2.,5.);
     47    parlist.AddToList(collarea);
    4348
    4449    tasklist.AddToList(&reader);   
     50
     51    MF filterhadrons("MHadronness.fHadronness<0.25");
     52    tasklist.AddToList(&filterhadrons);
     53    effi.SetFilter(&filterhadrons);
    4554
    4655    tasklist.AddToList(&effi);
     
    7584        return;
    7685
    77     TFile f(outname, "recreate");
    78     if (!f)
    79         return;
    80     parlist.FindObject("MHMcCT1CollectionArea")->Write();
     86    TFile f("area.root","recreate");
     87    collarea->GetHist()->Write();
     88    collarea->GetHAll()->Write();
     89    collarea->GetHSel()->Write();
     90    f.Close();
    8191}
    8292
Note: See TracChangeset for help on using the changeset viewer.