runAnalysisBC.C

Go to the documentation of this file.

// --- ROOT system ---
#include <TStopwatch.h>
#include <TROOT.h>

// --- ANALYSIS system ---
#include "BadChannelAna.h" //include when compile

void runAnalysisBC(Int_t nversion = -1, TString period = "LHC15n", TString train = "Train_603", TString trigger= "INT7", Int_t runNum= 245683, TString externalFile= "",TString listName="runList.txt",TString workDir=".")
{
  gROOT->ProcessLine("gErrorIgnoreLevel = kWarning;"); //..to supress a lot of standard output
    std::vector<Int_t> badcellsManual;

    //..15o block 1
    //badcellsManual.insert(badcellsManual.end(),{14655,14622,14640,14728,14726,14593,14599,14600,14645,14646,14759,14776});

  //..If nversion=-1
  //..this is detected as a run-by-run analysis
  if(nversion == -1)nversion=runNum; //..If you do the analysis run by run - this might be helpful

  TStopwatch* watch = new TStopwatch();
  watch->Start();

  BadChannelAna* Analysis;
  Analysis=new BadChannelAna(period,train,trigger,runNum,nversion,workDir,listName);

  //. . . . . . . . . . . . . . . . . . . . . . . .
  //. . Settings
  //. . . . . . . . . . . . . . . . . . . . . . . .
  Analysis->SetExternalMergedFile(externalFile);
    Analysis->AddManualMasking(badcellsManual);
  //Analysis->AddMaskSM(7);                     //..switch off entire SMs
  //Analysis->AddMaskSM(9);                     //..switch off entire SMs
  //Analysis->AddMaskSM(11);                    //..switch off entire SMs
  Analysis->SetLowerBound(1);                 //..If the Emin of the energy range (Emin-Emax) is higher than X GeV then dont apply a lower cut on the distribution
    //Analysis->SetStartEndCell(0,12288);       //..only EMCal
    //Analysis->SetStartEndCell(12288,17664);   //..only DCal
    //Analysis->SetQAChecks(1);                 //..1= Perform QA checks - takes a long time! Saves all good cells for cross check to pdf
  //Analysis->SetPrintOutput(1);              //..1= prints more information about excluded cells
  //Analysis->SetExternalBadMap("Version1OADB/LHC16s_INT7_Histograms_V1.root");

  //. . . . . . . . . . . . . . . . . . . . . . . .
  //. . Add different period analyses
  //. . . . . . . . . . . . . . . . . . . . . . . .
  //..the range of sigmas should be selected such
  //..that one does not cut into the natural fluctuation over the modules
  Analysis->AddPeriodAnalysis(2, 5.5,0.1,0.3);  // hits in cell in range Emin Emax
  Analysis->AddPeriodAnalysis(1, 4.0,0.1,0.3);  // energy/hit in range Emin Emax
  Analysis->AddPeriodAnalysis(2, 5.5,0.2,0.5);  // hits in cell range Emin Emax
  Analysis->AddPeriodAnalysis(1, 4.5,0.2,0.5);  // energy/hit in range Emin Emax
  Analysis->AddPeriodAnalysis(2, 5.5,0.5,1.0);  // hits in cell range Emin Emax
  Analysis->AddPeriodAnalysis(1, 4.5,0.5,1.0);  // energy/hit in range Emin Emax
  Analysis->AddPeriodAnalysis(2, 5.5,1.0,4.0);  // hits in cell range Emin Emax
    Analysis->AddPeriodAnalysis(1, 4.5,1.0,4.0);  // mean energy in range Emin Emax
  Analysis->AddPeriodAnalysis(2, 5.5,1.0,10.0); // hits in cell in range Emin Emax
  Analysis->AddPeriodAnalysis(1, 4.5,1.0,10.0); // energy/hit in range Emin Emax
  Analysis->AddPeriodAnalysis(2, 5.5,0.11,0.29);  // hits in cell in range Emin Emax


  //..special test for extra high energy fluctuations
  Analysis->AddPeriodAnalysis(1, 4.0,3.0,40.0); // energy/hit in cell in range Emin Emax
//  Analysis->AddPeriodAnalysis(2, 4.0,3.0,40.0); // hits in cell range Emin Emax
  //Analysis->AddPeriodAnalysis(1, 4.5,3.0,5.0);  // mean energy in range Emin Emax - cliff
  //Analysis->AddPeriodAnalysis(2, 5.5,3.0,5.0);  // hits in cell range Emin Emax   - cliff
  //Analysis->AddPeriodAnalysis(1, 4.5,5,20);     // energy/hit in range Emin Emax
  //Analysis->AddPeriodAnalysis(2, 5.5,5,20);     // hits in range Emin Emax

  //*test time stuff*/  Analysis->AddPeriodAnalysis(3, 6,-20,+20);// energy/hit in range Emin Emax

  //. . . . . . . . . . . . . . . . . . . . . . . .
  //. . Start the bad channel analysis
  //. . . . . . . . . . . . . . . . . . . . . . . .
  Bool_t mergeOnly=0;//.. =1 do only merge and filter
  Analysis->Run(mergeOnly);

  watch->Stop();
  cout<<"Finished BC analysis "<<watch->RealTime()/60<<" min"<<endl;
  watch->Continue();
}