#include <TObjArray.h>
#include <TClass.h>
#include <TKey.h>
-//#include "AliRecoParam.h"
+#include <TFile.h>
+#include <TFileMerger.h>
+#include <TAlienFile.h>
+//#include <TExec.h>
+#include <TSystem.h>
+#include <TGrid.h>
+#include <TGridResult.h>
+#include <Riostream.h>
+#include <TObjArray.h>
+#include <TClass.h>
+#include <TH1F.h>
+#include <TH2D.h>
+#include <TH2F.h>
+#include <TF1.h>
+#include <TStyle.h>
+#include <TCanvas.h>
+#include <TLegend.h>
+#include <TMath.h>
+#include "AliRecoParam.h"
#endif
-enum rawexpert_t{pattern,patternnorm,layer3norm,layer4norm,rawtot};
+enum rawexpert_t{pattern,patternnorm,layer3norm,layer4norm,layer3occ,layer4occ,rawtot};
enum nonrawexpert_t{layer3,layer4,nonrawtot};
-enum rpexpert_t{rlocaldistro,rlayer3norm,rlayer4norm,rmodpattern,rmodpatternnorm,rphilayer3,rphilayer4,rrlayer3,rrlayer4,rtot};
+enum rpexpert_t{rlocaldistro,rlayer3norm,rlayer4norm,rmodpattern,rmodpatternnorm,rphilayer3,rphilayer4,rrlayer3,rrlayer4,occlayer3,occlayer4,recraw2dlayer3,recraw2dlayer4,recrawlayer3,recrawlayer4,dedxlayer3,dedxlayer4,rtot};
enum rnonexpert_t{rphizlayer3,rphizlayer4,rglobalrz,rglobalxy,rchargelayer3,rchargelayer4,rtimelayer3,rtimelayer4,rlayer3,rlayer4,nrtot};
-enum canvsname_t{canvname0,canvname1,canvname2,canvname3,canvname4,canvname5,canvname6,canvname7,canvname8,canvname9,canvname10,canvname11,canvname12};
+enum canvsname_t{canvname0,canvname1,canvname2,canvname3,canvname4,canvname5,canvname6,canvname7,canvname8,canvname9,canvname10,canvname11,canvname12,canvname13,canvname14,canvname15,canvname16,canvname17,canvname18,canvname19,canvname20,canvname21,canvname22,canvname23};//
-void PlotQASDD(Char_t fileName[100]="File.QA.111333.2010.LHC09b.pass2.root",Char_t eventspecie[25]="LowMultiplicity")
+
+
+void PlotQASDD(Char_t fileName[100]="File.QA.111333.2010.LHC09b.pass2.root",Char_t eventspecie[25]="LowMultiplicity",Bool_t kDoEps=kFALSE)
{
- const TString Rawsexpertname[]={"SDDModPattern","SDDModPatternNORM","SDDphizL3NORM","SDDphizL4NORM"};//4
- const TString Rawsnonexpertname[]={"SDDphizL3","SDDphizL4"};//2
+ const TString Rawsexpertname[]={"SDDModPattern","SDDModPatternNORM","SDDphizL3","SDDphizL4","SDDL3_RelativeOccupancy","SDDL4_RelativeOccupancy"};//6
+ const TString Rawsnonexpertname[]={"SDDphizL3NORM","SDDphizL4NORM"};//2
+
+ const TString RecPointsexpertname[]={"SDDLocalCoordDistrib","SDDModPatternL3RP","SDDModPatternL4RP","SDDModPatternRP","SDDModPatternRPNORM","SDDphidistrib_Layer3","SDDphidistrib_Layer4","SDDrdistrib_Layer3","SDDrdistrib_Layer4","SDDL3_RelativeOccupancy","SDDL4_RelativeOccupancy","SDDL3_Rec2Raw_2D","SDDL4_Rec2Raw_2D","SDDL3_Rec2Raw","SDDL4_Rec2Raw","SDDL3_dedx","SDDL4_dedx"};//17
+ const TString RecPointsnonexpertname[]={"SDDGlobalCoordDistribL3PHIZ","SDDGlobalCoordDistribL4PHIZ","SDDGlobalCoordDistribRZ","SDDGlobalCoordDistribYX","SDDLay3TotCh","SDDLay4TotCh","SDDdrifttime_Layer3","SDDdrifttime_Layer4","SDDModPatternL3RPNORM","SDDModPatternL4RPNORM"};//10
- const TString RecPointsexpertname[]={"SDDLocalCoordDistrib","SDDModPatternL3RPNORM","SDDModPatternL4RPNORM","SDDModPatternRP","SDDModPatternRPNORM","SDDphidistrib_Layer3","SDDphidistrib_Layer4","SDDrdistrib_Layer3","SDDrdistrib_Layer4"};//9
- const TString RecPointsnonexpertname[]={"SDDGlobalCoordDistribL3PHIZ","SDDGlobalCoordDistribL4PHIZ","SDDGlobalCoordDistribRZ","SDDGlobalCoordDistribYX","SDDLay3TotCh","SDDLay4TotCh","SDDdrifttime_Layer3","SDDdrifttime_Layer4","SDDModPatternL3RP","SDDModPatternL4RP"};//1
+ const TString canvassavedname[]={"RawLayers","Rawpatterns","RawLayersNORM","RawLayersRelOcc","RawLayersRelOccBoth","RecPointsphiz","RecPointsGlobalDistributions","RecPointsCharge","RecPointsDriftTime","RecPointsLocalDistribution","RecPointsLayers","RecPointsLayersNORM","RecPointsLayersPatterns","RecPointsPhiDistribution","RecPointsRDistribution","RecPointsChargeBothLayers","RecPointsDriftTimeBothLayers","RecPointsRelLayersOcc","RecPointsLayersRelOccBoth","RecPointsLayersRecRaw2D","RecPointsLayersRecRaw","RecPointsLayersRecRaw","RecPointsLayersRecRawBoth","RecPointsLayersdedx","RecPointsLayersdedxBoth"};//15
+
+ Char_t trendfile[25]="FileQAtrend.root";
- const TString canvassavedname[]={"RawLayers","Rawpatterns","RawLayersNORM","RecPointsphiz","RecPointsGlobalDistributions","RecPointsCharge","RecPointsDriftTime","RecPointsLocalDistribution","RecPointsLayers","RecPointsLayersNORM","RecPointsLayersPatterns","RecPointsPhiDistribution","RecPointsRDistribution"};//13
gStyle->SetPalette(1);
- gStyle->SetOptStat(0);
- TFile mergedfile(fileName);
- if(mergedfile.GetKey("ITS")==0x0){
- printf("In this run ITS QA has not been executed.\n\nExit macro \n\n");
+ // gStyle->SetOptStat("ne");
+ if(gSystem->Exec(Form("ls %s >/dev/null 2>&1", fileName))!=0)
+ {
+ printf(" No file created --- Exit");
return;
}
+ TFile mergedfile(fileName);
+ if(mergedfile.GetKey("ITS")==0x0){
+ printf("In this run ITS QA has not been executed.\n\nExit macro \n\n");
+ return;
+ }
Char_t filepath[100];
Char_t namecanvas[50];
TDirectory *directory=NULL;
TH2D *historaw=NULL;
TH1F *historaw2=NULL;
TH2F *histodraw=NULL;
-
+
+ TH1F *histo2save[9];
+ for(Int_t i=0;i<9;i++){histo2save[i]=NULL;}
+
+ TF1 *f1 = new TF1("f1","1",-100000,100000); // this is the function defined to multiply histograms
+ Char_t layer[10];
+ float fmax=0.;
+ float fmaxold=0.;
+ float fmaxmargin;
+
+ Int_t cannum=0;
+ Int_t histo2savenumber[9]={4,5,6,7,7,8,3,5,1};
TCanvas *canvas1 = new TCanvas("canvas1","SDD QA Plot",1000,600);
-// canvas1->Divide(2,1);
-// TCanvas *canvas3 = new TCanvas("canvas3","RawLayersNORM",1000,600);
-// canvas3->Divide(2,1);
-// TCanvas *canvas2 = new TCanvas("canvas2","Rawpatterns",1000,600);
-// canvas2->Divide(1,2);
-
-// TCanvas *canvas4 = new TCanvas("canvas4","RecPointsphiz",1000,600);
-// canvas4->Divide(2,1);
-// TCanvas *canvas5 = new TCanvas("canvas5","RecPointsGlobalDistributions",1000,600);
-// canvas5->Divide(2,1);
-// TCanvas *canvas6 = new TCanvas("canvas6","RecPointsCharge",1000,600);
-// canvas6->Divide(1,2);
-// TCanvas *canvas7 = new TCanvas("canvas7","RecPointsDriftTime",1000,600);
-// canvas7->Divide(1,2);
-
-// TCanvas *canvas8 = new TCanvas("canvas8","RecPointsLocalDistribution",1000,600);
-// TCanvas *canvas9 = new TCanvas("canvas9","RecPointsLayers",1000,600);
-// canvas9->Divide(2,1);
+ // TLegend *legend=new TLegend(0.83,0.8,0.97,0.7);
+ TLegend *legend=new TLegend(0.81,0.895,0.95,0.995);
+ Bool_t updatecanvas[2];
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
-// TCanvas *canvas10 = new TCanvas("canvas10","RecPointsLayersNORM",1000,600);
-// canvas10->Divide(2,1);
-// TCanvas *canvas11 = new TCanvas("canvas11","RecPointsLayersPatterns",1000,600);
-// canvas11->Divide(1,2);
-// TCanvas *canvas12 = new TCanvas("canvas12","RecPointsPhiDistribution",1000,600);
-// canvas12->Divide(2,1);
-// TCanvas *canvas13 = new TCanvas("canvas13","RecPointsRDistribution",1000,600);
-// canvas13->Divide(2,1);
+ // -------------This part is to read the number of chunks that were merged
+ Float_t fChnknmbr=0.;
+ FILE * pChunkNumber = fopen ("ChunkNumber.txt","r");
+ fscanf(pChunkNumber, "%f", &fChnknmbr);
+ fclose (pChunkNumber);
+ gStyle->SetPalette(1);
+ float fCNinv=1./fChnknmbr;
+ // printf("\n====================>%f\n\n", fCNinv);
-
if(!mergedfile.IsOpen()){return;}else{printf("file is open\n");}
for(Int_t ispecie=0;ispecie<AliRecoParam::kNSpecies;ispecie++){
//__________________________________________________________________
printf("1\n");
canvas1->Clear();
canvas1->Divide(2,1);
- Bool_t updatecanvas[2];
+
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- for(Int_t iraws=0;iraws<nonrawtot;iraws++){//non expert raws
- canvas1->cd(iraws+1);
-
- sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsnonexpertname[iraws].Data());
- printf("histo name %s ",Rawsnonexpertname[iraws].Data());
+ for(Int_t iraws=2;iraws<4;iraws++){//non expert raws
+ canvas1->cd(iraws-1);
+
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[iraws].Data());
+ printf("histo name %s ",Rawsexpertname[iraws].Data());
historaw=(TH2D*)mergedfile.Get(histoname);
if(historaw){printf("...Found\n");historaw->DrawCopy("colz");}
- else{updatecanvas[iraws]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ else{updatecanvas[iraws-2]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw=NULL;
}//end for
-
+
if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
{
canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname0].Data());
- canvas1->SaveAs(namecanvas);
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
canvas1->Print("SDDQAPlot.ps");
}
//else{delete canvas1; canvas1=NULL;}
for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(1,2);
- printf("2 \n");
- for(Int_t inraws=0;inraws<2;inraws++){//non expert raws
- canvas1->cd(inraws+1);
- sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[inraws].Data());
- printf("histo name %s ",Rawsexpertname[inraws].Data());
- historaw2=(TH1F*)mergedfile.Get(histoname);
- if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[inraws]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
- historaw2=NULL;
- }//end for
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ canvas1->Clear();
+ canvas1->Divide(1,2);
+ printf("2 \n");
+ cannum++;
+ for(Int_t inraws=0;inraws<2;inraws++){//non expert raws
+ canvas1->cd(inraws+1);
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[inraws].Data());
+ printf("histo name %s ",Rawsexpertname[inraws].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+
+
+ // -----------------This is the part were I divide between the number of chunks to normalize the histogram----------------
+
+
+ if(historaw2)
{
- canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname1].Data());
+ printf("...Found\n");
+ if (inraws==1)
+ {
+ historaw2->Multiply(f1,fCNinv);
+
+ }
+ historaw2->DrawCopy();
+ }
+ else{updatecanvas[inraws]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
- }
- //else{delete canvas2; canvas2=NULL;}
-
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- printf("3\n");
- canvas1->Clear();
- canvas1->Divide(2,1);
- for(Int_t iraws=2;iraws<rawtot;iraws++){//non expert raws
- canvas1->cd(iraws-1);
-
- sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[iraws].Data());
- printf("histo name %s",Rawsexpertname[iraws].Data());
- historaw=(TH2D*)mergedfile.Get(histoname);
- if(historaw){printf("...Found\n");historaw->DrawCopy("colz");}else{updatecanvas[iraws-2]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
- historaw=NULL;
- }//end for
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas2; canvas2=NULL;}
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ printf("3\n");
+ cannum++;
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ for(Int_t iraws=0;iraws<2;iraws++){//non expert raws
+ canvas1->cd(iraws+1);
+
+ sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsnonexpertname[iraws].Data());
+ printf("histo name %s",Rawsnonexpertname[iraws].Data());
+ historaw=(TH2D*)mergedfile.Get(histoname);
+ if(historaw){
+ printf("...Found\n");
+ historaw->Multiply(f1,fCNinv);
+ historaw->DrawCopy("colz");
+ }
+ else
+ {updatecanvas[iraws]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas3; canvas3=NULL;}
+
+ //--------------------- new plots
+ printf("================= 4\n");
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ cannum++;
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ for(Int_t iraws=4;iraws<rawtot;iraws++){//non expert raws
+ canvas1->cd(iraws-3);
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[iraws].Data());
+ printf("histo name %s",Rawsexpertname[iraws].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ if(historaw2)
{
- canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname2].Data());
+ printf("...Found\n");
+ historaw2->DrawCopy();
+ }else
+ {
+ updatecanvas[iraws-4]=kFALSE;
+ printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");
+ }
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
- }
- //else{delete canvas3; canvas3=NULL;}
- }//end if directory
-
- //_______________________________________________________________________________
- //rec point
- sprintf(filepath,"ITS/RecPoints/%s",AliRecoParam::GetEventSpecieName(ispecie));
- printf("%s\n",filepath);
- directory2=(TDirectory*)mergedfile.Get(filepath);
- if(directory2){
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(2,1);
- printf("4\n");
- for(Int_t irp=0;irp<2;irp++){//non expert rec point
- canvas1->cd(irp+1);
-
- sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irp].Data());
- printf("histo name %s ",RecPointsnonexpertname[irp].Data());
- histodraw=(TH2F*)mergedfile.Get(histoname);
- if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[irp]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
- histodraw=NULL;
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas3; canvas3=NULL;}
+
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ legend->Clear();
+ printf("14\n");
+ cannum++;
+ for (Int_t i=4;i<rawtot;i++){
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ if(historaw2){
+ fmax=historaw2->GetMaximum();
+ if (i==4){
+ fmaxold=fmax;}
+ if (i!=4){
+ if(fmaxold<fmax){
+ fmaxold=fmax;
+ }
+ }
}
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ historaw2=NULL;
+ }
+ fmaxmargin=1.1*fmaxold;
+ for(Int_t irrpp=4;irrpp<rawtot;irrpp++){//non expert raws
+
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[irrpp].Data());
+ printf("histo name %s",Rawsexpertname[irrpp].Data());
+ sprintf(layer, "layer %d",irrpp-1);
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ gStyle->SetOptStat(0);
+
+ if(historaw2)
{
+ printf("...Found\n");
+ historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
+ historaw2->SetTitle("Raw Relative Occupancy");
+ if (irrpp==4) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
+ if (irrpp!=4) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
+ legend->AddEntry(historaw2,layer,"l");
+ legend->Draw();
canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname3].Data());
+ }
+ else{updatecanvas[irrpp-4]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+
+ }//end if directory
+
+
+
+ //_______________________________________________________________________________
+ //rec point
+ sprintf(filepath,"ITS/RecPoints/%s",AliRecoParam::GetEventSpecieName(ispecie));
+ printf("%s\n",filepath);
+ directory2=(TDirectory*)mergedfile.Get(filepath);
+ if(directory2){
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("4\n");
+ cannum++;
+ for(Int_t irp=0;irp<2;irp++){//non expert rec point
+ canvas1->cd(irp+1);
+
+ sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irp].Data());
+ printf("histo name %s ",RecPointsnonexpertname[irp].Data());
+ histodraw=(TH2F*)mergedfile.Get(histoname);
+ if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[irp]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ histodraw=NULL;
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
}
- // else{delete canvas4; canvas4=NULL;}
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ // else{delete canvas4; canvas4=NULL;}
+
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("5\n");
+ cannum++;
+ for(Int_t irrp=2;irrp<4;irrp++){//non expert raws
+ canvas1->cd(irrp-1);
+
+ sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrp].Data());
+ printf("histo name %s",RecPointsnonexpertname[irrp].Data() );
+ histodraw=(TH2F*)mergedfile.Get(histoname);
+ if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[irrp-2]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ histodraw=NULL;
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas5; canvas5=NULL;}
+
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(1,2);
+ printf("6\n");
+ cannum++;
+ for(Int_t irrpp=4;irrpp<6;irrpp++){//non expert raws
+ canvas1->cd(irrpp-3);
+
+ sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
+ printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+
+ if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[irrpp-4]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas6; canvas6=NULL;}
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(1,2);
+ printf("7\n");
+ cannum++;
+ for(Int_t irrpp=6;irrpp<8;irrpp++){//non expert raws
+ canvas1->cd(irrpp-5);
+ // printf("histo name\n");
+ sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
+ printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[irrpp-6]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas7; canvas7=NULL;}
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ printf("8\n");
+ cannum++;
+ //canvas1->Divide(2,1);
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[rlocaldistro].Data());
+ printf("histo name %s",RecPointsexpertname[rlocaldistro].Data());
+ histodraw=(TH2F*)mergedfile.Get(histoname);
+
+ if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[0]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ histodraw=NULL;
+ if(updatecanvas[0]==kTRUE){
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas8; canvas8=NULL;}
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("9\n");
+ cannum++;
+ for(Int_t i=1;i<3;i++)
+ {
+ canvas1->cd(i);
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(2,1);
- printf("5\n");
- for(Int_t irrp=2;irrp<4;irrp++){//non expert raws
- canvas1->cd(irrp-1);
-
- sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrp].Data());
- printf("histo name %s",RecPointsnonexpertname[irrp].Data() );
+ sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ printf("histo name %s",RecPointsexpertname[i].Data());
histodraw=(TH2F*)mergedfile.Get(histoname);
- if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[irrp-2]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+
+ if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[i-1]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
histodraw=NULL;
}
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
- {
- canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname4].Data());
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas9; canvas9=NULL;}
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("10\n");
+ cannum++;
+ for(Int_t i=8;i<10;i++)
+ {
+ canvas1->cd(i-7);
+
+ sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[i].Data());
+ histodraw=(TH2F*)mergedfile.Get(histoname);
+ printf("histo name %s",RecPointsnonexpertname[i].Data());
+
+ if(histodraw){
+ printf("...Found\n");histodraw->Multiply(f1,fCNinv);
+ histodraw->DrawCopy("colz");
+ }else{updatecanvas[i-8]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ histodraw=NULL;
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
}
- //else{delete canvas5; canvas5=NULL;}
-
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas10; canvas10=NULL;}
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(1,2);
- printf("6\n");
- for(Int_t irrpp=4;irrpp<6;irrpp++){//non expert raws
- canvas1->cd(irrpp-3);
-
- sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
- printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(1,2);
+ printf("11\n");
+ cannum++;
+ for(Int_t i=3;i<5;i++)
+ {
+ canvas1->cd(i-2);
+
+ sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ printf("histo name %s",RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
- if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[irrpp-4]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+
+ if(historaw2)
+ {
+ printf("...Found\n");
+ if (i==4){historaw2->Multiply(f1,fCNinv);}
+ historaw2->DrawCopy();
+ }
+ else{updatecanvas[i-3]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
- }//end for
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
- {
- canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname5].Data());
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
}
- //else{delete canvas6; canvas6=NULL;}
-
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(1,2);
- printf("7\n");
- for(Int_t irrpp=6;irrpp<8;irrpp++){//non expert raws
- canvas1->cd(irrpp-5);
- // printf("histo name\n");
- sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
- printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas11; canvas11=NULL;}
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("12\n");
+ cannum++;
+ for(Int_t i=5;i<7;i++)
+ {
+ canvas1->cd(i-4);
+
+ sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ printf("histo name %s",RecPointsexpertname[i].Data());
historaw2=(TH1F*)mergedfile.Get(histoname);
- if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[irrpp-6]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+
+ if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-5]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
historaw2=NULL;
- }//end for
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
- {
- canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname6].Data());
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
}
- //else{delete canvas7; canvas7=NULL;}
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //else{delete canvas12; canvas12=NULL;}
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- printf("8\n");
- //canvas1->Divide(2,1);
- sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[rlocaldistro].Data());
- printf("histo name %s",RecPointsexpertname[rlocaldistro].Data());
- histodraw=(TH2F*)mergedfile.Get(histoname);
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("13\n");
+ cannum++;
+ for(Int_t i=7;i<9;i++)
+ {
+ canvas1->cd(i-6);
+ sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ printf("histo name %s",RecPointsexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
- if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[0]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
- histodraw=NULL;
- if(updatecanvas[0]==kTRUE){
+ if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-7]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname7].Data());
- canvas1->SaveAs(namecanvas);
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
canvas1->Print("SDDQAPlot.ps");
}
- //else{delete canvas8; canvas8=NULL;}
-
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(2,1);
- printf("9\n");
- for(Int_t i=8;i<10;i++)
- {
- canvas1->cd(i-7);
-
- sprintf(histoname,"ITS/RecPoints/%s/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[i].Data());
- printf("histo name %s",RecPointsnonexpertname[i].Data());
- histodraw=(TH2F*)mergedfile.Get(histoname);
+
- if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[i-8]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
- histodraw=NULL;
+ //superimposed
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ legend->Clear();
+ printf("14\n");
+ cannum++;
+ for (Int_t i=4;i<6;i++){
+ sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ if(historaw2){
+ fmax=historaw2->GetMaximum();
+ if (i==4){
+ fmaxold=fmax;}
+ if (i!=4){
+ if(fmaxold<fmax){
+ fmaxold=fmax;
+ }
}
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ }
+ historaw2=NULL;
+ }
+ fmaxmargin=1.1*fmaxold;
+ for(Int_t irrpp=4;irrpp<6;irrpp++){//non expert raws
+
+ sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
+ printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
+ sprintf(layer, "layer %d",irrpp-1);
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ gStyle->SetOptStat(0);
+
+ if(historaw2)
{
+ printf("...Found\n");
+ historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
+ historaw2->SetTitle("Charge");
+ if (irrpp==4) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
+ if (irrpp!=4) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
+ legend->AddEntry(historaw2,layer,"l");
+ legend->Draw();
canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname8].Data());
+ }
+ else{updatecanvas[irrpp-4]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
- }
- //else{delete canvas9; canvas9=NULL;}
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(2,1);
- printf("10\n");
- for(Int_t i=1;i<3;i++)
- {
- canvas1->cd(i);
-
- sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
- histodraw=(TH2F*)mergedfile.Get(histoname);
- printf("histo name %s",RecPointsexpertname[i].Data());
-
- if(histodraw){printf("...Found\n");histodraw->DrawCopy("colz");}else{updatecanvas[i-1]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
- histodraw=NULL;
}
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ canvas1->Print("SDDQAPlot.ps");
+ }
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ legend->Clear();
+ printf("15\n");
+ cannum++;
+ for (Int_t i=6;i<8;i++){
+ sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ if(historaw2){
+ fmax=historaw2->GetMaximum();
+ if (i==6){
+ fmaxold=fmax;}
+ if (i!=6){
+ if(fmaxold<fmax){
+ fmaxold=fmax;}}}
+ fmaxmargin=1.1*fmaxold;
+ }
+ for(Int_t irrpp=6;irrpp<8;irrpp++){//non expert raws
+
+ sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[irrpp].Data());
+ printf("histo name %s",RecPointsnonexpertname[irrpp].Data());
+ sprintf(layer, "layer %d",irrpp-3);
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ gStyle->SetOptStat(0);
+
+ if(historaw2)
{
+ printf("...Found\n");
+ historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
+ historaw2->SetTitle("Drift Time");
+ if (irrpp==6) {historaw2->SetStats(0); historaw2->SetLineColor(2);historaw2->DrawCopy();}
+ if (irrpp!=6) {historaw2->SetStats(0); historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
+ legend->AddEntry(historaw2,layer,"l");
+ legend->Draw();
canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname9].Data());
+ }
+ else{updatecanvas[irrpp-6]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ printf("%s\n%s\n",historaw2->GetName(),historaw2->GetTitle());
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
}
- //else{delete canvas10; canvas10=NULL;}
+ canvas1->Print("SDDQAPlot.ps");
+ }
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(1,2);
- printf("11\n");
- for(Int_t i=3;i<5;i++)
- {
- canvas1->cd(i-2);
+ //------------------------------------------- new plot
- sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
- printf("histo name %s",RecPointsexpertname[i].Data());
- historaw2=(TH1F*)mergedfile.Get(histoname);
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("13\n");
+ cannum++;
+ for(Int_t i=9;i<11;i++)
+ {
+ canvas1->cd(i-8);
+ sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ printf("histo name %s",RecPointsexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
- if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-3]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
- historaw2=NULL;
+ if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-9]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+
+ //------------------------------------
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ legend->Clear();
+ printf("14\n");
+ cannum++;
+ for (Int_t i=9;i<11;i++){
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ if(historaw2){
+ fmax=historaw2->GetMaximum();
+ if (i==4){
+ fmaxold=fmax;}
+ if (i!=4){
+ if(fmaxold<fmax){
+ fmaxold=fmax;
+ }
}
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ }
+ historaw2=NULL;
+ }
+ fmaxmargin=1.1*fmaxold;
+ for(Int_t irrpp=9;irrpp<11;irrpp++){//non expert raws
+
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[irrpp].Data());
+ printf("histo name %s",RecPointsexpertname[irrpp].Data());
+ sprintf(layer, "layer %d",irrpp-6);
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ gStyle->SetOptStat(0);
+
+ if(historaw2)
{
+ printf("...Found\n");
+ historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
+ historaw2->SetTitle("RecPoint Relative Occupancy");
+ if (irrpp==9) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
+ if (irrpp!=9) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
+ legend->AddEntry(historaw2,layer,"l");
+ legend->Draw();
canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname10].Data());
+ }
+ else{updatecanvas[irrpp-9]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
}
- //else{delete canvas11; canvas11=NULL;}
+ canvas1->Print("SDDQAPlot.ps");
+ }
+
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("13\n");
+ cannum++;
+ for(Int_t i=11;i<13;i++)
+ {
+ canvas1->cd(i-10);
+ sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ printf("histo name %s",RecPointsexpertname[i].Data());
+ historaw=(TH2D*)mergedfile.Get(histoname);
+
+ if(historaw){printf("...Found\n");historaw->DrawCopy("colz");}else{updatecanvas[i-11]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(2,1);
- printf("12\n");
- for(Int_t i=5;i<7;i++)
- {
- canvas1->cd(i-4);
- sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
- printf("histo name %s",RecPointsexpertname[i].Data());
- historaw2=(TH1F*)mergedfile.Get(histoname);
+ //-----------------
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("13\n");
+ cannum++;
+ for(Int_t i=13;i<15;i++)
+ {
+ canvas1->cd(i-12);
+ sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ printf("histo name %s",RecPointsexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+
+ if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-13]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ //--------------------------------------------
- if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-5]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
- historaw2=NULL;
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ legend->Clear();
+ printf("14\n");
+ cannum++;
+ for (Int_t i=13;i<15;i++){
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ if(historaw2){
+ fmax=historaw2->GetMaximum();
+ if (i==13){
+ fmaxold=fmax;}
+ if (i!=13){
+ if(fmaxold<fmax){
+ fmaxold=fmax;
+ }
}
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ }
+ historaw2=NULL;
+ }
+ fmaxmargin=1.1*fmaxold;
+ for(Int_t irrpp=13;irrpp<15;irrpp++){//non expert raws
+
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[irrpp].Data());
+ printf("histo name %s",RecPointsexpertname[irrpp].Data());
+ sprintf(layer, "layer %d",irrpp-10);
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ gStyle->SetOptStat(0);
+
+ if(historaw2)
{
+ printf("...Found\n");
+ historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
+ historaw2->SetTitle("Rec2Raw Ratio");
+ if (irrpp==13) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
+ if (irrpp!=13) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
+ legend->AddEntry(historaw2,layer,"l");
+ legend->Draw();
canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname11].Data());
+ }
+ else{updatecanvas[irrpp-13]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
}
- //else{delete canvas12; canvas12=NULL;}
+ canvas1->Print("SDDQAPlot.ps");
+ }
- for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
- canvas1->Clear();
- canvas1->Divide(2,1);
- printf("13\n");
- for(Int_t i=7;i<9;i++)
- {
- canvas1->cd(i-6);
- sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
- printf("histo name %s",RecPointsexpertname[i].Data());
- historaw2=(TH1F*)mergedfile.Get(histoname);
+
+ //--------------------------------------
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ canvas1->Divide(2,1);
+ printf("13\n");
+ cannum++;
+ for(Int_t i=15;i<17;i++)
+ {
+ canvas1->cd(i-14);
+ sprintf(histoname,"ITS/RecPoints/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ printf("histo name %s",RecPointsexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+
+ if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-15]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
+ canvas1->SaveAs(namecanvas);
+ }
+ canvas1->Print("SDDQAPlot.ps");
+ }
+
- if(historaw2){printf("...Found\n");historaw2->DrawCopy();}else{updatecanvas[i-7]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
- historaw2=NULL;
+ //--------------------------------------------
+
+ for(Int_t i=0;i<2;i++){updatecanvas[i]=kTRUE;}
+ canvas1->Clear();
+ legend->Clear();
+ printf("14\n");
+ cannum++;
+ for (Int_t i=15;i<17;i++){
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[i].Data());
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ if(historaw2){
+ fmax=historaw2->GetMaximum();
+ if (i==15){
+ fmaxold=fmax;}
+ if (i!=15){
+ if(fmaxold<fmax){fmaxold=fmax;}
}
- if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ }
+ historaw2=NULL;
+ }
+ fmaxmargin=1.1*fmaxold;
+ for(Int_t irrpp=15;irrpp<17;irrpp++){//non expert raws
+
+ sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[irrpp].Data());
+ printf("histo name %s",RecPointsexpertname[irrpp].Data());
+ sprintf(layer, "layer %d",irrpp-12);
+ historaw2=(TH1F*)mergedfile.Get(histoname);
+ gStyle->SetOptStat(0);
+
+ if(historaw2)
{
+ printf("...Found\n");
+ historaw2->GetYaxis()->SetRangeUser(0,fmaxmargin);
+ historaw2->SetTitle("RecPoint dEdx");
+ if (irrpp==15) {historaw2->SetStats(0);historaw2->SetLineColor(2);historaw2->DrawCopy();}
+ if (irrpp!=15) {historaw2->SetStats(0);historaw2->SetLineColor(4);historaw2->DrawCopy("same");}
+ legend->AddEntry(historaw2,layer,"l");
+ legend->Draw();
canvas1->Update();
- sprintf(namecanvas,"%s.eps",canvassavedname[canvname12].Data());
+ }
+ else{updatecanvas[irrpp-15]=kFALSE;printf("...Not Found: the histogram or this QA has been done before the histograms was added to QA\n");}
+ historaw2=NULL;
+ }//end for
+ if(updatecanvas[0]==kTRUE||updatecanvas[1]==kTRUE)
+ {
+ canvas1->Update();
+ if(kDoEps){
+ sprintf(namecanvas,"%s.eps",canvassavedname[cannum].Data());
canvas1->SaveAs(namecanvas);
- canvas1->Print("SDDQAPlot.ps");
}
- canvas1->Print("SDDQAPlot.ps]");
- //else{delete canvas13; canvas13=NULL;}
- directory2=NULL;
- }//end directory
+ canvas1->Print("SDDQAPlot.ps");
+ }
+ canvas1->Print("SDDQAPlot.ps]");
+
+
+ for(Int_t isave=0;isave<9;isave++){
+ if(isave<4)sprintf(histoname,"%s/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsnonexpertname[histo2savenumber[isave]].Data());
+ if(isave>4&&isave<7)sprintf(histoname,"%s/Expert/%s_%s",filepath,AliRecoParam::GetEventSpecieName(ispecie),RecPointsexpertname[histo2savenumber[isave]].Data());
+ if(isave>7)sprintf(histoname,"ITS/Raws/%s/Expert/%s_%s",AliRecoParam::GetEventSpecieName(ispecie),AliRecoParam::GetEventSpecieName(ispecie),Rawsexpertname[histo2savenumber[isave]].Data());
+ printf("file2save name:\t %s\n",histoname);
+ histo2save[isave]=(TH1F*)mergedfile.Get(histoname);
+ }
+
+ TFile file2savefortrend(trendfile,"recreate");
+ file2savefortrend.cd();
+ for(Int_t iss=0;iss<9;iss++){printf("Saved %d\n",iss); histo2save[iss]->Write();}
+ file2savefortrend.Close();
+
+ //else{delete canvas13; canvas13=NULL;}
+ //directory2=NULL;
+ }//end directory
+
+
+
}//end for
+
+
+
delete directory;
directory=NULL;
delete directory2;
--- /dev/null
+#if !defined(__CINT__) || defined(__MAKECINT__)
+#include <TFile.h>
+#include <TFileMerger.h>
+#include <TAlienFile.h>
+//#include <TExec.h>
+#include <TSystem.h>
+#include <TGrid.h>
+#include <TGridResult.h>
+#include <Riostream.h>
+#include <TObjArray.h>
+#include <TClass.h>
+#include <TH1F.h>
+#include <TStyle.h>
+#include <TCanvas.h>
+#include <TLegend.h>
+#include <TMath.h>
+#endif
+void TrendQASDD( Char_t filelist[150]="LHC10b.txt", Char_t filename[20]="FileQAtrend",Char_t pass[7]="pass1",Int_t year=2010, Char_t period[7]="LHC10b", Bool_t kUseOriginalFile=kFALSE , Char_t especiename[50]="LowMultiplicity")
+{
+ Int_t m;
+ vector<Int_t> n;
+ //Char_t FileName[200];//="/home/msicilia/Desktop/createlists";
+ char FileName[150];
+ Char_t FileName1[150];
+ sprintf(FileName1,"%s",filelist);
+ FILE * pFile;
+ pFile = fopen (FileName1,"r");
+ while (fscanf (pFile,"%d\n",&m)==1){
+
+ if(kUseOriginalFile)sprintf(FileName,"run%d/%s/File.QA.%d.%s.%s.Run.%d.root",m,pass,year,period,pass,m);
+ else sprintf(FileName, "run%d/%s/%s.root", m,pass,filename);
+ if(gSystem->Exec(Form("ls %s >/dev/null 2>&1", FileName))==0){
+ TFile mergedfile(FileName);
+ if(kUseOriginalFile){
+ if(mergedfile.GetKey("ITS")==0x0){
+ printf("Run %d, In this run ITS QA has not been executed.-- Exit file\n",m);
+ continue;
+ }
+ else {n.push_back(m);} //this is to grow the vector
+ }else {n.push_back(m);}
+ }
+ else{printf("Run %d,: no file %s present -- Continue\n",m, FileName);continue;}
+ }
+ fclose (pFile);
+
+ Char_t filepath[200];
+ TH1F *histocharge;
+ TH1F *histochargetrend[2];
+
+ char name[100];
+ char title[100];
+
+ for(Int_t layer=3;layer<5;layer++)
+ {
+ sprintf(name,"histochargetrend_%d", layer);
+ sprintf(title,"Trend of the Charge of the Layer %d", layer);
+ histochargetrend[layer-3]=new TH1F(name,title,n.size(),0.,(float)n.size());
+ }
+
+ gStyle->SetOptStat(0);
+ TCanvas *canvas[6];
+
+ for(Int_t i=0;i<3;i++)
+ {
+ sprintf(name,"canvas%i",i);
+ canvas[i]=new TCanvas(name,name);
+ canvas[i]->Divide(2,1);
+ }
+
+ TLegend *legend =new TLegend(0.83,0.8,0.97,0.2);
+ TLegend *legend2=new TLegend(0.83,0.8,0.97,0.2);
+
+ Char_t runnmbr[20];
+ Float_t fmax=0;
+ Float_t fmaxold=0;
+ Float_t fmaxmargin=0;
+
+ Float_t fmaxtime=0;
+ Float_t fmaxoldtime=0;
+ Float_t fmaxmargintime=0;
+
+ Int_t j=0;
+ Int_t hh=0;
+ //Int_t kk=0;
+ char buffer [10];
+
+ printf("%d \n",n.size());
+
+ for(Int_t layer=3;layer<5;layer++){
+ for (Int_t irun=0; irun<(Int_t)(n.size());irun++)
+ {
+ if(kUseOriginalFile)sprintf(FileName,"run%d/%s/File.QA.%d.%s.%s.Run.%d.root",n[irun],pass,year,period,pass,n[irun]);
+ else sprintf(FileName, "run%d/%s/%s.root", n[irun],pass,filename);
+
+ TFile mergedfile(FileName);
+ if(kUseOriginalFile){
+ if(mergedfile.GetKey("ITS")==0x0){
+ printf("Run %d, In this run ITS QA has not been executed.-- Exit file\n",n[irun]);
+ continue;
+ }
+ }
+ if(kUseOriginalFile)sprintf(filepath,"ITS/RecPoints/%s/%s_SDDLay%dTotCh",especiename,especiename,layer);
+ else sprintf(filepath,"%s_SDDLay%dTotCh",especiename,layer);
+ histocharge=(TH1F*)(mergedfile.Get(filepath));
+
+ if(histocharge)
+ {
+ sprintf(buffer, "%d", n[irun]);
+ histochargetrend[layer-3]->Fill(buffer,histocharge->GetMean());
+ histochargetrend[layer-3]->SetMinimum(0.);
+ histochargetrend[layer-3]->SetMaximum(250.);
+ histochargetrend[layer-3]->GetXaxis()->SetTitle("# Run");
+ histochargetrend[layer-3]->GetYaxis()->SetTitle("keV");
+ fmax=histocharge->GetMaximum();
+ if (j==0) { fmaxold=fmax;}
+ if (j!=0) { if(fmaxold<fmax) { fmaxold=fmax;} }
+ j++;
+ }//end if histocharge
+
+ histocharge=NULL;
+
+ //drift time
+
+ if(kUseOriginalFile)sprintf(filepath,"ITS/RecPoints/%s/%s_SDDdrifttime_Layer%d",especiename,especiename,layer);
+ else sprintf(filepath,"%s_SDDdrifttime_Layer%d",especiename,layer);
+ histocharge=(TH1F*)(mergedfile.Get(filepath));
+
+ if(histocharge)
+ {
+ sprintf(buffer, "%d", n[irun]);
+ fmaxtime=histocharge->GetMaximum();
+ if (j==0) { fmaxoldtime=fmaxtime;}
+ if (j!=0) { if(fmaxoldtime<fmaxtime) { fmaxoldtime=fmaxtime;} }
+ j++;
+ }//end if histocharge
+
+ histocharge=NULL;
+
+ }//end for run
+ fmaxmargin=fmaxold;
+ fmaxmargintime=1.023*fmaxoldtime;
+ printf("Max Y Range is %4.1f time %4.1f\n",fmaxmargin,fmaxmargintime);
+
+ canvas[0]->SetFillColor(46);
+ canvas[0]->cd(layer-2)->SetFillColor(46);
+ histochargetrend[layer-3]->SetMarkerStyle(20);
+ histochargetrend[layer-3]->SetMarkerSize(1);
+ histochargetrend[layer-3]->SetMarkerColor(kGreen+4);
+ histochargetrend[layer-3]->SetFillColor(42);
+ canvas[0]->cd(layer-2)->SetFrameFillColor(kAzure-9);
+ histochargetrend[layer-3]->SetBarWidth(0.9);
+ histochargetrend[layer-3]->SetBarOffset(0.05);
+ histochargetrend[layer-3]->DrawCopy("P");
+
+ j=0;
+ hh=0;
+ Int_t color=0;
+ for (Int_t irun=0; irun< (Int_t)(n.size());irun++)
+ {
+
+ if(kUseOriginalFile) sprintf(FileName, "run%d/%s/File.QA.%d.%s.%s.Run.%d.root", n[irun],pass,year,period,pass,n[irun]);
+ else sprintf(FileName, "run%d/%s/%s.root", n[irun],pass,filename);
+
+ TFile mergedfile(FileName);
+ if(kUseOriginalFile){
+ if(mergedfile.GetKey("ITS")==0x0){
+ printf("run %d In this run ITS QA has not been executed.--- Exit file\n",n[irun]);
+ continue;
+ }
+ }
+ if(kUseOriginalFile)sprintf(filepath,"ITS/RecPoints/%s/%s_SDDLay%dTotCh",especiename,especiename,layer);
+ else sprintf(filepath,"%s_SDDLay%dTotCh",especiename,layer);
+ sprintf(runnmbr,"run%d",n[irun]);
+ printf("%s %s",filepath,runnmbr);
+ histocharge=(TH1F*)(mergedfile.Get(filepath));
+ if(histocharge)
+ {
+ color=2+j;
+ if (color>=12 && color<=27) { color=j+28; }
+ histocharge->SetLineColor(color);
+ if (j==0)
+ {
+ canvas[1]->cd(layer-2);
+
+ canvas[1]->cd(layer-2)->SetFrameFillColor(kGray+3);
+ histocharge->SetMaximum(fmaxmargin);
+ histocharge->GetYaxis()->SetTitle("");
+ histocharge->DrawNormalized();
+ histocharge->SetMaximum(fmaxmargin);
+ canvas[1]->cd(layer-2)->Update();
+ }
+ if (j!=0)
+ {
+
+ canvas[1]->cd(layer-2);
+ histocharge->DrawNormalized("same");
+ canvas[1]->cd(layer-2)->Update();
+ }
+ if(layer==3)legend->AddEntry(histocharge,runnmbr,"l");
+ canvas[1]->cd(layer-2);
+ legend->Draw();
+ canvas[1]->Update();
+ j++;
+ printf("...Found\n");
+
+ }//end if histocharge
+ else{printf("...Not Found....\n");}
+ histocharge=NULL;
+
+ //======================================================== Drift Time
+ if(kUseOriginalFile)sprintf(filepath,"ITS/RecPoints/%s/%s_SDDdrifttime_Layer%d",especiename,especiename,layer);
+ else sprintf(filepath,"%s_SDDdrifttime_Layer%d",especiename,layer);
+ sprintf(runnmbr,"run %d",n[irun]);
+ printf("%s %s",filepath,runnmbr);
+ histocharge=(TH1F*)(mergedfile.Get(filepath));
+ if(histocharge)
+ {
+ printf("...Found\n");
+ histocharge->SetLineColor(color);
+ if (hh==0)
+ {
+ canvas[2]->cd(layer-2);
+ canvas[2]->cd(layer-2)->SetFrameFillColor(kGray+3);
+ histocharge->SetMaximum(fmaxmargintime);
+ histocharge->GetYaxis()->SetTitle("");
+ histocharge->DrawNormalized();
+ histocharge->SetMaximum(fmaxmargintime);
+ }
+ if (hh!=0)
+ {
+ canvas[2]->cd(layer-2);
+ histocharge->DrawNormalized("same");
+ }
+ if(layer==3)legend2->AddEntry(histocharge,runnmbr,"l");
+ canvas[2]->cd(layer-2);
+ legend2->Draw();
+ canvas[2]->Update();
+ hh++;
+ }//end if histocharge
+ else{printf("...Not Found....\n");}
+ histocharge=NULL;
+ }//enf for
+ }//end for layer
+
+ TFile trendfile(Form("SDDQAtrend%s%s.root",period,pass),"recreate");
+ trendfile.cd();
+ for(Int_t ican=0;ican<3;ican++)canvas[ican]->Write();
+ trendfile.Close();
+ Char_t psfile[50];
+ sprintf(psfile,"SDDtrend%s%s.ps",period,pass);
+ canvas[0]->Print(Form("%s[",psfile));
+ for(Int_t ifile=0;ifile<3;ifile++){canvas[ifile]->Print(psfile);}
+
+ canvas[2]->Print(Form("%s]",psfile));
+
+
+ delete histocharge;
+
+
+
+
+}//end macro
git://git.uio.no / u / mrichter / AliRoot.git / commitdiff