ROOT logo

gSystem->Load("libSTAT.so");
TStatToolkit stat;
Int_t npoints;
Double_t chi2;
TVectorD vec;
TMatrixD mat;


TObjArray * array = AliTPCCalibViewerGUI::ShowGUI("laserTree.root");
AliTPCCalibViewerGUI * viewer = (AliTPCCalibViewerGUI*)array->At(0);
TTree * tree = viewer->GetViewer()->GetTree();
TFile fp("/data/calib/CalibTreePulser_run33834_Cside.root");
tree->AddFriend(treePulser,"P.");


tree->SetAlias("dt","(sector%36>30)*2+(sector<36)*0.3");
tree->SetAlias("T","T0_100_220.fElements-P..StandardTime0.fElements");
tree->SetAlias("Tm","T0_100_220_Median.fElements");
tree->SetAlias("Q","Q_100_220.fElements");
tree->SetAlias("Qm","Q_100_220_Median.fElements");

tree->SetAlias("Qcut","abs(Q/Qm-2)<1.4&&Q>6&&Q<200");
tree->SetAlias("Tcut","abs(T-Tm)<2");



TString strSector="";
{
  for (Int_t isec=54;isec<71;isec+=1){
    if (isec!=64) {
      strSector+="(sector==";
      strSector+=isec;
      strSector+=")++";
      strSector+="(lx.fElements-195.)*(sector==";
      strSector+=isec;
      strSector+=")++";
      strSector+="((lx.fElements-195)^2)*(sector==";
      strSector+=isec;
      strSector+=")++";
    }
  }
}


TCut cutA="Tcut&&Qcut&&sector%36>17";

TString *strFitG =stat.FitPlane(tree,"T+dt","gx.fElements++gy.fElements",cutA,chi2,npoints,vec,mat);

TString *strFitGL =stat.FitPlane(tree,"T+dt","lx.fElements++ly.fElements++gx.fElements++gy.fElements",cutA,chi2,npoints,vec,mat);

TString *strFitGL2 =stat.FitPlane(tree,"T+dt","lx.fElements++ly.fElements++gx.fElements++gy.fElements++lx.fElements^2++ly.fElements^2",cutA,chi2,npoints,vec,mat);

TString *strFitGLA =stat.FitPlane(tree,"T+dt",strSector+"lx.fElements++ly.fElements++gx.fElements++gy.fElements++lx.fElements^2++ly.fElements^2",cutA,chi2,npoints,vec,mat);


tree->SetAlias("tfitG",strFitG->Data())
tree->SetAlias("tfitGL",strFitGL->Data())
tree->SetAlias("tfitGL2",strFitGL2->Data())
tree->SetAlias("tfitGLA",strFitGLA->Data())


 AnalyzeLaser.C:1
 AnalyzeLaser.C:2
 AnalyzeLaser.C:3
 AnalyzeLaser.C:4
 AnalyzeLaser.C:5
 AnalyzeLaser.C:6
 AnalyzeLaser.C:7
 AnalyzeLaser.C:8
 AnalyzeLaser.C:9
 AnalyzeLaser.C:10
 AnalyzeLaser.C:11
 AnalyzeLaser.C:12
 AnalyzeLaser.C:13
 AnalyzeLaser.C:14
 AnalyzeLaser.C:15
 AnalyzeLaser.C:16
 AnalyzeLaser.C:17
 AnalyzeLaser.C:18
 AnalyzeLaser.C:19
 AnalyzeLaser.C:20
 AnalyzeLaser.C:21
 AnalyzeLaser.C:22
 AnalyzeLaser.C:23
 AnalyzeLaser.C:24
 AnalyzeLaser.C:25
 AnalyzeLaser.C:26
 AnalyzeLaser.C:27
 AnalyzeLaser.C:28
 AnalyzeLaser.C:29
 AnalyzeLaser.C:30
 AnalyzeLaser.C:31
 AnalyzeLaser.C:32
 AnalyzeLaser.C:33
 AnalyzeLaser.C:34
 AnalyzeLaser.C:35
 AnalyzeLaser.C:36
 AnalyzeLaser.C:37
 AnalyzeLaser.C:38
 AnalyzeLaser.C:39
 AnalyzeLaser.C:40
 AnalyzeLaser.C:41
 AnalyzeLaser.C:42
 AnalyzeLaser.C:43
 AnalyzeLaser.C:44
 AnalyzeLaser.C:45
 AnalyzeLaser.C:46
 AnalyzeLaser.C:47
 AnalyzeLaser.C:48
 AnalyzeLaser.C:49
 AnalyzeLaser.C:50
 AnalyzeLaser.C:51
 AnalyzeLaser.C:52
 AnalyzeLaser.C:53
 AnalyzeLaser.C:54
 AnalyzeLaser.C:55
 AnalyzeLaser.C:56
 AnalyzeLaser.C:57
 AnalyzeLaser.C:58
 AnalyzeLaser.C:59
 AnalyzeLaser.C:60
 AnalyzeLaser.C:61
 AnalyzeLaser.C:62
 AnalyzeLaser.C:63