Ignore:
Timestamp:
02/10/05 16:06:21 (20 years ago)
Author:
mazin
Message:
*** empty log message ***
Location:
trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses
Files:
1 added
4 edited

Legend:

Unmodified
Added
Removed

  • trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MFindSupercutsONOFF.cc

    r6232 r6355  
    319319   
    320320
     321    // Degree of polynomials used to fit the ON and the OFF data are retrieved from
     322    // the event loop
     323
     324    Int_t degreeON;
     325    Int_t degreeOFF;
     326
     327    MDataValue* DegreeONContainer = (MDataValue*) ONDataplistfcn->FindObject("DegreeON");
     328    degreeON = int(DegreeONContainer -> GetValue());
     329
     330    MDataValue* DegreeOFFContainer = (MDataValue*) OFFDataplistfcn->FindObject("DegreeOFF");
     331    degreeOFF = int(DegreeOFFContainer -> GetValue());
     332
     333    gLog << "fcnSupercuts : fDegreeON and fDegreeOFF retrieved from ONDataplistfcn and OFFDataplistfcn: " << endl
     334         << "fDegreeON, fDegreeOFF: = " << degreeON << " ," << degreeOFF << endl;
     335
     336
    321337    Double_t AlphaSig = 0.0;
    322338    MDataValue* AlphaSigContainer = (MDataValue*) ONDataplistfcn->FindObject("AlphaSigValue");
    … …  
    484500   // alpha min for bkg region in ON data
    485501    const Double_t alphamax = AlphaBkgMax; // alpha max for bkg region in ON data
    486     const Int_t    degreeON   =    2;
    487     // degree of polynomial used to fit ON data in Bkg  region
    488     const Int_t    degreeOFF   =    2;
    489     // degree of polynomial used to fit OFF data in ALL region
     502
    490503
    491504    Bool_t drawpoly;
    … …  
    523536   
    524537    TString DummyPs = ("Dummy");
     538 
    525539
    526540    const Bool_t rc = findsig.FindSigmaONOFF(&alphaHist,&alphaHistOFF,
    … …  
    788802    ////////////////////////////////////////////////////
    789803   
    790 
    791 
     804   
     805    // Degree of polynomials used to fit the ON OFF data is
     806    // set initially to 2.
     807
     808    fDegreeON = 2;
     809    fDegreeOFF = 2;
    792810
    793811
    … …  
    19311949
    19321950
    1933 
    1934 
    1935 
    1936    
    1937 
    1938 
    1939 
    1940    
    19411951    Double_t prob = whichfractiontrain/(whichfractiontrain + whichfractiontest);
    19421952
    … …  
    24922502    AlphaBkgMaxCont.SetName("AlphaBkgMaxValue");
    24932503   
     2504
     2505    MDataValue DegreeONCont (fDegreeON);
     2506    DegreeONCont.SetName("DegreeON");
     2507
     2508    MDataValue DegreeOFFCont (fDegreeOFF);
     2509    DegreeOFFCont.SetName("DegreeOFF");
     2510
     2511
     2512
     2513
    24942514   
    24952515    // Value of fUseFittedQuantities is stored in container
    … …  
    25932613    parlistfcn.AddToList(&AlphaBkgMinCont);
    25942614    parlistfcn.AddToList(&AlphaBkgMaxCont);
     2615    parlistfcn.AddToList(&DegreeONCont);
    25952616    parlistfcn.AddToList(&UseFittedQuantitiesCont);
    25962617    parlistfcn.AddToList(&UseNormFactorFromAlphaBkgCont);
    … …  
    26982719    parlistfcnOFF.AddToList(&alphaOFF);
    26992720
     2721    parlistfcnOFF.AddToList(&DegreeOFFCont);
    27002722   
    27012723
    … …  
    35213543   
    35223544    const Double_t alphamax = fAlphaBkgMax;
    3523     const Int_t    degree   =    2;
    3524     const Int_t    degreeOFF   =    2;
     3545//    const Int_t    degree   =    2;
     3546//    const Int_t    degreeOFF   =    2;
    35253547    const Bool_t   drawpoly  = kTRUE;
    35263548    const Bool_t   fitgauss  = kTRUE;
    … …  
    35743596                                       fNormFactorTest,
    35753597                                       alphamin, alphamax,
    3576                                        degree, degreeOFF,
     3598                                       fDegreeON, fDegreeOFF,
    35773599                                       alphasig, drawpolyTest,
    35783600                                       fitgaussTest, printTest,
    … …  
    36843706                           fNormFactorTest,
    36853707                           alphamin, alphamax,
    3686                            degree, degreeOFF,
     3708                           fDegreeON, fDegreeOFF,
    36873709                           alphasig, drawpoly, fitgauss, print,
    36883710                           saveplots, psfilename);
    … …  
    42084230   
    42094231    const Double_t alphamax = fAlphaBkgMax;
    4210     const Int_t    degree   =    2;
    4211     const Int_t    degreeOFF   =    2;
     4232//    const Int_t    degree   =    2;
     4233//    const Int_t    degreeOFF   =    2;
    42124234    const Bool_t   drawpoly  = kTRUE;
    42134235    const Bool_t   fitgauss  = kTRUE;
    … …  
    42624284                                      fNormFactorTrain,
    42634285                                      alphamin, alphamax,
    4264                                       degree, degreeOFF,
     4286                                      fDegreeON, fDegreeOFF,
    42654287                                      alphasig, drawpolyTest,
    42664288                                      fitgaussTest, printTest,
    … …  
    43654387                           fNormFactorTrain,
    43664388                           alphamin, alphamax,
    4367                            degree, degreeOFF,
     4389                           fDegreeON, fDegreeOFF,
    43684390                           alphasig, drawpoly, fitgauss, print,
    43694391                           saveplots, psfilename);

  • trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MFindSupercutsONOFF.h

    r6232 r6355  
    6666
    6767
    68 
    69 
     68    Int_t fDegreeON;
     69    Int_t fDegreeOFF;
    7070
    7171
    … …  
    317317  Bool_t CheckAlphaSigBkg();
    318318
    319 
    320 
     319  void SetDegreeON (Int_t d) {fDegreeON = d;}
     320  void SetDegreeOFF (Int_t d) {fDegreeOFF = d;}
    321321
    322322

  • trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MHFindSignificanceONOFF.cc

    r6232 r6355  
    547547    //*fLog << "MHFindSignificanceONOFF::FindSigma;" << endl;
    548548 
    549 
    550549    // Pointer to object TPostScript where plots will be stored
    551550    // is copied into member varialbe
    … …  
    12151214      // prepare fit of a polynomial :   (a0 + a1*x + a2*x**2 + a3*x**3 + ...)
    12161215
     1216
    12171217      TString funcname = "PolyOFF";
    12181218      Double_t xmin =   0.0;
    … …  
    17211721      //--------------------------------------------------
    17221722      // prepare fit of a polynomial :   (a0 + a1*x + a2*x**2 + a3*x**3 + ...)
    1723 
     1723 
    17241724      TString funcname = "Poly";
    17251725      Double_t xmin =   0.0;
    … …  
    32343234        fHist -> SetTitle("Alpha Plot");
    32353235        fHist-> SetTitleOffset(1.5, "Y");
     3236        fHist-> SetLineColor(kRed);
    32363237        fHist-> DrawCopy();
    32373238       
    … …  
    32543255        }
    32553256
    3256         fHistOFFNormalized->SetLineColor(kRed);
    3257         fHistOFFNormalized->DrawCopy("same");
     3257
     3258//      fHistOFFNormalized->SetLineColor(kRed);
     3259//      fHistOFFNormalized->DrawCopy("same");
    32583260
    32593261    }
    … …  
    32783280        *fLog << "MHFindSignificanceONOFF::DrawFit; fbackg = NULL" << endl;
    32793281
     3282    if (fHistOFFNormalized)
     3283      {
     3284      fHistOFFNormalized->SetFillStyle(3002);
     3285      fHistOFFNormalized->SetFillColor(9);
     3286      fHistOFFNormalized->DrawCopy("Psamebar");}
     3287
     3288
    32803289      if (fbackg)
    32813290      {
    32823291        // 10, 1 is white and solid
    3283         fbackg->SetLineColor(kRed);
     3292        //fbackg->SetLineColor(kRed);
     3293        fbackg->SetLineColor(196);
    32843294        fbackg->SetLineStyle(1);
    32853295        fbackg->SetLineWidth(4);
    3286         //fbackg->DrawCopy("same"); // I do not want to draw it... already too many things.
     3296        fbackg->DrawCopy("same"); // I do not want to draw it... already too many things.
    32873297      }
    32883298    }
    32893299
     3300//    fCanvas->SaveAs("plot.root");
     3301//    fCanvas->SaveAs("plot.ps");
    32903302
    32913303    //-------------------------------
    … …  
    34553467        if (!fPsFilenameString.IsNull())
    34563468        {
    3457             TString filename = (fPsFilenameString);
     3469            TString filename  = (fPsFilenameString);
     3470            TString filename2 = (fPsFilenameString);
    34583471            // Train or Test Sample is specified outside
    34593472            // class MHFindSignificanceONOFF, and included in
    34603473            // fPsFilenameString
    34613474           
    3462             filename += ("AlphaPlot.ps");
     3475            filename  += ("AlphaPlot.ps");
     3476            filename2 += ("AlphaPlot.C");
    34633477            cout << filename << endl;   
    34643478            fCanvas -> SaveAs(filename);
     3479            fCanvas -> SaveAs(filename2);
    34653480        }
    34663481       

  • trunk/MagicSoft/Mars/mtemp/mmpi/SupercutsONOFFClasses/MSupercutsCalcONOFF.cc

    r6232 r6355  
    108108    // as external variables.
    109109
    110     fSizeOffset = 3000; // photons
     110    fSizeOffset = 3000; // still in photons !
    111111    fDistOffset = 0.9; // degrees 
    112112    fThetaOffset = 20; // degrees   NOT USED FOR THE TIME BEING
Note: See TracChangeset for help on using the changeset viewer.