void AliAnalysisTaskPt::Terminate(Option_t *)
{
// Processing when the event loop is ended
-
+
+ AliInfo(Form(" *** %s Report:", GetName())) ;
+
TCanvas *c1 = new TCanvas("c1","Pt",10,10,310,310);
c1->SetFillColor(10);
c1->SetHighLightColor(10);
c1->cd(1)->SetLeftMargin(0.15);
c1->cd(1)->SetBottomMargin(0.15);
- c1->cd(1)->SetLogy();
+ if (fhPt->GetMaximum() > 0 )
+ c1->cd(1)->SetLogy();
fOutputContainer = (TObjArray*)GetOutputData(0);
fhPt = (TH1F*)fOutputContainer->At(0);
if (fhPt) fhPt->DrawCopy("E");
c1->Print("TracksPt.eps");
char line[1024] ;
- sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
gROOT->ProcessLine(line);
sprintf(line, ".!rm -fR *.eps");
gROOT->ProcessLine(line);
fhEMCALPhotons = (TH1I*)fOutputContainer->At(5);
fhEMCALInvariantMass = (TH1D*)fOutputContainer->At(6);
fhEMCALDigitsEvent = (TH1I*)fOutputContainer->At(7);
- printf("EMCALEnergy Mean : %5.3f , RMS : %5.3f \n", fhEMCALEnergy->GetMean(), fhEMCALEnergy->GetRMS() ) ;
- printf("EMCALDigits Mean : %5.3f , RMS : %5.3f \n", fhEMCALDigits->GetMean(), fhEMCALDigits->GetRMS() ) ;
- printf("EMCALRecParticles Mean : %5.3f , RMS : %5.3f \n", fhEMCALRecParticles->GetMean(), fhEMCALRecParticles->GetRMS() ) ;
- printf("EMCALPhotons Mean : %5.3f , RMS : %5.3f \n", fhEMCALPhotons->GetMean(), fhEMCALPhotons->GetRMS() ) ;
- printf("EMCALInvariantMass Mean : %5.3f , RMS : %5.3f \n", fhEMCALInvariantMass->GetMean(), fhEMCALInvariantMass->GetRMS() ) ;
- printf("EMCALDigitsEvent Mean : %5.3f , RMS : %5.3f \n", fhEMCALDigitsEvent->GetMean(), fhEMCALDigitsEvent->GetRMS() ) ;
+
+ AliInfo(Form(" *** %s Report:", GetName())) ;
+ printf(" EMCALEnergy Mean : %5.3f , RMS : %5.3f \n", fhEMCALEnergy->GetMean(), fhEMCALEnergy->GetRMS() ) ;
+ printf(" EMCALDigits Mean : %5.3f , RMS : %5.3f \n", fhEMCALDigits->GetMean(), fhEMCALDigits->GetRMS() ) ;
+ printf(" EMCALRecParticles Mean : %5.3f , RMS : %5.3f \n", fhEMCALRecParticles->GetMean(), fhEMCALRecParticles->GetRMS() ) ;
+ printf(" EMCALPhotons Mean : %5.3f , RMS : %5.3f \n", fhEMCALPhotons->GetMean(), fhEMCALPhotons->GetRMS() ) ;
+ printf(" EMCALInvariantMass Mean : %5.3f , RMS : %5.3f \n", fhEMCALInvariantMass->GetMean(), fhEMCALInvariantMass->GetRMS() ) ;
+ printf(" EMCALDigitsEvent Mean : %5.3f , RMS : %5.3f \n", fhEMCALDigitsEvent->GetMean(), fhEMCALDigitsEvent->GetRMS() ) ;
TCanvas * cEMCAL = new TCanvas("EMCAL", "EMCAL ESD Test", 400, 10, 600, 700);
cEMCAL->Divide(3, 2) ;
cEMCAL->cd(1) ;
- gPad->SetLogy();
+ if ( fhEMCALEnergy->GetMaximum() > 0. )
+ gPad->SetLogy();
fhEMCALEnergy->SetAxisRange(0, 25.);
fhEMCALEnergy->SetXTitle("Energy (GeV)");
fhEMCALEnergy->Draw();
cEMCAL->cd(2) ;
- gPad->SetLogy();
+ if ( fhEMCALDigits->GetMaximum() > 0. )
+ gPad->SetLogy();
fhEMCALDigits->SetAxisRange(0, 25.);
fhEMCALDigits->SetXTitle("DigitsPerCluster");
fhEMCALDigits->Draw();
cEMCAL->cd(3) ;
- gPad->SetLogy();
+ if ( fhEMCALRecParticles->GetMaximum() > 0. )
+ gPad->SetLogy();
fhEMCALRecParticles->SetAxisRange(0, 25.);
fhEMCALRecParticles->SetXTitle("RecParticles");
fhEMCALRecParticles->Draw();
cEMCAL->cd(4) ;
- gPad->SetLogy();
+ if ( fhEMCALPhotons->GetMaximum() > 0. )
+ gPad->SetLogy();
fhEMCALPhotons->SetAxisRange(0, 25.);
fhEMCALPhotons->SetXTitle("Photons");
fhEMCALPhotons->Draw();
fhEMCALInvariantMass->Draw();
cEMCAL->cd(6) ;
- gPad->SetLogy();
+ if ( fhEMCALDigitsEvent->GetMaximum() > 0. )
+ gPad->SetLogy();
fhEMCALDigitsEvent->SetAxisRange(0, 40.);
fhEMCALDigitsEvent->SetXTitle("DigitsPerEvent");
fhEMCALDigitsEvent->Draw();
cEMCAL->Print("EMCAL.eps");
char line[1024] ;
- sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
gROOT->ProcessLine(line);
sprintf(line, ".!rm -fR *.eps");
gROOT->ProcessLine(line);
fhHMPIDProb[2] = (TH1F*)fOutputContainer->At(6);
fhHMPIDProb[3] = (TH1F*)fOutputContainer->At(7);
fhHMPIDProb[4] = (TH1F*)fOutputContainer->At(8);
-
+
+ AliInfo(Form(" *** %s Report:", GetName())) ;
+
Float_t n = 1.292 ; //mean freon ref idx
TF1 * hHMPIDpPi = new TF1("RiPiTheo", "acos(sqrt(x*x+[0]*[0])/(x*[1]))", 1.2, 7) ;
hHMPIDpPi->SetLineWidth(1) ;
cHMPID->Print("HMPID.eps");
char line[1024] ;
- sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
gROOT->ProcessLine(line);
sprintf(line, ".!rm -fR *.eps");
gROOT->ProcessLine(line);
fhChi2match = new TH1F("hChi2match","Chi2 of trig/track matching",100, 0.,20.);
// create output container
- fOutputContainer = new TObjArray(10) ;
+ fOutputContainer = new TObjArray(12) ;
fOutputContainer->SetName(GetName()) ;
fOutputContainer->AddAt(fhMUONVertex, 0) ;
fOutputContainer->AddAt(fhMUONMult, 1) ;
effMatch = 100 * fnTrackTrig / ftracktot ;
fheffMatchT->Fill(effMatch);}
- AliInfo(Form("Terminate %s:", GetName())) ;
+ AliInfo(Form(" *** %s Report:", GetName())) ;
+
fOutputContainer = (TObjArray*)GetOutputData(0);
fhMUONVertex = (TH1F*)fOutputContainer->At(0);
fhMUONMult = (TH1F*)fOutputContainer->At(1);
fhUSHighpt=(TH1F*)fOutputContainer->At(7);
fhLSLowpt=(TH1F*)fOutputContainer->At(8);
fhLSHighpt=(TH1F*)fOutputContainer->At(9);
- //______________________________________________________________________
- printf("===================================================\n") ;
- printf("================ %s ESD SUMMARY ==============\n", GetName()) ;
- printf(" \n") ;
+
printf(" Total number of processed events %d \n", fnevents) ;
- printf("\n") ;
- printf("\n") ;
- printf("Table 1: \n") ;
- printf("===================================================\n") ;
- printf(" Global Trigger output Low pt High pt \n") ;
- printf(" number of Single :\t");
- printf("%i\t", fSLowpt) ;
+ printf(" \n") ;
+ printf(" \n") ;
+ printf(" Table 1: \n") ;
+ printf(" ===================================================\n") ;
+ printf(" Global Trigger output Low pt High pt \n") ;
+ printf(" number of Single :\t");
+ printf(" %i\t", fSLowpt) ;
printf("\n");
- printf(" number of UnlikeSign pair :\t");
- printf("%i\t%i\t", fUSLowpt, fUSHighpt) ;
+ printf(" number of UnlikeSign pair :\t");
+ printf(" %i\t%i\t", fUSLowpt, fUSHighpt) ;
printf("\n");
- printf(" number of LikeSign pair :\t");
- printf("%i\t%i\t", fLSLowpt, fLSHighpt) ;
+ printf(" number of LikeSign pair :\t");
+ printf(" %i\t%i\t", fLSLowpt, fLSHighpt) ;
printf("\n");
- printf("===================================================\n") ;
- printf("\n") ;
- printf("matching efficiency with the trigger for single tracks = %2d %% \n", effMatch);
+ printf(" matching efficiency with the trigger for single tracks = %2d %% \n", effMatch);
printf("\n") ;
- printf("===================================================\n") ;
TCanvas * cMUON1 = new TCanvas("cMUON1", "MUON ESD Vert & Mult", 400, 10, 600, 700) ;
cMUON1->Divide(1,2) ;
cMUON4->Print("MUON4.eps") ;
char line[1024] ;
- sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
gROOT->ProcessLine(line);
sprintf(line, ".!rm -fR *.eps");
gROOT->ProcessLine(line);
fhPHOSInvariantMass = (TH1D*)fOutputContainer->At(6);
fhPHOSDigitsEvent = (TH1I*)fOutputContainer->At(7);
- printf("PHOSEnergy Mean : %5.3f , RMS : %5.3f \n", fhPHOSEnergy->GetMean(), fhPHOSEnergy->GetRMS() ) ;
- printf("PHOSDigits Mean : %5.3f , RMS : %5.3f \n", fhPHOSDigits->GetMean(), fhPHOSDigits->GetRMS() ) ;
- printf("PHOSRecParticles Mean : %5.3f , RMS : %5.3f \n", fhPHOSRecParticles->GetMean(), fhPHOSRecParticles->GetRMS() ) ;
- printf("PHOSPhotons Mean : %5.3f , RMS : %5.3f \n", fhPHOSPhotons->GetMean(), fhPHOSPhotons->GetRMS() ) ;
- printf("PHOSInvariantMass Mean : %5.3f , RMS : %5.3f \n", fhPHOSInvariantMass->GetMean(), fhPHOSInvariantMass->GetRMS() ) ;
- printf("PHOSDigitsEvent Mean : %5.3f , RMS : %5.3f \n", fhPHOSDigitsEvent->GetMean(), fhPHOSDigitsEvent->GetRMS() ) ;
+ AliInfo(Form(" *** %s Report:", GetName())) ;
+ printf(" PHOSEnergy Mean : %5.3f , RMS : %5.3f \n", fhPHOSEnergy->GetMean(), fhPHOSEnergy->GetRMS() ) ;
+ printf(" PHOSDigits Mean : %5.3f , RMS : %5.3f \n", fhPHOSDigits->GetMean(), fhPHOSDigits->GetRMS() ) ;
+ printf(" PHOSRecParticles Mean : %5.3f , RMS : %5.3f \n", fhPHOSRecParticles->GetMean(), fhPHOSRecParticles->GetRMS() ) ;
+ printf(" PHOSPhotons Mean : %5.3f , RMS : %5.3f \n", fhPHOSPhotons->GetMean(), fhPHOSPhotons->GetRMS() ) ;
+ printf(" PHOSInvariantMass Mean : %5.3f , RMS : %5.3f \n", fhPHOSInvariantMass->GetMean(), fhPHOSInvariantMass->GetRMS() ) ;
+ printf(" PHOSDigitsEvent Mean : %5.3f , RMS : %5.3f \n", fhPHOSDigitsEvent->GetMean(), fhPHOSDigitsEvent->GetRMS() ) ;
TCanvas * cPHOS = new TCanvas("cPHOS", "PHOS ESD Test", 400, 10, 600, 700) ;
cPHOS->Divide(3, 2);
cPHOS->cd(1) ;
- gPad->SetLogy();
+ if ( fhPHOSEnergy->GetMaximum() > 0. )
+ gPad->SetLogy();
fhPHOSEnergy->SetAxisRange(0, 25.);
fhPHOSEnergy->SetLineColor(2);
fhPHOSEnergy->Draw();
fhPHOSDigits->Draw();
cPHOS->cd(3) ;
- gPad->SetLogy();
+ if ( fhPHOSRecParticles->GetMaximum() > 0. )
+ gPad->SetLogy();
fhPHOSRecParticles->SetAxisRange(0, 25.);
fhPHOSRecParticles->SetLineColor(2);
fhPHOSRecParticles->Draw();
cPHOS->cd(4) ;
- gPad->SetLogy();
+ if ( fhPHOSPhotons->GetMaximum() > 0. )
+ gPad->SetLogy();
fhPHOSPhotons->SetAxisRange(0,25.);
fhPHOSPhotons->SetLineColor(2);
fhPHOSPhotons->Draw();
fhPHOSInvariantMass->Draw();
cPHOS->cd(6) ;
- gPad->SetLogy();
+ if ( fhPHOSDigitsEvent->GetMaximum() > 0. )
+ gPad->SetLogy();
fhPHOSDigitsEvent->SetAxisRange(0,40.);
fhPHOSDigitsEvent->SetLineColor(2);
fhPHOSDigitsEvent->Draw();
cPHOS->Print("PHOS.eps");
char line[1024] ;
- sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
gROOT->ProcessLine(line);
sprintf(line, ".!rm -fR *.eps");
gROOT->ProcessLine(line);
fhPMDSM24 = (TH2F*)fOutputContainer->At(36);
fhPMDSM = (TH1F*)fOutputContainer->At(37);
+ AliInfo(Form(" *** %s Report:", GetName())) ;
+
gStyle->SetOptStat(110000);
gStyle->SetOptFit(1);
cPMD4->Print("CPVPREAdc.eps");
char line[1024] ;
- sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
gROOT->ProcessLine(line);
AliInfo(Form("!!! All the eps files are in %s.tar.gz !!! \n", GetName())) ;
fhT02 = (TH1F*)fOutputContainer->At(1);
fhT03 = (TH1F*)fOutputContainer->At(2);
+ AliInfo(Form(" *** %s Report:", GetName())) ;
+
Float_t mean = fhT02->GetMean();
printf ("mean time T0 ps %f\n", mean) ;
cTO1->Print("T0.eps");
char line[1024] ;
- sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
gROOT->ProcessLine(line);
sprintf(line, ".!rm -fR *.eps");
gROOT->ProcessLine(line);
fhTOFMatch->Fill(fracM);
PostData(0, fOutputContainer);
- AliInfo("Finishing event processing...") ;
+ // AliInfo("Finishing event processing...") ;
}
//______________________________________________________________________________
fhTOFeffMom->Draw();
fhTOFeffMom->Draw("histo,same");
- cTOFeff->Print("TOF_eff.gif");
cTOFeff->Print("TOF_eff.eps");
TCanvas * cTOFtime = new TCanvas("cTOFtime", "TOF measured Times ", 400, 30, 550, 630) ;
cTOFtime->Divide(1,2) ;
cTOFtime->cd(1);
- cTOFtime->GetPad(1)->SetLogy(1);
+ if ( fhTOFTime->GetMaximum() > 0 )
+ cTOFtime->GetPad(1)->SetLogy(1);
fhTOFTime->GetXaxis()->SetTitle("TOF time (ns)");
fhTOFTime->GetYaxis()->SetTitle("Entries");
fhTOFTime->SetFillColor(4);
fhTOFTime->Draw();
cTOFtime->cd(2);
- cTOFtime->GetPad(2)->SetLogy(1);
+ if ( fhTOFDeltaTime->GetMaximum() > 0 )
+ cTOFtime->GetPad(2)->SetLogy(1);
fhTOFDeltaTime->GetXaxis()->SetTitle("t^{TOF}-t^{exp}_{#pi} (ns)");
fhTOFDeltaTime->GetYaxis()->SetTitle("Entries");
fhTOFDeltaTime->SetFillColor(4);
fhTOFDeltaTime->Draw();
- cTOFtime->Print("TOF_time.gif");
cTOFtime->Print("TOF_time.eps");
TCanvas * cTOFpid = new TCanvas("cTOFpid", "TOF PID ", 400, 30, 550, 630) ;
cTOFpid->Divide(1,3) ;
cTOFpid->cd(1);
- cTOFpid->SetLogy(1);
+ if ( fhTOFMass->GetMaximum() > 0 )
+ cTOFpid->SetLogy(1);
fhTOFMass->GetXaxis()->SetTitle("Reconstructed Mass (GeV/c^{2})");
fhTOFMass->GetYaxis()->SetTitle("Entries");
fhTOFMass->SetFillColor(4);
fhTOFMass->Draw();
- cTOFpid->SetLogy(0);
+ if ( fhTOFMassVsMom->GetMaximum() > 0 )
+ cTOFpid->SetLogy(0);
cTOFpid->cd(2);
fhTOFMassVsMom->GetYaxis()->SetRange(0,400);
fhTOFMassVsMom->GetXaxis()->SetTitle("Reconstructed Mass (GeV/c^{2})");
texpr->Draw();
- cTOFpid->Print("TOF_pid.gif");
cTOFpid->Print("TOF_pid.eps");
// draw all
// some plots
- AliInfo("TOF QA Task: End of events loop");
fOutputContainer = (TObjArray*)GetOutputData(0);
fhTOFMatch = (TH1F*)fOutputContainer->At(0);
fhESDeffPhi = (TH1F*)fOutputContainer->At(1);
fhTOFIDSpecies = (TH1F*)fOutputContainer->At(15);
fhTOFMassVsMom = (TH2F*)fOutputContainer->At(16);
fhTOFMass = (TH1F*)fOutputContainer->At(17);
+
+ AliInfo(Form(" *** %s Report:", GetName())) ;
+
GetEfficiency();
// PostData(0, fOutputContainer);
DrawHistos() ;
+
+ char line[1024] ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
+ gROOT->ProcessLine(line);
+
+ AliInfo(Form("!!! All the eps files are in %s.tar.gz !!! \n", GetName())) ;
}
fOutputContainer->AddAt(fTrdSigMomPID[i], counter++);
}
- AliInfo(Form("Number of histograms = %d", counter));
+ //AliInfo(Form("Number of histograms = %d", counter));
}
void AliTRDQATask::Terminate(Option_t *)
{
// Processing when the event loop is ended
- AliInfo("TRD QA module");
fOutputContainer = (TObjArray*)GetOutputData(0);
int counter = 0;
fNTracks = (TH1D*)fOutputContainer->At(counter++);
}
// create efficiency histograms
+ AliInfo(Form(" *** %s Report:", GetName())) ;
CalculateEff();
-// PostData(0, fOutputContainer);
DrawESD() ;
DrawGeoESD() ;
//DrawConvESD() ;
DrawPidESD() ;
+
+ char line[1024] ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
+ gROOT->ProcessLine(line);
+
+ AliInfo(Form("!!! All the eps files are in %s.tar.gz !!! \n", GetName())) ;
+
}
//______________________________________________________________________________
{
// Makes a few plots
- AliInfo("Plotting....") ;
-
TCanvas * cTRD = new TCanvas("cTRD", "TRD ESD Test", 400, 10, 600, 700) ;
cTRD->Divide(6,3) ;
cTRD->cd(i+1) ;
// new TCanvas(names[i], names[nhist], 500, 300);
- gPad->SetLogy(logy[i]);
for(int j=0; j<nover[i]; j++) {
TH1D *hist = dynamic_cast<TH1D*>(gDirectory->FindObject(names[nhist++]));
if (!hist) continue;
-
+ if (hist->GetMaximum() > 0 )
+ gPad->SetLogy(logy[i]);
if (strstr(hist->GetName(), "eff")) {
hist->SetMarkerStyle(20);
hist->SetMinimum(0);
else hist->Draw("SAME");
}
}
- cTRD->Print("TRD_ESD.gif");
+ cTRD->Print("TRD_ESD.eps");
}
//______________________________________________________________________________
{
// Makes a few plots
- AliInfo("Plotting....") ;
-
TCanvas * cTRDGeo = new TCanvas("cTRDGeo", "TRD ESDGeo Test", 400, 10, 600, 700) ;
cTRDGeo->Divide(4,2) ;
//if (i<2) new TCanvas(names[i], names[i], 500, 300);
//else new TCanvas(names[i], names[i], 300, 900);
- gPad->SetLogy(logy[i]);
+ if (hist->GetMaximum() > 0 )
+ gPad->SetLogy(logy[i]);
if (strstr(opt[i],"colz")) gPad->SetRightMargin(0.1);
hist->Draw(opt[i]);
AliInfo(Form("%s\t%d", names[i], hist->GetEntries()));
}
- cTRDGeo->Print("TRD_Geo.gif");
+ cTRDGeo->Print("TRD_Geo.eps");
}
//______________________________________________________________________________
for(int i=0; i<nplots; i++) {
cTRDConv->cd(i+1) ;
//new TCanvas(names[i], names[i], 500, 300);
- gPad->SetLogy(logy[i]);
if (strstr(opt[i],"colz")) gPad->SetRightMargin(0.1);
for(int j=0; j<nover[i]; j++) {
TH1D *hist = dynamic_cast<TH1D*>(gDirectory->FindObject(names[nhist++]));
+ if ( hist->GetMaximum() > 0 )
+ gPad->SetLogy(logy[i]);
if (!j) hist->Draw(opt[i]);
else hist->Draw("same");
}
{
// Makes a few plots
- AliInfo("Plotting....") ;
TCanvas * cTRDPid = new TCanvas("cTRDPid", "TRD ESDPid Test", 400, 10, 600, 700) ;
cTRDPid->Divide(9,3) ;
if (!hist) continue;
//new TCanvas(names[i], names[i], 500, 300);
- gPad->SetLogy(logy[i]);
+ if ( hist->GetMaximum() > 0 )
+ gPad->SetLogy(logy[i]);
if (strstr(opt[i],"colz")) gPad->SetRightMargin(0.1);
if (strstr(names[i],"sigMom")) gPad->SetLogz(1);
hist->Draw(opt[i]);
AliInfo(Form("%s\t%d", names[i], hist->GetEntries()));
}
- cTRDPid->Print("TRD_Pid.gif");
+ cTRDPid->Print("TRD_Pid.eps");
}
fhVZEROMultA = (TH1I*)fOutputContainer->At(2);
fhVZEROMultC = (TH1I*)fOutputContainer->At(3);
- printf("V0A Multiplicity Mean : %5.3f , RMS : %5.3f \n",fhVZEROMultA->GetMean(),fhVZEROMultA->GetRMS());
- printf("V0C Multiplicity Mean : %5.3f , RMS : %5.3f \n",fhVZEROMultC->GetMean(),fhVZEROMultC->GetRMS());
+ AliInfo(Form(" *** %s Report:", GetName())) ;
+ printf(" V0A Multiplicity Mean : %5.3f , RMS : %5.3f \n",fhVZEROMultA->GetMean(),fhVZEROMultA->GetRMS());
+ printf(" V0C Multiplicity Mean : %5.3f , RMS : %5.3f \n",fhVZEROMultC->GetMean(),fhVZEROMultC->GetRMS());
TCanvas * c1 = new TCanvas("Number of PMs fired in V0A", "Number of PMs fired in V0A", 1);
fhVZERONbPMA->SetAxisRange(0, 99);
c4->Print("NumberV0CPMs.eps");
char line[1024] ;
- sprintf(line, ".!tar -zcvf %s.tar.gz *.eps", GetName()) ;
+ sprintf(line, ".!tar -zcf %s.tar.gz *.eps", GetName()) ;
gROOT->ProcessLine(line);
sprintf(line, ".!rm -fR *.eps");
gROOT->ProcessLine(line);
}
if ( strstr(pararchivename, "ESD") ) {
- gSystem->Load("libVMC.so");
+ //gSystem->Load("libVMC.so");
gSystem->Load("libESD.so");
- gSystem->Load("libRAliEn.so") ;
- gSystem->Load("libProof.so") ;
+ //gSystem->Load("libRAliEn.so") ;
+ //gSystem->Load("libProof.so") ;
}
printf("*** %s library loaded *** %s **\n", pararchivename);
if (! gIsAnalysisLoaded ) {
LoadLib("ESD") ;
LoadLib("ANALYSIS") ;
- printf("Include path = %s\n", gSystem->GetIncludePath()) ;
LoadLib("AnalysisCheck") ;
}
// definition of Tag cuts
const char * runCuts = 0x0 ;
- const char * evtCuts = 0x0 ; //"fEventTag.fNPHOSClustersMin == 1 && fEventTag.fNEMCALClustersMin == 1" ;
+ const char * evtCuts = 0x0 ;
+ const char * lhcCuts = 0x0 ;
+ const char * detCuts = 0x0 ;
+
+//"fEventTag.fNPHOSClustersMin == 1 && fEventTag.fNEMCALClustersMin == 1" ;
TString input = gSystem->Getenv("ANA_INPUT") ;
//create the ESD collection from the tag collection
input.ReplaceAll("tag?", "") ;
const char * collESD = "esdCollection.xml" ;
- ag->MakeEsdCollectionFromTagCollection(runCuts, evtCuts, input.Data(), collESD) ;
+ ag->MakeEsdCollectionFromTagCollection(runCuts, lhcCuts, detCuts, evtCuts, input.Data(), collESD) ;
sprintf(argument, "esd?%s", collESD) ;
} else if ( input.Contains("esd?") )
sprintf(argument, "%s", input.Data()) ;
TChain* analysisChain = new TChain("esdTree") ;
// input = "alien:///alice/cern.ch/user/a/aliprod/prod2006_2/output_pp/105/411/AliESDs.root" ;
// analysisChain->AddFile(input);
- input = "/home/schutz/group/schutz/work/analysis/QA/HEAD/RunV0/AliESDs.root" ;
+ input = "AliESDs.root" ;
analysisChain->AddFile(input);
ag->Process(analysisChain) ;
}