Split: removed dirs now in AliPhysics

[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / ChargedHadrons / dNdPt / macros / plots / readAliceInel.C

diff --git a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/macros/plots/readAliceInel.C b/PWGLF/SPECTRA/ChargedHadrons/dNdPt/macros/plots/readAliceInel.C
deleted file mode 100644 (file)
index 6213cce..0000000
--- a/PWGLF/SPECTRA/ChargedHadrons/dNdPt/macros/plots/readAliceInel.C
+++ /dev/null
@@ -1,130 +0,0 @@
-Int_t readAliceInel() 
-{
-// read ALICE INEL data from file
-//
-
-cout << endl;
-cout << "================================" << endl;
-cout << endl;
-cout << "read ALICE INEL data from file" <<endl;
-cout << "Number of bins: " << binsInelAlice <<endl;
-cout << "Filename:       " << filenameInelAlice <<endl;
-cout << endl;
-cout << "================================" << endl;
-cout << endl;
-
-ifstream fileInelAlice;
-fileInelAlice.open(filenameInelAlice);
-
-Int_t i = 0;
-while(!fileInelAlice.eof()) {
-    if(i == binsInelAlice) break;    
-    fileInelAlice >> centerPtInelAlice[i] >>  inelAlice[i] >> statInelAlice[i] >> systInelAlice[i];
-    //systInelAlice[i] = 0.15*inelAlice[i]+1e-7; // sys error has to be provided, currently 10%
-    
-    // the width of the pt bins, currently hardwired....
-    widthPtInelAlice[i] = 0.05;
-    if (centerPtInelAlice[i] > 1) widthPtInelAlice[i] = 0.1;
-    if (centerPtInelAlice[i] > 2) widthPtInelAlice[i] = 0.2;
-    if (centerPtInelAlice[i] > 4) widthPtInelAlice[i] = 0.5;
-    if (centerPtInelAlice[i] > 7) widthPtInelAlice[i] = 1.0;
-    
-    ptInelAlice[i]       = centerPtInelAlice[i];
-    errPtInelAlice[i]    = widthPtInelAlice[i] / 2.0;
-    lowPtInelAlice[i]    = centerPtInelAlice[i] - errPtInelAlice[i];
-    highPtInelAlice[i]   = centerPtInelAlice[i] + errPtInelAlice[i];      
-            
-    lowStatInelAlice[i]  = inelAlice[i] - statInelAlice[i];
-    highStatInelAlice[i] = inelAlice[i] + statInelAlice[i];
-    relStatInelAlice[i]  = statInelAlice[i] / inelAlice[i];    
-    lowSystInelAlice[i]  = inelAlice[i] - systInelAlice[i];
-    highSystInelAlice[i] = inelAlice[i] + systInelAlice[i];
-    relSystInelAlice[i]  = systInelAlice[i] / inelAlice[i];    
-    errInelAlice[i]      = systInelAlice[i] + statInelAlice[i];
-    lowErrInelAlice[i]   = inelAlice[i] - errInelAlice[i];
-    highErrInelAlice[i]  = inelAlice[i] + errInelAlice[i];
-    relErrInelAlice[i]   = errInelAlice[i] / inelAlice[i];
-    err2InelAlice[i]     = sqrt(systInelAlice[i]*systInelAlice[i] + statInelAlice[i]*statInelAlice[i]);
-    lowErr2InelAlice[i]  = inelAlice[i] - err2InelAlice[i];
-    highErr2InelAlice[i] = inelAlice[i] + errInelAlice[i];
-    relErr2InelAlice[i]  = err2InelAlice[i] / inelAlice[i];
-    
-    
-    
-    ptInel2PiPtAlice[i]       = ptInelAlice[i];
-    centerPtInel2PiPtAlice[i] = centerPtInelAlice[i];
-    widthPtInel2PiPtAlice[i]  = widthPtInelAlice[i];
-    errPtInel2PiPtAlice[i]    = errPtInelAlice[i];
-    lowPtInel2PiPtAlice[i]    = lowPtInelAlice[i];
-    highPtInel2PiPtAlice[i]   = highPtInelAlice[i];
-
-    inel2PiPtAlice[i]         = inelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    statInel2PiPtAlice[i]     = statInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    lowStatInel2PiPtAlice[i]  = lowStatInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    highStatInel2PiPtAlice[i] = highStatInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    relStatInel2PiPtAlice[i]  = relStatInelAlice[i];
-    systInel2PiPtAlice[i]     = systInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    lowSystInel2PiPtAlice[i]  = lowSystInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    highSystInel2PiPtAlice[i] = highSystInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    relSystInel2PiPtAlice[i]  = relSystInelAlice[i];
-    errInel2PiPtAlice[i]      = errInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    lowErrInel2PiPtAlice[i]   = lowErrInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    highErrInel2PiPtAlice[i]  = highErrInelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    relErrInel2PiPtAlice[i]   = relErrInelAlice[i];
-    err2Inel2PiPtAlice[i]     = err2InelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    lowErr2Inel2PiPtAlice[i]  = lowErr2InelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    highErr2Inel2PiPtAlice[i] = highErr2InelAlice[i]*centerPtInelAlice[i]*2*M_PI;
-    relErr2Inel2PiPtAlice[i]  = relErr2InelAlice[i];
-    
-        
-    
-    cout << "ptInelAlice[" << i << "]          = " << ptInelAlice[i] <<endl;
-    cout << "   centerPtInelAlice[" << i << "] = " << centerPtInelAlice[i] <<endl;
-    cout << "   widthPtInelAlice[" << i << "]  = " << widthPtInelAlice[i] <<endl;
-    cout << "   errPtInelAlice[" << i << "]    = " << errPtInelAlice[i] <<endl;
-    cout << "   lowPtInelAlice[" << i << "]    = " << lowPtInelAlice[i] <<endl;
-    cout << "   highPtInelAlice[" << i << "]   = " << highPtInelAlice[i] <<endl;
-    cout << "inelAlice[" << i << "]            = " << inelAlice[i] <<endl;
-    cout << "   statInelAlice[" << i << "]     = " << statInelAlice[i] <<endl;
-    cout << "   lowStatInelAlice[" << i << "]  = " << lowStatInelAlice[i] <<endl;
-    cout << "   highStatInelAlice[" << i << "] = " << highStatInelAlice[i] <<endl;
-    cout << "   relStatInelAlice[" << i << "]  = " << relStatInelAlice[i] <<endl;
-    cout << "   systInelAlice[" << i << "]     = " << systInelAlice[i] <<endl;
-    cout << "   lowSystInelAlice[" << i << "]  = " << lowSystInelAlice[i] <<endl;
-    cout << "   highSystInelAlice[" << i << "] = " << highSystInelAlice[i] <<endl;
-    cout << "   relSystInelAlice[" << i << "]  = " << relSystInelAlice[i] <<endl;
-    cout << "errInelAlice[" << i << "]         = " << errInelAlice[i] <<endl;
-    cout << "   lowErrInelAlice[" << i << "]   = " << lowErrInelAlice[i] <<endl;
-    cout << "   highErrInelAlice[" << i << "]  = " << highErrInelAlice[i] <<endl;
-    cout << "   relErrInelAlice[" << i << "]   = " << relErrInelAlice[i] <<endl;
-    cout << "err2InelAlice[" << i << "]        = " << err2InelAlice[i] <<endl;
-    cout << "   lowErr2InelAlice[" << i << "]  = " << lowErr2InelAlice[i] <<endl;
-    cout << "   highErr2InelAlice[" << i << "] = " << highErr2InelAlice[i] <<endl;
-    cout << "   relErr2InelAlice[" << i << "]  = " << relErr2InelAlice[i] <<endl;
-    cout << endl;
-    /*
-    ptInelAliceFit[i]      = ptInelAlice[i];
-    InelInelAliceFit[i]     = InelInelAlice[i]*ptInelAlice[i];
-    errInelInelAliceFit[i]  = errInelInelAlice[i];
-    widthPtInelAliceFit[i] = widthPtInelAlice[i];
-    errPtInelAliceFit[i]  = 0; //errPtInelAlice[i];
-    lowPtInelAliceFit[i]  = lowPtInelAlice[i];
-    highPtInelAliceFit[i] = highPtInelAlice[i];
-    */
-        
-   i++;
-} // while(!fileInelAlice.eof())
-fileInelAlice.close();
-//if (fileInelAlice) { delete fileInelAlice; }
-//fileInelAlice=0;
-
-cout << "================================" << endl;
-cout << endl;
-cout << "Finished reading ALICE INEL data" <<endl;
-cout << "Number of bins read: " << i <<endl;
-cout << endl;
-cout << "================================" << endl;
-cout << endl;
-
-return i;
-}
\ No newline at end of file