Changeset 14943 for fact/tools


Ignore:
Timestamp:
02/21/13 14:53:33 (12 years ago)
Author:
Jens Buss
Message:
disttribution histogramms added to hist list and saved to root file
File:
1 edited

Legend:

Unmodified
Added
Removed

  • fact/tools/rootmacros/PulseTemplates/FCalcPulseTemplate.C

    r14811 r14943  
    9696//----------------------------------------------------------------------------
    9797int FCalcPulseTemplate(
    98     TString     InRootFileName        = "20120802templ_0mV_entries.root",
    99 //    TString     InRootFileName        = "20120309_032.root",
    100 //    TString     InputPath           = "analysis/analysis/FPulseTemplate/20120309_017/Overlay/",
    101     TString     InputPath           = "",
    102 //    TString     InputPath           = "analysis/analysis/FPulseTemplate/20120309_032/Overlay/",
    103     TString     OutputRootFileName     = "test.root",
    104     TString     OutPutPath          = "",
    105     int         firstpixel          = 0,
    106     int         npixel              = 1,
    107     int         pixelSetSize        = 200,
    108     int         maxPulseOrder       = 1,
    109 //    TString     histoOptions        = "SRM",
    110     bool        ProduceGraphic      = true,
    111     bool        stats               = true,
    112     bool        saveResults         = false,
    113 //    bool        fitdata             = false,
    114     bool        debugPixel          = true,
    115 //    int         refresh_rate        = 500,      //refresh rate for canvases
    116     int         verbosityLevel      = 1        // different verbosity levels can be implemented here
     98
     99        // ------------{ Files }------------
     100//        TString     InRootFileName        = "20120802_0mV_newSelectAlg.root",
     101        TString     InRootFileName          = "20120802_0mV_Stacking_G11B-2AWW8.root",
     102//        TString     InRootFileName          = "20120802_200mV_Stacking_G13B-2AWW8.root",
     103        TString     InputPath               = "analysis/PulsTemplateAnalysis/20120802/",
     104        TString     OutputRootFileName      = "20120802_0mV_Extraction_G11B-2AWW8_model2.root",
     105//        TString     OutputRootFileName      = "20120802_200mV_Extraction_G13B-2AWW8_model0.root",
     106        TString     OutPutPath              = "analysis/PulsTemplateAnalysis/20120802/",
     107        //    TString     InRootFileName        = "20120309_032.root",
     108        //    TString     InRootFileName        = "20120309_017.root",
     109        //    TString     InputPath           = "analysis/analysis/FPulseTemplate/20120309_017/Overlay/",
     110        //    TString     InputPath           = "analysis/analysis/FPulseTemplate/20120309_032/Overlay/",
     111        // ------------{ data quantity }------------
     112        int     firstpixel              = 0,
     113        int     npixel                  = 1,
     114        int     pixelSetSize            = 200,
     115        int     maxPulseOrder           = 6,
     116        // ------------{ analysis parameters }------------
     117        float   meanGain                = 11.,
     118        float   meanBsl                 = -2.5,
     119        int     bootstrapIt             = 10,
     120        int     pulseModell             = 0,
     121        float   pixelOverlayXaxisLeft   = 0,
     122        float   pixelOverlayXaxisRight  = 400,
     123
     124        // ------------{ program options }------------
     125        //    TString     histoOptions        = "SRM",
     126
     127        // show graphical output
     128        bool        ProduceGraphic      = true,
     129
     130        // show stats in plots
     131        bool        stats               = true,
     132
     133        // save data
     134        bool        saveResults         = false,
     135        //    bool        fitdata             = false,
     136
     137        //pause after each pixel
     138        bool        debugPixel          = true,
     139
     140        //refresh rate for canvases
     141        //    int         refresh_rate        = 500,
     142
     143        TString fitOption = "Q",
     144
     145        // different verbosity levels can be set here
     146        int         verbosityLevel      = 1
    117147        )
    118148{
    … …  
    265295    TH1D hAmplitude[maxPulseOrder];
    266296    TH1D hRiseTime[maxPulseOrder];
     297    TList* histoList = new TList();
    267298
    268299    //Lopp over pulse order configure distribution histograms
    … …  
    272303        TString title;
    273304        title = "hBSL";
    274         title.Append(i);
     305        title += "_";
     306        title += i;
    275307        hBsl[i].SetNameTitle(title, title);
    276         hBsl[i].SetBins(400, -2-0.005, 2-0.005);
    277 
     308        hBsl[i].SetBins(800, -4-0.005, 4-0.005);
     309        histoList->Add( &hBsl[i] );
    278310        //Asymptote
    279311        title = "hAsymptote";
    280         title.Append(i);
     312        title += "_";
     313        title += i;
    281314        hAsymptote[i].SetNameTitle(title, title);
    282315        hAsymptote[i].SetBins(60, -10.5, 49.5);
     316        histoList->Add( &hAsymptote[i] );
    283317
    284318        //T0
    285319        title = "hT0";
    286         title.Append(i);
     320        title += "_";
     321        title += i;
    287322        hT0[i].SetNameTitle(title, title);
    288323        hT0[i].SetBins(80, 39.5, 119.5);
     324        histoList->Add( &hT0[i] );
    289325
    290326        //Tau1
    291327        title = "hTau1";
    292         title.Append(i);
     328        title += "_";
     329        title += i;
    293330        hTau1[i].SetNameTitle(title, title);
    294331        hTau1[i].SetBins(10, -10.5, 9.5);
     332        histoList->Add( &hTau1[i] );
    295333
    296334        //Tau2
    297335        title = "hTau2";
    298         title.Append(i);
     336        title += "_";
     337        title += i;
    299338        hTau2[i].SetNameTitle(title, title);
    300339        hTau2[i].SetBins(40, 19.5, 59.5);
     340        histoList->Add( &hTau2[i] );
    301341
    302342        //Integral
    303343        title = "hIntegral";
    304         title.Append(i);
     344        title += "_";
     345        title += i;
    305346        hIntegral[i].SetNameTitle(title, title);
    306347        hIntegral[i].SetBins(500, (i+1)*-0.5, (i+1)*499.5);
     348        histoList->Add( &hIntegral[i] );
    307349
    308350        //Amplitude
    309351        title = "hAmplitude";
    310         title.Append(i);
     352        title += "_";
     353        title += i;
    311354        hAmplitude[i].SetNameTitle(title, title);
    312355        hAmplitude[i].SetBins(400, (i+1)*4.95, (i+1)*24.95);
     356        histoList->Add( &hAmplitude[i] );
    313357
    314358        //hRiseTime
    315359        title = "hRiseTime";
    316         title.Append(i);
     360        title += "_";
     361        title += i;
    317362        hRiseTime[i].SetNameTitle(title, title);
    318363        hRiseTime[i].SetBins(20, -0, 10);
     364        histoList->Add( &hRiseTime[i] );
    319365    }
    320366
    … …  
    326372
    327373    //Lopp over pulse order configure distribution histograms
     374    float roi = 1024.;
    328375    for (int i = 0; i<maxPulseOrder; i++)
    329376    {
    … …  
    331378        TString title;
    332379        title = "hMedianErrors";
    333         title.Append(i);
     380        title += "_";
     381        title += i;
    334382        hMedianErrors[i].SetNameTitle(title, title);
    335         hMedianErrors[i].SetBins(300, -0.5, 300.5);
     383        hMedianErrors[i].SetBins(roi, -0.5, roi + 0.5);
    336384
    337385        //hMeanErrors
    338386        title = "hMeanErrors";
    339         title.Append(i);
     387        title += "_";
     388        title += i;
    340389        hMeanErrors[i].SetNameTitle(title, title);
    341         hMeanErrors[i].SetBins(300, -0.5, 300.5);
     390        hMeanErrors[i].SetBins(roi, -0.5, roi + 0.5);
    342391
    343392        //hMaxErrors
    344393        title = "hMaxErrors";
    345         title.Append(i);
     394        title += "_";
     395        title += i;
    346396        hMaxErrors[i].SetNameTitle(title, title);
    347         hMaxErrors[i].SetBins(300, -0.5, 300.5);
     397        hMaxErrors[i].SetBins(roi, -0.5, roi + 0.5);
    348398    }
    349399
    … …  
    355405    grPheVsAmpl.SetNameTitle("grPheVsAmpl", "Distribution of Pulse Amplitude vs. Photon equivalent");
    356406    grPheVsAmpl.SetMarkerStyle(21);
     407    histoList->Add( &grPheVsAmpl );
    357408
    358409    TGraphErrors grPheVsIntegral;
    359410    grPheVsIntegral.SetNameTitle("grPheVsIntegral", "Distribution of Pulse Integral vs. Photon equivalent");
    360411    grPheVsIntegral.SetMarkerStyle(21);
     412    histoList->Add( &grPheVsIntegral );
    361413
    362414    TGraphErrors grPheVsTau2;
    363415    grPheVsTau2.SetNameTitle("grPheVsTau2", "Distribution of Tau2 vs. Photon equivalent");
    364416    grPheVsTau2.SetMarkerStyle(21);
     417    histoList->Add( &grPheVsTau2 );
    365418//-----------------------------------------------------------------------------
    366419// Variables
    … …  
    462515                       stats,
    463516                       "L",
    464                        70,   //pixelOverlayXaxisLeft /TODO: get it from the root file
    465                        230,  //pixelOverlayXaxisRight /TODO: get it from the root file
    466                        -1,   //bslMean   /TODO: get it from the root file
    467                        9,    //gain mean /TODO: get it from the root file
     517                       pixelOverlayXaxisLeft,   //TODO: get it from the root file
     518                       pixelOverlayXaxisRight,  //TODO: get it from the root file
     519                       meanBsl,   //bslMean     /TODO: get it from the root file
     520                       meanGain,    //gain mean /TODO: get it from the root file
    468521                       inputRootFile,
    469522                       outputRootFile
    … …  
    472525            if (breakout)   break;
    473526
    474             //Preparing Camera
    475             if (first_pass)
    476             {
    477                 if (verbosityLevel > 0)
    478                 {
    479                     cout << endl << "...preparing camera" << endl;
    480                 }
    481 
    482                 wholeCamera = new PixelSum(
    483                             "AllPixel",
    484                             1440,
    485                             maxPulseOrder,
    486                             verbosityLevel,
    487                             stats,
    488                             "C",
    489                             pixel[pixelID]->mPixelOverlayXaxisLeft,
    490                             pixel[pixelID]->mPixelOverlayXaxisRight ,
    491                             pixel[pixelID]->mBSLMean ,
    492                             pixel[pixelID]->mGainMean ,
    493                             outputRootFile
    494                             );
    495 
    496                 first_pass = false;
    497             }
    498527            //-------------------------------------
    499528            // Histogramms of Maximas in Overlay Spectra
    … …  
    521550                //-------------------------------------
    522551
     552//                pixel[pixelID]->Normalize2Dhistos(pulse_order);
     553
    523554                //from Maximum Overlay
    524555                if (verbosityLevel > 2)
    … …  
    550581                            );
    551582
     583                pixel[pixelID]->SetRangeUser(pixelOverlayXaxisLeft, pixelOverlayXaxisRight, pulse_order);
     584
     585
    552586                PixelCsv.WritePointSet(
    553587                            pixel[pixelID],
    … …  
    618652
    619653//                pulses.cd(1);
    620                 Pulse maxMaxPulse("maxMaxPulse", pixel[pixelID]->hPixelMax[pulse_order], "Q", 1,pulse_order);
     654                Pulse maxMaxPulse("maxMaxPulse", pixel[pixelID]->hPixelMax[pulse_order], fitOption, pulseModell,pulse_order);
    621655                pulseFits[0] = &maxMaxPulse;
    622656
    623657//                pulses.cd(2);
    624                 Pulse maxMedianPulse("maxMedianPulse", pixel[pixelID]->hPixelMedian[pulse_order], "Q",1,pulse_order);
     658                Pulse maxMedianPulse("maxMedianPulse", pixel[pixelID]->hPixelMedian[pulse_order], fitOption,pulseModell,pulse_order);
    625659                pulseFits[1] = &maxMedianPulse;
    626660
    627661//                pulses.cd(3);
    628                 Pulse maxMeanPulse("maxMeanPulse", pixel[pixelID]->hPixelMean[pulse_order], "Q",1,pulse_order);
     662                Pulse maxMeanPulse("maxMeanPulse", pixel[pixelID]->hPixelMean[pulse_order], fitOption,pulseModell,pulse_order);
    629663                pulseFits[2] = &maxMeanPulse;
    630664
    631                 Pulse edgeMaxPulse("edgeMaxPulse", pixel[pixelID]->hPixelEdgeMax[pulse_order], "Q",1,pulse_order);
     665                Pulse edgeMaxPulse("edgeMaxPulse", pixel[pixelID]->hPixelEdgeMax[pulse_order], fitOption,pulseModell,pulse_order);
    632666                pulseFits[3] = &edgeMaxPulse;
    633667
    634                 Pulse edgeMedianPulse("edgeMedianPulse", pixel[pixelID]->hPixelEdgeMedian[pulse_order], "Q",1,pulse_order);
     668                Pulse edgeMedianPulse("edgeMedianPulse", pixel[pixelID]->hPixelEdgeMedian[pulse_order], fitOption,pulseModell,pulse_order);
    635669                pulseFits[4] = &edgeMedianPulse;
    636670
    637                 Pulse edgeMeanPulse("edgeMeanPulse", pixel[pixelID]->hPixelEdgeMean[pulse_order], "Q",1,pulse_order);
     671                Pulse edgeMeanPulse("edgeMeanPulse", pixel[pixelID]->hPixelEdgeMean[pulse_order], fitOption,pulseModell,pulse_order);
    638672                pulseFits[5] = &edgeMeanPulse;
    639673
    … …  
    754788//                            verbosityLevel
    755789//                            );
     790
     791                //Preparing Camera
     792                if (first_pass)
     793                {
     794                    if (verbosityLevel > 0)
     795                    {
     796                        cout << endl << "...preparing camera" << endl;
     797                    }
     798
     799                    wholeCamera = new PixelSum(
     800                                "AllPixel",
     801                                1440,
     802                                maxPulseOrder,
     803                                verbosityLevel,
     804                                stats,
     805                                "C",
     806                                pixel[pixelID]->mPixelOverlayXaxisLeft,
     807                                pixel[pixelID]->mPixelOverlayXaxisRight ,
     808                                pixel[pixelID]->mBSLMean ,
     809                                pixel[pixelID]->mGainMean ,
     810                                outputRootFile
     811                                );
     812
     813                    first_pass = false;
     814                }
     815
    756816                //-------------------------------------
    757817                // Fill Histogramms of Camera
    … …  
    897957                    "Maximum",
    898958                    pulse_order,
    899                     10,
     959                    bootstrapIt,
    900960                    verbosityLevel
    901961                    );
    … …  
    916976                    "Edge",
    917977                    pulse_order,
    918                     10,
     978                    bootstrapIt,
    919979                    verbosityLevel
    920980                    );
    … …  
    9591019
    9601020        pulses.cd(1);
    961         Pulse maxCamMaxPulse("maxMaxPulse", wholeCamera->hPixelMax[pulse_order], "Q", 1);
     1021        Pulse maxCamMaxPulse("maxMaxPulse", wholeCamera->hPixelMax[pulse_order], "Q", pulseModell);
    9621022        allPulseFits[0] = &maxCamMaxPulse;
    9631023
    9641024        pulses.cd(2);
    965         Pulse maxCamMedianPulse("maxMedianPulse", wholeCamera->hPixelMedian[pulse_order], "Q", 1);
     1025        Pulse maxCamMedianPulse("maxMedianPulse", wholeCamera->hPixelMedian[pulse_order], "Q", pulseModell);
    9661026        allPulseFits[1] = &maxCamMedianPulse;
    9671027
    9681028        pulses.cd(3);
    969         Pulse maxCamMeanPulse("maxMeanPulse", wholeCamera->hPixelMean[pulse_order],"Q", 1);
     1029        Pulse maxCamMeanPulse("maxMeanPulse", wholeCamera->hPixelMean[pulse_order],"Q", pulseModell);
    9701030        allPulseFits[2] = &maxCamMeanPulse;
    9711031
    972         Pulse edgeCamMaxPulse("edgeMaxPulse", wholeCamera->hPixelEdgeMax[pulse_order], "Q", 1);
     1032        Pulse edgeCamMaxPulse("edgeMaxPulse", wholeCamera->hPixelEdgeMax[pulse_order], "Q", pulseModell);
    9731033        allPulseFits[3] = &edgeCamMaxPulse;
    9741034
    975         Pulse edgeCamMedianPulse("edgeMedianPulse", wholeCamera->hPixelEdgeMedian[pulse_order],"Q", 1);
     1035        Pulse edgeCamMedianPulse("edgeMedianPulse", wholeCamera->hPixelEdgeMedian[pulse_order],"Q", pulseModell);
    9761036        allPulseFits[4] = &edgeCamMedianPulse;
    9771037
    978         Pulse edgeCamMeanPulse("edgeMeanPulse", wholeCamera->hPixelEdgeMean[pulse_order], "Q", 1);
     1038        Pulse edgeCamMeanPulse("edgeMeanPulse", wholeCamera->hPixelEdgeMean[pulse_order], "Q", pulseModell);
    9791039        allPulseFits[5] = &edgeCamMeanPulse;
    9801040
    … …  
    10611121                    saveResults
    10621122                    );
     1123
     1124        SaveList(
     1125                outFile,
     1126                    "",
     1127                histoList,
     1128                saveResults,
     1129                verbosityLevel
     1130                );
     1131
     1132
    10631133    }
    10641134
    … …  
    10901160    delete inputRootFile;
    10911161    delete outputRootFile;
     1162    delete histoList;
    10921163
    10931164    return( 0 );
Note: See TracChangeset for help on using the changeset viewer.