Ignore:
Timestamp:
05/18/12 23:35:09 (12 years ago)
Author:
Jens Buss
Message:
whitespace changes
File:
1 edited

Legend:

Unmodified
Added
Removed
  • fact/tools/rootmacros/PulseTemplates/FCalcPulseTemplate.C

    r13780 r13790  
    101101{
    102102
    103 
    104     InputPath = SetHostsPaths(true, InputPath );
    105     OutPutPath = SetHostsPaths(true, OutPutPath );
    106 
     103    InputPath       = SetHostsPaths(true, InputPath );
     104    OutPutPath      = SetHostsPaths(true, OutPutPath );
    107105
    108106//----------------------------------------------------------------------------
     
    114112    }
    115113
    116     TFile * inputRootFile   = OpenRootFile( InputPath, InRootFileName, verbosityLevel );
    117     TFile * outputRootFile  = OpenRootFile( OutPutPath, OutputRootFileName, verbosityLevel );
    118 //----------------------------------------------------------------------------
    119 //      global variable Settings
     114    TFile * inputRootFile
     115            = OpenRootFile( InputPath, InRootFileName, verbosityLevel );
     116
     117    TFile * outputRootFile
     118            = OpenRootFile( OutPutPath, OutputRootFileName, verbosityLevel );
     119//----------------------------------------------------------------------------
     120//      Define operation range
    120121//----------------------------------------------------------------------------
    121122    if ( npixel == -1 )
     
    129130        pixelSetSize = firstpixel +npixel;
    130131    }
     132
    131133//    float GainMean  = GainMean;  // this has to be extracted from root files
    132134//    float BSLMean   = BSLMean;    // this has to be extracted from root files
     
    155157
    156158    //Canvas Pad numbering
    157     int PixelCanvasFrameNrs[4] = {
     159    int PixelCanvasFrameNrs[4]      = {
    158160        1,  // Top left
    159161        2,  // Top right
     
    163165
    164166    //Canvas Pad numbering
    165     int AllpixelCanvasFrameNrs[4] = {
     167    int AllpixelCanvasFrameNrs[4]   = {
    166168        1,  // Top left
    167169        2,  // Top right
     
    174176
    175177        //Canvases
    176         cgpPixelPulses      = new TCanvas*[maxPulseOrder];
    177         cgpDistributions    = new TCanvas*[maxPulseOrder];
     178        cgpPixelPulses              = new TCanvas*[maxPulseOrder];
     179        cgpDistributions            = new TCanvas*[maxPulseOrder];
    178180
    179181        //TCanvas*    gpcDevelopment = NULL;
    180         TString cName   = "";
    181         TString cTitle  = "";
     182        TString cName               = "";
     183        TString cTitle              = "";
    182184
    183185        //naming of pulse canvases
     
    189191        {
    190192            cName   ="cgpDistributions";
    191             cName   += pulse_order;
    192             cTitle   ="Distributions of Pulses with Order of: ";
    193             cTitle   += pulse_order;
    194             cgpDistributions[pulse_order] = new TCanvas(cName,cTitle, 720,pulse_order*20,720,720);
     193            cName  += pulse_order;
     194
     195            cTitle  ="Distributions of Pulses with Order of: ";
     196            cTitle += pulse_order;
     197
     198            cgpDistributions[pulse_order]
     199                    = new TCanvas(cName,cTitle, 720,pulse_order*20,720,720);
    195200            cgpDistributions[pulse_order]->Divide(2, 2);
     201
    196202            cName   ="cgpPixelPulses";
    197             cName   += pulse_order;
    198             cTitle   ="Overlays of Pulses with Order of: ";
    199             cTitle   += pulse_order;
    200             cgpPixelPulses[pulse_order] = new TCanvas(cName,cTitle, 0,pulse_order*20,720,720);
     203            cName  += pulse_order;
     204
     205            cTitle  ="Overlays of Pulses with Order of: ";
     206            cTitle  += pulse_order;
     207
     208            cgpPixelPulses[pulse_order]
     209                     = new TCanvas(cName,cTitle, 0,pulse_order*20,720,720);
    201210            cgpPixelPulses[pulse_order]->Divide(2, 2);
    202211        }
    203 
    204 
    205         // Create (pointer to) Canvases, which are used in every run,
    206         // also in 'non-debug' runs
    207                 // Canvases only need if spike Debug, but I want to deklare
    208                 // the pointers anyway ...
    209 
    210 //        if (testmode)
    211 //        {
    212 //        //additional Test histograms
    213 //        cgpTestHistos = new TCanvas( "cgpTestHistos", "Test Histograms", 360, 420, 360, 360 );
    214 //        cgpTestHistos->Divide(2,0);
    215 //        }
    216212    }
    217213
     
    226222
    227223//----------------------------------------------------------------------------
    228 // Initialize Pixel
     224// Initialize Pixels
    229225//----------------------------------------------------------------------------
    230226    if (verbosityLevel > 0)
     
    233229    }
    234230
    235     Pixel** pixel = new Pixel*[NPIX];
     231    Pixel**     pixel      = new Pixel*[NPIX];
    236232
    237233    for (int i = 0 ; i < NPIX; i++)
     
    240236    }
    241237
    242     PixelSum* wholeCamera = NULL;
    243 
    244     bool first_pass = true;
     238    PixelSum*   wholeCamera = NULL;
     239
     240    bool        first_pass  = true;
     241
    245242//-------------------------------------
    246243// Loop over Pixel Sets
     
    252249        if (verbosityLevel >= 0)
    253250        {
    254             cout << "------------------------------------------------" << endl
     251            cout << "------------------------------------------------"
     252                 << endl
    255253                 << "...processing Pixel: "
    256254                 << firstPixelOfSet
    257255                 << " to Pixel: "
    258                  << firstPixelOfSet+pixelSetSize-1 << endl;
     256                 << firstPixelOfSet+pixelSetSize-1
     257                 << endl;
    259258        }
    260259
     
    263262        //--------------------------------------------------------------------
    264263        for ( int pixelID = firstPixelOfSet;
    265                 pixelID < firstPixelOfSet + pixelSetSize
    266                 && pixelID < firstpixel + npixel;
    267                 pixelID++ )
     264              pixelID < firstPixelOfSet + pixelSetSize
     265              && pixelID < firstpixel + npixel;
     266              pixelID++ )
    268267        {
    269268
    270269            if (verbosityLevel > 1)
    271270            {
    272              cout << "-------------------------------------" << endl
    273                   << "...processing Set from Pixel "
    274                   << firstPixelOfSet
    275                   << " to Pixel "
    276                   << firstPixelOfSet+pixelSetSize-1
    277                   << " Pixel: " << pixelID
    278                   << "/" << firstpixel + npixel -1 << endl;
     271                cout << "-------------------------------------"
     272                     << endl
     273                     << "...processing Set from Pixel "
     274                     << firstPixelOfSet
     275                     << " to Pixel "
     276                     << firstPixelOfSet+pixelSetSize-1
     277                     << " Pixel: " << pixelID
     278                     << "/" << firstpixel + npixel -1
     279                     << endl;
    279280            }
    280281
     
    286287                cout << "...creating pixel: " << pixelID << endl;
    287288            }
     289
    288290           pixel[pixelID] = new Pixel(
    289291                       pixelID,
     
    308310                    cout << endl << "...preparing camera" << endl;
    309311                }
     312
    310313                wholeCamera = new PixelSum(
    311314                            "AllPixel",
     
    321324                            outputRootFile
    322325                            );
     326
    323327                first_pass = false;
    324328            }
     
    333337                if (verbosityLevel > 2)
    334338                {
    335                     cout << "-------------------------------------" << endl
     339                    cout << "-------------------------------------"
     340                         << endl
    336341                         << "...processing Set from Pixel "
    337342                         << firstPixelOfSet
     
    340345                         << " Pixel: " << pixelID
    341346                         << "/" << firstpixel + npixel -1
    342                          << " Pulse-Order: " << pulse_order << endl;
     347                         << " Pulse-Order: " << pulse_order
     348                         << endl;
    343349                }
    344350
    345351                // Calculate Max Prop. Value of each slice
    346352                //-------------------------------------
     353
     354                //from Maximum Overlay
    347355                if (verbosityLevel > 2)
    348356                {
    349                     cout << "...extracting templates from Maximum Overlay " << endl;
     357                    cout << "...extracting templates from Maximum Overlay "
     358                         << endl;
    350359                }
    351                 //from Maximum Overlay
     360
    352361                ExtractPulseTemplate(
    353362                            pixel[pixelID],
     
    357366                            );
    358367
     368                //from Edge Overlay
    359369                if (verbosityLevel > 2)
    360370                {
    361                     cout << "...extracting templates from Edge Overlay " << endl;
     371                    cout << "...extracting templates from Edge Overlay "
     372                         << endl;
    362373                }
    363                 //from Edge Overlay
    364374                ExtractPulseTemplate(
    365375                            pixel[pixelID],
     
    384394                            verbosityLevel
    385395                            );
     396
    386397                if (ProduceGraphic)
    387398                {
    388                 pixel[pixelID]->DrawTemplateHistograms(
    389                             cgpPixelPulses,
    390                             PixelCanvasFrameNrs
    391                             );
    392 
    393                 pixel[pixelID]->DrawEdgeTemplateHistograms(
    394                             cgpPixelPulses,
    395                             PixelCanvasFrameNrs
    396                             );
     399                    pixel[pixelID]->DrawTemplateHistograms(
     400                                cgpPixelPulses,
     401                                PixelCanvasFrameNrs
     402                                );
     403
     404                    pixel[pixelID]->DrawEdgeTemplateHistograms(
     405                                cgpPixelPulses,
     406                                PixelCanvasFrameNrs
     407                                );
    397408                }
    398409                //-------------------------------------
    399410                // Fill Histogramms of Camera
    400411                //-------------------------------------
    401                 cout << "1" << endl;
     412
    402413                wholeCamera->hMaxOverlay[pulse_order]->Add(
    403414                            pixel[pixelID]->hMaxOverlay[pulse_order]
    404415                            );
    405                 cout << wholeCamera->hMaxOverlay[pulse_order]->GetName() << endl;
    406 cout << "2" << endl;
     416
    407417                wholeCamera->hMaxProfile[pulse_order]->Add(
    408418                            pixel[pixelID]->hMaxProfile[pulse_order]
    409419                            );
    410 cout << "3" << endl;
     420
    411421                wholeCamera->hEdgeOverlay[pulse_order]->Add(
    412422                            pixel[pixelID]->hEdgeOverlay[pulse_order]
    413423                            );
    414 cout << "4" << endl;
     424
    415425                wholeCamera->hEdgeProfile[pulse_order]->Add(
    416426                            pixel[pixelID]->hEdgeProfile[pulse_order]
    417427                            );
    418 cout << "5" << endl;
     428
    419429
    420430                if (verbosityLevel > 2)
    421431                {
    422                 cout << endl << "...End of pulseorder "
    423                     << pulse_order
    424                     << endl;
     432                cout << endl
     433                     << "...End of pulseorder "
     434                     << pulse_order
     435                     << endl;
    425436                }
    426437
     
    430441            if (ProduceGraphic)
    431442            {
    432             UpdateCanvases(
    433                         verbosityLevel,
    434                         MAX_PULS_ORDER,
    435                         false
    436                         );
     443                UpdateCanvases(
     444                            verbosityLevel,
     445                            MAX_PULS_ORDER,
     446                            false
     447                            );
    437448            }
    438449
     
    454465
    455466            //deleteCurrent Pixel from Heap
    456             delete pixel[pixelID];
    457             pixel[pixelID] = NULL;
     467            delete  pixel[pixelID];
     468            pixel[pixelID]  = NULL;
     469
    458470            if (verbosityLevel > 2)
    459471            {
    460             cout << endl << "...End of Pixel"
    461                 << endl << "------------------------------------------------"
    462                 << endl;
     472            cout << endl
     473                 << "...End of Pixel"
     474                 << endl
     475                 << "------------------------------------------------"
     476                 << endl;
    463477            }
    464478        }
     
    467481        if (verbosityLevel > 1)
    468482        {
    469         cout << endl << "...End of Loop over all Pixels of set"
    470             << endl << "------------------------------------------------"
    471             << endl;
     483        cout << endl
     484             << "...End of Loop over all Pixels of set"
     485             << endl
     486             << "------------------------------------------------"
     487             << endl;
    472488        }
    473489    }
     
    476492    if (verbosityLevel > 0)
    477493    {
    478     cout << endl << "...End of Loop over all Pixelsets"
    479         << endl << "------------------------------------------------"
    480         << endl;
     494    cout << endl
     495         << "...End of Loop over all Pixelsets"
     496         << endl
     497         << "------------------------------------------------"
     498         << endl;
    481499    }
    482500
    483501    delete[] pixel;
    484502    pixel = NULL;
     503
    485504//-------------------------------------
    486505// Draw All Pixel Histograms
     
    493512        if (verbosityLevel > 2)
    494513        {
    495             cout << "-------------------------------------" << endl
    496                  << "...processing Pulse-Order: " << pulse_order;
     514            cout << "-------------------------------------"
     515                 << endl
     516                 << "...processing Pulse-Order: "
     517                 << pulse_order;
    497518        }
    498519
    499 //         Calculate Max Prop. Value of each slice
    500 //        -------------------------------------
     520// Calculate Max Prop. Value of each slice
     521//-------------------------------------
    501522
    502523        //from Maximum Overlay
     
    507528                    verbosityLevel
    508529                    );
     530        WritePixelTemplateToCsv(
     531                    wholeCamera,
     532                    OutPutPath,
     533                    "Maximum",
     534                    pulse_order,
     535                    verbosityLevel
     536                    );
    509537
    510538        //from Edge Overlay
     
    515543                    verbosityLevel
    516544                    );
    517 
    518         WritePixelTemplateToCsv(
    519                     wholeCamera,
    520                     OutPutPath,
    521                     "Maximum",
    522                     pulse_order,
    523                     verbosityLevel
    524                     );
    525 
    526545        WritePixelTemplateToCsv(
    527546                    wholeCamera,
     
    543562    } //EOF: Draw All Pixel Histograms
    544563
     564
    545565//-------------------------------------
    546566// Save All Pixel Histograms
     
    560580//                verbosityLevel
    561581//                );
    562 
    563 
    564582
    565583//    if (ProduceGraphic)
     
    572590//    }
    573591
     592//-------------------------------------
     593// Delete Objects on Heap
     594//-------------------------------------
    574595    delete wholeCamera;
    575596    if (ProduceGraphic)
     
    579600    delete inputRootFile;
    580601    delete outputRootFile;
     602
    581603    return( 0 );
    582604}
     
    593615void
    594616DeletePixelCanvases(
    595         int maxPulseOrder,
    596         int verbosityLevel
     617        int         maxPulseOrder,
     618        int         verbosityLevel
    597619        )
    598620{
    599621    if (verbosityLevel > 2)
    600622    {
    601         cout << endl << "...delete pixel Canvases" << endl;
    602     }
    603     for (int pulse_order = 0; pulse_order < maxPulseOrder; pulse_order++ )
     623        cout << endl
     624             << "...delete pixel Canvases"
     625             << endl;
     626    }
     627
     628    for ( int pulse_order = 0; pulse_order < maxPulseOrder; pulse_order++ )
    604629    {
    605630        delete cgpPixelPulses[pulse_order];
    606         cgpPixelPulses[pulse_order] = NULL;
     631        cgpPixelPulses[pulse_order]     = NULL;
     632
    607633        delete cgpDistributions[pulse_order];
    608         cgpDistributions[pulse_order] = NULL;
    609     }
     634        cgpDistributions[pulse_order]   = NULL;
     635    }
     636
    610637    delete[] cgpPixelPulses;
    611     cgpPixelPulses = NULL;
     638    cgpPixelPulses      = NULL;
     639
    612640    delete[] cgpDistributions;
    613     cgpDistributions = NULL;
     641    cgpDistributions    = NULL;
    614642}
    615643
    616644void
    617645UpdateCanvases(
    618         int verbosityLevel,
    619         int max_pulse_order,
    620         bool testmode
     646        int         verbosityLevel,
     647        int         max_pulse_order,
     648        bool        testmode
    621649        )
    622650{
    623     if (verbosityLevel > 3) cout << endl << "...updating canvases" ;
     651    if (verbosityLevel > 3)
     652    {
     653        cout << endl << "...updating canvases" << endl;
     654    }
    624655    for (int pulse_order = 0; pulse_order < max_pulse_order; pulse_order++)
    625656    {
    626657        cgpPixelPulses[pulse_order]->Modified();
    627658        cgpPixelPulses[pulse_order]->Update();
     659
    628660        cgpDistributions[pulse_order]->Modified();
    629661        cgpDistributions[pulse_order]->Update();
     
    641673int main()
    642674{
    643 
    644 FCalcPulseTemplate();
    645 return 0;
    646 
     675    FCalcPulseTemplate();
     676    return 0;
    647677}
Note: See TracChangeset for help on using the changeset viewer.